Targeting the N-terminal domain of the RNA-binding protein of the SARS-CoV-2 with high affinity natural compounds to abrogate the protein-RNA interaction: a molecular dynamics study

Calcium-dependent lectin I from Pseudomonas aeruginosa (PA-IL) binds specifically to oligosaccharides presenting an alpha-galactose residue at their nonreducing end, such as the disaccharides alphaGal1-2betaGalOMe, alphaGal1-3betaGalOMe, and alphaGal1-4betaGalOMe. This provides a unique model for studying the effect of the glycosidic linkage of the ligands on structure and thermodynamics of the complexes by means of experimental and theoretical tools. The structural features of PA-IL in complex with the three disaccharides were established by docking and molecular dynamics simulations and compared with those observed in available crystal structures, including PA-IL.alphaGal1-2betaGalOMe complex, which was solved at 2.4 A resolution and reported herein. The role of a structural bridge water molecule in the binding site of PA-IL was also elucidated through molecular dynamics simulations and free energy calculations. This water molecule establishes three very stable hydrogen bonds with O6 of nonreducing galactose, oxygen from Pro-51 main chain, and nitrogen from Gln-53 main chain of the lectin binding site. Binding free energies for PA-IL in complex with the three disaccharides were investigated, and the results were compared with the experimental data determined by titration microcalorimetry. When the bridge water molecule was included in the free energy calculations, the simulations predicted the correct binding affinity trends with the 1-2-linked disaccharide presenting three times stronger affinity ligand than the other two. These results highlight the role of the water molecule in the binding site of PA-IL and indicate that it should be taken into account when designing glycoderivatives active against P. aeruginosa adhesion.

Molecular dynamics simulations and free energy calculation have been performed to study how the single-chain variable fragment (scFv) binds methamphetamine (METH) and amphetamine (AMP). The structures of the scFv:METH and the scFv:AMP complexes are analyzed by examining the time-dependence of their RMSDs, by analyzing the distance between some key atoms of the selected residues, and by comparing the averaged structures with their corresponding crystallographic structures. It is observed that binding an AMP to the scFv does not cause significant changes to the binding pocket of the scFv:ligand complex. The binding free energy of scFv:AMP without introducing an extra water into the binding pocket is much stronger than scFv:METH. This is against the first of the two scenarios postulated in the experimental work of Celikel et al. (Protein Science 18, 2336 (2009)). However, adding a water to the AMP (at the position of the methyl group of METH), the binding free energy of the scFv:AMP-H2O complex, is found to be significantly weaker than scFv:METH. This is consistent with the second of the two scenarios given by Celikel et al. Decomposition of the binding energy into ligand-residue pair interactions shows that two residues (Tyr175 and Tyr177) have nearly-zero interactions with AMP in the scFv:AMP-H2O complex, whereas their interactions with METH in the scFv:METH complex are as large as -0.8 and -0.74 kcal mol(-1). The insights gained from this study may be helpful in designing more potent antibodies in treating METH abuse.

4-(Phenylamino)-pyrrolo[2,1-f][1,2,4]triazines have been discovered as inhibitors of p38alpha. Experimental assays have proven that the configuration of alpha-Me-benzyl connected with amide at C6 is essential for the binding affinity. The S-configured inhibitor (11j) displays 80 times more potency than the R-configured one (11k). Here we investigated the mechanism how different configurations influence the binding affinity using molecular dynamics simulations, free energy calculations and free energy decomposition analysis. We found that the van der Waals interactions play the most important role in differentiating the activities between 11j and 11k with p38alpha. The difference of the van der Waals interactions is primarily determined by two residues, LEU108 and LEU167. Consequently stabilization of pyrrolo[2,1-f][1,2,4]triazine ring is important for the activities of inhibitors. Meanwhile we observed that the different configuration of the alpha-Me-benzyl group leads to the difference of binding between 11j and 11k. In conclusion, our work shows that it is feasible to analyze the chirality effect of inhibitors with different configurations by molecular dynamics simulations and free energy calculations, and provides useful information for drug design.

Bidentate inhibitors of protein tyrosine phosphatase 1B (PTP1B) are considered as a group of ideal inhibitors with high binding potential and high selectivity in treating type II diabetes. In this paper, the binding models of five bidentate inhibitors to PTP1B, TCPTP, and SHP-2 were investigated and compared by using molecular dynamics (MD) simulations and free energy calculations. The binding free energies were computed using the Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) methodology. The calculation results show that the predicted free energies of the complexes are well consistent with the experimental data. The Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) free energy decomposition analysis indicates that the residues ARG24, ARG254, and GLN262 in the second binding site of PTP1B are essential for the high selectivity of inhibitors. Furthermore, the residue PHE182 close to the active site is also important for the selectivity and the binding affinity of the inhibitors. According to our analysis, it can be concluded that in most cases the polarity of the portion of the inhibitor that binds to the second binding site of the protein is positive to the affinity of the inhibitors while negative to the selectivity of the inhibitors. We expect that the information we obtained here can help to develop potential PTP1B inhibitors with more promising specificity.

Computational simulation techniques have become increasingly powerful and widespread due to advancements in computing capabilities which provide valu- able insights into the behavior, properties, and interactions of nanoscale systems to accelerate development, optimize designs, and improve the safety and efficacy of nanomaterial-based applications. In the current study, adsorption equilibrium coefficients and free energy of some aromatic acids are calculated at different temperatures using molecular dynamics simulation along with free enegrgy calcu- lation to investigate the adsorption behavior of aromatic acid on graphene surface with temperature. The effect of temperature on the adsorption behavior of cer- tain aromatic acids -salicylic acid, benzoic acid, and phthalic acid- on carbon nanotubes was investigated using molecular dynamics simulation and free energy calculations at three different temperatures 260K, 300K, and 340K. Current anal- yses show that, in all temperature, salicylic acid has the highest adsorption and phthalic acid has the lowest (phthalic acid < benzoic acid < salicylic acid). It is discovered that the amount of carboxylic acids that adsorb on carbon nanotubes increases with temperature. Keywords: carbon nanotube, aromatic compounds, free energy calculation, adsorption equilibrium constant

We have studied the binding of two organic cations, an iminium (IM) and a guanidinium (GU), to a cyclophane host P4--4Na+, using molecular dynamics simulations and free energy calculations. A proper treatment of the long-range electrostatic forces is essential for the stability of these highly charged complexes, and a simple cutoff at 12 A results in an artifactual dissociation of the IM−P4--4Na+ complex. Since the host is highly aromatic and the guests cationic, cation−π interactions play an important role in the complex stability. In free energy calculations, using a simple additive force field, we calculate that the relative free energy of association of IM and GU binding to the host is 2.3 kcal/mol favoring IM, which is of the correct sign but 1.4 kcal/mol too small in magnitude. Differences in van der Waals interaction energies are mainly responsible for the different binding strengths, and the host adopts different shapes when accommodating IM compared to GU. To approximately estimate the contributi...

Molecular dynamics simulation and free energy calculation studies of kinase inhibitors binding to active and inactive conformations of VEGFR-2

ABTRACT The epidermal growth factor receptor (EGFR) kinase inhibitors Gefitinib, Erlotinib, Afatinib and Osimertinib have been approved for the treatments of non-small cell lung cancer patients harboring sensitive EGFR mutations, but resistance arises rapidly. To date all approved EGFR inhibitors are ATP-competitive inhibitors, highlighting the need for therapeutic agents with alternative mechanisms of action. Allosteric kinase inhibitors offer a promising new therapeutic strategy to ATP-competitive inhibitors. The mutant-selective allosteric EGFR inhibitors EAI045 exhibited higher potency for EGFRL858R&T790M compared to WT, which was also effective in EGFR-mutant models including those harboring the C797S mutation. However, it was not effective as a single-agent inhibitor, and require the co-administration of the anti-EGFR antibody Cetuximab. Further efforts produced a more potent analog JBJ-04-125-02, which can inhibit cell proliferation as a single-agent inhibitor. In the present study, molecular dynamics simulations and free energy calculations were performed and revealed the detailed inhibitory mechanism of JBJ-04-125-02 as more potent EGFR inhibitor. Moreover, the energy difference between HOMO and LUMO calculated by DFT implied the higher interaction of JBJ-04-125-02 than EAI045 in the active site of the EGFR. The identified key features obtained from the molecular modeling enabled us to design novel EGFR allosteric inhibitors. Communicated by Ramaswamy H. Sarma

Lung cancer is the leading cause of cancer death, and epidermal growth factor receptor (EGFR) kinase domain mutations are a common cause of non-small-cell lung cancer (NSCLC), a major subtype of lung cancers. Patients harboring most of these mutations respond well to the EGFR inhibitors Gefitinib and Erlotinib initially, but soon develop resistance to them due to the emergence of the gatekeeper mutation T790M. The new-generation inhibitors such as AZD9291, HM61713, CO-1686 and WZ4002 can overcome T790M through covalent binding to Cys 797, but ultimately lose their efficacy upon the emergence of the C797S mutation that abolishes the covalent bonding. Allosteric inhibitors EAI001 and EAI045 are a new type of EGFR inhibitors that bind to EGFR away from the ATP-binding site and not relying on Cys 797. In this study, molecular dynamics simulations and free energy calculations were carried out on EAI001 and EAI045 in complex with EGFR, revealing the detailed inhibitory mechanism of EAI001 and EAI045 as EGFR allosteric inhibitor, which was expected to provide a basis for rational drug design of the EGFR allosteric inhibitors.Communicated by Ramaswamy H. Sarma

4-Hydroxyphenylpyruvate dioxygenase is a relevant target in both pharmaceutical and agricultural research. We report on molecular dynamics simulations and free energy calculations on this enzyme, in complex with 12 inhibitors for which experimental affinities were determined. We applied the thermodynamic integration approach and the more efficient one-step perturbation. Even though simulations seem well converged and both methods show excellent agreement between them, the correlation with the experimental values remains poor. We investigate the effect of slight modifications on the charge distribution of these highly conjugated systems and find that accurate models can be obtained when using improved force field parameters. This study gives insight into the applicability of free energy methods and current limitations in force field parameterization.

Matrix metalloproteinases (MMP) are well-known biological targets implicated in tumour progression, homeostatic regulation, innate immunity, impaired delivery of pro-apoptotic ligands, and the release and cleavage of cell-surface receptors. Hence, the development of potent and selective inhibitors targeting these enzymes continues to be eagerly sought. In this paper, a number of alloxan-based compounds, initially conceived to bias other therapeutically relevant enzymes, were rationally modified and successfully repurposed to inhibit MMP-2 (also named gelatinase A) in the nanomolar range. Importantly, the alloxan core makes its debut as zinc binding group since it ensures a stable tetrahedral coordination of the catalytic zinc ion in concert with the three histidines of the HExxHxxGxxH metzincin signature motif, further stabilized by a hydrogen bond with the glutamate residue belonging to the same motif. The molecular decoration of the alloxan core with a biphenyl privileged structure allowed to sample the deep S1′ specificity pocket of MMP-2 and to relate the high affinity towards this enzyme with the chance of forming a hydrogen bond network with the backbone of Leu116 and Asn147 and the side chains of Tyr144, Thr145 and Arg149 at the bottom of the pocket. The effect of even slight structural changes in determining the interaction at the S1′ subsite of MMP-2 as well as the nature and strength of the binding is elucidated via molecular dynamics simulations and free energy calculations. Among the herein presented compounds, the highest affinity (pIC50 = 7.06) is found for BAM, a compound exhibiting also selectivity (>20) towards MMP-2, as compared to MMP-9, the other member of the gelatinases.

Designing small molecule inhibitors that are highly selective for Golgi alpha mannosidase II (GMII) over lysosomal alpha mannosidase (LM) remains crucial for the development of novel anticancer drugs targeting the N-glycosylation pathway. Studies have previously identified an unconserved ‘anchor site’ in GMII that represents an attractive target for achieving selectivity. In this study we conduct molecular dynamics simulations and free energy calculations of GMII and LM with their natural oligosaccharide substrates to investigate the potential of the anchor site. Our findings reveal that the corresponding N-acetylglucosamine residue remains tightly bound to the anchor site in GMII, helping stabilize overall substrate binding to GMII, while the lack of a similar conserved site in LM allows the substrate to fluctuate more freely. Our simulations also suggest that besides stabilizing the catalytic site mannose for cleavage, the anchor site residue may play a role in facilitating the release of the holding site mannose and its transition to the catalytic site in GMII. Subsequently, free energy calculations reveal that the anchor site N-acetylglucosamine contributes 4.107 kcal/mol to the free energy of binding in GMII, with a significantly smaller contribution of 1.035 kcal/mol in LM. This difference of 3.072 kcal/mol in the free energy of binding represents potential gains in selectivity that could be achieved by targeting the anchor site in GMII. Taken together, these findings provide further evidence on the potential of targeting the anchor site in the design of highly selective GMII inhibitors.

Human cytochrome P450 2E1 (CYP2E1) participates in the metabolism of over 2% of all the oral drugs. A hallmark peculiar feature of this enzyme is that it exhibits a pronounced negative cooperativity in substrate binding. However the mechanism by which the negative cooperativity occurs is unclear. Here, we performed molecular dynamics simulations and free energy calculations on human CYP2E1 to examine the structural differences between the substrate-free and the enzymes with one and two aniline molecules bound. Our results indicate that although the effector substrate does not bind in the active site cavity, it still can directly interact with the active site residues of human CYP2E1. The interaction of the effector substrate with the active site leads to a reorientation of active site residues, which thereby weakens the interactions of the active substrate with this site. We also identify a conserved residue T303 that plays a crucial role in the negative cooperative binding on the short-range effects. This residue is a key factor in the positioning of substrates and in proton delivery to the active site. Additionally, a long-range effect of the effector substrate is identified in which F478 is proposed to play a key role. As located in the interface between the active and effector sites, this residue structurally links the active and effector sites and is found to play a significant role in affecting substrate access and ligand positioning within the active site. In the negative cooperative binding, this residue can decrease the interactions of the active substrate with the active site by π-π stacking which then lowers the hydroxylation activity for the active substrate. These findings are in agreement with previous experimental observations and thus provide detailed atomistic insight into the poorly understood mechanism of the negative cooperativity in human CYP2E1.

The recent outbreak of novel “coronavirus disease 2019” (COVID-19) has spread rapidly worldwide, causing a global pandemic. In the present work, we have elucidated the mechanism of binding of two inhibitors, namely α-ketoamide and Z31792168, to SARS-CoV-2 main protease (Mpro or 3CLpro) by using all-atom molecular dynamics simulations and free energy calculations. We calculated the total binding free energy (ΔGbind) of both inhibitors and further decomposed ΔGbind into various forces governing the complex formation using the Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) method. Our calculations reveal that α-ketoamide is more potent (ΔGbind= − 9.05 kcal/mol) compared to Z31792168 (ΔGbind= − 3.25 kcal/mol) against COVID-19 3CLpro. The increase in ΔGbind for α-ketoamide relative to Z31792168 arises due to an increase in the favorable electrostatic and van der Waals interactions between the inhibitor and 3CLpro. Further, we have identified important residues controlling the 3CLpro-ligand binding from per-residue based decomposition of the binding free energy. Finally, we have compared ΔGbind of these two inhibitors with the anti-HIV retroviral drugs, such as lopinavir and darunavir. It is observed that α-ketoamide is more potent compared to lopinavir and darunavir. In the case of lopinavir, a decrease in van der Waals interactions is responsible for the lower binding affinity compared to α-ketoamide. On the other hand, in the case of darunavir, a decrease in the favorable intermolecular electrostatic and van der Waals interactions contributes to lower affinity compared to α-ketoamide. Our study might help in designing rational anti-coronaviral drugs targeting the SARS-CoV-2 main protease. Communicated by Ramaswamy H. Sarma

Binding of anti-apoptotic Bcl-2 with different BH3 peptides: A molecular dynamics study