Molecular Mechanics Generalized Born Surface Area Research Articles

Network pharmacology, molecular simulation, and binding free energy calculation-based investigation of Neosetophomone B revealed key targets for the treatment of cancer.

In the current study, Neosetophomone B (NSP-B) was investigated for its anti-cancerous potential using network pharmacology, quantum polarized ligand docking, molecular simulation, and binding free energy calculation. Using SwissTarget prediction, and Superpred, the molecular targets for NSP-B were predicted while cancer-associated genes were obtained from DisGeNet. Among the total predicted proteins, only 25 were reported to overlap with the disease-associated genes. A protein-protein interaction network was constructed by using Cytoscape and STRING databases. MCODE was used to detect the densely connected subnetworks which revealed three sub-clusters. Cytohubba predicted four targets, i.e., fibroblast growth factor , FGF20, FGF22, and FGF23 as hub genes. Molecular docking of NSP-B based on a quantum-polarized docking approach with FGF6, FGF20, FGF22, and FGF23 revealed stronger interactions with the key hotspot residues. Moreover, molecular simulation revealed a stable dynamic behavior, good structural packing, and residues' flexibility of each complex. Hydrogen bonding in each complex was also observed to be above the minimum. In addition, the binding free energy was calculated using the MM/GBSA (Molecular Mechanics/Generalized Born Surface Area) and MM/PBSA (Molecular Mechanics/Poisson-Boltzmann Surface Area) approaches. The total binding free energy calculated using the MM/GBSA approach revealed values of -36.85kcal/mol for the FGF6-NSP-B complex, -43.87kcal/mol for the FGF20-NSP-B complex, and -37.42kcal/mol for the FGF22-NSP-B complex, and -41.91kcal/mol for the FGF23-NSP-B complex. The total binding free energy calculated using the MM/PBSA approach showed values of -30.05kcal/mol for the FGF6-NSP-B complex, -39.62kcal/mol for the FGF20-NSP-B complex, -34.89kcal/mol for the FGF22-NSP-B complex, and -37.18kcal/mol for the FGF23-NSP-B complex. These findings underscore the promising potential of NSP-B against FGF6, FGF20, FGF22, and FGF23, which are reported to be essential for cancer signaling. These results significantly bolster the potential of NSP-B as a promising candidate for cancer therapy.

PHARMACOINFORMATICS ANALYSIS OF MORUS MACROURA FOR DRUG DISCOVERY AND DEVELOPMENT

Objective: Pharmacoinformatics is an innovative approach rapidly evolving in pharmaceutical research and drug development. This study focuses on analysing Morus macroura, a plant species with untapped pharmacological potential. This investigation aims to leverage pharmacoinformatics techniques to unveil the hidden potential of Morus macroura in drug discovery and development.
 Methods: The study includes analyses of protein-protein interactions, deep learning docking, adsorption tests, distribution, metabolism, excretion, molecular dynamics simulations and free energy calculation using Molecular Mechanics Generalized Born Surface Area (MMGBSA).
 Results: Nine active compounds were identified in Morus macroura, namely Andalasin A, Guangsangon K, Guangsangon L, Guangsangon M, Guangsangon N, Macrourone C, Mulberrofuran G, Mulberrofuran K, and Mulberroside C. These compounds exhibit protein-protein interaction activities against a cytochrome P450 monooxygenase that catalyses the conversion of C19 androgens. These plant compounds influence aromatase excess syndrome, deficiency, and ovarian dysgenesis. Regarding drug-likeness, Mulberroside C and Macrourone C demonstrated good absorption potential by adhering to Lipinski's rule of five. Deep learning docking simulations yielded affinity results of-9.62 kcal/mol for Guangsangon M,-10.44 kcal/mol for Macrourone C, and-10.99 kcal/mol for Guangsangon L. Subsequent molecular dynamics simulations indicated that Guangsangon L and Macrourone C remained stable during a 100 ns simulation.
 Conclusion: Morus macroura interacts with important proteins, particularly CYP19A1, which might influence health conditions like aromatase excess syndrome and ovarian dysgenesis. These findings provide potential paths for addressing specific health issues and advancing drug development. Molecular dynamics simulations indicated that Guangsangon L and Macrourone C remained stable during simulation.

The Interaction Mechanism of Picolinamide Fungicide Targeting on the Cytochrome bc1 Complex and Its Structural Modification.

Picolinamide fungicides, structurally related to UK-2A and antimycin-A, bind into the Qi-site in the bc1 complex. However, the detailed binding mode of picolinamide fungicides remains unknown. In the present study, antimycin-A and UK-2A were selected to study the binding mode of picolinamide inhibitors with four protonation states in the Qi-site by integrating molecular dynamics simulation, molecular docking, and molecular mechanics Generalized Born surface area (MM/GBSA) calculations. Subsequently, a series of new picolinamide derivatives were designed and synthesized to further understand the effects of substituents on the tail phenyl ring. The computational results indicated that the substituted aromatic rings in antimycin-A and UK-2A were the pharmacophore fragments and made the primary contribution when bound to a protein. Compound 9g-hydrolysis formed H-bonds with Hie201 and Ash228 and showed an IC50 value of 6.05 ± 0.24 μM against the porcine bc1 complex. Compound 9c, with a simpler chemical structure, showed higher control effects than florylpicoxamid against cucumber downy mildew and expanded the fungicidal spectrum of picolinamide fungicides. The structural and mechanistic insights obtained from the present study will provide a valuable clue for the future designing of new promising Qi-site inhibitors.

Bioactivity and molecular docking studies of selected plant compounds.

The shortcomings of synthetic pesticides are becoming increasingly evident. In the pursuit of looking for alternative, environmentally friendly pesticides, the potential use of two naturally occurring sesquiterpenes, α-bisabolol and bisabolol oxide A, as bactericides or acaricides, was investigated against three gram-positive and three gram-negative bacteria, and the eggs and larvae of the cattle tick Hyalomma scupense, respectively. The primary targeted site of synthetic insecticides is acetylcholinesterase (AChE), so the inhibitory activity of both sesquiterpenes was evaluated. Then, molecular docking and molecular mechanics generalized Born surface area (MM-GBSA)studies were performed to gain insight into their binding mode in the AChE active site. Bisabolol oxide A showed good activity against Staphylococcus aureus, Escherichia coli, and Salmonella enteritidis with inhibition zones ranging from 16 to 34 mm and minimal inhibitory concentration (CIM) of 3.125 mg mL-1. Ovicidal and larvicidal assays revealed promising activity of α-bisabolol and bisabolol oxide A against H. scupense, with bisabolol oxide A being more effective against larval ticks with lethal concentration (LC50) value of 0.78%. Both compounds also displayed potent anti-AChE activity with inhibition concentration (IC50) values of 37.09 and 28.14 μg mL-1, respectively. Furthermore, α-bisabolol and bisabolol oxide A exhibited good and comparable docking scores (-7.289 and -7.139 Kcal mol-1, respectively) and were found to accommodate in the active-site gorge of AChE via hydrogen bonding and hydrophobic interactions. Bisabolol oxide A and, to a lesser degree, α-bisabolol are active against bacteria and ectoparasites and may represent an economical and sustainable alternative to toxic synthetic pesticides to control pathogens. © 2024 Society of Chemical Industry.

The effect of phosphorylation on the conformational dynamics and allostery of the association of death-associated protein kinase with calmodulin

Protein phosphorylation plays an important role in the signal transduction and is capable of regulation of cell activity. The death-associated protein kinase 1 (DAPK1), as a Ser/Thr kinase, interacts with calmodulin (CaM) to regulate apoptotic and autophagic signaling. Autophosphorylation of DAPK1 at Ser308 located at the autoregulatory domain (ARD) blocks CaM binding and inhibits kinase catalytic activity. However, the mechanism underlying the influence of Ser308 phosphorylation (pS308) on the DAPK1 activity remains unclear. Here, we performed multiple, microsecond length molecular dynamics (MD) simulations, the molecular mechanics generalized Born/surface area (MM-GBSA) binding free energy calculations, principal component analysis, and dynamic cross-correlation analysis to unravel the conformational dynamics and allostery of the DAPK1 − CaM interaction triggered by the pS308 at the ARD. MD simulations showed that pS308 affected the conformational stability of the DAPK1 − CaM complex. Further energetic and structural exploration revealed that pS308 weakened the association of the phosphorylated DAPK1 to CaM, which lowered the susceptibility of DAPK1 to be activated by CaM. This result can provide mechanistic insights into the molecular underpinning through which the DAPK1 kinase activity is modulated by the auto-phosphorylation. Communicated by Ramaswamy H. Sarma

In silico and in vitro screening of pyrrole-based Hydrazide-Hydrazones as novel acetylcholinesterase inhibitors

Virtual screening is emerging as a highly applied technique and gained prominence as widely used method for the search and identification of potential hits, significantly reducing the time needed to discover novel and effective compounds compared to high-throughput screening. Recently, the superiority of simulations with multiple programs compared to a single software docking has been discussed. The aim of this work was to apply consensus docking, molecular mechanics/generalized Born surface area (MM/GBSA) free binding energy recalculations, and in vitro evaluations on an in-house dataset of recently synthesized pyrrole-based hydrazide-hydrazones in the search for novel acetylcholinesterase (AChE) inhibitors. Two licensed softwares – GOLD 5.3 and Glide, were employed for the virtual screenings, and several chemotherapeutic potential hits were identified. Furthermore, MM/GBSA free binding energy recalculations were provided to enhance the robustness of the in silico results. The MM/GBSA scores of the top ten pyrrole-based hydrazide-hydrazones were ranging from -60.44 to -70.93 Kcal/mol. Subsequent, in vitro evaluations of the top ranked compounds revealed that 12d exhibited the highest AChE inhibitory activity, with a 55% inhibition rate at a concentration of 10 μM. Moreover, this prominent pyrrole-based AChE inhibitor formed stable complex with the active site of the enzyme. Interactions with the active amino residues Tyr72 and Tyr286 indicated that 12d was located near the peripheral anionic site of the enzyme. Additionally, in silicoADME investigations using QikProp demonstrated that 12d possesses optimal pharmacokinetic properties. In conclusion, this study identified a novel pyrrole-based AChE inhibitor 12d through a combination of computational and experimental findings.

Binding behavior of receptor binding domain of the SARS-CoV-2 virus and ivermectin

The COVID-19 pandemic, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), sparked an international debate on effective ways to prevent and treat the virus. Specifically, there were many varying opinions on the use of ivermectin (IVM) throughout the world, with minimal research to support either side. IVM is an FDA-approved antiparasitic drug that was discovered in the 1970s and was found to show antiviral activity. The objective of this study is to examine the binding behavior and rates of association and dissociation between SARS-CoV-2 receptor binding domain (RBD), IVM, and their combination using aminopropylsilane (APS) biosensors as surrogates for the hydrophobic interaction between the viral protein and human angiotensin-converting enzyme 2 (ACE2) receptors to determine the potential of IVM as a repurposed drug for SARS-CoV-2 prevention and treatment. The IVM, RBD, and combination binding kinetics were analyzed using biolayer interferometry (BLI) and validated with multiple in silico techniques including protein–ligand docking, molecular dynamics simulation, molecular mechanics-generalized Born surface area (MM-GBSA), and principal component analysis (PCA). Our results suggest that with increasing IVM concentrations the association rate with the hydrophobic biosensor increases with a simultaneous decrease in dissociation. Significant kinetic changes to RBD, when combined with IVM, were found only at a concentration a thousand times the approved dosage with minimal changes found over a 35-min time period. Our study suggests that IVM is not an effective preventative or treatment method at the currently approved dosage.

Evaluation of novel Anti-SARS-CoV-2 compounds by targeting nucleoprotein and envelope protein through homology modeling, docking simulations, ADMET, and molecular dynamic simulations with the MM/GBSA calculation

The current prominent virus that induces severe acute respiratory syndrome is SARS-CoV-2. The incidence of COVID-19 cases is increasing, necessitating the immediate development of effective treatments. Our objective was to employ an in-silico approach to evaluate the effectiveness of conventional compounds against COVID-19's nucleoprotein and envelope protein. A docking simulation was performed on 9 compounds as SARS-coronavirus inhibitors using AMDock software. Anti-covid-19 activities were further evaluated for the compounds. Based on docking results, the binding affinity of "N-(4-carbamoylphenyl)-8-cyclopropyl-7-(naphthalen-1-ylmethyl)-5-oxo-2,3-dihydro-5H-thiazolo[3,2-a]pyridine-3-carboxamide,” also called compound 36 in this research, was found to be −8.8 ​kcal/mol for the modeled envelope protein and −7.3 ​kcal/mol for the template envelope protein, while −10.1 ​kcal/mol for the modeled nucleocapsid proteins (NP) and −8.7 ​kcal/mol for the template nucleocapsid proteins (NP) of SARS-coronavirus, respectively. The ligand and control drug (ritonavir) with high docking scores were subjected to pharmacological screening, molecular dynamic simulations, and Molecular Mechanics-generalized Born Surface Area (MM/GBSA) calculations. Furthermore, the jobs of pharmacokinetics were assessed, and the outcomes acquired show that the proposed compound 36 includes great oral bioavailability and a capacity to diffuse through various organic boundaries. The protein-ligand complexes were subjected to dynamic simulation analyses with a re-enactment time of 10 ns, likewise, their free binding energy was inspected operating the MM/GBSA approach. The docking (MD simulation) results acquired emphasize the pivotal residues answerable for the protein-ligand interaction, giving an understanding of the method of association. The MD simulation analysis verifies the structural stability of the selected complexes during the MD trajectory, with minor changes detected. The MM/GBSA data show that compound 36 has the lowest free energy of −12.498 ​kcal/mol for EP and −57.5185 ​kcal/mol for NP proteins of SARS-coronavirus, confirming the molecular docking result. As a result, the identified chemical can be used to develop a new family of antiviral medications against SARS-coronavirus-2.

In silico and in vitro evaluation of the anti-virulence potential of patuletin, a natural methoxy flavone, against Pseudomonas aeruginosa.

This study aimed to investigate the potential of patuletin, a rare natural flavonoid, as a virulence and LasR inhibitor against Pseudomonas aeruginosa. Various computational studies were utilized to explore the binding of Patuletin and LasR at a molecular level. Molecular docking revealed that Patuletin strongly interacted with the active pocket of LasR, with a high binding affinity value of -20.96 kcal/mol. Further molecular dynamics simulations, molecular mechanics generalized Born surface area (MM/GBSA), protein-ligand interaction profile (PLIP), and essential dynamics analyses confirmed the stability of the patuletin-LasR complex, and no significant structural changes were observed in the LasR protein upon binding. Key amino acids involved in binding were identified, along with a free energy value of -26.9 kcal/mol. In vitro assays were performed to assess patuletin's effects on P. aeruginosa. At a sub-inhibitory concentration (1/4 MIC), patuletin significantly reduced biofilm formation by 48% and 42%, decreased pyocyanin production by 24% and 14%, and decreased proteolytic activities by 42% and 20% in P. aeruginosa isolate ATCC 27853 (PA27853) and P. aeruginosa clinical isolate (PA1), respectively. In summary, this study demonstrated that patuletin effectively inhibited LasR activity in silico and attenuated virulence factors in vitro, including biofilm formation, pyocyanin production, and proteolytic activity. These findings suggest that patuletin holds promise as a potential therapeutic agent in combination with antibiotics to combat antibiotic-tolerant P. aeruginosa infections.

Monastrol disrupts KIFC1-ATP dynamics: Towards newer anticancer mechanism

Mitotic kinesins are eukaryotic proteins that play a vital role in cellular mitosis. Their over-expression in malignant rather than normal cells and their specific cellular role made kinesin inhibitors a promising cancer therapeutics. The realization of in-vitro antitumor activity of the investigational molecule-monastrol-was traced to its allosteric inhibition of the ATPase activity of the motor domain of Human Kensin-5 (Eg5) and prompted extensive efforts towards kinesin inhibitors. Many reports pointed to the prospect of Kinesin-like protein (KIFC1) as a novel anticancer target. Nevertheless, no verified KIFC1-inhibitor crystallized structure has been reported so far, confirming the scarcity of molecular studies devoted to deciphering claimed potentiality. The significant structural resemblance across the kinesins superfamily sparked our interest in investigating monastrol as a KIFC1 inhibitor.Time-scale findings of molecular dynamics simulation and molecular mechanics/generalized-born surface area (MM/GBSA) methods revealed that monastrol binds with high affinity to an allosteric pocket, inducing notable dynamical alterations evidenced by increased structural stability, rigidity, compactness, and overall folding tendency of KIFC1. Further, these conformational events negatively impacted the affinity of ATP to its site disrupting the essential ATP-KIFC1 dynamics leading to consequential loss of functionality. We concluded that monastrol can act as non-competitive inhibitor of KIFC1 and these findings speculate new anticancer mechanism contributing to the design of highly selective and novel cancer therapeutics.

Determining Structural Changes for Ligand Recognition between Human and Rat Phosphorylated BACE1 in Silico and Its Phosphorylation by GSK3β at Thr252 by in Vitro Studies.

Alzheimer's disease (AD) is a neurodegenerative disease affecting older adults. AD pathogenesis involves the production of the highly neurotoxic amyloid-β peptide 1-42 (Aβ1-42) from β-site amyloid precursor protein cleaving enzyme 1 (BACE1). The phosphorylation of BACE1 at Thr252 increases its enzymatic activity. This study examined the phosphorylation of BACE1 from human and rat BACE1 in silico through phosphorylation predictors. Besides, we explored how phosphorylation at various sites affected the BACE1 structure and its affinity with amyloid precursor protein (APP) and six BACE1 inhibitors. Additionally, we evaluated the phosphorylation of Thr252-BACE1 by glycogen synthase kinase 3 β (GSK3β) in vitro. The phosphorylation predictors showed that Thr252, Ser59, Tyr76, Ser71, and Ser83 could be phosphorylated. Also, Ser127 in rat BACE1 can be phosphorylated, but human BACE1 has a Gly at this position. Molecular dynamics simulations showed that Ser127 plays an important role in the open and closed BACE1 conformational structures. Docking studies and the molecular mechanics generalized Born surface area (MMGBSA) approach showed that human BACE1 phosphorylated at Thr252 and rat BACE1 phosphorylated at Ser71 have the best binding and free energy with APP, forming hydrogen bonds with Asp672. Importantly, inhibitors have a higher affinity for the phosphorylated rat BACE1 than for its human counterpart, which could explain their failure during clinical trials. Finally, in vitro experiments showed that GSK3β could phosphorylate BACE1. In conclusion, BACE1 phosphorylation influences the BACE1 conformation and its recognition of ligands and substrates. Thus, these features should be carefully considered in the design of BACE1 inhibitors.

Novichok Nerve Agents as Inhibitors of Acetylcholinesterase-In Silico Study of Their Non-Covalent Binding Affinity.

In silico studies were performed to assess the binding affinity of selected organophosphorus compounds toward the acetylcholinesterase enzyme (AChE). Quantum mechanical calculations, molecular docking, and molecular dynamics (MD) with molecular mechanics Generalized-Born surface area (MM/GBSA) were applied to assess quantitatively differences between the binding energies of acetylcholine (ACh; the natural agonist of AChE) and neurotoxic, synthetic correlatives (so-called "Novichoks", and selected compounds from the G- and V-series). Several additional quantitative descriptors like root-mean-square fluctuation (RMSF) and the solvent accessible surface area (SASA) were briefly discussed to give-to the best of our knowledge-the first quantitative in silico description of AChE-Novichok non-covalent binding process and thus facilitate the search for an efficient and effective treatment for Novichok intoxication and in a broader sense-intoxication with other warfare nerve agents as well.

Assessing the performance of MM/PBSA and MM/GBSA methods. 10. Prediction reliability of binding affinities and binding poses for RNA-ligand complexes.

Ribonucleic acid (RNA)-ligand interactions play a pivotal role in a wide spectrum of biological processes, ranging from protein biosynthesis to cellular reproduction. This recognition has prompted the broader acceptance of RNA as a viable candidate for drug targets. Delving into the atomic-scale understanding of RNA-ligand interactions holds paramount importance in unraveling intricate molecular mechanisms and further contributing to RNA-based drug discovery. Computational approaches, particularly molecular docking, offer an efficient way of predicting the interactions between RNA and small molecules. However, the accuracy and reliability of these predictions heavily depend on the performance of scoring functions (SFs). In contrast to the majority of SFs used in RNA-ligand docking, the end-point binding free energy calculation methods, such as molecular mechanics/generalized Born surface area (MM/GBSA) and molecular mechanics/Poisson Boltzmann surface area (MM/PBSA), stand as theoretically more rigorous approaches. Yet, the evaluation of their effectiveness in predicting both binding affinities and binding poses within RNA-ligand systems remains unexplored. This study first reported the performance of MM/PBSA and MM/GBSA with diverse solvation models, interior dielectric constants (εin) and force fields in the context of binding affinity prediction for 29 RNA-ligand complexes. MM/GBSA is based on short (5 ns) molecular dynamics (MD) simulations in an explicit solvent with the YIL force field; the GBGBn2 model with higher interior dielectric constant (εin = 12, 16 or 20) yields the best correlation (Rp = -0.513), which outperforms the best correlation (Rp = -0.317, rDock) offered by various docking programs. Then, the efficacy of MM/GBSA in identifying the near-native binding poses from the decoys was assessed based on 56 RNA-ligand complexes. However, it is evident that MM/GBSA has limitations in accurately predicting binding poses for RNA-ligand systems, particularly compared with notably proficient docking programs like rDock and PLANTS. The best top-1 success rate achieved by MM/GBSA rescoring is 39.3%, which falls below the best results given by docking programs (50%, PLNATS). This study represents the first evaluation of MM/PBSA and MM/GBSA for RNA-ligand systems and is expected to provide valuable insights into their successful application to RNA targets.

Exploring the most promising anti - Depressant drug targeting Microtubule Affinity Receptor Kinase 4 involved in Alzheimer's Disease through molecular docking and molecular dynamics simulation.

Alzheimer's Disease (AD) is the prevailing type of neurodegenerative illness, characterised by the accumulation of amyloid beta plaques. The symptoms associated with AD are memory loss, emotional variability, and a decline in cognitive functioning. To date, the pharmaceuticals currently accessible in the marketplace are limited to symptom management. According to several research, antidepressants have demonstrated potential efficacy in the management of AD. In this particular investigation, a total of 24 anti-depressant medications were selected as ligands, while the Microtubule Affinity Receptor Kinase 4 (MARK4) protein was chosen as the focal point of our study. The selection of MARK4 was based on its known involvement in the advancement of AD and other types of malignancies, rendering it a highly prospective target for therapeutic interventions. The initial step involved doing ADMET analysis, which was subsequently followed by molecular docking of 24 drugs. This was succeeded by molecular dynamics simulation and molecular mechanics generalised Born surface area (MMGBSA) calculations. Upon conducting molecular docking experiments, it has been determined that the binding affinities observed fall within the range of -5.5 kcal/mol to -9.0 kcal/mol. In this study, we selected six anti-depressant compounds (CID ID - 4184, 2771, 4205, 5533, 4543, and 2160) based on their binding affinities, which were determined to be -9.0, -8.7, -8.4, -8.3, -8.2, and -8.2, respectively. Molecular dynamics simulations were conducted for all six drugs, with donepezil serving as the control drug. Various analyses were performed, including basic analysis and post-trajectory analysis such as free energy landscape (FEL), polarizable continuum model (PCM), and MMGBSA calculations. Based on the findings from molecular dynamics simulations and the MMGBSA analysis, it can be inferred that citalopram and mirtazapine exhibit considerable potential as anti-depressant agents. Consequently, these compounds warrant further investigation through in vitro and in vivo investigations in the context of treating AD.

A theoretical study on the activity and selectivity of IDO/TDO inhibitors.

Indoleamine 2,3-dioxygenase 1 (IDO) is a tryptophan (Trp) metabolic enzyme along the kynurenine (NFK) pathway. Under pathological conditions, IDO overexpressed by tumor cells causes depletion of tryptophan and the accumulation of metabolic products, which inhibit the local immune response and form immune escape. Therefore, the suppression of IDO activity is one of the strategies for tumor immunotherapy, and drug design for this target has been the focus of research for more than two decades. Apart from IDO, tryptophan dioxygenase (TDO) of the same family can also catalyze the same biochemical reaction in the human body, but it has different tissue distribution and substrate selectivity from IDO. Based on the principle of drug design with high potency and low cross-reactivity to specific targets, in this subject, the activity and selectivity of IDO and TDO toward small molecular inhibitors were studied from the perspective of thermodynamics and kinetics. The aim was to elucidate the structural requirements for achieving favorable biological activity and selectivity of IDO and TDO inhibitors. Specifically, the interactions of inhibitors from eight families with IDO and TDO were initially investigated through molecular docking and molecular dynamics simulations, and the thermodynamic data for binding of inhibitors were predicted by the molecular mechanics/generalized Born surface area (MM/GBSA) method. Secondly, we explored the free energy landscape of JKloops, the kinetic control element of IDO/TDO, using temperature replica exchange molecular dynamics (T-REMD) simulations and elucidated the connection between the rules of IDO/TDO conformational changes and the inhibitor selectivity mechanism. Furthermore, the binding and dissociation processes of the C1 inhibitor (NLG919) were simulated by the adaptive steering molecular dynamics (ASMD) method, which not only addressed the possible stable, metastable, and transition states for C1 inhibitor-IDO/TDO interactions, but also accurately predicted kinetic data for C1 inhibitor binding and dissociation. In conclusion, we have constructed a complete process from enzyme (IDO/TDO) conformational activation to inhibitor binding/dissociation and used the thermodynamic and kinetic data of each link as clues to verify the control mechanism of IDO/TDO on inhibitor selectivity. This is of great significance for us to understand the design principles of tumor immunotherapy drugs and to avoid drug resistance of immunotherapy drugs.

Insights into the chirality-dependent recognition of Danshensu Bingpian Zhi stereoisomers with PPARγ.

Peroxisome proliferator-activated receptor γ (PPARγ), a nuclear receptor involved in metabolic processes, inflammation, and energy balance, represents a promising therapeutic target for cardiovascular diseases. Danshensu Bingpian Zhi (DBZ), a chiral compound derived from traditional Chinese medicine, exhibits potential as a PPARγ agonist. Using an ensemble-based docking approach, molecular dynamics (MD) simulations, and the molecular mechanics generalized born surface area (MM/GBSA) methods, we explored the binding modes and energetics of DBZ stereoisomers with the PPARγ ligand-binding domain (LBD). The results indicated that the right-handed stereoisomer (DBZR) binds like a full agonist, while the left-handed stereoisomer (DBZS) binds as a partial agonist with stronger binding energies (ΔGbind), indicating a robust interaction with PPARγ. Both the stereoisomers stabilize the β-sheet region of PPARγ-LBD, potentially protecting Ser245 from phosphorylation by Cdk5, a process implicated in atherosclerosis. Principal component analysis (PCA) and dynamic cross-correlation matrices (DCCM) revealed the complex structural dynamics within the Ω loop, β-sheet, and AF-2 region of PPARγ-LBD upon ligand binding, which may contribute to the unique binding mode and efficacy of DBZS. These findings provide insights into the molecular recognition of PPARγ-LBD by DBZ stereoisomers and their impact on the conformational dynamics of PPARγ, highlighting the therapeutic potential of DBZ and the significance of chirality in drug design.

Accelerated molecular dynamics study of the interaction mechanism between small molecule inhibitors and phosphoglycerate mutase 1.

In 2020, cancer-related deaths reached 9.96 million globally, of which China accounted for 3 million, ranking first in the world. Phosphoglycerate mutase 1 (PGAM1) is a key metabolic enzyme in glycolysis, catalysing the conversion of 3-phosphoglycerate to 2-phosphoglycerate. Based on the excellent anticancer activity of PGMI-004A and HKB99, new small molecules with an anthraquinone core were synthesised to inhibit tumour growth. Developing small molecules with an anthraquinone core targeting PGAM1 may be an effective strategy for treating cancer. In this study, accelerated molecular dynamics (aMD) simulation, dynamic cross-correlation map (DCCM) calculation, principal component analysis (PCA) and free energy landscape (FEL) analysis were used to analyse conformational changes of PGAM1 caused by binding of inhibitors 8KX, 9HU and HKB. DCCM calculations and PCA showed that inhibitor binding significantly affected the kinetic behaviour of PGAM1 and conformational rearrangement of PGAM1. The binding ability and mechanism of 8KX, 9HU and HKB to PGAM1 were studied using the molecular mechanics generalised Born surface area (MM-GBSA) method. The results showed that compared with 8KX, the binding ability of 9HU and HKB to PGAM1 was enhanced by sulphonamide reversal and aminocarboxyl trifluoromethyl substitution. There were several hydrophobic interactions between inhibitors and PGAM1, providing significant contributions for inhibitor binding. Calculation of residue-based free energy decomposition revealed that F22, R90, Y92, L95, V112, W115, R116, V121, P123, P124, R191 and M206 were key residues of the PGAM1-inhibitor interaction and could be used as effective targets for designing drugs that inhibit the activity of PGAM1.

The potentiality of isoflavones from Derris robusta (DC.) Benth. against α-glucosidase.

The four recently described isoflavones, derrisrobustones A-D, and eight known isoflavones extracted from Derris robusta (DC.) Benth. (Leguminosae) were identified as potential isoflavones for α-glucosidase inhibitory activity. However, the inhibitory ability of each isoflavone was quite different and the structural and dynamical properties of the isoflavones in the binding with α-glucosidase have not been studied. In this study, the molecular dynamics (MD) simulation and the molecular mechanics generalized Born surface area (MM/GBSA) methods were applied to gain insight into isoflavones bound with the α-glucosidase active site. The results found that the isoflavones with a hydroxy substituent could be stabilized in the α-glucosidase active site with low values of the binding free energy, especially derrisrobustone D and derrubone. The isoflavones, except for derrubone, were bound with α-glucosidase by means of electrostatic force effectively.

Investigating the Potency of Erythrina‒Derived Flavonoids as Cholinesterase Inhibitors and Free Radical Scavengers Through in silico Approach: Implications for Alzheimer's Disease Therapy.

This study aimed to evaluate the potency of 471 flavonoids from the genus Erythrina as potential acetylcholinesterase (AChE) inhibitors and free radical scavengers through computational studies to develop Alzheimer's disease (AD) therapies from natural products. A total of 471 flavonoids from the genus Erythrina were subjected to molecular docking against AChE, followed by toxicity screening. The potential AChE inhibitors with the least toxic profile were subjected to further investigation through molecular dynamics (MD) simulations, density functional theory (DFT) study, and in silico pharmacokinetic predictions. A combination of molecular docking and in silico toxicity screening led to the identification of 2(S)‒5,7‒dihydroxy‒5'‒methoxy‒[2'',2''‒(3''‒hydroxy)‒dimethylpyrano]‒(5'',6'':3',4') flavanone (89) and Abyssinoflavanone IV (83) as potential AChE inhibitors. These compounds had stable binding to AchE and exhibited lower Root Mean Square Deviation (RMSD) values compared to the apo state of AChE. In addition, Molecular Mechanics Generalized Born Surface Area (MMGBSA) analysis revealed that the binding energies of 89 and 83 were significantly lower compared to acetylcholine, the natural substrate of AChE. Based on DFT study, these compounds exhibited a higher energy in the highest occupied molecular orbital (EHOMO) and lower electron affinity (EA) than Quercetin. This indicated that 89 and 83 could be potential radical scavengers through their electron-donating activity. Although this study primarily relied on computational methods, the results showed the dual functionality of compounds 89 and 83 as both potential AChE inhibitors and free radical scavengers. Further investigation in wet laboratory experiments is required to validate their therapeutic potential for AD.

Unearthing phytochemicals as natural inhibitors for pantothenate synthetase in Mycobacterium tuberculosis: A computational approach.

Pantothenate synthetase protein plays a pivotal role in the biosynthesis of coenzyme A (CoA), which is a crucial molecule involved in a number of cellular processes including the metabolism of fatty acid, energy production, and the synthesis of various biomolecules, which is necessary for the survival of Mycobacterium tuberculosis (Mtb). Therefore, inhibiting this protein could disrupt CoA synthesis, leading to the impairment of vital metabolic processes within the bacterium, ultimately inhibiting its growth and survival. This study employed molecular docking, structure-based virtual screening, and molecular dynamics (MD) simulation to identify promising phytochemical compounds targeting pantothenate synthetase for tuberculosis (TB) treatment. Among 239 compounds, the top three (rutin, sesamin, and catechin gallate) were selected, with binding energy values ranging from -11 to -10.3kcal/mol, and the selected complexes showed RMSD (<3Å) for 100 ns MD simulation time. Furthermore, molecular mechanics generalized Born surface area (MM/GBSA) binding free energy calculations affirmed the stability of these three selected phytochemicals with binding energy ranges from -82.24 ± 9.35 to -66.83 ± 4.5kcal/mol. Hence, these identified natural plant-derived compounds as potential inhibitors of pantothenate synthetase could be used to inhibit TB infection in humans.