Insights into drug resistance of mutations D30N and I50V to HIV-1 protease inhibitor TMC-114: Free energy calculation and molecular dynamic simulation

The single mutations D30N and I50V are considered as the key residue mutations of the HIV-1 protease drug resistance to inhibitors in clinical use. In this work, molecular dynamics (MD) simulations combined with the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method have been performed to investigate the drug-resistant mechanisms of D30N and I50V to an inhibitor TMC-114. The analyses of absolute binding free energies using the separate trajectory approach suggests that the decrease in the van der Waals energy and electrostatic energy in the gas phase results in the drug resistance of D30N to TMC-114, while for I50V, the decrease in the electrostatic energy mainly drive its drug resistance to TMC-114. Detailed binding free energies between TMC-114 and individual protein residues are computed by using a per-residue basis decomposition method, which provides insights into the inhibitor-protein binding mechanism and also explains the drug-resistant mechanisms of mutations D30N and I50V to TMC-114. The study shows that the loss of the hydrogen bond between TMC-114 and the side chain of Asn30′ is the main driving force of the resistance of D30N to TMC-114, and in the case of I50V, the increase in the polar solvation energies between TMC-114 and two residues Val50′ and Asp30′ definitively drives the resistance of I50V to TMC-114. We expect that this work can provide some helpful insights into the nature of mutational effect and aid the future design of better inhibitors.     

Evaluating the potency of HIV-1 protease drugs to combat resistance

HIV-1 protease has been an important drug target for the antiretroviral treatment of HIV infection. The efficacy of protease drugs is impaired by the rapid emergence of resistant virus strains. Understanding the molecular basis and evaluating the potency of an inhibitor to combat resistance are no doubt important in AIDS therapy. In this study, we first identified residues that have significant contributions to binding with six substrates using molecular dynamics simulations and Molecular Mechanics Generalized Born Surface Area calculations. Among the critical residues, Asp25, Gly27, Ala28, Asp29, and Gly49 are well conserved, with which the potent drugs should form strong interactions. We then calculated the contribution of each residue to binding with eight FDA approved drugs. We analyzed the conservation of each protease residue and also compared the interaction between the HIV protease and individual residues of the drugs and substrates. Our analyses showed that resistant mutations usually occur at less conserved residues forming more favorable interactions with drugs than with substrates. To quantitatively integrate the binding free energy and conservation information, we defined an empirical parameter called free energy/variability (FV) value, which is the product of the contribution of a single residue to the binding free energy and the sequence variability at that position. As a validation, the FV value was shown to identify single resistant mutations with an accuracy of 88%. Finally, we evaluated the potency of a newly approved drug, darunavir, to combat resistance and predicted that darunavir is more potent than amprenavir but may be susceptible to mutations on Val32 and Ile84.     

Determination of the protonation state of the Asp dyad: conventional molecular dynamics versus thermodynamic integration

The protonation state of the Asp dyad is important as it can reveal enzymatic mechanisms, and the information this provides can be used in the development of drugs for proteins such as memapsin 2 (BACE-1), HIV-1 protease, and rennin. Conventional molecular dynamics (MD) simulations have been successfully used to determine the preferred protonation state of the Asp dyad. In the present work, we demonstrate that the results obtained from conventional MD simulations can be greatly influenced by the particular force field applied or the values used for control parameters. In principle, free-energy changes between possible protonation states can be used to determine the protonation state. We show that protonation state prediction by the thermodynamic integration (TI) method is insensitive to force field version or to the cutoff for calculating nonbonded interactions (a control parameter). In the present study, the protonation state of the Asp dyad predicted by TI calculations was the same regardless of the force field and cutoff value applied. Contrary to the intuition that conventional MD is more efficient, our results clearly show that the TI method is actually more efficient and more reliable for determining the protonation state of the Asp dyad.     

Conformational changes in HIV-1 proteinase: effect of protonation of the active center on conformation of HIV-1 proteinase in water

At weak acidic pH, where HIV-1 proteinase is most stable and active, its catalytic Asp 25/25' dyad shares one proton. At a physiological pH the dyad is deprotonated, however, 2 ns molecular dynamics simulations of the HIV-1 protease with monoprotonated and deprotonated Asp25/25' dyad is performed, in order to investigate the influence of Asp25/25' protonation state on the proteinase dynamics. For net charge neutralization the 4 Cl- ions were included. In case of deprotonated active site the significant tertiary structure deviation of HIV-1 PR structure from crystal structure is observed, while in the monoprotonated one the tertiary structure fluctuates near starting structure. Possible mechanism of the influence of the Asp25/25' protonation state on proteinase dynamics is discussed.     

HIV-1跨膜蛋白gp41与N-取代吡咯衍生物的结合模式研究

采用分子对接,分子动力学(MD)模拟和分子力学/泊松-波尔兹曼溶剂可有面积方法与分子力学/广义伯恩溶剂可及面积方法(MM-PBSA/MM-GBSA),预测两种N-取代吡咯衍生物与HIV-1 跨膜蛋白gp41疏水口袋的结合模式与作用机理．分子对接采用多种受体构象,并从结果中选取几种可能的结合模式进行MD 模拟,然后通过MM-PBSA计算结合能的方法识别最优的结合模式. MM-PBSA计算结果表明,范德华相互作用是结合的主要驱动力,而极性相互作用决定了配体在结合过程中的取向．进一步的结合能分解显示,配体的羧基与gp41残基Arg579的静电相互作用对结合有重要贡献．上述工作为进一步优化N-取代吡咯衍生物类的HIV-1融合抑制剂建立了良好的理论基础．     

Probing the electronic structures and properties of neutral and anionic ScSi n (0,−1) (n = 1–6) clusters using ccCA-TM and G4 theory

We apply molecular docking, molecular dynamics (MD) simulation, and binding free energy calculation to investigate and reveal the binding mechanism between five xanthine inhibitors and DPP-4. The electrostatic and van der Waals interactions of the five inhibitors with DPP-4 are analyzed and discussed. The computed binding free energies using MM-PBSA method are in qualitatively agreement with experimental inhibitory potency of five inhibitors. The hydrogen bonds of inhibitors with Ser630 and Asp663 can stabilize the inhibitors in binding sites. The van der Waals interactions, especially the key contacts with His740, Asn710, Trp629, and Tyr666 have larger contributions to the binding free energy and play important roles in distinguishing the variant bioactivity of five inhibitors.     

Molecular dynamics simulations and MM–PBSA calculations of the lectin from snowdrop (Galanthus nivalis)

Galanthus nivalis agglutinin (GNA), a mannose-specific lectin from snowdrop bulbs, is a member of the monocot mannose-specific lectin family and exhibits antiviral activity toward HIV. In the present study, molecular dynamics (MD) simulations were performed to study the interaction between GNA and its carbohydrate ligand over a specific time span. By analysis of the secondary structures, it was observed that the GNA conformation maintains rather stable along the trajectories and the high fluctuations were only centered on the carbohydrate recognition domains. Our MD simulations also reproduced most of the hydrogen bonds observed in the x-ray crystal structure. Furthermore, the obtained MD trajectories were used to estimate the binding free energy of the complex using the molecular mechanics/Poisson Boltzmann surface area (MM-PBSA) method. It was revealed by the inspection of the binding free energy components that the major contributions to the complex stability arose from electrostatic interactions.     

Litchi chinensis-derived terpenoid as anti-HIV-1 protease agent: structural design from molecular dynamics simulations

The molecular structures of the binding between human immunodeficiency virus-1 protease (HIV-1PR) and various inhibitors including existing extensive natural products extracts have been investigated for anti-HIV drug development. In this study, the binding of HIV-1PR and a terpenoid from Litchi chinensis extracts (3-oxotrirucalla-7,24-dien-21-oic acid) was investigated in order to clarify the inhibition effectiveness of this compound. Molecular dynamics (MD) simulations of HIV-1PR complex with 3-oxotrirucalla-7,24-dien-21-oic acid were performed including water molecules. The MD simulation results indicated the formation of hydrogen bonds between the oxygen atoms of the inhibitor and the catalytic aspartates, which are commonly found in inhibitors–protease complexes. On the other hand, no hydrogen bonding of this particular inhibitor to the flap region was found. In addition, the radial distribution function of water oxygens around the catalytic carboxylate nitrogens of Asp29 and Asp30 suggests that at least one or two water molecules are in the active site region whereas direct interaction of the inhibitor was found for catalytic carboxylate oxygen of Asp25. The results of this simulation, in comparison with the structures of other HIV-PR inhibitor complexes, could lead to a better understanding of the activity of 3-oxotrirucalla-7,24-dien-21-oic acid.     

Investigation into mechanism of orotidine 5′-monophosphate decarboxylase enzyme by MM-PBSA/MM-GBSA and molecular docking

In this work, we studied the binding affinity of orotidine 5′-monophosphate (OMP) and 6-hydroxy-UMP (BMP) for Saccharomyces cerevisiae orotidine 5′-monophosphate decarboxylase (OMPDC) enzyme by using Molecular Mechanics-Poisson–Boltzmann Surface Area (MM-PBSA) and the Molecular Mechanics-Generalised Born Surface Area (MM-GBSA) calculations. In all simulations, Asp91, which is an important residue in the enzyme active site, was considered in both anionic (present in the native form of the enzyme) and neutral states. A series of 10-ns molecular dynamics simulations were performed for the four OMPDC–ligand complexes, two ligand-free enzymes and two free ligands, followed by MM-PBSA/MM-GBSA calculations on the collected snapshots, and molecular docking calculations using the free enzymes and ligands. The results of MM-PBSA/MM-GBSA calculations indicate that all of the OMPDC–ligand complexes form favourable systems in water, which is in agreement with corresponding experimental data. The results of the MM-PBSA and molecular docking methods also showed that OMPDC–BMP complexes, transition state analogue and inhibitor of the OMPDC enzyme have the highest binding affinities. The fact that in the native anionic state BMP shows a higher binding affinity compared with the substrate suggests the contribution of a transition state stabilisation mechanism in the debatable catalytic mechanism of the OMPDC enzyme.     

Application of free energy simulations to the binding of a transition-state-analogue inhibitor to HIV protease.

Free energy simulations (slow-change method) have been used to estimate quantitatively the ratio of the binding constants of (S) and (R) isomers of a novel HIV protease inhibitor, JG365. As a starting geometry, we used the X-ray crystallographic structure of a complex of HIV protease and JG365 provided by A. Wlodawer. According to our results the (S) configuration, i.e. the form previously identified experimentally, binds considerably more tightly to the protease (delta delta G degrees = 2.9 kcal/mol). When the (S) inhibitor is bound, there is a very strong preference for protonation of the Asp125 (rather than the Asp25) residue of the protease. This study is the first to apply a new method for quantitatively assessing the precision of free energies calculated by the slow-change method.     

Characterization of the binding of angiotensin II receptor blockers to human serum albumin using docking and molecular dynamics simulation

Human serum albumin (HSA), the most abundant protein found in blood plasma, transports many drugs and ligands in the circulatory system. The drug binding ability of HSA strongly influences free drug concentrations in plasma, and is directly related to the effectiveness of clinical therapy. In current work, binding of HSA to angiotensin II receptor blockers (ARBs) are investigated using docking and molecular dynamics (MD) simulations. Docking results demonstrate that the main HSA–ARB binding site is subdomain IIIA of HSA. Simulation results reveal clearly how HSA binds with valsartan and telmisartan. Interestingly, electrostatic interactions appear to be more important than hydrophobic interactions in stabilizing binding of valsartan to HSA, and vice versa for HSA–telmisartan. The molecular distance between HSA Trp214 (donor) and the drug (acceptor) can be measured by fluorescence resonance energy transfer (FRET) in experimental studies. The average distances between Trp-214 and ARBs are estimated here based on our MD simulations, which could be valuable to future FRET studies. This work will be useful in the design of new ARB drugs with desired HSA binding affinity.     

Moracin M from Morus alba L. is a natural phosphodiesterase-4 inhibitor

Phosphodiesterase-4 (PDE4) has been identified to be a promising target for treatment of asthma. Moracin M extracted from Chinese herbal drug 'Sang-Bai-Pi' (Morus alba L.) was studied for the inhibitory affinity towards PDE4. It inhibited PDE4D2, PDE4B2, PDE5A1, and PDE9A2 with the IC(50) values of 2.9, 4.5, >40, and >100 μM, respectively. Our molecular docking and 8ns molecular dynamics (MD) simulations demonstrated that moracin M forms three hydrogen bonds with Gln369, Asn321, and Asp318 in the active site and stacks against Phe372. In addition, comparative kinetics analysis of its analog moracin C was carried out to qualitatively validate their inhibitory potency as predicted by the binding free energy calculations after MD simulations.     

A molecular dynamics study of the structural stability of HIV-1 protease under physiological conditions: the role of Na+ ions in stabilizing the active site

HIV-1 protease is most active under weakly acidic conditions (pH 3.5-6.5), when the catalytic Asp25 and Asp25' residues share 1 proton. At neutral pH, this proton is lost and the stability of the structure is reduced. Here we present an investigation of the effect of pH on the dynamics of HIV-1 protease using MD simulation techniques. MD simulations of the solvated HIV-1 protease with the Asp25/25' residues monoprotonated and deprotonated have been performed. In addition we investigated the effect of the inclusion of Na(+) and Cl(-) ions to mimic physiological salt conditions. The simulations of the monoprotonated form and deprotonated form including Na(+) show very similar behavior. In both cases the protein remained stable in the compact, "self-blocked" conformation in which the active site is blocked by the tips of the flaps. In the deprotonated system a Na(+) ion binds tightly to the catalytic dyad shielding the repulsion between the COO(-) groups. Ab initio calculations also suggest the geometry of the active site with the Na(+) bound closely resembles that of the monoprotonated case. In the simulations of the deprotonated form (without Na(+) ions), a water molecule bound between the Asp25 Asp25' side-chains. This disrupted the dimerization interface and eventually led to a fully open conformation.     

Computational analysis of aspartic protease plasmepsin II complexed with EH58 inhibitor: a QM/MM MD study

Plasmepsin (PM) II is one of four enzymes in the food vacuole of Plasmodium falciparum. It has become an attractive target for combating malaria through research regarding its importance in the P. falciparum metabolism and life cycle, making it the target of choice for structure-based drug design. This paper reports the results of hybrid quantum mechanics / molecular mechanics (QM/MM) molecular dynamics (MD) simulations employed to study the details of the interactions established between PM II and N-(3-{(2-benzo[1, 3]dioxol-5-yl-ethyl)[3-(1-methyl-3-oxo-1,3-dihydro-isoindol-2-yl) propionyl]-amino}-1-benzyl-2-(hydroxyl-propyl)-4-benzyloxy-3,5dimethoxy-benzamide (EH58), a well-known potent inhibitor for this enzyme. Electrostatic binding free energy and energy terms decomposition have been computed for PM II complexed with the EH58 inhibitor. The results reveal that there is a strong interaction between Asp34, Val78, Ser79, Tyr192 and Asp214 residues and the EH58 inhibitor. In addition, we have computed the potential of the mean force (PMF) profile in order to assign the protonation state of the two catalytic aspartates in PM II-EH58 complex. The results indicate that the protonation of Asp214 favors a stable active site structure, which is consistent with our electrostatic binding free energy calculation and with previous published works.     

Insights into the interactions between HIV-1 integrase and human LEDGF/p75 by molecular dynamics simulation and free energy calculation

Human cellular protein LEDGF/p75 (lens epithelium-derived growth factor) is an important binding partner of human immunodeficiency virus type 1 (HIV-1) integrase (IN). Without LEDGF/p75, HIV-1 can not complete its life cycle. To study the detailed interactions between LEDGF/p75 and HIV-1 IN, and then obtain the hotspots at the binding interface, 13 ns molecular dynamics simulations were carried out here. One-hundred snapshots extracted from the last 4 ns trajectories were used for calculation of binding free energy and decomposition of the energy by residue. First, the structural changes and their dynamic interactions were investigated focused on the production stage. And then, the free energy was discussed. On the basis of the above results, it could be suggested that residues Gln168, Glu170, and Thr174 in chain A of IN, Thr125, and Trp131 in chain B of IN as well as Ile365, Asp366, Phe406, and Val408 in LEDGF/p75 were responsible for their binding. These results might be helpful for discovery and design of small molecules to interrupt the interaction between HIV-1 IN and LEDGF/p75.     

Molecular dynamics simulation of interleukin-2 and its complex and determination of the binding free energy

Interleukin-2 (IL-2) protein belongs to the signal modulator cytokine's family and therefore it is prevalent for immunological responses. It has been identified as a centrally important potential drug target for the inhibition of protein-protein interactions; so as to suppress the immunological responses associated with autoimmune, inflammatory and immunological diseases, and cancer. In the present work, we have performed two independent 100?ns of molecular dynamics (MD) simulations on the apo IL-2 protein and its ligand-bound complex (with a potent inhibitor FRG), to study the effect of inhibitor binding on the dynamics and stability of the protein. The calculation of binding free energy via post-processing end state method of Molecular Mechanics Poisson Boltzmann Surface Area (MM-PBSA) and Molecular Mechanics Generalised Born Surface Area (MM-GBSA) has inferred a good correlation in accordance with the already reported experimental data, demonstrating that the free energy of binding calculated by the two methods has no significant difference. The investigation of individual components of free energy revealed that the association of IL-2 protein with FRG ligand is primarily driven by the van der Waals energy contribution that represents the non-polar/hydrophobic energy contribution as dominant in this case of ligand binding.     

Calculation of binding energy using BLYP/MM for the HIV-1 integrase complexed with the S-1360 and two analogues

Integrase (IN) is one of the three human immunodeficiency virus type 1 (HIV-1) enzymes essential for effective viral replication. S-1360 is a potent and selective inhibitor of HIV-1 IN. In this work, we have carried out molecular dynamics (MD) simulations using a hybrid Quantum Mechanics/Molecular Mechanics (QM/MM) approach, to determine the protein-ligand interaction energy for S-1360 and two analogues. Analysis of the MD trajectories reveals that the strongest protein-inhibitor interactions, observed in the three studied complexes, are established with Lys-159 residue and Mg(2+) cation. Calculations of binding energy using BLYP/MM level of theory reveal that there is a direct relationship between this theoretical computed property and the experimental determined anti-HIV activity.     

Binding free energies for nicotine analogs inhibiting cytochrome P450 2A6 by a combined use of molecular dynamics simulations and QM/MM-PBSA calculations

Molecular dynamics (MD) simulations and hybrid quantum mechanical/molecular mechanical (QM/MM) calculations have been performed to explore the dynamic behaviors of cytochrome P450 2A6 (CYP2A6) binding with nicotine analogs (that are typical inhibitors) and to calculate their binding free energies in combination with Poisson–Boltzmann surface area (PBSA) calculations. The combined MD simulations and QM/MM-PBSA calculations reveal that the most important structural parameters affecting the CYP2A6-inhibitor binding affinity are two crucial internuclear distances, that is, the distance between the heme iron atom of CYP2A6 and the coordinating atom of the inhibitor, and the hydrogen-bonding distance between the N297 side chain of CYP2A6 and the pyridine nitrogen of the inhibitor. The combined MD simulations and QM/MM-PBSA calculations have led to dynamic CYP2A6-inhibitor binding structures that are consistent with the observed dynamic behaviors and structural features of CYP2A6-inhibitor binding, and led to the binding free energies that are in good agreement with the experimentally-derived binding free energies. The agreement between the calculated binding free energies and the experimentally-derived binding free energies suggests that the combined MD and QM/MM-PBSA approach may be used as a valuable tool to accurately predict the CYP2A6-inhibitor binding affinities in future computational design of new, potent and selective CYP2A6 inhibitors.     

Molecular dynamics (MD) simulations for the prediction of chiral discrimination of N-acetylphenylalanine enantiomers by cyclomaltoheptaose (beta-cyclodextrin, beta-CD) based on the MM-PBSA (molecular mechanics-Poisson-Boltzmann surface area) approach

Molecular dynamics (MD) simulations were performed for the prediction of chiral discrimination of N-acetylphenylalanine enantiomers by cyclomaltoheptaose (beta-cyclodextrin, beta-CD). Binding free energies and various conformational properties were obtained using by the MM-PBSA (molecular mechanics Poisson-Boltzmann/surface area) approach. The calculated relative difference (DeltaDeltabinding) of binding free energy was in fine agreement with the experimentally determined value. The difference of rotameric distributions of guest N-acetylphenylalanine enantiomers complexed with the host, beta-CD, was observed after the conformational analyses, suggesting that the conformational changes of guest captured within host cavity would be a decisive factor for enantiodifferentiation at a molecular level.