The folding free-energy surface of HIV-1 protease: insights into the thermodynamic basis for resistance to inhibitors

Rheostats, toggles, and neutrals, Oh my! A new framework for understanding how amino acid changes modulate protein function

Advances in personalized medicine and protein engineering require accurately predicting outcomes of amino acid substitutions. Many algorithms correctly predict that evolutionarily-conserved positions show "toggle" substitution phenotypes, which is defined when a few substitutions at that position retain function. In contrast, predictions often fail for substitutions at the less-studied "rheostat" positions, which are defined when different amino acid substitutions at a position sample at least half of the possible functional range. This review describes efforts to understand the impact and significance of rheostat positions: (1) They have been observed in globular soluble, integral membrane, and intrinsically disordered proteins; within single proteins, their prevalence can be up to 40%. (2) Substitutions at rheostat positions can have biological consequences and ∼10% of substitutions gain function. (3) Although both rheostat and "neutral" (defined when all substitutions exhibit wild-type function) positions are nonconserved, the two classes have different evolutionary signatures. (4) Some rheostat positions have pleiotropic effects on function, simultaneously modulating multiple parameters (e.g., altering both affinity and allosteric coupling). (5) In structural studies, substitutions at rheostat positions appear to cause only local perturbations; the overall conformations appear unchanged. (6) Measured functional changes show promising correlations with predicted changes in protein dynamics; the emergent properties of predicted, dynamically coupled amino acid networks might explain some of the complex functional outcomes observed when substituting rheostat positions. Overall, rheostat positions provide unique opportunities for using single substitutions to tune protein function. Future studies of these positions will yield important insights into the protein sequence/function relationship.

The denatured state of HIV-1 protease under native conditions

The denatured state of several proteins has been shown to display transient structures that are relevant for folding, stability, and aggregation. To detect them by nuclear magnetic resonance (NMR) spectroscopy, the denatured state must be stabilized by chemical agents or changes in temperature. This makes the environment different from that experienced in biologically relevant processes. Using high‐resolution heteronuclear NMR spectroscopy, we have characterized several denatured states of a monomeric variant of HIV‐1 protease, which is natively structured in water, induced by different concentrations of urea, guanidinium chloride, and acetic acid. We have extrapolated the chemical shifts and the relaxation parameters to the denaturant‐free denatured state at native conditions, showing that they converge to the same values. Subsequently, we characterized the conformational properties of this biologically relevant denatured state under native conditions by advanced molecular dynamics simulations and validated the results by comparison to experimental data. We show that the denatured state of HIV‐1 protease under native conditions displays rich patterns of transient native and non‐native structures, which could be of relevance to its guidance through a complex folding process.

A synergy of activity, stability, and inhibitor-interaction of HIV-1 protease mutants evolved under drug-pressure

A clinically-relevant, drug-resistant mutant of HIV-1 protease (PR), termed Flap+_(I54V) and containing L10I, G48V, I54V and V82A mutations, is known to produce significant changes in the entropy and enthalpy balance of drug-PR interactions, compared to wild-type PR. A similar mutant, Flap+_(I54A) , which evolves from Flap+_(I54V) and contains the single change at residue 54 relative to Flap+_(I54V) , does not. Yet, how Flap+_(I54A) behaves in solution is not known. To understand the molecular basis of V54A evolution, we compared nuclear magnetic resonance (NMR) spectroscopy, fluorescence spectroscopy, isothermal titration calorimetry, and enzymatic assay data from four PR proteins: PR (pWT), Flap+_(I54V) , Flap+_(I54A) , and Flap+_(I54) , a control mutant that contains only L10I, G48V and V82A mutations. Our data consistently show that selection to the smaller side chain at residue 54, not only decreases inhibitor affinity, but also restores the catalytic activity.

Diverse Folding Pathways of HIV-1 Protease Monomer on a Rugged Energy Landscape

The modern energy landscape theory of protein folding predicts multiple folding pathways connecting a myriad of unfolded conformations and a well-defined folded state. However, direct experimental observation of heterogeneous folding pathways is difficult. Naturally evolved proteins typically exhibit a smooth folding energy landscape for fast and efficient folding by avoiding unfavorable kinetic traps. In this case, rapid fluctuations between unfolded conformations result in apparent two-state behavior and make different pathways indistinguishable. However, the landscape roughness can be different, depending on the selection pressures during evolution. Here, we characterize the unusually rugged folding energy landscape of human immunodeficiency virus-1 protease monomer using single-molecule Förster resonance energy transfer spectroscopy. Our data show that fluctuations between unfolded conformations are slow, which enables the experimental observation of heterogeneous folding pathways as predicted by the landscape theory. Although the landscape ruggedness is sensitive to the mutations and fluorophore locations, the folding rate is similar for various protease constructs. The natural evolution of the protease to have a rugged energy landscape likely results from intrinsic pressures to maintain robust folding when human immunodeficiency virus-1 mutates frequently, which is essential for its survival.

Unified understanding of folding and binding mechanisms of globular and intrinsically disordered proteins

Extensive experimental and theoretical studies have advanced our understanding of the mechanisms of folding and binding of globular proteins, and coupled folding and binding of intrinsically disordered proteins (IDPs). The forces responsible for conformational changes and binding are common in both proteins; however, these mechanisms have been separately discussed. Here, we attempt to integrate the mechanisms of coupled folding and binding of IDPs, folding of small and multi-subdomain proteins, folding of multimeric proteins, and ligand binding of globular proteins in terms of conformational selection and induced-fit mechanisms as well as the nucleation–condensation mechanism that is intermediate between them. Accumulating evidence has shown that both the rate of conformational change and apparent rate of binding between interacting elements can determine reaction mechanisms. Coupled folding and binding of IDPs occurs mainly by induced-fit because of the slow folding in the free form, while ligand binding of globular proteins occurs mainly by conformational selection because of rapid conformational change. Protein folding can be regarded as the binding of intramolecular segments accompanied by secondary structure formation. Multi-subdomain proteins fold mainly by the induced-fit (hydrophobic collapse) mechanism, as the connection of interacting segments enhances the binding (compaction) rate. Fewer hydrophobic residues in small proteins reduce the intramolecular binding rate, resulting in the nucleation–condensation mechanism. Thus, the folding and binding of globular proteins and IDPs obey the same general principle, suggesting that the coarse-grained, statistical mechanical model of protein folding is promising for a unified theoretical description of all mechanisms.

Evolution under Drug Pressure Remodels the Folding Free-Energy Landscape of Mature HIV-1 Protease

Using high-pressure NMR spectroscopy and differential scanning calorimetry, we investigate the folding landscape of the mature HIV-1 protease homodimer. The cooperativity of unfolding was measured in the absence or presence of a symmetric active site inhibitor for the optimized wild type protease (PR), its inactive variant PRD25N, and an extremely multidrug-resistant mutant, PR20. The individual fit of the pressure denaturation profiles gives rise to first order, ∆GNMR, and second order, ∆VNMR (the derivative of ∆GNMR with pressure); apparent thermodynamic parameters for each amide proton considered. Heterogeneity in the apparent ∆VNMR values reflects departure from an ideal cooperative unfolding transition. The narrow to broad distribution of ∆VNMR spanning the extremes from inhibitor-free PR20D25N to PR-DMP323 complex, and distinctively for PRD25N-DMP323 complex, indicated large variations in folding cooperativity. Consistent with this data, the shape of thermal unfolding transitions varies from asymmetric for PR to nearly symmetric for PR20, as dimer-inhibitor ternary complexes. Lack of structural cooperativity was observed between regions located close to the active site, including the hinge and tip of the glycine-rich flaps, and the rest of the protein. These results strongly suggest that inhibitor binding drastically decreases the cooperativity of unfolding by trapping the closed flap conformation in a deep energy minimum. To evade this conformational trap, PR20 evolves exhibiting a smoother folding landscape with nearly an ideal two-state (cooperative) unfolding transition. This study highlights the malleability of retroviral protease folding pathways by illustrating how the selection of mutations under drug pressure remodels the free-energy landscape as a primary mechanism.

HIV-1 Protease Dimerization Dynamics Reveals a Transient Druggable Binding Pocket at the Interface

The binding mechanism of HIV-1 protease monomers leading to the catalytically competent dimeric enzyme has been investigated by means of state-of-the-art atomistic simulations. The emerging picture allows a deeper understanding of experimental observations and reveals that water molecules trapped at the interface have an important role in slowing down the kinetics of the association process. Unexpectedly, a cryptic binding pocket is identified at the interface of the complex, corresponding to a partially bound dimer that lacks enzymatic function. The pocket has a transient nature with a lifetime longer than 1 μs, and it displays very favorable druggability features. Docking as well as MM-GBSA free-energy calculations further support the possibility to target the new binding site by means of inhibitors able to prevent the complete dimerization by capturing the inactive conformation. This discovery could open the way to the rational design of a new class of anti-HIV drugs.

The maturation of HIV-1 protease precursor studied by discrete molecular dynamics

The equilibrium properties of a HIV-1-protease precursor are studied by means of an efficient molecular dynamics scheme, which allows for the simulation of the folding of the protein monomers and their dimerization into an active form and compare them with those of the mature protein. The results of the model provide, with atomic detail, an overall account of several experimental findings, including the NMR conformation of the mature dimer, the calorimetric properties of the system, the effects of the precursor tail on the dimerization constant, the secondary chemical shifts of the monomer, and the paramagnetic relaxation enhancement data associated with the conformations of the precursor. It is found that although the mature protein can dimerize in a unique, single way, the precursor populates several dimeric conformations in which monomers are always native-like, but their binding can be non-native.

Gag-Pol processing during HIV-1 virion maturation: a systems biology approach

Proteolytic processing of Gag and Gag-Pol polyproteins by the viral protease (PR) is crucial for the production of infectious HIV-1, and inhibitors of the viral PR are an integral part of current antiretroviral therapy. The process has several layers of complexity (multiple cleavage sites and substrates; multiple enzyme forms; PR auto-processing), which calls for a systems level approach to identify key vulnerabilities and optimal treatment strategies. Here we present the first full reaction kinetics model of proteolytic processing by HIV-1 PR, taking into account all canonical cleavage sites within Gag and Gag-Pol, intermediate products and enzyme forms, enzyme dimerization, the initial auto-cleavage of full-length Gag-Pol as well as self-cleavage of PR. The model allows us to identify the rate limiting step of virion maturation and the parameters with the strongest effect on maturation kinetics. Using the modelling framework, we predict interactions and compensatory potential between individual cleavage rates and drugs, characterize the time course of the process, explain the steep dose response curves associated with PR inhibitors and gain new insights into drug action. While the results of the model are subject to limitations arising from the simplifying assumptions used and from the uncertainties in the parameter estimates, the developed framework provides an extendable open-access platform to incorporate new data and hypotheses in the future.

Thermodynamics of strongly allosteric inhibition: a model study of HIV-1 protease

Protein inhibitors that shift the thermodynamic equilibrium towards a denatured state escape, in general, the straightforward framework of competitive or allosteric inhibitors. The equilibrium properties of peptides which compete with the folding, or more precisely destabilize the native state, of the human immunodeficiency virus (HIV)-1 protease monomer are studied within a structure-based model. The effect of peptides that disrupt the hydrophobic core of the protein can still be summarized in terms of an inhibition constant, which depends on the thermal stability of the protein. The state of the protein denatured by such a peptide is more structured than its intrinsic denatured state, but displays the same degree of compactness. Peptides that target less buried regions of the protein are less efficient and display a more complex thermodynamics that cannot be captured in a simple way.

Protein intrinsic disorder as a flexible armor and a weapon of HIV-1

Many proteins and protein regions are disordered in their native, biologically active states. These proteins/regions are abundant in different organisms and carry out important biological functions that complement the functional repertoire of ordered proteins. Viruses, with their highly compact genomes, small proteomes, and high adaptability for fast change in their biological and physical environment utilize many of the advantages of intrinsic disorder. In fact, viral proteins are generally rich in intrinsic disorder, and intrinsically disordered regions are commonly used by viruses to invade the host organisms, to hijack various host systems, and to help viruses in accommodation to their hostile habitats and to manage their economic usage of genetic material. In this review, we focus on the structural peculiarities of HIV-1 proteins, on the abundance of intrinsic disorder in viral proteins, and on the role of intrinsic disorder in their functions.

The H2A-H2B dimeric kinetic intermediate is stabilized by widespread hydrophobic burial with few fully native interactions

The H2A-H2B histone heterodimer folds via monomeric and dimeric kinetic intermediates. Within ∼5 ms, the H2A and H2B polypeptides associate in a nearly diffusion limited reaction to form a dimeric ensemble, denoted I₂ and I₂*, the latter being a subpopulation characterized by a higher content of nonnative structure (NNS). The I₂ ensemble folds to the native heterodimer, N₂, through an observable, first-order kinetic phase. To determine the regions of structure in the I₂ ensemble, we characterized 26 Ala mutants of buried hydrophobic residues, spanning the three helices of the canonical histone folds of H2A and H2B and the H2B C-terminal helix. All but one targeted residue contributed significantly to the stability of I₂, the transition state and N₂; however, only residues in the hydrophobic core of the dimer interface perturbed the I₂* population. Destabilization of I₂* correlated with slower folding rates, implying that NNS is not a kinetic trap but rather accelerates folding. The pattern of Φ values indicated that residues forming intramolecular interactions in the peripheral helices contributed similar stability to I₂ and N₂, but residues involved in intermolecular interactions in the hydrophobic core are only partially folded in I₂. These findings suggest a dimerize-then-rearrange model. Residues throughout the histone fold contribute to the stability of I₂, but after the rapid dimerization reaction, the hydrophobic core of the dimer interface has few fully native interactions. In the transition state leading to N₂, more native-like interactions are developed and nonnative interactions are rearranged.

The evolution of catalytic function in the HIV-1 protease

The evolution of species is a complex phenomenon based on the optimization of a multidimensional function referred to as fitness. At the level of biomolecular evolution, the fitness function can be reduced to include physiochemical properties relevant to the biological function of a particular molecule. In this work, questions involving the physical-chemical mechanisms underlying the evolution of HIV-1 protease are addressed through molecular simulation and subsequent analysis of thermodynamic properties related to the activity of the enzyme. Specifically, the impact of 40 single amino acid mutations on the binding affinity toward the matrix/capsid (MA/CA) substrate and corresponding transition state intermediate has been characterized using a molecular mechanics Poisson-Boltzmann surface area approach. We demonstrate that this approach is capable of extracting statistically significant information relevant to experimentally determined catalytic activity. Further, no correlation was observed between the effect of mutations on substrate and transition state binding, suggesting independent evolutionary pathways toward optimizing substrate specificity and catalytic activity. In addition, a detailed analysis of calculated binding affinity data suggests that ground-state destabilization (reduced binding affinity for the substrate) could be a contributing factor in the evolutionary optimization of HIV-1 protease. A numerical model is developed to demonstrate that ground-state destabilization is a valid mechanism for activity optimization given the high concentrations of substrate experienced by the functional enzyme in vivo.

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.

Insights into the mechanism of drug resistance: X-ray structure analysis of G48V/C95F tethered HIV-1 protease dimer/saquinavir complex

The mutation G48V in HIV-1 protease is a major resistance mutation against the drug saquinavir. Recently, G48V mutation is found to co-exist with the mutation C95F in AIDS patients treated with saquinavir. We report here the three-dimensional crystal structure of G48V/C95F tethered HIV-1 protease/saquinavir complex. The structure indicates following as the possible causes of drug resistance: (1) loss of direct van der Waals interactions between saquinavir and enzyme residues PHE-53 and PRO-1081, (2) loss of water-mediated hydrogen bonds between the carbonyl oxygen atoms in saquinavir and amide nitrogen atoms of flap residues 50 and 1050, (3) changes in inter-monomer interactions, which could affect the energetics of domain movements associated with inhibitor-binding, and (4) significant reduction in the stability of the mutant dimer. The present structure also provides a rationale for the clinical observation that the resistance mutations C95F/G48V/V82A occur as a cluster in AIDS patients.

Metal-free ALS variants of dimeric human Cu,Zn-superoxide dismutase have enhanced populations of monomeric species

Amino acid replacements at dozens of positions in the dimeric protein human, Cu,Zn superoxide dismutase (SOD1) can cause amyotrophic lateral sclerosis (ALS). Although it has long been hypothesized that these mutations might enhance the populations of marginally-stable aggregation-prone species responsible for cellular toxicity, there has been little quantitative evidence to support this notion. Perturbations of the folding free energy landscapes of metal-free versions of five ALS-inducing variants, A4V, L38V, G93A, L106V and S134N SOD1, were determined with a global analysis of kinetic and thermodynamic folding data for dimeric and stable monomeric versions of these variants. Utilizing this global analysis approach, the perturbations on the global stability in response to mutation can be partitioned between the monomer folding and association steps, and the effects of mutation on the populations of the folded and unfolded monomeric states can be determined. The 2- to 10-fold increase in the population of the folded monomeric state for A4V, L38V and L106V and the 80- to 480-fold increase in the population of the unfolded monomeric states for all but S134N would dramatically increase their propensity for aggregation through high-order nucleation reactions. The wild-type-like populations of these states for the metal-binding region S134N variant suggest that even wild-type SOD1 may also be prone to aggregation in the absence of metals.

Interactions across the interface contribute the stability of homodimeric 3alpha-hydroxysteroid dehydrogenase/carbonyl reductase

The dimerization of 3alpha-hydroxysteroid dehydrogenase/carbonyl reductase was studied by interrupting the salt bridge interactions between D249 and R167 in the dimeric interface. Substitution of alanine, lysine and serine for D249 decreased catalytic efficiency 30, 1400 and 1.4-fold, and lowered the melting temperature 6.9, 5.4 and 7.6 degrees C, respectively. The mutated enzymes have the dimeric species but the equilibrium between monomer and dimer for these mutants varies from each other, implying that these residues might contribute differently to the dimer stability. Thermal and urea-induced unfolding profiles for wild-type and mutant enzymes appeared as a two-state transition and three-state transition, respectively. In addition, mutation on D249 breaks the salt bridges and causes different effects on the loss of enzymatic activity for D249A, D249K and D249S mutants in the urea-induced unfolding profiles. Hence, D249 at the dimeric interface in 3alpha-HSD/CR is essential for conformational stability, oligomeric integrity and enzymatic activity.