Curling of flap tips in HIV-1 protease as a mechanism for substrate entry and tolerance of drug resistance

BACKGROUND: The human immunodeficiency virus type 1 (HIV-1) protease is an essential viral protein that is a major drug target in the fight against Acquired Immune Deficiency Syndrome (AIDS). Access to the active site of this homodimeric enzyme is gained when two large flaps, one from each monomer, open. The flap movements are therefore central to the function of the enzyme, yet determining how these flaps move at an atomic level has not been experimentally possible. RESULTS: In the present study, we observe the flaps of HIV-1 protease completely opening during a 10 ns solvated molecular dynamics simulation starting from the unliganded crystal structure. This movement is on the time scale observed by Nuclear Magnetic Resonance (NMR) relaxation data. The highly flexible tips of the flaps, with the sequence Gly-Gly-Ile-Gly-Gly, are seen curling back into the protein and thereby burying many hydrophobic residues. CONCLUSIONS: This curled-in conformational change has never been previously described. Previous models of this movement, with the flaps as rigid levers, are not consistent with the experimental data. The residues that participate in this hydrophobic cluster as a result of the conformational change are highly sensitive to mutation and often contribute to drug resistance when they do change. However, several of these residues are not part of the active site cavity, and their essential role in causing drug resistance could possibly be rationalized if this conformational change actually occurs. Trapping HIV-1 protease in this inactive conformation would provide a unique opportunity for future drug design.     

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.     

Characterization of two hydrophobic methyl clusters in HIV-1 protease by NMR spin relaxation in solution

Nearly 50 % of the amino acid residues of HIV-1 protease contain methyl side-chains, most of which appear to be organized into two clusters: the inner cluster that nearly surrounds the active site and the outer cluster that contains the hydrophobic core which stabilizes the inhibitor-free protease structure. NMR relaxation experiments sensitive to motions of methyl groups on the sub-nanosecond and the milli-microsecond time-scales revealed flexible methyl groups in residues that link the two clusters, the methyl groups of L10, L23, V75, and L76. We hypothesize that flexibility at the junctions of these clusters allows the protease to minimize conformational changes upon drug-binding. The two-methyl cluster motif appears to be a common structural feature among retroviral proteases and may play a similar role throughout this family of enzymes.     

Mutation patterns and structural correlates in human immunodeficiency virus type 1 protease following different protease inhibitor treatments

Although many human immunodeficiency virus type 1 (HIV-1)-infected persons are treated with multiple protease inhibitors in combination or in succession, mutation patterns of protease isolates from these persons have not been characterized. We collected and analyzed 2,244 subtype B HIV-1 isolates from 1,919 persons with different protease inhibitor experiences: 1,004 isolates from untreated persons, 637 isolates from persons who received one protease inhibitor, and 603 isolates from persons receiving two or more protease inhibitors. The median number of protease mutations per isolate increased from 4 in untreated persons to 12 in persons who had received four or more protease inhibitors. Mutations at 45 of the 99 amino acid positions in the protease-including 22 not previously associated with drug resistance-were significantly associated with protease inhibitor treatment. Mutations at 17 of the remaining 99 positions were polymorphic but not associated with drug treatment. Pairs and clusters of correlated (covarying) mutations were significantly more likely to occur in treated than in untreated persons: 115 versus 23 pairs and 30 versus 2 clusters, respectively. Of the 115 statistically significant pairs of covarying residues in the treated isolates, 59 were within 8 A of each other-many more than would be expected by chance. In summary, nearly one-half of HIV-1 protease positions are under selective drug pressure, including many residues not previously associated with drug resistance. Structural factors appear to be responsible for the high frequency of covariation among many of the protease residues. The presence of mutational clusters provides insight into the complex mutational patterns required for HIV-1 protease inhibitor resistance.     

Conserved hydrogen bonds and water molecules in MDR HIV-1 protease substrate complexes

The success of highly active antiretroviral therapy (HAART) in anti-HIV therapy is severely compromised by the rapidly developing drug resistance. HIV-1 protease inhibitors, part of HAART, are losing their potency and efficacy in inhibiting the target. Multi-drug resistant (MDR) 769 HIV-1 protease (resistant mutations at residues 10, 36, 46, 54, 62, 63, 71, 82, 84, 90) was selected for the present study to understand the binding to its natural substrates. The nine crystal structures of MDR769 HIV-1 protease substrate hepta-peptide complexes were analyzed in order to reveal the conserved structural elements for the purpose of drug design against MDR HIV-1 protease. Our structural studies demonstrated that highly conserved hydrogen bonds between the protease and substrate peptides, together with the conserved crystallographic water molecules, played a crucial role in the substrate recognition, substrate stabilization and protease stabilization. In addition, the absence of the key flap-ligand bridging water molecule might imply a different catalytic mechanism of MDR769 HIV-1 protease compared to that of wild type (WT) HIV-1 protease.     

Multiple routes and milestones in the folding of HIV-1 protease monomer

Proteins fold on a time scale incompatible with a mechanism of random search in conformational space thus indicating that somehow they are guided to the native state through a funneled energetic landscape. At the same time the heterogeneous kinetics suggests the existence of several different folding routes. Here we propose a scenario for the folding mechanism of the monomer of HIV-1 protease in which multiple pathways and milestone events coexist. A variety of computational approaches supports this picture. These include very long all-atom molecular dynamics simulations in explicit solvent, an analysis of the network of clusters found in multiple high-temperature unfolding simulations and a complete characterization of free-energy surfaces carried out using a structure-based potential at atomistic resolution and a combination of metadynamics and parallel tempering. Our results confirm that the monomer in solution is stable toward unfolding and show that at least two unfolding pathways exist. In our scenario, the formation of a hydrophobic core is a milestone in the folding process which must occur along all the routes that lead this protein towards its native state. Furthermore, the ensemble of folding pathways proposed here substantiates a rational drug design strategy based on inhibiting the folding of HIV-1 protease.     

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.     

Why is Substrate Peptide Binding Unsusceptible to Multidrug-Resistant Mutations in HIV-1 Protease? A Structural and Energetic Analysis

Understanding of the molecular mechanism and biological implication underlying the difference in binding of substrate peptides and small-molecule inhibitors to multidrug-resistant mutants of HIV-1 protease would help to develop new anti-HIV agents combating drug resistance. Here, an integration of rigorous quantum mechanics/molecular mechanics (QM/MM) analysis and empirical Poisson–Boltzmann/surface area (PB/SA) model is described to investigate the structural basis and energetic property of wild-type HIV-1 protease and its mutants in recognizing and binding with a wide variety of ligands, including the peptides derived from its cognate cleavage sites and the cleavage site variants as well as a number of FDA-approved protease inhibitors, attempting to explain why is substrate binding unsusceptible to most observed HIV-1 protease mutations. A preliminary test study demonstrates that the combined QM/MM–PB/SA scheme is able to effectively reproduce the relative ligand binding energy changes upon protease single- and double-mutations, albeit the absolute values appear to be different significantly between the calculated and experimental results. With the QM/MM–PB/SA calculations a complete mutation energy map of HIV-1 protease–ligand interactions is created, which unravels distinct affinity pictures of wild-type substrates, substrate variants and, particularly, the protease inhibitors bound to HIV-1 protease mutants, suggesting that, on the one hand, the evaluation pressure under anti-HIV chemotherapies addresses site-directed protease mutations that impair and undermine the intermolecular interactions specific to inhibitors but not substrates; on the other hand, co-evaluation of protease and its substrate peptides provides a more effective mechanism to avoid therapeutic surveillance. Further, nonbonded interaction analysis and computational alanine scanning reveal 12 key residues that is critical for substrate binding, from which the Asn25, Gly27, Ala28, Asp29 and Pro81 are identified that have not yet been found to cause drug resistance and hence would be the promising sites targeted by new protease inhibitors.     

Closing of the flaps of HIV-1 protease induced by substrate binding: a model of a flap closing mechanism in retroviral aspartic proteases

The active site of aspartic proteases is covered by one or more flaps, which control access to the active site and play a significant role in the binding of the substrate. An extensive conformational change of the flaps takes place upon binding of substrate to the active site. A long molecular dynamics simulation was performed on the complex consisting of a peptide (CA-p2) from a natural substrate cleavage site of the gag/pol polyprotein placed in the active site of HIV-1 protease (PR) with an open flap conformation. During the simulation, the substrate induced the closing of the flaps into the closed conformation in an asymmetrical way through a hydrophobic intermediate state cluster. The nature of the residues of HIV-1 PR identified to be important in the flap closing mechanism is conserved across known structures of retroviral aspartic proteases family. The flap closing mechanism described in HIV-1 PR is proposed to be a general model for flap closing in retroviral aspartic proteases.     

Insights into saquinavir resistance in the G48V HIV-1 protease: quantum calculations and molecular dynamic simulations

The spread of acquired immune deficiency syndrome has increasingly become a great concern owing largely to the failure of chemotherapies. The G48V is considered the key signature residue mutation of HIV-1 protease developing with saquinavir therapy. Molecular dynamics simulations of the wild-type and the G48V HIV-1 protease complexed with saquinavir were carried out to explore structure and interactions of the drug resistance. The molecular dynamics results combined with the quantum-based and molecular mechanics Poisson-Boltzmann surface area calculations indicated a monoprotonation took place on D25, one of the triad active site residues. The inhibitor binding of the triad residues and its interaction energy in the mutant were similar to those in the wild-type. The overall structure of both complexes is almost identical. However, the steric conflict of the substituted valine results in the conformational change of the P2 subsite and the disruption of hydrogen bonding between the -NH of the P2 subsite and the backbone -CO of the mutated residue. The magnitude of interaction energy changes was comparable to the experimental K(i) data. The designing for a new drug should consider a reduction of steric repulsion on P2 to enhance the activity toward this mutant strain.     

Predicting drug resistance of the HIV-1 protease using molecular interaction energy components

Drug resistance significantly impairs the efficacy of AIDS therapy. Therefore, precise prediction of resistant viral mutants is particularly useful for developing effective drugs and designing therapeutic regimen. In this study, we applied a structure-based computational approach to predict mutants of the HIV-1 protease resistant to the seven FDA approved drugs. We analyzed the energetic pattern of the protease-drug interaction by calculating the molecular interaction energy components (MIECs) between the drug and the protease residues. Support vector machines (SVMs) were trained on MIECs to classify protease mutants into resistant and nonresistant categories. The high prediction accuracies for the test sets of cross-validations suggested that the MIECs successfully characterized the interaction interface between drugs and the HIV-1 protease. We conducted a proof-of-concept study on a newly approved drug, darunavir (TMC114), on which no drug resistance data were available in the public domain. Compared with amprenavir, our analysis suggested that darunavir might be more potent to combat drug resistance. To quantitatively estimate binding affinities of drugs and study the contributions of protease residues to causing resistance, linear regression models were trained on MIECs using partial least squares (PLS). The MIEC-PLS models also achieved satisfactory prediction accuracy. Analysis of the fitting coefficients of MIECs in the regression model revealed the important resistance mutations and shed light into understanding the mechanisms of these mutations to cause resistance. Our study demonstrated the advantages of characterizing the protease-drug interaction using MIECs. We believe that MIEC-SVM and MIEC-PLS can help design new agents or combination of therapeutic regimens to counter HIV-1 protease resistant strains.     

Mechanism of action of aspartic proteinases. V. Conformational characteristics of fragments of substrate-binding sites in rhizopuspepsin and HIV-1 proteinase

The conformational states of side chains of catalytic Asp residues in active sites of HIV-1 protease and rhizopuspepsin in the potential field of free enzymes were studied by using theoretical conformational analysis. Structural factors that stabilize the conformation of these residues in free enzymes were revealed. Methods of molecular mechanics were used to estimate the stabilization energy of the Met46-Phe53 labile fragments of HIV-1 protease in the potential field of their nearest surrounding amino acid residues for the conformations characteristic of the free protein and similar to that of the protein in enzyme-inhibitor complexes. In solution, the conformational state of the fragments of the free enzyme was concluded to be similar to that observed in the enzyme complex with the ligand and different from that determined by X-ray diffraction analysis. This difference was ascribed to the effect of crystal packing.     

Genetic evidence that interhelical packing interactions in the gp41 core are critical for transition of the human immunodeficiency virus type 1 envelope glycoprotein to the fusion-active state

The envelope glycoprotein complex (gp120-gp41) of human immunodeficiency virus type 1 (HIV-1) promotes the fusion of viral and cellular membranes through formation of the fusion-active six-helix bundle in the gp41 ectodomain. This gp41 core structure consists of three C-terminal helices packed in an antiparallel manner into hydrophobic grooves on the surface of the N-terminal trimeric coiled coil. Alanine mutations that destabilize the N- and C-terminal interhelical packing interactions also reduce viral infectivity. Here we show that viruses bearing these mutations exhibit a marked potentiation of inhibition by peptides that make up the gp41 core. By contrast, these viruses are unchanged in their sensitivities to soluble CD4, the CXCR4 coreceptor ligand SDF-1alpha, and human anti-HIV immunoglobulin, reagents that impact the initial, receptor-induced conformational changes in the envelope glycoprotein. Our results support the notion that these alanine mutations specifically affect the conformational transition to the fusion-active gp41 structure. The mutations also increase viral sensitivity to the gp41-directed monoclonal antibody 2F5, suggesting that this broadly neutralizing antibody may also interfere with this transition. The conformational activation of the HIV-1 envelope glycoprotein likely represents a viable target for vaccine and antiviral drug development.     

Comparative analysis of the sequences and structures of HIV-1 and HIV-2 proteases

The different isolates available for HIV-1 and HIV-2 were compared for the region of the protease (PR) sequence, and the variations in amino acids were analyzed with respect to the crystal structure of HIV-1 PR with inhibitor. Based on the extensive homology (39 identical out of 99 residues), models were built of the HIV-2 PR complexed with two different aspartic protease inhibitors, acetylpepstatin and a renin inhibitor, H-261. Comparison of the HIV-1 PR crystal structure and the HIV-2 PR model structure and the analysis of the changes found in different isolates showed that correlated substitutions occur in the hydrophobic interior of the molecule and at surface residues involved in ionic or hydrogen bond interactions. The substrate binding residues of HIV-1 and HIV-2 PRs show conservative substitutions of four residues. The difference in affinity of HIV-1 and HIV-2 PRs for the two inhibitors appears to be due in part to the change of Val 32 in HIV-1 PR to Ile in HIV-2 PR.     

Small-molecule inhibition of human immunodeficiency virus type 1 replication by specific targeting of the final step of virion maturation

Despite the effectiveness of currently available human immunodeficiency virus type 1 (HIV-1) therapies, a continuing need exists for new drugs to treat HIV-1 infection. We investigated the mechanism by which 3-O-[3',3'-dimethylsuccinyl]-betulinic acid (DSB) inhibits HIV-1 replication. DSB functions at a late stage of the virus life cycle but does not inhibit the HIV-1 protease in vitro or interfere with virus assembly or release. DSB specifically delays the cleavage of Gag between the capsid (CA) and p2, resulting in delayed formation of the mature viral core and reduced HIV-1 infectivity. Replication of simian immunodeficiency virus (SIV) was resistant to DSB; however, a chimeric SIV carrying CA-p2 sequences from HIV-1 was inhibited by the drug, indicating that susceptibility to DSB maps to the CA-p2 region of the HIV-1 Gag protein. A single point mutation at the CA-p2 cleavage site of HIV-1 conferred strong resistance to DSB, confirming the target of the drug. HIV-1 strains that are resistant to a variety of protease inhibitors were sensitive to DSB. These findings indicate that DSB specifically protects the CA-p2 cleavage site from processing by the viral protease during virion maturation, thereby revealing a novel mechanism for pharmacologic inhibition of HIV-1 replication.     

Functional tolerance of Streptomyces subtilisin inhibitor toward conformational and stability changes caused by single-point mutations in the hydrophobic core

Single amino acid mutations of Met103 in the hydrophobic core of a serine protease inhibitor, Streptomyces subtilisin inhibitor, caused little change in the inhibitory activity, as measured by the inhibitor constant, although some altered the thermodynamic stability of the protein considerably. (1)H NMR investigations showed that the conformational stress caused by the replacement of Met103 with Gly, Ala, Val, and Ile, namely, the effects of the cavities generated by replacements with smaller side-chains and of the steric distortions generated by beta-branched side-chains, caused considerable changes in the structural arrangement of the side-chains within the core. However, these structural changes were absorbed within the hydrophobic core, without distorting the structure of the reactive site essential for the protein function. These results provide an excellent example of the conformational flexibility of a protein core and the degree of its tolerance of an amino acid replacement. The results also reveal the crucially designed structural relationship between the core of the inhibitor and the enzyme-binding segment with the reactive site in a serine protease inhibitor.     

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.     

C-terminal residues of mature human T-lymphotropic virus type 1 protease are critical for dimerization and catalytic activity

HTLV-1 [HTLV (human T-cell lymphotrophic virus) type 1] is associated with a number of human diseases. HTLV-1 protease is essential for virus replication, and similarly to HIV-1 protease, it is a potential target for chemotherapy. The primary sequence of HTLV-1 protease is substantially longer compared with that of HIV-1 protease, and the role of the ten C-terminal residues is controversial. We have expressed C-terminally-truncated forms of HTLV-1 protease with and without N-terminal His tags. Removal of five of the C-terminal residues caused a 4-40-fold decrease in specificity constants, whereas the removal of an additional five C-terminal residues rendered the protease completely inactive. The addition of the N-terminal His tag dramatically decreased the activity of HTLV-1 protease forms. Pull-down experiments carried out with His-tagged forms, gel-filtration experiments and dimerization assays provided the first unequivocal experimental results for the role of the C-terminal residues in dimerization of the enzyme. There is a hydrophobic tunnel on the surface of HTLV-1 protease close to the C-terminal ends that is absent in the HIV-1 protease. This hydrophobic tunnel can accommodate the extra C-terminal residues of HTLV-1 protease, which was predicted to stabilize the dimer of the full-length enzyme and provides an alternative target site for protease inhibition.     

Accessory mutations balance the marginal stability of the HIV-1 protease in drug resistance

The HIV-1 protease is a major target of inhibitor drugs in AIDS therapies. The therapies are impaired by mutations of the HIV-1 protease that can lead to resistance to protease inhibitors. These mutations are classified into major mutations, which usually occur first and clearly reduce the susceptibility to protease inhibitors, and minor, accessory mutations that occur later and individually do not substantially affect the susceptibility to inhibitors. Major mutations are predominantly located in the active site of the HIV-1 protease and can directly interfere with inhibitor binding. Minor mutations, in contrast, are typically located distal to the active site. A central question is how these distal mutations contribute to resistance development. In this article, we present a systematic computational investigation of stability changes caused by major and minor mutations of the HIV-1 protease. As most small single-domain proteins, the HIV-1 protease is only marginally stable. Mutations that destabilize the folded, active state of the protease therefore can shift the conformational equilibrium towards the unfolded, inactive state. We find that the most frequent major mutations destabilize the HIV-1 protease, whereas roughly half of the frequent minor mutations are stabilizing. An analysis of protease sequences from patients in treatment indicates that the stabilizing minor mutations are frequently correlated with destabilizing major mutations, and that highly resistant HIV-1 proteases exhibit significant fractions of stabilizing mutations. Our results thus indicate a central role of minor mutations in balancing the marginal stability of the protease against the destabilization induced by the most frequent major mutations.