Conserved hydrogen bonds and water molecules in MDR HIV-1 protease substrate complexes

The success of highly active antiretroviral therapy (HAART) in anti-HIV therapy is severely compromised by the rapidly developing drug resistance. HIV-1 protease inhibitors, part of HAART, are losing their potency and efficacy in inhibiting the target. Multi-drug resistant (MDR) 769 HIV-1 protease (resistant mutations at residues 10, 36, 46, 54, 62, 63, 71, 82, 84, 90) was selected for the present study to understand the binding to its natural substrates. The nine crystal structures of MDR769 HIV-1 protease substrate hepta-peptide complexes were analyzed in order to reveal the conserved structural elements for the purpose of drug design against MDR HIV-1 protease. Our structural studies demonstrated that highly conserved hydrogen bonds between the protease and substrate peptides, together with the conserved crystallographic water molecules, played a crucial role in the substrate recognition, substrate stabilization and protease stabilization. In addition, the absence of the key flap-ligand bridging water molecule might imply a different catalytic mechanism of MDR769 HIV-1 protease compared to that of wild type (WT) HIV-1 protease.     

Pentacycloundecane derived hydroxy acid peptides: a new class of irreversible non-scissile ether bridged type isoster as potential HIV-1 wild type C-SA protease inhibitors

Novel peptides incorporating the PCU derived hydroxy acid (5-hydroxy-4-oxahexacyclo[5.4.1.0(2,6).0(3,10).0(5,9).0(8,11)]dodecane) were synthesized and their activity against the resistance-prone wild type C-South African (C-SA) HIV-protease is reported. The attachment of peptides and peptoids to the PCU derived hydroxy acid resulted in a series of structurally diverse promising HIV-1 protease inhibitors. Amongst the nine novel compounds, 16, 17, 20 and 23 gave IC(50) values ranging from 0.6 to 5.0 μM against the wild type C-SA HIV-1 protease enzyme. Docking studies and molecular dynamic (MD) simulations have been carried out in order to understand the binding mode of the PCU moiety at the active site of the HIV protease enzyme. A conserved hydrogen bonding pattern between the PCU derived hydroxy ether and the active site residues, ASP25/ASP25', was observed in all active compounds.     

Extreme Entropy-Enthalpy Compensation in a Drug-Resistant Variant of HIV-1 Protease

The development of HIV-1 protease inhibitors has been the historic paradigm of rational structure-based drug design, where structural and thermodynamic analyses have assisted in the discovery of novel inhibitors. While the total enthalpy and entropy change upon binding determine the affinity, often the thermodynamics are considered in terms of inhibitor properties only. In the current study, profound changes are observed in the binding thermodynamics of a drug-resistant variant compared to wild-type HIV-1 protease, irrespective of the inhibitor bound. This variant (Flap+) has a combination of flap and active site mutations and exhibits extremely large entropy-enthalpy compensation compared to wild-type protease, 5-15 kcal/mol, while losing only 1-3 kcal/mol in total binding free energy for any of six FDA-approved inhibitors. Although entropy-enthalpy compensation has been previously observed for a variety of systems, never have changes of this magnitude been reported. The co-crystal structures of Flap+ protease with four of the inhibitors were determined and compared with complexes of both the wild-type protease and another drug-resistant variant that does not exhibit this energetic compensation. Structural changes conserved across the Flap+ complexes, which are more pronounced for the flaps covering the active site, likely contribute to the thermodynamic compensation. The finding that drug-resistant mutations can profoundly modulate the relative thermodynamic properties of a therapeutic target independent of the inhibitor presents a new challenge for rational drug design.     

Design and synthesis of highly potent HIV-1 protease inhibitors with novel isosorbide-derived P2 ligands

The design, synthesis, and biological evaluation of a series of six HIV-1 protease inhibitors incorporating isosorbide moiety as novel P2 ligands are described. All the compounds are very potent HIV-1 protease inhibitors with IC₅₀ values in the nanomolar or picomolar ranges (0.05–0.43 nM). Molecular docking studies revealed the formation of an extensive hydrogen-bonding network between the inhibitor and the active site. Particularly, the isosorbide-derived P2 ligand is involved in strong hydrogen bonding interactions with the backbone atoms.     

Strength of hydrogen bond network takes crucial roles in the dissociation process of inhibitors from the HIV-1 protease binding pocket

To understand the underlying mechanisms of significant differences in dissociation rate constant among different inhibitors for HIV-1 protease, we performed steered molecular dynamics (SMD) simulations to analyze the entire dissociation processes of inhibitors from the binding pocket of protease at atomistic details. We found that the strength of hydrogen bond network between inhibitor and the protease takes crucial roles in the dissociation process. We showed that the hydrogen bond network in the cyclic urea inhibitors AHA001/XK263 is less stable than that of the approved inhibitor ABT538 because of their large differences in the structures of the networks. In the cyclic urea inhibitor bound complex, the hydrogen bonds often distribute at the flap tips and the active site. In contrast, there are additional accessorial hydrogen bonds formed at the lateral sides of the flaps and the active site in the ABT538 bound complex, which take crucial roles in stabilizing the hydrogen bond network. In addition, the water molecule W301 also plays important roles in stabilizing the hydrogen bond network through its flexible movement by acting as a collision buffer and helping the rebinding of hydrogen bonds at the flap tips. Because of its high stability, the hydrogen bond network of ABT538 complex can work together with the hydrophobic clusters to resist the dissociation, resulting in much lower dissociation rate constant than those of cyclic urea inhibitor complexes. This study may provide useful guidelines for design of novel potent inhibitors with optimized interactions.     

Role of hydroxyl group and R/S configuration of isostere in binding properties of HIV-1 protease inhibitors.

The crystal structure of the complex between human immunodeficiency virus type 1 (HIV-1) protease and a peptidomimetic inhibitor of ethyleneamine type has been refined to R factor of 0.178 with diffraction limit 2.5 A. The peptidomimetic inhibitor Boc-Phe-Psi[CH2CH2NH]-Phe-Glu-Phe-NH2 (denoted here as OE) contains the ethyleneamine replacement of the scissile peptide bond. The inhibitor lacks the hydroxyl group which is believed to mimic tetrahedral transition state of proteolytic reaction and thus is suspected to be necessary for good properties of peptidomimetic HIV-1 protease inhibitors. Despite the missing hydroxyl group the inhibition constant of OE is 1.53 nm and it remains in the nanomolar range also towards several available mutants of HIV-1 protease. The inhibitor was found in the active site of protease in an extended conformation with a unique hydrogen bond pattern different from hydroxyethylene and hydroxyethylamine inhibitors. The isostere nitrogen forms a hydrogen bond to one catalytic aspartate only. The other aspartate forms two weak hydrogen bridges to the ethylene group of the isostere. A comparison with other inhibitors of this series containing isostere hydroxyl group in R or S configuration shows different ways of accommodation of inhibitor in the active site. Special attention is devoted to intermolecular contacts between neighbouring dimers responsible for mutual protein adhesion and for a special conformation of Met46 and Phe53 side chains not expected for free protein in water solution.     

Crystal structures of multidrug-resistant HIV-1 protease in complex with two potent anti-malarial compounds

Two potent inhibitors (compounds 1 and 2) of malarial aspartyl protease, plasmepsin-II, were evaluated against wild type (NL4-3) and multidrug-resistant clinical isolate 769 (MDR) variants of human immunodeficiency virus type-1 (HIV-1) aspartyl protease. Enzyme inhibition assays showed that both 1 and 2 have better potency against NL4-3 than against MDR protease. Crystal structures of MDR protease in complex with 1 and 2 were solved and analyzed. Crystallographic analysis revealed that the MDR protease exhibits a typical wide-open conformation of the flaps (Gly48 to Gly52) causing an overall expansion in the active site cavity, which, in turn caused unstable binding of the inhibitors. Due to the expansion of the active site cavity, both compounds showed loss of direct contacts with the MDR protease compared to the docking models of NL4-3. Multiple water molecules showed a rich network of hydrogen bonds contributing to the stability of the ligand binding in the distorted binding pockets of the MDR protease in both crystal structures. Docking analysis of 1 and 2 showed a decrease in the binding affinity for both compounds against MDR supporting our structure-function studies. Thus, compounds 1 and 2 show promising inhibitory activity against HIV-1 protease variants and hence are good candidates for further development to enhance their potency against NL4-3 as well as MDR HIV-1 protease variants.     

Comprehensive bioinformatic analysis of the specificity of human immunodeficiency virus type 1 protease

Rapidly developing viral resistance to licensed human immunodeficiency virus type 1 (HIV-1) protease inhibitors is an increasing problem in the treatment of HIV-infected individuals and AIDS patients. A rational design of more effective protease inhibitors and discovery of potential biological substrates for the HIV-1 protease require accurate models for protease cleavage specificity. In this study, several popular bioinformatic machine learning methods, including support vector machines and artificial neural networks, were used to analyze the specificity of the HIV-1 protease. A new, extensive data set (746 peptides that have been experimentally tested for cleavage by the HIV-1 protease) was compiled, and the data were used to construct different classifiers that predicted whether the protease would cleave a given peptide substrate or not. The best predictor was a nonlinear predictor using two physicochemical parameters (hydrophobicity, or alternatively polarity, and size) for the amino acids, indicating that these properties are the key features recognized by the HIV-1 protease. The present in silico study provides new and important insights into the workings of the HIV-1 protease at the molecular level, supporting the recent hypothesis that the protease primarily recognizes a conformation rather than a specific amino acid sequence. Furthermore, we demonstrate that the presence of 1 to 2 lysine residues near the cleavage site of octameric peptide substrates seems to prevent cleavage efficiently, suggesting that this positively charged amino acid plays an important role in hindering the activity of the HIV-1 protease.     

Accessory mutations balance the marginal stability of the HIV-1 protease in drug resistance

The HIV-1 protease is a major target of inhibitor drugs in AIDS therapies. The therapies are impaired by mutations of the HIV-1 protease that can lead to resistance to protease inhibitors. These mutations are classified into major mutations, which usually occur first and clearly reduce the susceptibility to protease inhibitors, and minor, accessory mutations that occur later and individually do not substantially affect the susceptibility to inhibitors. Major mutations are predominantly located in the active site of the HIV-1 protease and can directly interfere with inhibitor binding. Minor mutations, in contrast, are typically located distal to the active site. A central question is how these distal mutations contribute to resistance development. In this article, we present a systematic computational investigation of stability changes caused by major and minor mutations of the HIV-1 protease. As most small single-domain proteins, the HIV-1 protease is only marginally stable. Mutations that destabilize the folded, active state of the protease therefore can shift the conformational equilibrium towards the unfolded, inactive state. We find that the most frequent major mutations destabilize the HIV-1 protease, whereas roughly half of the frequent minor mutations are stabilizing. An analysis of protease sequences from patients in treatment indicates that the stabilizing minor mutations are frequently correlated with destabilizing major mutations, and that highly resistant HIV-1 proteases exhibit significant fractions of stabilizing mutations. Our results thus indicate a central role of minor mutations in balancing the marginal stability of the protease against the destabilization induced by the most frequent major mutations.     

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.     

Effects of HIV-1 protease on cellular functions and their potential applications in antiretroviral therapy

ABSTRACT: Human Immunodeficiency Virus Type 1 (HIV-1) protease inhibitors (PIs) are the most potent class of drugs in antiretroviral therapies. However, viral drug resistance to PIs could emerge rapidly thus reducing the effectiveness of those drugs. Of note, all current FDA-approved PIs are competitive inhibitors, i.e., inhibitors that compete with substrates for the active enzymatic site. This common inhibitory approach increases the likelihood of developing drug resistant HIV-1 strains that are resistant to many or all current PIs. Hence, new PIs that move away from the current target of the active enzymatic site are needed. Specifically, allosteric inhibitors, inhibitors that block HIV-1 protease active site, should be sought. Another common feature of current PIs is they were all developed based on the structure-based design. Drugs derived from a structure-based strategy may generate target specific and potent inhibitors. However, this type of drug design can only target one site at a time and drugs discovered by this method are often associated with strong side effects such as cellular toxicity, limiting its number of target choices, efficacy, and applicability. In contrast, a cell-based system may provide a useful alternative strategy that can overcome many of the inherited shortcomings associated with structure-based drug designs. For example, allosteric PIs can be sought using a cell-based system without considering the site or mechanism of inhibition. In addition, a cell-based system can eliminate those PIs that have strong cytotoxic effect. Most importantly, a simple, economical, and easy-to-maintained eukaryotic cellular system such as yeast will allow us to search for potential PIs in a large-scaled high throughput screening (HTS) system, thus increasing the chance of success. Based on our many years of experience in using fission yeast as a model system to study HIV-1 Vpr, we propose the use of fission yeast as a possible surrogate system to study the effects of HIV-1 protease on cellular functions and to explore its utility as a HTS system to search for new PIs to battle HIV-1 strains resistant to the current PI drugs.     

Nine Crystal Structures Determine the Substrate Envelope of the MDR HIV-1 Protease

Under drug selection pressure, emerging mutations render HIV-1 protease drug resistant, leading to the therapy failure in anti-HIV treatment. It is known that nine substrate cleavage site peptides bind to wild type (WT) HIV-1 protease in a conserved pattern. However, how the multidrug-resistant (MDR) HIV-1 protease binds to the substrate cleavage site peptides is yet to be determined. MDR769 HIV-1 protease (resistant mutations at residues 10, 36, 46, 54, 62, 63, 71, 82, 84, and 90) was selected for present study to understand the binding to its natural substrates. MDR769 HIV-1 protease was co-crystallized with nine substrate cleavage site hepta-peptides. Crystallographic studies show that MDR769 HIV-1 protease has an expanded substrate envelope with wide open flaps. Furthermore, ligand binding energy calculations indicate weaker binding in MDR769 HIV-1 protease-substrate complexes. These results help in designing the next generation of HIV-1 protease inhibitors by targeting the MDR HIV-1 protease.     

A major role for a set of non-active site mutations in the development of HIV-1 protease drug resistance

A major problem in the chemotherapy of HIV-1 infection is the appearance of drug resistance. In the case of HIV-1 protease inhibitors, resistance originates from mutations in the protease molecule that lower the affinity of inhibitors while still maintaining a viable enzymatic profile. Drug resistance mutations can be classified as active site or non-active site mutations depending on their location within the protease molecule. Active site mutations directly affect drug/target interactions, and their action can be readily understood in structural terms. Non-active site mutations influence binding from distal locations, and their mechanism of action is not immediately apparent. In this paper, we have characterized a mutant form of the HIV-1 protease, ANAM-11, identified in clinical isolates from HIV-1 infected patients treated with protease inhibitors. This mutant protease contains 11 mutations, 10 of which are located outside the active site (L10I/M36I/S37D/M46I/R57K/L63P/A71V/G73S/L90M/I93L) and 1 within the active site (I84V). ANAM-11 lowers the binding affinity of indinavir, nelfinavir, saquinavir, and ritonavir by factors of 4000, 3300, 5800, and 80000, respectively. Surprisingly, most of the loss in inhibitor affinity is due to the non-active site mutations as demonstrated by additional experiments performed with a protease containing only the 10 non-active site mutations (NAM-10) and another containing only the active site mutation (A-1). Kinetic analysis with two different substrates yielded comparable catalytic efficiencies for A-1, ANAM-11, NAM-10, and the wild-type protease. These studies demonstrate that non-active site mutations can be the primary source of resistance and that their role is not necessarily limited to compensate deleterious effects of active site mutations. Analysis of the structural stability of the proteases by differential scanning calorimetry reveals that ANAM-11 and NAM-10 are structurally more stable than the wild-type protease while A-1 is less stable. Together, the binding and structural thermodynamic results suggest that the non-active site mutants affect inhibitor binding by altering the geometry of the binding site cavity through the accumulation of mutations within the core of the protease molecule.     

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.     

Catalytic Water Co-Existing with a Product Peptide in the Active Site of HIV-1 Protease Revealed by X-Ray Structure Analysis

Background

It is known that HIV-1 protease is an important target for design of antiviral compounds in the treatment of Acquired Immuno Deficiency Syndrome (AIDS). In this context, understanding the catalytic mechanism of the enzyme is of crucial importance as transition state structure directs inhibitor design. Most mechanistic proposals invoke nucleophilic attack on the scissile peptide bond by a water molecule. But such a water molecule coexisting with any ligand in the active site has not been found so far in the crystal structures.

Principal Findings

We report here the first observation of the coexistence in the active site, of a water molecule WAT1, along with the carboxyl terminal product (Q product) peptide. The product peptide has been generated in situ through cleavage of the full-length substrate. The N-terminal product (P product) has diffused out and is replaced by a set of water molecules while the Q product is still held in the active site through hydrogen bonds. The position of WAT1, which hydrogen bonds to both the catalytic aspartates, is different from when there is no substrate bound in the active site. We propose WAT1 to be the position from where catalytic water attacks the scissile peptide bond. Comparison of structures of HIV-1 protease complexed with the same oligopeptide substrate, but at pH 2.0 and at pH 7.0 shows interesting changes in the conformation and hydrogen bonding interactions from the catalytic aspartates.

Conclusions/Significance

The structure is suggestive of the repositioning, during substrate binding, of the catalytic water for activation and subsequent nucleophilic attack. The structure could be a snap shot of the enzyme active site primed for the next round of catalysis. This structure further suggests that to achieve the goal of designing inhibitors mimicking the transition-state, the hydrogen-bonding pattern between WAT1 and the enzyme should be replicated.     

Design and synthesis of potent HIV-1 protease inhibitors incorporating hydroxyprolinamides as novel P2 ligands

A series of new HIV-1 protease inhibitors with the hydroxyethylamine core and different hydroxyprolinamide P2 ligands were designed and synthesized. Variation of substitutions at the P2 significantly affected the enzyme inhibitory potency of the inhibitors. Compounds 2a and 2d showed excellent enzyme inhibitory activity with IC₅₀ values in the nanomolar range. An active site binding model for inhibitors 2a and 2d was suggested based upon the computational-docking results of the ligand with HIV-1 protease. This model offers molecular insights regarding ligand-binding site interactions of the hydroxyprolinamide-derived novel P2-ligand.     

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.     

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.     

X-ray crystallographic structure of ABT-378 (lopinavir) bound to HIV-1 protease

The crystal structure of ABT-378 (lopinavir), bound to the active site of HIV-1 protease is described. A comparison with crystal structures of ritonavir, A-78791, and BILA-2450 shows some analogous features with previous reported compounds. A cyclic urea unit in the P(2) position of ABT-378 is novel and makes two bidentate hydrogen bonds with Asp 29 of HIV-1 protease. In addition, a previously unreported shift in the Gly 48 carbonyl position is observed. A discussion of the structural features responsible for its high potency against wild-type HIV protease is given along with an analysis of the effect of active site mutations on potency in in vitro assays.