Dynamic flaps in HIV-1 protease adopt unique ordering at different stages in the catalytic cycle

The flexibility of HIV-1 protease flaps is known to be essential for the enzymatic activity. Here we attempt to capture a multitude of conformations of the free and substrate-bound HIV-1 protease that differ drastically in their flap arrangements. The substrate binding process suggests the opening of active site gate in conjunction with a reversal of flap tip ordering, from the native semiopen state. The reversed-flap, open-gated enzyme readily transforms to a closed conformation after proper placement of the substrate into the binding cleft. After substrate processing, the closed state protease which possessed opposite flap ordering relative to the semiopen state, encounters another flap reversal via a second open conformation that facilitates the evolution of native semiopen state of correct flap ordering. The complicated transitional pathway, comprising of many high and low energy states, is explored by combining standard and activated molecular dynamics (MD) simulation techniques. The study not only complements the existing findings from X-ray, NMR, EPR, and MD studies but also provides a wealth of detailed information that could help the structure-based drug design process.     

Adaptability and flexibility of HIV-1 protease.

Even though more than 200 three-dimensional structures of HIV-1 protease complexed to a variety of inhibitors are available in the Protein Data Bank; very few structures of unliganded protein have been determined. We have recently solved structures of unliganded HIV-1 protease tethered dimer mutants to resolutions of 1.9 A and 2.1 A, and have found that the flaps assume closed-flap conformation even in the absence of any bound ligand. We report comparison of the unliganded closed-flap structure with structures of HIV-1 protease inhibitor complexes with a view to accurately identifying structural changes that the ligand can induce on binding to HIV-1 protease in the crystal. These studies reveal that the least flexible region present in the active site of HIV-1 protease need not also be the least adaptable to external stress, thus highlighting the conceptual difference between flexibility and adaptability of proteins in general.     

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.     

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.     

A molecular dynamics study of the structural stability of HIV-1 protease under physiological conditions: the role of Na+ ions in stabilizing the active site

HIV-1 protease is most active under weakly acidic conditions (pH 3.5-6.5), when the catalytic Asp25 and Asp25' residues share 1 proton. At neutral pH, this proton is lost and the stability of the structure is reduced. Here we present an investigation of the effect of pH on the dynamics of HIV-1 protease using MD simulation techniques. MD simulations of the solvated HIV-1 protease with the Asp25/25' residues monoprotonated and deprotonated have been performed. In addition we investigated the effect of the inclusion of Na(+) and Cl(-) ions to mimic physiological salt conditions. The simulations of the monoprotonated form and deprotonated form including Na(+) show very similar behavior. In both cases the protein remained stable in the compact, "self-blocked" conformation in which the active site is blocked by the tips of the flaps. In the deprotonated system a Na(+) ion binds tightly to the catalytic dyad shielding the repulsion between the COO(-) groups. Ab initio calculations also suggest the geometry of the active site with the Na(+) bound closely resembles that of the monoprotonated case. In the simulations of the deprotonated form (without Na(+) ions), a water molecule bound between the Asp25 Asp25' side-chains. This disrupted the dimerization interface and eventually led to a fully open conformation.     

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.     

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.     

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.     

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.     

Rapid structural fluctuations of the free HIV protease flaps in solution: Relationship to crystal structures and comparison with predictions of dynamics calculations

Crystal structures have shown that the HIV-1 protease flaps, domains that control access to the active site, are closed when the active site is occupied by a ligand. Although flap structures ranging from closed to semi-open are observed in the free protease, crystal structures reveal that even the semi-open flaps block access to the active site, indicating that the flaps are mobile in solution. The goals of this paper are to characterize the secondary structure and fast (sub-ns) dynamics of the flaps of the free protease in solution, to relate these results to X-ray structures and to compare them with predictions of dynamics calculations. To this end we have obtained nearly complete backbone and many sidechain signal assignments of a fully active free-protease construct that is stabilized against autoproteolysis by three point mutations. The secondary structure of this protein was characterized using the chemical shift index, measurements of (3h)J(NC') couplings across hydrogen bonds, and NOESY connectivities. Analysis of these measurements indicates that the protease secondary structure becomes irregular near the flap tips, residues 49-53. Model-free analysis of (15)N relaxation parameters, T(1), T(2) (T(1rho)) and (15)N-[(1)H] NOE, shows that residues in the flap tips are flexible on the sub-ns time scale, in contrast with previous observations on the inhibitor-bound protease. These results are compared with theoretical predictions of flap dynamics and the possible biological significance of the sub-ns time scale dynamics of the flap tips is discussed.     

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.     

Comparison of NMR and MD N–H bond order parameters: example of HIV-1 protease

It has been suggested that the cause of disagreements between molecular dynamics (MD) and NMR N–H bond order parameters is the fact that the NMR order parameter is determined for different amino acid residues at different time intervals, while the MD one is derived for all residues from the same MD trajectory of the same time interval. Therefore, it has been proposed for correct comparison with NMR data to calculate the MD order parameter for different amino acid residues separately for trajectory ranges close to NMR correlation time. The MD simulation of the human immunodeficiency virus type-1 protease (HIV-1 PR) with monoprotonated active centre was performed for verification of the proposition. It has been shown that the protease in aqueous solution adopts a set of conformations, which are intermediate between semiopen and closed ones. The calculated MD N–H bond order parameters are in agreement with literature NMR data in confidence interval limits.     

Restrained molecular dynamics simulations of HIV-1 protease: the first step in validating a new target for drug design

To test the anticorrelated relationship that was recently displayed in conventional molecular dynamics (MD) simulations, several different restrained MD simulations on a wild type and on the V82F/I84V drug-resistant mutant of HIV-1 protease were performed. This anticorrelated relationship refers to the observation that compression of the peripheral ear-to-cheek region of HIV protease (i.e., the elbow of the flap to the fulcrum and the cantilever) occurred as the active site flaps were opening, and, conversely, expansion of that ear-to-cheek region occurred as both flaps were closing. An additional examination of this anticorrelated relationship was necessary to determine whether it can be harnessed in a useful manner. Consequently, six different MD experiments were performed that incorporated pairwise distance restraints in that ear-to-cheek region (i.e., the distance between the alpha-carbons of Gly40 and Gln61 was restrained to either 7.7 or 10.5 A, in both monomers). Pushing the backbones of the ear and the cheek regions away from each other slightly did force the flaps that guard the active site to remain closed in both the wild type and the mutant systems-even though there were no ligands in the active sites. Thus, these restrained MD simulations provided evidence that the anticorrelated relationship can be exploited to affect the dynamic behavior of the flaps that guard the active site of HIV-1 protease. These simulations supported our hypothesis of the mechanism governing flap motion, and they are the first step towards validating that peripheral surface as a new target for drug design.     

A molecular dynamics study of the 41-56 beta-hairpin from B1 domain of protein G.

The structural and dynamical behavior of the 41-56 beta-hairpin from the protein G B1 domain (GB1) has been studied at different temperatures using molecular dynamics (MD) simulations in an aqueous environment. The purpose of these simulations is to establish the stability of this hairpin in view of its possible role as a nucleation site for protein folding. The conformation of the peptide in the crystallographic structure of the protein GB1 (native conformation) was lost in all simulations. The new equilibrium conformations are stable for several nanoseconds at 300K (>10 ns), 350 K (>6.5 ns), and even at 450 K (up to 2.5 ns). The new structures have very similar hairpin-like conformations with properties in agreement with available experimental nuclear Overhauser effect (NOE) data. The stability of the structure in the hydrophobic core region during the simulations is consistent with the experimental data and provides further evidence for the role played by hydrophobic interactions in hairpin structures. Essential dynamics analysis shows that the dynamics of the peptide at different temperatures spans basically the same essential subspace. The main equilibrium motions in this subspace involve large fluctuations of the residues in the turn and ends regions. Of the six interchain hydrogen bonds, the inner four remain stable during the simulations. The space spanned by the first two eigenvectors, as sampled at 450 K, includes almost all of the 47 different hairpin structures found in the database. Finally, analysis of the hydration of the 300 K average conformations shows that the hydration sites observed in the native conformation are still well hydrated in the equilibrium MD ensemble.     

HIV-1跨膜蛋白gp41与N-取代吡咯衍生物的结合模式研究

采用分子对接,分子动力学(MD)模拟和分子力学/泊松-波尔兹曼溶剂可有面积方法与分子力学/广义伯恩溶剂可及面积方法(MM-PBSA/MM-GBSA),预测两种N-取代吡咯衍生物与HIV-1 跨膜蛋白gp41疏水口袋的结合模式与作用机理．分子对接采用多种受体构象,并从结果中选取几种可能的结合模式进行MD 模拟,然后通过MM-PBSA计算结合能的方法识别最优的结合模式. MM-PBSA计算结果表明,范德华相互作用是结合的主要驱动力,而极性相互作用决定了配体在结合过程中的取向．进一步的结合能分解显示,配体的羧基与gp41残基Arg579的静电相互作用对结合有重要贡献．上述工作为进一步优化N-取代吡咯衍生物类的HIV-1融合抑制剂建立了良好的理论基础．     

Comparison of inhibitor binding in HIV-1 protease and in non-viral aspartic proteases: the role of the flap

The crystal structure of HIV-1 protease with an inhibitor has been compared with the structures of non-viral aspartic proteases complexed with inhibitors. In the dimeric HIV-1 protease, two 4-stranded beta-sheets are formed by half of the inhibitor, residues 27-29, and the flap from each monomer. In the monomeric non-viral enzyme the single flap does not form a beta-sheet with an inhibitor. The HIV-1 protease shows more interactions with a longer peptide inhibitor than are observed in non-viral aspartic protease-inhibitor complexes. This, and the large movement of the flaps, restricts the conformation of the protease cleavage sites in the retroviral polyprotein precursor.     

Mechanism of drug resistance revealed by the crystal structure of the unliganded HIV-1 protease with F53L mutation

Mutations in HIV-1 protease (PR) that produce resistance to antiviral PR inhibitors are a major problem in AIDS therapy. The mutation F53L arising from antiretroviral therapy was introduced into the flexible flap region of the wild-type PR to study its effect and potential role in developing drug resistance. Compared to wild-type PR, PR(F53L) showed lower (15%) catalytic efficiency, 20-fold weaker inhibition by the clinical drug indinavir, and reduced dimer stability, while the inhibition constants of two peptide analog inhibitors were slightly lower than those for PR. The crystal structure of PR(F53L) was determined in the unliganded form at 1.35 Angstrom resolution in space group P4(1)2(1)2. The tips of the flaps in PR(F53L) had a wider separation than in unliganded wild-type PR, probably due to the absence of hydrophobic interactions of the side-chains of Phe53 and Ile50'. The changes in interactions between the flaps agreed with the reduced stability of PR(F53L) relative to wild-type PR. The altered flap interactions in the unliganded form of PR(F53L) suggest a distinct mechanism for drug resistance, which has not been observed in other common drug-resistant mutants.     

A comparative analysis of the equilibrium dynamics of a designed protein inferred from NMR,X‐ray,and computations

A detailed analysis of high‐resolution structural data and computationally predicted dynamics was carried out for a designed sugar‐binding protein. The mean‐square deviations in the positions of residues derived from nuclear magnetic resonance (NMR) models and those inferred from X‐ray crystallographic B‐factors for two different crystal forms were compared with the predictions based on the Gaussian Network Model (GNM) and the results from molecular dynamics (MD) simulations. GNM systematically yielded a higher correlation than MD, with experimental data, suggesting that the lack of atomistic details in the coarse‐grained GNM is more than compensated for by the mathematically exact evaluation of fluctuations using the native contacts topology. Evidence is provided that particular loop motions are curtailed by intermolecular contacts in the crystal environment causing a discrepancy between theory and experiments. Interestingly, the information conveyed by X‐ray crystallography becomes more consistent with NMR models and computational predictions when ensembles of X‐ray models are considered. Less precise (broadly distributed) ensembles indeed appear to describe the accessible conformational space under native state conditions better than B‐factors. Our results highlight the importance of using multiple conformations obtained by alternative experimental methods, and analyzing results from both coarse‐grained models and atomic simulations, for accurate assessment of motions accessible to proteins under native state conditions. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

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.     

iGNM: a database of protein functional motions based on Gaussian Network Model

MOTIVATION: The knowledge of protein structure is not sufficient for understanding and controlling its function. Function is a dynamic property. Although protein structural information has been rapidly accumulating in databases, little effort has been invested to date toward systematically characterizing protein dynamics. The recent success of analytical methods based on elastic network models, and in particular the Gaussian Network Model (GNM), permits us to perform a high-throughput analysis of the collective dynamics of proteins. RESULTS: We computed the GNM dynamics for 20 058 structures from the Protein Data Bank, and generated information on the equilibrium dynamics at the level of individual residues. The results are stored on a web-based system called iGNM and configured so as to permit the users to visualize or download the results through a standard web browser using a simple search engine. Static and animated images for describing the conformational mobility of proteins over a broad range of normal modes are accessible, along with an online calculation engine available for newly deposited structures. A case study of the dynamics of 20 non-homologous hydrolases is presented to illustrate the utility of the iGNM database for identifying key residues that control the cooperative motions and revealing the connection between collective dynamics and catalytic activity.