Some insights into mechanism for binding and drug resistance of wild type and I50V V82A and I84V mutations in HIV-1 protease with GRL-98065 inhibitor from molecular dynamic simulations

Decoding drug resistant mechanism of V32I, I50V and I84V mutations of HIV-1 protease on amprenavir binding by using molecular dynamics simulations and MM-GBSA calculations

Mutations V32I, I50V and I84V in the HIV-1 protease (PR) induce drug resistance towards drug amprenavir (APV). Multiple short molecular dynamics (MSMD) simulations and molecular mechanics generalized Born surface area (MM-GBSA) method were utilized to investigate drug-resistant mechanism of V32I, I50V and I84V towards APV. Dynamic information arising from MSMD simulations suggest that V32I, I50V and I84V highly affect structural flexibility, motion modes and conformational behaviours of two flaps in the PR. Binding free energies calculated by MM-GBSA method suggest that the decrease in binding enthalpy and the increase in binding entropy induced by mutations V32I, I50V and I84V are responsible for drug resistance of the mutated PRs on APV. The energetic contributions of separate residues on binding of APV to the PR show that V32I, I50V and I84V highly disturb the interactions of two flaps with APV and mostly drive the decrease in binding ability of APV to the PR. Thus, the conformational changes of two flaps in the PR caused by V32I, I50V and I84V play key roles in drug resistance of three mutated PR towards APV. This study can provide useful dynamics information for the design of potent inhibitors relieving drug resistance.

Decoding molecular mechanism underlying binding of drugs to HIV-1 protease with molecular dynamics simulations and MM-GBSA calculations

HIV-1 protease (PR) is thought to be efficient targets of anti-AIDS drug design. Molecular dynamics (MD) simulations and multiple post-processing analysis technologies were applied to decipher molecular mechanism underlying binding of three drugs Lopinavir (LPV), Nelfinavir (NFV) and Atazanavir (ATV) to the PR. Binding free energies calculated by molecular mechanics generalized Born surface area (MM-GBSA) suggest that compensation between binding enthalpy and entropy plays a vital role in binding of drugs to PR. Dynamics analyses show that binding of LPV, NFV and ATV highly affects structural flexibility, motion modes and dynamics behaviour of the PR, especially for two flaps. Computational alanine scanning and interaction network analysis verify that although three drugs have structural difference, they share similar binding modes to the PR and common interaction clusters with the PR. The current findings also confirm that residues located interaction clusters, such as Asp25/Asp25', Gly27/Gly27', Ala28/Ala28', Asp29, Ile47/Ile47', Gly49/Gly49', Ile50/Ile50', Val82/Val82' and Ile84/Ile84, can be used as efficient targets of clinically available inhibitors towards the PR.

Insights into effect of the Asp25/Asp25' protonation states on binding of inhibitors Amprenavir and MKP97 to HIV-1 protease using molecular dynamics simulations and MM-GBSA calculations

The protonation states of two aspartic acids in the catalytic strands of HIV-1 protease (PR) remarkably affect bindings of inhibitors to PR. It is requisite for the design of potent inhibitors towards PR to investigate the influences of Asp25/Asp25' protonated states on dynamics behaviour of PR and binding mechanism of inhibitors to PR. In this work, molecular dynamics (MD) simulations, MM-GBSA method and principal component (PC) analysis were coupled to explore the effect of Asp25/Asp25' protonation states on conformational changes of PR and bindings of Amprenavir and MKP97 to PR. The results show that the Asp25/Asp25' protonation states exert different impacts on structural fluctuations, flexibility and motion modes of PR. Dynamics analysis verifies that Asp25/Asp25' protonated states highly affect conformational dynamics of two flaps in PR. The binding free energy calculations results suggest that the Asp25/Asp25' protonated states obviously strengthen bindings of inhibitors to PR compared to the non-protonation state. Calculations of residue-based free energy decomposition indicate that the Asp25/Asp25' protonation not only disturbs the interaction network of inhibitors with PR but also stabilizes bindings of inhibitors to PR by cancelling the electrostatic repulsive interaction. Therefore, special attentions should be paid to the Asp25/Asp25' protonation in the design of potent inhibitors towards PR.

Non-active site mutants of HIV-1 protease influence resistance and sensitisation towards protease inhibitors

Background

HIV-1 can develop resistance to antiretroviral drugs, mainly through mutations within the target regions of the drugs. In HIV-1 protease, a majority of resistance-associated mutations that develop in response to therapy with protease inhibitors are found in the protease’s active site that serves also as a binding pocket for the protease inhibitors, thus directly impacting the protease-inhibitor interactions. Some resistance-associated mutations, however, are found in more distant regions, and the exact mechanisms how these mutations affect protease-inhibitor interactions are unclear. Furthermore, some of these mutations, e.g. N88S and L76V, do not only induce resistance to the currently administered drugs, but contrarily induce sensitivity towards other drugs. In this study, mutations N88S and L76V, along with three other resistance-associated mutations, M46I, I50L, and I84V, are analysed by means of molecular dynamics simulations to investigate their role in complexes of the protease with different inhibitors and in different background sequence contexts.

Results

Using these simulations for alchemical calculations to estimate the effects of mutations M46I, I50L, I84V, N88S, and L76V on binding free energies shows they are in general in line with the mutations’ effect on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$IC_{50}$$\end{document}IC50 values. For the primary mutation L76V, however, the presence of a background mutation M46I in our analysis influences whether the unfavourable effect of L76V on inhibitor binding is sufficient to outweigh the accompanying reduction in catalytic activity of the protease. Finally, we show that L76V and N88S changes the hydrogen bond stability of these residues with residues D30/K45 and D30/T31/T74, respectively.

Conclusions

We demonstrate that estimating the effect of both binding pocket and distant mutations on inhibitor binding free energy using alchemical calculations can reproduce their effect on the experimentally measured \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$IC_{50}$$\end{document}IC50 values. We show that distant site mutations L76V and N88S affect the hydrogen bond network in the protease’s active site, which offers an explanation for the indirect effect of these mutations on inhibitor binding. This work thus provides valuable insights on interplay between primary and background mutations and mechanisms how they affect inhibitor binding.

Exploration of the effects of sequence variations between HIV-1 and HIV-2 proteases on their three-dimensional structures

HIV protease inhibitors (PIs) approved by the FDA (US Food and Drug Administration) are a major class of antiretroviral. HIV-2 protease (PR2) is naturally resistant to most of them as PIs were designed for HIV-1 protease (PR1). In this study, we explored the impact of amino-acid substitutions between PR1 and PR2 on the structure of protease (PR) by comparing the structural variability of 13 regions using 24 PR1 and PR2 structures complexed with diverse ligands. Our analyses confirmed structural rigidity of the catalytic region and highlighted the important role of three regions in the conservation of the catalytic region conformation. Surprisingly, we showed that the flap region, corresponding to a flexible region, exhibits similar conformations in PR1 and PR2. Furthermore, we identified regions exhibiting different conformations in PR1 and PR2, which could be explained by the intrinsic flexibility of these regions, by crystal packing, or by PR1 and PR2 substitutions. Some substitutions induce structural changes in the R2 and R4 regions that could have an impact on the properties of PI-binding site and could thus modify PI binding mode. Substitutions involved in structural changes in the elbow region could alter the flexibility of the PR2 flap regions relative to PR1, and thus play a role in the transition from the semi-open form to the closed form, and have an impact on ligand binding. These results improve the understanding of the impact of sequence variations between PR1 and PR2 on the natural resistance of HIV-2 to commercially available PIs.Communicated by Ramaswamy H. Sarma.

Insight into binding mechanisms of inhibitors MKP56, MKP73, MKP86, and MKP97 to HIV-1 protease by using molecular dynamics simulation

HIV-1 protease (PR) has been a significant target for design of potent inhibitors curing acquired immunodeficiency syndrome. Molecular dynamics simulations coupled with molecular mechanics Poisson-Boltzmann surface area method were performed to study interaction modes of four inhibitors MKP56, MKP73, MKP86, and MKP97 with PR. The results suggest that the main force controlling interactions of inhibitors with PR should be contributed by van der Waals interactions between inhibitors and PR. The cross-correlation analyses based on MD trajectories show that inhibitor binding produces significant effect on the flap dynamics of PR. Hydrogen bond analyses indicate that inhibitors can form stable hydrogen bonding interactions with the residues from the catalytic strands of PR. The contributions of separate residues to inhibitor bindings are evaluated by using residue-based free energy decomposition method and the results demonstrate that the CH-π and CH-CH interactions between the hydrophobic groups of inhibitors with residues drive the associations of inhibitors with PR. We expect that this study can provide a significant theoretical aid for design of potent inhibitors targeting PR.

Effect of allosteric molecules on structure and drug affinity of HIV-1 protease by molecular dynamics simulations

Recent experiments show that small molecules can bind onto the allosteric sites of HIV-1 protease (PR), which provides a starting point for developing allosteric inhibitors. However, the knowledge of the effect of such binding on the structural dynamics and binding free energy of the active site inhibitor and PR is still lacking. Here, we report 200ns long molecular dynamics simulation results to gain insight into the influences of two allosteric molecules (1H-indole-6-carboxylic acid, 1F1 and 2-methylcyclohexano, 4D9). The simulations demonstrate that both allosteric molecules change the PR conformation and stabilize the structures of PR and the inhibitor; the residues of the flaps are sensitive to the allosteric molecules and the flexibility of the residues is pronouncedly suppressed; the additions of the small molecules to the allosteric sites strengthen the binding affinities of 3TL-PR by about 12-15kal/mol in the binding free energy, which mainly arises from electrostatic term. Interestingly, it is found that the action mechanisms of 1F1 and 4D9 are different, the former behaviors like a doorman that keeps the inhibitor from escape and makes the flaps (door) partially open; the latter is like a wedge that expands the allosteric space and meanwhile closes the flaps. Our data provide a theoretical support for designing the allosteric inhibitor.

Computational Studies of a Mechanism for Binding and Drug Resistance in the Wild Type and Four Mutations of HIV-1 Protease with a GRL-0519 Inhibitor

Drug resistance of mutations in HIV-1 protease (PR) is the most severe challenge to the long-term efficacy of HIV-1 PR inhibitor in highly active antiretroviral therapy. To elucidate the molecular mechanism of drug resistance associated with mutations (D30N, I50V, I54M, and V82A) and inhibitor (GRL-0519) complexes, we have performed five molecular dynamics (MD) simulations and calculated the binding free energies using the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method. The ranking of calculated binding free energies is in accordance with the experimental data. The free energy spectra of each residue and inhibitor interaction for all complexes show a similar binding model. Analysis based on the MD trajectories and contribution of each residues show that groups R2 and R3 mainly contribute van der Waals energies, while groups R1 and R4 contribute electrostatic interaction by hydrogen bonds. The drug resistance of D30N can be attributed to the decline in binding affinity of residues 28 and 29. The size of Val50 is smaller than Ile50 causes the residue to move, especially in chain A. The stable hydrophobic core, including the side chain of Ile54 in the wild type (WT) complex, became unstable in I54M because the side chain of Met54 is flexible with two alternative conformations. The binding affinity of Ala82 in V82A decreases relative to Val82 in WT. The present study could provide important guidance for the design of a potent new drug resisting the mutation inhibitors.

Drug resistance mechanisms of three mutations V32I, I47V and V82I in HIV-1 protease toward inhibitors probed by molecular dynamics simulations and binding free energy predictions

Molecular dynamics simulation and binding free energy calculations were used to probe drug resistance of HIV-1 protease mutations toward inhibitors.

Probing Origin of Binding Difference of inhibitors to MDM2 and MDMX by Polarizable Molecular Dynamics Simulation and QM/MM-GBSA Calculation

Binding abilities of current inhibitors to MDMX are weaker than to MDM2. Polarizable molecular dynamics simulations (MD) followed by Quantum mechanics/molecular mechanics generalized Born surface area (QM//MM-GBSA) calculations were performed to investigate the binding difference of inhibitors to MDM2 and MDMX. The predicted binding free energies not only agree well with the experimental results, but also show that the decrease in van der Walls interactions of inhibitors with MDMX relative to MDM2 is a main factor of weaker bindings of inhibitors to MDMX. The analyses of dihedral angles based on MD trajectories suggest that the closed conformation formed by the residues M53 and Y99 in MDMX leads to a potential steric clash with inhibitors and prevents inhibitors from arriving in the deep of MDMX binding cleft, which reduces the van der Waals contacts of inhibitors with M53, V92, P95 and L98. The calculated results using the residue-based free energy decomposition method further prove that the interaction strength of inhibitors with M53, V92, P95 and L98 from MDMX are obviously reduced compared to MDM2. We expect that this study can provide significant theoretical guidance for designs of potent dual inhibitors to block the p53-MDM2/MDMX interactions.

A comparative study of trypsin specificity based on QM/MM molecular dynamics simulation and QM/MM GBSA calculation

Hydrogen bonding and polar interactions play a key role in identification of protein-inhibitor binding specificity. Quantum mechanics/molecular mechanics molecular dynamics (QM/MM MD) simulations combined with DFT and semi-empirical Hamiltonian (AM1d, RM1, PM3, and PM6) methods were performed to study the hydrogen bonding and polar interactions of two inhibitors BEN and BEN1 with trypsin. The results show that the accuracy of treating the hydrogen bonding and polar interactions using QM/MM MD simulation of PM6 can reach the one obtained by the DFT QM/MM MD simulation. Quantum mechanics/molecular mechanics generalized Born surface area (QM/MM-GBSA) method was applied to calculate binding affinities of inhibitors to trypsin and the results suggest that the accuracy of binding affinity prediction can be significantly affected by the accurate treatment of the hydrogen bonding and polar interactions. In addition, the calculated results also reveal the binding specificity of trypsin: (1) the amidinium groups of two inhibitors generate favorable salt bridge interaction with Asp189 and form hydrogen bonding interactions with Ser190 and Gly214, (2) the phenyl of inhibitors can produce favorable van der Waals interactions with the residues His58, Cys191, Gln192, Trp211, Gly212, and Cys215. This systematic and comparative study can provide guidance for the choice of QM/MM MD methods and the designs of new potent inhibitors targeting trypsin.

Exploring the drug resistance of V32I and M46L mutant HIV-1 protease to inhibitor TMC114: flap dynamics and binding mechanism

Inhibitors of HIV-1 protease (HIV-1-pr) generally only bind to the active site of the protease. However, for some mutants such as V32I and M46L the TMC114 can bind not only to the active cavity but also to the groove of the flexible flaps. Although the second binding site suggests the higher efficiency of the drug against HIV-1-pr, the drug resistance in HIV-1-pr due to mutations cannot be ignored, which prompts us to investigate the molecular mechanisms of drug resistance and behavior of double bound TMC114 (2T) to HIV-1-pr. The conformational dynamics of HIV-1-pr and the binding of TMC114 to the WT, V32I and M46L mutants were investigated with all-atom molecular dynamic (MD) simulation. The 20 ns MD simulation shows many fascinating effects of the inhibitor binding to the WT and mutant proteases. MM-PBSA calculations explain the binding free energies unfavorable for the M46L and V32I mutants as compared to the WT. For the single binding (1T) the less binding affinity can be attributed to the entropic loss for both V32I-1T and M46L-1T. Although the second binding of TMC114 with flap does increase binding energy for the mutants (V32I-2T and M46L-2T), the considerable entropy loss results in the lower binding Gibbs free energies. Thus, binding of TMC114 in the flap region does not help much in the total gain in binding affinity of the system, which was verified from this study and thereby validating experiments.

Revealing origin of decrease in potency of darunavir and amprenavir against HIV-2 relative to HIV-1 protease by molecular dynamics simulations

Clinical inhibitors Darunavir (DRV) and Amprenavir (APV) are less effective on HIV-2 protease (PR2) than on HIV-1 protease (PR1). To identify molecular basis associated with the lower inhibition, molecular dynamics (MD) simulations and molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) calculations were performed to investigate the effectiveness of the PR1 inhibitors DRV and APV against PR1/PR2. The rank of predicted binding free energies agrees with the experimental determined one. Moreover, our results show that two inhibitors bind less strongly to PR2 than to PR1, again in agreement with the experimental findings. The decrease in binding free energies for PR2 relative to PR1 is found to arise from the reduction of the van der Waals interactions induced by the structural adjustment of the triple mutant V32I, I47V and V82I. This result is further supported by the difference between the van der Waals interactions of inhibitors with each residue in PR2 and in PR1. The results from the principle component analysis suggest that inhibitor binding tends to make the flaps of PR2 close and the one of PR1 open. We expect that this study can theoretically provide significant guidance and dynamics information for the design of potent dual inhibitors targeting PR1/PR2.

Molecular dynamic simulations reveal the mechanism of binding between xanthine inhibitors and DPP-4

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.

Interchain hydrophobic clustering promotes rigidity in HIV-1 protease flap dynamics: new insights from molecular dynamics

The dynamics of HIV-1 protease (HIV-pr), a drug target for HIV infection, has been studied extensively by both computational and experimental methods. The flap dynamics of HIV-pr is considered to be more important for better ligand binding and enzymatic actions. Moreover, it has been demonstrated that the drug-induced mutations can change the flap dynamics of HIV-pr affecting the binding affinity of the ligands. Therefore, detailed understanding of flap dynamics is essential for designing better inhibitors. Previous computational investigations observed significant variation in the flap opening in nanosecond time scale indicating that the dynamics is highly sensitive to the simulation protocols. To understand the sensitivity of the flap dynamics on the force field and simulation protocol, molecular dynamics simulations of HIV-pr have been performed with two different AMBER force fields, ff99 and ff02. Two different trajectories (20 ns each) were obtained using the ff99 and ff02 force field. The results showed polarizable force field (ff02) make the flap tighter than the nonpolarizable force field (ff99). Some polar interactions and hydrogen bonds involving flap residues were found to be stronger with ff02 force field. The formation of interchain hydrophobic cluster (between flap tip of one chain and active site wall of another chain) was found to be dominant in the semi-open structures obtained from the simulations irrespective of the force field. It is proposed that an inhibitor, which will promote this interchain hydrophobic clustering, may make the flaps more rigid, and presumably the effect of mutation would be small on ligand binding.

Ligand selectivity of estrogen receptors by a molecular dynamics study

Estrogen receptors α (ERα) and β (ERβ) have different physiological functions and expression levels in different tissues. ERα and ERβ are highly homologous and have only two residue substitutions in the binding pocket. This high similarity at the active site stimulates the interests for discovering subtype selective ligands. In this study, molecular dynamics (MD) simulations combined with molecular mechanics generalized Born surface area (MM-GBSA) method have been carried out to analyze the basis of selectivity of three ligands (659, 818 and 041). The calculated binding free energies show that all the ligands bind more tightly to ERβ than to ERα. The dominant free energy components of selectivity for 659 are similar to that for 041, but different from that for 818. The decompositions of free energy contributions and structural analysis imply that there are eight residues primarily contributing to the selectivity for 659, five residues for 041, as well as two residues for 818. The structural analysis implies that two residue substitutions in binding packet cause the position of 659 in ERβ-659 complex to shift relative to that in ERα-659 complex and also cause the conformational changes of other residues in the binding pocket. The higher selectivity for 041 is mainly caused by three residues, Ile373 (Met421), His475 (His524) and Leu476 (Leu525).

Molecular dynamic simulations give insight into the mechanism of binding between 2-aminothiazole inhibitors and CDK5

Molecular docking, molecular dynamics (MD) simulations, and binding free energy analysis were performed to reveal differences in the binding affinities between five 2-aminothiazole inhibitors and CDK5. The hydrogen bonding and hydrophobic interactions between inhibitors and adjacent residues are analyzed and discussed. The rank of calculated binding free energies using the MM-PBSA method is consistent with experimental result. The results illustrate that hydrogen bonds with Cys83 favor inhibitor binding. The van der Waals interactions, especially the important contact with Ile10, dominate in the binding free energy and play a crucial role in distinguishing the different bioactivity of the five inhibitors.

Molecular dynamics and free energy studies on the carboxypeptidases complexed with peptide/small molecular inhibitor: mechanism for drug resistance

As one potent plant protease inhibitor, potato carboxypeptidase inhibitor (PCI) can competitively inhibit insect digestive metallocarboxypeptidases (MCPs) through interfering with its digestive system that causes amino acid deficiencies and leading to serious developmental delay and mortality. However, this effective biological pest control is significantly impaired by the PCI-resistant insect MCPs. Therefore, deep understanding of the resistant mechanism of insect MCPs is particularly necessary for designing new durable pest control regimen and developing effective pesticides. In this study, the binding of PCI and small molecular inhibitor THI to insect PCI-sensitive/-resistant MCPs and human MCP was investigated by docking, molecular dynamics (MD) simulations and thermodynamic analysis. The structural analysis from MD simulations indicates that the PCI-resistant mechanism of CPBHz is mainly dominated by the Trp277A, which changes the conformation of β8-α9 loop and therefore narrow the access to the active site of CPBHz, prohibiting the entrance of the C termini tail of PCI. Additionally, the insertion of Gly247A weakens the stabilization of CPBHz and PCI through disrupting the hydrogen bond formation with its surrounding residues. Furthermore, the predicted binding free energies gives explanation of structure affinity relationship of PCI and THI with MCPs and suggest that the electrostatic energy is the main contribution term affecting the difference in binding affinities. Finally, the decomposition analysis of binding free energies infers that the key residues Glu72, Arg127, Ile247/Leu247 and Glu270 are critical for the binding of PCI/THI to MCPs.

Interaction of I50V mutant and I50L/A71V double mutant HIV-protease with inhibitor TMC114 (darunavir): molecular dynamics simulation and binding free energy studies

In the present work, the binding of inhibitor TMC114 (darunavir) to wild-type (WT), single (I50V) as well as double (I50L/A71V) mutant HIV-proteases (HIV-pr) was investigated with all-atom molecular dynamics (MD) simulations as well as molecular mechanic-Poisson-Boltzmann surface area (MM-PBSA) calculation. For both the apo and complexed HIV-pr, many intriguing effects due to double mutant, I50L/A71V, are observed. For example, the flap-flap distance and the distance from the active site to the flap residues in the apo I50L/A71V-HIV-pr are smaller than those of WT- and I50V-HIV-pr, probably making the active site smaller in volume and closer movement of flaps. For the complexed HIV-pr with TMC114, the double mutant I50L/A71V shows a less curling of the flap tips and less flexibility than WT and the single mutant I50V. As for the other previous studies, the present results also show that the single mutant I50V decreases the binding affinity of I50V-HIV-pr to TMC, resulting in a drug resistance; whereas the double mutant I50L/A71V increases the binding affinity, and as a result of the stronger binding, the I50L/A71V may be well adapted by the TMC114. The energy decomposition analysis suggests that the increase of the binding for the double mutant I50L/A71V-HIV-pr can be mainly attributed to the increase in electrostatic energy by -5.52 kacl/mol and van der Waals by -0.42 kcal/mol, which are canceled out in part by the increase of polar solvation energy of 1.99 kcal/mol. The I50L/A71V mutant directly increases the binding affinity by approximately -0.88 (Ile50 to Leu50) and -0.90 (Ile50' to Leu50') kcal/mol, accounting 45% for the total gain of the binding affinity. Besides the direct effects from the residues Leu50 and Leu50', the residue Gly49' increases the binding affinity of I50L/A71V-HIV-pr to the inhibitor by -0.74 kcal/mol, to which the electrostatic interaction of Leu50's backbone contributes by -1.23 kcal/mol. Another two residues Ile84 and Ile47' also increase the binding affinity by -0.22 and -0.29 kcal/mol, respectively, which can be mainly attributed to van der Waals terms (ΔT(vdw) = -0.21 and -0.39 kcal/mol).

Insights into scFv:drug binding using the molecular dynamics simulation and free energy calculation

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

A computational analysis of the binding model of MDM2 with inhibitors

It is a new and promising strategy for anticancer drug design to block the MDM2-p53 interaction using a non-peptide small-molecule inhibitor. We carry out molecular dynamics simulations to study the binding of a set of six non-peptide small-molecule inhibitors with the MDM2. The relative binding free energies calculated using molecular mechanics Poisson-Boltzmann surface area method produce a good correlation with experimentally determined results. The study shows that the van der Waals energies are the largest component of the binding free energy for each complex, which indicates that the affinities of these inhibitors for MDM2 are dominated by shape complementarity. The A-ligands and the B-ligands are the same except for the conformation of 2,2-dimethylbutane group. The quantum mechanics and the binding free energies calculation also show the B-ligands are the more possible conformation of ligands. Detailed binding free energies between inhibitors and individual protein residues are calculated to provide insights into the inhibitor-protein binding model through interpretation of the structural and energetic results from the simulations. The study shows that G1, G2 and G3 group mimic the Phe19, Trp23 and Leu26 residues in p53 and their interactions with MDM2, but the binding model of G4 group differs from the original design strategy to mimic Leu22 residue in p53.