Molecular dynamics of HIV-1 protease.

Molecular dynamics simulations have been carried out based on the GROMOS force field on the aspartyl protease (PR) of the human immunodeficiency virus HIV-1. The principal simulation treats the HIV-1 PR dimer and 6990 water molecules in a hexagonal prism cell under periodic boundary conditions and was carried out for a trajectory of 100 psec. Corresponding in vacuo simulations, i.e., treating the isolated protein without solvent, were carried out to study the influence of solvent on the simulation. The results indicate that including waters explicitly in the simulation results in a model considerably closer to the crystal structure than when solvent is neglected. Detailed conformational and helicoidal analysis was performed on the solvated form to determine the exact nature of the dynamical model and the exact points of agreement and disagreement with the crystal structure. The calculated dynamical model was further elucidated by means of studies of the time evolution of the cross-correlation coefficients for atomic displacements of the atoms comprising the protein backbone. The cross-correlation analysis revealed significant aspects of structure originating uniquely in the dynamical motions of the molecule. In particular, an unanticipated through-space, domain-domain correlation was found between the mobile flap region covering the active site and a remote regions of the structure, which collectively act somewhat like a molecular cantilever. The significance of these results is discussed with respect to the inactivation of the protease by site-specific mutagenesis, and in the design of inhibitors.     

Molecular dynamics studies on HIV-1 protease: a comparison of the flap motions between wild type protease and the M46I/G51D double mutant

The emergence of drug-resistant mutants of HIV-1 is a tragic effect associated with conventional long-treatment therapies against acquired immunodeficiency syndrome. These mutations frequently involve the aspartic protease encoded by the virus; knowledge of the molecular mechanisms underlying the conformational changes of HIV-1 protease mutants may be useful in developing more effective and longer lasting treatment regimes. The flap regions of the protease are the target of a particular type of mutations occurring far from the active site. These mutations modify the affinity for both substrate and ligands, thus conferring resistance. In this work, molecular dynamics simulations were performed on a native wild type HIV-1 protease and on the drug-resistant M46I/G51D double mutant. The simulation was carried out for a time of 3.5 ns using the GROMOS96 force field, with implementation of the SPC216 explicit solvation model. The results show that the flaps may exist in an ensemble of conformations between a “closed” and an “open” conformation. The behaviour of the flap tips during simulations is different between the native enzyme and the mutant. The mutation pattern leads to stabilization of the flaps in a semi-open configuration.     

HIV-1 protease substrate binding and product release pathways explored with coarse-grained molecular dynamics

We analyze the encounter of a peptide substrate with the native HIV-1 protease, the mechanism of substrate incorporation in the binding cleft, and the dissociation of products after substrate hydrolysis. To account for the substrate, we extend a coarse-grained model force field, which we previously developed to study the flap opening dynamics of HIV-1 protease on a microsecond timescale. Molecular and Langevin dynamics simulations show that the flaps need to open for the peptide to bind and that the protease interaction with the substrate influences the flap opening frequency and interval. On the other hand, release of the products does not require flap opening because they can slide out from the binding cleft to the sides of the enzyme. Our data show that in the protease-substrate complex the highest fluctuations correspond to the 17- and 39-turns and the substrate motion is anticorrelated with the 39-turn. Moreover, the active site residues and the flap tips move in phase with the peptide. We suggest some mechanistic principles for how the flexibility of the protein may be involved in ligand binding and release.     

Purification of HIV-1 wild-type protease and characterization of proteolytically inactive HIV-1 protease mutants by pepstatin A affinity chromatography

Recombinant wild-type protease of human immunodeficiency virus, type [(HIV-1) expressed in E. coli was purified by pepstatin A affinity chromatography. An 88-fold purification was achieved giving a protease preparation with a specific enzymatic activity of approximately 3700 pmol/min/μg. Two proteolytically inactive HIV-1 mutant proteases (Arg-87 → Lys; Asn-88 → Glu) were found to bind to pepstatin A agarose, and they were purified as the wild-type protease. A third mutant protease (Arg-87 → Glu) was apparently unable to bind to pepstatin A under similar conditions. Binding to pepstatin A indicates the binding ability of the substrate binding site and the ability to form dimers. These features may be used to purify and to characterize other mutated HIV-1 proteases.     

Comprehensive mutagenesis of HIV-1 protease: a computational geometry approach

A computational geometry technique based on Delaunay tessellation of protein structure, represented by C(alpha) atoms, is used to study effects of single residue mutations on sequence-structure compatibility in HIV-1 protease. Profiles of residue scores derived from the four-body statistical potential are constructed for all 1881 mutants of the HIV-1 protease monomer and compared with the profile of the wild-type protein. The profiles for an isolated monomer of HIV-1 protease and the identical monomer in a dimeric state with an inhibitor are analyzed to elucidate changes to structural stability. Protease residues shown to undergo the greatest impact are those forming the dimer interface and flap region, as well as those known to be involved in inhibitor binding.     

Role of conformational fluctuations in the enzymatic reaction of HIV-1 protease

The emergence of compensatory drug-resistant mutations in HIV-1 protease challenges the common view of the reaction mechanism of this enzyme. Here, we address this issue by performing classical and ab initio molecular dynamics simulations (MD) on a complex between the enzyme and a peptide substrate. The classical MD calculation reveals large-scale protein motions involving the flaps and the cantilever. These motions modulate the conformational properties of the substrate at the cleavage site. The ab initio calculations show in turn that substrate motion modulates the activation free energy barrier of the enzymatic reaction dramatically. Thus, the catalytic power of the enzyme does not arise from the presence of a pre-organized active site but from the protein mechanical fluctuations. The implications of this finding for the emergence of drug-resistance are discussed.     

Analysis of the structure of HIV-1 protease complexed with a hexapeptide inhibitor. Part II: Molecular dynamic studies of the active site region

Six models of the catalytic site of HIV-1 protease complexed with a reduced peptide inhibitor, MVT-101, were investigated. These studies focused on the details of protonation of the active site, its total net charge and hydrogen bonding pattern, which was consistent with both the observed coplanar configuration of the acidic groups of the catalytic aspartates (Asp-25 and Asp-125) and the observed binding mode of the inhibitor. Molecular dynamic simulations using AMBER 4.0 indicated that the active site should be neutral. The planarity of the aspartate dyad may be due to the formation of two hydrogen bonds: one between the inner O_δ1oxygen atoms of the two catalytic aspartates and another between the O_δ2atom of Asp-125 and the nitrogen atom of the reduced peptide bond of the bound inhibitor. This would require two additional protonations, either of both aspartates, or of one Asp and the amido nitrogen atom of Nle-204. Our results favor the Asp-inhibitor protonation but the other one is not excluded. Implications of these findings for the mechanism of enzymatic catalysis are discussed. Dynamic properties of the hydrogen bond network in the active site and an analysis of the interaction energy between the inhibitor and the protease are presented. © 1997 Wiley-Liss, Inc.     

Molecular dynamics simulations on HIV-1 Tat

Molecular dynamics simulations are used to investigate dynamics and intramolecular interactions of the HIV-1 transactivator (Tat) in aqueous solution. The calculations are based on the AMBER force field with particle mesh Ewald treatment for long-range electrostatics. The Tat structure exhibits a large flexibility, consistent with its absence of secondary structure elements. From an analysis of the correlation matrix and of electrostatic interactions we suggest that segments expressed by the two exons (amino acids 1-72 and 73-86, respectively) exhibit rather separated dynamic and energetic properties. We also identify intramolecular interactions of importance for structure stabilization. In particular, significant electrostatic interactions are recognized between the N-terminus and the basic domain of the protein, consistent with site-directed mutagenesis performed in this work.     

Insights into the functional role of protonation states in the HIV-1 protease-BEA369 complex: molecular dynamics simulations and free energy calculations

The molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method combined with molecular dynamics (MD) simulations were used to investigate the functional role of protonation in human immunodeficiency virus type 1 (HIV-1) protease complexed with the inhibitor BEA369. Our results demonstrate that protonation of two aspartic acids (Asp25/Asp25′) has a strong influence on the dynamics behavior of the complex, the binding free energy of BEA369, and inhibitor–residue interactions. Relative binding free energies calculated using the MM-PBSA method show that protonation of Asp25 results in the strongest binding of BEA369 to HIV-1 protease. Inhibitor–residue interactions computed by the theory of free energy decomposition also indicate that protonation of Asp25 has the most favorable effect on binding of BEA369. In addition, hydrogen-bond analysis based on the trajectories of the MD simulations shows that protonation of Asp25 strongly influences the water-mediated link of a conserved water molecule, Wat301. We expect that the results of this study will contribute significantly to binding calculations for BEA369, and to the design of high affinity inhibitors.     

Relation between sequence and structure of HIV-1 protease inhibitor complexes: a model system for the analysis of protein flexibility.

The flexibility of different regions of HIV-1 protease was examined by using a database consisting of 73 X-ray structures that differ in terms of sequence, ligands or both. The root-mean-square differences of the backbone for the set of structures were shown to have the same variation with residue number as those obtained from molecular dynamics simulations, normal mode analyses and X-ray B-factors. This supports the idea that observed structural changes provide a measure of the inherent flexibility of the protein, although specific interactions between the protease and the ligand play a secondary role. The results suggest that the potential energy surface of the HIV-1 protease is characterized by many local minima with small energetic differences, some of which are sampled by the different X-ray structures of the HIV-1 protease complexes. Interdomain correlated motions were calculated from the structural fluctuations and the results were also in agreement with molecular dynamics simulations and normal mode analyses. Implications of the results for the drug-resistance engendered by mutations are discussed briefly.     

HIV-1 protease inhibition potential of functionalized polyoxometalates

Polyoxometalates (POMs) are interesting biomedical agents due to their versatile anticancer and antiviral properties, such as remarkable anti-HIV activity. Although POMs are tunable and easily accessible inorganic drug prototypes in principle, their full potential can only be tapped by enhancing their biocompatibility, for example, through organic functionalization. We have therefore investigated the HIV-1 protease inhibition potential of functionalized Keggin- and Dawson-type POMs with organic side chains. Their inhibitory performance was furthermore compared to other POM types, and the buffer dependence of the results is discussed. In addition, chemical shift mapping NMR experiments were performed to exclude POM-substrate interactions. Whereas the introduction of organic side chains into POMs is a promising approach in principle, the influence of secondary effects on the reaction system also merits detailed investigation.     

Fighting a virus with a virus: a dynamic model for HIV-1 therapy

A mathematical model examined a potential therapy for controlling viral infections using genetically modified viruses. The control of the infection is an indirect effect of the selective elimination by an engineered virus of infected cells that are the source of the pathogens. Therefore, this engineered virus could greatly compensate for a dysfunctional immune system compromised by AIDS. In vitro studies using engineered viruses have been shown to decrease the HIV-1 load about 1000-fold. However, the efficacy of this potential treatment for reducing the viral load in AIDS patients is unknown. The present model studied the interactions among the HIV-1 virus, its main host cell (activated CD4+ T cells), and a therapeutic engineered virus in an in vivo context; and it examined the conditions for controlling the pathogen. This model predicted a significant drop in the HIV-1 load, but the treatment does not eradicate HIV. A basic estimation using a currently engineered virus indicated an HIV-1 load reduction of 92% and a recovery of host cells to 17% of their normal level. Greater success (98% HIV reduction, 44% host cells recovery) is expected as more competent engineered viruses are designed. These results suggest that therapy using viruses could be an alternative to extend the survival of AIDS patients.     

High-performance liquid chromatographic assay of a new HIV-1 protease inhibitor, LB71350, in the plasma of dogs

A reliable reversed-phase high-performance liquid chromatographic method has been developed for the determination of LB71350 in the plasma of dogs. The analyte was deproteinized with 1.5 volumes of methanol and 0.5 volumes of 10% zinc sulfate, and the supernatant was injected into a 5-μm Capcell Pak C₁₈ column (150×4.6 mm I.D.). The mobile phase was a stepwise gradient mixture of acetonitrile and 0.2% triethylamine–HCl with a flow-rate of 1 ml/min and detection at UV 245 nm. The proportion of acetonitrile was kept at 52% for the first 6 min, increased to 100% for the next 0.5 min, kept at 100% for the next 2 min, decreased to 52% for the next 0.5 min, and finally kept at 52% for the next 7 min. The retention time of LB71350 was 6.9 min. The calibration was linear over the concentration range of 0.1–100 mg/l for dog plasma (r>0.997) and the limit of quantitation was 0.1 mg/l using 0.1 ml plasma. The quality control samples were reproducible with acceptable accuracy and precision at 0.1, 1, 10 and 100 mg/l concentrations. The within-day recovery (n=5) was 90.2–93.9%, the between-day recovery (n=5) was 89.5–93.5%, and the absolute between-day recovery (n=5) was 77–81%. The within-day precision (n=5) and between-day precision (n=5) were 2.59–5.82% and 3.17–4.55%, respectively. No interferences from endogenous substances were observed. Taken together, the above HPLC assay method by deproteinization and UV detection was suitable for the determination of LB71350 in the preclinical pharmacokinetics.     

Molecular dynamics studies of the transmembrane domain of gp41 from HIV-1

Helix-helix interactions in the putative three-helix bundle formation of the gp41 transmembrane (TM) domain may contribute to the process of virus-cell membrane fusion in HIV-1 infection. In this study, molecular dynamics is used to analyze and compare the conformations of monomeric and trimeric forms of the TM domain in various solvent systems over the course of 4 to 23-ns simulations. The trimeric bundles of the TM domain were stable as helices and remained associated in a hydrated POPE lipid bilayer for the duration of the 23-ns simulation. Several stable inter-chain hydrogen bonds, mostly among the three deprotonated arginine residues located at the center of each of the three TM domains, formed in a right-handed bundle embedded in the lipid bilayer. No such bonds were observed when the bundle was left-handed or when the central arginine residue in each of the three TM helices was replaced with isoleucine (R_I mutant), suggesting that the central arginine residues may play an essential role in maintaining the integrity of the three-helix bundle. These observations suggest that formation of the three-helix bundle of the TM domain may play a role in the trimerization of gp41, thought to occur during the virus-cell membrane fusion process.     

HIV-1 dynamics at different time scales under antiretroviral therapy

We exploit a model that considers three compartments: blood plasma (BP), lymphoid tissue-interstitial spaces (LT-IS), and follicular dendritic cells (FDC), for the HIV-1 dynamics under the application of highly active antiretroviral therapy (HAART) which allowed us to unravel distinct viral dynamics occurring in short- (2 days), middle- (21 days), and long-term (183 days) time scales. The different time scales are determined by the viral clearance rate, the ratio of productively infected CD4⁺ T cells to chronically infected cells, and the dissociation rate of HIV-1 complexes from FDC. This generates a scenario in which, after an initial transient stage, the viral BP dynamics decouples and becomes governed by the lymphoid tissue (LT) dynamics; in a later stage, a new decoupling occurs in which the LT-IS dynamics is slaved to that of the FDC dynamics. We observed an initial increase in the viremia after HAART in a patient who did not receive protease inhibitors (PI). By means of the above-mentioned model we were able to highlight the relevant parameters which need to be estimated at three different time scales after HAART.     

Molecular dynamics simulations applied to the study of subtypes of HIV-1 protease common to Brazil, Africa, and Asia

Africa accounts for the majority of HIV-1 infections worldwide caused mainly by the A and C viral subtypes rather than B subtype, which prevails in the United States and Western Europe. In Brazil, B subtype is the major subtype, but F, C, and A also circulate. These non-B subtypes present polymorphisms, and some of them occur at sites that have been associated with drug resistance, including the HIV-1 protease (PR), one important drug target. Here, we report a Molecular Dynamics study of the B and non-B PR complexed with the inhibitor ritonavir to delineate the behavior of each subtype. We compare root mean squared deviation, binding free energy by linear interaction energy approach, hydrogen bonds, and intermolecular contact surface area between inhibitor and PR. From our results, we can provide a basis to understand the molecular mechanism of drug resistance in non-B subtypes. In this sense, we found a decrease of approx 4 kcal/mol in ΔG of binding between B and non-B subtypes. This corresponds to the loss of one hydrogen bond, which is in agreement with our H-bond analysis. Previous experimental affinity studies reported analogous results with inhibition constant values for non-B PR.     

Structural and Dynamical Properties of a Full-length HIV-1 Integrase: Molecular Dynamics Simulations

Abstract

The structural and dynamical properties of the complete full-length structure of HIV-1 integrase were investigated using Molecular Dynamics approach. Simulations were carried out for the three systems, core domain only (CORE), full-length structure without (FULL) and with a Mg²⁺ (FULL+ION) in its active site, aimed to investigate the difference in the molecular properties of the full-length models due to their different construction procedures as well as the effects of the two ends, C- and N-terminal, on those properties in the core domain. The full-length structure was prepared from the two experimental structures of two-domain fragment. The following properties were observed to differ significantly from the previous reports: (i) relative topology formed by an angle between the three domains; (ii) the cavity size defined by the catalytic triad, Asp64, Asp116, and Glul52; (iii) distances and solvation of the Mg²⁺; and (iv) conformation of the catalytic residues. In addition, the presence of the two terminal domains decreases the mobility of the central core domain significantly.     

The L33F darunavir resistance mutation acts as a molecular anchor reducing the flexibility of the HIV-1 protease 30s and 80s loops

HIV-1 protease (PR) is a 99 amino acid protein responsible for proteolytic processing of the viral polyprotein – an essential step in the HIV-1 life cycle. Drug resistance mutations in PR that are selected during antiretroviral therapy lead to reduced efficacy of protease inhibitors (PI) including darunavir (DRV). To identify the structural mechanisms associated with the DRV resistance mutation L33F, we performed X-ray crystallographic studies with a multi-drug resistant HIV-1 protease isolate that contains the L33F mutation (MDR769 L33F). In contrast to other PR L33F DRV complexes, the structure of MDR769 L33F complexed with DRV reported here displays the protease flaps in an open conformation. The L33F mutation increases noncovalent interactions in the hydrophobic pocket of the PR compared to the wild-type (WT) structure. As a result, L33F appears to act as a molecular anchor, reducing the flexibility of the 30s loop (residues 29–35) and the 80s loop (residues 79–84). Molecular anchoring of the 30s and 80s loops leaves an open S1/S1′ subsite and distorts the conserved hydrogen-bonding network of DRV. These findings are consistent with previous reports despite structural differences with regards to flap conformation.