Investigation on the mechanism for the binding and drug resistance of wild type and mutations of G86 residue in HIV-1 protease complexed with Darunavir by molecular dynamic simulation and free energy calculation

Interaction of rice glutelin with soybean 7S globulin formed co-assemblies with improved functional properties

Rice protein (RP) molecules tend to aggregate, limiting their functional properties and posing a significant challenge to the intensive processing of rice products. This study explored the interaction between soybean 7S globulin (7S) and rice glutelin (RG) to improve the structure and properties of RP. The results show that the tertiary structure of the RG-7S co-assemblies undergoes a certain degree of extension, increases the α-helix content, and reduces the β-sheet content in the secondary structure. Molecular dynamics provide further verification that hydrogen bonding and hydrophobic interactions are the main drivers of these conformational changes. When the RG-7S ratio was 1:1.500 (g/g), the solubility increased 20-fold, and the emulsifying activity index and foaming capacity increased by 3.33 and 1.89 times, respectively. This study confirms that 7S co-assembly with RG enhances the functional properties of RG, demonstrating the potential of this strategy for application in the food industry.

Molecular-Level Strategy from Bottom-Up to Acquire High-Efficiency Antimicrobial Peptides

The escalating threat of multidrug-resistant pathogens necessitates efficient antimicrobial development, in which antimicrobial peptides (AMPs) have been extensively studied due to their broad-spectrum antibiotic activity. Combining sum frequency generation and molecular dynamics simulation, we rationally designed AMPs by identifying two structural principles. Simply increasing the number of basic amino acids does not reliably improve antimicrobial efficacy; terminal (N' and C') phenylalanine residues can enhance AMP's membrane interfacial activity via the hydrophobic effect. After two rounds of sequence optimization, among the derivatives we designed, an artificial AMP named GF demonstrated superior membrane binding (especially insertion) and spatial conformation stability. In vitro and in vivo evaluations revealed GF's potential broad-spectrum efficacy against common bacteria and drug-resistant bacteria. Notably, GF exhibited enhanced antimicrobial potency over conventional antibiotics at lower concentrations. Our study established a bottom-up (mechanism-driven) design framework and provided a template for developing precision antimicrobials against resistant infections.

Isolation, structural characterization of natural chondroitin sulfate oligosaccharides and their binding study with anti-angiogenic factors

Drugs that inhibit tumor angiogenesis, promote vascular normalization and improve the tumor microenvironment. However, their application is limited by adaptive or compensatory resistance. Chondroitin sulfate (CS) regulates numerous proteins including pro-angiogenic growth factors, for whom binding affinity depends on sulfation of CS. In this study, we aimed to determine how sulfation of natural tetrasaccharides and hexasaccharides of CS affected binding to the vascular endothelial growth factor (VEGF-A) and fibroblast growth factor 2 (FGF-2). Twenty-eight CS oligosaccharide isomers were obtained by preparative HPLC, tagged with the AEAB fluorescent linker, and identified using an improved chemical derivatization strategy combined with tandem mass spectrometry. CS oligosaccharide microarrays revealed that VEGF-A and FGF-2 bound preferentially to highly sulfated CS, and the GalNAc(4S)GlcA(2S)GalNAc(6S) sequence was found to be indispensable for binding to these proteins. By integrating glycan microarrays with computational modeling, this study revealed the relationship between the structure of CS and its interactions with pro-angiogenic factors. The degree and the specific sulfation patterns on CS should be taken into account when designing anti-angiogenic drugs.

Identification of FDA Approved Drugs with Antiviral Activity against SARS-CoV-2: A Tale from structure-based drug repurposing to host-cell mechanistic investigation

The continuing heavy toll of the COVID-19 pandemic necessitates development of therapeutic options. We adopted structure-based drug repurposing to screen FDA-approved drugs for inhibitory effects against main protease enzyme (Mpro) substrate-binding pocket of SARS-CoV-2 for non-covalent and covalent binding. Top candidates were screened against infectious SARS-CoV-2 in a cell-based viral replication assay. Promising candidates included atovaquone, mebendazole, ouabain, dronedarone, and entacapone, although atovaquone and mebendazole were the only two candidates with IC50s that fall within their therapeutic plasma concentration. Additionally, we performed Mpro assays on the top hits, which demonstrated inhibition of Mpro by dronedarone (IC50 18 µM), mebendazole (IC50 19 µM) and entacapone (IC50 9 µM). Atovaquone showed only modest Mpro inhibition, and thus we explored other potential mechanisms. Although atovaquone is Dihydroorotate dehydrogenase (DHODH) inhibitor, we did not observe inhibition of DHODH at the respective SARS-CoV-2 IC50. Metabolomic profiling of atovaquone treated cells showed dysregulation of purine metabolism pathway metabolite, showing that ecto-5′-nucleotidase (NT5E) is downregulated by atovaquone at concentrations equivalent to its antiviral IC50. Atovaquone and mebendazole are promising candidates targeting SARS-CoV-2, however atovaquone did not significantly inhibit Mpro at therapeutically meaningful concentrations but may inhibit SARS-CoV-2 viral replication by targeting host purine metabolism.

In Silico Study of the Mechanisms Underlying the Action of the Snake Natriuretic-Like Peptide Lebetin 2 during Cardiac Ischemia

Lebetin 2 (L2), a natriuretic-like peptide (NP), exerts potent cardioprotection in myocardial infarction (MI), with stronger effects than B-type natriuretic peptide (BNP). To determine the molecular mechanisms underlying its cardioprotection effect, we used molecular modeling, molecular docking and molecular dynamics (MD) simulation to describe the binding mode, key interaction residues as well as mechanistic insights into L2 interaction with NP receptors (NPRs). L2 binding affinity was determined for human, rat, mouse and chicken NPRs, and the stability of receptor-ligand complexes ascertained during 100 ns-long MD simulations. We found that L2 exhibited higher affinity for all human NPRs compared to BNP, with a rank preference for NPR-A > NPR-C > NPR-B. Moreover, L2 affinity for human NPR-A and NPR-C was higher in other species. Both docking and MD studies revealed that the NPR-C-L2 interaction was stronger in all species compared to BNP. Due to its higher affinity to human receptors, L2 could be used as a therapeutic approach in MI patients. Moreover, the stronger interaction of L2 with NPR-C could highlight a new L2 signaling pathway that would explain its additional effects during cardiac ischemia. Thus, L2 is a promising candidate for drug design toward novel compounds with high potency, affinity and stability.

PLAS-5k: Dataset of Protein-Ligand Affinities from Molecular Dynamics for Machine Learning Applications

Computational methods and recently modern machine learning methods have played a key role in structure-based drug design. Though several benchmarking datasets are available for machine learning applications in virtual screening, accurate prediction of binding affinity for a protein-ligand complex remains a major challenge. New datasets that allow for the development of models for predicting binding affinities better than the state-of-the-art scoring functions are important. For the first time, we have developed a dataset, PLAS-5k comprised of 5000 protein-ligand complexes chosen from PDB database. The dataset consists of binding affinities along with energy components like electrostatic, van der Waals, polar and non-polar solvation energy calculated from molecular dynamics simulations using MMPBSA (Molecular Mechanics Poisson-Boltzmann Surface Area) method. The calculated binding affinities outperformed docking scores and showed a good correlation with the available experimental values. The availability of energy components may enable optimization of desired components during machine learning-based drug design. Further, OnionNet model has been retrained on PLAS-5k dataset and is provided as a baseline for the prediction of binding affinities.

Identification of SARS-CoV-2 surface therapeutic targets and drugs using molecular modeling methods for inhibition of the virus entry

Although COVID-19 emerged as a major concern to public health around the world, no licensed medication has been found as of yet to efficiently stop the virus spread and treat the infection. The SARS-CoV-2 entry into the host cell is driven by the direct interaction of the S1 domain with the ACE-2 receptor followed by conformational changes in the S2 domain, as a result of which fusion peptide is inserted into the target cell membrane, and the fusion process is mediated by the specific interactions between the heptad repeats 1 and 2 (HR1 and HR2) that form the six-helical bundle. Since blocking this interaction between HRs stops virus fusion and prevents its subsequent replication, the HRs inhibitors can be used as anti-COVID drugs. The initial drug selection is based on existing molecular databases to screen for molecules that may have a therapeutic effect on coronavirus. Based on these premises, we chose two approved drugs to investigate their interactions with the HRs (based on docking methods). To this end, molecular dynamics simulations and molecular docking were carried out to investigate the changes in the structure of the SARS-CoV-2 spike protein. Our results revealed, cefpiramide has the highest affinity to S protein, thereby revealing its potential to become an anti-COVID-19 clinical medicine. Therefore, this study offers new ways to re-use existing drugs to combat SARS-CoV-2 infection.

Role of ACE2 receptor and the landscape of treatment options from convalescent plasma therapy to the drug repurposing in COVID-19

Since the first case reports in Wuhan, China, the SARS-CoV-2 has caused a pandemic and took lives of > 8,35,000 people globally. This single-stranded RNA virus uses Angiotensin-converting enzyme 2 (ACE2) as a receptor for entry into the host cell. Overexpression of ACE2 is mainly observed in hypertensive, diabetic and heart patients that make them prone to SARS-CoV-2 infection. Mitigations strategies were opted globally by the governments to minimize transmission of SARS-CoV-2 via the implementation of social distancing norms, wearing the facemasks, and spreading awareness using digital platforms. The lack of an approved drug treatment regimen, and non-availability of a vaccine, collectively posed a challenge for mankind to fight against the SARS-CoV-2 pandemic. In this scenario, repurposing of existing drugs and old treatment options like convalescent plasma therapy can be one of the potential alternatives to treat the disease. The drug repurposing provides a selection of drugs based on the scientific rationale and with a shorter cycle of clinical trials, while plasma isolated from COVID-19 recovered patients can be a good source of neutralizing antibody to provide passive immunity. In this review, we provide in-depth analysis on these two approaches currently opted all around the world to treat COVID-19 patients. For this, we used "Boolean Operators" such as AND, OR & NOT to search relevant research articles/reviews from the PUBMED for the repurposed drugs and the convalescent plasma in the COVID-19 treatment. The repurposed drugs like Chloroquine and Hydroxychloroquine, Tenofovir, Remdesivir, Ribavirin, Darunavir, Oseltamivir, Arbidol (Umifenovir), Favipiravir, Anakinra, and Baricitinib are already being used in clinical trials to treat the COVID-19 patients. These drugs have been approved for a different indication and belong to a diverse category such as anti-malarial/anti-parasitic, anti-retroviral/anti-viral, anti-cancer, or against rheumatoid arthritis. Although, the vaccine would be an ideal option for providing active immunity against the SARS-CoV-2, but considering the current situation, drug repurposing and convalescent plasma therapy and repurposed drugs are the most viable option against SARS-CoV-2.

Role of ACE2 receptor and the landscape of treatment options from convalescent plasma therapy to the drug repurposing in COVID-19

Since the first case reports in Wuhan, China, the SARS-CoV-2 has caused a pandemic and took lives of > 8,35,000 people globally. This single-stranded RNA virus uses Angiotensin-converting enzyme 2 (ACE2) as a receptor for entry into the host cell. Overexpression of ACE2 is mainly observed in hypertensive, diabetic and heart patients that make them prone to SARS-CoV-2 infection. Mitigations strategies were opted globally by the governments to minimize transmission of SARS-CoV-2 via the implementation of social distancing norms, wearing the facemasks, and spreading awareness using digital platforms. The lack of an approved drug treatment regimen, and non-availability of a vaccine, collectively posed a challenge for mankind to fight against the SARS-CoV-2 pandemic. In this scenario, repurposing of existing drugs and old treatment options like convalescent plasma therapy can be one of the potential alternatives to treat the disease. The drug repurposing provides a selection of drugs based on the scientific rationale and with a shorter cycle of clinical trials, while plasma isolated from COVID-19 recovered patients can be a good source of neutralizing antibody to provide passive immunity. In this review, we provide in-depth analysis on these two approaches currently opted all around the world to treat COVID-19 patients. For this, we used "Boolean Operators" such as AND, OR & NOT to search relevant research articles/reviews from the PUBMED for the repurposed drugs and the convalescent plasma in the COVID-19 treatment. The repurposed drugs like Chloroquine and Hydroxychloroquine, Tenofovir, Remdesivir, Ribavirin, Darunavir, Oseltamivir, Arbidol (Umifenovir), Favipiravir, Anakinra, and Baricitinib are already being used in clinical trials to treat the COVID-19 patients. These drugs have been approved for a different indication and belong to a diverse category such as anti-malarial/anti-parasitic, anti-retroviral/anti-viral, anti-cancer, or against rheumatoid arthritis. Although, the vaccine would be an ideal option for providing active immunity against the SARS-CoV-2, but considering the current situation, drug repurposing and convalescent plasma therapy and repurposed drugs are the most viable option against SARS-CoV-2.

Revealing binding selectivity of ligands toward murine double minute 2 and murine double minute X based on molecular dynamics simulations and binding free energy calculations

It is well known that the interactions of p53 with murine double minute 2 and murine double minute X, namely MDM2 and MDMX, have been significant targets of efficient anti-cancer drug design. In this study, molecular dynamics (MD) simulations, principal component (PC) analysis and binding free energy calculations are combined to recognize binding selectivity of three ligands to MDM2 and MDMX. The binding free energies were estimated by using molecular mechanics generalized Born surface area (MM-GBSA) method and the obtained results display that the increase in the binding enthalpy of three ligands to MDM2 relative to MDMX mainly drives the binding selectivity of them toward MDM2 and MDMX. The information obtained from PC analysis shows that the associations of ligands exert important impacts on internal dynamics of MDM2 and MDMX. Meanwhile, the calculations of residue-based free energy decomposition not only identify the hot interaction spots of ligands with MDM2 and MDMX, but also show the residues (L54, M53), (Y67, Y66), (V93, V92), (H96, P95), (I99, I98) and (Y100, Y99) in (MDM2, MDMX) are responsible for most contributions to the binding selectivity of three ligands toward MDM2 and MDMX. It is believed that this work can provide useful information for design of highly selective and dual inhibitors targeting MDM2 and MDMX.Communicated by Ramaswamy H. Sarma.

Theoretical Study on Zearalenol Compounds Binding with Wild Type Zearalenone Hydrolase and V153H Mutant

Zearalenone hydrolase (ZHD) is the only reported α/β-hydrolase that can detoxify zearalenone (ZEN). ZHD has demonstrated its potential as a treatment for ZEN contamination that will not result in damage to cereal crops. Recent researches have shown that the V153H mutant ZHD increased the specific activity against α-ZOL, but decreased its specific activity to β-ZOL. To understand whyV153H mutation showed catalytic specificity for α-ZOL, four molecular dynamics simulations combining with protein network analysis for wild type ZHD α-ZOL, ZHD β-ZOL, V153H α-ZOL, and V153H β-ZOL complexes were performed using Gromacs software. Our theoretical results indicated that the V153H mutant could cause a conformational switch at the cap domain (residues Gly161–Thr190) to affect the relative position catalytic residue (H242). Protein network analysis illustrated that the V153H mutation enhanced the communication with the whole protein and residues with high betweenness in the four complexes, which were primarily assembled in the cap domain and residues Met241 to Tyr245 regions. In addition, the existence of α-ZOL binding to V153H mutation enlarged the distance from the O_AE atom in α-ZOL to the NE2 atom in His242, which prompted the side chain of H242 to the position with catalytic activity, thereby increasing the activity of V153H on the α-ZOL. Furthermore, α-ZOL could easily form a right attack angle and attack distance in the ZHD and α-ZOL complex to guarantee catalytic reaction. The alanine scanning results indicated that modifications of the residues in the cap domain produced significant changes in the binding affinity for α-ZOL and β-ZOL. Our results may provide useful theoretical evidence for the mechanism underlying the catalytic specificity of ZHD.

Recent Developments and Applications of the MMPBSA Method

The Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) approach has been widely applied as an efficient and reliable free energy simulation method to model molecular recognition, such as for protein-ligand binding interactions. In this review, we focus on recent developments and applications of the MMPBSA method. The methodology review covers solvation terms, the entropy term, extensions to membrane proteins and high-speed screening, and new automation toolkits. Recent applications in various important biomedical and chemical fields are also reviewed. We conclude with a few future directions aimed at making MMPBSA a more robust and efficient method.

Effect of polarization on HIV-1protease and fluoro-substituted inhibitors binding energies by large scale molecular dynamics simulations

Molecular dynamics simulations in explicit water are carried out to study the binding of six inhibitors to HIV-1 protease (PR) for up to 700 ns using the standard AMBER force field and polarized protein-specific charge (PPC). PPC is derived from quantum mechanical calculation for protein in solution and therefore it includes electronic polarization effect. Our results show that in all six systems, the bridging water W301 drifts away from the binding pocket in AMBER simulation. However, it is very stable in all six complexes systems using PPC. Especially, intra-protease, protease-inhibitor hydrogen bonds are dynamic stabilized in MD simulation. The computed binding free energies of six complexes have a significantly linear correlation with those experiment values and the correlation coefficient is found to be 0.91 in PPC simulation. However, the result from AMBER simulation shows a weaker correlation with the correlation coefficient of −0.51 due to the lack of polarization effect. Detailed binding interactions of W301, inhibitors with PR are further analyzed and discussed. The present study provides important information to quantitative understanding the interaction mechanism of PR-inhibitor and PR-W301 and these data also emphasizes the importance of both the electronic polarization and the bridging water molecule in predicting precisely binding affinities.

Effects of drug-resistant mutations on the dynamic properties of HIV-1 protease and inhibition by Amprenavir and Darunavir

Molecular dynamics simulations are performed to investigate the dynamic properties of wild-type HIV-1 protease and its two multi-drug-resistant variants (Flap + (L10I/G48V/I54V/V82A) and Act (V82T/I84V)) as well as their binding with APV and DRV inhibitors. The hydrophobic interactions between flap and 80 s (80’s) loop residues (mainly I50-I84’ and I50’-I84) play an important role in maintaining the closed conformation of HIV-1 protease. The double mutation in Act variant weakens the hydrophobic interactions, leading to the transition from closed to semi-open conformation of apo Act. APV or DRV binds with HIV-1 protease via both hydrophobic and hydrogen bonding interactions. The hydrophobic interactions from the inhibitor is aimed to the residues of I50 (I50’), I84 (I84’), and V82 (V82’) which create hydrophobic core clusters to further stabilize the closed conformation of flaps, and the hydrogen bonding interactions are mainly focused with the active site of HIV-1 protease. The combined change in the two kinds of protease-inhibitor interactions is correlated with the observed resistance mutations. The present study sheds light on the microscopic mechanism underlying the mutation effects on the dynamics of HIV-1 protease and the inhibition by APV and DRV, providing useful information to the design of more potent and effective HIV-1 protease inhibitors.

Investigate the binding of catechins to trypsin using docking and molecular dynamics simulation

To explore the inhibitory mechanism of catechins for digestive enzymes, we investigated the binding mode of catechins to a typical digestive enzyme-trypsin and analyzed the structure-activity relationship of catechins, using an integration of molecular docking, molecular dynamics simulation and binding free energy calculation. We found that catechins with different structures bound to a conservative pocket S1 of trypsin, which is comprised of residues 189–195, 214–220 and 225–228. In the trypsin-catechin complexes, Asp189 by forming strong hydrogen bonding, and Gln192, Trp215 and Gly216 through hydrophobic interactions, all significantly contribute to the binding of catechins. The number and the position of hydroxyl and aromatic groups, the structure of stereoisomers, and the orientation of catechins in the binding pocket S1 of trypsin all affect the binding affinity. The binding affinity is in the order of Epigallocatechin gallate (EGCG) > Epicatechin gallate (ECG) > Epicatechin (EC) > Epigallocatechin (EGC), and 2R-3R EGCG shows the strongest binding affinity out of other stereoisomers. Meanwhile, the synergic conformational changes of residues and catechins were also analyzed. These findings will be helpful in understanding the knowledge of interactions between catechins and trypsin and referable for the design of novel polyphenol based functional food and nutriceutical formulas.