Interactions of a novel inhibitor from an extremophilic Bacillus sp. with HIV-1 protease: implications for the mechanism of inactivation

The active site cleft of the HIV-1 protease (PR) is bound by two identical conformationally mobile loops known as flaps, which are important for substrate binding and catalysis. The present article reports, for the first time, an HIV-1 PR inhibitor, ATBI, from an extremophilic Bacillus sp. The inhibitor is found to be a hydrophilic peptide with Mr of 1147, and an amino acid sequence of Ala-Gly-Lys-Lys-Asp-Asp-Asp-Asp-Pro-Pro-Glu. Sequence homology exhibited no similarity with the reported peptidic inhibitors of HIV-1 PR. Investigation of the kinetics of the enzyme-inhibitor interactions revealed that ATBI is a noncompetitive and tight binding inhibitor with the IC(50) and K(i) values 18.0 and 17.8 nm, respectively. The binding of the inhibitor with the enzyme and the subsequent induction of the localized conformational changes in the flap region of the HIV-1 PR were monitored by exploiting the intrinsic fluorescence of the surface exposed Trp-42 residues, which are present at the proximity of the flaps. We have demonstrated by fluorescence and circular dichroism studies that ATBI binds in the active site of the HIV-1 PR and thereby leads to the inactivation of the enzyme. Based on our results, we propose that the inactivation is due to the reorganization of the flaps impairing its flexibility leading toward inaccessibility of the substrate to the active site of the enzyme.     

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

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

Mechanism of substrate recognition by drug-resistant human immunodeficiency virus type 1 protease variants revealed by a novel structural intermediate

Human immunodeficiency virus type 1 (HIV-1) protease processes and cleaves the Gag and Gag-Pol polyproteins, allowing viral maturation, and therefore is an important target for antiviral therapy. Ligand binding occurs when the flaps open, allowing access to the active site. This flexibility in flap geometry makes trapping and crystallizing structural intermediates in substrate binding challenging. In this study, we report two crystal structures of two HIV-1 protease variants bound with their corresponding nucleocapsid-p1 variant. One of the flaps in each of these structures exhibits an unusual "intermediate" conformation. Analysis of the flap-intermediate and flap-closed crystal structures reveals that the intermonomer flap movements may be asynchronous and that the flap which wraps over the P3 to P1 (P3-P1) residues of the substrate might close first. This is consistent with our hypothesis that the P3-P1 region is crucial for substrate recognition. The intermediate conformation is conserved in both the wild-type and drug-resistant variants. The structural differences between the variants are evident only when the flaps are closed. Thus, a plausible structural model for the adaptability of HIV-1 protease to recognize substrates in the presence of drug-resistant mutations has been proposed.     

Effect of the active site D25N mutation on the structure, stability, and ligand binding of the mature HIV-1 protease

All aspartic proteases, including retroviral proteases, share the triplet DTG critical for the active site geometry and catalytic function. These residues interact closely in the active, dimeric structure of HIV-1 protease (PR). We have systematically assessed the effect of the D25N mutation on the structure and stability of the mature PR monomer and dimer. The D25N mutation (PR(D25N)) increases the equilibrium dimer dissociation constant by a factor >100-fold (1.3 +/- 0.09 microm) relative to PR. In the absence of inhibitor, NMR studies reveal clear structural differences between PR and PR(D25N) in the relatively mobile P1 loop (residues 79-83) and flap regions, and differential scanning calorimetric analyses show that the mutation lowers the stabilities of both the monomer and dimer folds by 5 and 7.3 degrees C, respectively. Only minimal differences are observed in high resolution crystal structures of PR(D25N) complexed to darunavir (DRV), a potent clinical inhibitor, or a non-hydrolyzable substrate analogue, Ac-Thr-Ile-Nle-r-Nle-Gln-Arg-NH(2) (RPB), as compared with PR.DRV and PR.RPB complexes. Although complexation with RPB stabilizes both dimers, the effect on their T(m) is smaller for PR(D25N) (6.2 degrees C) than for PR (8.7 degrees C). The T(m) of PR(D25N).DRV increases by only 3 degrees C relative to free PR(D25N), as compared with a 22 degrees C increase for PR.DRV, and the mutation increases the ligand dissociation constant of PR(D25N).DRV by a factor of approximately 10(6) relative to PR.DRV. These results suggest that interactions mediated by the catalytic Asp residues make a major contribution to the tight binding of DRV to PR.     

Molecular dynamics studies on HIV-1 protease: a comparison of the flap motions between wild type protease and the M46I/G51D double mutant

The emergence of drug-resistant mutants of HIV-1 is a tragic effect associated with conventional long-treatment therapies against acquired immunodeficiency syndrome. These mutations frequently involve the aspartic protease encoded by the virus; knowledge of the molecular mechanisms underlying the conformational changes of HIV-1 protease mutants may be useful in developing more effective and longer lasting treatment regimes. The flap regions of the protease are the target of a particular type of mutations occurring far from the active site. These mutations modify the affinity for both substrate and ligands, thus conferring resistance. In this work, molecular dynamics simulations were performed on a native wild type HIV-1 protease and on the drug-resistant M46I/G51D double mutant. The simulation was carried out for a time of 3.5 ns using the GROMOS96 force field, with implementation of the SPC216 explicit solvation model. The results show that the flaps may exist in an ensemble of conformations between a “closed” and an “open” conformation. The behaviour of the flap tips during simulations is different between the native enzyme and the mutant. The mutation pattern leads to stabilization of the flaps in a semi-open configuration.     

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.     

Insight into the structural similarity between HIV protease and secreted aspartic protease-2 and binding mode analysis of HIV-Candida albicans inhibitors

The analysis of the structural similarity between Candida albicans Sap2 and HIV-1 aspartic proteases by molecular modeling gave insight into the common requirements for inhibition of both targets. Structure superimposition of Sap2 and HIV-1 protease confirmed the similarity between their active sites and flap regions. HIV-1 protease inhibitors herein investigated can fit the active site of Sap2, adopting very similar ligand-backbone conformations. In particular, key anchoring sites consisting of Gly85 in Sap2 and Ile50 in HIV-1 protease, both belonging to their corresponding flap regions, were found as elements of a similar binding-mode interaction. The knowledge of the molecular basis for binding to both Sap2 and HIV-1 proteases may ultimately lead to the development of single inhibitor acting on both targets.     

Structural and evolutionary relationships between retroviral and eucaryotic aspartic proteinases.

Three-dimensional crystal structures of the homologous retroviral proteinases from Rous sarcoma virus (RSV PR) and from human immunodeficiency virus (HIV-1 PR) are to a large extent similar and bear close resemblance to the six known structures of the bilobal fungal and mammalian aspartic proteinases. Systematic three-dimensional structural superpositions were carried out between the retroviral and the eucaryotic aspartic proteinases. Both retroviral enzymes were found to be similarly related to their fungal and mammalian counterparts. The most strongly conserved parts correspond to those regions in the N- and C-domains of the eucaryotic enzymes that are related by the interdomain dyad and consist of a combination of secondary structural elements that form the psi loop-alpha helix motif at the active sites of the aspartic proteinases. The retroviral proteinase monomer exhibits nearly the same degree of structural equivalence to the N- and C-domains of the eucaryotic enzymes. In light of the deduced structural relationships between HIV-1 PR and RSV PR, sequence alignments were performed for a number of retroviral proteinases from different subfamilies. There are three highly conserved amino acid sequence stretches, of which two belong to the psi loop-alpha helix motif, that bear moderate sequence similarity with the eucaryotic enzymes. The third conserved sequence stretch among the retroviral proteinases belongs to the flap and bears no resemblance to the flap sequences in the eucaryotic enzymes. The interdomain antiparallel beta sheet in the cellular enzymes differs from the intersubunit beta sheet in the number, arrangement, and directionality of strands, suggesting the possibility of convergent evolution.(ABSTRACT TRUNCATED AT 250 WORDS)     

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

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

HIV-1 protease with 20 mutations exhibits extreme resistance to clinical inhibitors through coordinated structural rearrangements

The escape mutant of HIV-1 protease (PR) containing 20 mutations (PR20) undergoes efficient polyprotein processing even in the presence of clinical protease inhibitors (PIs). PR20 shows >3 orders of magnitude decreased affinity for PIs darunavir (DRV) and saquinavir (SQV) relative to PR. Crystal structures of PR20 crystallized with yttrium, substrate analogue p2-NC, DRV, and SQV reveal three distinct conformations of the flexible flaps and diminished interactions with inhibitors through the combination of multiple mutations. PR20 with yttrium at the active site exhibits widely separated flaps lacking the usual intersubunit contacts seen in other inhibitor-free dimers. Mutations of residues 35-37 in the hinge loop eliminate interactions and perturb the flap conformation. Crystals of PR20/p2-NC contain one uninhibited dimer with one very open flap and one closed flap and a second inhibitor-bound dimer in the closed form showing six fewer hydrogen bonds with the substrate analogue relative to wild-type PR. PR20 complexes with PIs exhibit expanded S2/S2' pockets and fewer PI interactions arising from coordinated effects of mutations throughout the structure, in agreement with the strikingly reduced affinity. In particular, insertion of the large aromatic side chains of L10F and L33F alters intersubunit interactions and widens the PI binding site through a network of hydrophobic contacts. The two very open conformations of PR20 as well as the expanded binding site of the inhibitor-bound closed form suggest possible approaches for modifying inhibitors to target extreme drug-resistant HIV.     

HIV-1 protease molecular dynamics of a wild-type and of the V82F/I84V mutant: possible contributions to drug resistance and a potential new target site for drugs

The protease from type 1 human immunodeficiency virus (HIV-1) is a critical drug target against which many therapeutically useful inhibitors have been developed; however, the set of viral strains in the population has been shifting to become more drug-resistant. Because indirect effects are contributing to drug resistance, an examination of the dynamic structures of a wild-type and a mutant could be insightful. Consequently, this study examined structural properties sampled during 22 nsec, all atom molecular dynamics (MD) simulations (in explicit water) of both a wild-type and the drug-resistant V82F/I84V mutant of HIV-1 protease. The V82F/I84V mutation significantly decreases the binding affinity of all HIV-1 protease inhibitors currently used clinically. Simulations have shown that the curling of the tips of the active site flaps immediately results in flap opening. In the 22-nsec MD simulations presented here, more frequent and more rapid curling of the mutant's active site flap tips was observed. The mutant protease's flaps also opened farther than the wild-type's flaps did and displayed more flexibility. This suggests that the effect of the mutations on the equilibrium between the semiopen and closed conformations could be one aspect of the mechanism of drug resistance for this mutant. In addition, correlated fluctuations in the active site and periphery were noted that point to a possible binding site for allosteric inhibitors.     

Computational design of novel fullerene analogues as potential HIV-1 PR inhibitors: Analysis of the binding interactions between fullerene inhibitors and HIV-1 PR residues using 3D QSAR, molecular docking and molecular dynamics simulations

A series of experimentally reported as well as computationally designed monoadducts and bisadducts of [60]fullerene analogues have been used in order to analyze the binding interactions between fullerene based inhibitors and HIV-1 PR employing docking studies. MD simulations of ligand-free and the inhibitor bound HIV-1 PR systems complemented the above studies and provided proper input structure of HIV-1 PR in docking simulations. The obtained results revealed a different orientation of the beta-hairpin flaps at these two systems. In inhibitor bound system, the flaps of the enzyme are pulled in toward the bottom of the active site (the closed form) while, in ligand-free system flaps shifted away from the dual Asp25 catalytic site and this system adopts a semi-open form. The structural analysis of these systems at catalytic and flexible flap regions of the HIV-1 PR through the simulation, assisted in understanding the structural preferences of these regions, as well as, the adopted orientations of fullerene derivatives within the active site of the enzyme. Five different combinations of steroelectronic fields of 3D QSAR/CoMSIA models were obtained from the set of biologically evaluated and computationally designed fullerene derivatives (training set=43, test set=6) in order to predict novel compounds with improved inhibition effect. The best 3D QSAR/CoMSIA model yielded a cross validated r(2) value of 0.739 and a non-cross validated r(2) value of 0.993. The derived model indicated the importance of steric (42.6%), electrostatic (12.7%), H-bond donor (16.7%) and H-bond acceptor (28.0%) contributions. The derived contour plots together with de novo drug design were then used as pilot models for proposing the novel analogues with enhanced binding affinities. Such structures may trigger the interest of medicinal chemists for novel HIV-1 PR inhibitors possessing higher bioactivity.     

Studies on adaptability of binding residues and flap region of TMC-114 resistance HIV-1 protease mutants

Drug resistant mutations have severely restricted the success of HIV therapy. These mutations frequently involve the aspartic protease encoded by the virus. Knowledge of the molecular mechanisms underlying the conformational changes of HIV-1 protease mutants may be useful in developing more effective and longer lasting treatment regimes. The flap regions of the protease are the target of a particular type of mutations occurring far from the active site, which are able to produce significant resistance against the anti-HIV drug TMC-114. We provide insight into the molecular basis of TMC-114 resistance major flap mutations (I50V and I54M) in HIV-1 protease. It reports the shape complementarity and receptor-ligand interaction analysis supported by unrestrained all-atom molecular dynamics simulations of wild and major flap mutants of HIV-1 protease that sample large conformational changes of the flaps and active site binding residues. Both resistant flap mutants showed less atomic interaction toward TMC-114 and more structural deviation compared to wild HIV-protease. It is due to increasing flexibility at TMC-114 binding cavity and deviation of binding residues in 3-D space. Distortion in binding cavity and deviation in binding residues are the result of alteration in hydrogen bonding. Flap region also exhibited similar behaviour due to changes in number of hydrogen bonds during simulations.     

HIV-1 protease substrate binding and product release pathways explored with coarse-grained molecular dynamics

We analyze the encounter of a peptide substrate with the native HIV-1 protease, the mechanism of substrate incorporation in the binding cleft, and the dissociation of products after substrate hydrolysis. To account for the substrate, we extend a coarse-grained model force field, which we previously developed to study the flap opening dynamics of HIV-1 protease on a microsecond timescale. Molecular and Langevin dynamics simulations show that the flaps need to open for the peptide to bind and that the protease interaction with the substrate influences the flap opening frequency and interval. On the other hand, release of the products does not require flap opening because they can slide out from the binding cleft to the sides of the enzyme. Our data show that in the protease-substrate complex the highest fluctuations correspond to the 17- and 39-turns and the substrate motion is anticorrelated with the 39-turn. Moreover, the active site residues and the flap tips move in phase with the peptide. We suggest some mechanistic principles for how the flexibility of the protein may be involved in ligand binding and release.     

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.     

HIV protease (HIV PR) inhibitor structure-activity-selectivity, and active site molecular modeling of high affinity Leu [CH(OH)CH2]Val modified viral and nonviral substrate analogs.

This report details the structure-activity relationships of the HIV gag substrate analog Val-Ser-Gln-Asn-Leu psi[CH(OH)CH2]Val-Ile-Val (U-85548E), an inhibitor exhibiting subnanomolar affinity towards HIV type-1 aspartic proteinase (HIV-1 PR). Our data show that the P1-P2' tripeptidyl sequence provides the minimal chemical determinant for HIV-1 PR binding. We describe the structure-activity properties of Leu psi[CH(OH)CH2]Val substitution in other peptidyl ligands of nonviral substrate origin (e.g., angiotensinogen, insulin and pepstatin). Furthermore, the aspartic proteinase selectivities of a few key compounds are summarized relative to evaluation against human renin, human pepsin, and the fungal enzyme, rhizopuspepsin. These studies have led to the rational design of nanomolar potent inhibitors of both HIV-1 and HIV-2 PR. Finally, a 2.5 A resolution X-ray crystallographic structure of U-85548E complexed to synthetic HIV-1 PR dimer (Jaskolski et al., Biochemistry 30, 1600 [1991]) provided a 3-D picture of the inhibitor bound to the enzyme active site, and we performed computer-assisted molecular modeling studies to explore the possible binding modes of the above series of Leu psi[CH(OH)CH2]Val substituted HIV-1 PR inhibitors.     

Mutations that alter the activity of the Rous sarcoma virus protease.

Mutations designed by analysis of the Rous sarcoma virus (RSV) and human immunodeficiency virus (HIV)-1 protease (PR) crystal structures were introduced into 1) the substrate binding pocket, 2) the substrate enclosing "flaps," and 3) surface loops of RSV PR. Each mutant PR was expressed in Escherichia coli. Changes in activity were detected by following cleavage of a truncated (NC-PR) precursor polypeptide in E. coli and cleavage of synthetic peptide substrates representing RSV and HIV-1 PR cleavage sites in vitro. Mutations in the substrate binding pocket exchanged amino acid residues located close to the substrate in the HIV-1 PR for structurally equivalent residues in the RSV PR. Changing histidine 65 to glycine (H65G) gave an inactive enzyme, while a double mutant R105P,G106V, as well as the triple mutant, H65G,R105P,G106V, produced enzymes which showed significant activity toward a substrate that represented a HIV-1 cleavage site. Mutating the catalytic aspartate (D37S) or an adjacent conserved alanine to threonine (A40T), produced inactive enzymes. In contrast, the substitution A40S was active, but showed a reduced rate of catalysis. Mutations in the flaps of conserved glycines (G69L, G70L) produced inactive PRs. Two extended RSV PR surface loops were shortened to the size found in HIV-1 PR and resulted in drastically reduced activity. These results have confirmed some of the basic predictions made from structural models but have also revealed unexpected roles and interactions in the protein.     

Probing the S1/S1' substrate binding pocket geometry of HIV-1 protease with modified aspartic acid analogues

Aspartates 25 and 125, the active site residues of HIV-1 protease, participate functionally in proteolysis by what is believed to be a general acid-general base mechanism. However, the structural role that these residues may play in the formation and maintenance of the neighboring S1/S1' substrate binding pockets remains largely unstudied. Because the active site aspartic acids are essential for catalysis, alteration of these residues to any other naturally occurring amino acid by conventional site-directed mutagenesis renders the protease inactive, and hence impossible to characterize functionally. To investigate whether Asp-25 and Asp-125 may also play a structural role that influences substrate processing, a series of active site protease mutants has been produced in a cell-free protein synthesizing system via readthrough of mRNA nonsense (UAG) codons by chemically misacylated suppressor tRNAs. The suppressor tRNAs were activated with the unnatural aspartic acid analogues erythro-beta-methylaspartic acid, threo-beta-methylaspartic acid, or beta,beta-dimethylaspartic acid. On the basis of the specific activity measurements of the mutants that were produced, the introduction of the beta-methyl moiety was found to alter protease function to varying extents depending upon its orientation. While a beta-methyl group in the erythro orientation was the least deleterious to the specific activity of the protease, a beta-methyl group in the threo orientation, present in the modified proteins containing threo-beta-methylaspartate and beta,beta-dimethylaspartate, resulted in specific activities between 0 and 45% of that of the wild type depending upon the substrate and the substituted active site position. Titration studies of pH versus specific activity and inactivation studies, using an aspartyl protease specific suicide inhibitor, demonstrated that the mutant proteases maintained bell-shaped pH profiles, as well as suicide-inhibitor susceptibilities that are characteristic of aspartyl proteases. A molecular dynamics simulation of the beta-substituted aspartates in position 25 of HIV-1 protease indicated that the threo-beta-methyl moiety may partially obstruct the adjacent S1' binding pocket, and also cause reorganization within the pocket, especially with regard to residues Val-82 and Ile-84. This finding, in conjunction with the biochemical studies, suggests that the active site aspartate residues are in proximity to the S1/S1' binding pocket and may be spatially influenced by the residues presented in these pockets upon substrate binding. It thus seems possible that the catalytic residues cooperatively interact with the residues that constitute the S1/S1' binding pockets and can be repositioned during substrate binding to orient the active site carboxylates with respect to the scissile amide bond, a process that likely affects the facility of proteolysis.     

A solution NMR study of the binding kinetics and the internal dynamics of an HIV-1 protease-substrate complex

NMR studies of the binding of a substrate to an inactive HIV-1 protease construct, containing an active site mutation PR(D25N), are reported. Substrate titration measurements monitored by HSQC spectra and a (15)N-edited NOESY experiment show that the chromogenic substrate analog of the capsid/p2 cleavage site binds to PR(D25N) with an equilibrium dissociation constant, K(D), of 0.27 +/- 0.05 mM, and upper limits of the association and dissociation rate constants, 2 x 10(4) M(-1)s(-1) and 10 s(-1), respectively, at 20 degrees C, pH 5.8. This association rate constant is not in the diffusion limit, suggesting that association is controlled by a rare event, such as opening of the protease flaps. Analysis of (15)N relaxation experiments reveals a slight reduction of S(2) values in the flap region, indicating a small increase in the amplitude of internal motion on the sub-nsec timescale. In addition, several residues in the flap region are mobile on the conformational exchange timescale, msec-microsec. Flap dynamics of the protease-substrate complex are compared with those of protease-inhibitor complexes, and the implications of these results for substrate-binding models are discussed.     

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

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