Crystal molecular dynamics simulations to speed up MM/PB(GB)SA evaluation of binding free energies of di-mannose deoxy analogs with P51G-m4-Cyanovirin-N

An insight into the binding mechanism of Viprinin and its morpholine and piperidine derivatives with HIV-1 Vpr: molecular dynamics simulation, principal component analysis and binding free energy calculation study

HIV-1 Vpr is an accessory protein responsible for a plethora of functions inside the host cell to promote viral pathogenesis. One of the major functions of Vpr is the G₂ cell cycle arrest. Among several small molecule inhibitors, Viprinin, a coumarin derivative, has been shown to specifically inhibit the G₂ cell cycle arrest activity of Vpr thus making it an excellent choice for a lead molecule to design antiretroviral drug. But the exact mechanism of binding of the Viprinin and its two potent derivatives with Vpr is still not understood. In this study with combined molecular docking, molecular dynamics simulation, Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) method, Principal component analysis and Umbrella sampling simulation, we have explored the binding mechanism of Viprinin and its two derivatives with Vpr. MM-PBSA and Umbrella sampling calculations suggest that Viprinin and ViprininD1 have higher binding energy than ViprininD2. Molecular dynamics simulation shows that the ligands are not very stable inside the initial binding pocket and various hydrophobic interactions are responsible to hold the ligands with Vpr. Vpr backbone Principle Component Analysis (PCA) shows various unique essential motions of Vpr bound with Viprinin and its two derivatives. This study may give detailed insight of the mode of binding of the specified compounds at atomic scale and provide valuable information about the possibility of using these compounds as a potent Vpr inhibitor. Communicated by Ramaswamy H. Sarma.

In silico study on miRNA regulation and NSs protein interactome characterization of the SFTS virus

Severe fever with thrombocytopenia syndrome causing virus i.e. SFTS virus has increased in the last few years. The underlying cause and mechanism of disease progression and development of symptoms is not well known. Many viruses including Hepatitis B, Hepatitis C, HIV-1, Herpes virus, Dengue virus and many others have been seen to regulate their functions at the miRNA level. This study aimed to find out those cellular miRNAs, which can be mimicked or antagonized by the viral genome and analyze the effect of these miRNAs on various gene functions. Investigations in this study suggest a correlation between miRNA regulation with the disease symptoms and progression. By exhaustive literature survey we have tried to identify the interacting partners of the Non Structural S (NSs) protein and characterized the protein-protein interactions. The binding interface that can serve as target for therapeutic studies involving the interfacial residues was analyzed. This study would serve as an avenue to design therapeutics making use of not only protein-protein interactions but also miRNA based regulation as well.

Luteolin Binds Streptolysin O Toxin and Inhibits Its Hemolytic Effects and Cytotoxicity

Group A streptococcus (GAS, Streptococcus pyogenes) is a common pathogen that can cause a variety of human diseases. Streptolysin O (SLO) is an exotoxin produced by GAS. It is a pore-forming toxin (PFT) that exhibits high in vivo toxicity. SLO enables GAS to evade phagocytosis and clearance by neutrophils, induces eukaryotic cell lysis, and activates inflammatory bodies. Luteolin is a natural compound that is produced by a wide range of plant species, and recent studies have shown that luteolin can inhibit the growth and alter the morphological of GAS. Here, we reported that luteolin can weaken the cytotoxicity and hemolytic activity of SLO in vitro. Briefly, luteolin bound SLO with high affinity, inhibited its dissolution of erythrocytes, affected its conformational stability and inhibited the formation of oligomers. To further verify the protective effect of luteolin, we used an in vitro SLO-induced human laryngeal carcinoma epithelial type-2 cells (HEp-2) model. Notably, our results showed luteolin protected HEp-2 cells from SLO induced cytotoxicity and changed in cell membrane permeability. In addition, we explored the role of luteolin in protecting mice from GAS-mediated injury using an aerosolized lung delivery model, and our results indicate that luteolin increases murine survival rate following inoculation with a lethal dose of GAS, and that survival was also associated with decreased pathological damage to lung tissue. Our results suggest that luteolin may be a novel drug candidate for the treatment of GAS infection.

Integrated Approach to Interaction Studies of Pyrene Derivatives with Bovine Serum Albumin: Insights from Theory and Experiment

This work aimed to investigate the interaction of bovine serum albumin with newly synthesized potent new pyrene derivatives (PS1 and PS2), which might prove useful to have a better antibacterial character as found for similar compounds in the previous report [Low et al. Bioconjugate Chemistry2014, 12, 2269−2284]. However, to date, binding studies with plasma protein are still unknown. Steady-state fluorescence spectroscopy and lifetime fluorescence studies show that the static interaction binding mode and binding constants of PS1 and PS2 are 7.39 and 7.81 [K_b × 10⁵ (M^–1)], respectively. The experimental results suggest that hydrophobic forces play a crucial role in interacting pyrene derivatives with BSA protein. To verify this, molecular docking and molecular dynamics simulations were performed to predict the nature of the interaction and the dynamic behavior of the two compounds in the BSA complex, PS1 and PS2, under physiological conditions of pH = 7.1. In addition, the free energies of binding for the BSA-PS1 and BSA-PS2 complexes were estimated at 300 K based on the molecular mechanics of the Poisson–Boltzmann surface (MMPBSA) with the Gromacs package. PS2 was found to have a higher binding affinity than PS1. To determine the behavior of the orbital transitions in the ground state geometry, we found that both compounds have similar orbital transitions from HOMO–LUMO via π → π* and HOMO–1–LUMO+1 via n → π*, which was included in the FMO analysis. A cytotoxicity study was performed to determine the toxicity of the compounds. Based on the MD study, the stability of the compounds with BSA and the dynamic binding modes were further revealed, as well as the nature of the binding force components involved and the important residues involved in the binding process. From the binding energy analysis, it can be assumed that PS2 may be more active than PS1.

Designing of nanobodies against Dengue virus Capsid: a computational affinity maturation approach

Dengue virus, an arbovirus, is one of the most prevalent diseases in the tropical environment and leads to huge number of casualties every year. No therapeutics are available till date against the viral disease and the only medications provide symptomatic relief. In this study, we have focused on utilizing conventional nanobodies and repurposing them for Dengue. Computationally affinity matured, best binding nanobodies tagged with constant antibody regions, could be proposed as therapeutics. These could also be applied for drug delivery purposes due to their high specificity against the viral Capsid. Another application of these nanobodies has been thought to utilize them for diagnostic purposes, to use the nanobodies for viral detection from patient samples at the earliest stage using ELISA. This study may open a new avenue for immunologic study in foreseeable future with the usage of the same molecules for multiple purposes. HighlightsNatural nanobodies against viruses were modified for use against Dengue virus Capsid conserved regions.Computational affinity maturation was performed making use of change in binding affinities upon mutating various residues in the complementary determining regions.Docking studies performed to inspect the docking groove, interface analysis and energy calculations.MM/GBSA calculations done to calculate binding free energy of the complex to determine stability of the complex.Communicated by Ramaswamy H. Sarma.

Exploring the toxic effects and mechanism of methoxylated polybrominated diphenyl ethers (MeO-PBDEs) on thyroxine-binding globulin (TBG): Synergy between spectroscopic and computations

The interaction mechanism between thyroxine-binding globulin (TBG) and three methoxylated polybrominated diphenyl ethers (MeO-PBDEs) was analyzed by steady-state fluorescence, ultraviolet-visible (UV-visible) spectroscopy, circular dichroism (CD), molecular docking and molecular dynamics simulation methods. The results of the molecular docking technique revealed that 2'-MeO-BDE-3, 5-MeO-BDE-47, and 3-MeO-BDE-100 combined with TBG at the active site. The steady-state fluorescence spectra displayed that MeO-PBDEs quenched the endogenous fluorescence of TBG through static quenching mechanism, and complex formation between MeO-PBDEs and TBG was further indicated by UV-vis spectroscopy. The thermodynamic quantities showed that the binding process is spontaneous, and the major forces responsible for the binding are hydrogen bonding and hydrophobic interactions, which are consistent with the results of molecular docking to a certain extent. The results of CD confirmed that the secondary structure of TBG was changed after combining with MeO-PBDEs. The dynamic simulation results illustrated that the protein structure is more compact and changes in the secondary structure of TBG after binding to MeO-PBDEs. Additionally, we also utilized the molecular mechanics/Poisson-Boltzmann surface area (MM-PBSA) method to analyze the binding free energy of TBG and MeO-PBDEs. The results suggest that van der Waals force plays an essential role in the combination.

Identification and characterization of novel RdRp and Nsp15 inhibitors for SARS-COV2 using computational approach

The World Health Organization has declared COVID-19 as a global health emergency. COVID-19 is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and highlights an urgent need for therapeutics. Here, we have employed a series of computer-aided drug repurposing campaign to discover inhibitors of RNA dependent RNA polymerase (RdRp) and Nsp15/EndoU. Subsequently, MD simulation has been performed to observe dynamic behavior of identified leads at the active site of RdRp and Nsp15. We successfully identified novel lead molecule such as Alectinib for RdRp while Naldemedine and Ergotamine for NSP15. These lead molecules were accommodated in the active site of the enzyme and stabilized by the networks of the hydrogen bond, pi type and hydrophobic interaction with key residues of either target. Interestingly, identified compounds show molecular mimicry in terms of molecular interactions with key residues of RdRp and Nsp15 essential for catalysis and substrate interaction. Previously, Alectinib, Naldemedine and Ergotamine were used as drug in different diseases might be repurposed against selected protein targets of COVID19. Finally, we propose that the identified inhibitors represent a novel lead molecule to design a more effective inhibitor to stop the progress of pathogen.

Communicated by Ramaswamy H. Sarma

Computational Simulations Identified Two Candidate Inhibitors of Cdk5/p25 to Abrogate Tau-associated Neurological Disorders

Deregulation of Cdk5 is a hallmark in neurodegenerative diseases and its complex with p25 forms Cdk5/p25, thereby causes severe neuropathological insults. Cdk5/p25 abnormally phosphorylates tau protein, and induces tau-associated neurofibrillary tangles in neurological disorders. Therefore, the pharmacological inhibition of Cdk5/p25 alleviates tau-associated neurological disorders. Herein, computational simulations probed two candidate inhibitors of Cdk5/p25. Structure-based pharmacophore investigated the essential complementary chemical features of ATP-binding site of Cdk5 in complex with roscovitine. Resultant pharmacophore harbored polar interactions with Cys83 and Asp86 residues and non-polar interactions with Ile10, Phe80, and Lys133 residues of Cdk5. The chemical space of selected pharmacophore was comprised of two hydrogen bond donors, one hydrogen bond acceptor, and three hydrophobic features. Decoy test validation of pharmacophore obtained highest Guner-Henry score (0.88) and enrichment factor score (7.23). The screening of natural product drug-like databases by validated pharmacophore retrieved 1126 compounds as candidate inhibitors of Cdk5/p25. The docking of candidate inhibitors filtered 10 molecules with docking score >80.00 and established polar and non-polar interactions with the ATP-binding site residues of Cdk5/p25. Finally, molecular dynamics simulation and binding free energy analyses identified two candidate inhibitors of Cdk5/p25. During 30 ns simulation, the candidate inhibitors established <3.0 Å root mean square deviation and stable hydrogen bond interactions with the ATP-binding site residues of Cdk5/p25. The final candidate inhibitors obtained lowest binding free energies of -122.18 kJ/mol and - 117.26 kJ/mol with Cdk5/p25. Overall, we recommend two natural product candidate inhibitors to target the pharmacological inhibition of Cdk5/p25 in tau-associated neurological disorders.

Investigation of novel chemical scaffolds targeting prolyl oligopeptidase for neurological therapeutics

Prolyl oligopeptidase (POP) is a potential therapeutic target for treatment of several neurological disorders and α-synucleinopathies including Parkinson's disease. Most of the known POP inhibitors failed in the clinical trials due to poor pharmacokinetic properties and blood-brain impermeability. Therefore, a training set of 30 structurally diverse compounds with a wide range of inhibitory activity against POP was used to generate a quantitative pharmacophore model, Hypo 3, to identify potential POP inhibitors with desirable drug-like properties. Validations through test set, cost analysis, and Fisher's randomization methods proved that Hypo 3 accurately predicted the known inhibitors among inactive compounds. Hypo 3 was employed as 3D query for virtual screening on an in-house drug-like chemical database containing compounds with good brain permeability and ADMET parameters. Database screening with Hypo 3 resulted in 99 compounds that were narrowed down to 21 compounds through molecular docking. Among them, five compounds were identified in our earlier studies, while two compounds showed in vitro POP inhibition. The current study proposed new 16 virtually screened compounds as potential inhibitors against POP that possess Gold docking score in the range of 64.61-75.74 and Chemscore of -32.25 to -38.35. Furthermore, the top scoring four hit compounds were subjected to molecular dynamics simulations to reveal their appropriate binding modes and assessing binding free energies. The hit compounds interacted with POP effectively via hydrogen bonds with important active site residues along with hydrophobic interactions. Moreover, the hit compounds had key inter-molecular interactions and better binding free energies as compared to the reference inhibitor. A potential new hydrogen bond interaction was discovered between Hit 2 with the Arg252 residue of POP. To conclude, we propose four hit compounds with new structural scaffolds against POP for the lead development of POP-based therapeutics for neurological disorders.

Effects of Tyr555 and Trp678 on the processivity of cellobiohydrolase A from Ruminiclostridium thermocellum: A simulation study

Cellobiohydrolase A from Ruminiclostridium thermocellum (Cbh9A) is a processive exoglucanase from family 9 and is an important cellobiohydrolase that hydrolyzes cello-oligosaccharide into cellobiose. Residues Tyr555 and Trp678 considerably affect catalytic activity, but their mechanisms are still unknown. To investigate how the Tyr555 and Trp678 affect the processivity of Cbh9A, conventional molecular dynamics, steered molecular dynamics, and free energy calculation were performed to simulate the processive process of wild type (WT)-Cbh9A, Y555S mutant, and W678G mutant. Analysis of simulation results suggests that the binding free energies between the substrate and WT-Cbh9A are lower than those of Y555S and W678G mutants. The pull forces and energy barrier in Y555S and W678G mutants also reduced significantly during the steered molecular dynamics (SMD) simulation compared with that of the WT-Cbh9A. And the potential mean force calculations showed that the pulling energy barrier of Y555S and W678G mutants is much lower than that of WT-Cbh9A.

Accuracy Comparison of Generalized Born Models in the Calculation of Electrostatic Binding Free Energies

The need for accurate yet efficient representation of the aqueous environment in biomolecular modeling has led to the development of a variety of generalized Born (GB) implicit solvent models. While many studies have focused on the accuracy of available GB models in predicting solvation free energies, a systematic assessment of the quality of these models in binding free energy calculations, crucial for rational drug design, has not been undertaken. Here, we evaluate the accuracies of eight common GB flavors (GB-HCT, GB-OBC, GB-neck2, GBNSR6, GBSW, GBMV1, GBMV2, and GBMV3), available in major molecular dynamics packages, in predicting the electrostatic binding free energies ( ΔΔ G_el) for a diverse set of 60 biomolecular complexes belonging to four main classes: protein-protein, protein-drug, RNA-peptide, and small complexes. The GB flavors are examined in terms of their ability to reproduce the results from the Poisson-Boltzmann (PB) model, commonly used as accuracy reference in this context. We show that the agreement with the PB of ΔΔ G_el estimates varies widely between different GB models and also across different types of biomolecular complexes, with R² correlations ranging from 0.3772 to 0.9986. A surface-based "R6" GB model recently implemented in AMBER shows the closest overall agreement with reference PB ( R² = 0.9949, RMSD = 8.75 kcal/mol). The RNA-peptide and protein-drug complex sets appear to be most challenging for all but one model, as indicated by the large deviations from the PB in ΔΔ G_el. Small neutral complexes present the least challenge for most of the GB models tested. The quantitative demonstration of the strengths and weaknesses of the GB models across the diverse complex types provided here can be used as a guide for practical computations and future development efforts.

Exploration of binding and inhibition mechanism of a small molecule inhibitor of influenza virus H1N1 hemagglutinin by molecular dynamics simulation

Influenza viruses are a major public health threat worldwide. The influenza hemagglutinin (HA) plays an essential role in the virus life cycle. Due to the high conservation of the HA stem region, it has become an especially attractive target for inhibitors for therapeutics. In this study, molecular simulation was applied to study the mechanism of a small molecule inhibitor (MBX2329) of influenza HA. Behaviors of the small molecule under neutral and acidic conditions were investigated, and an interesting dynamic binding mechanism was found. The results suggested that the binding of the inhibitor with HA under neutral conditions facilitates only its intake, while it interacts with HA under acidic conditions using a different mechanism at a new binding site. After a series of experiments, we believe that binding of the inhibitor can prevent the release of HA1 from HA2, further maintaining the rigidity of the HA2 loop and stabilizing the distance between the long helix and short helices. The investigated residues in the new binding site show high conservation, implying that the new binding pocket has the potential to be an effective drug target. The results of this study will provide a theoretical basis for the mechanism of new influenza virus inhibitors.

Insights into structure and activity of natural compound inhibitors of pneumolysin

Pneumolysin is the one of the major virulence factor of the bacterium Streptococcus pneumoniae. In previous report, it is shown that β-sitosterol, a natural compound without antimicrobial activity, is a potent antagonist of pneumolysin. Here, two new pneumolysin natural compound inhibitors, with differential activity, were discovered via haemolysis assay. To explore the key factor of the conformation for the inhibition activity, the interactions between five natural compound inhibitors with differential activity and pneumolysin were reported using molecular modelling, the potential of mean force profiles. Interestingly, it is found that incorporation of the single bond (C22-C23-C24-C25) to replace the double bond (hydrocarbon sidechain) improved the anti-haemolytic activity. In view of the molecular modelling, binding of the five inhibitors to the conserved loop region (Val372, Leu460, and Tyr461) of the cholesterol binding sites led to stable complex systems, which was consistent with the result of β-sitosterol. Owing to the single bond (C22-C23-C24-C25), campesterol and brassicasterol could form strong interactions with Val372 and show higher anti-haemolytic activity, which indicated that the single bond (C22-C23-C24-C25) in inhibitors was required for the anti-haemolytic activity. Overall, the current molecular modelling work provides a starting point for the development of rational design and higher activity pneumolysin inhibitors.

Effect of the Crystal Environment on Side-Chain Conformational Dynamics in Cyanovirin-N Investigated through Crystal and Solution Molecular Dynamics Simulations

Side chains in protein crystal structures are essential for understanding biochemical processes such as catalysis and molecular recognition. However, crystal packing could influence side-chain conformation and dynamics, thus complicating functional interpretations of available experimental structures. Here we investigate the effect of crystal packing on side-chain conformational dynamics with crystal and solution molecular dynamics simulations using Cyanovirin-N as a model system. Side-chain ensembles for solvent-exposed residues obtained from simulation largely reflect the conformations observed in the X-ray structure. This agreement is most striking for crystal-contacting residues during crystal simulation. Given the high level of correspondence between our simulations and the X-ray data, we compare side-chain ensembles in solution and crystal simulations. We observe large decreases in conformational entropy in the crystal for several long, polar and contacting residues on the protein surface. Such cases agree well with the average loss in conformational entropy per residue upon protein folding and are accompanied by a change in side-chain conformation. This finding supports the application of surface engineering to facilitate crystallization. Our simulation-based approach demonstrated here with Cyanovirin-N establishes a framework for quantitatively comparing side-chain ensembles in solution and in the crystal across a larger set of proteins to elucidate the effect of the crystal environment on protein conformations.

The opening/closure of the P-loop and hinge of BCR-ABL1 decodes the low/high bioactivities of dasatinib and axitinib

The open/closed conformations reveal the low/high bioactivities of the ligands.

Investigation of an "alternate water supply system" in enzymatic hydrolysis in the processive endocellulase Cel7A from Rasamsonia emersonii by molecular dynamics simulation

Cel7A from Rasamsonia emersonii is one of the processive endocellulases classified under family 7 glycoside hydrolase. Molecular dynamics simulations were carried out to obtain the optimized sliding and hydrolyzing conformations, in which the reducing ends of sugar chains are located on different sites. Hydrogen bonds are investigated to clarify the interactions between protein and substrate in either conformation. Nine hydrogen bonding interactions are identified in the sliding conformation, and six similar interactions are also found correspondingly in the hydrolyzing conformation. In addition, four strong hydrophobic interactions are also determined. The domain cross-correlation map analysis shows movement correlation of protein including autocorrelation between residues. The root mean square fluctuations analysis represents the various flexibilities of different fragment in the two conformations. Comparing the two conformations reveals the water-supply mechanism of selective hydrolysis of cellulose in Cel7A. The mechanism can be described as follow. When the reducing end of substrate slides from the unhydrolyzing site (sliding conformation) to the hydrolyzing site (hydrolyzing conformation), His225 is pushed down and rotated, the rotation leads to the movement of Glu209 with the interstrand hydrogen bonding in β-sheet. It further makes Asp211 close to the hydrolysis center and provides a water molecule bounding on its carboxyl in the previous unhydrolyzing site. After the hydrolysis takes place and the product is excluded from the enzyme, the Asp211 comes back to its initial position. In summary, Asp211 acts as an elevator to transport outer water molecules into the hydrolysis site for every other glycosidic bond.

Structural Basis of Fullerene Derivatives as Novel Potent Inhibitors of Protein Tyrosine Phosphatase 1B: Insight into the Inhibitory Mechanism through Molecular Modeling Studies

Protein tyrosine phosphatase 1B (PTP1B) has become an outstanding target for the treatment of diabetes and obesity. Recent research has demonstrated that some fullerene derivatives serve as a new nanoscale-class of potent inhibitors of PTP1B, but the specific mechanism remains unclear. Several molecular modeling methods (molecular docking, molecular dynamics simulations, and molecular mechanics/generalized Born surface area calculations) were integrated to provide insight into the binding mode and inhibitory mechanism of the new class of fullerene inhibitors. The results reveal that PTP1B with an open WPD loop is more susceptible to the combination with the fullerene inhibitor because of their comparable shapes and sizes. When the WPD loop fluctuates to the open conformation, the inhibitor falls into the active pocket and induces conformational rotation of the WPD loop. This rotation is closely related to the reduction of the catalytic activity of PTP1B. In addition, it is suggested that compound 1, like compound 2, is a competitive inhibitor since it blocks the active site to prevent the binding of the substrate. The high binding affinity of fullerene-based compounds and the transition of the WPD loop, caused by the specific structural property of the hydrophobic fullerene core and the appended polar groups, make these fullerene derivatives efficient competitive inhibitors. The theoretical results provide useful clues for further investigation of the noval inhibitors of PTP1B at the nanoscale.

Molecular interactions of UvrB protein and DNA from Helicobacter pylori: Insight into a molecular modeling approach

Helicobacter pylori (H. pylori) persevere in the human stomach, an environment in which they encounter many DNA-damaging conditions, including gastric acidity. The pathogenicity of H. pylori is enhanced by its well-developed DNA repair mechanism, thought of as 'machinery,' such as nucleotide excision repair (NER). NER involves multi-enzymatic excinuclease proteins (UvrABC endonuclease), which repair damaged DNA in a sequential manner. UvrB is the central component in prokaryotic NER, essential for damage recognition. Therefore, molecular modeling studies of UvrB protein from H. pylori are carried out with homology modeling and molecular dynamics (MD) simulations. The results reveal that the predicted structure is bound to a DNA hairpin with 3-bp stem, an 11-nucleotide loop, and 3-nt 3' overhang. In addition, a mutation of the Y96A variant indicates reduction in the binding affinity for DNA. Free-energy calculations demonstrate the stability of the complex and help identify key residues in various interactions based on residue decomposition analysis. Stability comparative studies between wild type and mutant protein-DNA complexes indicate that the former is relatively more stable than the mutant form. This predicted model could also be useful in designing new inhibitors for UvrB protein, as well as preventing the pathogenesis of H. pylori.

Exploration of Novel Inhibitors for Bruton's Tyrosine Kinase by 3D QSAR Modeling and Molecular Dynamics Simulation

Bruton’s tyrosine kinase (BTK) is a cytoplasmic, non-receptor tyrosine kinase which is expressed in most of the hematopoietic cells and plays an important role in many cellular signaling pathways. B cell malignancies are dependent on BCR signaling, thus making BTK an efficient therapeutic target. Over the last few years, significant efforts have been made in order to develop BTK inhibitors to treat B-cell malignancies, and autoimmunity or allergy/hypersensitivity but limited success has been achieved. Here in this study, 3D QSAR pharmacophore models were generated for Btk based on known IC₅₀ values and experimental energy scores with extensive validations. The five features pharmacophore model, Hypo1, includes one hydrogen bond acceptor lipid, one hydrogen bond donor, and three hydrophobic features, which has the highest correlation coefficient (0.98), cost difference (112.87), and low RMS (1.68). It was further validated by the Fisher’s randomization method and test set. The well validated Hypo1 was used as a 3D query to search novel Btk inhibitors with different chemical scaffold using high throughput virtual screening technique. The screened compounds were further sorted by applying ADMET properties, Lipinski’s rule of five and molecular docking studies to refine the retrieved hits. Furthermore, molecular dynamic simulation was employed to study the stability of docked conformation and to investigate the binding interactions in detail. Several important hydrogen bonds with Btk were revealed, which includes the gatekeeper residues Glu475 and Met 477 at the hinge region. Overall, this study suggests that the proposed hits may be more effective inhibitors for cancer and autoimmune therapy.

A Rigid Hinge Region Is Necessary for High-Affinity Binding of Dimannose to Cyanovirin and Associated Constructs

Mutations in the hinge region of cyanovirin-N (CVN) dictate its preferential oligomerization state. Constructs with the Pro51Gly mutation preferentially exist as monomers, whereas wild-type cyanovirin can form domain-swapped dimers under certain conditions. Because the hinge region is an integral part of the high-affinity binding site of CVN, we investigated whether this mutation affects the shape, flexibility, and binding affinity of domain B for dimannose. Our studies indicate that the capability of monomeric wild-type CVN to resist mechanical perturbations is enhanced when compared to that of constructs in which the hinge region is more flexible. Our computational results also show that enhanced flexibility leads to blocking of the binding site by allowing different rotational isomeric states of Asn53. Moreover, at higher temperatures, this observed flexibility leads to an interaction between Asn53 and Asn42, further hindering access to the binding site. On the basis of these results, we predicted that binding affinity for dimannose would be more favorable for cyanovirin constructs containing a wild-type hinge region, whereas affinity would be impaired in the case of mutants containing Pro51Gly. Experimental characterization by isothermal titration calorimetry of a set of cyanovirin mutants confirms this hypothesis. Those possessing the Pro51Gly mutation are consistently inferior binders.

Novel chemical scaffolds of the tumor marker AKR1B10 inhibitors discovered by 3D QSAR pharmacophore modeling

AIM

Recent evidence suggests that aldo-keto reductase family 1 B10 (AKR1B10) may be a potential diagnostic or prognostic marker of human tumors, and that AKR1B10 inhibitors offer a promising choice for treatment of many types of human cancers. The aim of this study was to identify novel chemical scaffolds of AKR1B10 inhibitors using in silico approaches.

METHODS

The 3D QSAR pharmacophore models were generated using HypoGen. A validated pharmacophore model was selected for virtual screening of 4 chemical databases. The best mapped compounds were assessed for their drug-like properties. The binding orientations of the resulting compounds were predicted by molecular docking. Density functional theory calculations were carried out using B3LYP. The stability of the protein-ligand complexes and the final binding modes of the hit compounds were analyzed using 10 ns molecular dynamics (MD) simulations.

RESULTS

The best pharmacophore model (Hypo 1) showed the highest correlation coefficient (0.979), lowest total cost (102.89) and least RMSD value (0.59). Hypo 1 consisted of one hydrogen-bond acceptor, one hydrogen-bond donor, one ring aromatic and one hydrophobic feature. This model was validated by Fischer's randomization and 40 test set compounds. Virtual screening of chemical databases and the docking studies resulted in 30 representative compounds. Frontier orbital analysis confirmed that only 3 compounds had sufficiently low energy band gaps. MD simulations revealed the binding modes of the 3 hit compounds: all of them showed a large number of hydrogen bonds and hydrophobic interactions with the active site and specificity pocket residues of AKR1B10.

CONCLUSION

Three compounds with new structural scaffolds have been identified, which have stronger binding affinities for AKR1B10 than known inhibitors.

Pentacycloundecane lactam vs lactone norstatine type protease HIV inhibitors: binding energy calculations and DFT study

Background

Novel pentacycloundecane (PCU)-lactone-CO-EAIS peptide inhibitors were designed, synthesized, and evaluated against wild-type C-South African (C-SA) HIV-1 protease. Three compounds are reported herein, two of which displayed IC₅₀ values of less than 1.00 μM. A comparative MM-PB(GB)SA binding free energy of solvation values of PCU-lactam and lactone models and their enantiomers as well as the PCU-lactam-NH-EAIS and lactone-CO-EAIS peptide inhibitors and their corresponding diastereomers complexed with South African HIV protease (C-SA) was performed. This will enable us to rationalize the considerable difference between inhibitory concentration (IC₅₀) of PCU-lactam-NH-EAIS and PCU-lactone-CO-EAIS peptides.

Results

The PCU-lactam model exhibited more negative calculated binding free energies of solvation than the PCU-lactone model. The same trend was observed for the PCU-peptide inhibitors, which correspond to the experimental activities for the PCU-lactam-NH-EAIS peptide (IC₅₀ = 0.076 μM) and the PCU-lactone-CO-EAIS peptide inhibitors (IC₅₀ = 0.850 μM). Furthermore, a density functional theory (DFT) study on the natural atomic charges of the nitrogen and oxygen atoms of the three PCU-lactam, PCU-lactim and PCU-lactone models were performed using natural bond orbital (NBO) analysis. Electrostatic potential maps were also used to visualize the electron density around electron-rich regions. The asymmetry parameter (η) and quadrupole coupling constant (χ) values of the nitrogen and oxygen nuclei of the model compounds were calculated at the same level of theory. Electronic molecular properties including polarizability and electric dipole moments were also calculated and compared. The Gibbs theoretical free solvation energies of solvation (∆G_solv) were also considered.

Conclusions

A general trend is observed that the lactam species appears to have a larger negative charge distribution around the heteroatoms, larger quadrupole constant, dipole moment and better solvation energy, in comparison to the PCU-lactone model. It can be argued that these characteristics will ensure better eletronic interaction between the lactam and the receptor, corresponding to the observed HIV protease activities in terms of experimental IC₅₀ data.

Electronic supplementary material

The online version of this article (doi:10.1186/s12929-015-0115-5) contains supplementary material, which is available to authorized users.

Insight into the interactive residues between two domains of human somatic Angiotensin-converting enzyme and Angiotensin II by MM-PBSA calculation and steered molecular dynamics simulation

Angiotensin-converting enzyme (ACE), a membrane-bound zinc metallopeptidase, catalyzes the formation of Angiotensin-II (AngII) and the deactivation of bradykinin in the renin-angiotensin-aldosterone and kallikrein-kinin systems. As a hydrolysis product of ACE, AngII is regarded as an inhibitor and displays stronger competitive inhibition in the C-domain than the N-domain of ACE. However, the AngII binding differences between the two domains and the mechanisms behind AngII dissociation from the C-domain are rarely explored. In this work, molecular docking, Molecular Mechanics/Poisson-Boltzmann Surface Area calculation, and steered molecular dynamics (SMD) are applied to explore the structures and interactions in the binding or unbinding of AngII with the two domains of human somatic ACE. Calculated free energy values suggest that the C-domain-AngII complex is more stable than the N-domain-AngII complex, consistent with available experimental data. SMD simulation results imply that electrostatic interaction is dominant in the dissociation of AngII from the C-domain. Moreover, Gln106, Asp121, Glu123, and Tyr213 may be the key residues in the unbinding pathway of AngII. The simulation results in our work provide insights into the interactions between the two domains of ACE and its natural peptide inhibitor AngII at a molecular level. Moreover, the results provide theoretical clues for the design of new inhibitors.

A flexible docking scheme efficiently captures the energetics of glycan-cyanovirin binding

Cyanovirin-N (CVN), a cyanobacterial lectin, exemplifies a class of antiviral agents that inhibit HIV by binding to the highly glycosylated envelope protein gp120. Here, we investigate the energetics of glycan recognition using a computationally inexpensive flexible docking approach, backbone perturbation docking (BP-Dock). We benchmarked our method using two mutants of CVN: P51G-m4-CVN, which binds dimannose with high affinity through domain B, and CVN((mutDB)), in which binding to domain B has been abolished through mutation of five polar residues to small nonpolar side chains. We investigated the energetic contribution of these polar residues along with the additional position 53 by docking dimannose to single-point CVN mutant models. Analysis of the docking simulations indicated that the E41A/G and T57A mutations led to a significant decrease in binding energy scores due to rearrangements of the hydrogen-bond network that reverberated throughout the binding cavity. N42A decreased the binding score to a level comparable to that of CVN((mutDB)) by affecting the integrity of the local protein structure. In contrast, N53S resulted in a high binding energy score, similar to P51G-m4-CVN. Experimental characterization of the five mutants by NMR spectroscopy confirmed the binding affinity pattern predicted by BP-Dock. Despite their mostly conserved fold and stability, E41A, E41G, and T57A displayed dissociation constants in the millimolar range. N53S showed a binding constant in the low micromolar range, similar to that observed for P51G-m4-CVN. No binding was observed for N42A. Our results show that BP-Dock is a useful tool for rapidly screening the relative binding affinity pattern of in silico-designed mutants compared with wild-type, supporting its use to design novel mutants with enhanced binding properties.

Practical Aspects of Free-Energy Calculations: A Review

Free-energy calculations in the framework of classical molecular dynamics simulations are nowadays used in a wide range of research areas including solvation thermodynamics, molecular recognition, and protein folding. The basic components of a free-energy calculation, that is, a suitable model Hamiltonian, a sampling protocol, and an estimator for the free energy, are independent of the specific application. However, the attention that one has to pay to these components depends considerably on the specific application. Here, we review six different areas of application and discuss the relative importance of the three main components to provide the reader with an organigram and to make nonexperts aware of the many pitfalls present in free energy calculations.

The antiviral lectin cyanovirin-N: probing multivalency and glycan recognition through experimental and computational approaches

CVN (cyanovirin-N), a small lectin isolated from cyanobacteria, exemplifies a novel class of anti-HIV agents that act by binding to the highly glycosylated envelope protein gp120 (glycoprotein 120), resulting in inhibition of the crucial viral entry step. In the present review, we summarize recent work in our laboratory and others towards determining the crucial role of multivalency in the antiviral activity, and we discuss features that contribute to the high specificity and affinity for the glycan ligand observed in CVN. An integrated approach that encompasses structural determination, mutagenesis analysis and computational work holds particular promise to clarify aspects of the interactions between CVN and glycans.

The role of Glu41 in the binding of dimannose to P51G-m4-CVN

The carbohydrate binding protein, Cyanovirin-N, obtained from cyanobacteria, consists of high-affinity and low-affinity binding domains. To avoid the formation of a domain swapped structure in solution and also to better focus on the binding of carbohydrates at the high-affinity site, the Ghirlanda group (Biochemistry, 46, 2007, 9199-9207) engineered the P51G-m4-CVN mutant which does not dimerize nor binds at the low-affinity site. This mutant provides an excellent starting point for the experimental and computational study of further transformations to enhance binding at the high-affinity site as well as to retool this site for the possible binding of different sugars. However, before such endeavors are pursued, detailed understanding of apparently key interactions both present in wild-type and P51G-m4-CVN at the high-affinity site must be derived and controversies about the importance of certain residues must be resolved. One such interaction is that of Glu41, a charged residue in intimate contact with 2'OH of dimannose at the nonreducing end. We do so computationally by performing two mutations using the thermodynamic integration formalism in explicit solvent. Mutations of P51G-m4-CVN Glu41 to Ala41 and Gly41 reveal that whereas the loss of Coulomb interactions result in a free energy penalty of about 2.1 kcal/mol, this is significantly compensated by favorable contributions to the Lennard-Jones portion of the transformation, resulting in almost no change in the free energy of binding. At least in terms of free energetics, and in the case of this particular CVN mutant, Glu41 does not appear to be as important as previously thought. This is not because of lack of extensive hydrogen bonding with the ligand but instead because of other compensating factors.

Packing interface energetics in different crystal forms of the λ Cro dimer

Variation among crystal structures of the λ Cro dimer highlights conformational flexibility. The structures range from a wild type closed to a mutant fully open conformation, but it is unclear if each represents a stable solution state or if one may be the result of crystal packing. Here we use molecular dynamics (MD) simulation to investigate the energetics of crystal packing interfaces and the influence of site-directed mutagenesis on them in order to examine the effect of crystal packing on wild type and mutant Cro dimer conformation. Replica exchange MD of mutant Cro in solution shows that the observed conformational differences between the wild type and mutant protein are not the direct consequence of mutation. Instead, simulation of Cro in different crystal environments reveals that mutation affects the stability of crystal forms. Molecular Mechanics Poisson-Boltzmann Surface Area binding energy calculations reveal the detailed energetics of packing interfaces. Packing interfaces can have diverse properties in strength, energetic components, and some are stronger than the biological dimer interface. Further analysis shows that mutation can strengthen packing interfaces by as much as ∼5 kcal/mol in either crystal environment. Thus, in the case of Cro, mutation provides an additional energetic contribution during crystal formation that may stabilize a fully open higher energy state. Moreover, the effect of mutation in the lattice can extend to packing interfaces not involving mutation sites. Our results provide insight into possible models for the effect of crystallization on Cro conformational dynamics and emphasize careful consideration of protein crystal structures.

Insight into the dynamic interaction between different flavonoids and bovine serum albumin using molecular dynamics simulations and free energy calculations

In this study, the binding of Bovine serum albumin (BSA) with three flavonoids, kaempferol-3-O-a-L-rhamnopyranosyl-(1-3)-a-L-rhamnopyranosyl-(1-6)-b-D-galacto- pyranoside (drug 1),kaempfol-7-O-rhamnosyl-3-O-rutinoside (drug 2)andkaempferide-7-O-(4"-O-acetylrhamnosyl)-3-O-ruti- noside (drug 3) is investigated by molecular docking, molecular dynamics (MD) simulation, and binding free energy calculation. The free energies are consistent with available experimental results and suggest that the binding site of BSA-drug1 is more stable than those of BSA-drug2 and BSA-drug3. The energy decomposition analysis is performed and reveals that the electrostatic interactions play an important role in the stabilization of the binding site of BSA-drug1 while the van der Waals interactions contribute largely to stabilization of the binding site of BSA-drug2 and BSA-drug3. The key residues stabilizing the binding sites of BSA-drug1, BSA-drug2 and BSA-drug3 are identified based on the residue decomposition analysis.

Free energy calculations of protein-ligand interactions

In the calculation of free energies of binding for protein-ligand complexes, we distinguish endpoint methods, methods involving alchemical modifications and methods that physically displace the ligand from the protein. Most methodological advances seem to come from a clever combination of multiple existing methods to enhance the sampling or to utilize specific advantages of various approaches. The coupling parameters common in thermodynamic integration and in Hamiltonian replica exchange are for instance combined to yield replica exchange thermodynamic integration. As new methods mostly aim to improve efficiency or to attain more complete sampling, there are good prospects to understand and tackle the sampling problem better and to shift the focus towards the scoring problem in the context of more robust and accurate force fields.