Free Energy Decomposition Analysis Research Articles

Exploring protein-berberine interactions via molecular dynamics and MM/GBSA calculations

Structural modification on berberine has garnered considerable attention among medicinal chemists as a means to enhance its therapeutic potential. The dissection of macrobiomolecule-berberine interactions is pivotal for guiding such modifications. In this study, four diverse crystal complexes of protein-berberine were selected, and simulated by molecular dynamics, in which each complex was conducted five times. Subsequently, the interactions of protein-berberine were collected, and analyzed by in-house scripts. Our findings indicated that the four oxygen atoms, namely O15, O16, O17, and O18, actively participated in the formation of hydrogen bonds and water bridges; the carbon atoms within the D and E rings were able to form hydrophobic interactions more than other carbon atoms; the positively charged N7 atom predominantly participated in the formation of cation-π interactions; π-π interactions were prevalent in the A, B, and D rings, with their frequency influenced by the presence of the residues Trp, Tyr, His, and Phe in the binding pocket. In the free energy decomposition analysis, the terms ΔG_Lipo and ΔG_vdW emerged as the most significant contributors to the overall energy, with ΔG_vdW accounting for over 50% of the energy contribution. Notably, the contribution of ΔG_Packing to the free energy was substantially greater than that of ΔG_Hbond, which contributed less than 1%. These insights into molecular interactions provide valuable guidance for the structural modification of berberine, such as site selection of structural modification, group optimization with the improvement of van der Waals and hydrophobic interactions, and the trade-off between hydrogen bonding and π-packing interactions in terms of interaction weights.

Pharmacophore-Assisted Covalent Docking Identifies a Potential Covalent Inhibitor for Drug-Resistant Genotype 3 Variants of Hepatitis C Viral NS3/4A Serine Protease.

The emergence of drug-resistance-inducing mutations in Hepatitis C virus (HCV) coupled with genotypic heterogeneity has made targeting NS3/4A serine protease difficult. In this work, we investigated the mutagenic variations in the binding pocket of Genotype 3 (G3) HCV NS3/4A and evaluated ligands for efficacious inhibition. We report mutations at 14 positions within the ligand-binding residues of HCV NS3/4A, including H57R and S139P within the catalytic triad. We then modelled each mutational variant for pharmacophore-based virtual screening (PBVS) followed by covalent docking towards identifying a potential covalent inhibitor, i.e., cpd-217. The binding stability of cpd-217 was then supported by molecular dynamic simulation followed by MM/GBSA binding free energy calculation. The free energy decomposition analysis indicated that the resistant mutants alter the HCV NS3/4A-ligand interaction, resulting in unbalanced energy distribution within the binding site, leading to drug resistance. Cpd-217 was identified as interacting with all NS3/4A G3 variants with significant covalent docking scores. In conclusion, cpd-217 emerges as a potential inhibitor of HCV NS3/4A G3 variants that warrants further in vitro and in vivo studies. This study provides a theoretical foundation for drug design and development targeting HCV G3 NS3/4A.

Prediction of the binding mechanism of a selective DNA methyltransferase 3A inhibitor by molecular simulation

DNA methylation is an epigenetic mechanism that introduces a methyl group at the C5 position of cytosine. This reaction is catalyzed by DNA methyltransferases (DNMTs) and is essential for the regulation of gene transcription. The DNMT1 and DNMT3A or -3B family proteins are known targets for the inhibition of DNA hypermethylation in cancer cells. A selective non-nucleoside DNMT3A inhibitor was developed that mimics S-adenosyl-l-methionine and deoxycytidine; however, the mechanism of selectivity is unclear because the inhibitor–protein complex structure determination is absent. Therefore, we performed docking and molecular dynamics simulations to predict the structure of the complex formed by the association between DNMT3A and the selective inhibitor. Our simulations, binding free energy decomposition analysis, structural isoform comparison, and residue scanning showed that Arg688 of DNMT3A is involved in the interaction with this inhibitor, as evidenced by its significant contribution to the binding free energy. The presence of Asn1192 at the corresponding residues in DNMT1 results in a loss of affinity for the inhibitor, suggesting that the interactions mediated by Arg688 in DNMT3A are essential for selectivity. Our findings can be applied in the design of DNMT-selective inhibitors and methylation-specific drug optimization procedures.

BINDING AFFINITIES OF SANGGENON DERIVATIVES AS PTP1B INHIBITORS; USING MOLECULAR DYNAMICS AND FREE ENERGY CALCULATIONS

Recently, protein tyrosine phosphatase 1B (PTP1B) inhibitors have become the frontier as possible targeting for anti-cancer and antidiabetic drugs. The contemporary observe represents a pc assisted version to investigate the importance of precise residues within the binding web site of PTP1B with numerous Sanggenon derivatives remoted from nature. Molecular dynamics (MD) simulations were performed to estimate the dynamics of the complexes, and absolute binding unfastened energies have been calculated with exclusive additives, and carried out through the usage of the Molecular Mechanics-Poisson-Boltzmann floor region (MM-PB/SA) and Generalized Born surface vicinity (MM-GB/SA) strategies. The effects show that the expected free energies of the complexes are normally constant with the available experimental statistics. MM/GBSA free energy decomposition analysis shows that the residues Asp29, Arg24, Met258, and , Arg254 in the second active site in PTP1B are crucial for the excessive selectivity of the inhibitors

Decoding molecular mechanism of species-selective targeting of fungal versus human HSP90 using multiple replica molecular dynamics simulations and binding free energy calculations

As a highly evolutionarily conserved molecular chaperone, heat shock protein (HSP90), plays an important role in virulence traits, representing a therapeutic target for the treatment of fungal infections. The close evolutionary relationship between fungi and their human hosts poses a key challenge for the development of selective antifungal agents. In this work, molecular docking, multiple replica microsecond-based molecular dynamics (MD) simulations, and binding free energy calculations were performed to decode molecular mechanism of species-selective targeting of fungal versus human HSP90 triggered by the compound A11. MD simulations reveal that binding of compound A11 to human HSP90 nucleotide-binding domain (NBD) leads to obvious conformational changes relative to fungal HSP90 NBD. Binding free energy calculations show that the binding of compound A11 to fungal HSP90 NBD is stronger than that to human HSP90 NBD. Per residue-based free energy decomposition analysis was used to evaluate the inhibitor − residue interaction profile. The results efficiently identify the hot spot residues that play vital roles in favorable binding of compound A11 to fungal HSP90 NBD. This study is expected to provide a useful guidance for the development of selective inhibitors toward fungal HSP90. Communicated by Ramaswamy H. Sarma

Molecular docking and molecular dynamics simulations reveal the clinical resistance of the gatekeeper mutation V564F of FGFR2 against Infigratinib

Fibroblast growth factor receptor 2 (FGFR2), as a transmembrane receptor tyrosine kinase, is implicated in a plethora of human cancers, including intrahepatic cholangiocarcinomas, breast cancers, and non-small cell lung cancer. The clinically relevant V564F gatekeeper mutation conferred resistance to current FGFR2 drug − Infigratinib. In this study, the protein − ligand interactions between FGFR2 kinase domain (wild-type and V564F) and Infigratinib were compared through an integrated computational method. The multiple, large-scale molecular dynamics (MD) simulations, together with dynamic cross-correlation analysis and binding free energy calculations suggested that the resistant mutation may not trigger the conformational changes of the FGFR2 kinase domain. The simulation results also indicated that the driving force to decrease the binding affinity of Infigratinib to the FGFR2 V564F variant derived from the difference in the protein − ligand hydrogen bonding interactions. Moreover, the per-residue free energy decomposition analysis revealed that the reduced contributions from several residues in the ATP-binding site of FGFR2, especially Glu565 and Ala567 located at the kinase hinge domain, were the key residues responsible for the loss of binding affinity of Infigratinib to the V564F variant. This study may offer useful information for the design of novel selective kinase inhibitors targeting FGFR2.

Interface-Based Design of High-Affinity Affibody Ligands for the Purification of RBD from Spike Proteins.

The outbreak of coronavirus disease 2019 (COVID-19) has sparked an urgent demand for advanced diagnosis and vaccination worldwide. The discovery of high-affinity ligands is of great significance for vaccine and diagnostic reagent manufacturing. Targeting the receptor binding domain (RBD) from the spike protein of severe acute respiratory syndrome-coronavirus 2, an interface at the outer surface of helices on the Z domain from protein A was introduced to construct a virtual library for the screening of ZRBD affibody ligands. Molecular docking was performed using HADDOCK software, and three potential ZRBD affibodies, ZRBD-02, ZRBD-04, and ZRBD-07, were obtained. Molecular dynamics (MD) simulation verified that the binding of ZRBD affibodies to RBD was driven by electrostatic interactions. Per-residue free energy decomposition analysis further substantiated that four residues with negative-charge characteristics on helix α1 of the Z domain participated in this process. Binding affinity analysis by microscale thermophoresis showed that ZRBD affibodies had high affinity for RBD binding, and the lowest dissociation constant was 36.3 nmol/L for ZRBD-07 among the three potential ZRBD affibodies. Herein, ZRBD-02 and ZRBD-07 affibodies were selected for chromatographic verifications after being coupled to thiol-activated Sepharose 6 Fast Flow (SepFF) gel. Chromatographic experiments showed that RBD could bind on both ZRBD SepFF gels and was eluted by 0.1 mol/L NaOH. Moreover, the ZRBD-07 SepFF gel had a higher affinity for RBD. This research provided a new idea for the design of affibody ligands and validated the potential of affibody ligands in the application of RBD purification from complex feedstock.

Investigating the novel-binding site of RPA2 on Menin and predicting the effect of point mutation of Menin through protein–protein interactions

Protein–protein interactions (PPIs) play a critical role in all biological processes. Menin is tumor suppressor protein, mutated in multiple endocrine neoplasia type 1 syndrome and has been shown to interact with multiple transcription factors including (RPA2) subunit of replication protein A (RPA). RPA2, heterotrimeric protein required for DNA repair, recombination and replication. However, it’s still remains unclear the specific amino acid residues that have been involved in Menin-RPA2 interaction. Thus, accurately predicting the specific amino acid involved in interaction and effects of MEN1 mutations on biological systems is of great interests. The experimental approaches for identifying amino acids in menin-RPA2 interactions are expensive, time-consuming, and challenging. This study leverages computational tools, free energy decomposition and configurational entropy scheme to annotate the menin-RPA2 interaction and effect on menin point mutation, thereby proposing a viable model of menin-RPA2 interaction. The menin–RPA2 interaction pattern was calculated on the basis of different 3D structures of menin and RPA2 complexes, constructed using homology modeling and docking strategy, generating three best-fit models: Model 8 (− 74.89 kJ/mol), Model 28 (− 92.04 kJ/mol) and Model 9 (− 100.4 kJ/mol). The molecular dynamic (MD) was performed for 200 ns and binding free energies and energy decomposition analysis were calculated using Molecular Mechanics Poisson–Boltzmann Surface Area (MM/PBSA) in GROMACS. From binding free energy change, model 8 of Menin-RPA2 exhibited most negative binding energy of − 205.624 kJ/mol, followed by model 28 of Menin-RPA2 with − 177.382 kJ/mol. After S606F point mutation in Menin, increase of BFE (ΔGbind) by − 34.09 kJ/mol in Model 8 of mutant Menin-RPA2 occurs. Interestingly, we found a significant reduction of BFE (ΔGbind) and configurational entropy by − 97.54 kJ/mol and − 2618 kJ/mol in mutant model 28 as compared the o wild type. Collectively, this is the first study to highlight the configurational entropy of protein–protein interactions thereby strengthening the prediction of two significant important interaction sites in menin for the binding of RPA2. These predicted sites could be vulnerable for structural alternation in terms of binding free energy and configurational entropy after missense mutation in menin.

A Free Energy Decomposition Analysis: Insight into Binding Thermodynamics from Absolutely Localized Molecular Orbitals.

Traditional energy decomposition analysis (EDA) methods can provide an interpretive decomposition of non-covalent electronic binding energies. However, by definition, they neglect entropic effects and nuclear contributions to the enthalpy. With the objective of revealing the chemical origins of trends in free energies of binding, we introduce the concept of a Gibbs decomposition analysis (GDA) by coupling the absolutely localized molecular orbital treatment of electrons in non-covalent interactions with the simplest possible quantum rigid rotor-harmonic oscillator treatment of nuclear motion at finite temperature. The resulting pilot GDA is employed to decompose enthalpic and entropic contributions to the free energy of association of the water dimer, fluoride-water dimer, and water binding to an open metal site in the metal-organic framework Cu(I)-MFU-4l. The results show enthalpic trends that generally track the electronic binding energy and entropic trends that reveal the increasing price of loss of translational and rotational degrees freedom with temperature.

Demystifying Solvent Effects on Diels–Alder Reactions in Pure and Mixed Solvents: A Combined Electronic DFT and QM/MM Study

Although solvent effects on chemical reactions in pure solvents are well elucidated, both reaction mechanisms and solvent effects in mixed solvents remain poorly understood. Herein, the free energy profiles and the solute–solvent energy pair distribution of the Diels–Alder reaction in pure water, acetonitrile (ACN), carbon tetrachloride (CTC), and their mixtures ACN-water and ACN-CTC are calculated using the quantum and molecular mechanical (QM/MM) Monte Carlo method. The results show that the presence of solvent can significantly reduce the reaction free energy barrier, and the barrier in the mixed solvent is intermediate between the two pure solvents. Additionally, the free energy decomposition analysis demonstrated that the QM region dominates the reaction free energy. Electronic DFT calculations are employed to delineate the interactions between reactants and solvent molecules in pure and mixed solvents, with a particular focus on the role of electrostatic potential characteristics and an independent gradient model based on Hirshfeld partition. By rigorously analyzing these interactions, this study demystifies the mechanism underlying the phenomenon of preferential solvation and competition for reacting sites in mixed solvents, which alters the microenvironment of the solute.

Insights Into Resveratrol as an Inhibitor Against Aβ1-42 Peptide Aggregation: A Molecular Dynamics Simulation Study

Background: Resveratrol (RSV), a polyphenolic compound, is reported to have anti-aggregation properties against Amyloid-beta peptides. It is, therefore, significant to understand the mechanism of inhibition of Aβ1-42 peptide aggregation by the RSV at the molecular level. We have used Molecular docking along with Molecular dynamics (MD) simulation techniques to address the role of RSV in the inhibition of Aβ1-42 pep-tide aggregation. Objective: To understand the role of Resveratrol on the Aβ1-42 peptide aggregation. Method: In this computational study, we have docked the RSV to Aβ1-42 peptide using Molecular Docking software and then performed MD simulation for the Aβ1-42 peptide monomer Aβ1-42 peptide-RSV complex using the AMBER force field. From the analysis of MD trajectories, we obtained salient structural features and determined the Binding Free Energy(BFE) and Per-residue Energy Decomposition Analysis (PRED) using MM-PBSA/GBSA method. Results: The secondary structure and the conformational analysis obtained from MD trajectories show that the binding of RSV with the Aβ1-42 peptide monomer causes an increase in the helical content in the structure of the Aβ1-42 peptide. The BFE and PRED results show a high binding affinity (GBtotal=-11.07 kcal mol-1; PBtotal= -1.82 kcal mol-1) of RSV with Aβ1-42 peptide. Also, we found the RSV to interact with cru-cial residues (Asp 23 and Lys 28) of the Aβ1-42 peptide. These residues play a signif-icant role in facilitating the formation of toxic amyloid oligomers and amyloid fibrils. The salt bridge interaction between these residues D23–K28 was found to be destabi-lized in the Aβ1-42 peptide when it is complexed with RSV. Conclusion: In summary, it can be concluded that Resveratrol greatly aids the preven-tion of Aβ1-42 peptide aggregation. Therefore, it can be considered a possible drug candidate for therapeutic strategies for Alzheimer’s disease.

Molecular insights on PS-PLA1 lipase activity of human ABHD16B

The human alpha beta hydrolase domain (ABHD) proteins are ubiquitous and regulate the cellular lipids' anabolic and catabolic processes. The structural aspects for specific biochemical function of many ABHD proteins related to physiological disorders and its link to pathological conditions remain unknown. Here putative human ABHD16B protein was overexpressed in Saccharomyces cerevisiae for its biological activity. In-vitro enzymatic assay of the recombinant ABHD16B protein with fluorescently tagged glycerophospholipids revealed that the PLA1 activity is observed with phosphatidylserine (PS). In addition, it efficiently hydrolyzed monoacylglycerol over triacylglycerols. Further, molecular dynamic simulations and per residue binding free energy decomposition analysis revealed that the origin of PS-specific PLA1 activity of ABHD16B is due to the electrostatic interaction of the PS head group with K8, R319, and E178, which led to having the hydrogen bond interaction of sn-1 acyl chain ester to the catalytic site residues. Site-directed mutagenesis of the 245GXSXG249 motif of ABHD16B reduced the maximal lipase activity of PS and MAG. In summary, these results revealed that ABHD16B plays a vital role in PS selectivity that in turn, controls the specific subcellular pools of 2-LPS metabolism in the tissues at low pH.

Phosphorylation of Tau R2 Repeat Destabilizes Its Binding to Microtubules: A Molecular Dynamics Simulation Study.

Phosphorylation, the most popular post-translational modification of tau protein, plays an important role in regulating tau physiological functions. However, aberrant phosphorylation attenuates the binding affinity of tau to a microtubule (MT), resulting in MT destabilization followed by accumulation of neurofibrillary tangles in the brain. There are in total 85 potential phosphorylation sites in a full-length tau protein, and about half of them are abnormal as they occur in tau of Alzheimer's disease (AD) brain only. In this work, we investigated the impact of abnormal Ser289, Ser293, and Ser289/Ser293 phosphorylation on tau R2-MT binding and the conformation of tau R2 using molecular dynamics simulation. We found that the phosphorylation significantly affected R2-MT interaction and reduced the binding affinity of tau R2 peptides to MTs. Free energy decomposition analysis suggested that the post-translational modified residues themselves made a significant contribution to destabilize tau repeat R2-MT binding. Therefore, the phosphorylation may attenuate the binding affinity of tau to MTs. Additionally, the phosphorylation also enhanced helix-coil transition of monomeric R2 peptides, which may result in the acceleration of tau aggregation. Since these phosphorylated sites have not been examined in previous experimental studies, our finding through all-atom molecular dynamics simulations and free energy analysis can inspire experimental scientists to investigate the impact of the phosphorylation on MT binding and aggregation of full-length tau and the pathological roles of the phosphorylation at those sites in AD development through in vitro/in vivo assays.

Evaluation of interactions between the hepatitis C virus NS3/4A and sulfonamidobenzamide based molecules using molecular docking, molecular dynamics simulations and binding free energy calculations.

The Hepatitis C Virus (HCV) NS3/4A is an attractive target for the treatment of Hepatitis C infection. Herein, we present an investigation of HCV NS3/4A inhibitors based on a sulfonamidobenzamide scaffold. Inhibitor interactions with HCV NS3/4A were explored by molecular docking, molecular dynamics simulations, and MM/PBSA binding free energy calculations. All of the inhibitors adopt similar molecular docking poses in the catalytic site of the protease that are stabilized by hydrogen bond interactions with G137 and the catalytic S139, which are known to be important for potency and binding stability. The quantitative assessments of binding free energies from MM/PBSA correlate well with the experimental results, with a high coefficient of determination, R2 of 0.92. Binding free energy decomposition analyses elucidate the different contributions of Q41, F43, H57, R109, K136, G137, S138, S139, A156, M485, and Q526 in binding different inhibitors. The importance of these sidechain contributions was further confirmed by computational alanine scanning mutagenesis. In addition, the sidechains of K136 and S139 show crucial but distinct contributions to inhibitor binding with HCV NS3/4A. The structural basis of the potency has been elucidated, demonstrating the importance of the R155 sidechain conformation. This extensive exploration of binding energies and interactions between these compounds and HCV NS3/4A at the atomic level should benefit future antiviral drug design.

The effect of genotype variation and M423 resistance mutations to the binding of phosphonomidate-based inhibitor IDX17119 with the thumb-II domain of Hepatitis C virus RdRp: an integrated molecular dynamics and binding free energy study

ABSTRACT The IDX17119 is a phosphonomidate-based pan-genotypic non-nucleoside inhibitor of Hepatitis C Virus (HCV) NS5B RdRp thumb-II domain. IDX17119 exhibits significant antiviral activity against HCV genotype (GT) 1a, 2a, 3a and 4a. The molecular docking analysis was performed to understand the binding and intermolecular interactions of the IDX17119 molecule at the thumb-II domain of 1a, 2a, 3a and 4a GT NS5B RdRps. The molecular dynamics simulation has been performed to determine the stability of the ligand–protein complexes. We found GT 1a complex is relatively stable compared with other drug-receptor complexes. The binding free energy calculations and free energy decomposition analysis were carried out for all the complexes, the result revealed that the increase in the unfavourable polar contribution affects the binding of IDX17119 with GT 2a NS5B RdRp compared with all other complexes. In addition to this, the effect of mutation at M423 site (M423T/V) alters the conformation and activity of the IDX17119 molecule.

Engineering human JMJD2A tudor domains for an improved understanding of histone peptide recognition.

JMJD2A is a histone lysine demethylase which recognizes and demethylates H3K9me3 and H3K36me3 residues and is overexpressed in various cancers. It utilizes a tandem tudor domain to facilitate its own recruitment to histone sites, recognizing various di- and tri-methyl lysine residues with moderate affinity. In this study, we successfully engineered the tudor domain of JMJD2A to specifically bind to H4K20me3 with a 20-fold increase of affinity and improved selectivity. To reveal the molecular basis, we performed molecular dynamics and free energy decomposition analysis on the human JMJD2A tandem tudor domains bound to H4K20me2, H4K20me3, and H3K23me3 peptides to uncover the residues and conformational changes important for the enhanced binding affinity and selectivity toward H4K20me2/3. These analyses revealed new insights into understanding chromatin reader domains recognizing histone modifications and improving binding affinity and selectivity of these domains. Furthermore, we showed that the tight binding of JMJD2A to H4K20me2/3 is not sufficient to improve the efficiency of CRISPR-CAS9 mediated homology directed repair (HDR), suggesting a complicated relationship between JMJD2A and the DNA damage response beyond binding affinity toward the H4K20me2/3 mark.

Studies on the antibacterial activities and molecular mechanism of GyrB inhibitors by 3D-QSAR, molecular docking and molecular dynamics simulation

DNA gyrase B (GyrB) plays a critical role in DNA replication, repair, recombination, and transcription and has become an attractive target for a number of antibacterial agents. In the present work, to better comprehend the structure–activity relationship and mechanism of action of novel series of GyrB inhibitors, a theoretical study of three-dimensional-quantitative structure–activity relationship (3D-QSAR), molecular docking, molecular dynamics (MD) simulations, and binding free energy analysis were performed. The results showed that CoMFA (R 2 cv = 0.591, R 2 pred = 0.7698) and CoMSIA (R 2 cv = 0.629, R 2 pred = 0.6848) models possessed robust stability and predictability. The derived contour maps of steric, electrostatic, and hydrogen bond donor field further displayed the modified information of these inhibitors. Molecular docking was further conducted to provide the initial binding conformation for MD simulations. Subsequently, MD simulations were applied to explore the key amino acids and binding modes at the active site. The binding free energy decomposition analysis further indicates that the residues Ile54, Glu55, Arg83, Ala85, Val86 and Thr128 are essential for the high selectivity of inhibitors. Overall, these results would serve as a significant guideline for the discovery and design of novel GyrB inhibitors.

Deciphering the Mechanism of Binding Selectivity of Chlorofluoroacetamide-Based Covalent Inhibitors toward L858R/T790M Resistance Mutation.

Covalent modification of the oncogenic mutant epidermal growth factor receptor (EGFR) by small molecules is an efficient strategy for achieving an enhanced and sustained pharmacological effect in the treatment of non-small-cell lung cancer. NSP-037 (18), an irreversible inhibitor of the L858R/T790M double-mutant EGFR (EGFRDM) using α-chlorofluoroacetamide (CFA) as a novel warhead, has seven times the inhibition selectivity for EGFRDM over the wild type (EGFRWT), as compared to clinically approved osimertinib (7). Here, we employ multiple computational approaches to elucidate the mechanism underlining this improved selectivity, as well as the effect of CFA on the selectivity enhancement of inhibitor 18 over 7. We find that EGFRDM undergoes significantly larger conformational changes than EGFRWT upon binding to 18. The conformational stability of the diamine side chain and the CFA motif of 18 in the orthosteric site of EGFRDM is identified as key for the disparate binding mechanism and inhibitory prowess of 18 with respect to EGFRWT and EGFRDM and 18's higher selectivity than 7. The binding free energy of the 18-bound complexes is -6.38 kcal/mol greater than that of the 7-bound complexes, explaining the difference in selectivity of these inhibitors. Further, free energy decomposition analysis indicates that the electrostatic contribution of key residues plays an important role in the 18-bound complexes. QM/MM calculations show that the most favored mechanism for the Cys797 alkylation reaction is the direct displacement mechanism through a CFA-based inhibitor, producing a reaction with the lowest energy barrier and most stable product.

Simulation analysis of selective alanine mutation effect on stability of human prion protein

Prion diseases are neurodegenerative disorders caused by spongiform degeneration of the brain. Understanding the fundamental mechanism of prion protein aggregation caused by mutations is very crucial to resolve the pathology of prion diseases. To help understand the roles of individual residues on the stability of the human prion protein, the computational method of free energy simulations based on atomistic molecular dynamics trajectories is applied to Phe175 → Ala, Val180 → Ala, and Val209 → Ala mutations of the human prion protein. The simulations show that all three alanine mutations destabilize the human prion protein. The calculated free energy change differences, ΔΔG, for the Phe175 → Ala, Val180 → Ala, and Val209 → Ala mutations are in good agreement with the experimental values. The significant destabilizing effects on the mutants relative to the wild-type protein arise from van der Waals terms. Furthermore, our free energy decomposition analysis shows that the major contribution to destabilizing the V180A and V209A mutants relative to the wild-type protein is originated from van der Waals interactions from residues near the mutation sites. In contrast, the contribution to destabilizing the F175A mutant is mainly caused by van der Waals interactions from residues near and far away from the mutation site. Our results show that the free energy simulation with a thermodynamic integration approach for selected alanine scanning mutations is beneficial for understanding the detailed mechanism of human prion protein destabilization, specific residues' role, and the hydrophobic effect on protein stability. Communicated by Ramaswamy H. Sarma

Modelling the DFT structural and reactivity study of feverfew and evaluation of its potential antiviral activity against COVID-19 using molecular docking and MD simulations.

The unavailability of a proper drug against SARS-CoV-2 infections and the emergence of various variants created a global crisis. In the present work, we have studied the antiviral behavior of feverfew plant in treating COVID-19. We have reported a systematic in silico study with the antiviral effects of various phytoconstituents Borneol (C10H18O), Camphene (C10H16), Camphor (C10H16O), Alpha-thujene (C10H16), Eugenol (C10H14O), Carvacrol (C10H14O) and Parthenolide (C15H20O3) of feverfew on the viral protein of SARS-CoV-2. Parthenolide shows the best binding affinity with both main protease (Mpro) and papain-like protease (PLpro). The molecular electrostatic potential and Mulliken atomic charges of the Parthenolide molecule shows the high chemical reactivity of the molecule. The docking of Parthenolide with PLpro give score of −8.0 kcal/mol that validates the good binding of Parthenolide molecule with PLpro. This complex was further considered for molecular dynamics simulations. The binding energy of the complex seems to range in between −3.85 to −11.07 kcal/mol that is high enough to validate the stability of the complex. Free energy decomposition analysis have been also performed to understand the contribution of residues that reside into the binding site. Good binding affinity and reactivity response suggested that Parthenolide can be used as a promising drug against the COVID-19.Graphical Supplementary InformationThe online version contains supplementary material available at 10.1007/s11696-022-02067-6.