Nine crystal structures determine the substrate envelope of the MDR HIV-1 protease

Decoding molecular mechanism underlying binding of drugs to HIV-1 protease with molecular dynamics simulations and MM-GBSA calculations

HIV-1 protease (PR) is thought to be efficient targets of anti-AIDS drug design. Molecular dynamics (MD) simulations and multiple post-processing analysis technologies were applied to decipher molecular mechanism underlying binding of three drugs Lopinavir (LPV), Nelfinavir (NFV) and Atazanavir (ATV) to the PR. Binding free energies calculated by molecular mechanics generalized Born surface area (MM-GBSA) suggest that compensation between binding enthalpy and entropy plays a vital role in binding of drugs to PR. Dynamics analyses show that binding of LPV, NFV and ATV highly affects structural flexibility, motion modes and dynamics behaviour of the PR, especially for two flaps. Computational alanine scanning and interaction network analysis verify that although three drugs have structural difference, they share similar binding modes to the PR and common interaction clusters with the PR. The current findings also confirm that residues located interaction clusters, such as Asp25/Asp25', Gly27/Gly27', Ala28/Ala28', Asp29, Ile47/Ile47', Gly49/Gly49', Ile50/Ile50', Val82/Val82' and Ile84/Ile84, can be used as efficient targets of clinically available inhibitors towards the PR.

Exploring the concerted mechanistic pathway for HIV-1 PR-substrate revealed by umbrella sampling simulation

HIV-1 protease (HIV-1 PR) is an essential enzyme for the replication process of its virus, and therefore considered an important target for the development of drugs against the acquired immunodeficiency syndrome (AIDS). Our previous study shows that the catalytic mechanism of subtype B/C-SA HIV-1 PR follows a one-step concerted acyclic hydrolysis reaction process using a two-layered ONIOM B3LYP/6-31++G(d,p) method. This present work is aimed at exploring the proposed mechanism of the proteolysis catalyzed by HIV-1 PR and to ensure our proposed mechanism is not an artefact of a single theoretical technique. Hence, we present umbrella sampling method that is suitable for calculating potential mean force (PMF) for non-covalent ligand/substrate-enzyme association/dissociation interactions which provide thermodynamic details for molecular recognition. The free activation energy results were computed in terms of PMF analysis within the hybrid QM(DFTB)/MM approach. The theoretical findings suggest that the proposed mechanism corresponds in principle with experimental data. Given our observations, we suggest that the QM/MM MD method can be used as a reliable computational technique to rationalize lead compounds against specific targets such as the HIV-1 protease.

Concerted hydrolysis mechanism of HIV-1 natural substrate against subtypes B and C-SA PR: insight through molecular dynamics and hybrid QM/MM studies

It is well known that understanding the catalytic mechanism of HIV-1 PR is the rationale on which its inhibitors were developed; therefore, a better understanding of the mechanism of natural substrate hydrolysis is important. Herein, the reaction mechanism of HIV-1 natural substrates with subtypes B and common mutant in South Africa (subtype C-SA) protease were studied through transition state modelling, using a general acid-general base (GA-GB) one-step concerted process. The activation free energies of enzyme-substrate complexes were compared based on their rate of hydrolysis using a two-layered ONIOM (B3LYP/6-31++G(d,p):AMBER) method. We expanded our computational model to obtain a better understanding of the mechanism of hydrolysis as well as how the enzyme recognises or chooses the cleavage site of the scissile bonds. Using this model, a potential substrate-based inhibitor could be developed with better potency. The calculated activation energies of natural substrates in our previous study correlated well with experimental data. A similar trend was observed for the Gag and Gag-Pol natural substrates in the present work for both enzyme complexes except for the PR-RT substrate. Natural bond orbital (NBO) analysis was also applied to determine the extent of charge transfer within the QM part of both enzymes considered and the PR-RT natural substrate. The result of this study shows that the method can be utilized as a dependable computational technique to rationalize lead compounds against specific targets.

A substrate selected by phage display exhibits enhanced side-chain hydrogen bonding to HIV-1 protease

Crystal structures of inactive variants of HIV-1 protease bound to peptides have revealed how the enzyme recognizes its endogenous substrates. The best of the known substrates is, however, a nonnatural substrate that was identified by directed evolution. The crystal structure of the complex between this substrate and the D25N variant of the protease is reported at a resolution of 1.1 Å. The structure has several unprecedented features, especially the formation of additional hydrogen bonds between the enzyme and the substrate. This work expands the understanding of molecular recognition by HIV-1 protease and informs the design of new substrates and inhibitors.

The Molecular Origin of Enthalpy/Entropy Compensation in Biomolecular Recognition

Biomolecular recognition can be stubborn; changes in the structures of associating molecules, or the environments in which they associate, often yield compensating changes in enthalpies and entropies of binding and no net change in affinities. This phenomenon—termed enthalpy/entropy (H/S) compensation—hinders efforts in biomolecular design, and its incidence—often a surprise to experimentalists—makes interactions between biomolecules difficult to predict. Although characterizing H/S compensation requires experimental care, it is unquestionably a real phenomenon that has, from an engineering perspective, useful physical origins. Studying H/S compensation can help illuminate the still-murky roles of water and dynamics in biomolecular recognition and self-assembly. This review summarizes known sources of H/ S compensation (real and perceived) and lays out a conceptual framework for understanding and dissecting—and, perhaps, avoiding or exploiting—this phenomenon in biophysical systems.

GRL-09510, a Unique P2-Crown-Tetrahydrofuranylurethane -Containing HIV-1 Protease Inhibitor, Maintains Its Favorable Antiviral Activity against Highly-Drug-Resistant HIV-1 Variants in vitro

We report that GRL-09510, a novel HIV-1 protease inhibitor (PI) containing a newly-generated P2-crown-tetrahydrofuranylurethane (Crwn-THF), a P2′-methoxybenzene, and a sulfonamide isostere, is highly active against laboratory and primary clinical HIV-1 isolates (EC₅₀: 0.0014–0.0028 μM) with minimal cytotoxicity (CC₅₀: 39.0 μM). Similarly, GRL-09510 efficiently blocked the replication of HIV-1_NL4-3 variants, which were capable of propagating at high-concentrations of atazanavir, lopinavir, and amprenavir (APV). GRL-09510 was also potent against multi-drug-resistant clinical HIV-1 variants and HIV-2_ROD. Under the selection condition, where HIV-1_NL4-3 rapidly acquired significant resistance to APV, an integrase inhibitor raltegravir, and a GRL-09510 congener (GRL-09610), no variants highly resistant against GRL-09510 emerged over long-term in vitro passage of the virus. Crystallographic analysis demonstrated that the Crwn-THF moiety of GRL-09510 forms strong hydrogen-bond-interactions with HIV-1 protease (PR) active-site amino acids and is bulkier with a larger contact surface, making greater van der Waals contacts with PR than the bis-THF moiety of darunavir. The present data demonstrate that GRL-09510 has favorable features for treating patients infected with wild-type and/or multi-drug-resistant HIV-1 variants, that the newly generated P2-Crwn-THF moiety confers highly desirable anti-HIV-1 potency. The use of the novel Crwn-THF moiety sheds lights in the design of novel PIs.

HIV Genome-Wide Protein Associations: a Review of 30 Years of Research

The HIV genome encodes a small number of viral proteins (i.e., 16), invariably establishing cooperative associations among HIV proteins and between HIV and host proteins, to invade host cells and hijack their internal machineries. As a known example, the HIV envelope glycoprotein GP120 is closely associated with GP41 for viral entry. From a genome-wide perspective, a hypothesis can be worked out to determine whether 16 HIV proteins could develop 120 possible pairwise associations either by physical interactions or by functional associations mediated via HIV or host molecules. Here, we present the first systematic review of experimental evidence on HIV genome-wide protein associations using a large body of publications accumulated over the past 3 decades. Of 120 possible pairwise associations between 16 HIV proteins, at least 34 physical interactions and 17 functional associations have been identified. To achieve efficient viral replication and infection, HIV protein associations play essential roles (e.g., cleavage, inhibition, and activation) during the HIV life cycle. In either a dispensable or an indispensable manner, each HIV protein collaborates with another viral protein to accomplish specific activities that precisely take place at the proper stages of the HIV life cycle. In addition, HIV genome-wide protein associations have an impact on anti-HIV inhibitors due to the extensive cross talk between drug-inhibited proteins and other HIV proteins. Overall, this study presents for the first time a comprehensive overview of HIV genome-wide protein associations, highlighting meticulous collaborations between all viral proteins during the HIV life cycle.

HIV Genome-Wide Protein Associations: a Review of 30 Years of Research

The HIV genome encodes a small number of viral proteins (i.e., 16), invariably establishing cooperative associations among HIV proteins and between HIV and host proteins, to invade host cells and hijack their internal machineries. As a known example, the HIV envelope glycoprotein GP120 is closely associated with GP41 for viral entry. From a genome-wide perspective, a hypothesis can be worked out to determine whether 16 HIV proteins could develop 120 possible pairwise associations either by physical interactions or by functional associations mediated via HIV or host molecules. Here, we present the first systematic review of experimental evidence on HIV genome-wide protein associations using a large body of publications accumulated over the past 3 decades. Of 120 possible pairwise associations between 16 HIV proteins, at least 34 physical interactions and 17 functional associations have been identified. To achieve efficient viral replication and infection, HIV protein associations play essential roles (e.g., cleavage, inhibition, and activation) during the HIV life cycle. In either a dispensable or an indispensable manner, each HIV protein collaborates with another viral protein to accomplish specific activities that precisely take place at the proper stages of the HIV life cycle. In addition, HIV genome-wide protein associations have an impact on anti-HIV inhibitors due to the extensive cross talk between drug-inhibited proteins and other HIV proteins. Overall, this study presents for the first time a comprehensive overview of HIV genome-wide protein associations, highlighting meticulous collaborations between all viral proteins during the HIV life cycle.

Anti-HIV Drug Discovery and Development: Current Innovations and Future Trends

The early effectiveness of combinatorial antiretroviral therapy (cART) in the treatment of HIV infection has been compromised to some extent by rapid development of multidrug-resistant HIV strains, poor bioavailability, and cumulative toxicities, and so there is a need for alternative strategies of antiretroviral drug discovery and additional therapeutic agents with novel action modes or targets. From this perspective, we first review current strategies of antiretroviral drug discovery and optimization, with the aid of selected examples from the recent literature. We highlight the development of phosphate ester-based prodrugs as a means to improve the aqueous solubility of HIV inhibitors, and the introduction of the substrate envelope hypothesis as a new approach for overcoming HIV drug resistance. Finally, we discuss future directions for research, including opportunities for exploitation of novel antiretroviral targets, and the strategy of activation of latent HIV reservoirs as a means to eradicate the virus.

The role of mutations at codons 32, 47, 54, and 90 in HIV-1 protease flap dynamics

Treatment of Human Immunodeficiency Virus remains challenging due to the emergence of drug resistant strains under the selective pressure produced by standard anti-retroviral therapy. To explore the structural mechanisms of drug resistance, we performed 40 ns molecular dynamics simulations on three multi-drug resistant HIV-1 protease clinical isolates from patients attending an infectious diseases clinic in Detroit, MI. We identify a novel structural role for the I47V, V32I, I54M and L90M major resistance mutations in flap opening and closure of MDR-PR isolates. Our studies suggest I47V is involved in flap opening and the interaction between I47V and V32I tethers the flaps to the active site. Also, I54M and L90M may be responsible for asymmetric movement of the protease flaps. These findings can be utilized to improve drug design strategies against MDR HIV-1 PR variants.

Ligand modifications to reduce the relative resistance of multi-drug resistant HIV-1 protease

Proper proteolytic processing of the HIV-1 Gag/Pol polyprotein is required for HIV infection and viral replication. This feature has made HIV-1 protease an attractive target for antiretroviral drug design for the treatment of HIV-1 infected patients. To examine the role of the P1 and P1'positions of the substrate in inhibitory efficacy of multi-drug resistant HIV-1 protease 769 (MDR 769), we performed a series of structure-function studies. Using the original CA/p2 cleavage site sequence, we generated heptapeptides containing one reduced peptide bond with an L to F and A to F double mutation at P1 and P1' (F-r-F), and an A to F at P1' (L-r-F) resulting in P1/P1' modified ligands. Here, we present an analysis of co-crystal structures of CA/p2 F-r-F, and CA/p2 L-r-F in complex with MDR 769. To examine conformational changes in the complex structure, molecular dynamic (MD) simulations were performed with MDR769-ligand complexes. MD trajectories show the isobutyl group of both the lopinavir analog and the CA/p2 L-r-F substrate cause a conformational change of in the active site of MDR 769. IC50 measurements suggest the non identical P1/P1' ligands (CA/p2 L-r-F and lopinavir analog) are more effective against MDR proteases as opposed to identical P1/P1'ligands. Our results suggest that a non identical P1/P1'composition may be more favorable for the inhibition of MDR 769 as they induce conformational changes in the active site of the enzyme resulting in disruption of the two-fold symmetry of the protease, thus, stabilizing the inhibitor in the active site.

Crystallographic study of multi-drug resistant HIV-1 protease lopinavir complex: mechanism of drug recognition and resistance

Lopinavir (LPV) is a second generation HIV-1 protease inhibitor. Drug resistance has rapidly emerged against LPV since its US FDA approval on September 15, 2000. Mutations at residues 32I, L33F, 46I, 47A, I54V, V82A, I84V, and L90M render the protease drug resistant against LPV. We report the crystal structure of a clinical isolate multi-drug resistant (MDR) 769 HIV-1 protease (resistant mutations at residues 10, 36, 46, 54, 62, 63, 71, 82, 84, and 90) complexed with LPV and the in vitro enzymatic IC50 of LPV against MDR 769. The structural and functional studies demonstrate significant drug resistance of MDR 769 against LPV, arising from reduced interactions between LPV and the protease target.

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.

Correlating conformational shift induction with altered inhibitor potency in a multidrug resistant HIV-1 protease variant

Inhibitor-induced conformational ensemble shifts in a multidrug resistant HIV-1 protease variant, MDR769, are characterized by site-directed spin labeling double electron-electron resonance spectroscopy. For MDR769 compared to the native enzyme, changes in inhibitor IC(50) values are related to a parameter defined as |ΔC|, which is the relative change in the inhibitor-induced shift to the closed state. Specifically, a linear correlation is found between |ΔC| and the magnitude of the change in IC(50), provided that inhibitor binding is not too weak. Moreover, inhibitors that exhibit MDR769 resistance no longer induce a strong shift to a closed conformational ensemble as seen previously in the native enzyme.

From base pair to bedside: molecular simulation and the translation of genomics to personalized medicine

Despite the promises made that genomic sequencing would transform therapy by introducing a new era of personalized medicine, relatively few tangible breakthroughs have been made. This has led to the recognition that complex interactions at multiple spatial, temporal, and organizational levels may often combine to produce disease. Understanding this complexity requires that existing and future models are used and interpreted within a framework that incorporates knowledge derived from investigations at multiple levels of biological function. It also requires a computational infrastructure capable of dealing with the vast quantities of data generated by genomic approaches. In this review, we discuss the use of molecular modeling to generate quantitative and qualitative insights at the smallest scales of the systems biology hierarchy, how it can play an important role in the development of a systems understanding of disease and in the application of such knowledge to help discover new therapies and target existing ones on a personal level.

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.

Context surrounding processing sites is crucial in determining cleavage rate of a subset of processing sites in HIV-1 Gag and Gag-Pro-Pol polyprotein precursors by viral protease

Processing of the human immunodeficiency virus type 1 (HIV-1) Gag and Gag-Pro-Pol polyproteins by the HIV-1 protease (PR) is essential for the production of infectious particles. However, the determinants governing the rates of processing of these substrates are not clearly understood. We studied the effect of substrate context on processing by utilizing a novel protease assay in which a substrate containing HIV-1 matrix (MA) and the N-terminal domain of capsid (CA) is labeled with a FlAsH (fluorescein arsenical hairpin) reagent. When the seven cleavage sites within the Gag and Gag-Pro-Pol polyproteins were placed at the MA/CA site, the rates of cleavage changed dramatically compared with that of the cognate sites in the natural context reported previously. The rate of processing was affected the most for three sites: CA/spacer peptide 1 (SP1) (≈10-fold increase), SP1/nucleocapsid (NC) (≈10-30-fold decrease), and SP2/p6 (≈30-fold decrease). One of two multidrug-resistant (MDR) PR variants altered the pattern of processing rates significantly. Cleavage sites within the Pro-Pol region were cleaved in a context-independent manner, suggesting for these sites that the sequence itself was the determinant of rate. In addition, a chimera consisting of SP1/NC P4-P1 and MA/CA P1'-P4' residues (ATIM↓PIVQ) abolished processing by wild type and MDR proteases, and the reciprocal chimera consisting of MA/CA P4-P1 and SP1/NC P1'-4' (SQNY↓IQKG) was cleaved only by one of the MDR proteases. These results suggest that complex substrate interactions both beyond the active site of the enzyme and across the scissile bond contribute to defining the rate of processing by the HIV-1 PR.

The higher barrier of darunavir and tipranavir resistance for HIV-1 protease

Darunavir and tipranavir are two inhibitors that are active against multi-drug resistant (MDR) HIV-1 protease variants. In this study, the invitro inhibitory efficacy was tested against a MDR HIV-1 protease variant, MDR 769 82T, containing the drug resistance mutations of 46L/54V/82T/84V/90M. Crystallographic and enzymatic studies were performed to examine the mechanism of resistance and the relative maintenance of potency. The key findings are as follows: (i) The MDR protease exhibits decreased susceptibility to all nine HIV-1 protease inhibitors approved by the US Food and Drug Administration (FDA), among which darunavir and tipranavir are the most potent; (ii) the threonine 82 mutation on the protease greatly enhances drug resistance by altering the hydrophobicity of the binding pocket; (iii) darunavir or tipranavir binding facilitates closure of the wide-open flaps of the MDR protease; and (iv) the remaining potency of tipranavir may be preserved by stabilizing the flaps in the inhibitor-protease complex while darunavir maintains its potency by preserving protein main chain hydrogen bonds with the flexible P2 group. These results could provide new insights into drug design strategies to overcome multi-drug resistance of HIV-1 protease variants.