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

Design and Development of Highly Potent HIV-1 Protease Inhibitors with a Crown-Like Oxotricyclic Core as the P2-Ligand To Combat Multidrug-Resistant HIV Variants

Design, synthesis, and evaluation of a new class of exceptionally potent HIV-1 protease inhibitors are reported. Inhibitor 5 displayed superior antiviral activity and drug-resistance profiles. In fact, this inhibitor showed several orders of magnitude improved antiviral activity over the FDA approved drug darunavir. This inhibitor incorporates an unprecedented 6-5-5 ring-fused crown-like tetrahydropyranofuran as the P2 ligand and an aminobenzothiazole as the P2' ligand with the (R)-hydroxyethylsulfonamide isostere. The crown-like P2 ligand for this inhibitor has been synthesized efficiently in an optically active form using a chiral Diels-Alder catalyst providing a key intermediate in high enantiomeric purity. Two high resolution X-ray structures of inhibitor-bound HIV-1 protease revealed extensive interactions with the backbone atoms of HIV-1 protease and provided molecular insight into the binding properties of these new inhibitors.

Application of computational methods for anticancer drug discovery, design, and optimization

Developing a novel drug is a complex, risky, expensive and time-consuming venture. It is estimated that the conventional drug discovery process ending with a new medicine ready for the market can take up to 15 years and more than a billion USD. Fortunately, this scenario has recently changed with the arrival of new approaches. Many novel technologies and methodologies have been developed to increase the efficiency of the drug discovery process, and computational methodologies have become a crucial component of many drug discovery programs. From hit identification to lead optimization, techniques such as ligand- or structure-based virtual screening are widely used in many discovery efforts. It is the case for designing potential anticancer drugs and drug candidates, where these computational approaches have had a major impact over the years and have provided fruitful insights into the field of cancer. In this paper, we review the concept of rational design presenting some of the most representative examples of molecules identified by means of it. Key principles are illustrated through case studies including specifically successful achievements in the field of anticancer drug design to demonstrate that research advances, with the aid of in silico drug design, have the potential to create novel anticancer drugs.

Application of computational methods for anticancer drug discovery, design, and optimization

Developing a novel drug is a complex, risky, expensive and time-consuming venture. It is estimated that the conventional drug discovery process ending with a new medicine ready for the market can take up to 15 years and more than a billion USD. Fortunately, this scenario has recently changed with the arrival of new approaches. Many novel technologies and methodologies have been developed to increase the efficiency of the drug discovery process, and computational methodologies have become a crucial component of many drug discovery programs. From hit identification to lead optimization, techniques such as ligand- or structure-based virtual screening are widely used in many discovery efforts. It is the case for designing potential anticancer drugs and drug candidates, where these computational approaches have had a major impact over the years and have provided fruitful insights into the field of cancer. In this paper, we review the concept of rational design presenting some of the most representative examples of molecules identified by means of it. Key principles are illustrated through case studies including specifically successful achievements in the field of anticancer drug design to demonstrate that research advances, with the aid of in silico drug design, have the potential to create novel anticancer drugs.

Investigational protease inhibitors as antiretroviral therapies

INTRODUCTION

Highly Active Antiretroviral Therapy (HAART) has tremendously improved the life expectancy of the HIV-infected population over the past three decades. Protease inhibitors have been one of the major classes of drugs in HAART regimens that are effective in treating HIV. However, the emergence of resistance and cross-resistance against protease inhibitors encourages researchers to develop new PIs with broad-spectrum activity, as well as novel means of enhancing the efficacy of existing PIs.

AREAS COVERED

In this article we discuss recent advances in HIV protease inhibitor (PI) development, focusing on both investigational and experimental agents. We also include a section on pharmacokinetic booster drugs for improved bioavailability of protease inhibitors. Further, we discuss novel drug delivery systems using a variety of nanocarriers for the delivery of PIs across the blood-brain barrier to treat the HIV in the brain.

EXPERT OPINION

We discuss our opinion on the promises and challenges on the development of novel investigational and experimental PIs that are less toxic and more effective in combating drug-resistance. Further, we discuss the future of novel nanocarriers that have been developed to deliver PIs to the brain cells. Although these are promising findings, many challenges need to be overcome prior to making them a viable option.

Probing Lipophilic Adamantyl Group as the P1-Ligand for HIV-1 Protease Inhibitors: Design, Synthesis, Protein X-ray Structural Studies, and Biological Evaluation

A series of potent HIV-1 protease inhibitors with a lipophilic adamantyl P1 ligand have been designed, synthesized, and evaluated. We have developed an enantioselective synthesis of adamantane-derived hydroxyethylamine isosteres utilizing Sharpless asymmetric epoxidation as the key step. Various inhibitors incorporating P1-adamantylmethyl in combination with P2 ligands such as 3-(R)-THF, 3-(S)-THF, bis-THF, and THF-THP were examined. The S1' pocket was also probed with phenyl and phenylmethyl ligands. Inhibitor 15d, with an isobutyl P1' ligand and a bis-THF P2 ligand, proved to be the most potent of the series. The cLogP value of inhibitor 15d is improved compared to inhibitor 2 with a phenylmethyl P1-ligand. X-ray structural studies of 15d, 15h, and 15i with HIV-1 protease complexes revealed molecular insight into the inhibitor-protein interaction.

Discovery of HIV Type 1 Aspartic Protease Hit Compounds through Combined Computational Approaches

A combination of computational techniques and inhibition assay experiments was employed to identify hit compounds from commercial libraries with enhanced inhibitory potency against HIV type 1 aspartic protease (HIV PR). Extensive virtual screening with the aid of reliable pharmacophore models yielded five candidate protease inhibitors. Subsequent molecular dynamics and molecular mechanics Poisson-Boltzmann surface area free-energy calculations for the five ligand-HIV PR complexes suggested a high stability of the systems through hydrogen-bond interactions between the ligands and the protease's flaps (Ile50/50'), as well as interactions with residues of the active site (Asp25/25'/29/29'/30/30'). Binding-energy calculations for the three most promising compounds yielded values between -5 and -10 kcal mol(-1) and suggested that van der Waals interactions contribute most favorably to the total energy. The predicted binding-energy values were verified by in vitro inhibition assays, which showed promising results in the high nanomolar range. These results provide structural considerations that may guide further hit-to-lead optimization toward improved anti-HIV drugs.

Recent Progress in the Development of HIV-1 Protease Inhibitors for the Treatment of HIV/AIDS

HIV-1 protease inhibitors continue to play an important role in the treatment of HIV/AIDS, transforming this deadly ailment into a more manageable chronic infection. Over the years, intensive research has led to a variety of approved protease inhibitors for the treatment of HIV/AIDS. In this review, we outline current drug design and medicinal chemistry efforts toward the development of next-generation protease inhibitors beyond the currently approved drugs.

Virological failure in patients with HIV-1 subtype C receiving antiretroviral therapy: an analysis of a prospective national cohort in Sweden

BACKGROUND

People with HIV-1 in low-income and middle-income countries increasingly need second-line regimens with boosted protease inhibitors. However, data are scarce for treatment response in patients with HIV-1 subtype C (HIV-1C), which is predominant in these regions. We aimed to examine factors associated with virological failure in patients in a standardised national health-care setting.

METHODS

We analysed data for participants in InfCare HIV, a prospective national cohort that includes more than 99% of people with HIV in Sweden. We extracted data for the cohort from the InfCare HIV database on Jan 14, 2015. Baseline was initiation of antiretroviral therapy. We used logistic regression to assess factors associated with primary virological failure (failure to suppress HIV-1 within 9 months) in patients with HIV-1B and HIV-1C and calculated odds ratios (OR) for failure. We also used Cox regression models to calculate hazard ratios (HR) for time-to-secondary virological failure (detectable viral load after initial virological suppression). We did homology-based molecular modelling to assess docking.

FINDINGS

We included 1077 patients with HIV-1B and 596 with HIV-1C. In multivariate regression analysis, pre-therapy higher viral load (OR 1·82, 95% CI 1·49-2·21; p<0·0001), subtype C infection (1·75, 1·06-2·88; p=0·028), and boosted protease inhibitor-based regimens (1·55, 1·45-2·11; p=0·004) were associated with increased risk of primary virological failure. Individuals with HIV-1C who were given therapy with boosted protease inhibitors had earlier time-to-secondary virological failure than did those with HIV-1B given similar regimens (adjusted HR 1·92, 95% CI 1·30-2·83; p=0·002). Molecular modelling suggested lower affinity for protease inhibitors to HIV-1C protease than to HIV-1B.

INTERPRETATION

Our findings suggest an increased risk of virological failure in patients with HIV-1C, especially in those on boosted protease inhibitor-based regimens. Future studies should further dissect the biochemical and viral mechanisms of resistance to protease inhibitors in patients with non-B subtypes of HIV-1, including clinical studies to assess the efficacy of boosted protease inhibitor-based regimens in low-income and middle income countries.

FUNDING

Karolinska Institutet Research Foundation, Swedish Research Council, Stockholm County Council, Swedish Physicians against AIDS, US National Institutes of Health, University of Missouri.

Long-Range Electrostatics-Induced Two-Proton Transfer Captured by Neutron Crystallography in an Enzyme Catalytic Site

Neutron crystallography was used to directly locate two protons before and after a pH-induced two-proton transfer between catalytic aspartic acid residues and the hydroxy group of the bound clinical drug darunavir, located in the catalytic site of enzyme HIV-1 protease. The two-proton transfer is triggered by electrostatic effects arising from protonation state changes of surface residues far from the active site. The mechanism and pH effect are supported by quantum mechanics/molecular mechanics (QM/MM) calculations. The low-pH proton configuration in the catalytic site is deemed critical for the catalytic action of this enzyme and may apply more generally to other aspartic proteases. Neutrons therefore represent a superb probe to obtain structural details for proton transfer reactions in biological systems at a truly atomic level.

Design, synthesis, biological evaluation and X-ray structural studies of HIV-1 protease inhibitors containing substituted fused-tetrahydropyranyl tetrahydrofuran as P2-ligands

Design, synthesis, biological and X-ray crystallographic studies of a series of potent HIV-1 protease inhibitors are described. Various polar functionalities have been incorporated on the tetrahydropyranyl-tetrahydrofuran-derived P2 ligand to interact with the backbone atoms in the S2-subsite. The majority of the inhibitors showed very potent enzyme inhibitory and antiviral activity. Two high-resolution X-ray structures of 30b- and 30j-bound HIV-1 protease provide insight into ligand-binding site interactions. In particular, the polar functionalities on the P2-ligand appear to form unique hydrogen bonds with Gly48 amide NH and amide carbonyl groups in the flap region.

Molecular mechanisms and design principles for promiscuous inhibitors to avoid drug resistance: lessons learned from HIV-1 protease inhibition

Molecular recognition is central to biology and ranges from highly selective to broadly promiscuous. The ability to modulate specificity at will is particularly important for drug development, and discovery of mechanisms contributing to binding specificity is crucial for our basic understanding of biology and for applications in health care. In this study, we used computational molecular design to create a large dataset of diverse small molecules with a range of binding specificities. We then performed structural, energetic, and statistical analysis on the dataset to study molecular mechanisms of achieving specificity goals. The work was done in the context of HIV‐1 protease inhibition and the molecular designs targeted a panel of wild‐type and drug‐resistant mutant HIV‐1 protease structures. The analysis focused on mechanisms for promiscuous binding to bind robustly even to resistance mutants. Broadly binding inhibitors tended to be smaller in size, more flexible in chemical structure, and more hydrophobic in nature compared to highly selective ones. Furthermore, structural and energetic analyses illustrated mechanisms by which flexible inhibitors achieved binding; we found ligand conformational adaptation near mutation sites and structural plasticity in targets through torsional flips of asymmetric functional groups to form alternative, compensatory packing interactions or hydrogen bonds. As no inhibitor bound to all variants, we designed small cocktails of inhibitors to do so and discovered that they often jointly covered the target set through mechanistic complementarity. Furthermore, using structural plasticity observed in experiments, and potentially in simulations, is suggested to be a viable means of designing adaptive inhibitors that are promiscuous binders. Proteins 2015; 83:351–372. © 2014 Wiley Periodicals, Inc.

Design of HIV-1 Protease Inhibitors with Amino-bis-tetrahydrofuran Derivatives as P2-Ligands to Enhance Backbone-Binding Interactions: Synthesis, Biological Evaluation, and Protein-Ligand X-ray Studies

Structure-based design, synthesis, and biological evaluation of a series of very potent HIV-1 protease inhibitors are described. In an effort to improve backbone ligand-binding site interactions, we have incorporated basic-amines at the C4 position of the bis-tetrahydrofuran (bis-THF) ring. We speculated that these substituents would make hydrogen bonding interactions in the flap region of HIV-1 protease. Synthesis of these inhibitors was performed diastereoselectively. A number of inhibitors displayed very potent enzyme inhibitory and antiviral activity. Inhibitors 25f, 25i, and 25j were evaluated against a number of highly-PI-resistant HIV-1 strains, and they exhibited improved antiviral activity over darunavir. Two high resolution X-ray structures of 25f- and 25g-bound HIV-1 protease revealed unique hydrogen bonding interactions with the backbone carbonyl group of Gly48 as well as with the backbone NH of Gly48 in the flap region of the enzyme active site. These ligand-binding site interactions are possibly responsible for their potent activity.

Conformational variation of an extreme drug resistant mutant of HIV protease

Molecular mechanisms leading to high level drug resistance have been analyzed for the clinical variant of HIV-1 protease bearing 20 mutations (PR20); which has several orders of magnitude worse affinity for tested drugs. Two crystal structures of ligand-free PR20 with the D25N mutation of the catalytic aspartate (PR20D25N) revealed three dimers with different flap conformations. The diverse conformations of PR20D25N included a dimer with one flap in a unique "tucked" conformation; directed into the active site. Analysis of molecular dynamics (MD) simulations of the ligand-free PR20 and wild-type enzymes showed that the mutations in PR20 alter the correlated interactions between two monomers in the dimer. The two flaps tend to fluctuate more independently in PR20 than in the wild type enzyme. Combining the results of structural analysis by X-ray crystallography and MD simulations; unusual flap conformations and weakly correlated inter-subunit motions may contribute to the high level resistance of PR20.

Fluorogenic Assay for Inhibitors of HIV-1 Protease with Sub-picomolar Affinity

A fluorogenic substrate for HIV-1 protease was designed and used as the basis for a hypersensitive assay. The substrate exhibits a k_cat of 7.4 s⁻¹, K_M of 15 μM, and an increase in fluorescence intensity of 104-fold upon cleavage, thus providing sensitivity that is unmatched in a continuous assay of HIV-1 protease. These properties enabled the enzyme concentration in an activity assay to be reduced to 25 pM, which is close to the K_d value of the protease dimer. By fitting inhibition data to Morrison’s equation, K_i values of amprenavir, darunavir, and tipranavir were determined to be 135, 10, and 82 pM, respectively. This assay, which is capable of measuring K_i values as low as 0.25 pM, is well-suited for characterizing the next generation of HIV-1 protease inhibitors.

Structure-based design of potent HIV-1 protease inhibitors with modified P1-biphenyl ligands: synthesis, biological evaluation, and enzyme-inhibitor X-ray structural studies

We report the design, synthesis, X-ray structural studies, and biological evaluation of a novel series of HIV-1 protease inhibitors. We designed a variety of functionalized biphenyl derivatives to make enhanced van der Waals interactions in the S1 subsite of HIV-1 protease. These biphenyl derivatives were conveniently synthesized using a Suzuki-Miyaura cross-coupling reaction as the key step. We examined the potential of these functionalized biphenyl-derived P1 ligands in combination with 3-(S)-tetrahydrofuranyl urethane and bis-tetrahydrofuranyl urethane as the P2 ligands. Inhibitor 21e, with a 2-methoxy-1,1'-biphenyl derivative as P1 ligand and bis-THF as the P2 ligand, displayed the most potent enzyme inhibitory and antiviral activity. This inhibitor also exhibited potent activity against a panel of multidrug-resistant HIV-1 variants. A high resolution X-ray crystal structure of related Boc-derivative 17a-bound HIV-1 protease provided important molecular insight into the ligand-binding site interactions of the biphenyl core in the S1 subsite of HIV-1 protease.

Structure-based design, synthesis, X-ray studies, and biological evaluation of novel HIV-1 protease inhibitors containing isophthalamide-derived P2-ligands

We describe the design, synthesis and biological evaluation of a series of novel HIV-1 protease inhibitors bearing isophthalamide derivatives as the P2-P3 ligands. We have investigated a range of acyclic and heterocyclic amides as the extended P2-P3 ligands. These inhibitors displayed good to excellent HIV-1 protease inhibitory activity. Also, a number of inhibitors showed very good antiviral activity in MT cells. Compound 5n has shown an enzyme Ki of 0.17 nM and antiviral IC50 of 14 nM. An X-ray crystal structure of inhibitor 5o-bound to HIV-1 protease was determined at 1.11Å resolution. This structure revealed important molecular insight into the inhibitor-HIV-1 protease interactions in the active site.

Disubstituted Bis-THF Moieties as New P2 Ligands in Nonpeptidal HIV-1 Protease Inhibitors (II)

A series of darunavir analogues featuring a substituted bis-THF ring as P2 ligand have been synthesized and evaluated. Very high affinity protease inhibitors (PIs) with an interesting activity on wild-type HIV and a panel of multi-PI resistant HIV-1 mutants containing clinically observed, primary mutations were identified using a cell-based assay. Crystal structure analysis was conducted on a number of PI analogues in complex with HIV-1 protease.

Organic carbamates in drug design and medicinal chemistry

The carbamate group is a key structural motif in many approved drugs and prodrugs. There is an increasing use of carbamates in medicinal chemistry and many derivatives are specifically designed to make drug-target interactions through their carbamate moiety. In this Perspective, we present properties and stabilities of carbamates, reagents and chemical methodologies for the synthesis of carbamates, and recent applications of carbamates in drug design and medicinal chemistry.

HIV protease inhibitors: a review of molecular selectivity and toxicity

Highly active antiretroviral therapy (HAART) is recognized as the most effective treatment method for AIDS, and protease inhibitors play a very important role in HAART. However, poor bioavailability and unbearable toxicity are their common disadvantages. Thus, the development of safer and potentially promising protease inhibitors is eagerly needed. In this review, we introduced the chemical characteristics and associated side effects of HIV protease inhibitors, as well as the possible off-target mechanisms causing the side effects. From the chemical structures of HIV protease inhibitors and their possible off-target molecules, we could obtain hints for optimizing the molecular selectivity of the inhibitors, to provide help in the design of new compounds with enhanced bioavailability and reduced side effects.

Exploring the drug resistance of V32I and M46L mutant HIV-1 protease to inhibitor TMC114: flap dynamics and binding mechanism

Inhibitors of HIV-1 protease (HIV-1-pr) generally only bind to the active site of the protease. However, for some mutants such as V32I and M46L the TMC114 can bind not only to the active cavity but also to the groove of the flexible flaps. Although the second binding site suggests the higher efficiency of the drug against HIV-1-pr, the drug resistance in HIV-1-pr due to mutations cannot be ignored, which prompts us to investigate the molecular mechanisms of drug resistance and behavior of double bound TMC114 (2T) to HIV-1-pr. The conformational dynamics of HIV-1-pr and the binding of TMC114 to the WT, V32I and M46L mutants were investigated with all-atom molecular dynamic (MD) simulation. The 20 ns MD simulation shows many fascinating effects of the inhibitor binding to the WT and mutant proteases. MM-PBSA calculations explain the binding free energies unfavorable for the M46L and V32I mutants as compared to the WT. For the single binding (1T) the less binding affinity can be attributed to the entropic loss for both V32I-1T and M46L-1T. Although the second binding of TMC114 with flap does increase binding energy for the mutants (V32I-2T and M46L-2T), the considerable entropy loss results in the lower binding Gibbs free energies. Thus, binding of TMC114 in the flap region does not help much in the total gain in binding affinity of the system, which was verified from this study and thereby validating experiments.

Design of gem-difluoro-bis-tetrahydrofuran as P2 ligand for HIV-1 protease inhibitors to improve brain penetration: synthesis, X-ray studies, and biological evaluation

The structure-based design, synthesis, biological evaluation, and X-ray structural studies of fluorine-containing HIV-1 protease inhibitors are described. The synthesis of both enantiomers of the gem-difluoro-bis-THF ligands was carried out in a stereoselective manner using a Reformatskii-Claisen reaction as the key step. Optically active ligands were converted into protease inhibitors. Two of these inhibitors, (3R,3aS,6aS)-4,4-difluorohexahydrofuro[2,3-b]furan-3-yl(2S,3R)-3-hydroxy-4-((N-isobutyl-4-methoxyphenyl)sulfonamido)-1-phenylbutan-2-yl) carbamate (3) and (3R,3aS,6aS)-4,4-difluorohexahydrofuro[2,3-b]furan-3-yl(2S,3R)-3-hydroxy-4-((N-isobutyl-4-aminophenyl)sulfonamido)phenylbutan-2-yl) carbamate (4), exhibited HIV-1 protease inhibitory Ki values in the picomolar range. Both 3 and 4 showed very potent antiviral activity, with respective EC50 values of 0.8 and 3.1 nM against the laboratory strain HIV-1LAI . The two inhibitors exhibited better lipophilicity profiles than darunavir, and also showed much improved blood-brain barrier permeability in an in vitro model. A high-resolution X-ray structure of inhibitor 4 in complex with HIV-1 protease was determined, revealing that the fluorinated ligand makes extensive interactions with the S2 subsite of HIV-1 protease, including hydrogen bonding interactions with the protease backbone atoms. Moreover, both fluorine atoms on the bis-THF ligand formed strong interactions with the flap Gly 48 carbonyl oxygen atom.

Design and synthesis of potent macrocyclic HIV-1 protease inhibitors involving P1-P2 ligands

A series of potent macrocyclic HIV-1 protease inhibitors have been designed and synthesized. The compounds incorporated 16- to 19-membered macrocyclic rings between a nelfinavir-like P2 ligand and a tyrosine side chain containing a hydroxyethylamine sulfonamide isostere. All cyclic inhibitors are more potent than their corresponding acyclic counterparts. Saturated derivatives showed slight reduction of potency compared to the respective unsaturated derivatives. Compound containing a 16-membered ring as the P1-P2 ligand showed the most potent enzyme inhibitory and antiviral activity.

Placental transfer of darunavir in an ex vivo human cotyledon perfusion model

Placental transfer of the HIV protease inhibitor darunavir was investigated in 5 term human cotyledons perfused with darunavir (1,000 ng/ml) in the maternal to fetal direction. The mean (± the standard deviation [SD]) fetal transfer rate (FTR) (fetal/maternal concentration at steady state from 30 to 90 min) was 15.0%±2.1%, and the mean (±SD) clearance index (darunavir FTR/antipyrine FTR) was 40.3%±5.8%. This shows that darunavir crosses the placenta at a relatively low rate, resulting in fetal exposure.

Systematic molecular dynamics, MM-PBSA, and ab initio approaches to the saquinavir resistance mechanism in HIV-1 PR due to 11 double and multiple mutations

Mutations in the human immunodeficiency virus (HIV) enable virus replication even when appropriate antiretroviral therapy is followed, thus leading to the emergence of drug resistance. In a previous work, we systematically examined seven single mutations that are associated with saquinavir (SQV) resistance in HIV-1 protease (Tzoupis, H.; Leonis, G.; Mavromoustakos, T.; Papadopoulos, M. G. J. Chem. Theory Comput. 2013, 9, 1754-1764). Herein, we extend our analysis, which includes seven double (G48V-V82A, L10I-G48V, G48V-L90M, I84V-L90M, L10I-V82A, L10I-L63P, A71V-G73S) and four multiple (L10I-L63P-A71V, L10I-G48V-V82A, G73S-I84V-L90M, L10I-L63P-A71V-G73S-I84V-L90M) SQV-HIV-1 PR mutant complexes, in an attempt to generalize our findings and formulate the main elements of the SQV resistance mechanism in the protease. On the basis of molecular dynamics (MD), molecular mechanics Poisson-Boltzmann surface area (MM-PBSA), and ab initio computational approaches, we identified specific features that constitute the HIV-1 PR mechanism of resistance at the molecular level: the low flexibility of SQV in the binding cavity and the preservation of hydrogen bonding (HB) and van der Waals interactions between SQV and several active-site (Gly27/27', Asp29/29'/30/30', especially Asp25/25') and flap (Ile50/50', Gly48/48') residues of the protease contribute significantly to efficient binding. The total enthalpy loss in all mutants is mostly due to the loss in enthalpy of the active-site region. Furthermore, it was observed that mutation accumulation may induce stabilization to SQV and to the flaps through enhanced HB interactions that lead to improved inhibition (e.g., accumulation of mutations in complexes containing L10I, G48V, L63P, I84V, or L90M single mutations). It was also concluded that permanent flap closure is obtained independently of mutations and SQV binding is mostly driven by van der Waals, nonpolar, and exchange-energy contributions. Importantly, it was indicated that the optimal positioning of SQV and the structure of the binding cavity are tightly coupled, since small changes in geometry may affect the binding energy greatly. The results of our theoretical approaches are in agreement with experimental evidence and provide a reliable description of SQV resistance in HIV-1 PR.

Prediction of HIV drug resistance from genotype with encoded three-dimensional protein structure

BACKGROUND

Drug resistance has become a severe challenge for treatment of HIV infections. Mutations accumulate in the HIV genome and make certain drugs ineffective. Prediction of resistance from genotype data is a valuable guide in choice of drugs for effective therapy.

RESULTS

In order to improve the computational prediction of resistance from genotype data we have developed a unified encoding of the protein sequence and three-dimensional protein structure of the drug target for classification and regression analysis. The method was tested on genotype-resistance data for mutants of HIV protease and reverse transcriptase. Our graph based sequence-structure approach gives high accuracy with a new sparse dictionary classification method, as well as support vector machine and artificial neural networks classifiers. Cross-validated regression analysis with the sparse dictionary gave excellent correlation between predicted and observed resistance.

CONCLUSION

The approach of encoding the protein structure and sequence as a 210-dimensional vector, based on Delaunay triangulation, has promise as an accurate method for predicting resistance from sequence for drugs inhibiting HIV protease and reverse transcriptase.

Structures of darunavir-resistant HIV-1 protease mutant reveal atypical binding of darunavir to wide open flaps

The molecular basis for high resistance to clinical inhibitors of HIV-1 protease (PR) was examined for the variant designated PR_P51 that was selected for resistance to darunavir (DRV). High resolution crystal structures of PR_P51 with the active site D25N mutation revealed a ligand-free form and an inhibitor-bound form showing a unique binding site and orientation for DRV. This inactivating mutation is known to increase the dimer dissociation constant and decrease DRV affinity of PR. The PR_P51-D25N dimers were in the open conformation with widely separated flaps, as reported for other highly resistant variants. PR_P51-D25N dimer bound two DRV molecules and showed larger separation of 8.7 Å between the closest atoms of the two flaps compared with 4.4 Å for the ligand-free structure of this mutant. The ligand-free structure, however, lacked van der Waals contacts between Ile50 and Pro81′ from the other subunit in the dimer, unlike the majority of PR structures. DRV is bound inside the active site cavity; however, the inhibitor is oriented almost perpendicular to its typical position and exhibits only 2 direct hydrogen bond and two water-mediated interactions with atoms of PR_P51-D25N compared with 11 hydrogen bond interactions seen for DRV bound in the typical position in wild-type enzyme. The atypical location of DRV may provide opportunities for design of novel inhibitors targeting the open conformation of PR drug-resistant mutants.

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

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

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

Understanding HIV-1 protease autoprocessing for novel therapeutic development

In the infected cell, HIV-1 protease (PR) is initially synthesized as part of the GagPol polyprotein. PR autoprocessing is a virus-specific process by which the PR domain embedded in the precursor catalyzes proteolytic reactions responsible for liberation of free mature PRs, which then recognize and cleave at least ten different peptide sequences in the Gag and GagPol polyproteins. Despite extensive structure and function studies of the mature PRs as well as the successful development of ten US FDA-approved catalytic-site inhibitors, the precursor autoprocessing mechanism remains an intriguing yet-to-be-solved puzzle. This article discusses current understanding of the autoprocessing mechanism, in an effort to prompt the development of novel anti-HIV drugs that selectively target precursor autoprocessing.

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

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

Highly potent HIV-1 protease inhibitors with novel tricyclic P2 ligands: design, synthesis, and protein-ligand X-ray studies

The design, synthesis, and biological evaluation of a series of HIV-1 protease inhibitors incorporating stereochemically defined fused tricyclic P2 ligands are described. Various substituent effects were investigated to maximize the ligand-binding site interactions in the protease active site. Inhibitors 16a and 16f showed excellent enzyme inhibitory and antiviral activity, although the incorporation of sulfone functionality resulted in a decrease in potency. Both inhibitors 16a and 16f maintained activity against a panel of multidrug resistant HIV-1 variants. A high-resolution X-ray crystal structure of 16a-bound HIV-1 protease revealed important molecular insights into the ligand-binding site interactions, which may account for the inhibitor's potent antiviral activity and excellent resistance profiles.

Novel P2 tris-tetrahydrofuran group in antiviral compound 1 (GRL-0519) fills the S2 binding pocket of selected mutants of HIV-1 protease

GRL-0519 (1) is a potent antiviral inhibitor of HIV-1 protease (PR) possessing tris-tetrahydrofuran (tris-THF) at P2. The high resolution X-ray crystal structures of inhibitor 1 in complexes with single substitution mutants PR(R8Q), PR(D30N), PR(I50V), PR(I54M), and PR(V82A) were analyzed in relation to kinetic data. The smaller valine side chain in PR(I50V) eliminated hydrophobic interactions with inhibitor and the other subunit consistent with 60-fold worse inhibition. Asn30 in PR(D30N) showed altered interactions with neighboring residues and 18-fold worse inhibition. Mutations V82A and I54M showed compensating structural changes consistent with 6-7-fold lower inhibition. Gln8 in PR(R8Q) replaced the ionic interactions of wild type Arg8 with hydrogen bond interactions without changing the inhibition significantly. The carbonyl oxygen of Gly48 showed two alternative conformations in all structures likely due to the snug fit of the large tris-THF group in the S2 subsite in agreement with high antiviral efficacy of 1 on resistant virus.

Binding of single walled carbon nanotube to WT and mutant HIV-1 proteases: analysis of flap dynamics and binding mechanism

Most of the currently treated HIV-1 protease (HIV-PR) inhibitors have been prone to suffer from the mutations associated drug resistance. Therefore, it is necessary to search for potent alternatives against the drug resistance. In the current study we have tested the single-walled carbon nanotube (SWCNT) as an inhibitor in wild type (WT) as well as in three primary mutants (I50V(PR), V82A(PR) and I84V(PR)) of the HIV-1-PR through docking the SWCNT in the active site region, and then performed all-atom MD simulations for the complexes. The conformational dynamics of HIV-PR with a 20 ns trajectory reveals that the SWCNT can effectively bind to the HIV-1-PR active site and regulate the flap dynamics such as maintaining the flap-flap closed. To gain an insight into the binding affinity, we also performed the MM-PBSA based binding free energy calculations for the four HIV-PR/SWCNT complexes. It was observed that, although the binding between the SWCNT and the HIV-PR decreases due to the mutations, the SWCNTs bind to the HIV-PRs 3-5 folds stronger than the most potent HIV-1-PR inhibitor, TMC114. Remarkably, the significant interactions with binding energy higher than 1kcal/mol focus on the flap and active regions, which favors closing flap-flap and deactivating the active residues of the HIV-PR. The flap dynamics and binding strength information for HIV-PR and SWCNTs can help design SWCNT-based HIV-1-PR inhibitors.

Structure-based methods for predicting target mutation-induced drug resistance and rational drug design to overcome the problem

Drug resistance has become one of the biggest challenges in drug discovery and/or development and has attracted great research interests worldwide. During the past decade, computational strategies have been developed to predict target mutation-induced drug resistance. Meanwhile, various molecular design strategies, including targeting protein backbone, targeting highly conserved residues and dual/multiple targeting, have been used to design novel inhibitors for combating the drug resistance. In this article we review recent advances in development of computational methods for target mutation-induced drug resistance prediction and strategies for rational design of novel inhibitors that could be effective against the possible drug-resistant mutants of the target.

Computational studies of darunavir into HIV-1 protease and DMPC bilayer: necessary conditions for effective binding and the role of the flaps

Human immunodeficiency virus type 1 protease (HIV-1 PR) is one of the main targets toward AIDS therapy. We have selected the potent drug darunavir and a weak inhibitor (fullerene analog) as HIV-1 PR substrates to compare protease's conformational features upon binding. Molecular dynamics (MD), molecular mechanics Poisson-Boltzmann surface area (MM-PBSA), and quantum-mechanical (QM) calculations indicated the importance of the stability of HIV-1 PR flaps toward effective binding: a weak inhibitor may induce flexibility to the flaps, which convert between closed and semiopen states. A water molecule in the darunavir-HIV-1 PR complex bridged the two flap tips of the protease through hydrogen bonding (HB) interactions in a stable structure, a feature that was not observed for the fullerene-HIV-1 PR complex. Additionally, despite that van der Waals interactions and nonpolar contribution to solvation favored permanent fullerene entrapment into the cavity, these interactions alone were not sufficient for effective binding; enhanced electrostatic interactions as observed in the darunavir-complex were the crucial component of the binding energy. An alternative pathway to the usual way of a ligand to access the cavity was also observed for both compounds. Each ligand entered the binding cavity through an opening between the one flap of the protease and a neighboring loop. This suggested that access to the cavity is not necessarily regulated by flap opening. Darunavir exerts its biological action inside the cell, after crossing the membrane barrier. Thus, we also initiated a study on the interactions between darunavir and the DMPC bilayer to reveal that the drug was accommodated inside the bilayer in conformations that resembled its structure into HIV-1 PR, being stabilized via HBs with the lipids and water molecules.

Enamino-oxindole HIV protease inhibitors

We have designed and synthesized a series of HIV protease inhibitors (PIs) with enamino-oxindole substituents optimized to interact with the S2' subsite of the HIV protease binding pocket. Several of these inhibitors have sub-nanomolar K(i) and antiviral IC(50) in the low nM range against WT HIV and against a panel of multi-drug resistant (MDR) strains.

idTarget: a web server for identifying protein targets of small chemical molecules with robust scoring functions and a divide-and-conquer docking approach

Identification of possible protein targets of small chemical molecules is an important step for unravelling their underlying causes of actions at the molecular level. To this end, we construct a web server, idTarget, which can predict possible binding targets of a small chemical molecule via a divide-and-conquer docking approach, in combination with our recently developed scoring functions based on robust regression analysis and quantum chemical charge models. Affinity profiles of the protein targets are used to provide the confidence levels of prediction. The divide-and-conquer docking approach uses adaptively constructed small overlapping grids to constrain the searching space, thereby achieving better docking efficiency. Unlike previous approaches that screen against a specific class of targets or a limited number of targets, idTarget screen against nearly all protein structures deposited in the Protein Data Bank (PDB). We show that idTarget is able to reproduce known off-targets of drugs or drug-like compounds, and the suggested new targets could be prioritized for further investigation. idTarget is freely available as a web-based server at http://idtarget.rcas.sinica.edu.tw.

Potent antiviral HIV-1 protease inhibitor GRL-02031 adapts to the structures of drug resistant mutants with its P1'-pyrrolidinone ring

GRL-02031 (1) is an HIV-1 protease (PR) inhibitor containing a novel P1' (R)-aminomethyl-2-pyrrolidinone group. Crystal structures at resolutions of 1.25-1.55 Å were analyzed for complexes of 1 with the PR containing major drug resistant mutations, PR(I47V), PR(L76V), PR(V82A), and PR(N88D). Mutations of I47V and V82A alter residues in the inhibitor-binding site, while L76V and N88D are distal mutations having no direct contact with the inhibitor. Substitution of a smaller amino acid in PR(I47V) and PR(L76V) and the altered charge of PR(N88D) are associated with significant local structural changes compared to the wild-type PR(WT), while substitution of alanine in PR(V82A) increases the size of the S1' subsite. The P1' pyrrolidinone group of 1 accommodates to these local changes by assuming two different conformations. Overall, the conformation and interactions of 1 with PR mutants resemble those of PR(WT) with similar inhibition constants in good agreement with the antiviral potency on multidrug resistant HIV-1.

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.