First Assessment of Acquired HIV-1 Drug Resistance and Mutation Patterns in Suriname.

HIV drug resistance testing is fundamental in clinical patient management, but data on HIV-1 drug-resistant mutations (DRMs) is scarce in the Caribbean and in Suriname limited to one survey on transmitted resistance. The aim of this study was to address this gap, to gain insight in acquired HIV drug resistance (ADR) prevalence and mutation patterns, and to improve HIV-1 treatment outcome of people living with HIV (PLHIV) in Suriname. A prospective cross-sectional study was conducted from July 2018 through January 2019 among treatment-experienced PLHIV (n = 72), with either treatment failure or antiretroviral therapy restart. Genotypic drug resistance testing was performed and DRM impact on drug effectiveness was examined. Genotypic drug resistance testing revealed 97.2% HIV-1 subtype B, 2.8% B/D recombinants and a ADR prevalence of 63.2% in treatment failure patients, with a predominance of nonnucleoside reverse transcriptase inhibitor (NNRTI) mutations. The most common DRMs were M184V (23.6%) and K103N (18.8%). A high level of non-DRM polymorphisms was observed in both the reverse transcriptase (RT) and protease (PR) gene. Interesting deviations from the existing mutation datasets were noted at position E248 and R83 of the RT gene and L63 and V77 in the PR gene. Full susceptibility to all examined drugs was 54.2%, while high-level drug resistance was estimated at 37.5%, which seems promising for treatment outcomes for PLHIV in Suriname, although cross-resistance to next-generation NNRTIs was already estimated for nearly a quarter of the patients. The meager 2.9% of PR DRMs rendered protease inhibitors as an effective rescue HIV-1 treatment.

HIV drug resistance following a decade of the free antiretroviral therapy programme in India: A review.

OBJECTIVE: The objective of this review was to assess the burden of HIV drug resistance mutations (DRM) in Indian adults exposed to first-line antiretroviral therapy (ART) as per national guidelines. METHODS: An advanced search of the published literature on HIV drug resistance in India was performed in the PubMed and Scopus databases. Data pertaining to age, sex, CD4 count, viral load, and prevalence of nucleoside reverse transcriptase inhibitor (NRTI)/non-nucleoside reverse transcriptase inhibitor (NNRTI) DRM were extracted from each publication. Year-wise Indian HIV-1 reverse transcriptase (RT) sequences were retrieved from the Los Alamos HIV database and mutation analyses were performed. A time trend analysis of the proportion of sequences showing NRTI resistance mutations among individuals exposed to first-line ART was conducted. RESULTS: Overall, 23 studies (1046 unique RT sequences) were identified indicating a prevalence of drug resistance to NRTI and NNRTI. The proportion of RT sequences with any DRM, any NRTI DRM, and any NNRTI DRM was 78.39%, 68.83%, and 73.13%, respectively. The temporal trend analysis of individual DRM from sequences retrieved during 2004-2014 indicated a rising trend in K65R mutations (p=0.013). CONCLUSIONS: Although the overall burden of resistance against first-line ART agents remained steady over the study decade, periodic monitoring is essential. There is the need to develop an HIV-1 subtype C-specific resistance database in India.

Investigation of the solvent-dependent photolysis of a nonnucleoside reverse-transcriptase inhibitor, antiviral agent efavirenz.

This study sought to investigate the solvent-dependency on the photolysis of efavirenz to gain insight into the photoprocesses involved. The primary mechanisms were firstly the excited-state intramolecular proton transfer (i.e. phototautomerization), which generated the imidic acid phototautomer observed as [M-H]- quasimolecular ion at m/z 314.0070 in the high-performance liquid chromatography-electrospray ionization-time-of-flight mass spectrometry in the negative mode. Secondly, the photoinduced α-cleavage with the loss of a carbonyl group occurred (i.e. photodecarbonylation) to form the photoproduct at m/z 286.0395. The ultraviolet-visible spectra illustrated a large, hyperchromic, and slight bathochromic effect in both the π→π* and n→π* electronic transitions. The largest bathochromic effect was prevalent in the chloroform solvent, i.e. chloroform (π* = 0.58; ß = 0.00; α = 0.44) > methanol (π* = 0.60; ß = 0.66; α = 0.98) > acetonitrile (π* = 0.75; ß = 0.40; α = 0.19). This is due to the significant interaction of the amino group with the excited carbonyl moiety which is attributed to intramolecular phototautomerization resulting in a larger energy shift of the electronic state. A plausible explanation is due to the hydrogen bond donor ability of the polar methanol and nonpolar chloroform solvents, which stabilized the polarized imidic acid phototautomer by means of hydrogen bonding interactions, as opposed to the aprotic acetonitrile which exhibits no hydrogen bonding interactions. The study would form the basis for further photolytic analyses and syntheses to generate a plethora of novel photoproducts with anti-HIV activity based on the biologically active benzoxazinone framework of efavirenz.

Performance of an in-house human immunodeficiency virus type 1 genotyping system for assessment of drug resistance in Cuba.

As commercial human immunodeficiency virus type 1 drug resistance assays are expensive, they are not commonly used in resource-limited settings. Hence, a more affordable in-house procedure was set up taking into account the specific epidemiological and economic circumstances of Cuba. The performance characteristics of the in-house assay were evaluated using clinical samples with various subtypes and resistance patterns. The lower limit of amplification was determined on dilutions series of 20 clinical isolates and ranged from 84 to 529 RNA copies/mL. For the assessment of trueness, 14 clinical samples were analyzed and the ViroSeq HIV-1 Genotyping System v2.0 was used as the reference standard. The mean nucleotide sequence identity between the two assays was 98.7% ± 1.0. Additionally, 99.0% of the amino acids at drug resistance positions were identical. The sensitivity and specificity in detecting drug resistance mutations was respectively 94.1% and 99.5%. Only few discordances in drug resistance interpretation patterns were observed. The repeatability and reproducibility were evaluated using 10 clinical samples with 3 replicates per sample. The in-house test was very precise as nucleotide sequence identity among paired nucleotide sequences ranged from 98.7% to 99.9%. The acceptance criteria were met by the in-house test for all performance characteristics, demonstrating a high degree of accuracy. Subsequently, the applicability in routine clinical practice was evaluated on 380 plasma samples. The amplification success rate was 91% and good quality consensus sequences encoding the entire protease and the first 335 codons in reverse transcriptase could be obtained for 99% of the successful amplicons. The reagent cost per sample using the in-house procedure was around € 80 per genotyping attempt. Overall, the in-house assay provided good results, was feasible with equipment and reagents available in Cuba and was half as expensive as commercial assays.

Molecular dynamics and Monte Carlo simulations for protein-ligand binding and inhibitor design.

BACKGROUND: Non-nucleoside inhibitors of HIV reverse transcriptase are an important component of treatment against HIV infection. Novel inhibitors are sought that increase potency against variants that contain the Tyr181Cys mutation. METHODS: Molecular dynamics based free energy perturbation simulations have been run to study factors that contribute to protein-ligand binding, and the results are compared with those from previous Monte Carlo based simulations and activity data. RESULTS: Predictions of protein-ligand binding modes are very consistent for the two simulation methods; the accord is attributed to the use of an enhanced sampling protocol. The Tyr181Cys binding pocket supports large, hydrophobic substituents, which is in good agreement with experiment. CONCLUSIONS: Although some discrepancies exist between the results of the two simulation methods and experiment, free energy perturbation simulations can be used to rapidly test small molecules for gains in binding affinity. GENERAL SIGNIFICANCE: Free energy perturbation methods show promise in providing fast, reliable and accurate data that can be used to complement experiment in lead optimization projects. This article is part of a Special Issue entitled "Recent developments of molecular dynamics".

QSAR models for HEPT derivates as NNRTI inhibitors based on Monte Carlo method.

A series of 107 1-[(2-hydroxyethoxy)-methyl]-6-(phenylthio) thymine (HEPT) with anti-HIV-1 activity as a non-nucleoside reverse transcriptase inhibitor (NNRTI) has been studied. Monte Carlo method has been used as a tool to build up the quantitative structure-activity relationships (QSAR) for anti-HIV-1 activity. The QSAR models were calculated with the representation of the molecular structure by simplified molecular input-line entry system and by the molecular graph. Three various splits into training and test set were examined. Statistical quality of all build models is very good. Best calculated model had following statistical parameters: for training set r(2) = 0.8818, q(2) = 0.8774 and r(2) = 0.9360, q(2) = 0.9243 for test set. Structural indicators (alerts) for increase and decrease of the IC50 are defined. Using defined structural alerts computer aided design of new potential anti-HIV-1 HEPT derivates is presented.

The high cost of fidelity.

The notoriously low fidelity of HIV-1 replication is largely responsible for the virus's rapid mutation rate, facilitating escape from immune or drug control. The error-prone activity of the viral reverse transcriptase (RT) is predicted to be the most influential mechanism for generating mutations. The low fidelity of RT has been successfully exploited by nucleoside and nucleotide analogue reverse transcriptase inhibitors (NRTIs) that halt viral replication upon incorporation. Consequently, drug-resistant strains have arisen in which the viral RT has an increased fidelity of replication, thus reducing analogue incorporation. Higher fidelity, however, impacts on viral fitness. The appearance of compensatory mutations in combination with higher fidelity NRTI resistance mutations and the subsequent reversion of NRTI-resistant mutations upon cessation of antiretroviral treatment lend support to the notion that higher fidelity exacts a fitness cost. Potential mechanisms for reduced viral fitness are a smaller pool of mutant strains available to respond to immune or drug pressure, slower rates of replication, and a limitation to the dNTP tropism of the virus. Unraveling the relationship between replication fidelity and fitness should lead to a greater understanding of the evolution and control of HIV.

Picomolar inhibitors of HIV reverse transcriptase featuring bicyclic replacement of a cyanovinylphenyl group.

Members of the catechol diether class are highly potent non-nucleoside inhibitors of HIV-1 reverse transcriptase (NNRTIs). The most active compounds yield EC50 values below 0.5 nM in assays using human T-cells infected by wild-type HIV-1. However, these compounds such as rilpivirine, the most recently FDA-approved NNRTI, bear a cyanovinylphenyl (CVP) group. This is an uncommon substructure in drugs that gives reactivity concerns. In the present work, computer simulations were used to design bicyclic replacements for the CVP group. The predicted viability of a 2-cyanoindolizinyl alternative was confirmed experimentally and provided compounds with 0.4 nM activity against the wild-type virus. The compounds also performed well with EC50 values of 10 nM against the challenging HIV-1 variant that contains the Lys103Asn/Tyr181Cys double mutation in the RT enzyme. Indolyl and benzofuranyl analogues were also investigated; the most potent compounds in these cases have EC50 values toward wild-type HIV-1 near 10 nM and high-nanomolar activities toward the double-variant. The structural expectations from the modeling were much enhanced by obtaining an X-ray crystal structure at 2.88 Å resolution for the complex of the parent 2-cyanoindolizine 10b and HIV-1 RT. The aqueous solubilities of the most potent indolizine analogues were also measured to be ~40 µg/mL, which is similar to that for the approved drug efavirenz and ~1000-fold greater than for rilpivirine.

Emergence of HIV drug resistance during first- and second-line antiretroviral therapy in resource-limited settings.

INTRODUCTION: Antiretroviral therapy (ART) in resource-limited settings has expanded in the last decade, reaching >8 million individuals and reducing AIDS mortality and morbidity. Continued success of ART programs will require understanding the emergence of HIV drug resistance patterns among individuals in whom treatment has failed and managing ART from both an individual and public health perspective. We review data on the emergence of HIV drug resistance among individuals in whom first-line therapy has failed and clinical and resistance outcomes of those receiving second-line therapy in resource-limited settings. RESULTS: Resistance surveys among patients initiating first-line nonnucleoside reverse-transcriptase inhibitor (NNRTI)-based therapy suggest that 76%-90% of living patients achieve HIV RNA suppression by 12 months after ART initiation. Among patients with detectable HIV RNA at 12 months, HIV drug resistance, primarily due to M184V and NNRTI mutations, has been identified in 60%-72%, although the antiretroviral activity of proposed second-line regimens has been preserved. Complex mutation patterns, including thymidine-analog mutations, K65R, and multinucleoside mutations, are prevalent among cases of treatment failure identified by clinical or immunologic methods. Approximately 22% of patients receiving second-line therapy do not achieve HIV RNA suppression by 6 months, with poor adherence, rather than HIV drug resistance, driving most failures. Major protease inhibitor resistance at the time of second-line failure ranges from 0% to 50%, but studies are limited. CONCLUSIONS: Resistance of HIV to first-line therapy is predictable at 12 months when evaluated by means of HIV RNA monitoring and, when detected, largely preserves second-line therapy options. Optimizing adherence, performing resistance surveillance, and improving treatment monitoring are critical for long-term prevention of drug resistance.

Cost-effectiveness of the third-agent class in treatment-naive human immunodeficiency virus-infected patients in Portugal.

INTRODUCTION: Current Portuguese HIV treatment guidelines recommend initiating antiretroviral therapy with a regimen composed of two Nucleoside Reverse Transcriptase Inhibitors plus one Non-nucleoside Reverse Transcriptase Inhibitor (2NRTI+NNRTI) or two Nucleoside Reverse Transcriptase Inhibitors plus one boosted protease inhibitor (2NRTI+PI/r). Given the lower daily cost of NNRTI as the third agent when compared to the average daily costs of PI/r, it is relevant to estimate the long term impact of each treatment option in the Portuguese context. METHODS: We developed a microsimulation discrete events model for cost-effectiveness analysis of HIV treatment, simulating individual paths from ART initiation to death. Four driving forces determine the course of events: CD4+ cell count, viral load, resistance and adherence. Distributions of time to event are conditional to individuals' characteristics and past history. Time to event was modeled using parametric survival analysis using Stata 11®. Disease progression was structured according to therapy lines and the model was parameterized with cohort Portuguese observational data. All resources were valued at 2009 prices. The National Health Service's perspective was assumed considering a lifetime horizon and a 5% annual discount rate. RESULTS: In this analysis, initiating therapy with two Nucleoside Reverse Transcriptase Inhibitors plus one Non-nucleoside Reverse Transcriptase Inhibitor reduces the average number of switches by 17%, saves 19.573€ per individual and increases life expectancy by 1.7 months showing to be a dominant strategy in 57% of the simulations when compared to two Nucleoside Reverse Transcriptase Inhibitors plus one boosted protease inhibitor. CONCLUSION: This study suggests that, when clinically valid, initiating therapy with two Nucleoside Reverse Transcriptase Inhibitors plus one Non-nucleoside Reverse Transcriptase Inhibitor is a cost-saving strategy and equally effective when compared to two Nucleoside Reverse Transcriptase Inhibitors plus one boosted protease inhibitor as the first regimen.

Monte Carlo sampling and multivariate adaptive regression splines as tools for QSAR modelling of HIV-1 reverse transcriptase inhibitors.

The present work focuses on the development of an interpretable quantitative structure-activity relationship (QSAR) model for predicting the anti-HIV activities of 67 thiazolylthiourea derivatives. This set of molecules has been proposed as potent HIV-1 reverse transcriptase inhibitors (RT-INs). The molecules were encoded to a diverse set of molecular descriptors, spanning different physical and chemical properties. Monte Carlo (MC) sampling and multivariate adaptive regression spline (MARS) techniques were used to select the most important descriptors and to predict the activity of the molecules. The most important descriptor was found to be the aspherisity index. The analysis of variance (ANOVA) and interpretable spline equations showed that the geometrical shape of the molecules has considerable effect on their activities. It seems that the linear molecules are more active than symmetric top compounds. The final MARS model derived displayed a good predictive ability judging from the determination coefficient corresponding to the leave multiple out (LMO) cross-validation technique, i.e. r (2 )= 0.828 (M = 12) and r (2 )= 0.813 (M = 20). The results of this work showed that the developed spline model is robust, has a good predictive power, and can then be used as a reliable tool for designing novel HIV-1 RT-INs.

Transmitted HIV drug resistance among drug-naive subjects recently infected with HIV in Mexico City: a World Health Organization survey to classify resistance and to field test two alternative patient enrollment methods.

In 2004, the World Health Organization performed a survey to assess transmitted drug resistance in Mexico City among drug-naive persons with newly diagnosed human immunodeficiency virus (HIV) infection and likely to be recently infected who were attending 3 voluntary counseling and testing sites. A parallel study comparing 2 alternative methods of enrolling survey participant was conducted in 9 voluntary counseling and testing sites in central Mexico. In study arm 1, subject information, consent and blood specimens were obtained during the HIV diagnostic testing visit. In study arm 2, consent and blood specimens were obtained at the return visit, only from those who were HIV infected. This survey classified nonnucleoside reverse-transcriptase inhibitor and nucleoside reverse-transcriptase inhibitor transmitted drug resistance as <5% and 5%-15%, respectively. Arm 2 yielded major advantages in cost and workload, with no evidence of increased sampling bias.

Rilpivirine: a second-generation nonnucleoside reverse transcriptase inhibitor.

PURPOSE: The pharmacology, pharmacokinetics, drug interactions, clinical efficacy, adverse effects, dosage and administration, and place in therapy of rilpivirine are reviewed. SUMMARY: Rilpivirine is a human immunodeficiency virus type 1 (HIV-1)-specific nonnucleoside reverse transcriptase inhibitor (NNRTI) indicated for use in combination with other antiretroviral agents in adult patients not previously treated with antiretroviral therapy. Rilpivirine is primarily metabolized by the cytochrome P-450 (CYP)3A isoenzyme system. Therefore, providers should be cautious when administering drugs that are inhibitors or inducers of this pathway. Coadministration with CYP 3A inhibitors may lead to increased concentrations of rilpivirine, thereby increasing the risk of adverse effects. Coadministration with inducers of CYP3A isoenzymes or drugs that increase gastric pH may lead to decreased concentrations of rilpivirine, thus promoting virological failure or resistance to rilpivirine. Two Phase III, randomized, double-blind, double-dummy, active-controlled trials compared rilpivirine with efavirenz in HIV-infected adults not previously treated with an antiretroviral. The investigators concluded that rilpivirine, when combined with two nucleoside or nucleotide reverse transcriptase inhibitors, was noninferior to efavirenz for reaching the endpoint of confirmed virological response (HIV-1 RNA level of <50 copies/mL) in adults with HIV infection not previously treated with antiretroviral therapy. The most commonly reported adverse effects included depression, insomnia, headache, and rash. Rilpivirine is administered as a single 25-mg tablet given once daily in combination with other antiretroviral drugs in order to optimize efficacy and reduce resistance. CONCLUSION: Rilpivirine is a viable NNRTI for HIV-infected patients who have not previously received antiretroviral therapy.

Calculation of direct antiretroviral treatment costs and potential cost savings by using generics in the German HIV ClinSurv cohort.

UNLABELLED: BACKGROUND/AIM OF THE STUDY: The study aimed to determine the cost impacts of antiretroviral drugs by analysing a long-term follow-up of direct costs for combined antiretroviral therapy, cART, -regimens in the nationwide long-term observational multi-centre German HIV ClinSurv Cohort. The second aim was to develop potential cost saving strategies by modelling different treatment scenarios. METHODS: Antiretroviral regimens (ART) from 10,190 HIV-infected patients from 11 participating ClinSurv study centres have been investigated since 1996. Biannual data cART-initiation, cART-changes, surrogate markers, clinical events and the Centre of Disease Control- (CDC)-stage of HIV disease are reported. Treatment duration was calculated on a daily basis via the documented dates for the beginning and end of each antiretroviral drug treatment. Prices were calculated for each individual regimen based on actual office sales prices of the branded pharmaceuticals distributed by the license holder including German taxes. RESULTS: During the 13-year follow-up period, 21,387,427 treatment days were covered. Cumulative direct costs for antiretroviral drugs of €812,877,356 were determined according to an average of €42.08 per day (€7.52 to € 217.70). Since cART is widely used in Germany, the costs for an entire regimen increased by 13.5%. Regimens are more expensive in the advanced stages of HIV disease. The potential for cost savings was calculated using non-nucleotide-reverse-transcriptase-inhibitor, NNRTI, more frequently instead of ritonavir-boosted protease inhibitor, PI/r, in first line therapy. This calculation revealed cumulative savings of 10.9% to 19.8% of daily treatment costs (50% and 90% substitution of PI/r, respectively). Substituting certain branded drugs by generic drugs showed potential cost savings of between 1.6% and 31.8%. CONCLUSIONS: Analysis of the data of this nationwide study reflects disease-specific health services research and will give insights into the cost impacts of antiretroviral therapy, and might allow a more rational allocation of resources within the German health care system.

Replicative fitness costs of nonnucleoside reverse transcriptase inhibitor drug resistance mutations on HIV subtype C.

Single-dose nevirapine (NVP) is quite effective in preventing transmission of the human immunodeficiency virus (HIV) from mother to child; however, many women develop resistance to NVP in this setting. Comparing outcomes of clinical studies reveals an increased amount of resistance in subtype C relative to that in other subtypes. This study investigates how nonnucleoside reverse transcriptase inhibitor (NNRTI) drug resistance mutations of subtype C affect replication capacity. The 103N, 106A, 106M, 181C, 188C, 188L, and 190A drug resistance mutations were placed in a reverse transcriptase (RT) that matches the consensus subtype C sequence as well as the HXB2 RT, as a subtype B reference. The replicative fitness of each mutant was compared with that of the wild type in a head-to-head competition assay. The 106A mutant of subtype C would not grow in the competition assay, making it the weakest virus tested. The effect of the 106M mutation was weaker than those of the 181C and 188C mutations in the consensus C RT, but in subtype B, this difference was not seen. To see if the 106A mutation in a different subtype C background would have a different replicative profile, the same NNRTI resistance mutations were added to the MJ4 RT, a reference subtype C molecular clone. In the context of MJ4 RT, the 106A mutant was not the only mutant that showed poor replicative fitness; the 106M, 188C, and 190A mutants also failed to replicate. These results suggest that NNRTIs may be a cost-effective alternative for salvage therapy if deleterious mutations are present in a subtype C setting.

Evolutionary dynamics of complex networks of HIV drug-resistant strains: the case of San Francisco.

Over the past two decades, HIV resistance to antiretroviral drugs (ARVs) has risen to high levels in the wealthier countries of the world, which are able to afford widespread treatment. We have gained insights into the evolution and transmission dynamics of ARV resistance by designing a biologically complex multistrain network model. With this model, we traced the evolutionary history of ARV resistance in San Francisco and predict its future dynamics. By using classification and regression trees, we identified the key immunologic, virologic, and treatment factors that increase ARV resistance. Our modeling shows that 60% of the currently circulating ARV-resistant strains in San Francisco are capable of causing self-sustaining epidemics, because each individual infected with one of these strains can cause, on average, more than one new resistant infection. It is possible that a new wave of ARV-resistant strains that pose a substantial threat to global public health is emerging.

Effect of prophylactic lamivudine for chemotherapy-associated hepatitis B reactivation in lymphoma: a meta-analysis of published clinical trials and a decision tree addressing prolonged prophylaxis and maintenance.

Lamivudine prophylaxis is an effective strategy in HbSAg-positive patients receiving cancer chemotherapy. Recent data indicate that a lamividune-prophylaxis strategy results in a decrease of hepatitis B virus (HBV) reactivation rates, though its effect on HBV-mortality remains equivocal. This report evaluates the benefits from this strategy among lymphoma patients and develops a management approach for patients with prolonged immunosuppression. A Medline search was conducted to retrieve published trials on HBsAg-positive lymphoma patients receiving prophylactic lamivudine during chemotherapy. Basic inclusion criterion was to report HBV-reactivation rates with and without lamivudine prophylaxis. A meta-analysis of the risk of HBV-reactivation and HBV-related mortality was conducted, and the pooled effect was calculated as risk ratio (RR). We found that lamivudine prophylaxis is associated with a significant reduction in hepatitis B virus reactivation (RR 0.21, 95%CI 0.13-0.35) and a trend in reducing HBV-related mortality (RR 0.68, 95%CI 0.19-2.49). In order to study the long-term effects of anti-HBV prophylaxis when prolonged immunosuppression is needed, we used our findings to model a decision tree. Overall survival was the main outcome used in the analysis. Rituximab maintenance in B-cell lymphomas was used as a paradigm of prolonged immunosuppression. We found that extended anti-HBV prophylaxis can improve survival rates by 2.4% in HBsAg-positive patients. If 1,000 HBsAg-positive lymphoma patients receive prophylaxis, one will die from hepatitis B virus reactivation versus 25/1,000 if no prophylaxis is administered. This effect is probably mediated through a reduction of hepatitis B virus reactivation and HBV-related mortality. The ideal antiviral agent needs to be determined.

Efficient drug lead discovery and optimization.

During the 1980s, advances in the abilities to perform computer simulations of chemical and biomolecular systems and to calculate free energy changes led to the expectation that such methodology would soon show great utility for guiding molecular design. Important potential applications included design of selective receptors, catalysts, and regulators of biological function including enzyme inhibitors. This time also saw the rise of high-throughput screening and combinatorial chemistry along with complementary computational methods for de novo design and virtual screening including docking. These technologies appeared poised to deliver diverse lead compounds for any biological target. As with many technological advances, realization of the expectations required significant additional effort and time. However, as summarized here, striking success has now been achieved for computer-aided drug lead generation and optimization. De novo design using both molecular growing and docking are illustrated for lead generation, and lead optimization features free energy perturbation calculations in conjunction with Monte Carlo statistical mechanics simulations for protein-inhibitor complexes in aqueous solution. The specific applications are to the discovery of non-nucleoside inhibitors of HIV reverse transcriptase (HIV-RT) and inhibitors of the binding of the proinflammatory cytokine MIF to its receptor CD74. A standard protocol is presented that includes scans for possible additions of small substituents to a molecular core, interchange of heterocycles, and focused optimization of substituents at one site. Initial leads with activities at low-micromolar concentrations have been advanced rapidly to low-nanomolar inhibitors.

Optimization of azoles as anti-human immunodeficiency virus agents guided by free-energy calculations.

Efficient optimization of an inactive 2-anilinyl-5-benzyloxadiazole core has been guided by free energy perturbation (FEP) calculations to provide potent non-nucleoside inhibitors of human immunodeficiency virus (HIV) reverse transcriptase (NNRTIs). An FEP "chlorine scan" was performed to identify the most promising sites for substitution of aryl hydrogens. This yielded NNRTIs 8 and 10 with activities (EC50) of 820 and 310 nM for protection of human T-cells from infection by wild-type HIV-1. FEP calculations for additional substituent modifications and change of the core heterocycle readily led to oxazoles 28 and 29, which were confirmed as highly potent anti-HIV agents with activities in the 10-20 nM range. The designed compounds were also monitored for possession of desirable pharmacological properties by use of additional computational tools. Overall, the trends predicted by the FEP calculations were well borne out by the assay results. FEP-guided lead optimization is confirmed as a valuable tool for molecular design including drug discovery; chlorine scans are particularly attractive since they are both straightforward to perform and highly informative.