Using Gaussian accelerated molecular dynamics combined with Markov state models to explore the mechanism of action of new oral inhibitors on Complex I.

In this study, our focus was on investigating H-1,2,3-triazole derivative HP661 as a novel and highly efficient oral OXPHOS inhibitor, with its molecular-level inhibitory mechanism not yet fully understood. We selected the ND1, NDUFS2, and NDUFS7 subunits of Mitochondrial Complex I as the receptor proteins and established three systems for comparative analysis: protein-IACS-010759, protein-lead compound 10, and protein-HP661. Through extensive analysis involving 500 ns Gaussian molecular dynamics simulations, we gained insights into these systems. Additionally, we constructed a Markov State Models to examine changes in secondary structures during the motion processes. The research findings suggest that the inhibitor HP661 enhances the extensibility and hydrophilicity of the receptor protein. Furthermore, HP661 induces the unwinding of the α-helical structure in the region of residues 726-730. Notably, key roles were identified for Met37, Phe53, and Pro212 in the binding of various inhibitors. In conclusion, we delved into the potential molecular mechanisms of triazole derivative HP661 in inhibiting Complex I. These research outcomes provide crucial information for a deeper understanding of the mechanisms underlying OXPHOS inhibition, offering valuable theoretical support for drug development and disease treatment design.

Elucidating the structural basis of vitamin B12 derivatives as novel potent inhibitors of PTP1B: Insights from inhibitory mechanisms using Gaussian accelerated molecular dynamics (GaMD) and in vitro study.

Vitamin B12 is a group of biologically active cobalamin compounds. In this study, we investigated the inhibitory effects of methylcobalamin (MeCbl) and hydroxocobalamin acetate (OHCbl Acetate) on protein tyrosine phosphatase 1B (PTP1B). MeCbl and OHCbl Acetate exhibited an IC50 of approximately 58.390 ± 2.811 µM and 8.998 ± 0.587 µM, respectively. The Ki values of MeCbl and OHCbl Acetate were 25.01 µM and 4.04 µM respectively. To elucidate the inhibition mechanism, we conducted a 500 ns Gaussian accelerated molecular dynamics (GaMD) simulation. Utilizing PCA and tICA, we constructed Markov state models (MSM) to examine secondary structure changes during motion. Our findings revealed that the α-helix at residues 37-42 remained the most stable in the PTP1B-OHCbl Acetate system. Furthermore, upon binding of OHCbl Acetate or MeCbl, the WPD loop of PTP1B moved inward to the active pocket, forming a closed conformation and potentially obstructs substrate entry. Protein-ligand interaction analysis and MM-PBSA showed that OHCbl Acetate exhibited lower binding free energy and engaged in more residue interactions with PTP1B. In summary, our study confirmed the substantial inhibitory activity of OHCbl Acetate against PTP1B, with its inhibitory potency notably surpassing that of MeCbl. We demonstrated potential molecular mechanisms of OHCbl Acetate inhibiting PTP1B.

Building a Kokumi Database and Machine Learning-Based Prediction: A Systematic Computational Study on Kokumi Analysis.

Kokumi is a subtle sensation characterized by a sense of fullness, continuity, and thickness. Traditional methods of taste discovery and analysis, including those of kokumi, have been labor-intensive and costly, thus necessitating the emergence of computational methods as critical strategies in molecular taste analysis and prediction. In this study, we undertook a comprehensive analysis, prediction, and screening of the kokumi compounds. We categorized 285 kokumi compounds from a previously unreleased kokumi database into five groups based on their molecular characteristics. Moreover, we predicted kokumi/non-kokumi and multi-flavor compositions using six structure-taste relationship models: MLP-E3FP, MLP-PLIF, MLP-RDKFP, SVM-RDKFP, RF-RDKFP, and WeaveGNN feature of Atoms and Bonds. These six predictors exhibited diverse performance levels across two different models. For kokumi/non-kokumi prediction, the WeaveGNN model showed an exceptional predictive AUC value (0.94), outperforming the other models (0.87, 0.90, 0.89, 0.92, and 0.78). For multi-flavor prediction, the MLP-E3FP model demonstrated a higher predictive AUC and MCC value (0.94 and 0.74) than the others (0.73 and 0.33; 0.92 and 0.70; 0.95 and 0.73; 0.94 and 0.64; and 0.88 and 0.69). This data highlights the model's proficiency in accurately predicting kokumi molecules. As a result, we sourced kokumi active compounds through a high-throughput screening of over 100 million molecules, further refined by toxicity and similarity screening. Lastly, we launched a web platform, KokumiPD (https://www.kokumipd.com/), offering a comprehensive kokumi database and online prediction services for users.

Probing the mechanisms of hydrazide-based HDAC inhibitors binding to HDAC3 using Gaussian accelerated molecular dynamics (GaMD) simulations.

Histone deacetylases (HDACs) have emerged as promising targets for anticancer drug development. They regulate gene expression by removing acetyl groups from lysine residues on histone tails, leading to chromatin condensation. A hydrazide-based HDAC inhibitor, N-(4-(2-Propylhydrazine-1-carbonyl)benzyl)-1H-indole-2-carboxamide (11h), has been reported to exhibit significant in vivo antitumor activity. In comparison to the lead compound N-(4-(2-Propylhydrazine-1-carbonyl)benzyl)cinnamamide (17), compound 11h demonstrates 2- to 5-fold higher HDAC inhibition and cell-based antitumor activity. However, the inhibitory mechanism of 11h remains insufficiently explored. In this study, we conducted 500 ns Gaussian Accelerated Molecular Dynamics (GaMD) simulations on Histone deacetylase 3 (HDAC3) and two complex systems (HDAC3-17 and HDAC3-11h). Our findings revealed that upon inhibitor binding, the active pocket volume of HDAC3 undergone alterations, and the movement of the L6-loop toward the active site impeded substrate entry. Moreover, we observed a destabilization of the α-helix in the aa75-89 region of HDAC3 compared to the two complex systems, indicating partial unwinding. Notably, 11h exhibited a closer proximity of its carbonyl oxygen to the active pocket's Zn2+ metal compared to 17, increasing the likelihood of coordination with the Zn2+ metal. The analysis of protein-ligand interactions highlighted a greater number of hydrogen bonds and other interactions between 11h and the receptor protein when compared to 17, underscoring the stronger binding of 11h to HDAC3. In conclusion, our study provided theoretical insights into the inhibitory mechanism of hydrazide-based HDAC inhibitors on HDAC3, thereby contributing to the development of improved drug targets for cancer therapy.Communicated by Ramaswamy H. Sarma.

Molecular mechanisms of metal ions in regulating the catalytic efficiency of D-psicose 3-epimerase revealed by multiple short molecular dynamic simulations and free energy predictions.

Communicated by Ramaswamy H. Sarma.

Computational Study on the Effect of Inactivating/Activating Mutations on the Inhibition of MEK1 by Trametinib.

Activation of the mitogen-activated protein kinase (MAPK) signaling pathway regulated by human MAP kinase 1 (MEK1) is associated with the carcinogenesis and progression of numerous cancers. In addition, two active mutations (P124S and E203K) have been reported to enhance the activity of MEK1, thereby eventually leading to the tumorigenesis of cancer. Trametinib is an MEK1 inhibitor for treating EML4-ALK-positive, EGFR-activated, and KRAS-mutant lung cancers. Therefore, in this study, molecular docking and molecular dynamic (MD) simulations were performed to explore the effects of inactive/active mutations (A52V/P124S and E203K) on the conformational changes of MEK1 and the changes in the interaction of MEK1 with trametinib. Moreover, steered molecular dynamic (SMD) simulations were further utilized to compare the dissociation processes of trametinib from the wild-type (WT) MEK1 and two active mutants (P124S and E203K). As a result, trametinib had stronger interactions with the non-active MEK1 (WT and A52V mutant) than the two active mutants (P124S and E203K). Moreover, two active mutants may make the allosteric channel of MEK1 wider and shorter than that of the non-active types (WT and A52V mutant). Hence, trametinib could dissociate from the active mutants (P124S and E203K) more easily compared with the WT MEK1. In summary, our theoretical results demonstrated that the active mutations may attenuate the inhibitory effects of MEK inhibitor (trametinib) on MEK1, which could be crucial clues for future anti-cancer treatment.