Revisiting monomeric HIV-1 protease. Characterization and redesign for improved properties

Evolution of drug resistance drives destabilization of flap region dynamics in HIV-1 protease

The HIV-1 protease is one of several common key targets of combination drug therapies for human immunodeficiency virus infection and acquired immunodeficiency syndrome. During the progression of the disease, some individual patients acquire drug resistance due to mutational hotspots on the viral proteins targeted by combination drug therapies. It has recently been discovered that drug-resistant mutations accumulate on the "flap region" of the HIV-1 protease, which is a critical dynamic region involved in nonspecific polypeptide binding during invasion and infection of the host cell. In this study, we utilize machine learning-assisted comparative molecular dynamics, conducted at single amino acid site resolution, to investigate the dynamic changes that occur during functional dimerization and drug binding of wild-type and common drug-resistant versions of the main protease. We also use a multiagent machine learning model to identify conserved dynamics of the HIV-1 main protease that are preserved across simian and feline protease orthologs. We find that a key conserved functional site in the flap region, a solvent-exposed isoleucine (Ile50) that controls flap dynamics is functionally targeted by drug resistance mutations, leading to amplified molecular dynamics affecting the functional ability of the flap region to hold the drugs. We conclude that better long-term patient outcomes may be achieved by designing drugs that target protease regions that are less dependent upon single sites with large functional binding effects.

Relative domain orientation of the L289K HIV-1 reverse transcriptase monomer

HIV-1 reverse transcriptase (RT) is a heterodimer comprised p66 and p51 subunits (p66/p51). Several single amino acid substitutions in RT, including L289K, decrease p66/p51 dimer affinity, and reduce enzymatic functioning. Here, small-angle X-ray scattering (SAXS) with proton paramagnetic relaxation enhancement (PRE), 19 F site-specific NMR, and size exclusion chromatography (SEC) were performed for the p66 monomer with the L289K mutation, p66L289K . NMR and SAXS experiments clearly elucidated that the thumb and RNH domains in the monomer do not rigidly interact with each other but are spatially close to the RNH domain. Based on this structural model of the monomer, p66L289K and p51 were predicted to form a heterodimer while p66 and p51L289K not. We tested this hypothesis by SEC analysis of p66 and p51 containing L289K in different combinations and clearly demonstrated that L289K substitution in the p51 subunit, but not in the p66 subunit, reduces p66/p51 formation. Based on the derived monomer model and the importance of the inter-subunit RNH-thumb domain interaction in p66/p51, validated by SEC, the mechanism of p66 homodimer formation was discussed.

Understanding the mechanism of HIV-1 protease inhibition by monoclonal antibodies

HIV-1 protease is an essential enzyme in the life cycle of human immunodeficiency virus (HIV) and hence is one of the most important targets for antiviral drug design. Although there are ten FDA approved drugs against HIV protease (PR), their long term usage elicits mutations leading to drug resistance. As a result, novel therapeutic approaches are being explored including synthetic antibodies. Recently, a murine monoclonal antibody, mAB1696 (mAB) was reported to inhibit PR by preventing dimerization. Crystallographic data could reveal only six protease residues that interact with mAB. The present study employs a range of computational techniques, starting from protein-protein docking to all-atomic molecular dynamics simulations to generate plausible 3D structures of PR-mAB complex. Results show that mAB interacts very strongly with several PR dimer interface residues, such as Gln7, Arg8 (N-terminal), Cys95, Leu97 (C-terminal), Thr26, Gly27 (active site), Gly49, Ile50 (flap), apart from its interactions with the PR epitope region, Pro1-Trp6 (N-terminal). These observations support the hypothesis that binding of mAB prevents the dimerization of PR. The interactions and binding conformations identified in this study could form the basis for designing allosteric inhibitors preventing the dimerization of HIV-1 Protease.

A Single Amino Acid Substitution at the HIV-1 Protease Termini Dimer Interface Significantly Reduces Viral Particles Processing Efficiency

The dimeric form of HIV-1 protease (PR) is required for its full proteolytic activity. The stability of the dimer primarily depends on the termini interface, with N-terminal residues 1-4 of one monomer encountering C-terminal residues 96-99 of another. We made an alanine substitution for valine 3 (V3) or leucine 97 (L97) at the termini dimer interface and tested their proteolytic activity. We found that an alanine substitution for L97 (PRL97A) completely inhibited the proteolytic activity of the PR. However, an alanine substitution for V3 (PRV3A) partially impaired the proteolytic activity. We then introduced two forced-dimerization systems involving nucleocapsid (NC) replacement or the addition of 1-2 leucine zippers to determine whether the proteolytic activity of dimer-defective PRs could be restored. We found that two forced-dimerization systems compensated for the defect in PRV3A, but not in PRL97A. This implies that PRV3A and PRL97A potentially impair the PR via different mechanisms or cause defects in PR activity to different extents. These novel findings will likely serve as a foundation for developing new PR inhibitors for treating drug-resistant HIV-1 infections in the future.

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.

Inhibition of human endogenous retrovirus-K by antiretroviral drugs

Background

Human endogenous retroviruses (HERVs) are genomic sequences of retroviral origin which were believed to be integrated into germline chromosomes millions of years ago and account for nearly 8% of the human genome. Although mostly defective and inactive, some of the HERVs may be activated under certain physiological and pathological conditions. While no drugs are designed specifically targeting HERVs, there are a panel of antiretroviral drugs designed against the human immunodeficiency virus and approved by the Federal Drug Administration (FDA).

Results

We determined if these antiretroviral drugs may also be effective in inhibiting HERVs. We constructed a plasmid with consensus HERV-K sequence for testing the effect of antiretroviral drugs on HERV-K. We first determined the effects of nucleoside and non-nucleotide reverse transcriptase (RT) inhibitors on HERV-K by product enhanced reverse transcription assay. We found that all RT inhibitors could significantly inhibit HERV-K RT activity. To determine the effects of antiretroviral drugs on HERV-K infection and viral production, we pseudotyped HERV-K with VSV-G and used the pseudotyped HERV-K virus to infect HeLa cells. HERV-K production was measured by quantitative real time polymerase chain reaction. We found that RT inhibitors Abacavir and Zidovudine, and integrase inhibitor Raltegravir could effectively block HERV-K infection and production. However, protease inhibitors were not as effective as RT and integrase inhibitors.

Conclusions

In summary, we identified several FDA approved antiretroviral drugs that can effectively inhibit HERV-K. These antiretrovirals may open new prospects for studying HERV-K pathophysiology and potentially for exploring treatment of diseases in which HERV-K has been implicated.

Binding of Clinical Inhibitors to a Model Precursor of a Rationally Selected Multidrug Resistant HIV-1 Protease Is Significantly Weaker Than That to the Released Mature Enzyme

We have systematically validated the activity and inhibition of a HIV-1 protease (PR) variant bearing 17 mutations (PR(S17)), selected to represent high resistance by machine learning on genotype-phenotype data. Three of five mutations in PR(S17) correlating with major drug resistance, M46L, G48V, and V82S, and five of 11 natural variations differ from the mutations in two clinically derived extreme mutants, PR20 and PR22 bearing 19 and 22 mutations, respectively. PR(S17), which forms a stable dimer (<10 nM), is ∼10- and 2-fold less efficient in processing the Gag polyprotein than the wild type and PR20, respectively, but maintains the same cleavage order. Isolation of a model precursor of PR(S17) flanked by the 56-amino acid transframe region (TFP-p6pol) at its N-terminus, which is impossible upon expression of an analogous PR20 precursor, allowed systematic comparison of inhibition of TFP-p6pol-PR(S17) and mature PR(S17). Resistance of PR(S17) to eight protease inhibitors (PIs) relative to PR (Ki) increases by 1.5-5 orders of magnitude from 0.01 to 8.4 μM. Amprenavir, darunavir, atazanavir, and lopinavir, the most effective of the eight PIs, inhibit precursor autoprocessing at the p6pol/PR site with IC50 values ranging from ∼7.5 to 60 μM. Thus, this process, crucial for stable dimer formation, shows inhibition ∼200-800-fold weaker than that of the mature PR(S17). TFP/p6pol cleavage, which occurs faster, is inhibited even more weakly by all PIs except darunavir (IC50 = 15 μM); amprenavir shows a 2-fold increase in IC50 (∼15 μM), and atazanavir and lopinavir show increased IC50 values of >42 and >70 μM, respectively.

Dynamical Network of HIV-1 Protease Mutants Reveals the Mechanism of Drug Resistance and Unhindered Activity

HIV-1 protease variants resist drugs by active and non-active-site mutations. The active-site mutations, which are the primary or first set of mutations, hamper the stability of the enzyme and resist the drugs minimally. As a result, secondary mutations that not only increase protein stability for unhindered catalytic activity but also resist drugs very effectively arise. While the mechanism of drug resistance of the active-site mutations is through modulating the active-site pocket volume, the mechanism of drug resistance of the non-active-site mutations is unclear. Moreover, how these allosteric mutations, which are 8-21 Å distant, communicate to the active site for drug efflux is completely unexplored. Results from molecular dynamics simulations suggest that the primary mechanism of drug resistance of the secondary mutations involves opening of the flexible protease flaps. Results from both residue- and community-based network analyses reveal that this precise action of protease is accomplished by the presence of robust communication paths between the mutational sites and the functionally relevant regions: active site and flaps. While the communication is more direct in the wild type, it traverses across multiple intermediate residues in mutants, leading to weak signaling and unregulated motions of flaps. The global integrity of the protease network is, however, maintained through the neighboring residues, which exhibit high degrees of conservation, consistent with clinical data and mutagenesis studies.

Drug resistance mechanisms of three mutations V32I, I47V and V82I in HIV-1 protease toward inhibitors probed by molecular dynamics simulations and binding free energy predictions

Molecular dynamics simulation and binding free energy calculations were used to probe drug resistance of HIV-1 protease mutations toward inhibitors.

Mutations Proximal to Sites of Autoproteolysis and the α-Helix That Co-evolve under Drug Pressure Modulate the Autoprocessing and Vitality of HIV-1 Protease

N-Terminal self-cleavage (autoprocessing) of the HIV-1 protease precursor is crucial for liberating the active dimer. Under drug pressure, evolving mutations are predicted to modulate autoprocessing, and the reduced catalytic activity of the mature protease (PR) is likely compensated by enhanced conformational/dimer stability and reduced susceptibility to self-degradation (autoproteolysis). One such highly evolved, multidrug resistant protease, PR20, bears 19 mutations contiguous to sites of autoproteolysis in retroviral proteases, namely clusters 1-3 comprising residues 30-37, 60-67, and 88-95, respectively, accounting for 11 of the 19 mutations. By systematically replacing corresponding clusters in PR with those of PR20, and vice versa, we assess their influence on the properties mentioned above and observe no strict correlation. A 10-35-fold decrease in the cleavage efficiency of peptide substrates by PR20, relative to PR, is reflected by an only ∼4-fold decrease in the rate of Gag processing with no change in cleavage order. Importantly, optimal N-terminal autoprocessing requires all 19 PR20 mutations as evaluated in vitro using the model precursor TFR-PR20 in which PR is flanked by the transframe region. Substituting PR20 cluster 3 into TFR-PR (TFR-PR(PR20-3)) requires the presence of PR20 cluster 1 and/or 2 for autoprocessing. In accordance, substituting PR clusters 1 and 2 into TFR-PR20 affects the rate of autoprocessing more drastically (>300-fold) compared to that of TFR-PR(PR20-3) because of the cumulative effect of eight noncluster mutations present in TFR-PR20(PR-12). Overall, these studies imply that drug resistance involves a complex synchronized selection of mutations modulating all of the properties mentioned above governing PR regulation and function.

Dimerization of HIV-1 protease occurs through two steps relating to the mechanism of protease dimerization inhibition by darunavir

Dimerization of HIV-1 protease (PR) subunits is an essential process for PR's acquisition of proteolytic activity, which plays a critical role in the maturation of HIV-1. Recombinant wild-type PR (PR(WT)) proved to dimerize, as examined with electrospray ionization mass spectrometry; however, two active site interface PR mutants (PR(T26A) and PR(R87K)) remained monomeric. On the other hand, two termini interface PR mutants (PR(1-C95A) and PR(97/99)) took both monomeric and dimeric forms. Differential scanning fluorimetry indicated that PR(1-C95A) and PR(97/99) dimers were substantially less stable than PR(WT) dimers. These data indicate that intermolecular interactions of two monomers occur first at the active site interface, generating unstable or transient dimers, and interactions at the termini interface subsequently occur, generating stable dimers. Darunavir (DRV), an HIV-1 protease inhibitor, inhibits not only proteolytic activity but also PR dimerization. DRV bound to protease monomers in a one-to-one molar ratio, inhibiting the first step of PR dimerization, whereas conventional protease inhibitors (such as saquinavir) that inhibit enzymatic activity but not dimerization failed to bind to monomers. DRV also bound to mutant PRs containing the transframe region-added PR (TFR-PR(D25N) and TFR-PR(D25N-7AA)), whereas saquinavir did not bind to TFR-PR(D25N) or TFR-PR(D25N-7AA). Notably, DRV failed to bind to mutant PR containing four amino acid substitutions (V32I, L33F, I54M, and I84V) that confer resistance to DRV on HIV-1. To our knowledge, the present report represents the first demonstration of the two-step PR dimerization dynamics and the mechanism of dimerization inhibition by DRV, which should help design further, more potent novel PIs.

Drug Resistance Mutations Alter Dynamics of Inhibitor-Bound HIV-1 Protease

Under the selective pressure of therapy, HIV-1 protease mutants resistant to inhibitors evolve to confer drug resistance. Such mutations can impact both the dynamics and structures of the bound and unbound forms of the enzyme. Flap+ is a multidrug-resistant variant of HIV-1 protease with a combination of primary and secondary resistance mutations (L10I, G48V, I54V, V82A) and a strikingly altered thermodynamic profile for darunavir (DRV) binding relative to the wild-type protease. We elucidated the impact of these mutations on protein dynamics in the DRV-bound state using molecular dynamics simulations and NMR relaxation experiments. Both methods concur in that the conformational ensemble and dynamics of protease are impacted by the drug resistance mutations in Flap+ variant. Surprisingly this change in ensemble dynamics is different from that observed in the unliganded form of the same variant (Cai, Y. et al. J. Chem. Theory Comput.2012, 8, 3452–3462). Our comparative analysis of both inhibitor-free and bound states presents a comprehensive picture of the altered dynamics in drug-resistant mutant HIV-1 protease and underlies the importance of incorporating dynamic analysis of the whole system, including the unliganded state, into revealing drug resistance mechanisms.

Investigation on the mechanism for the binding and drug resistance of wild type and mutations of G86 residue in HIV-1 protease complexed with Darunavir by molecular dynamic simulation and free energy calculation

Residue Gly86 is considered as the highly conversed residue in the HIV-1 protease. In our work, the detailed binding free energies for the wild-type (WT) and mutated proteases binding to the TMC-114 are estimated to investigate the protein-inhibitor binding and drug resistance mechanism by molecule dynamic simulations and molecular mechanics Poisson Boltzmann surface area (MM-PBSA) method. The binding affinities between the mutants and inhibitor are different than that in the wild-type complex and the major resistance to Darunavir (DRV) of G86A and G86S originate from the electrostatic energy and entropy, respectively. Furthermore, free energy decomposition analysis for the WT and mutated complexes on the basis of per-residue indicates that the mutagenesis influences the energy contribution of the residue located at three regions: active site region (residue 24-32), the flap region, and the region around the mutated residue G86 (residue 79-88), especially the flap region. Finally, further hydrogen bonds and structure analysis are carried out to detect the relationship between the energy and conformation. In all, the G86 mutations change the flap region's conformation. The experimental results are in good agreement with available results.

Enhanced stability of monomer fold correlates with extreme drug resistance of HIV-1 protease

During treatment, mutations in HIV-1 protease (PR) are selected rapidly that confer resistance by decreasing affinity to clinical protease inhibitors (PIs). As these unique drug resistance mutations can compromise the fitness of the virus to replicate, mutations that restore conformational stability and activity while retaining drug resistance are selected on further evolution. Here we identify several compensating mechanisms by which an extreme drug-resistant mutant bearing 20 mutations (PR20) with >5-fold increased Kd and >4000-fold decreased affinity to the PI darunavir functions. (1) PR20 cleaves, albeit poorly, Gag polyprotein substrates essential for viral maturation. (2) PR20 dimer, which exhibits distinctly enhanced thermal stability, has highly attenuated autoproteolysis, thus likely prolonging its lifetime in vivo. (3) The enhanced stability of PR20 results from stabilization of the monomer fold. Both monomeric PR20(T26A) and dimeric PR20 exhibit Tm values 6-7.5 °C higher than those for their PR counterparts. Two specific mutations in PR20, L33F and L63P at sites of autoproteolysis, increase the Tm of monomeric PR(T26A) by ~8 °C, similar to PR20(T26A). However, without other compensatory mutations as seen in PR20, L33F and L63P substitutions, together, neither restrict autoproteolysis nor significantly reduce binding affinity to darunavir. To determine whether dimer stability contributes to binding affinity for inhibitors, we examined single-chain dimers of PR and PR(D25N) in which the corresponding identical monomer units were covalently linked by GGSSG sequence. Linking of the subunits did not appreciably change the ΔTm on inhibitor binding; thus stabilization by tethering appears to have little direct effect on enhancing inhibitor affinity.

Differential Flap Dynamics in Wild-type and a Drug Resistant Variant of HIV-1 Protease Revealed by Molecular Dynamics and NMR Relaxation

In the rapidly evolving disease of HIV drug resistance readily emerges, nullifying the effectiveness of therapy. Drug resistance has been extensively studied in HIV-1 protease where resistance occurs when the balance between enzyme inhibition and substrate recognition and turn-over is perturbed to favor catalytic activity. Mutations which confer drug resistance can impact the dynamics and structure of both the bound and unbound forms of the enzyme. Flap+ is a multi-drug-resistant variant of HIV-1 protease with a combination of mutations at the edge of the active site, within the active site, and in the flaps (L10I, G48V, I54V, V82A). The impact of these mutations on the dynamics in the unliganded form in comparison with the wild-type protease was elucidated with Molecular Dynamic simulations and NMR relaxation experiments. The comparative analyses from both methods concur in showing that the enzyme's dynamics are impacted by the drug resistance mutations in Flap+ protease. These alterations in the enzyme dynamics, particularly within the flaps, likely modulate the balance between substrate turn-over and drug binding, thereby conferring drug resistance.

How Does Darunavir Prevent HIV-1 Protease Dimerization?

The drug Darunavir (DRV) is a potent inhibitor of HIV-1 protease (PR), a homodimeric essential enzyme of the AIDS virus. Recent experimental data suggest that DRV is able to prevent dimerization of HIV-1 PR, which, together with its high affinity for the mature enzyme, has been linked to the high genetic barrier to the development of viral resistance. The mechanism of dimerization inhibition and the binding mode(s) of DRV to monomeric HIV-1 PR are unknown. Here, multiple molecular dynamics simulations with explicit solvent (for a total of 11 μs with the CHARMM force field and 1 μs with the Amber force field) show that the monomer of HIV-1 PR is structurally stable and reveal a major binding mode of DRV. This binding mode is stabilized by favorable interactions between the apolar groups of DRV and the hydrophobic residues Ile32, Ile47, Ile50, Ile54, Pro79, Val82, and Ile84. The binding mode to monomeric HIV-1 PR identified by molecular dynamics is different from the two binding modes observed in the crystal structure of the complex with dimeric HIV-1 PR. As an example, there are no interactions between DRV and the catalytic Asp25 in the binding mode to monomeric HIV-1 PR revelead by the simulations. In contrast, the simulations show extensive and stable interactions between DRV and the flap (residues 46-55), which are likely to sterically hinder the formation of the flap interface as observed in the dimeric structure. Which of the two mechanisms of inhibition (dimerization inhibition by association with the flap or binding to the active site of the mature enzyme) dominates might depend on the HIV-1 PR mutations, and it is likely that dimerization inhibition is predominant for multiple mutations at the active site in multidrug resistant strains.

Detailed atomistic analysis of the HIV-1 protease interface

HIV-1 protease is a very attractive target for the development of new anti-HIV drugs and has been extensively studied over the past decades. In this study, we present a detailed atomic level characterization of the dimer interface in the enzyme HIV-1 protease through computational alanine scanning mutagenesis and molecular dynamics simulations. In addition to a full mapping of the amino acid residues present at the subunit interface, in terms of the corresponding energetic contribution for dimer formation and of their classification as hot spots, warm spots, and null spots, we trace a dynamic analysis of the subunit interacting and solvent accessible surface areas and of the most important hydrogen bonds between subunits. The results presented illustrate the high energetic importance for dimer formation of a small set of five amino acid residue pairs at the subunit interface-Leu5, Ile50, Arg87, Leu97, and Phe99-and provide important clues on the most important structural and energetic determinants for dimer formation. In addition, the results presented suggest several key targets at the subunit interface for the development of new molecules that aim to inhibit HIV-1 protease (PR) activity through blocking the formation of the fully active PR homodimeric form, providing important clues for drug design.

Globular structure of a human immunodeficiency virus-1 protease (1DIFA dimer) in an effective solvent medium by a Monte Carlo simulation

A coarse-grained model is used to study the structure and dynamics of a human immunodeficiency virus-1 protease (1DIFA dimer) consisting of 198 residues in an effective solvent medium on a cubic lattice by Monte Carlo simulations for a range of interaction strengths. Energy and mobility profiles of residues are found to depend on the interaction strength and exhibit remarkable segmental symmetries in two monomers. Lowest energy residues such as Arg(41) and Arg(140) (most electrostatic and polar) are not the least mobile; despite the higher energy, the hydrophobic residues (Ile, Leu, and Val) are least mobile and form the core by pinning down the local segments for the globular structure. Variations in the gyration radius (R(g)) and energy (E(c)) of the protein show nonmonotonic dependence on the interaction strength with the smallest R(g) around the largest value of E(c). Pinning of the conformations by the hydrophobic residues at high interaction strength seems to provide seed for the protein chain to collapse.

Structural and energetic analysis on the complexes of clinically isolated subtype C HIV-1 proteases and approved inhibitors by molecular dynamics simulation

HIV-1 has a large genetic diversity. Subtype B HIV-1 is commonly found in patients in developed countries. In contrast, an increasing number of patients are infected with the non-B subtype viruses, especially with subtype C HIV-1, in developing countries. It remains to be clarified how mutations or polymorphisms in non-B subtype HIV-1 influence the efficacy of the approved inhibitors. In this study, we have performed molecular dynamics simulations on clinically isolated subtype C HIV-1 proteases in complex with three kinds of approved inhibitors. From the structural and energetic viewpoints, we identified the polymorphisms influencing on the binding of the inhibitors. The effect of the V82I mutation on the association with chemicals and the reason for rare appearance of the D30N mutation in subtype C HIV-1 were discussed in terms of the change of geometry of the residues in HIV-1 protease.

Highly conserved glycine 86 and arginine 87 residues contribute differently to the structure and activity of the mature HIV-1 protease

The structural and functional role of conserved residue G86 in HIV-1 protease (PR) was investigated by NMR and crystallographic analyses of substitution mutations of glycine to alanine and serine (PR(G86A) and PR(G86S)). While PR(G86S) had undetectable catalytic activity, PR(G86A) exhibited approximately 6000-fold lower catalytic activity than PR. (1)H-(15)N NMR correlation spectra revealed that PR(G86A) and PR(G86S) are dimeric, exhibiting dimer dissociation constants (K(d)) of approximately 0.5 and approximately 3.2 muM, respectively, which are significantly lower than that seen for PR with R87K mutation (K(d) > 1 mM). Thus, the G86 mutants, despite being partially dimeric under the assay conditions, are defective in catalyzing substrate hydrolysis. NMR spectra revealed no changes in the chemical shifts even in the presence of excess substrate, indicating very poor binding of the substrate. Both NMR chemical shift data and crystal structures of PR(G86A) and PR(G86S) in the presence of active-site inhibitors indicated high structural similarity to previously described PR/inhibitor complexes, except for specific perturbations within the active site loop and around the mutation site. The crystal structures in the presence of the inhibitor showed that the region around residue 86 was connected to the active site by a conserved network of hydrogen bonds, and the two regions moved further apart in the mutants. Overall, in contrast to the role of R87 in contributing significantly to the dimer stability of PR, G86 is likely to play an important role in maintaining the correct geometry of the active site loop in the PR dimer for substrate binding and hydrolysis. Proteins 2010. (c) 2009 Wiley-Liss, Inc.

Revealing the dimer dissociation and existence of a folded monomer of the mature HIV-2 protease

Purification and in vitro protein-folding schemes were developed to produce monodisperse samples of the mature wild-type HIV-2 protease (PR2), enabling a comprehensive set of biochemical and biophysical studies to assess the dissociation of the dimeric protease. An E37K substitution in PR2 significantly retards autoproteolytic cleavage during expression. Furthermore, it permits convenient measurement of the dimer dissociation of PR2(E37K) (elevated K(d) approximately 20 nM) by enzyme kinetics. Differential scanning calorimetry reveals a T(m) of 60.5 for PR2 as compared with 65.7 degrees C for HIV-1 protease (PR1). Consistent with weaker binding of the clinical inhibitor darunavir (DRV) to PR2, the T(m) of PR2 increases by 14.8 degrees C in the presence of DRV as compared with 22.4 degrees C for PR1. Dimer interface mutations, such as a T26A substitution in the active site (PR2(T26A)) or a deletion of the C-terminal residues 96-99 (PR2(1-95)), drastically increase the K(d) (>10(5)-fold). PR2(T26A) and PR2(1-95) consist predominantly of folded monomers, as determined by nuclear magnetic resonance (NMR) and size-exclusion chromatography coupled with multiangle light scattering and refractive index measurements (SMR), whereas wild-type PR2 and its active-site mutant PR2(D25N) are folded dimers. Addition of twofold excess active-site inhibitor promotes dimerization of PR2(T26A) but not of PR2(1-95), indicating that subunit interactions involving the C-terminal residues are crucial for dimer formation. Use of SMR and NMR with PR2 facilitates probing for potential inhibitors that restrict protein folding and/or dimerization and, thus, may provide insights for the future design of inhibitors to circumvent drug resistance.

Interactions of different inhibitors with active-site aspartyl residues of HIV-1 protease and possible relevance to pepsin

The importance of the active site region aspartyl residues 25 and 29 of the mature HIV-1 protease (PR) for the binding of five clinical and three experimental protease inhibitors [symmetric cyclic urea inhibitor DMP323, nonhydrolyzable substrate analog (RPB) and the generic aspartic protease inhibitor acetyl-pepstatin (Ac-PEP)] was assessed by differential scanning calorimetry. DeltaT(m) values, defined as the difference in T(m) for a given protein in the presence and absence of inhibitor, for PR with DRV, ATV, SQV, RTV, APV, DMP323, RPB, and Ac-PEP are 22.4, 20.8, 19.3, 15.6, 14.3, 14.7, 8.7, and 6.5 degrees C, respectively. Binding of APV and Ac-PEP is most sensitive to the D25N mutation, as shown by DeltaT(m) ratios [DeltaT(m)(PR)/DeltaT(m)(PR(D25N))] of 35.8 and 16.3, respectively, whereas binding of DMP323 and RPB (DeltaT(m) ratios of 1-2) is least affected. Binding of the substrate-like inhibitors RPB and Ac-PEP is nearly abolished (DeltaT(m)(PR)/DeltaT(m)(PR(D29N)) > or = 44) by the D29N mutation, whereas this mutation only moderately affects binding of the smaller inhibitors (DeltaT(m) ratios of 1.4-2.2). Of the nine FDA-approved clinical HIV-1 protease inhibitors screened, APV, RTV, and DRV competitively inhibit porcine pepsin with K(i) values of 0.3, 0.6, and 2.14 microM, respectively. DSC results were consistent with this relatively weak binding of APV (DeltaT(m) 2.7 degrees C) compared with the tight binding of Ac-PEP (DeltaT(m) > or = 17 degrees C). Comparison of superimposed structures of the PR/APV complex with those of PR/Ac-PEP and pepsin/pepstatin A complexes suggests a role for Asp215, Asp32, and Ser219 in pepsin, equivalent to Asp25, Asp25', and Asp29 in PR in the binding and stabilization of the pepsin/APV complex.

Residue energy and mobility in sequence to global structure and dynamics of a HIV-1 protease (1DIFA) by a coarse-grained Monte Carlo simulation

Energy, mobility, and structural profiles of residues in a specific sequence of human immunodeficiency virus (HIV)-1 protease chain and its global conformation and dynamics are studied by a coarse-grained computer simulation model on a cubic lattice. HIV-1 protease is described by a chain of 99 residues (nodes) in a specific sequence (1DIFA) with N- and C-terminals on the lattice, where empty lattice sites represent an effective solvent medium. Internal structures of the residues are ignored but their specificities are captured via an interaction (epsilon(ij)) matrix (residue-residue, residue-solvent) of the coefficient (fepsilon(ij)) of the Lennard-Jones potential. Simulations are performed for a range of interaction strength (f) with the solvent-residue interaction describing the quality of the solvent. Snapshots of the protein show considerable changes in the conformation of the protein on varying the interaction. From the mobility and energy profiles of the residues, it is possible to identify the active (and not so active) segments of the protein and consequently their role in proteolysis. Contrary to interaction thermodynamics, the hydrophobic residues possess higher configurational energy and lower mobility while the electrostatic and polar residues are more mobile despite their lower interaction energy. Segments of hydrophobic core residues, crucial for the structural evolution of the protein are identified-some of which are consistent with recent molecular dynamics simulation in context to possible clinical observations. Global energy and radius of gyration of the protein exhibit nonmonotonic dependence on the interaction strength (f) with opposite trends, e.g., rapid transition into globular structure with higher energy. Variations of the rms displacement of the protein and that of a tracer residue, Gly(49), with the time steps show how they slow down on increasing the interaction strength.

Analysis and characterization of dimerization inhibition of a multi-drug-resistant human immunodeficiency virus type 1 protease using a novel size-exclusion chromatographic approach

Active-site inhibitors of HIV-1 PR (protease) block viral replication by preventing viral maturation. However, HIV-1 often develops resistance to active-site inhibitors through multiple mutations in PR and therefore recent efforts have focused on inhibiting PR dimerization as an alternative approach. Dimerization inhibitors have been identified using kinetic analysis, but additional characterization of the effect of these inhibitors on PR by physical methods has been difficult. In the present study, we identified a PR(MDR) (multi-drug-resistant HIV-1 PR) that was highly resistant to autoproteolysis. Using this PR and a novel size-exclusion chromatographic approach that incorporated fluorescence and MS detection, we were able to demonstrate inhibition of dimerization using P27 (peptide 27), a peptide dimerization inhibitor of PR previously identified on the basis of kinetic analysis. Incubation of PR(MDR) with P27, or other dimerization inhibitors, led to a dose- and time-dependent formation of PR monomers based on the change in elution time by size exclusion and its similar elution time to engineered forms of monomeric PR, namely PR(T26A) and glutathionylated PR. In contrast, incubation of PR(MDR) with a potent active-site inhibitor did not change the elution time for the PR(MDR) dimer. The monomeric PR induced by P27 had fluorescent characteristics which were consistent with unfolded PR. Structure-activity studies identified the active regions of P27 and experiments were performed to examine the effect of other dimerization inhibitors on PR. The present study is the first characterization of dimerization inhibition of PR(MDR), a prime target for these inhibitors, using a novel size-exclusion chromatographic approach.

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

Spontaneous mutations at numerous sites distant from the active site of human immunodeficiency virus type 1 protease enable resistance to inhibitors while retaining enzymatic activity. As a benchmark for probing the effects of these mutations on the conformational adaptability of this dimeric beta-barrel protein, the folding free-energy surface of a pseudo-wild-type variant, HIV-PR(*), was determined by a combination of equilibrium and kinetic experiments on the urea-induced unfolding/refolding reactions. The equilibrium unfolding reaction was well described by a two-state model involving only the native dimeric form and the unfolded monomer. The global analysis of the kinetic folding mechanism reveals the presence of a fully folded monomeric intermediate that associates to form the native dimeric structure. Independent analysis of a stable monomeric version of the protease demonstrated that a small-amplitude fluorescence phase in refolding and unfolding, not included in the global analysis of the dimeric protein, reflects the presence of a transient intermediate in the monomer folding reaction. The partially folded and fully folded monomers are only marginally stable with respect to the unfolded state, and the dimerization reaction provides a modest driving force at micromolar concentrations of protein. The thermodynamic properties of this system are such that mutations can readily shift the equilibrium from the dimeric native state towards weakly folded states that have a lower affinity for inhibitors but that could be induced to bind to their target proteolytic sites. Presumably, subsequent secondary mutations increase the stability of the native dimeric state in these variants and, thereby, optimize the catalytic properties of the resistant human immunodeficiency virus type 1 protease.

Analysis of correlations between energy and residue fluctuations in native proteins and determination of specific sites for binding

The Gaussian network model is used to derive the correlations between energy and residue fluctuations in native proteins. Residues are identified that respond strongly to energy fluctuations and that display correlations with the remaining residues of the protein at the highest modes. We postulate that these residues are located at specific sites for drug binding. We test the validity of this postulate on a data set of 33 structurally distinct proteins in the unbound state. Detailed results are presented for drug binding to the HIV protease.

Novel macromolecular inhibitors of human immunodeficiency virus-1 protease

An intracellularly expressed defective human immunodeficiency virus type-1 (HIV-1) protease (PR) monomer could function as a dominant-negative inhibitor of the enzyme that requires dimerization for activity. Based on in silico studies, two mutant PRs harboring hydrophilic mutations, capable of forming favorable inter- and intra-subunit interactions, were selected: PR(RE) containing Asp25Arg and Gly49Glu mutations, and PR(RER) containing an additional Ile50Arg mutation. The mutants were expressed and tested by PR assays, nuclear magnetic resonance (NMR) and cell culture experiments. The mutant PRs showed dose-dependent inhibition of the wild-type PR in a fluorescent microtiter plate PR assay. Furthermore, both mutants were retained by hexahistidine-tagged wild-type HIV-1 PR immobilized on nickel-chelate affinity resin. For the first time, heterodimerization between wild-type and dominant-negative mutant PRs were also demonstrated by NMR spectroscopy. (1)H-(15)N Heteronuclear Single Quantum Coherence NMR spectra showed that although PR(RE) has a high tendency to aggregate, PR(RER) exists mainly as a folded monomer at 25-35 microM concentration, but in the presence of wild-type PR in a ratio of 1:1, heterodimerization occurs with both mutants. While the recombinant virus containing the PR(RE) sequence showed only very low level of expression, expression of the viral proteins of the virus with the PR(RER) sequence was comparable with that of the wild-type. In cell culture experiments, infectivity of viral particles containing PR(RER) protein was reduced by 82%, at mutant to wild-type infective DNA ratio of 2:1.

The solution structure of the simian foamy virus protease reveals a monomeric protein

In contrast to orthoretroviruses, foamy viruses (FVs) express their Pol polyprotein from a separate pol-specific transcript. Only the integrase domain is cleaved off, leading to a protease-reverse transcriptase (PR-RT) protein. We purified the separate PR domain (PRshort) of simian FV from macaques by expressing the recombinant gene in Escherichia coli. Sedimentation analyses and size exclusion chromatography indicate that PRshort is a stable monomer in solution. This allowed us to determine the structure of the PRshort monomer using 1426 experimental restraints derived from NMR spectroscopy. The superposition of 20 conformers resulted in a backbone atom rmsd of 0.55 A for residues Gln8-Leu93. Although the overall folds are similar, the macaque simian FV PRshort reveals significant differences in the dimerization interface relative to other retroviral PRs, such as HIV-1 (human immunodeficiency virus type 1) PR, which appear to be rather stable dimers. Especially the flap region and the N- and C-termini of PRshort are highly flexible. Neglecting these regions, the backbone atom rmsd drops to 0.32 A, highlighting the good definition of the central part of the protein. To exclude that the monomeric state of PRshort is due to cleaving off the RT, we purified the complete PR-RT and performed size exclusion chromatography. Our data show that PR-RT is also monomeric. We thus conclude adoption of a monomeric state of PR-RT to be a regulatory mechanism to inhibit PR activity before virus assembly in order to reduce packaging problems. Dimerization might therefore be triggered by additional viral or cellular factors.

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 beta-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.

Potent Inhibition of HIV-1 Replication by Novel Non-peptidyl Small Molecule Inhibitors of Protease Dimerization

Dimerization of HIV-1 protease subunits is essential for its proteolytic activity, which plays a critical role in HIV-1 replication. Hence, the inhibition of protease dimerization represents a unique target for potential intervention of HIV-1. We developed an intermolecular fluorescence resonance energy transfer-based HIV-1-expression assay employing cyan and yellow fluorescent protein-tagged protease monomers. Using this assay, we identified non-peptidyl small molecule inhibitors of protease dimerization. These inhibitors, including darunavir and two experimental protease inhibitors, blocked protease dimerization at concentrations of as low as 0.01 microm and blocked HIV-1 replication with IC(50) values of 0.0002-0.48 microm. These agents also inhibited the proteolytic activity of mature protease. Other approved anti-HIV-1 agents examined except tipranavir, a CCR5 inhibitor, and soluble CD4 failed to block the dimerization event. Once protease monomers dimerize to become mature protease, mature protease is not dissociated by this dimerization inhibition mechanism, suggesting that these agents block dimerization at the nascent stage of protease maturation. The proteolytic activity of mature protease that managed to undergo dimerization despite the presence of these agents is likely to be inhibited by the same agents acting as conventional protease inhibitors. Such a dual inhibition mechanism should lead to highly potent inhibition of HIV-1.

Mutational and structural studies aimed at characterizing the monomer of HIV-1 protease and its precursor

An experimental protocol for folding the mature human immunodeficiency virus-1 (HIV-1) protease is presented that facilitates NMR studies at a low protein concentration of approximately 20 micoM. Under these conditions, NMR spectra show that the mature protease lacking its terminal beta-sheet residues 1-4 and 96-99 (PR(5-95)) exhibits a stable monomer fold spanning the region 10-90 that is similar to that of the single subunit of the wild-type dimer and the dimer bearing a D25N mutation (PR(D25N)). Urea-induced unfolding monitored both by changes in (1)H-(15)N heteronuclear single quantum correlation spectra and by protein fluorescence indicates that although PR(5-95) monomer displays a transition profile similar to that of the PR(D25N) dimer (50% unfolded (U(50)) = approximately 1.9 M), extending the protease with 4 residues (SFNF) of its N-terminally flanking sequence in the Gag-Pol precursor ((SFNF)PR(D25N)) decreases the stability of the fold (U(50) = approximately 1.5 M). Assigned backbone chemical shifts were used to elucidate differences in the stability of the PR(T26A) (U(50) = 2.5 M) and (SFNF)PR(D25N) monomers and compared with PR(D25N/T26A) monomer. Discernible differences in the backbone chemical shifts were observed for N-terminal protease residues 3-6 of (SFNF)PR(D25N) that may relate to the increase in the equilibrium dissociation constant (K(d)) and the very low catalytic activity of the protease prior to its autoprocessing at its N terminus from the Gag-Pol precursor.

Inherent chaperone-like activity of aspartic proteases reveals a distant evolutionary relation to double-psi barrel domains of AAA-ATPases

Chaperones and proteases share the ability to interact with unfolded proteins. Here we show that enzymatically inactive forms of the aspartic proteases HIV-1 protease and pepsin have inherent chaperone-like activity and can prevent the aggregation of denatured substrate proteins. In contrast to proteolysis, which requires dimeric enzymes, chaperone-like activity could be observed also with monomeric domains. The involvement of the active site cleft in the chaperone-like function was demonstrated by the inhibitory effect of peptide substrate inhibitors. The high structural similarity between aspartic proteases and the N-terminal double-psi barrels of Cdc48-like proteins, which are involved in the unfolding and dissociation of proteins, suggests that they share a common ancestor. The latent chaperone-like activity in aspartic proteases can be seen as a relic that has further evolved to serve substrate binding in the context of proteolytic activity.

The Role of the S-S Bridge in Retroviral Protease Function and Virion Maturation

Retroviral proteases are translated as a part of Gag-related polyproteins, and are released and activated during particle release. Mason-Pfizer monkey virus (M-PMV) Gag polyproteins assemble into immature capsids within the cytoplasm of the host cells; however, their processing occurs only after transport to the plasma membrane and subsequent release. Thus, the activity of M-PMV protease is expected to be highly regulated during the replication cycle. It has been proposed that reversible oxidation of protease cysteine residues might be responsible for such regulation. We show that cysteine residues in M-PMV protease can form an intramolecular S-S bridge. The disulfide bridge shifts the monomer/dimer equilibrium in favor of the dimer, and increases the proteolytic activity significantly. To investigate the role of this disulfide bridge in virus maturation and replication, we engineered an M-PMV clone in which both protease cysteine residues were replaced by alanine (M-PMV(PRC7A/C106A)). Surprisingly, the cysteine residues were dispensable for Gag polyprotein processing within the virus, indicating that even low levels of protease activity are sufficient for polyprotein processing during maturation. However, the long-term infectivity of M-PMV(PRC7A/C106A) was noticeably compromised. These results show clearly that the proposed redox mechanism does not rely solely on the formation of the stabilizing S-S bridge in the protease. Thus, in addition to the protease disulfide bridge, reversible oxidation of cysteine and/or methionine residues in other domains of the Gag polyprotein or in related cellular proteins must be involved in the regulation of maturation.

Disruption of the HIV-1 protease dimer with interface peptides: Structural studies using NMR spectroscopy combined with [2-13C]-Trp selective labeling

HIV-1 protease (HIV-1 PR), which is encoded by retroviruses, is required for the processing of gag and pol polyprotein precursors, hence it is essential for the production of infectious viral particles. In vitro inhibition of the enzyme results in the production of progeny virions that are immature and noninfectious, suggesting its potential as a therapeutic target for AIDS. Although a number of potent protease inhibitor drugs are now available, the onset of resistance to these agents due to mutations in HIV-1 PR has created an urgent need for new means of HIV-1 PR inhibition. Whereas enzymes are usually inactivated by blocking of the active site, the structure of dimeric HIV-1 PR allows an alternative inhibitory mechanism. Since the active site is formed by two half-enzymes, which are connected by a four-stranded antiparallel beta-sheet involving the N- and C- termini of both monomers, enzyme activity can be abolished by reagents targeting the dimer interface in a region relatively free of mutations would interfere with formation or stability of the functional HIV-1 PR dimer. This strategy has been explored by several groups who targeted the four-stranded antiparallel beta-sheet that contributes close to 75% of the dimerization energy. Interface peptides corresponding to native monomer N- or C-termini of several of their mimetics demonstrated, mainly on the basis of kinetic analyses, to act as dimerization inhibitors. However, to the best of our knowledge, neither X-ray crystallography nor NMR structural studies of the enzyme-inhibitor complex have been performed to date. In this article we report a structural study of the dimerization inhibition of HIV-1 PR by NMR using selective Trp side chain labeling.

Computational mutagenesis studies of protein structure‐function correlations

Topological scores, measures of sequence-structure compatibility, are calculated for all 1,881 single point mutants of the human immunodeficiency virus (HIV)-1 protease using a four-body statistical potential function based on Delaunay tessellation of protein structure. Comparison of the mutant topological score data with experimental data from alanine scan studies specifically on the dimer interface residues supports previous findings that 1) L97 and F99 contribute greatly to the Gibbs energy of HIV-1 protease dimerization, 2) Q2 and T4 contribute the least toward the Gibbs energy, and 3) C-terminal residues are more sensitive to mutations than those at the N-terminus. For a more comprehensive treatment of the relationship between protease structure and function, mutant topological scores are compared with the activity levels for a set of 536 experimentally synthesized protease mutants, and a significant correlation is observed. Finally, this structure-function correlation is similarly identified by examining model systems consisting of 2,015 single point mutants of bacteriophage T4 lysozyme as well as 366 single point mutants of HIV-1 reverse transcriptase and is hypothesized to be a property generally applicable to all proteins.

A folding inhibitor of the HIV-1 protease

Because the human immunodeficiency virus type 1 protease (HIV-1-PR) is an essential enzyme in the viral life cycle, its inhibition can control AIDS. The folding of single-domain proteins, like each of the monomers forming the HIV-1-PR homodimer, is controlled by local elementary structures (LES, folding units stabilized by strongly interacting, highly conserved, as a rule hydrophobic, amino acids). These LES have evolved over myriad generations to recognize and strongly attract each other, so as to make the protein fold fast and be stable in its native conformation. Consequently, peptides displaying a sequence identical to those segments of the monomers associated with LES are expected to act as competitive inhibitors and thus destabilize the native structure of the enzyme. These inhibitors are unlikely to lead to escape mutants as they bind to the protease monomers through highly conserved amino acids, which play an essential role in the folding process. The properties of one of the most promising inhibitors of the folding of the HIV-1-PR monomers found among these peptides are demonstrated with the help of spectrophotometric assays and circular dichroism spectroscopy.

Folding Regulates Autoprocessing of HIV-1 Protease Precursor

Autoprocessing of HIV-1 protease (PR) precursors is a crucial step in the generation of the mature protease. Very little is known regarding the molecular mechanism and regulation of this important process in the viral life cycle. In this context we report here the first and complete residue level investigations on the structural and folding characteristics of the 17-kDa precursor TFR-PR-C(nn) (161 residues) of HIV-1 protease. The precursor shows autoprocessing activity indicating that the solution has a certain population of the folded active dimer. Removal of the 5-residue extension, C(nn) at the C-terminal of PR enhanced the activity to some extent. However, NMR structural characterization of the precursor containing a mutation, D25N in the PR at pH 5.2 and 32 degrees C under different conditions of partial and complete denaturation by urea, indicate that the precursor has a high tendency to be unfolded. The major population in the ensemble displays some weak folding propensities in both the TFR and the PR regions, and many of these in the PR region are the non-native type. As both D25N mutant and wild-type PR are known to fold efficiently to the same native dimeric form, we infer that TFR cleavage enables removal of the non-native type of preferences in the PR domain to cause constructive folding of the protein. These results indicate that intrinsic structural and folding preferences in the precursor would have important regulatory roles in the autoprocessing reaction and generation of the mature enzyme.

The folding and dimerization of HIV-1 protease: evidence for a stable monomer from simulations

HIV-1 protease (PR) is a major drug target in combating AIDS, as it plays a key role in maturation and replication of the virus. Six FDA-approved drugs are currently in clinical use, all designed to inhibit enzyme activity by blocking the active site, which exists only in the dimer. An alternative inhibition mode would be required to overcome the emergence of drug-resistance through the accumulation of mutations. This might involve inhibiting the formation of the dimer itself. Here, the folding of HIV-1 PR dimer is studied with several simulation models appropriate for folding mechanism studies. Simulations with an off-lattice Gō-model, which corresponds to a perfectly funneled energy landscape, indicate that the enzyme is formed by association of structured monomers. All-atom molecular dynamics simulations strongly support the stability of an isolated monomer. The conjunction of results from a model that focuses on the protein topology and a detailed all-atom force-field model suggests, in contradiction to some reported equilibrium denaturation experiments, that monomer folding and dimerization are decoupled. The simulation result is, however, in agreement with the recent NMR detection of folded monomers of HIV-1 PR mutants with a destabilized interface. Accordingly, the design of dimerization inhibitors should not focus only on the flexible N and C termini that constitute most of the dimer interface, but also on other structured regions of the monomer. In particular, the relatively high phi values for residues 23-35 and 79-87 in both the folding and binding transition states, together with their proximity to the interface, highlight them as good targets for inhibitor design.

Kinetics of the dimerization of retroviral proteases: the "fireman's grip" and dimerization

All retroviral proteases belong to the family of aspartic proteases. They are active as homodimers, each unit contributing one catalytic aspartate to the active site dyad. An important feature of all aspartic proteases is a conserved complex scaffold of hydrogen bonds supporting the active site, called the "fireman's grip," which involves the hydroxyl groups of two threonine (serine) residues in the active site Asp-Thr(Ser)-Gly triplets. It was shown previously that the fireman's grip is indispensable for the dimer stability of HIV protease. The retroviral proteases harboring Ser in their active site triplet are less active and, under natural conditions, are expressed in higher enzyme/substrate ratio than those having Asp-Thr-Gly triplet. To analyze whether this observation can be attributed to the different influence of Thr or Ser on dimerization, we prepared two pairs of the wild-type and mutant proteases from HIV and myeloblastosis-associated virus harboring either Ser or Thr in their Asp-Thr(Ser)-Gly triplet. The equilibrium dimerization constants differed by an order of magnitude within the relevant pairs. The proteases with Thr in their active site triplets were found to be approximately 10 times more thermodynamically stable. The dimer association contributes to this difference more than does the dissociation. We propose that the fireman's grip might be important in the initial phases of dimer formation to help properly orientate the two subunits of a retroviral protease. The methyl group of threonine might contribute significantly to fixing such an intermediate conformation.

Three-dimensional structure of a monomeric form of a retroviral protease

The assembly of Mason-Pfizer monkey virus Gag polyproteins into immature capsids and their cleavage by the encoded protease are temporally and spatially separated processes, making the virus a particularly useful model for investigation of protease activation. Here we present a high resolution NMR structure of a fully folded monomer of a 12 kDa M-PMV protease (wt 12 PR) and of a Cys7Ala/Asp26Asn/Cys106Ala mutant (12 PR(D26N/C7A/C106A)). The overall structures of both wt 12 PR and 12 PR(D26N/C7A/C106A) follow the conservative structural motif of other retroviral proteases. The most prominent difference from the canonical fold of retroviral proteases is the absence of the interfacial beta-sheet, which leads to the loss of the principal force stabilizing the dimer of M-PMV PR. The monomer-dimer equilibrium can be shifted in favor of the dimer by adding a substrate or an inhibitor, partially compensating for the missing role of the beta-sheet. We also show that cysteines C7 and C106 play a crucial role in stabilizing the dimer and consequently increasing the proteolytic activity of M-PMV PR. This is consistent with the role of reversible oxidative modification of the cysteine residues in the regulation of the maturation of assembled M-PMV capsids in the cytoplasm.

Solution structure of the mature HIV-1 protease monomer: insight into the tertiary fold and stability of a precursor

We present the first solution structure of the HIV-1 protease monomer spanning the region Phe1-Ala95 (PR1-95). Except for the terminal regions (residues 1-10 and 91-95) that are disordered, the tertiary fold of the remainder of the protease is essentially identical to that of the individual subunit of the dimer. In the monomer, the side chains of buried residues stabilizing the active site interface in the dimer, such as Asp25, Asp29, and Arg87, are now exposed to solvent. The flap dynamics in the monomer are similar to that of the free protease dimer. We also show that the protease domain of an optimized precursor flanked by 56 amino acids of the N-terminal transframe region is predominantly monomeric, exhibiting a tertiary fold that is quite similar to that of PR1-95 structure. This explains the very low catalytic activity observed for the protease prior to its maturation at its N terminus as compared with the mature protease, which is an active stable dimer under identical conditions. Adding as few as 2 amino acids to the N terminus of the mature protease significantly increases its dissociation into monomers. Knowledge of the protease monomer structure and critical features of its dimerization may aid in the screening and design of compounds that target the protease prior to its maturation from the Gag-Pol precursor.