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.     

Insight into the folding inhibition of the HIV-1 protease by a small peptide

It has recently been shown that the highly protected segments 24-34 (S2) and 83-93 (S8) of each of the two 99-mers of human immunodeficiency virus type 1 protease play an essential role in the folding of the monomers, giving rise to the so-called (postcritical) folding nucleus ((FN) minimum condensation unit ensuring folding) when they dock. This scenario received further support from model calculations that demonstrated that the peptide p-S8, displaying an amino acid sequence identical to the corresponding (83-93) segment of the monomer, can be used to interfere with the formation of the FN and eventually to inhibit folding by docking the fragment 24-34. Experiments in vitro and in cells infected with ex vivo wild-type and multiresistant HIV isolates confirm that the inhibition power of p-S8 is robust. On the other hand, there is no direct evidence demonstrating the validity of the proposed mechanism of inhibition associated with p-S8. To shed light on this question and to provide the basis for the design of a molecule mimetic to p-S8, to be used as lead of an eventual drug against AIDS, we study, in this paper, with the help of all-atom simulations in explicit solvent and the novel method of metadynamics combined with parallel tempering: a), the free energy and the equilibrium structure of each of the peptides p-S2 and p-S8; b), the details of the docking mechanism of the two peptides and the free energy associated with this process. Whereas p-S8 is found to be well structured, p-S2 is rather flexible, wrapping itself around p-S8 to give rise to the FN, which is stabilized by three particular hydrogen bonds.     

Disruption of the HIV-1 protease dimer with interface peptides: structural studies using NMR spectroscopy combined with [2-(13)C]-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.     

Litchi chinensis-derived terpenoid as anti-HIV-1 protease agent: structural design from molecular dynamics simulations

The molecular structures of the binding between human immunodeficiency virus-1 protease (HIV-1PR) and various inhibitors including existing extensive natural products extracts have been investigated for anti-HIV drug development. In this study, the binding of HIV-1PR and a terpenoid from Litchi chinensis extracts (3-oxotrirucalla-7,24-dien-21-oic acid) was investigated in order to clarify the inhibition effectiveness of this compound. Molecular dynamics (MD) simulations of HIV-1PR complex with 3-oxotrirucalla-7,24-dien-21-oic acid were performed including water molecules. The MD simulation results indicated the formation of hydrogen bonds between the oxygen atoms of the inhibitor and the catalytic aspartates, which are commonly found in inhibitors–protease complexes. On the other hand, no hydrogen bonding of this particular inhibitor to the flap region was found. In addition, the radial distribution function of water oxygens around the catalytic carboxylate nitrogens of Asp29 and Asp30 suggests that at least one or two water molecules are in the active site region whereas direct interaction of the inhibitor was found for catalytic carboxylate oxygen of Asp25. The results of this simulation, in comparison with the structures of other HIV-PR inhibitor complexes, could lead to a better understanding of the activity of 3-oxotrirucalla-7,24-dien-21-oic acid.     

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.     

Structural evidence for effectiveness of darunavir and two related antiviral inhibitors against HIV-2 protease

No drug has been targeted specifically for HIV-2 (human immunodeficiency virus type 2) infection despite its increasing prevalence worldwide. The antiviral HIV-1 (human immunodeficiency virus type 1) protease (PR) inhibitor darunavir and the chemically related GRL98065 and GRL06579A were designed with the same chemical scaffold and different substituents at P2 and P2′ to optimize polar interactions for HIV-1 PR (PR1). These inhibitors are also effective antiviral agents for HIV-2-infected cells. Therefore, crystal structures of HIV-2 PR (PR2) complexes with the three inhibitors have been solved at 1.2-Å resolution to analyze the molecular basis for their antiviral potency. Unusually, the crystals were grown in imidazole and zinc acetate buffer, which formed interactions with the PR2 and the inhibitors. Overall, the structures were very similar to the corresponding inhibitor complexes of PR1 with an RMSD of 1.1 Å on main-chain atoms. Most hydrogen-bond and weaker C-H…O interactions with inhibitors were conserved in the PR2 and PR1 complexes, except for small changes in interactions with water or disordered side chains. Small differences were observed in the hydrophobic contacts for the darunavir complexes, in agreement with relative inhibition of the two PRs. These near-atomic-resolution crystal structures verify the inhibitor potency for PR1 and PR2 and will provide the basis for the development of antiviral inhibitors targeting PR2.     

Insights from atomic-resolution X-ray structures of chemically synthesized HIV-1 protease in complex with inhibitors

The human immunodeficiency virus 1 (HIV-1) protease (PR) is an aspartyl protease essential for HIV-1 viral infectivity. HIV-1 PR has one catalytic site formed by the homodimeric enzyme. We chemically synthesized fully active HIV-1 PR using modern ligation methods. When complexed with the classic substrate-derived inhibitors JG-365 and MVT-101, the synthetic HIV-1 PR formed crystals that diffracted to 1.04- and 1.2-A resolution, respectively. These atomic-resolution structures revealed additional structural details of the HIV-1 PR's interactions with its active site ligands. Heptapeptide inhibitor JG-365, which has a hydroxyethylamine moiety in place of the scissile bond, binds in two equivalent antiparallel orientations within the catalytic groove, whereas the reduced isostere hexapeptide MVT-101 binds in a single orientation. When JG-365 was converted into the natural peptide substrate for molecular dynamic simulations, we found putative catalytically competent reactant states for both lytic water and direct nucleophilic attack mechanisms. Moreover, free energy perturbation calculations indicated that the insertion of catalytic water into the catalytic site is an energetically favorable process.     

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.     

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.     

High resolution crystal structures of HIV-1 protease with a potent non-peptide inhibitor (UIC-94017) active against multi-drug-resistant clinical strains

The compound UIC-94017 (TMC-114) is a second-generation HIV protease inhibitor with improved pharmacokinetics that is chemically related to the clinical inhibitor amprenavir. UIC-94017 is a broad-spectrum potent inhibitor active against HIV-1 clinical isolates with minimal cytotoxicity. We have determined the high-resolution crystal structures of UIC-94017 in complexes with wild-type HIV-1 protease (PR) and mutant proteases PR(V82A) and PR(I84V) that are common in drug-resistant HIV. The structures were refined at resolutions of 1.10-1.53A. The crystal structures of PR and PR(I84V) with UIC-94017 ternary complexes show that the inhibitor binds to the protease in two overlapping positions, while the PR(V82A) complex had one ordered inhibitor. In all three structures, UIC-94017 forms hydrogen bonds with the conserved main-chain atoms of Asp29 and Asp30 of the protease. These interactions are proposed to be critical for the potency of this compound against HIV isolates that are resistant to multiple protease inhibitors. Other small differences were observed in the interactions of the mutants with UIC-94017 as compared to PR. PR(V82A) showed differences in the position of the main-chain atoms of residue 82 compared to PR structure that better accommodated the inhibitor. Finally, the 1.10A resolution structure of PR(V82A) with UIC-94017 showed an unusual distribution of electron density for the catalytic aspartate residues, which is discussed in relation to the reaction mechanism.     

Structural studies of FIV and HIV-1 proteases complexed with an efficient inhibitor of FIV protease

Three forms of feline immunodeficiency virus protease (FIV PR), the wild type (wt) and two single point mutants, V59I and Q99V, as well as human immunodeficiency virus type 1 protease (HIV-1 PR), were cocrystallized with the C2-symmetric inhibitor, TL-3. The mutants of FIV PR were designed to replace residues involved in enzyme-ligand interactions by the corresponding HIV-1 PR residues at the structurally equivalent position. TL-3 shows decreased (improved) inhibition constants with these FIV PR mutants relative to wt FIV PR. Despite similar modes of binding of the inhibitor to all PRs (from P3 to P3'), small differences are evident in the conformation of the Phe side chains of TL-3 at the P1 and P1' positions in the complexes with the mutated FIV PRs. The differences mimick the observed binding of TL-3 in HIV-1 PR and correlate with a significant improvement in the inhibition constants of TL-3 with the two mutant FIV PRs. Large differences between the HIV-1 and FIV PR complexes are evident in the binding modes of the carboxybenzyl groups of TL-3 at P4 and P4'. In HIV-1 PR:TL-3, these groups bind over the flap region, whereas in the FIV PR complexes, the rings are located along the major axis of the active site. A significant difference in the location of the flaps in this region of the HIV-1 and FIV PRs correlates with the observed conformational changes in the binding mode of the peptidomimetic inhibitor at the P4 and P4' positions. These findings provide a structural explanation of the observed Ki values for TL-3 with the different PRs and will further assist in the development of improved inhibitors.     

Dimerization inhibitors of HIV-1 protease

By targeting the highly conserved antiparallel beta-sheet formed by the interdigitation of the N- and C-terminal strands of each monomer, dimerization inhibitors of HIV-1 protease may be useful to overcome the drug resistance observed with current active-site directed antiproteases. Sequestration of the monomer by the inhibitor (or disruption of the dimer interface) prevents the correct assembly of the inactive monomers to active enzyme. Strategies for the design of drugs targeting the dimer interface are described. Various dimerization inhibitors are reported including N- and C-terminal mimetics, lipopeptides and cross-linked interface peptides.     

Design and synthesis of new inhibitors of HIV-1 protease dimerization with conformationally constrained templates.

To inhibit the HIV-1 protease dimerization necessary to exhibit enzymatic activity, we synthesized and evaluated a new beta-sheet peptide (compound 1), containing 4-(2-aminoethyl)-6-dibenzofuranpropionic acid as a conformationally restricted linker and a non-peptidic beta-strand mimetic, 2-[3-([2-[(9-fluorenylmethoxy)carbonyl]hydrazino]carbonyl)-4-methoxyanilino]-2-oxoacetic acid (Fmoc-Hao-OH, compound 2). Kinetic analysis showed that compound 1 inhibited the dimerization of HIV-1 protease by a dissociative mechanism with a K(id) value of 5.4 microM at 37 degrees C (pH 5.0). However, compound 2 showed a small shift in the slope of the lines in the Zhang-Poorman plot (K(id)=9.1 microM), suggesting that compound 2 inhibits the dimerization of HIV-1 PR not only through a dissociative mechanism but also through an active-site directed mechanism partly. This is the first study of a non-peptidic inhibitor of HIV-1 protease dimerization.     

Free energy calculations on dimer stability of the HIV protease using molecular dynamics and a continuum solvent model

Dimerization of HIV-I protease (HIV PR) monomers is an essential prerequisite for viral proteolytic activity and the subsequent generation of infectious virus particles. Disrupting dimerization of the enzyme can inhibit its activity. We have calculated the relative binding free energies between different dimers of the HIV protease using molecular dynamics and a continuum model, which we call MM/PBSA. We examined the dominant negative inhibition of the HIV PR by a mutated form of the protease and found relative dimerization free energies of homo- and hetero-dimerization consistent with experimental data. We also developed a rapid screening method, which was called the virtual mutagenesis method to consider other mutations which might stabilize non-wild-type heterodimers. Using this approach, we considered the mutations near the dimer interface which might cause dominant negative inhibition of the HIV PR. The rapid method we developed can be used in studying any ligand-protein and protein-protein interaction, in order to identify mutations that can enhance the binding affinities of the complex.     

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

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

Molecular dynamics simulations of HIV-1 protease monomer: Assembly of N-terminus and C-terminus into beta-sheet in water solution

HIV-1 protease (HIV-PR) consists of two identical subunits that are united together through a four-stranded antiparallel beta-sheet formed of the peptide termini of each monomer. Since the active site exists only in the dimer, a strategy that is attracting more and more attention in inhibitor design and which may overcome the serious drug resistance caused by competitive inhibitors is to block the peptide termini of the monomer, thereby interfering with formation of the active dimer. In the present work, we performed several extensive molecular dynamics (MD) simulations of the HIV-PR monomer in water to illustrate its solvated conformation and dynamics behavior. We found that the peptide termini usually assembled into beta-sheet after several nanoseconds' simulation, and became much less flexible. This beta-sheet is stabilized by intramolecular interactions and is not easily disaggregated under the present MD simulation conditions. This transformation may be an important transition during the relaxing and equilibrating of the HIV-PR monomer in aqueous solution, and the terminal beta-sheet may be one of the major conformations of the solvated HIV-PR monomer termini in water. This work may provide new insights into the dynamics behavior and dimerization mechanism of HIV-PR, and more significantly, offer a more rational receptor model for the design and discovery of novel dimerization inhibitors than crystalline structures.     

Yeast two-hybrid assay for examining human immunodeficiency virus protease heterodimer formation with dominant-negative inhibitors and multidrug-resistant variants

The yeast two-hybrid assay was used to study the dimerization of engineered and naturally occurring variants of human immunodeficiency virus (HIV) protease (PR) monomers. Defective monomers that were previously shown to exhibit a dominant-negative (D-N) effect in cultured mammalian cells were tested for their ability to interact in the two-hybrid assay. Similarly, monomers with dimer-interface substitutions and monomers harboring in vivo selected mutations that confer multidrug resistance (mdr) in an AIDS patient were tested for interaction in yeast. Dimer formation between wt monomers with catalytic aspartates was not detected in yeast, whereas the dimerization of PR monomers harboring the acid active site substitution D25N was readily demonstrated. The use of inactive monomers harboring the D25N substitution as a genetic background for studying additional HIV PR mutations allowed for the probing of interactions between monomers with mdr-associated mutations with those based on the HIV-1 HXB2R sequence. The HTLVIII/HIV-1 HXB2R clone has been the basis for a large number of HIV-related plasmids, primers, antibodies, and other specific reagents throughout the HIV research community. The results of our assay suggest that HXB2R-based D-N PR inhibitors associate with variant monomers based on the recently obtained nucleotide sequence from an AIDS patient with a multidrug-resistant virus. These results further encourage the use of D-N PR inhibitors as antiviral agents which may complement existing small-molecule combination therapies.     

Molecular dynamics simulations of the first steps of the reaction catalyzed by HIV-1 protease

The mechanism of the first steps of the reaction catalyzed by HIV-1 protease was studied through molecular dynamics simulations. The potential energy surface in the active site was generated using the approximate valence bond method. The approximate valence bond (AVB) method was parameterized based on density functional calculations. The surrounding protein and explicit water environment was modeled with conventional, classical force field. The calculations were performed based on HIV-1 protease complexed with the MVT-101 inhibitor that was modified to a model substrate. The protonation state of the catalytic aspartates was determined theoretically. Possible reaction mechanisms involving the lytic water molecule are accounted for in this study. The modeled steps include the dissociation of the lytic water molecule and proton transfer onto Asp-125, the nucleophilic attack followed by a proton transfer onto peptide nitrogen. The simulations show that in the active site most preferable energetically are structures consisting of ionized or polarized molecular fragments that are not accounted for in conventional molecular dynamics. The mobility of the lytic water molecule, the dynamics of the hydrogen bond network, and the conformation of the aspartates in the active center were analyzed.     

Illustration of HIV-1 protease folding through a molten-globule-like intermediate using an experimental model that implicates alpha-crystallin and calcium ions

The folding of HIV-1 protease to its active form involves the coordination of structure formation and dimerization, which follows a hierarchy consisting of folding nuclei spanning from the active site, hinge region, and dimerization domain. However, the biochemical characteristics of the folding intermediates of this protein remain to be elucidated. In an experimental model, the denaturation of the tethered dimer of HIV-1 protease by guanidine hydrochloride revealed an alternative conformation resembling the molten-globule state. The molten-globule state binds to the molecular chaperone alpha-crystallin and prevents its aggregation; however, the chaperone alone failed to reconstitute HIV-1 protease into its active form. Calcium ion assisted in the release of active enzyme from the chaperone complex. Alpha-crystallin, a member of the small heat-shock protein, assists proteins to fold correctly; however, the underlying principle of signals responsible for chaperone-mediated protein folding remains enigmatic. X-ray photoelectron spectroscopy has been employed to provide the evidence of calcium binding to alpha-crystallin and to decipher the effect of calcium binding on the chaperone-mediated refolding of HIV-1 protease. On the basis of our spectroscopic data, we propose that calcium ions interact with the carboxyl groups of the surface-exposed acidic amino acids of alpha-crystallin bringing electrostatic interference, which plays a pivotal role in inducing conformational changes in the chaperone responsible for the release of the active enzyme.     

HIV-1 protease folding and the design of drugs which do not create resistance

Human immunodeficiency virus type 1 (HIV-1) protease (PR) plays an essential role in the life cycle of the virus. Consequently, its inhibition can control acquired immunodeficiency syndrome (AIDS). Any pharmacological treatment targeting the active site of the protease is known to generate escape mutants. On the other hand, if a drug targets a site crucial for the correct folding of the protease, mutations affecting this region would denaturate the protein and thus will not be expressed. We review the progress in our understanding of the folding of the protease, which has been instrumental in the design of a (non-conventional) folding inhibitor. The transferability of these results to other proteins testify to the universality of the folding-inhibition scenario for the design of leads of drugs which are unlikely to generate resistance.