Cytochrome c (CytC), a one-electron carrier, transfers electrons from complex bc1 to cytochrome c oxidase (CcO) in the electron-transport chain. Electrostatic interaction with the partners, complex bc1 and CcO, is ensured by a lysine cluster near the heme forming the Universal Binding Site (UBS). We constructed three mutant variants of mitochondrial CytC with one (2Mut), four (5Mut), and five (8Mut) Lys->Glu substitutions in the UBS and some compensating Glu->Lys substitutions at the periphery of the UBS for charge compensation. All mutants showed a 4-6 times increased peroxidase activity and accelerated binding of cyanide to the ferric heme of CytC. In contrast, decomposition of the cyanide complex with ferrous CytC, as monitored by magnetic circular dichroism spectroscopy, was slower in mutants compared to WT. Molecular dynamic simulations revealed the increase in the fluctuations of Cα atoms of individual residues of mutant CytC compared to WT, especially in the Ω-loop (70-85), which can cause destabilization of the Fe…S(Met80) coordination link, facilitation of the binding of exogenous ligands cyanide and peroxide, and an increase in peroxidase activity. It was found that only one substitution K72E is enough to induce all these changes, indicating the significance of K72 and the Ω-loop (70-85) for the structure and physiology of mitochondrial CytC. In this work, we also propose using a ferro-ferricyanide buffer as a substrate to monitor the peroxidase activity of CytC. This new approach allows us to determine the rate of peroxidase activity at moderate (200 µM) concentrations of H2O2 and avoid complications of radical formation during the reaction.

β-Glucuronidase, a crucial enzyme in drug metabolism and detoxification, represents a promising target for therapeutic intervention due to its potential to modulate drug pharmacokinetics and enhance therapeutic efficacy. Herein, we assessed the inhibitory potential of phytochemicals from Hibiscus trionum against β-glucuronidase. Grossamide and grossamide K emerged as the most potent β-glucuronidase inhibitors with IC50 values of 0.73 ± 0.03 and 1.24 ± 0.03 μM, respectively. The investigated alkaloids effectively inhibited β-glucuronidase-catalyzed PNPG hydrolysis through a noncompetitive inhibition mode, whereas steppogenin displayed a mixed inhibition mechanism. Molecular docking analyses highlighted grossamide and grossamide K as inhibitors with the lowest binding free energy, all compounds successfully docked into the same main binding site occupied by the reference drug Epigallocatechin gallate (EGCG). We explored the interaction dynamics of isolated compounds with β-glucuronidase through a 200 ns molecular dynamics (MD) simulation. Analysis of various MD parameters revealed that grossamide and grossamide K maintained stable trajectories and demonstrated significant energy stabilization upon binding to β-glucuronidase. Additionally, these compounds exhibited the lowest average interaction energies with the target enzyme. The MM/PBSA calculations further supported these findings, showing the lowest binding free energies for grossamide and grossamide K. These computational results are consistent with experimental data, suggesting that grossamide and grossamide K could be potent inhibitors of β-glucuronidase.

Unraveling the intricacies of β-glucuronidase inhibition is pivotal for developing effective strategies in applications specific to gastrointestinal health and drug metabolism. Our study investigated the efficacy of some Hibiscus trionum phytochemicals as β-glucuronidase inhibitors. The results showed that cleomiscosin A and mansonone H emerged as the most potent inhibitors, with IC50 values of 3.97 ± 0.35 μM and 10.32 ± 1.85 μM, respectively. Mechanistic analysis of β-glucuronidase inhibition indicated that cleomiscosin A and the reference drug EGCG displayed a mixed inhibition mode against β-glucuronidase, while mansonone H exhibited noncompetitive inhibition against β-glucuronidase. Docking studies revealed that cleomiscosin A and mansonone H exhibited the lowest binding affinities, occupying the same site as EGCG, and engaged significant key residues in their binding mechanisms. Using a 30 ns molecular dynamics (MD) simulation, we explored the interaction dynamics of isolated compounds with β-glucuronidase. Analysis of various MD parameters showed that cleomiscosin A and mansonone H exhibited consistent trajectories and significant energy stabilization with β-glucuronidase. These computational insights complemented experimental findings, underscoring the potential of cleomiscosin A and mansonone H as β-glucuronidase inhibitors.

SHP2 belongs to a cytoplasmic non-receptor protein tyrosine phosphatase class. It plays a critical role in the development of various cancers, such as gastric cancer, leukemia, and breast cancer. Thus, SHP2 has gained the interest of researchers as a potential target for inhibiting tumor cell proliferation in SHP2-dependent cancers. This study employed pharmacophore-based virtual screening, molecular docking, molecular dynamic (MD) simulations, MM/PBSA, and principal component analysis (PCA), followed by ADME prediction. We selected three potential hits from a collective database of more than one million chemical compounds. The stability of these selected hit-protein complexes was analyzed using 500 ns MD simulations and binding free energy calculations. The identified hits Lig_1, Lig_6, and Lig_14 demonstrated binding free energies of -161.49 kJ/mol, -151.28 kJ/mol, and -107.13 kJ/mol, respectively, compared to the reference molecule (SHP099) with a ΔG of -71.48 kJ/mol. Our results showed that the identified compounds could be used as promising candidates for selective SHP2 allosteric inhibition in cancer.

Alzheimer\'s disease (AD) poses a significant global health challenge, necessitating the exploration of novel therapeutic strategies. Fyn Tyrosine Kinase has emerged as a key player in AD pathogenesis, making it an attractive target for drug development. This study focuses on investigating the potential of Papaveroline as a drug candidate for AD by targeting Fyn Tyrosine Kinase. The research employed high-throughput virtual screening and QSAR analysis were conducted to identify compounds with optimal drug-like properties, emphasizing adherence to ADMET parameters for further evaluation. Molecular dynamics simulations to analyze the binding interactions between Papaveroline and Staurosporine with Fyn Tyrosine Kinase over a 200-ns period. The study revealed detailed insights into the binding mechanisms and stability of the Papaveroline-Fyn complex, showcasing the compound\'s potential as an inhibitor of Fyn Tyrosine Kinase. Comparative analysis with natural compounds and a reference compound highlighted Papaveroline\'s unique characteristics and promising therapeutic implications for AD treatment. Overall, the findings underscore Papaveroline\'s potential as a valuable drug candidate for targeting Fyn Tyrosine Kinase in AD therapy, offering new avenues for drug discovery in neurodegenerative diseases. This study contributes to advancing our understanding of molecular interactions in AD pathogenesis and paves the way for further research and development in this critical area.

Prolyl oligopeptidase (POP) is a compelling therapeutic target associated with aging and neurodegenerative disorders due to its pivotal role in neuropeptide processing. Despite initial promise demonstrated by early-stage POP inhibitors, their progress in clinical trials has been halted at Phase I or II. This impediment has prompted the pursuit of novel inhibitors. The current study seeks to contribute to the identification of efficacious POP inhibitors through the design, synthesis, and comprehensive evaluation (both in vitro and in silico) of thiazolyl thiourea derivatives (5a-r). In vitro experimentation exhibited that the compounds displayed significant higher potency as POP inhibitors. Compound 5e demonstrated an IC50 value of 16.47 ± 0.54 μM, representing a remarkable potency. A meticulous examination of the structure-activity relationship indicated that halogen and methoxy substituents were the most efficacious. In silico investigations delved into induced fit docking, pharmacokinetics, and molecular dynamics simulations to elucidate the intricate interactions, orientation, and conformational changes of these compounds within the active site of the enzyme. Moreover, our pharmacokinetic assessments confirmed that the majority of the synthesized compounds possess attributes conducive to potential drug development.

Alpha-fetoprotein (AFP) is a glycoprotein primarily expressed during embryogenesis, with declining levels postnatally. Elevated AFP levels correlate with pathological conditions such as liver fibrosis, cirrhosis, and hepatocellular carcinoma (HCC). Recent investigations underscore AFP\'s intracellular role in HCC progression, wherein it forms complexes with proteins like Phosphatase and tensin homolog (PTEN), Caspase 3 (CASP3), and Retinoic acid receptors and Retinoid X receptors (RAR/RXR). RAR and RXR regulate gene expression linked to cell death and tumorigenesis in normal physiology. AFP impedes RAR/RXR dimerization, nuclear translocation, and function, promoting gene expression favoring cancer progression in HCC that provoked us to target AFP as a drug candidate. Despite extensive studies, inhibitors targeting AFP to disrupt complex formation and activities remain scarce. In this study, employing protein-protein docking, amino acid residues involved in AFP-RARβ interaction were identified, guiding the definition of AFP\'s active site for potential inhibitor screening. Currently, kinase inhibitors play a significant role in cancer treatment and, the present study explores the potential of repurposing FDA-approved protein kinase inhibitors to target AFP. Molecular docking with kinase inhibitors revealed Lapatinib as a candidate drug of the AFP-RARβ complex. Molecular dynamics simulations and binding energy calculations, employing Mechanic/Poisson-Boltzmann Surface Area (MM-PBSA), confirmed Lapatinib\'s stability with AFP. The study suggests Lapatinib\'s potential in disrupting the AFP-RARβ complex, providing a promising avenue for treating molecularly stratified AFP-positive HCC or its early stages.

Tumor suppressor cylindromatosis protein (CYLD) regulates NF-κB and JNK signaling pathways by cleaving K63-linked poly-ubiquitin chain from its substrate molecules and thus preventing the progression of tumorigenesis and metastasis of the cancer cells. Mutations in CYLD can cause aberrant structure and abnormal functionality leading to tumor formation. In this study, we utilized several computational tools such as PANTHER, PROVEAN, PredictSNP, PolyPhen-2, PhD-SNP, PON-P2, and SIFT to find out deleterious nsSNPs. We also highlighted the damaging impact of those deleterious nsSNPs on the structure and function of the CYLD utilizing ConSurf, I-Mutant, SDM, Phyre2, HOPE, Swiss-PdbViewer, and Mutation 3D. We shortlisted 18 high-risk nsSNPs from a total of 446 nsSNPs recorded in the NCBI database. Based on the conservation profile, stability status, and structural impact analysis, we finalized 13 nsSNPs. Molecular docking analysis and molecular dynamic simulation concluded the study with the findings of two significant nsSNPs (R830K, H827R) which have a remarkable impact on binding affinity, RMSD, RMSF, radius of gyration, and hydrogen bond formation during CYLD-ubiquitin interaction. The principal component analysis compared native and two mutants R830K and H827R of CYLD that signify structural and energy profile fluctuations during molecular dynamic (MD) simulation. Finally, the protein-protein interaction network showed CYLD interacts with 20 proteins involved in several biological pathways that mutations can impair. Considering all these in silico analyses, our study recommended conducting large-scale association studies of nsSNPs of CYLD with cancer as well as designing precise medications against diseases associated with these polymorphisms.

Isoprenoid glycerol dialkyl glycerol tetraethers (GDGTs), characteristic membrane lipids of archaea, are widely used in ecological and geochemical studies, especially for paleoenvironmental reconstruction. Glycerol monoalkyl glycerol tetraethers (GMGTs, also known as H-GDGTs), a unique variant of GDGTs, have covalent bonds linking the two alkyl chains. Despite some studies suggesting a link between GMGTs and high temperatures, the reliability and mechanisms remain unclear. Using molecular dynamics simulations, we elucidated the mechanism connecting GMGTs to high temperatures. Our findings show that H-bridging linkages reduce the distance between alkyl chains, leading to thicker and denser membranes with lower fluidity and permeability. The diffusion coefficient of GMGTs decreased by approximately 35 % compared to GDGTs, indicating their role as a archaeal high-temperature adaptation. This study provides a mechanistic basis for using archaeal GMGTs in geochemical studies and enhances confidence in their use for paleotemperature reconstruction.

A cell\'s ability to survive and to evade cancer is contingent on its ability to retain genomic integrity, which can be seriously compromised when nucleic acid phosphodiester bonds are disrupted. DNA Ligase 1 (LIG1) plays a key role in genome maintenance by sealing single-stranded nicks that are produced during DNA replication and repair. Autosomal recessive mutations in a limited number of individuals have been previously described for this gene. Here we report a homozygous LIG1 mutation (p.A624T), affecting a universally conserved residue, in a patient presenting with leukopenia, neutropenia, lymphopenia, pan-hypogammaglobulinemia, and diminished in vitro response to mitogen stimulation. Patient fibroblasts expressed normal levels of LIG1 protein but exhibited impaired growth, poor viability, high baseline levels of gamma-H2AX foci, and an enhanced susceptibility to DNA-damaging agents. The mutation reduced LIG1 activity by lowering its affinity for magnesium 2.5-fold. Remarkably, it also increased LIG1 fidelity > 50-fold against 3\' end 8-Oxoguanine mismatches, exhibiting a marked reduction in its ability to process such nicks. This is expected to yield increased ss- and dsDNA breaks. Molecular dynamic simulations, and Residue Interaction Network studies, predicted an allosteric effect for this mutation on the protein loops associated with the LIG1 high-fidelity magnesium, as well as on DNA binding within the adenylation domain. These dual alterations of suppressed activity and enhanced fidelity, arising from a single mutation, underscore the mechanistic picture of how a LIG1 defect can lead to severe immunological disease.