Binding of novel fullerene inhibitors to HIV-1 protease: insight through molecular dynamics and molecular mechanics Poisson–Boltzmann surface area calculations

The objectives of this study include the design of a series of novel fullerene-based inhibitors for HIV-1 protease (HIV-1 PR), by employing two strategies that can also be applied to the design of inhibitors for any other target. Additionally, the interactions which contribute to the observed exceptionally high binding free energies were analyzed. In particular, we investigated: (1) hydrogen bonding (H-bond) interactions between specific fullerene derivatives and the protease, (2) the regions of HIV-1 PR that play a significant role in binding, (3) protease changes upon binding and (4) various contributions to the binding free energy, in order to identify the most significant of them. This study has been performed by employing a docking technique, two 3D-QSAR models, molecular dynamics (MD) simulations and the molecular mechanics Poisson–Boltzmann surface area (MM–PBSA) method. Our computed binding free energies are in satisfactory agreement with the experimental results. The suitability of specific fullerene derivatives as drug candidates was further enhanced, after ADMET (absorption, distribution, metabolism, excretion and toxicity) properties have been estimated to be promising. The outcomes of this study revealed important protein–ligand interaction patterns that may lead towards the development of novel, potent HIV-1 PR inhibitors.     

Solution structure of HIV-1 protease flaps probed by comparison of molecular dynamics simulation ensembles and EPR experiments

The introduction of multidrug treatment regimens has dramatically prolonged the progression and survival of AIDS patients. However, the success of the long-term treatment has been hindered by strains of HIV that are increasingly resistant to inhibitors of targets such as HIV protease (HIV PR). Therefore, the need for a thorough understanding of the structure and dynamics of HIV PR and how these are altered in resistant mutants is crucial for the design of more effective treatments. Crystal structures of unbound HIV PR show significant heterogeneity and often have extensive crystal packing interactions. Recent site-directed spin labeling (SDSL) and double electron-electron resonance (DEER) spectroscopy studies characterized flap conformations in HIV-1 protease in an inhibited and uninhibited form and distinguished the extent of flap opening in an unbound form. However, the correlation between EPR-measured interspin distances and structural/dynamic features of the flaps has not been established. In this report, we link EPR-based data and 900 ns of MD simulation in explicit water to gain insight into the ensemble of conformations sampled by HIV PR flaps in solution, both in the presence and in the absence of an FDA-approved HIV PR inhibitor.     

Insights into the dynamics of HIV-1 protease: a kinetic network model constructed from atomistic simulations

The conformational dynamics in the flaps of HIV-1 protease plays a crucial role in the mechanism of substrate binding. We develop a kinetic network model, constructed from detailed atomistic simulations, to determine the kinetic mechanisms of the conformational transitions in HIV-1 PR. To overcome the time scale limitation of conventional molecular dynamics (MD) simulations, our method combines replica exchange MD with transition path theory (TPT) to study the diversity and temperature dependence of the pathways connecting functionally important states of the protease. At low temperatures the large-scale flap opening is dominated by a small number of paths; at elevated temperatures the transition occurs through many structurally heterogeneous routes. The expanded conformation in the crystal structure 1TW7 is found to closely mimic a key intermediate in the flap-opening pathways at low temperature. We investigated the different transition mechanisms between the semi-open and closed forms. The calculated relaxation times reveal fast semi-open ? closed transitions, and infrequently the flaps fully open. The ligand binding rate predicted from this kinetic model increases by 38-fold from 285 to 309 K, which is in general agreement with experiments. To our knowledge, this is the first application of a network model constructed from atomistic simulations together with TPT to analyze conformational changes between different functional states of a natively folded protein.     

HIV-1 protease flaps spontaneously close to the correct structure in simulations following manual placement of an inhibitor into the open state

We report unrestrained, all-atom molecular dynamics simulations of HIV-1 protease (HIV-PR) with a continuum solvent model that reproducibly sample closing of the active site flaps following manual placement of a cyclic urea inhibitor into the substrate binding site of the open protease. The open form was obtained from the unbound, semi-open HIV-PR crystal structure, which we recently reported (Hornak, V.; et al. Proc. Natl. Acad. Sci. U.S.A. 2006, 103, 915-920.) to have spontaneously opened during unrestrained dynamics. In those simulations, the transiently open flaps always returned to the semi-open form that is observed in all crystal structures of the free protease. Here, we show that manual docking of the inhibitor reproducibly induces spontaneous conversion to the closed form as seen in all inhibitor-bound HIV-PR crystal structures. These simulations reproduced not only the greater degree of flap closure, but also the striking difference in flap "handedness" between bound and free enzyme. In most of the simulations, the final structures were highly accurate. Root-mean-square deviations (RMSD) from the crystal structure of the complex were approximately 1.5 A (averaged over the last 100 ps) for the inhibitor and each flap despite initial RMSD of 2-5 A for the inhibitors and 6-11 A for the flaps. Key hydrogen bonds were formed between the flap tips and between flaps and inhibitor that match those seen in the crystal structure. The results demonstrate that all-atom simulations have the ability to significantly improve poorly docked ligand conformations and reproduce large-scale receptor conformational changes that occur upon binding.     

Combining Molecular Dynamic Information and an Aspherical-Atom Data Bank in the Evaluation of the Electrostatic Interaction Energy in Multimeric Protein-Ligand Complex: A Case Study for HIV-1 Protease

Computational analysis of protein–ligand interactions is of crucial importance for drug discovery. Assessment of ligand binding energy allows us to have a glimpse of the potential of a small organic molecule to be a ligand to the binding site of a protein target. Available scoring functions, such as in docking programs, all rely on equations that sum each type of protein–ligand interactions in order to predict the binding affinity. Most of the scoring functions consider electrostatic interactions involving the protein and the ligand. Electrostatic interactions constitute one of the most important part of total interactions between macromolecules. Unlike dispersion forces, they are highly directional and therefore dominate the nature of molecular packing in crystals and in biological complexes and contribute significantly to differences in inhibition strength among related enzyme inhibitors. In this study, complexes of HIV-1 protease with inhibitor molecules (JE-2147 and darunavir) were analyzed by using charge densities from the transferable aspherical-atom University at Buffalo Databank (UBDB). Moreover, we analyzed the electrostatic interaction energy for an ensemble of structures, using molecular dynamic simulations to highlight the main features of electrostatic interactions important for binding affinity.     

Molecular dynamic and free energy studies of primary resistance mutations in HIV-1 protease-ritonavir complexes

To understand the basis of drug resistance of the HIV-1 protease, molecular dynamic (MD) and free energy calculations of the wild-type and three primary resistance mutants, V82F, I84V, and V82F/I84V, of HIV-1 protease complexed with ritonavir were carried out. Analysis of the MD trajectories revealed overall structures of the protein and the hydrogen bonding of the catalytic residues to ritonavir were similar in all four complexes. Substantial differences were also found near the catalytic binding domain, of which the double mutant complex has the greatest impact on conformational changes of the protein and the inhibitor. The tip of the HIV-1 protease flap of the double mutant has the greater degree of opening with respect to that of the others. Additionally, the phenyl ring of Phe82 moves away from the binding pocket S1', and the conformational change of ritonavir subsite P1' consequently affects the cavity size of the protein and the conformational energy of the inhibitor. Calculations of binding free energy using the solvent continuum model were able to reproduce the same trend of the experimental inhibition constant. The results show that the resistance mutants require hydrophobic residues to maintain the interactions in the binding pocket. Changes of the cavity volume correlate well with free energy penalties due to the mutation and are responsible for the loss of drug susceptibility.     

Quantum study of HIV-1 protease-bridge water interaction

We present a fully quantum mechanical calculation for binding interaction between HIV-1 protease (PR) and the water molecule W301 which bridges the flaps of the protease with the inhibitors of PR. The quantum calculation is made possible by applying a recently developed molecular fractionation with conjugate caps (MFCC) method which divides a protein molecule into capped amino acid-based fragments and their conjugate caps. These individual fragments are properly treated to preserve the chemical property of bonds that are cut. Ab initio methods at HF, B3LYP, and MP2 levels with a fixed basis set 6-31+G* have been employed in the present calculation. The MFCC calculation produces a quantum mechanical interaction "map" representing interactions between individual residues of PR and W301. This enables a detailed quantitative analysis on binding of W301 to specific residues of PR at quantum mechanical level.     

Prediction and analysis of binding affinities for chemically diverse HIV-1 PR inhibitors by the modified SAFE_p approach

One of the biggest challenges in the "in silico" screening of enzyme ligands is to have a protocol that could predict the ligand binding free energies. In our group we have developed a very simple screening function (referred to as solvent accessibility free energy of binding predictor, SAFE_p) which we have applied previously to the study of peptidic HIV-1 protease (HIV-1 PR) inhibitors and later to cyclic urea type HIV-1 PR inhibitors. In this work, we have extended the SAFE_p protocol to a chemically diverse set of HIV-1 PR inhibitors with binding constants that differ by several orders of magnitude. The resulting function is able to reproduce the ranking and in many cases the value of the inhibitor binding affinities for the HIV-1 PR, with accuracy comparable with that of costlier protocols. We also demonstrate that the binding pocket SAFE_p analysis can contribute to the understanding of the physical forces that participate in ligand binding. The analysis tools afforded by our protocol have allowed us to identify an induced fit phenomena mediated by the inhibitor and have demonstrated that larger fragments do not necessarily contribute the most to the binding free energy, an outcome partially brought about by the substantial role the desolvation penalty plays in the energetics of binding. Finally, we have revisited the effect of the Asp dyad protonation state on the predicted binding affinities.     

Drug-resistant molecular mechanism of CRF01_AE HIV-1 protease due to V82F mutation

Human immunodeficiency virus type 1 protease (HIV-1 PR) is one of the major targets of anti-AIDS drug discovery. The circulating recombinant form 01 A/E (CRF01_AE, abbreviated AE) subtype is one of the most common HIV-1 subtypes, which is infecting more humans and is expanding rapidly throughout the world. It is, therefore, necessary to develop inhibitors against subtype AE HIV-1 PR. In this work, we have performed computer simulation of subtype AE HIV-1 PR with the drugs lopinavir (LPV) and nelfinavir (NFV), and examined the mechanism of resistance of the V82F mutation of this protease against LPV both structurally and energetically. The V82F mutation at the active site results in a conformational change of 79′s loop region and displacement of LPV from its proper binding site, and these changes lead to rotation of the side-chains of residues D25 and I50′. Consequently, the conformation of the binding cavity is deformed asymmetrically and some interactions between PR and LPV are destroyed. Additionally, by comparing the interactive mechanisms of LPV and NFV with HIV-1 PR we discovered that the presence of a dodecahydroisoquinoline ring at the P1′ subsite, a [2-(2,6-dimethylphenoxy)acetyl]amino group at the P2′ subsite, and an N2 atom at the P2 subsite could improve the binding affinity of the drug with AE HIV-1 PR. These findings are helpful for promising drug design. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Conformational analysis of TMC114, a novel HIV-1 protease inhibitor

TMC114, a potent novel HIV-1 protease inhibitor, remains active against a broad spectrum of mutant viruses. In order to bind to a variety of mutants, the compound needs to make strong, preferably backbone, interactions and have enough conformational flexibility to adapt to the changing geometry of the active site. The conformational analysis of TMC114 in the gas phase yielded 43 conformers in which five types of intramolecular H-bond interactions could be observed. All 43 conformers were subject to both rigid and flexible ligand docking in the wild-type and a triple mutant (L63P/V82T/I84V) of HIV-1 protease. The largest binding energy was calculated for the conformations that are close to the conformation observed in the X-ray complexes of TMC114 and HIV-1 protease.     

An insight to the molecular interactions of the FDA approved HIV PR drugs against L38L↑N↑L PR mutant

The aspartate protease of the human immune deficiency type-1 virus (HIV-1) has become a crucial antiviral target in which many useful antiretroviral inhibitors have been developed. However, it seems the emergence of new HIV-1 PR mutations enhances drug resistance, hence, the available FDA approved drugs show less activity towards the protease. A mutation and insertion designated L38L↑N↑L PR was recently reported from subtype of C-SA HIV-1. An integrated two-layered ONIOM (QM:MM) method was employed in this study to examine the binding affinities of the nine HIV PR inhibitors against this mutant. The computed binding free energies as well as experimental data revealed a reduced inhibitory activity towards the L38L↑N↑L PR in comparison with subtype C-SA HIV-1 PR. This observation suggests that the insertion and mutations significantly affect the binding affinities or characteristics of the HIV PIs and/or parent PR. The same trend for the computational binding free energies was observed for eight of the nine inhibitors with respect to the experimental binding free energies. The outcome of this study shows that ONIOM method can be used as a reliable computational approach to rationalize lead compounds against specific targets. The nature of the intermolecular interactions in terms of the host–guest hydrogen bond interactions is discussed using the atoms in molecules (AIM) analysis. Natural bond orbital analysis was also used to determine the extent of charge transfer between the QM region of the L38L↑N↑L PR enzyme and FDA approved drugs. AIM analysis showed that the interaction between the QM region of the L38L↑N↑L PR and FDA approved drugs are electrostatic dominant, the bond stability computed from the NBO analysis supports the results from the AIM application. Future studies will focus on the improvement of the computational model by considering explicit water molecules in the active pocket. We believe that this approach has the potential to provide information that will aid in the design of much improved HIV-1 PR antiviral drugs.     

Fast peptidyl cis-trans isomerization within the flexible Gly-rich flaps of HIV-1 protease

The catalytic aspartyl protease of the HIV-1 virus is a homodimer with two flaps that control access to the active site and are known to be flexible. However, knowledge of the atomistic mechanism of the flexibility is lacking. We show that the Gly-Gly omega-bond in the glycine-rich flap tips undergoes fast cis-trans isomerization on the microsecond to millisecond time scale rather than in the usual seconds. Further study reveals that the unexpectedly fast isomerization is a direct consequence of the beta-hairpin loop structure of the flap tips, which appears to be counterintuitive. After loop formation of a linear peptide containing the Gly-Gly motif, the rate of isomerization is shown to increase by many orders of magnitude.     

Interpretation of protein/ligand crystal structure using QM/MM calculations: case of HIV-1 protease/metallacarborane complex

Deltahedral metallacarborane compounds have recently been discovered as potent, specific, stable, and nontoxic inhibitors of HIV-1 protease (PR), the major target for AIDS therapy. The 2.15 A-resolution X-ray structure has exhibited a nonsymmetrical binding of the parental compound [Co(3+)-(C2B9H11)2](-) (GB-18) into PR dimer and a symmetrical arrangement in the crystal of two PR dimer complexes into a tetramer. In order to explore structural and energetic details of the inhibitor binding, quantum mechanics coupled with molecular mechanics approach was utilized. Realizing the close positioning of anionic inhibitors in the active site cavity, the possibility of an exchange of structural water molecules Wat50 and Wat128 by Na+ counterions was studied. The energy profiles for the rotation of the GB-18 molecules along their longitudinal axes in complex with PR were calculated. The results show that two Na+ counterions are present in the active site cavity and provide energetically favorable and unfavorable positions for carbon atoms within the carborane cages. Eighty-one rotamer combinations of four molecules of GB-18 bound to PR out of 4 x 10(5) are predicted to be highly populated. These results lay ground for further calculations of interaction energies between GB-18 and amino acids of PR active site and will make it possible to interpret computationally the binding of similar metallacarborane molecules to PR as well as to resistant PR variants. Moreover, this computational tool will allow the design of new, more potent metallacarborane-based HIV-1 protease inhibitors.     

Structure-based QSAR analysis of a set of 4-hydroxy-5,6-dihydropyrones as inhibitors of HIV-1 protease: an application of the receptor-dependent (RD) 4D-QSAR formalism

Receptor-dependent (RD) 4D-QSAR models were constructed for a set of 39 4-hydroxy-5,6-dihydropyrone analogue HIV-1 protease inhibitors. The receptor model used in this QSAR analysis was derived from the HIV-1 protease (PDB ID ) crystal structure. The bound ligand in the active site of the enzyme, also a 4-hydroxy-5,6-dihydropyrone analogue, was used as the reference ligand for docking the data set compounds. The optimized RD 4D-QSAR models are not only statistically significant (r(2) = 0.86, q(2) = 0.80 for four- and greater-term models) but also possess reasonable predictivity based on test set predictions. The proposed "active" conformations of the docked analogues in the active site of the enzyme are consistent in overall molecular shape with those suggested from crystallographic studies. Moreover, the RD 4D-QSAR models also "capture" the existence of specific induced-fit interactions between the enzyme active site and each specific inhibitor. Hydrophobic interactions, steric shape requirements, and hydrogen bonding of the 4-hydroxy-5,6-dihydropyrone analogues with the HIV-1 protease binding site model dominate the RD 4D-QSAR models in a manner again consistent with experimental conclusions. Some possible hypotheses for the development of new lead HIV-1 protease inhibitors can be inferred from the RD 4D-QSAR models.     

基于“底物包膜”假说筛选新型HIV-1蛋白酶抑制剂

基于“底物包膜”假说, 以现有HIV-1蛋白酶抑制剂Darunavir为模板构建药效团模型并对中药化学数据库进行搜索; 采用分子对接方法进一步考察化合物与HIV-1蛋白酶结合情况及其与“底物包膜”符合程度, 优先选出两个化合物Annomonicin和去乙酰蟾蜍它灵; 应用分子动力学方法对这两个化合物进行动力学模拟, 观察它们与蛋白酶结合的复合物在动力学过程中的稳定性并计算其结合自由能, 综合评价筛选结果, 最终确定化合物Annomonicin具有更潜在的深入研究价值.     

The open structure of a multi-drug-resistant HIV-1 protease is stabilized by crystal packing contacts

The introduction of HIV-1 protease (HIV-PR) inhibitors has led to a dramatic increase in patient survival; however, these gains are threatened by the emergence of multi-drug-resistant strains. Design of inhibitors that overcome resistance would be greatly facilitated by deeper insight into the mechanistic events associated with binding of substrates and inhibitors, as well as an understanding of the effects of resistance mutations on the structure and dynamic behavior of HIV-PR. We previously reported a series of simulations that provide a model for HIV-PR dynamics, with spontaneous conversions between the bound and unbound crystal forms upon addition or removal of an inhibitor. Importantly, the unbound protease transiently sampled a third fully open state that permits entry to the active site, unlike both crystallographic forms. Recently, a crystal structure of unbound HIV-PR was reported for the MDR 769 isolate (PDB: 1TW7); unlike all previous experimental structures, the binding pocket is open. It is suggested that drug resistance in this strain arises at least in part from the inability of inhibitors to induce closing. We carried out simulations of the MDR 769 HIV-PR mutant and observed that the reported structure is unstable in solution and rapidly adopts the semi-open conformation observed for the unbound wild-type protease in solution. Further analysis suggests that the wide-open structure observed for MDR 769 arises not from sequence variation, but instead is an artifact from crystal packing. Thus, despite being the first experimental structure to reveal flap opening sufficient for substrate access to the active site, this structure may not be directly relevant to studies of inhibitor entry or to the cause of HIV-PR drug resistance.     

Structure, dynamics and solvation of HIV-1 protease/saquinavir complex in aqueous solution and their contributions to drug resistance: molecular dynamic simulations

As it is known that the understanding of the basic properties of the enzyme/inhibitor complex leads directly to enhancing the capability in drug designing and drug discovery. Molecular dynamics simulations have been performed to examine detailed information on the structure and dynamical properties of the HIV-1 PR complexed with saquinavir in the three protonated states, monoprotonates at Asp25 (Mono-25) and Asp25'(Mono-25') and diprotonate (Di-Pro) at both Asp25 and Asp25'. The obtained results support clinical data which reveal that Ile84 and Gly48 are two of the most frequent residues where mutation toward a protease inhibitor takes place. In contrast to the Ile84 mutation due to high displacement of Ile84 in the presence of saquinavir, source of the Gly48 mutation was observed to be due to the limited space in the HIV-1 PR pocket. The Gly48 was, on one side, found to form strong hydrogen bonds with saquinavir, while on the other side this residue was repelled by the hydrophobic Phe53 residue. In terms of inhibitor/enzyme binding, interactions between saquinavir and a catalytic triad of the HIV-1 PR were calculated using the ab initio method. The results show an order of the binding energy of Mono25相似文献   

Resistant mechanism against nelfinavir of human immunodeficiency virus type 1 proteases

Inhibitors against human immunodeficiency virus type-1 (HIV-1) proteases are finely effective for anti-HIV-1 treatments. However, the therapeutic efficacy is reduced by the rapid emergence of inhibitor-resistant variants of the protease. Among patients who failed in the inhibitor nelfinavir (NFV) treatment, D30N, N88D, and L90M mutations of HIV-1 protease are often observed. Despite the serious clinical problem, it is not clear how these mutations, especially nonactive site mutations N88D and L90M, affect the affinity of NFV or why they cause the resistance to NFV. In this study, we executed molecular dynamics simulations of the NFV-bound proteases in the wild-type and D30N, N88D, D30N/N88D, and L90M mutants. Our simulations clarified the conformational change at the active site of the protease and the change of the affinity with NFV for all of these mutations, even though the 88th and 90th residues are not located in the NFV-bound cavity and not able to directly interact with NFV. D30N mutation causes the disappearance of the hydrogen bond between the m-phenol group of NFV and the 30th residue. N88D mutation alters the active site conformation slightly and induces a favorable hydrophobic contact. L90M mutation dramatically changes the conformation at the flap region and leads to an unfavorable distortion of the binding pocket of the protease, although 90M is largely far apart from the flap region. Furthermore, the changes of binding energies of the mutants from the wild-type protease are shown to be correlated with the mutant resistivity previously reported by the phenotypic experiments.     

Computational simulations of HIV-1 proteases--multi-drug resistance due to nonactive site mutation L90M

Human immunodeficiency virus type 1 protease (HIV-1 PR) is one of the proteins that currently available anti-HIV-1 drugs target. Inhibitors of HIV-1 PR have become available, and they have lowered the rate of mortality from acquired immune deficiency syndrome (AIDS) in advanced countries. However, the rate of emergence of drug-resistant HIV-1 variants is quite high because of their short retroviral life cycle and their high mutation rate. Serious drug-resistant mutations against HIV-1 PR inhibitors (PIs) frequently appear at the active site of PR. Exceptionally, some other mutations such as L90M cause drug resistance, although these appear at nonactive sites. The mechanism of resistance due to nonactive site mutations is difficult to explain. In this study, we carried out computational simulations of L90M PR in complex with each of three kinds of inhibitors and one typical substrate, and we clarified the mechanism of resistance. The L90M mutation causes changes in interaction between the side chain atoms of the 90th residue and the main chain atoms of the 25th residue, and a slight dislocation of the 25th residue causes rotation of the side chain at the 84th residue. The rotation of the 84th residue leads to displacement of the inhibitor from the appropriate binding location, resulting in a collision with the flap or loop region. The difference in levels of resistance to the three inhibitors has been explained from energetic and structural viewpoints, which provides the suggestion for promising drugs keeping its efficacy even for the L90M mutant.