Molecular Dynamics Free Energy Research Articles

Crystal structures of coronaviral main proteases in complex with the non-covalent inhibitor X77

Main proteases (Mpros) are a class of conserved cysteine hydrolases among coronaviruses and play a crucial role in viral replication. Therefore, Mpros are ideal targets for the development of pan-coronavirus drugs. X77, previously developed against SARS-CoV Mpro, was repurposed as a non-covalent tight binder inhibitor against SARS-CoV-2 Mpro during COVID-19 pandemic. Many novel inhibitors with favorable efficacy have been discovered using X77 as a reference, suggesting that the structure of X77 could be a valuable scaffold for drug design. However, the broad-spectrum performance of X77 and underlying mechanism remain less understood. Here, we reported the crystal structures of Mpros from SARS-CoV-2, SARS-CoV, and MERS-CoV, and several Mpro mutants from SARS-CoV-2 variants bound to X77. A detailed analysis of these structures reveals key structural determinants essential for interaction and elucidates the binding modes of X77 to different coronaviral Mpros. The potencies of X77 against these investigated Mpros were further evaluated through molecular dynamic simulation and binding free energy calculation. These data provide molecular insights into broad-spectrum inhibition against coronaviral Mpros by X77 and the similarities and differences of X77 when bound to various Mpros, which will promote X77-based design of novel antivirals with broad-spectrum efficacy against different coronaviruses and SARS-CoV-2 variants.

Computational investigation of nanoemulsion-released tocopherol for treating cardiovascular diseases: Insights from molecular dynamic simulation and deep learning prediction

This research investigated the potential of tocopherol released from nanoemulsion as a therapeutic agent for cardiovascular diseases (CVD) through comprehensive in-silico site-specific molecular docking analysis and AI- driven Molecular Dynamic Simulation, between tocopherol and seven CVD-related proteins (Angiotensin Converting Enzyme (ACE) (1O8a), Human Angiotensin Receptor (4YAY), P38 Mitogen-Activated Protein Kinase (MAPK) (4DLI), HMG-CoA Reductase (1HW9), β-1-Adrenergic Receptor (2YCW), Human C-reactive Protein (1BO9), and Cyclooxygenase-2 (1CX2)). Our findings revealed significant interactions. α-Tocopherol exhibited strong binding affinities, particularly with the 2YCW protein, evidenced by a Vina score of -11.9 ± 0.86 kcal/mol and an inhibition constant of 1.4 ± 0.12 µM. Extensive hydrophobic interactions and significant hydrogen bonds were identified as critical for the stability of these complexes. Furthermore, Molecular Dynamics (MD) simulations, RMSD, Rg, RMSF, and Free Energy Landscape (FEL) profiles, and a deep convolutional neural network (CNN) used for predicting interaction sites from protein sequences, illustrated that the protein-ligand complexes attained dynamic configuration stability with RMSD values converging around 0.3 nm, underscoring the robustness of tocopherol-protein interactions. The pathway analysis highlights the possibility of tocopherol's involvement in interacting with several key biological pathways, Dilated Cardiomyopathy, Arrhythmogenic Right Ventricular Cardiomyopathy, Hypertrophic Cardiomyopathy, and Apoptosis.This study revealed significant insights into the interactions of Nanoemulsion-released tocopherol with key proteins involved in the development of cardiovascular diseases, aimed at improving cardiovascular health. However, the findings of the study require meticulous experimental study and clinical trials to provide accurate results and validate tocopherol's therapeutic efficacy.

H and CO Co-Induced Roughening of Cu Surface in CO2 Electroreduction Conditions.

The dynamic restructuring of Cu has been observed under electrochemical conditions, and it has been hypothesized to underlie the unique reactivity of Cu toward CO2 electroreduction. Roughening is one of the key surface phenomena for Cu activation, whereby numerous atomic vacancies and adatoms form. However, the atomic structure of such surface motifs in the presence of relevant adsorbates has remained elusive. Here, we explore the chemical space of Cu surface restructuring under coverage of CO and H in realistic electroreduction conditions, by combining grand canonical DFT and global optimization techniques, from which we construct a potential-dependent grand canonical ensemble representation. The regime of intermediate and mixed CO and H coverage─where structures exhibit some elevated surface Cu─is thermodynamically unfavorable yet kinetically inevitable. Therefore, we develop a quasi-kinetic Monte Carlo simulation to track the system's evolution during a simulated cathodic scan. We reveal the evolution path of the system across coverage space and identify the accessible metastable structures formed along the way. Chemical bonding analysis is performed on the metastable structures with elevated Cu*CO species to understand their formation mechanism. By molecular dynamics simulations and free energy calculations, the surface chemistry of the Cu*CO species is explored, and we identify plausible mechanisms via which the Cu*CO species may diffuse or dimerize. This work provides rich atomistic insights into the phenomenon of surface roughening and the structure of involved species. It also features generalizable methods to explore the chemical space of restructuring surfaces with mixed adsorbates and their nonequilibrium evolution.

Direct Interaction of Minocycline to p47phox Contributes to its Attenuation of TNF-α-Mediated Neuronal PC12 Cell Death: Experimental and Simulation Validation.

Kaempferol illustrates an example of attempting to discover new treatments for neurodegeneration by investigating the potential efficacy of natural products. Despite the identification of several molecular targets for this bio-active compound, the precise underlying pathways are not well elucidated yet. Recently, it has been shown through pulldown assay that kaemferol directly interacts with p47phox, the organizer subunit of NADPH oxidase-2 (NOX2) complex. Hence, in this study, we used homology modelling, computational docking, mutation analysis, molecular dynamics simulations and free energy calculations to determine how kaempferol interacts with p47phox. Firstly, 3D structure of p47phox was generated using x-ray structures of its domains. Then, it was docked with kaempferol, and finally 100-ns molecular dynamics (MD) simulations were performed and the global properties like root-mean square deviation (RMSD) and root-mean square fluctuations (RMSF) were calculated. Literature survey and computational analysis of key interacting amino acid residues of p47phox provided insights into a possible binding site for kaempferol, approximately around Trp193 and Cys196 located within the N-terminal SH3 domain of p47phox. Moreover, free energy calculations indicated that in silico substitution of Trp193 and Cys196 with arginine and alanine, respectively, results in less favorable interaction corroborating their importance in binding with kaempferol. Taken together, these findings suggest that kaempferol directly attaches to N-SH3 domain p47 phox, with a subsequent diminution of p47phox protein-protein interaction and possibly attenuation of NOX2 complex assembly, which reduces reactive oxygen species (ROS) generation. These observations will be beneficial for researchers exploring neuroprotection and for the development of p47phox inhibitors.

Phloretin prolongs lifespan of Caenorhabditis elegans via inhibition of NDUFS1 and NDUFS6 at mitochondrial complex Ⅰ

Phloretin has been widely perceived as an antioxidant. However, the bioavailability of phloretin in vivo is generally far too low to elicit a direct antioxidant effect by scavenging reactive oxygen species (ROS). Here we showed that administration of phloretin of apple polyphenols extended lifespan of Caenorhabditis elegans and promoted fitness. Specially phloretin enhanced the survival rates of nematodes under oxidants in an inverted U-shaped dose-response manner. The lifespan-extending effects of phloretin were mediated by ROS via mitochondrial complex I inhibition. The increase of ROS stimulated p38 MAPK/PMK-1 as well as transcription factors of NRF2/SKN-1 and FOXO/DAF-16. Consistent with the involvement of NRF2/SKN-1 and FOXO/DAF-16 in lifespan-extending effects, activities of superoxide dismutase (SOD) and catalase (CAT) were enhanced by phloretin. The exogenous application of antioxidants butylated hydroxyanisole and N-acetylcysteine abolished the increase of ROS, the enhancement of SOD and CAT activities, and the lifespan extending effects of phloretin. Meanwhile, with the inhibition of mitochondrial complex I, ATP was instantly decreased. Both energy sensors of AMPK/AAK-2 and SIRT1/SIR-2.1 were involved in the lifespan extension by phloretin. Transcriptomic, real-time qPCR and molecular docking analyses demonstrated that the binding of phloretin at complex I located at NDUFS1/NUO-5, NDUFS2/GAS-1, and NDUFS6/NDUF-6. The molecular dynamic simulation and binding free energy calculations showed that phloretin had high binding affinities towards NDUFS1 (−7.21 kcal/mol) and NDUFS6 (−7.02 kcal/mol). Collectively, our findings suggested phloretin had effects of life expectancy enhancement and fitness promotion via redox regulations in vivo. NDUFS1/NUO-5 and NDUFS6/NDUF-6 might be new targets in the lifespan and wellness regulations.

Minimally modified human blood coagulation factor X to bypass direct factor Xa inhibitors

BackgroundDirect oral factor (F)Xa inhibitors are widely used as alternatives to conventional vitamin K antagonists in managing venous thromboembolism and nonvalvular atrial fibrillation. Unfortunately, bleeding-related adverse events remain a major concern in clinical practice. In case of bleeding or emergency surgery, rapid-onset reversal agents may be required to counteract the anticoagulant activity. ObjectivesThe ability of FXa variants to bypass the direct oral FXa inhibitors was assessed. MethodsHuman FXa variants were generated through substitution of phenylalanine 174 (F174) for either alanine, isoleucine, or serine. FXa variants were stably expressed in HEK293 cells and purified to homogeneity using ion-exchange chromatography. ResultsF174-substituted human FX variants demonstrated efficacy in restoring thrombin generation in plasma containing direct FXa inhibitors (apixaban, rivaroxaban, edoxaban). Their ability to bypass the anticoagulant effects stems from a significantly reduced sensitivity for the direct FXa inhibitors due to a decrease in binding affinity determined using molecular dynamics simulations and free energy computation. Furthermore, F174 modification resulted in a partial loss of inhibition by tissue factor pathway inhibitor, enhancing the procoagulant effect of F174-substituted FX. Consequently, the F174A- and F174S-substituted FX variants effectively counteracted the effects of 2 widely used anticoagulants, apixaban and rivaroxaban, in plasma of atrial fibrillation and venous thromboembolism patients. ConclusionThese human FX variants have the potential to serve as a rescue reversal strategy to overcome the effect of direct FXa inhibitors in case of life-threatening bleeding events or emergency surgical interventions.

Unveiling Therapeutic Avenues for Crohn's Disease Management: Exploring Inhibitors for Adherent-Invasive Escherichia coli Propanediol Dehydratase.

Inflammatory Bowel Disease (IBD) encompasses a group of chronic disorders distinguished by inflammation of the gastrointestinal tract. Among these, Crohn's Disease (CD) stands out as a complex and impactful condition due to challenges for both diagnosis and management, making it a cynosure of research. In CD, there is the predominance of proinflammatory bacteria, including the Adherentinvasive Escherichia coli (AIEC) with virulence-associated metabolic enzyme Propanediol Dehydratase (pduC), which has been identified as a therapeutic target for the management of CD. Herein, molecular modeling techniques, including molecular docking, Molecular Mechanics with Generalized Born and Surface Area (MMGBSA), drug-likeness, and pharmacokinetics profiling, were utilized to probe the potentials of eighty antibacterial compounds to serve as inhibitors of pduC. The results of this study led to the identification of five compounds with promising potentials; the results of the molecular docking simulation revealed the compounds as possessing better binding affinities for the target compared to the standard drug (sulfasalazine), while Lipinski's rule of five-based assessment of their drug-likeness properties revealed them as potential oral drugs. MMGBSA free energy calculation and Molecular Dynamics (MD) simulation of the complexes formed a sequel to molecular docking, revealing the compounds as stable binders in the active site of the protein. Ultimately, the results of this study have revealed five compounds to possess the potential to serve as inhibitors of pduC of AIEC. However, experimental studies are still needed to validate the findings of this study.

Prediction of deleterious non-synonymous SNPs of TMPRSS2 protein combined with Molecular Dynamics Simulations and free energy analysis to identify the potential peptide substrates against SARS-CoV-2.

Globally the SARS-CoV-2 viral infection demands for the new drugs, the TMPRSS2 target plays a vital role in facilitating the virus entry. The aim of the present study is to identify the potential peptide substrate from the Anti-viral database against TMPRSS2 of SARS-CoV-2. The compound screening and variation analysis were performed using molecular docking analysis and online tools such as PROVEAN and SNAP2 server, respectively. The re-docked crystal structure peptide substrate exhibits -128.151 kcal/mol whereas the RRKK peptide substrate shows -134.158 kcal/mol. Further, the selected compounds were proceeded with Molecular Dynamics Simulation, it explores the stability of the complex by revealing the hotspot residues (His296 and Ser441) were active for nucleophilic attack against TMPRSS2. The average Binding Free Energy values computed through MM/GBSA for RRKK, Camostat, and Crystal Structure were shown -69.9278 kcal/mol, -64.5983 kcal/mol, and -63.9755 kcal/mol, respectively against TMPRSS2. The 'rate of acylation' emerges as an indicator for RRKK's efficacy, it maintains the distance of 3.2 Å with Ser441 resembles, whilst its -NH backbone stabilizes at 2.5 Å 'Michaelis Complex' which leads to prevent the entry of SARS-CoV-2 to human cells. The sequence variation analysis explores that the V160 and G6 substitutions are essential to emphasize the uncover possibilities for the ongoing drug discovery research. Therefore, the identified peptide substrate found to be potent against SARS-CoV-2 and these results will be valuable for ongoing drug discovery research.Communicated by Ramaswamy H. Sarma.

Exploring the Super-Relaxed State of Human Cardiac β-Myosin by Molecular Dynamics Simulations.

Human β-cardiac myosin plays a critical role in generating the mechanical forces necessary for cardiac muscle contraction. This process relies on a delicate dynamic equilibrium between the disordered relaxed state (DRX) and the super-relaxed state (SRX) of myosin. Disruptions in this equilibrium due to mutations can lead to heart diseases. However, the structural characteristics of SRX and the molecular mechanisms underlying pathogenic mutations have remained elusive. To bridge this gap, we conducted molecular dynamics simulations and free energy calculations to explore the conformational changes in myosin. Our findings indicate that the size of the phosphate-binding pocket can serve as a valuable metric for characterizing the transition from the DRX to SRX state. Importantly, we established a global dynamic coupling network within the myosin motor head at the residue level, elucidating how the pathogenic mutation E483K impacts the equilibrium between SRX and DRX through allosteric effects. Our work illuminates molecular details of SRX and offers valuable insights into disease treatment through the regulation of SRX.

Iso-mukaadial acetate and ursolic acid acetate bind to Plasmodium Falciparum heat shock protein 70: towards targeting parasite protein folding pathway.

Plasmodium falciparum is the most lethal malaria parasite. P. falciparum Hsp70 (PfHsp70) is an essential molecular chaperone (facilitates protein folding) and is deemed a prospective antimalarial drug target. The present study investigates the binding capabilities of select plant derivatives, iso-mukaadial acetate (IMA) and ursolic acid acetate (UAA), against P. falciparum using an in silico docking approach. The interaction between the ligands and PfHsp70 was evaluated using plasmon resonance (SPR) analysis. Molecular docking, binding free energy analysis and molecular dynamics simulations were conducted towards understanding the mechanisms by which the compounds bind to PfHsp70. The molecular docking results revealed ligand flexibilities, conformations and positions of key amino acid residues and protein-ligand interactions as crucial factors accounting for selective inhibition of Hsp70. The simulation results also suggest protein-ligand van der Waals forces as the driving force guiding the interaction of these compounds with PfHsp70. Of the two compounds, UAA and IMA bound to PfHsp70 within the micromolar range based on surface plasmon resonance (SPR) based binding assay. Our findings pave way for future rational design of new selective compounds targeting PfHsp70.

Lomitapide as a Potential Estrogen Receptor Inhibitor: A Computational Drug Repurposing Study

Objective: Estrogen receptor (ER) inhibitors have significant therapeutic potential for hormone-dependent cancers and related disorders. Tamoxifen, a well-known selective estrogen receptor modulator, has been widely used as adjuvant therapy for estrogen receptor-positive breast cancer. However, tamoxifen may exhibit a tendency to develop resistance with prolonged usage and particularly elevate the risk of uterine cancer. Therefore, there is a need for the discovery and development of new ER modulators or inhibitors. In this study, we identified potential estrogen receptor inhibitors through computational drug repositioning.
 
 Methods: A set of 2048 compounds, encompassing FDA-approved drugs and active metabolites, were subjected to molecular docking, molecular dynamics simulations, and free energy calculations to evaluate their interaction with estrogen receptor α (ERα).
 
 Results: Among the compounds evaluated, conivaptan, atogepant, and lomitapide exhibited the highest affinities for ERα. Lomitapide displayed a superior docking score (-12 kcal/mol) compared to the established ER inhibitor, tamoxifen (-10 kcal/mol). Further investigation using molecular dynamics simulations and free energy calculations disclosed lomitapide's heightened binding affinity of -380.727 kJ/mol, surpassing tamoxifen's binding affinity of -352.029 kJ/mol.
 
 Conclusion: This comprehensive computational exploration underscores lomitapide's potential as a compelling candidate with an envisaged stronger estrogen receptor affinity than the acknowledged standard, tamoxifen. To validate lomitapide's promise as a novel ER inhibitor, essential in vitro and in vivo studies are suggested. These investigations will provide essential insights into lomitapide's reposition in addressing the challenges tied to hormone-dependent cancers and associated maladies.

In silico validation of hyaluronic acid – drug conjugates based targeted drug delivery for the treatment of COVID-19

The impact of COVID-19 urges scientists to develop targeted drug delivery to manage Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) viral infections with a fast recovery rate. The aim of the study is to develop Hyaluronic Acid (HA) drug conjugates of viral drugs to target two important enzymes (Mpro and PLpro) of SARS-CoV-2. Three antiviral drugs, namely Dexamethasone (DEX), Favipiravir (FAV), and Remdesivir (REM), were chosen for HA conjugation due to their reactive functional groups. Free forms of drugs (DEX, FAV, REM) and HA drug conjugates (HA-DEX, HA-FAV, HA-REM, HA-RHA, HA-RHE) were validated against Mpro (PDB ID 6LU7) and PLpro (PDB 7LLZ), which play an essential role in the replication and reproduction of the SARS-CoV-2 virus. The results of the present study revealed that HA-drug conjugates possess higher binding affinity and the best docking score towards the Mpro and PLpro target proteins of SARS-CoV-2 than their free forms of drugs. ADMET screening resulted that HA-drug conjugates exhibited better pharmacokinetic profiles than their pure forms of drugs. Further, molecular dynamic simulation studies, essential dynamics and free energy landscape analyses show that HA antiviral drug conjugates possess good trajectories and energy status, with the PLpro target protein (PDB 7LLZ) of SARS-CoV-2 through long-distance (500 ns) simulation screening. The research work recorded the best drug candidate for Cell-Targeted Drug Delivery (CTDD) for SARS-CoV-2-infected cells through hyaluronic acid conjugates of antiviral drugs. Communicated by Ramaswamy H. Sarma

Computational analysis followed by invitro studies to explore cytokines (TNF-α, IL-6 and IL-1β) inhibition potential of the new natural molecule polonilignan.

Targeting pro-inflammatory cytokines and their production is found to be of therapeutic benefit for the regulation of inflammation in various chronic autoimmune diseases. Our continued efforts to discover small molecular-weight pro-inflammatory cytokine inhibitors resulted in identifying a novel natural lignan molecule named polonilignan, isolated from the culture broth extract of an endophytic fungus Penicillium polonicum. An in silico study (molecular docking, ADME predictions, binding free energy calculation and molecular dynamics simulation) of the polonilignan over the pro-inflammatory cytokines proteins TNF-α, IL-6 and IL-1β was performed using Schrodinger LLC software to understand the binding interactions, drug-like properties, and stability of the ligand-protein complex. Further, in-vitro testing of inhibition of TNF-α, IL-6 and IL-1β by polonilignan was carried out using ELISA and RT-PCR on LPS-induced RAW 264.7 cell lines along with the testing of nitrite production effect (Griess assay) and cytotoxicity (MTT) analysis. Under the computational study, polonilignan revealed good docking scores, binding interactions, and stability under MDS and desirable in silico ADME results over the proteins TNF-α, IL-1β and IL-6. Poloniligan showed significant inhibition of IL-1β, IL-6 and TNF-α with IC50 values of 2.01 μM, 6.59 μM and 42.10 μM, respectively. Also, it reduced the translocation of the NF-κB subunit p65 to the nucleus (confocal microscopy). The mRNA expression levels of pro-inflammatory markers IL-1β, TNF-α and IL-6 levels were lowered significantly (p < .001) by the compound, and the diminution was higher with IL-1β. Further, the lignan was non-cytotoxic and effective in attenuating nitrite release (IC50 48.56 μM). Thus, polonilignan has been identified as a new pan-cytokine and NO inhibitor, it is recommended to optimise a method for the synthesis of this small molecular weight lignan and explore its pharmacokinetic characteristics, toxicity and therapeutic effect under various chronic inflammatory disease models.

Cyclization Scaffolding for Improved Vaccine Immunogen Stability: Application to Tau Protein in Alzheimer's Disease.

Effective scaffolding of immunogens is crucial for generating conformationally selective antibodies through active immunization, particularly in the treatment of protein misfolding diseases such as Alzheimer's and Parkinson's disease. Previous computational work has revealed that a disorder-prone region of the tau protein, when in a stacked form, is predicted to structurally resemble a small, soluble protofibril, having conformational properties similar to those of experimental in vitro tau oligomers. Such an oligomeric structural mimic has the potential to serve as a vaccine immunogen design for Alzheimer's disease. In this study, we developed a cyclization scaffolding method in Rosetta, in which multiple cyclic peptides are stacked into a protofibril. Cyclization results in significant stabilization of protofibril-like structures by constraining the conformational space. Applying this method to the disorder-prone region of the tau fibril, we evaluated the metastability of the cyclized tau immunogen using molecular dynamics simulations, and we identified sequences of two cyclic constructs having high metastability in the protofibril. We then assessed their thermodynamic stability by computing the free energy required to separate a distal chain from the rest of the stacked structure. Our computational results, based on molecular dynamics simulations and free energy calculations, demonstrate that two cyclized constructs, cyclo-(VKSEKLDFKDRVQSKIFyN) and cyclo-(VKSEKLDFKDRVQSKIYvG) (lowercase letters indicate d-form amino acids), possess significantly increased thermodynamic stability in the protofibril over an uncyclized linear construct VKSEKLDFKDRVQSKI. The cyclization scaffolding approach proposed here holds promise as a means to effectively design immunogens for protein misfolding diseases, particularly those involving liposome-conjugated peptide constructs.

Molecular Mechanisms Underlying Medium-Chain Free Fatty Acid-Regulated Activity of the Phospholipase PlaF from Pseudomonas aeruginosa.

PlaF is a membrane-bound phospholipase A1 from Pseudomonas aeruginosa that is involved in remodeling membrane glycerophospholipids (GPLs) and modulating virulence-associated signaling and metabolic pathways. Previously, we identified the role of medium-chain free fatty acids (FFAs) in inhibiting PlaF activity and promoting homodimerization, yet the underlying molecular mechanism remained elusive. Here, we used unbiased and biased molecular dynamics simulations and free energy computations to assess how PlaF interacts with FFAs localized in the water milieu surrounding the bilayer or within the bilayer and how these interactions regulate PlaF activity. Medium-chain FFAs localized in the upper bilayer leaflet can stabilize inactive dimeric PlaF, likely through interactions with charged surface residues, as has been experimentally validated. Potential of mean force (PMF) computations indicate that membrane-bound FFAs may facilitate the activation of monomeric PlaF by lowering the activation barrier for changing into a tilted, active configuration. We estimated that the coupled equilibria of PlaF monomerization-dimerization and tilting at the physiological concentration of PlaF lead to the majority of PlaF forming inactive dimers when in a cell membrane loaded with decanoic acid (C10). This is in agreement with a suggested in vivo product feedback loop and gas chromatography-mass spectrometry profiling results, indicating that PlaF catalyzes the release of C10 from P. aeruginosa membranes. Additionally, we found that C10 in the water milieu can access the catalytic site of active monomeric PlaF, contributing to the competitive component of C10-mediated PlaF inhibition. Our study provides mechanistic insights into how medium-chain FFAs may regulate the activity of PlaF, a potential bacterial drug target.

Harnessing Brazilian biodiversity database: identification of flavonoids as potential inhibitors of SARS-CoV-2 main protease using computational approaches and all-atom molecular dynamics simulation.

SARS-CoV-2 (Severe Acute Respiratory Syndrome Coronavirus 2) is the etiological agent responsible for the global outbreak of COVID-19 (Coronavirus Disease 2019). The main protease of SARS-CoV-2, Mpro, is a key enzyme that plays a vital role in mediating viral replication and transcription. In this study, a comprehensive computational approach was employed to investigate the binding affinity, selectivity, and stability of natural product candidates as potential new antivirals acting on the viral polyprotein processing mediated by SARS-CoV-2 Mpro. A library of 288 flavonoids extracted from Brazilian biodiversity was screened to select potential Mpro inhibitors. An initial filter based on Lipinski's rule of five was applied, and 204 compounds that did not violate any of the Lipinski rules were selected. The compounds were then docked into the active site of Mpro using the GOLD program, and the poses were subsequently re-scored using MM-GBSA (Molecular Mechanics Generalized Born Surface Area) binding free energy calculations performed by AmberTools23. The top five flavonoids with the best MM-GBSA binding free energy values were selected for analysis of their interactions with the active site residues of the protein. Next, we conducted a toxicity and drug-likeness analysis, and non-toxic compounds were subjected to molecular dynamics simulation and free energy calculation using the MM-PBSA (Molecular Mechanics Poisson-Boltzmann Surface Area) method. It was observed that the five selected flavonoids had lower MM-GBSA binding free energy with Mpro than the co-crystal ligand. Furthermore, these compounds also formed hydrogen bonds with two important residues, Cys145 and Glu166, in the active site of Mpro. Two compounds that passed the drug-likeness filter showed stable conformations during the molecular dynamics simulations. Among these, NuBBE_867 exhibited the best MM-PBSA binding free energy value compared to the crystallographic inhibitor. Therefore, this study suggests that NuBBE_867 could be a potential inhibitor against the main protease of SARS-CoV-2 and may be further examined to confirm our results.

An Efficient Approach to the Accurate Prediction of Mutational Effects in Antigen Binding to the MHC1.

The major histocompatibility complex (MHC) can recognize and bind to external peptides to generate effective immune responses by presenting the peptides to T cells. Therefore, understanding the binding modes of peptide-MHC complexes (pMHC) and predicting the binding affinity of pMHCs play a crucial role in the rational design of peptide vaccines. In this study, we employed molecular dynamics (MD) simulations and free energy calculations with an Alanine Scanning with Generalized Born and Interaction Entropy (ASGBIE) method to investigate the protein-peptide interaction between HLA-A*02:01 and the G9209 peptide derived from the melanoma antigen gp100. The energy contribution of individual residue was calculated using alanine scanning, and hotspots on both the MHC and the peptides were identified. Our study shows that the pMHC binding is dominated by the van der Waals interactions. Furthermore, we optimized the ASGBIE method, achieving a Pearson correlation coefficient of 0.91 between predicted and experimental binding affinity for mutated antigens. This represents a significant improvement over the conventional MM/GBSA method, which yields a Pearson correlation coefficient of 0.22. The computational protocol developed in this study can be applied to the computational screening of antigens for the MHC1 as well as other protein-peptide binding systems.

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.&#x0D; 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).&#x0D; 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.&#x0D; 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.

Potts Hamiltonian Models and Molecular Dynamics Free Energy Simulations for Predicting the Impact of Mutations on Protein Kinase Stability.

Single-point mutations in kinase proteins can affect their stability and fitness, and computational analysis of these effects can provide insights into the relationships among protein sequence, structure, and function for this enzyme family. To assess the impact of mutations on protein stability, we used a sequence-based Potts Hamiltonian model trained on a kinase family multiple-sequence alignment (MSA) to calculate the statistical energy (fitness) effects of mutations and compared these against relative folding free energies (ΔΔGs) calculated from all-atom molecular dynamics free energy perturbation (FEP) simulations in explicit solvent. The fitness effects of mutations in the Potts model (ΔEs) showed good agreement with experimental thermostability data (Pearson r = 0.68), similar to the correlation we observed with ΔΔGs predicted from structure-based relative FEP simulations. Recognizing the possible advantages of using Potts models to rapidly estimate protein stability effects of kinase mutations seen in cancer genomics data, we used the Potts statistical energy model to estimate the stability effects of 65 conservative and nonconservative mutations across three distinct kinases (Wee1, Abl1, and Cdc7) with somatic mutations reported in the Genomic Data Commons (GDC) database. The ΔEs of these mutations calculated from the Potts model are consistent with the corresponding ΔΔGs from FEP simulations (Pearson ratio of 0.72). The agreement between these methods suggests that the Potts model may be used as a sequence-based tool for high-throughput screening of mutational effects as part of a computational pipeline for predicting the stability effects of mutations. We also demonstrate how the scalability of the fitness-based Potts model calculations permits analyses that are not easily accessed using FEP simulations. To this end, we employed site-saturation mutagenesis in the Potts model in order to investigate the relative stability effects of mutations seen in different cancer evolutionary scenarios. We used this approach to analyze the effects of drug pressure in Abl kinase by contrasting the relative fitness penalties of somatic mutations seen in miscellaneous cancer types with those calculated for mutations associated with cancer drug resistance. We observed that, in contrast to somatic mutations of Abl seen in various tumors that appear to have evolved neutrally, cancer mutations that evolved under drug pressure in Abl-targeted therapies tend to preserve enzyme stability.

Interaction of Asn297-Linked Glycan Ligands with the Fc Fragment of the Immunoglobulin Class G1: A Molecular Dynamics Simulation Study.

The allosteric modulation of the homodimeric H10-03-6 protein to glycan ligands L1 and L2, and the STAB19 protein to glycan ligands L3 and L4, respectively, has been studied by molecular dynamics simulations and free energy calculations. The results revealed that the STAB19 protein has a significantly higher affinity for L3 (-11.38 ± 2.32 kcal/mol) than that for L4 (-5.51 ± 1.92 kcal/mol). However, the combination of the H10-03-6 protein with glycan L2 (1.23 ± 6.19 kcal/mol) is energetically unfavorable compared with that of L1 (-13.96 ± 0.35 kcal/mol). Further, the binding of glycan ligands L3 and L4 to STAB19 would result in the significant closure of the two CH2 domains of the STAB19 conformation with the decrease of the centroid distances between the two CH2 domains compared with the H10-03-6/L1/L2 complex. The CH2 domain closure of STAB19 relates directly to the formation of new hydrogen bonds and hydrophobic interactions between the residues Ser239, Val240, Asp265, Glu293, Asn297, Thr299, Ser337, Asp376, Thr393, Pro395, and Pro396 in STAB19 and glycan ligands L3 and L4, which suggests that these key residues would contribute to the specific regulation of STAB19 to L3 and L4. In addition, the distance analysis revealed that the EF loop in the H10-03-6/L1/L2 model presents a high flexibility and partial disorder compared with the stabilized STAB19/L3/L4 complex. These results will be helpful in understanding the specific regulation through the asymmetric structural characteristics in the CH2 and CH3 domains of the H10-03-6 and STAB19 proteins.