Role of invariant Thr80 in human immunodeficiency virus type 1 protease structure, function, and viral infectivity

Sequence variability associated with human immunodeficiency virus type 1 (HIV-1) is useful for inferring structural and/or functional constraints at specific residues within the viral protease. Positions that are invariant even in the presence of drug selection define critically important residues for protease function. While the importance of conserved active-site residues is easily understood, the role of other invariant residues is not. This work focuses on invariant Thr80 at the apex of the P1 loop of HIV-1, HIV-2, and simian immunodeficiency virus protease. In a previous study, we postulated, on the basis of a molecular dynamics simulation of the unliganded protease, that Thr80 may play a role in the mobility of the flaps of protease. In the present study, both experimental and computational methods were used to study the role of Thr80 in HIV protease. Three protease variants (T80V, T80N, and T80S) were examined for changes in structure, dynamics, enzymatic activity, affinity for protease inhibitors, and viral infectivity. While all three variants were structurally similar to the wild type, only T80S was functionally similar. Both T80V and T80N had decreased the affinity for saquinavir. T80V significantly decreased the ability of the enzyme to cleave a peptide substrate but maintained infectivity, while T80N abolished both activity and viral infectivity. Additionally, T80N decreased the conformational flexibility of the flap region, as observed by simulations of molecular dynamics. Taken together, these data indicate that HIV-1 protease functions best when residue 80 is a small polar residue and that mutations to other amino acids significantly impair enzyme function, possibly by affecting the flexibility of the flap domain.     

HIV-1 protease molecular dynamics of a wild-type and of the V82F/I84V mutant: possible contributions to drug resistance and a potential new target site for drugs

The protease from type 1 human immunodeficiency virus (HIV-1) is a critical drug target against which many therapeutically useful inhibitors have been developed; however, the set of viral strains in the population has been shifting to become more drug-resistant. Because indirect effects are contributing to drug resistance, an examination of the dynamic structures of a wild-type and a mutant could be insightful. Consequently, this study examined structural properties sampled during 22 nsec, all atom molecular dynamics (MD) simulations (in explicit water) of both a wild-type and the drug-resistant V82F/I84V mutant of HIV-1 protease. The V82F/I84V mutation significantly decreases the binding affinity of all HIV-1 protease inhibitors currently used clinically. Simulations have shown that the curling of the tips of the active site flaps immediately results in flap opening. In the 22-nsec MD simulations presented here, more frequent and more rapid curling of the mutant's active site flap tips was observed. The mutant protease's flaps also opened farther than the wild-type's flaps did and displayed more flexibility. This suggests that the effect of the mutations on the equilibrium between the semiopen and closed conformations could be one aspect of the mechanism of drug resistance for this mutant. In addition, correlated fluctuations in the active site and periphery were noted that point to a possible binding site for allosteric inhibitors.     

Close correspondence between the motions from principal component analysis of multiple HIV-1 protease structures and elastic network modes

The large number of available HIV-1 protease structures provides a remarkable sampling of conformations of the different conformational states, which can be viewed as direct structural information about the dynamics of the HIV-1 protease. After structure matching, we apply principal component analysis (PCA) to obtain the important apparent motions for both bound and unbound structures. There are significant similarities between the first few key motions and the first few low-frequency normal modes calculated from a static representative structure with an elastic network model (ENM), strongly suggesting that the variations among the observed structures and the corresponding conformational changes are facilitated by the low-frequency, global motions intrinsic to the structure. Similarities are also found when the approach is applied to an NMR ensemble, as well as to molecular dynamics (MD) trajectories. Thus, a sufficiently large number of experimental structures can directly provide important information about protein dynamics, but ENM can also provide similar sampling of conformations.     

Studies on adaptability of binding residues and flap region of TMC-114 resistance HIV-1 protease mutants

Drug resistant mutations have severely restricted the success of HIV therapy. 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, which are able to produce significant resistance against the anti-HIV drug TMC-114. We provide insight into the molecular basis of TMC-114 resistance major flap mutations (I50V and I54M) in HIV-1 protease. It reports the shape complementarity and receptor-ligand interaction analysis supported by unrestrained all-atom molecular dynamics simulations of wild and major flap mutants of HIV-1 protease that sample large conformational changes of the flaps and active site binding residues. Both resistant flap mutants showed less atomic interaction toward TMC-114 and more structural deviation compared to wild HIV-protease. It is due to increasing flexibility at TMC-114 binding cavity and deviation of binding residues in 3-D space. Distortion in binding cavity and deviation in binding residues are the result of alteration in hydrogen bonding. Flap region also exhibited similar behaviour due to changes in number of hydrogen bonds during simulations.     

Studies on flexibility and binding affinity of Asp25 of HIV-1 protease mutants

We have investigated and highlighted the behavior of binding residue, Asp25 by computational analysis, which play an important role in understanding docking process with drug molecule, Ritonavir (Norvir^®) and the flexibility nature of the Human Immunodeficiency Virus-1 (HIV-1) protease enzyme. It is well known that Ritonavir is a potent and a selective HIV-1 protease inhibitor. Molecular dockings were performed in order to gain insights regarding the binding mode of this inhibitor. In our analysis, we observed Ritonavir had different rank orders of scores against different mutant of this enzyme. Asp25 of the enzyme was found to be the active site for all the mutants. The results clearly suggest that Ritonavir is not able to appropriately bind at the active site of each HIV-1 protease mutant due to RMSD difference of the amino acid (Asp) at the position 25 of all mutants. These findings support the concept that 3D space of active site is a qualitative assessment for binding affinity of inhibitor with an enzyme. The investigation on the flexibility nature of Asp25 by normal mode analysis, show that binding residue posses less flexibility due to its solvation potential. The overall analysis of our study brings clarity to the binding behavior with respect to the different mutants with Ritonavir on the basis RMSD and also on the flexible nature of HIV-1 protease enzyme with respect to Asp25 position.     

Atomistic simulations of the HIV-1 protease folding inhibition

Biochemical experiments have recently revealed that the p-S8 peptide, with an amino-acid sequence identical to the conserved fragment 83-93 (S8) of the HIV-1 protease, can inhibit catalytic activity of the enzyme by interfering with protease folding and dimerization. In this study, we introduce a hierarchical modeling approach for understanding the molecular basis of the HIV-1 protease folding inhibition. Coarse-grained molecular docking simulations of the flexible p-S8 peptide with the ensembles of HIV-1 protease monomers have revealed structurally different complexes of the p-S8 peptide, which can be formed by targeting the conserved segment 24-34 (S2) of the folding nucleus (folding inhibition) and by interacting with the antiparallel termini β-sheet region (dimerization inhibition). All-atom molecular dynamics simulations of the inhibitor complexes with the HIV-1 PR monomer have been independently carried out for the predicted folding and dimerization binding modes of the p-S8 peptide, confirming the thermodynamic stability of these complexes. Binding free-energy calculations of the p-S8 peptide and its active analogs are then performed using molecular dynamics trajectories of the peptide complexes with the HIV-1 PR monomers. The results of this study have provided a plausible molecular model for the inhibitor intervention with the HIV-1 PR folding and dimerization and have accurately reproduced the experimental inhibition profiles of the active folding inhibitors.     

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

Abstract

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     

Hydrophobic sliding: a possible mechanism for drug resistance in human immunodeficiency virus type 1 protease

Hydrophobic residues outside the active site of HIV-1 protease frequently mutate in patients undergoing protease inhibitor therapy; however, the mechanism by which these mutations confer drug resistance is not understood. From analysis of molecular dynamics simulations, 19 core hydrophobic residues appear to facilitate the conformational changes that occur in HIV-1 protease. The hydrophobic core residues slide by each other, exchanging one hydrophobic van der Waal contact for another, with little energy penalty, while maintaining many structurally important hydrogen bonds. Such hydrophobic sliding may represent a general mechanism by which proteins undergo conformational changes. Mutation of these residues in HIV-1 protease would alter the packing of the hydrophobic core, affecting the conformational flexibility of the protease. Therefore these residues impact the dynamic balance between processing substrates and binding inhibitors, and thus contribute to drug resistance.     

Cooperative fluctuations of unliganded and substrate-bound HIV-1 protease: a structure-based analysis on a variety of conformations from crystallography and molecular dynamics simulations

The dynamics of HIV-1 protease, both in unliganded and substrate-bound forms have been analyzed by using an analytical method, Gaussian network model (GNM). The method is applied to different conformations accessible to the protein backbone in the native state, observed in crystal structures and snapshots from fully atomistic molecular dynamics (MD) simulation trajectories. The modes of motion obtained from GNM on different conformations of HIV-1 protease are conserved throughout the MD simulations. The flaps and 40's loop of the unliganded HIV-1 protease structure are identified as the most mobile regions. However, in the liganded structure these flaps lose mobility, and terminal regions of the monomers become more flexible. Analysis of the fast modes shows that residues important for stability are in the same regions of all the structures examined. Among these, Gly86 appears to be a key residue for stability. The contribution of residues in the active site region and flaps to the stability is more pronounced in the substrate-bound form than in the unliganded form. The convergence of modes in GNM to similar regions of HIV-1 protease, regardless of the conformation of the protein, supports the robustness of GNM as a potentially useful and predictive tool.     

Comparative analysis of HIV-1 Tat variants

HIV-1 Tat protein is a crucial element for viral replication; therefore, its inhibition might be exploited against the AIDS infection. To gain insights on the natural variability of this protein, we present a comparative investigation on the relationship between the primary sequences and the experimentally available three-dimensional structures from the HIV-1 Tat variants Z2, BRU, and MAL. Our computational tools include sequence conservation algorithms, structural analysis, electrostatic modeling, and molecular dynamics (MD) simulations. We find that two regions located between residues 10-18 and 41-52 display the highest primary sequence conservation, while the most conserved region among the available structures corresponds approximately to the segment between positions approximately 44 and 50. Furthermore, in spite of their large structural divergence, Tat variants share a common mode for long-range intramolecular interactions. Finally, the flexibility of the Z2, BRU, and MAL variants, as emerging from multinanosecond MD simulations, is rather similar. Based on this work, we conclude that the turnlike region between amino acids 44 and 50 is structurally most conserved, emerging as an important motif for pharmaceutical targeting aimed toward inhibiting Tat action.     

Efficiency of a second-generation HIV-1 protease inhibitor studied by molecular dynamics and absolute binding free energy calculations

A subnanomolar inhibitor of human immunodeficiency virus type 1 (HIV-1) protease, designated QF34, potently inhibits the wild-type and drug-resistant enzyme. To explain its broad activity, the binding of QF34 to the wild-type HIV-1 protease is investigated by molecular dynamics simulations and compared to the binding of two inhibitors that are used clinically, saquinavir (SQV) and indinavir (IDV). Analysis of the flexibility of protease residues and inhibitor segments in the complex reveals that segments of QF34 were more mobile during the dynamics studies than the segments of SQV and IDV. The dynamics of hydrogen bonding show that QF34 forms a larger number of stable hydrogen bonds than the two inhibitors that are used clinically. Absolute binding free energies were calculated with molecular mechanics-generalized Born surface area (MM-GBSA) methodology using three protocols. The most consistent results were obtained using the single-trajectory approach, due to cancellation of errors and inadequate sampling in the separate-trajectory protocols. For all three inhibitors, energy components in favor of binding include van der Waals and electrostatic terms, whereas polar solvation and entropy terms oppose binding. Decomposition of binding energies reveals that more protease residues contribute significantly to the binding of QF34 than to the binding of SQV and IDV. Moreover, contributions from protease main chains and side chains are balanced in the case of QF34 (52:48 ratio, respectively), whereas side chain contributions prevail in both SQV and IDV (main-chain:side-chain ratios of 41:59 and 45:55, respectively). The presented results help explain the ability of QF34 to inhibit multiple resistant mutants and should be considered in the design of broad-specificity second-generation HIV-1 protease inhibitors.     

Influence of Structural Symmetry on Protein Dynamics

Structural symmetry in homooligomeric proteins has intrigued many researchers over the past several decades. However, the implication of protein symmetry is still not well understood. In this study, we performed molecular dynamics (MD) simulations of two forms of trp RNA binding attenuation protein (TRAP), the wild-type 11-mer and an engineered 12-mer, having two different levels of circular symmetry. The results of the simulations showed that the inter-subunit fluctuations in the 11-mer TRAP were significantly smaller than the fluctuations in the 12-mer TRAP while the internal fluctuations were larger in the 11-mer than in the 12-mer. These differences in thermal fluctuations were interpreted by normal mode analysis and group theory. For the 12-mer TRAP, the wave nodes of the normal modes existed at the flexible interface between the subunits, while the 11-mer TRAP had its nodes within the subunits. The principal components derived from the MD simulations showed similar mode structures. These results demonstrated that the structural symmetry was an important determinant of protein dynamics in circularly symmetric homooligomeric proteins.     

The effect of calciums on molecular motions of proteinase K

The native serine protease proteinase K binds two calcium cations. It has been reported that Ca²⁺ removal decreased the enzyme’s thermal stability and to some extent the substrate affinity, but has discrepant effects on catalytic activity of the enzyme. Molecular dynamics simulations were performed on the Ca²⁺-bound and Ca²⁺-free proteases to investigate the mechanism by which the calciums affect the structural stability, molecular motions, and catalytic activity of proteinase K. Very similar structural properties were observed between these two forms of proteinase K during simulations; and several long-lived hydrogen bonds and salt bridges common to both forms of proteinase K were found to be crucial in maintaining the local conformations around these two Ca²⁺ sites. Although Ca²⁺ removal enhanced the overall flexibility of proteinase K, the flexibility in a limited number of segments surrounding the substrate-binding pockets decreased. The largest differences in the equilibrium structures of the two simulations indicate that, upon the removal of Ca²⁺, the large concerted motion originating from the Ca1 site can transmit to the substrate-binding regions but not to the catalytic triad residues. In conjunction with the large overlap of the essential subspaces between the two simulations, these results not only provide insight into the dynamics of the underlying molecular mechanism responsible for the unchanged enzymatic activity as well as the decreased thermal stability and substrate affinity of proteinase K upon Ca²⁺ removal, but also complement the experimentally determined structural and biochemical data.     

A diverse view of protein dynamics from NMR studies of HIV-1 protease flaps

Internal motion in proteins fulfills a multitude of roles in biological processes. NMR spectroscopy has been applied to elucidate protein dynamics at the atomic level, albeit at a low resolution, and is often complemented by molecular dynamics simulation. However, it is critical to justify the consistency between simulation results and conclusions often drawn from multiple experiments in which uncertainties arising from assumed motional models may not be explicitly evaluated. To understand the role of the flaps of HIV-1 protease dimer in substrate recognition and protease function, many molecular dynamics simulations have been performed. The simulations have resulted in various proposed models of the flap dynamics, some of which are more consistent than others with our working model previously derived from experiments. However, using the working model to discriminate among the simulation results is not straightforward because the working model was derived from a combination of NMR experiments and crystal structure data. In this study, we use the NMR chemical shifts and relaxation data of the protease "monomer" rather than structural data to narrow down the possible conformations of the flaps of the "dimer". For the first time, we show that the tips of the flaps in the unliganded protease dimer interact with each other in solution. Accordingly, we discuss the consistency of the simulations with the model derived from all experimental data.     

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.     

Molecular dynamics studies of the wild-type and double mutant HIV-1 integrase complexed with the 5CITEP inhibitor: mechanism for inhibition and drug resistance

The human immunodeficiency virus type 1 (HIV-1) integrase (IN) is an essential enzyme in the life cycle of the virus and is an attractive target for the development of new drugs useful in acquired immunodeficiency syndrome multidrug therapy. Starting from the crystal structure of the 5CITEP inhibitor bound to the active site in the catalytic domain of the HIV-1 IN, two different molecular dynamics simulations in water have been carried out. In the first simulation the wild-type IN was used, whereas in the second one the double mutation T66I/M154I, described to lead to drug resistance, was introduced in the protein. Compelling differences have been observed in these two structures during analyses of the molecular dynamics trajectories, particularly in the inhibitor binding modes and in the conformational flexibility of the loop (residues 138-149) located near the three catalytic residues in the active site (Asp(64), Asp(116), Glu(152)). Because the conformational flexibility of this region is important for efficient biological activity and its behavior is quite different in the two models, we suggest a hypothetical mechanism for the inhibition and drug resistance of HIV-1 IN. These results can be useful for the rational design of more potent and selective integrase inhibitors and may allow for the design of inhibitors that will be more robust against known resistance mutations.     

Theoretical studies of relaxation of a monomeric subunit of HIV-1 protease in water using molecular dynamics

The dynamic behavior of one 99-residue subunit of the dimeric aspartyl protease of HIV-1 was studied in a 160 psec molecular dynamics simulation at 300 K in water. The crystal structure of one of the identical subunits of the dimer was the starting point, with the aqueous phase modeled by 4,331 explicit waters in a restrained spherical droplet Analysis of the simulations showed that the monomer displayed considerable flexibility in the interfacial portions of the flap (the region which folds over the substrate), the N- and C-0termini, and, to a lesser extent, the active site. The flap undergoes significant motion as an independent rigid finger, but without the cantilever previously reported hi a simulation of the dimer. The N-terminus displayed the greatest fluctuational disorder whereas the C-terminus exhibited the greatest root mean square movement from the crystal structure. The central core of the monomer had a heavy-atom root mean square deviation from the initial structure of about 3.0 Å during the latter half of the simulation. Although this is larger than the 1.6 Å found for comparable simulations of typical globular proteins, the general features of the tertiary structure were preserved over the course of the simulation. Overall, these results indicate that the relaxed structure obtained in these simulations may provide a better model for the tertiary structure of the solvated HIV-1 protease monomer than the subunit conformation seen in the X-ray crystallographic structure of the dimer. Except in the flap region, the design of compounds intended to interfere with dimerization should take this relaxation and the flexibility of the solvated monomer, especially at the termini, into account. © 1993 Wiley-Liss, Inc.