Numeric analysis of reversibility of classic movement equations and constructive criteria of estimating quality of molecular dynamic simulations

Interaction of the renin inhibitor aliskiren with the SARS-CoV-2 main protease: a molecular docking study

The renin protein is an upstream enzymatic regulator of the renin-aldosterone-angiotensin system (RAAS) essential for the maintenance of blood pressure. The angiotensin-converting enzyme-2 (ACE2) is a major component of the RAAS and a cell surface receptor exploited by the SARS-CoV-2 virus to enter host cells. A recent molecular modeling study has revealed that the direct renin peptide inhibitor remikiren can bind to the catalytic site of SARS-CoV-2 main protease (M^pro). By analogy, we postulated that the non-peptidic drug aliskiren, a more potent renin inhibitor than remikiren and a drug routinely used to treat hypertension, may also be able to interact with M^pro. An in silico comparison of the binding of the two compounds to M^pro indicates that aliskiren (ΔE = -75.9 kcal/mol) can form stable complexes with the main viral protease, binding to the active site, as remikiren (ΔE = -83.2 kcal/mol). The comparison with a panoply of 30 references compounds (mainly antiviral drugs) indicated that remikiren is a potent M^pro binder comparable to drugs like glecaprevir and pibrentasvir (ΔE = -96.5 kcal/mol). The energy of interaction (ΔE) of aliskiren with M^pro is about 10% lower than with remikiren, comparable to that calculated with drugs like velpatasvir and sofosbuvir. A model is proposed to define the drug binding site, with the best binders (including remikiren) penetrating deeply into the site, whereas the less potent binders (including aliskiren) interact more superficially with the protein.Communicated by Ramaswamy H. Sarma.

A comparative study of different docking methodologies to assess the protein-ligand interaction for the E. coli MurB enzyme

We have investigated the active site of E. coli MurB using the Quantum Mechanics/Molecular Mechanics (QM/MM) methodology. The docking of three novel series of 4-thiazolidinone derivatives has been performed using two methods: rigid docking and flexible docking (Induced Fit Docking: IFD). The results have been compared to understand the conformational aspects of the enzyme. The docking results from rigid docking show that the ligands with highly negative ΔG_bind have poor docking scores. In addition, the value of the regression coefficient (R) obtained on correlating the ΔG_bind and the experimental pMIC values is insignificant. On keeping the protein flexible, there is a remarkable improvement in both the docking score and ΔG_bind, along with a good value of R (0.64). Two important residues, Tyr254 and Try190 are found to be highly displaced during the flexible docking and hence their role in effective ligand binding has been confirmed. Thus, comparing the two methodologies, IFD has emerged as the more appropriate one for studying the E. coli MurB enzyme. To further substantiate the findings, MD studies over a time period of 20 ns have been performed on the IFD-LIII j and Rigid/XP-LIII j complexes and the results shows the former complex to be more stable, with lower average RMSD and higher average ΔG_bind.Communicated by Ramaswamy H. Sarma.

Molecular docking study of GSK-3β interaction with nomilin, kihadanin B, and related limonoids and triterpenes with a furyl-δ-lactone core

Glycogen synthase kinase-3β (GSK-3β) is a target enzyme considered for the treatment of multiple human diseases, from neurodegenerative pathologies to viral infections and cancers. Numerous inhibitors of GSK-3β have been discovered but thus far only a few have reached clinical trials and only one drug, tideglusib (1), has been registered. Natural products targeting GSK-3β have been identified, including the two anticancer limonoids obacunone (5) and gedunin (4), both presenting a furyl-δ-lactone core. To help identifying novel GSK-3β ligands, we have performed a molecular docking study with 15 complementary natural products bearing a furyl-δ-lactone unit (such as limonin (6) and kihadanins A (8) and B (9)) or a closely related structure (such as cedrelone (10) and nimbolide (11)). The formation of GSK-3β-binding complexes for those natural products was compared to reference GSK-3β ATP-competitive inhibitors LY2090314 (3) and AR-A014418 (2). Our in silico analysis led to the identification of two new GSK-3β-binding natural products: kihadanin B (9) and nomilin (7). The latter surpassed the reference compounds in terms of calculated empirical energy of interaction (ΔE). Nomilin (7) can possibly bind to the active site of GSK-3β, notably via the furyl-δ-lactone core and its 1-acetyl group, implicated in the protein interaction. Compound structure-binding relationships are discussed. The study should help the discovery of novel natural products targeting GSK-3β.

Molecular docking study of britannin binding to PD-L1 and related anticancer pseudoguaianolide sesquiterpene lactones

The pseudoguaianolide-type sesquiterpene lactone (SL) britannin (BRT), found in different Inula species, has been characterized as a potent anticancer agent acting via modulation of the transcription factor NFkB and the Nrf2-Keap1 signaling pathway. In addition, a BRT-induced down-regulation of the immune checkpoint PD-L1 (programmed cell death ligand 1) expressed on cancer cells has been evidenced. Here we have performed a docking analysis of the direct binding of BRT to the PD-L1 protein, both in its monomeric and dimeric state. BRT appears to form stable complexes with PD-L1, with a preference for the dimeric form, binding at the interface of the two monomers. The calculated empirical energy of interaction (ΔE) value reaches -63.1 kcal/mol for the BRT-PD-L1 dimer complex, not far from the value calculated with the reference PD-L1 ligand BMS-202 (ΔE = -73.4 kcal/mol) under identical conditions. We also studied the potential PD-L1 dimer binding of 15 pseudoguaianolide sesquiterpene lactones analogues to BRT, including helenalin, gaillardin, bigelovin, coronopilin, and others. The docking analysis predicted that the SL chamissonolide (CHM) can also form equally stable complexes with PD-L1 dimer (ΔE = -64.8 kcal/mol). Preliminary compound structure-PD-L1 binding relationships have been delineated. This computational study supports the proposed interaction of BRT with PD-L1 and provides a guidance to the design of novel PD-L1 binders incorporating a SL-like tricyclic core unit.

Advances in Molecular Dynamics Simulations and Enhanced Sampling Methods for the Study of Protein Systems

Molecular dynamics (MD) simulation is a rigorous theoretical tool that when used efficiently could provide reliable answers to questions pertaining to the structure-function relationship of proteins. Data collated from protein dynamics can be translated into useful statistics that can be exploited to sieve thermodynamics and kinetics crucial for the elucidation of mechanisms responsible for the modulation of biological processes such as protein-ligand binding and protein-protein association. Continuous modernization of simulation tools enables accurate prediction and characterization of the aforementioned mechanisms and these qualities are highly beneficial for the expedition of drug development when effectively applied to structure-based drug design (SBDD). In this review, current all-atom MD simulation methods, with focus on enhanced sampling techniques, utilized to examine protein structure, dynamics, and functions are discussed. This review will pivot around computer calculations of protein-ligand and protein-protein systems with applications to SBDD. In addition, we will also be highlighting limitations faced by current simulation tools as well as the improvements that have been made to ameliorate their efficiency.