Crystal structures of 8-Cl and 9-Cl TIBO complexed with wild-type HIV-1 RT and 8-Cl TIBO complexed with the Tyr181Cys HIV-1 RT drug-resistant mutant

Kinetic coevolutionary models predict the temporal emergence of HIV-1 resistance mutations under drug selection pressure

Drug resistance in HIV type 1 (HIV-1) is a pervasive problem that affects the lives of millions of people worldwide. Although records of drug-resistant mutations (DRMs) have been extensively tabulated within public repositories, our understanding of the evolutionary kinetics of DRMs and how they evolve together remains limited. Epistasis, the interaction between a DRM and other residues in HIV-1 protein sequences, is key to the temporal evolution of drug resistance. We use a Potts sequence-covariation statistical-energy model of HIV-1 protein fitness under drug selection pressure, which captures epistatic interactions between all positions, combined with kinetic Monte-Carlo simulations of sequence evolutionary trajectories, to explore the acquisition of DRMs as they arise in an ensemble of drug-naive patient protein sequences. We follow the time course of 52 DRMs in the enzymes protease, RT, and integrase, the primary targets of antiretroviral therapy. The rates at which DRMs emerge are highly correlated with their observed acquisition rates reported in the literature when drug pressure is applied. This result highlights the central role of epistasis in determining the kinetics governing DRM emergence. Whereas rapidly acquired DRMs begin to accumulate as soon as drug pressure is applied, slowly acquired DRMs are contingent on accessory mutations that appear only after prolonged drug pressure. We provide a foundation for using computational methods to determine the temporal evolution of drug resistance using Potts statistical potentials, which can be used to gain mechanistic insights into drug resistance pathways in HIV-1 and other infectious agents.

Selected Milestones in Antiviral Drug Development

This review article will describe the (wide) variety of approaches that I envisaged to develop a specific therapy for viral infections: (i) interferon and its inducers, (ii) HSV, VZV and CMV inhibitors, (iii) NRTIs (nucleoside reverse transcriptase inhibitors), NtRTIs (nucleotide reverse transcriptase inhibitors) and NNRTIs (non-nucleoside reverse transcriptase inhibitors) as HIV inhibitors, (iv) NtRTIs as HBV inhibitors, and finally, (v) the transition of an HIV inhibitor to a stem cell mobilizer, as exemplified by AMD-3100 (Mozobil®).

Machine intelligence-guided selection of optimized inhibitor for human immunodeficiency virus (HIV) from natural products

The human immunodeficiency virus (HIV) connects to the cluster of differentiation (CD4) and any of the entry co-receptors (CCR5 and CXCR4); followed by unloading the viral genome, reverse transcriptase, and integrase enzymes within the host cell. The co-receptors facilitate the entry of virus and vital enzymes, leading to replication and pre-maturation of viral particles within the host. The protease enzyme transforms the immature viral vesicles into the mature virion. The pivotal role of co-receptors and enzymes in homeostasis and growth makes the crucial target for anti-HIV drug discovery, and the availability of X-ray crystal structures is an asset. Here, we used the machine intelligence-driven framework (A-HIOT) to identify and optimize target-based potential hit molecules for five significant protein targets from the ZINC15 database (natural products dataset). Following validation with dynamic motion behavior analysis and molecular dynamics simulation, the optimized hits were evaluated using in silico ADMET filtration. Furthermore, three molecules were screened, optimized, and validated: ZINC00005328058 for CCR5 and protease, ZINC000254014855 for CXCR4 and integrase, and ZINC000000538471 for reverse transcriptase. In clinical trials, the ZINC000254014855 and ZINC000254014855 were passed in primary screens for vif-HIV-1, and we reported the specific receptor as well as interactions. As a result, the validated molecules may be investigated further in experimental studies targeting specific receptors in order to design and synergize an anti-HIV regimen.

Elucidation of Antiviral and Antioxidant Potential of C-Phycocyanin against HIV-1 Infection through In Silico and In Vitro Approaches

Antiretroviral therapy is the single existing therapy for patients infected with HIV; however, it has drawbacks in terms of toxicity and resistance. Thus, there is a continuous need to explore safe and efficacious anti-retroviral agents. C-Phycocyanin (C-PC) is a phycobiliprotein, which has been known for various biological properties; however, its effect on HIV-1 replication needs revelation. This study aimed to identify the inhibitory effects of C-PC on HIV-1 using in vitro and in silico approaches and to assess its role in the generation of mitochondrial reactive oxygen species (ROS) during HIV-1 infection. In vitro anti-HIV-1 activity of C-PC was assessed on TZM-bl cells through luciferase gene assay against four different clades of HIV-1 strains in a dose-dependent manner. Results were confirmed in PBMCs, using the HIV-1 p24 antigen assay. Strong associations between C-PC and HIV-1 proteins were observed through in silico molecular simulation-based interactions, and the in vitro mechanistic study confirmed its target by inhibition of reverse transcriptase and protease enzymes. Additionally, the generation of mitochondrial ROS was detected by the MitoSOX and DCF-DA probe through confocal microscopy. Furthermore, our results confirmed that C-PC treatment notably subdued the fluorescence in the presence of the virus, thus reduction of ROS and the activation of caspase-3/7 in HIV-1-infected cells. Overall, our study suggests C-PC as a potent and broad in vitro antiviral and antioxidant agent against HIV-1 infection.

Insights into HIV-1 Reverse Transcriptase (RT) Inhibition and Drug Resistance from Thirty Years of Structural Studies

The enzyme reverse transcriptase (RT) plays a central role in the life cycle of human immunodeficiency virus (HIV), and RT has been an important drug target. Elucidations of the RT structures trapping and detailing the enzyme at various functional and conformational states by X-ray crystallography have been instrumental for understanding RT activities, inhibition, and drug resistance. The structures have contributed to anti-HIV drug development. Currently, two classes of RT inhibitors are in clinical use. These are nucleoside/nucleotide reverse transcriptase inhibitors (NRTIs) and non-nucleoside reverse transcriptase inhibitors (NNRTIs). However, the error-prone viral replication generates variants that frequently develop resistance to the available drugs, thus warranting a continued effort to seek more effective treatment options. RT also provides multiple additional potential druggable sites. Recently, the use of single-particle cryogenic electron microscopy (cryo-EM) enabled obtaining structures of NNRTI-inhibited HIV-1 RT/dsRNA initiation and RT/dsDNA elongation complexes that were unsuccessful by X-ray crystallography. The cryo-EM platform for the structural study of RT has been established to aid drug design. In this article, we review the roles of structural biology in understanding and targeting HIV RT in the past three decades and the recent structural insights of RT, using cryo-EM.

Review on fluorinated nucleoside/non-nucleoside FDA-approved antiviral drugs

FDA-approved antiviral agents represent an important class that has attracted attention in recent years to combat current and future threats of viral pandemics. Fluorine ameliorates the electronic, lipophilic and steric problems of drugs. Additionally, fluorine can prolong drug activity and improve metabolic stability, thereby, modifying their pharmacodynamic and pharmacokinetic character. Herein, we summarized the fluorinated FDA-approved antiviral agents, dealing with biological aspects, mechanisms of action, and synthetic pathways.

FDA-approved antiviral agents represent an important class that has attracted attention in recent years to combat current and future threats of viral pandemics.

Avoiding Drug Resistance in HIV Reverse Transcriptase

HIV reverse transcriptase (RT) is an enzyme that plays a major role in the replication cycle of HIV and has been a key target of anti-HIV drug development efforts. Because of the high genetic diversity of the virus, mutations in RT can impart resistance to various RT inhibitors. As the prevalence of drug resistance mutations is on the rise, it is necessary to design strategies that will lead to drugs less susceptible to resistance. Here we provide an in-depth review of HIV reverse transcriptase, current RT inhibitors, novel RT inhibitors, and mechanisms of drug resistance. We also present novel strategies that can be useful to overcome RT's ability to escape therapies through drug resistance. While resistance may not be completely avoidable, designing drugs based on the strategies and principles discussed in this review could decrease the prevalence of drug resistance.

Structure-activity relationship, molecular docking, and molecular dynamic studies of diterpenes from marine natural products with anti-HIV activity

HIV-1 infection is a global epidemic whose treatment is limited majorly by viral resistance and adverse effects. Natural products from algae have been studied for many years, including antiviral, being an alternative to anti-HIV drug design. Since the isolation of natural products can be a hurdle, molecular modeling is an important tool to study these compounds. Herein, structure-activity relationship, molecular docking, and molecular dynamic studies were performed to direct the studies of ten marine natural products with anti-HIV activity. In the structure-activity relationship, descriptors were identified associating the anti-HIV activity of five diterpenes with possible action on the reverse transcriptase allosteric site. These diterpenes were evaluated by molecular docking, and it was identified that only dolabelladienetriol interacted in the allosteric site. Molecular dynamics suggested that the dolabelladienetriol might interfere with the viral RNA binding to HIV-1 RT by inducing a conformational change of the enzyme. Also, in silico ADMET simulations predicts that the dolabelladienetriol present a high potential to be successfully developed as a drug. Thus, applying in silico approaches was possible to suggest potential anti-HIV compounds derived from marine natural products.Communicated by Ramaswamy H. Sarma.

The Impact of HIV-1 Drug Escape on the Global Treatment Landscape

The rising prevalence of HIV drug resistance (HIVDR) could threaten gains made in combating the HIV epidemic and compromise the 90-90-90 target proposed by United Nations Programme on HIV/AIDS (UNAIDS) to have achieved virological suppression in 90% of all persons receiving antiretroviral therapy (ART) by the year 2020. HIVDR has implications for the persistence of HIV, the selection of current and future ART drug regimens, and strategies of vaccine and cure development. Focusing on drug classes that are in clinical use, this Review critically summarizes what is known about the mechanisms the virus utilizes to escape drug control. Armed with this knowledge, strategies to limit the expansion of HIVDR are proposed.

Polymerases of Coronaviruses

Coronaviruses (CoVs) are a major group of viruses known to be responsible for wide spectrum of diseases in multiple species. The CoVs affecting human population are referred to as human coronaviruses (HCoVs). They lead to multiple respiratory diseases, such as common cold, pneumonia, bronchitis, severe acute respiratory syndrome, and Middle East respiratory syndrome. CoVs are RNA viruses that require RNA-dependent RNA polymerases (RdRPs) for various steps in their life cycle. Action of RdRP is needed in several steps in the life cycle of CoVs and thus RdRPs constitute potential targets for drugs and other therapeutic interventions for the treatment of diseases caused by CoVs. The chapter therefore presents a detailed discussion on the structure and functions of CoV polymerases and the development of their potential inhibitors.

Torsional flexibility of undecorated catechol diether compound as potent NNRTI targeting HIV-1 reverse transcriptase

Conformational adaptation of non-nucleoside reverse transcriptase inhibitor (NNRTI) via torsional flexibility is found to be very significant for targeting human immunodeficiency virus type 1 reverse transcriptase (HIV-1 RT) mutants. Catechol diether derivative including flexible torsions is new potent NNRTI with picomolar activity. Moreover, this derivative also reveals the good solubility, low toxicity and potent inhibition for HIV-1 mutants. In this study, torsional flexibility of an undecorated catechol diether compound in the binding pocket of wild type and mutants (Y181C and K103N/Y181C) HIV-1 RT is investigated by using QM/MM calculations. From the results, the uracil ring is found to exhibit more flexibility in the NNIBP. On the contrary, potential energy surfaces show that high energy is encountered by changing of the corresponding torsion of the cyanovinyl aryl ring indicating the limitation for torsional flexibility. For pointing out the key interaction for the binding, the residual interaction energies are performed by means of QM calculations. Important attractive interactions through hydrogen bonds between the inhibitor and K102, K/N103, V106, and Y188 are observed. The catechol ring is proposed to be modified in order to strengthen interactions with surrounding amino acids. The results may help for the designing of new potent NNRTIs.

Emerging reverse transcriptase inhibitors for HIV-1 infection

INTRODUCTION

There are 36.7 million people living with HIV with 20.9 million having access to antiretroviral therapy (ART). Nucleoside/nucleotide reverse transcriptase inhibitors (NRTIs) remain the 'backbone' of ART. However, the currently available nine NRTIs and five non-nucleoside reverse transcriptase inhibitors (NNRTIs) have significant side effects and resistance profiles. Areas covered: We summarize the mechanisms of resistance and other limitations of the existing NRTIs/NNRTIs. GS-9131, MK-8591, Elsulfavirine and Doravirine are four new agents that are furthest along in development. Expert opinion: ART development has evolved with several new promising agents. Longer-acting agents, like MK-8591 are extremely attractive to enhance drug adherence and patient satisfaction. Doravirine offers an NNRTI effective against common mutations that has fewer side effects, limitations on dosing and drug interactions. GS-9131 is very potent and active against a variety of NRTI mutants but it is too early in its development to understand its full risks and benefits. Finally, Elsulfavirine has a long half-life and preliminary data suggests fewer side effects than the most commonly used NNRTI, efavirenz. Each of these new agents shows promise and potential to improve ART in the future. The newer generation of reverse transcriptase inhibitors have longer half-lives, more favorable adverse effect profiles, and fewer drug interactions.

Molecular insight on the binding of NNRTI to K103N mutated HIV-1 RT: molecular dynamics simulations and dynamic pharmacophore analysis

Regardless of advances in anti-HIV therapy, HIV infection remains an immense challenge due to the rapid onset of mutation instigating drug resistance. Rilpivirine is a second generation di-aryl pyrimidine (DAPY) derivative, known to effectively inhibit wild-type (WT) as well as various mutant HIV-1 reverse transcriptase (RT). In this study, a cumulative 240 ns of molecular dynamic (MD) simulations of WT HIV-1 RT and its corresponding K103N mutated form, complexed with rilpivirine, were performed in solution. Conformational analysis of the NNRTI inside the binding pocket (NNIBP) revealed the ability of rilpivirine to adopt different conformations, which is possibly the reason for its reasonable activity against mutant HIV-1 RT. Binding free energy (MM-PB/GB SA) calculations of rilpivirine with mutant HIV-1 RT are in agreement with experimental data. The dynamics of interaction patterns were investigated based on the MD simulations using dynophores, a novel approach for MD-based ligand-target interaction mapping. The results from this interaction profile analysis suggest an alternate interaction between the linker N atom of rilpivirine and Lys 101, potentially providing the stability for ligand binding. PCA analysis and per residue fluctuation has highlighted the significant role of flexible thumb and finger sub-domains of RT in its biological activity. This study investigated the underlying reason for rilpivirine's improved inhibitory profile against mutant RT, which could be helpful to understand the molecular basis of HIV-1 RT drug resistance and design novel NNRTIs with improved drug resistance tolerance.

CHARMM-GUI ligand reader and modeler for CHARMM force field generation of small molecules

Reading ligand structures into any simulation program is often nontrivial and time consuming, especially when the force field parameters and/or structure files of the corresponding molecules are not available. To address this problem, we have developed Ligand Reader & Modeler in CHARMM-GUI. Users can upload ligand structure information in various forms (using PDB ID, ligand ID, SMILES, MOL/MOL2/SDF file, or PDB/mmCIF file), and the uploaded structure is displayed on a sketchpad for verification and further modification. Based on the displayed structure, Ligand Reader & Modeler generates the ligand force field parameters and necessary structure files by searching for the ligand in the CHARMM force field library or using the CHARMM general force field (CGenFF). In addition, users can define chemical substitution sites and draw substituents in each site on the sketchpad to generate a set of combinatorial structure files and corresponding force field parameters for throughput or alchemical free energy simulations. Finally, the output from Ligand Reader & Modeler can be used in other CHARMM-GUI modules to build a protein-ligand simulation system for all supported simulation programs, such as CHARMM, NAMD, GROMACS, AMBER, GENESIS, LAMMPS, Desmond, OpenMM, and CHARMM/OpenMM. Ligand Reader & Modeler is available as a functional module of CHARMM-GUI at http://www.charmm-gui.org/input/ligandrm. © 2017 Wiley Periodicals, Inc.

N-Pyrrylarylsulfones with High Therapeutic Potential

This review illustrates the various studies made to investigate the activity of N-pyrrylarylsulfone containing compounds as potential antiviral, anticancer and SNC drugs. A number of synthetic approaches to obtain tetracyclic, tricyclic and non-cyclic compounds, and their biological activity with regard to structure-activity relationships (SARs) have been reviewed. The literature reviewed here may provide useful information on the potential of N-pyrrylarylsulfone pharmacophore as well as suggest concepts for the design and synthesis of new N-pyrrylarylsulfone based agents.

Mechanistic Study of Common Non-Nucleoside Reverse Transcriptase Inhibitor-Resistant Mutations with K103N and Y181C Substitutions

Non-nucleoside reverse transcriptase inhibitors (NNRTIs) are a mainstay of therapy for human immunodeficiency type 1 virus (HIV-1) infections. However, their effectiveness can be hampered by the emergence of resistant mutations. To aid in designing effective NNRTIs against the resistant mutants, it is important to understand the resistance mechanism of the mutations. Here, we investigate the mechanism of the two most prevalent NNRTI-associated mutations with K103N or Y181C substitution. Virus and reverse transcriptase (RT) with K103N/Y188F, K103A, or K103E substitutions and with Y181F, Y188F, or Y181F/Y188F substitutions were employed to study the resistance mechanism of the K103N and Y181C mutants, respectively. Results showed that the virus and RT with K103N/Y188F substitutions displayed similar resistance levels to the virus and RT with K103N substitution versus NNRTIs. Virus and RT containing Y181F, Y188F, or Y181F/Y188F substitution exhibited either enhanced or similar susceptibility to NNRTIs compared with the wild type (WT) virus. These results suggest that the hydrogen bond between N103 and Y188 may not play an important role in the resistance of the K103N variant to NNRTIs. Furthermore, the results from the studies with the Y181 or Y188 variant provide the direct evidence that aromatic π-π stacking plays a crucial role in the binding of NNRTIs to RT.

Conformational Plasticity of the NNRTI-Binding Pocket in HIV-1 Reverse Transcriptase: A Fluorine Nuclear Magnetic Resonance Study

HIV-1 reverse transcriptase (RT) is a major drug target in the treatment of HIV-1 infection. RT inhibitors currently in use include non-nucleoside, allosteric RT inhibitors (NNRTIs), which bind to a hydrophobic pocket, distinct from the enzyme's active site. We investigated RT-NNRTI interactions by solution (19)F nuclear magnetic resonance (NMR), using singly (19)F-labeled RT proteins. Comparison of (19)F chemical shifts of fluorinated RT and drug-resistant variants revealed that the fluorine resonance is a sensitive probe for identifying mutation-induced changes in the enzyme. Our data show that in the unliganded enzyme, the NNRTI-binding pocket is highly plastic and not locked into a single conformation. Upon inhibitor binding, the binding pocket becomes rigidified. In the inhibitor-bound state, the (19)F signal of RT is similar to that of drug-resistant mutant enzymes, distinct from what is observed for the free state. Our results demonstrate the power of (19)F NMR spectroscopy to characterize conformational properties using selectively (19)F-labeled protein.

Approved Antiviral Drugs over the Past 50 Years

Since the first antiviral drug, idoxuridine, was approved in 1963, 90 antiviral drugs categorized into 13 functional groups have been formally approved for the treatment of the following 9 human infectious diseases: (i) HIV infections (protease inhibitors, integrase inhibitors, entry inhibitors, nucleoside reverse transcriptase inhibitors, nonnucleoside reverse transcriptase inhibitors, and acyclic nucleoside phosphonate analogues), (ii) hepatitis B virus (HBV) infections (lamivudine, interferons, nucleoside analogues, and acyclic nucleoside phosphonate analogues), (iii) hepatitis C virus (HCV) infections (ribavirin, interferons, NS3/4A protease inhibitors, NS5A inhibitors, and NS5B polymerase inhibitors), (iv) herpesvirus infections (5-substituted 2'-deoxyuridine analogues, entry inhibitors, nucleoside analogues, pyrophosphate analogues, and acyclic guanosine analogues), (v) influenza virus infections (ribavirin, matrix 2 protein inhibitors, RNA polymerase inhibitors, and neuraminidase inhibitors), (vi) human cytomegalovirus infections (acyclic guanosine analogues, acyclic nucleoside phosphonate analogues, pyrophosphate analogues, and oligonucleotides), (vii) varicella-zoster virus infections (acyclic guanosine analogues, nucleoside analogues, 5-substituted 2'-deoxyuridine analogues, and antibodies), (viii) respiratory syncytial virus infections (ribavirin and antibodies), and (ix) external anogenital warts caused by human papillomavirus infections (imiquimod, sinecatechins, and podofilox). Here, we present for the first time a comprehensive overview of antiviral drugs approved over the past 50 years, shedding light on the development of effective antiviral treatments against current and emerging infectious diseases worldwide.

Relating the Structure of HIV-1 Reverse Transcriptase to Its Processing Step

Abstract By treating an enzyme as a coarse-grained uniform block of material, utilizing only the α-Carbon positions, the normal modes of motion can be obtained. For reverse transcriptase the slower of these motions are suggestive of being involved in the processing step, where the RNA or DNA strand is copied onto a new DNA strand at a polymerase site, and the RNA strand is subsequently cut up at the distant Ribonuclease H site. The slowest mode of motion involves hinge bending about a site midway between the polymerase and Ribonuclease H sites, suggesting that it can push or pull the RNA strand between these two sites. Pulling the nucleic acid strand would require tight binding to the RNase H site. The next slowest mode involves a hinge that opens and closes the protein like a clamp, which could facilitate the release of the nucleic acids for their step-wise progression. The third mode could rotate the substrate. An overall description of the step-wise processing step would involve close coordination among these steps. Results suggest that the smaller p51 subunit serves only as ballast to support the various modes of motion involving the different parts of the p66 subunit.

Reverse Transcription of Retroviruses and LTR Retrotransposons

The enzyme reverse transcriptase (RT) was discovered in retroviruses almost 50 years ago. The demonstration that other types of viruses, and what are now called retrotransposons, also replicated using an enzyme that could copy RNA into DNA came a few years later. The intensity of the research in both the process of reverse transcription and the enzyme RT was greatly stimulated by the recognition, in the mid-1980s, that human immunodeficiency virus (HIV) was a retrovirus and by the fact that the first successful anti-HIV drug, azidothymidine (AZT), is a substrate for RT. Although AZT monotherapy is a thing of the past, the most commonly prescribed, and most successful, combination therapies still involve one or both of the two major classes of anti-RT drugs. Although the basic mechanics of reverse transcription were worked out many years ago, and the first high-resolution structures of HIV RT are now more than 20 years old, we still have much to learn, particularly about the roles played by the host and viral factors that make the process of reverse transcription much more efficient in the cell than in the test tube. Moreover, we are only now beginning to understand how various host factors that are part of the innate immunity system interact with the process of reverse transcription to protect the host-cell genome, the host cell, and the whole host, from retroviral infection, and from unwanted retrotransposition.

Exploring isoxazole and carboxamide derivatives as potential non-nucleoside reverse transcriptase inhibitors

Nonnucleoside reverse transciptase inhibitors (NNRTI) are a class of drug molecules with a specific target of HIV-1 reverse transcriptase (RT). In the present work, we evaluated a set of selected oxazole and carboxamide derivatives to identify potential pharmacophoric features using molecular docking approach. The docking approach employed has been validated by enrichment factor calculation at top 1% (EF1%). It shows a considerable improvement in EF1%value compared to earlier reported study carried out on specific dataset of ligands and decoys for RT, in the directory of useful decoys (DUD). The carboxamide derivatives show better activity as NNRT inhibitors than oxazole derivatives. From this study, four pharmacophoric groups including a triazine ring, an aniline substituent, a benzyl amide moiety and a trimethylphenoxy substituent have been recognized and used for designing new NNRT inhibitors. Newly designed molecules show significant enhancement in docking scores over the native ligand, parent and other training set molecules. In addition, some functional groups have also been identified to assist in improving the activity of these pharmacophores. Thus a nitrile group, an amide and fluoro substitution turn out to be an important requisite for NNRT potential inhibitors.

Rilpivirine analogs potently inhibit drug-resistant HIV-1 mutants

Background

Nonnucleoside reverse transcriptase inhibitors (NNRTIs) are a class of antiretroviral compounds that bind in an allosteric binding pocket in HIV-1 RT, located about 10 Å from the polymerase active site. Binding of an NNRTI causes structural changes that perturb the alignment of the primer terminus and polymerase active site, preventing viral DNA synthesis. Rilpivirine (RPV) is the most recent NNRTI approved by the FDA, but like all other HIV-1 drugs, suboptimal treatment can lead to the development of resistance. To generate better compounds that could be added to the current HIV-1 drug armamentarium, we have developed several RPV analogs to combat viral variants that are resistant to the available NNRTIs.

Results

Using a single-round infection assay, we identified several RPV analogs that potently inhibited a broad panel of NNRTI resistant mutants. Additionally, we determined that several resistant mutants selected by either RPV or Doravirine (DOR) caused only a small increase in susceptibility to the most promising RPV analogs.

Conclusions

The antiviral data suggested that there are RPV analogs that could be candidates for further development as NNRTIs, and one of the most promising compounds was modeled in the NNRTI binding pocket. This model can be used to explain why this compound is broadly effective against the panel of NNRTI resistance mutants.

Electronic supplementary material

The online version of this article (doi:10.1186/s12977-016-0244-2) contains supplementary material, which is available to authorized users.

Focus on Chirality of HIV-1 Non-Nucleoside Reverse Transcriptase Inhibitors

Chiral HIV-1 non-nucleoside reverse transcriptase inhibitors (NNRTIs) are of great interest since one enantiomer is often more potent than the corresponding counterpart against the HIV-1 wild type (WT) and the HIV-1 drug resistant mutant strains. This review exemplifies the various studies made to investigate the effect of chirality on the antiretroviral activity of top HIV-1 NNRTI compounds, such as nevirapine (NVP), efavirenz (EFV), alkynyl- and alkenylquinazolinone DuPont compounds (DPC), diarylpyrimidine (DAPY), dihydroalkyloxybenzyloxopyrimidine (DABO), phenethylthiazolylthiourea (PETT), indolylarylsulfone (IAS), arylphosphoindole (API) and trifluoromethylated indole (TFMI) The chiral separation, the enantiosynthesis, along with the biological properties of these HIV-1 NNRTIs, are discussed.

In-silico Studies Show Potent Inhibition of HIV-1 Reverse Transcriptase Activity by a Herbal Drug

Acquired immunodeficiency syndrome (AIDS) is a life threatening disease of the human immune system caused by human immunodeficiency virus (HIV). Effective inhibition of reverse transcriptase activity is a prominent, clinically viable approach for the treatment of AIDS. Few non-nucleoside reverse transcriptase inhibitors (NNRTIs) have been approved by the United States Food and Drug Administration (US FDA) as drugs for AIDS. In order to enhance therapeutic options against AIDS we examined novel herbal compounds of 4-thiazolidinone and its derivatives that are known to have remarkable antiviral potency. Our molecular docking and simulation experiments have identified one such herbal molecule known as (5E)-3-(2-aminoethyl)-5-benzylidene-1, 3-thiazolidine-2,4-dione that may bind HIV-1RT with high affinity to cause noncompetitive inhibition. Results are also compared with other US FDA approved drugs. Long de novo simulations and docking study suggest that the ligand (5E)-3-(2-aminoethyl)-5-benzylidene-1, 3-thiazolidine-2,4-dione (CID: 1656714) has strong binding interactions with Asp113, Asp110, Asp185 and Asp186 amino acids, all of which belong to one or the other catalytic pockets of HIV-1RT. It is expected that these interactions could be critical in the inhibitory activity of the HIV-1RT. Therefore, this study provides an evidence for consideration of (5E)-3-(2-aminoethyl)-5-benzylidene-1, 3-thiazolidine-2,4-dione as a valuable natural molecule in the treatment and prevention of HIV-associated disorders.

Structure-based evaluation of non-nucleoside inhibitors with improved potency and solubility that target HIV reverse transcriptase variants

The development of novel non-nucleoside inhibitors (NNRTIs) with activity against variants of HIV reverse transcriptase (RT) is crucial for overcoming treatment failure. The NNRTIs bind in an allosteric pocket in RT ∼10 Å away from the active site. Earlier analogues of the catechol diether compound series have picomolar activity against HIV strains with wild-type RT but lose potency against variants with single Y181C and double K103N/Y181C mutations. As guided by structure-based and computational studies, removal of the 5-Cl substitution of compound 1 on the catechol aryl ring system led to a new analogue compound 2 that maintains greater potency against Y181C and K103N/Y181C variants and better solubility (510 μg/mL). Crystal structures were determined for wild-type, Y181C, and K103N/Y181C RT in complex with both compounds 1 and 2 to understand the structural basis for these findings. Comparison of the structures reveals that the Y181C mutation destabilizes the binding mode of compound 1 and disrupts the interactions with residues in the pocket. Compound 2 maintains the same conformation in wild-type and mutant structures, in addition to several interactions with the NNRTI binding pocket. Comparison of the six crystal structures will assist in the understanding of compound binding modes and future optimization of the catechol diether series.

Altered viral fitness and drug susceptibility in HIV-1 carrying mutations that confer resistance to nonnucleoside reverse transcriptase and integrase strand transfer inhibitors

Nonnucleoside reverse transcriptase (RT) inhibitors (NNRTI) and integrase (IN) strand transfer inhibitors (INSTI) are key components of antiretroviral regimens. To explore potential interactions between NNRTI and INSTI resistance mutations, we investigated the combined effects of these mutations on drug susceptibility and fitness of human immunodeficiency virus type 1 (HIV-1). In the absence of drug, single-mutant viruses were less fit than the wild type; viruses carrying multiple mutations were less fit than single-mutant viruses. These findings were explained in part by the observation that mutant viruses carrying NNRTI plus INSTI resistance mutations had reduced amounts of virion-associated RT and/or IN protein. In the presence of efavirenz (EFV), a virus carrying RT-K103N together with IN-G140S and IN-Q148H (here termed IN-G140S/Q148H) mutations was fitter than a virus with a RT-K103N mutation alone. Similarly, in the presence of EFV, the RT-E138K plus IN-G140S/Q148H mutant virus was fitter than one with the RT-E138K mutation alone. No effect of INSTI resistance mutations on the fitness of RT-Y181C mutant viruses was observed. Conversely, RT-E138K and -Y181C mutations improved the fitness of the IN-G140S/Q148H mutant virus in the presence of raltegravir (RAL); the RT-K103N mutation had no effect. The NNRTI resistance mutations had no effect on RAL susceptibility. Likewise, the IN-G140S/Q148H mutations had no effect on EFV or RPV susceptibility. However, both the RT-K103N plus IN-G140S/Q148H and the RT-E138K plus IN-G140S/Q148H mutant viruses had significantly greater fold increases in 50% inhibitory concentration (IC50) of EFV than viruses carrying a single NNRTI mutation. Likewise, the RT-E138K plus IN-G140S/Q148H mutant virus had significantly greater fold increases in RAL IC50 than that of the IN-G140S/Q148H mutant virus. These results suggest that interactions between RT and IN mutations are important for NNRTI and INSTI resistance and viral fitness.

IMPORTANCE

Nonnucleoside reverse transcriptase inhibitors and integrase inhibitors are used to treat infection with HIV-1. Mutations that confer resistance to these drugs reduce the ability of HIV-1 to reproduce (that is, they decrease viral fitness). It is known that reverse transcriptase and integrase interact and that some mutations can disrupt their interaction, which is necessary for proper functioning of these two enzymes. To determine whether resistance mutations in these enzymes interact, we investigated their effects on drug sensitivity and viral fitness. Although individual drug resistance mutations usually reduced viral fitness, certain combinations of mutations increased fitness. When present in certain combinations, some integrase inhibitor resistance mutations increased resistance to nonnucleoside reverse transcriptase inhibitors and vice versa. Because these drugs are sometimes used together in the treatment of HIV-1 infection, these interactions could make viruses more resistant to both drugs, further limiting their clinical benefit.

Arylsulfone-based HIV-1 non-nucleoside reverse transcriptase inhibitors

HIV-1 non-nucleoside reverse transcriptase inhibitors (NNRTIs) represent one of the most significant classes of drugs for the treatment of AIDS/HIV infection. Over the past two decades several potent arylsulfone-based HIV-1 NNRTIs and related analogs have been developed. This review provides an essential overview of the structure–activity relationships of the arylsulfone-based HIV-1 NNRTIs. Furthermore, structural information useful for the design and development of new sulfur containing NNRTIs with enhanced antiretroviral activity against HIV-1 wild type and clinically relevant drug resistant HIV-1 mutant strains will be discussed.

Biophysical Insights into the Inhibitory Mechanism of Non-Nucleoside HIV-1 Reverse Transcriptase Inhibitors

HIV-1 reverse transcriptase (RT) plays a central role in HIV infection. Current United States Federal Drug Administration (USFDA)-approved antiretroviral therapies can include one of five approved non-nucleoside RT inhibitors (NNRTIs), which are potent inhibitors of RT activity. Despite their crucial clinical role in treating and preventing HIV-1 infection, their mechanism of action remains elusive. In this review, we introduce RT and highlight major advances from experimental and computational biophysical experiments toward an understanding of RT function and the inhibitory mechanism(s) of NNRTIs.

Investigating the mutation resistance of nonnucleoside inhibitors of HIV-RT using multiple microsecond atomistic simulations

Inhibiting HIV reverse transcriptase through the use of nonnucleoside reverse transcriptase inhibitors (NNRTIs) has become an essential component in drug regimens for the treatment of HIV. Older NNRTIs, such as nevirapine, are structurally rigid, exhibiting decreased inhibitory function on development of common mutations in the NNRTI-binding pocket, which is located around 10 Å from the catalytically active binding site. The newer generation of drugs, such as rilpivirine, are more flexible and resistant to binding pocket mutations but the mechanism by which they actually inhibit protein function and avoid mutations is not well-understood. To this end, we have performed 2-2.4 µs simulations with explicit solvent in an isobaric-isothermal ensemble of six different systems: apo wild-type, apo K103N/Y181C mutant, nevirapine-bound wild-type, nevirapine-bound mutant, rilpivirine-bound wild type, and rilpivirine-bound mutant. Analysis of protein conformations, principal components of motion, and mutual information between residues points to an inhibitory mechanism in which the primer grip stretches away from the catalytic triad of aspartic acids necessary for polymerization of HIV-encoding DNA, but is still unable to reveal a specific structural mechanism behind mutation resistance.

The interplay of structure and dynamics: insights from a survey of HIV-1 reverse transcriptase crystal structures

HIV-1 reverse transcriptase (RT) is a critical drug target for HIV treatment, and understanding the exact mechanisms of its function and inhibition would significantly accelerate the development of new anti-HIV drugs. It is well known that structure plays a critical role in protein function, but for RT, structural information has proven to be insufficient-despite enormous effort-to explain the mechanism of inhibition and drug resistance of non-nucleoside RT inhibitors. We hypothesize that the missing link is dynamics, information about the motions of the system. However, many of the techniques that give the best information about dynamics, such as solution nuclear magnetic resonance and molecular dynamics simulations, cannot be easily applied to a protein as large as RT. As an alternative, we combine elastic network modeling with simultaneous hierarchical clustering of structural and dynamic data. We present an extensive survey of the dynamics of RT bound to a variety of ligands and with a number of mutations, revealing a novel mechanism for drug resistance to non-nucleoside RT inhibitors. Hydrophobic core mutations restore active-state motion to multiple functionally significant regions of HIV-1 RT. This model arises out of a combination of structural and dynamic information, rather than exclusively from one or the other.

Targeted resequencing of HIV variants by microarray thermodynamics

Within a single infected individual, a virus population can have a high genomic variability. In the case of HIV, several mutations can be present even in a small genomic window of 20–30 nucleotides. For diagnostics purposes, it is often needed to resequence genomic subsets where crucial mutations are known to occur. In this article, we address this issue using DNA microarrays and inputs from hybridization thermodynamics. Hybridization signals from multiple probes are analysed, including strong signals from perfectly matching (PM) probes and a large amount of weaker cross-hybridization signals from mismatching (MM) probes. The latter are crucial in the data analysis. Seven coded clinical samples (HIV-1) are analyzed, and the microarray results are in full concordance with Sanger sequencing data. Moreover, the thermodynamic analysis of microarray signals resolves inherent ambiguities in Sanger data of mixed samples and provides additional clinically relevant information. These results show the reliability and added value of DNA microarrays for point-of-care diagnostic purposes.

Dancing with chemical formulae of antivirals: a personal account

A chemical structure is a joy forever, and this is how I perceived the chemical structures of a number of antiviral compounds with which I have been personally acquainted over the past 3 decades: (1) amino acid esters of acyclovir (i.e. valaciclovir); (2) 5-substituted 2'-deoxyuridines (i.e. brivudin); (3) 2',3'-dideoxynucleoside analogues (i.e. stavudine); (4) acyclic nucleoside phosphonates (ANPs) (i.e. cidofovir, adefovir); (5) tenofovir disoproxil fumarate (TDF) and drug combinations therewith; (6) tenofovir alafenamide (TAF, GS-7340), a new phosphonoamidate prodrug of tenofovir; (7) pro-prodrugs of PMEG (i.e. GS-9191 and GS-9219); (8) new ANPs: O-DAPy and 5-aza-C phosphonates; (9) non-nucleoside reverse transcriptase inhibitors (NNRTIs): HEPT and TIBO derivatives; and (10) bicyclam derivatives (i.e. AMD3100).

Antiviral drug development--success and failure: a personal perspective with a Japanese connection

At the 25th International Conference on Antiviral Research, I received a special recognition for my contribution to the International Society of Antiviral Research over a period of 25 years (from 1987 until 2012). This review follows the theme of my presentation at that event, which comprised 10 reminiscences, all with a Japanese connection concerning the success, or otherwise, in the clinical development of: double- and single-stranded polynucleotides; suramin, a polysulfonate; dextran sulfate, a polysulfate; brivudin; BVaraU; 2',3'-dideoxynucleoside analogues; HEPT; adefovir and tenofovir; CXCR4 antagonists; and elvitegravir.

HIV-1 reverse transcriptase and antiviral drug resistance. Part 2

Structures of RT and its complexes combined with biochemical and clinical data help in illuminating the molecular mechanisms of different drug-resistance mutations. The NRTI drugs that are used in combinations have different primary mutation sites. RT mutations that confer resistance to one drug can be hypersensitive to another RT drug. Structure of an RT-DNA-nevirapine complex revealed how NNRTI binding forbids RT from forming a polymerase competent complex. Collective knowledge about various mechanisms of drug resistance by RT has broader implications for understanding and targeting drug resistance in general. In Part 1, we discussed the role of RT in developing HIV-1 drug resistance, structural and functional states of RT, and the nucleoside/nucleotide analog (NRTI) and non-nucleoside (NNRTI) drugs used in treating HIV-1 infections. In this part, we discuss structural understanding of various mechanisms by which RT confers antiviral drug resistance.

Molecular basis of human immunodeficiency virus type 1 drug resistance: overview and recent developments

The introduction of potent combination therapies in the mid-90s had a tremendous effect on AIDS mortality. However, drug resistance has been a major factor contributing to antiretroviral therapy failure. Currently, there are 26 drugs approved for treating human immunodeficiency virus (HIV) infections, although some of them are no longer prescribed. Most of the available antiretroviral drugs target HIV genome replication (i.e. reverse transcriptase inhibitors) and viral maturation (i.e. viral protease inhibitors). Other drugs in clinical use include a viral coreceptor antagonist (maraviroc), a fusion inhibitor (enfuvirtide) and two viral integrase inhibitors (raltegravir and elvitegravir). Elvitegravir and the nonnucleoside reverse transcriptase inhibitor rilpivirine have been the most recent additions to the antiretroviral drug armamentarium. An overview of the molecular mechanisms involved in antiretroviral drug resistance and the role of drug resistance-associated mutations was previously presented (Menéndez-Arias, L., 2010. Molecular basis of human immunodeficiency virus drug resistance: an update. Antiviral Res. 85, 210-231). This article provides now an updated review that covers currently approved drugs, new experimental agents (e.g. neutralizing antibodies) and selected drugs in preclinical or early clinical development (e.g. experimental integrase inhibitors). Special attention is dedicated to recent research on resistance to reverse transcriptase and integrase inhibitors. In addition, recently discovered interactions between HIV and host proteins and novel strategies to block HIV assembly or viral entry emerge as promising alternatives for the development of effective antiretroviral treatments.

The nucleoside reverse transcriptase inhibitors, nonnucleoside reverse transcriptase inhibitors, and protease inhibitors in the treatment of HIV infections (AIDS)

The majority of the drugs currently used for the treatment of HIV infections (AIDS) belong to either of the following three classes: nucleoside reverse transcriptase inhibitors (NRTIs), nonnucleoside reverse transcriptase inhibitors (NNRTIs), and protease inhibitors (PIs). At present, there are 7 NRTIs, 5 NNRTIs, and 10 PIs approved for clinical use. They are discussed from the following viewpoints: (i) chemical formulae; (ii) mechanism of action; (iii) drug combinations; (iv) clinical aspects; (v) preexposure prophylaxis; (vi) prevention of mother-to-child transmission; (vii) their use in children; (viii) toxicity; (ix) adherence (compliance); (x) resistance; (xi) new NRTIs, NNRTIs, or PIs in (pre)clinical development; and (xii) the prospects for a "cure" of the disease.

Bifunctional inhibition of HIV-1 reverse transcriptase: a first step in designing a bifunctional triphosphate

The onset of resistance to approved anti-AIDS drugs by HIV necessitates the search for novel inhibitors of HIV-1 reverse transcriptase (RT). Developing single molecular agents concurrently occupying the nucleoside and nonnucleoside binding sites in RT is an intriguing idea but the proof of concept has so far been elusive. As a first step, we describe molecular modeling to guide focused chemical syntheses of conjugates having nucleoside (d4T) and nonnucleoside (TIBO) moieties tethered by a flexible polyethylene glycol (PEG) linker. A triphosphate of d4T-6PEG-TIBO conjugate was successfully synthesized that is recognized as a substrate by HIV-1 RT and incorporated into a double-stranded DNA.

Biochemical mechanism of HIV-1 resistance to rilpivirine

Rilpivirine (RPV) is a second generation nonnucleoside reverse transcriptase (RT) inhibitor (NNRTI) that efficiently inhibits HIV-1 resistant to first generation NNRTIs. Virological failure during therapy with RPV and emtricitabine is associated with the appearance of E138K and M184I mutations in RT. Here we investigate the biochemical mechanism of RT inhibition and resistance to RPV. We used two transient kinetics approaches (quench-flow and stopped-flow) to determine how subunit-specific mutations in RT p66 or p51 affect association and dissociation of RPV to RT as well as their impact on binding of dNTP and DNA and the catalytic incorporation of nucleotide. We compared WT with four subunit-specific RT mutants, p66(M184I)/p51(WT), p66(E138K)/p51(E138K), p66(E138K/M184I)/p51(E138K), and p66(M184I)/p51(E138K). Ile-184 in p66 (p66(184I)) decreased the catalytic efficiency of RT (k(pol)/K(d)(.dNTP)), primarily through a decrease in dNTP binding (K(d)(.dNTP)). Lys-138 either in both subunits or in p51 alone abrogated the negative effect of p66(184I) by restoring dNTP binding. Furthermore, p51(138K) reduced RPV susceptibility by altering the ratio of RPV dissociation to RPV association, resulting in a net reduction in RPV equilibrium binding affinity (K(d)(.RPV) = k(off.RPV)/k(on.RPV)). Quantum mechanics/molecular mechanics hybrid molecular modeling revealed that p51(E138K) affects access to the RPV binding site by disrupting the salt bridge between p51(E138) and p66(K101). p66(184I) caused repositioning of the Tyr-183 active site residue and decreased the efficiency of RT, whereas the addition of p51(138K) restored Tyr-183 to a WT-like conformation, thus abrogating the Ile-184-induced functional defects.

The structural biology of HIV-1: mechanistic and therapeutic insights

Three-dimensional molecular structures can provide detailed information on biological mechanisms and, for cases in which the molecular function affects human health, can significantly aid in the development of therapeutic interventions. For almost 25 years, key components of the lentivirus HIV-1, including the envelope glycoproteins, the capsid and the replication enzymes reverse transcriptase, integrase and protease, have been scrutinized to near atomic-scale resolution. Moreover, structural analyses of the interactions between viral and host cell components have yielded key insights into the mechanisms of viral entry, chromosomal integration, transcription and egress from cells. Here, we review recent advances in HIV-1 structural biology, focusing on the molecular mechanisms of viral replication and on the development of new therapeutics.

HIV-1 reverse transcriptase complex with DNA and nevirapine reveals non-nucleoside inhibition mechanism

Combinations of nucleoside and nonnucleoside inhibitors (NNRTIs) of HIV-1 reverse transcriptase (RT) are widely used in anti-AIDS therapies. Five NNRTIs including nevirapine are clinical drugs; however, the molecular mechanism of inhibition by NNRTIs is not clear. We determined the crystal structures of RT–DNA–nevirapine, RT–DNA, and RT–DNA–AZT-triphosphate complexes at 2.85, 2.70, and 2.80 Å, respectively. The RT–DNA complex in the crystal could bind nevirapine or AZT-triphosphate; however, not both. Binding of nevirapine led to opening of the NNRTI-binding pocket. The pocket formation caused shifting of the 3’-end of DNA primer by ~5.5 Å away from its polymerase active site position. Nucleic acid interactions with fingers and palm subdomains were reduced, the dNTP-binding pocket was distorted, and the thumb opened up. The structures elucidate complementary roles of nucleoside and nonnucleoside inhibitors in inhibiting RT.

CONFORMATIONAL ANALYSIS OF NEVIRAPINE IN SOLUTIONS BASED ON NMR SPECTROSCOPY AND QUANTUM CHEMICAL CALCULATIONS

The conformational analysis of HIV-1 Reverse Transcriptase Inhibitor, nevirapine, 11-cyclopropyl-5,-11dihydro-4-methyl-6H-dipyrido[3,2-b2′,3′-e][1,4]diazepin-6-one, was investigated using ab initio and density functional theory calculations. The fully optimized structures and rotational potential energies of the nitrogen and carbon bonds in the cyclopropyl ring (C15-N11-C17-C19, α) were examined in detail. Geometries obtained from all applied calculations show similarities to the complex structure with HIV-1 reverse transcriptase. To obtain more information on the structure, conformational minima of nevirapine, optimized at the B3LYP/6-31G** level, were calculated for the 1H, 13C, and 15N-NMR chemical shifts at the B3LYP/6-311++G** level using the GIAO approach in DMSO and chloroform IEFPCM solvation models. The calculated 1H, 13C-NMR chemical shifts agree well with the experimental data, which indicates that the geometry of nevirapine in solution is similar to that of the molecule in the inhibition complex. Solvation free energies (ΔG sol ) of nevirapine in DMSO and chloroform were also obtained.

Mutagenesis of human immunodeficiency virus reverse transcriptase p51 subunit defines residues contributing to vinylogous urea inhibition of ribonuclease H activity

The vinylogous urea, NSC727447, was proposed to allosterically inhibit ribonuclease H (RNase H) activity of human immunodeficiency virus type 1 reverse transcriptase (HIV-1 RT) by interacting with the thumb subdomain of its non-catalytic p51 subunit. Proximity of the p51 thumb to the p66 RNase H domain implied that inhibitor binding altered active site geometry, whereas protein footprinting suggested a contribution from α-helix I residues Cys-280 and Lys-281. To more thoroughly characterize the vinylogous urea binding site, horizontal alanine scanning mutagenesis between p51 residues Lys-275 and Thr-286 (comprising α-helix I and portions of the neighboring αH/αI and αI/αJ connecting loops) was combined with a limited vertical scan of Cys-280. A contribution from Cys-280 was strengthened by our observation that all substitutions at this position rendered selectively mutated, reconstituted p66/p51 heterodimers ∼45-fold less sensitive to inhibition. An ∼19-fold reduced IC(50) for p51 mutant T286A coupled with a 2-8-fold increased IC(50) when intervening residues were substituted supports our original proposal of p51 α-helix I as the vinylogous urea binding site. In contrast to these allosteric inhibitors, mutant enzymes retained equivalent sensitivity to the natural product α-hydroxytropolone inhibitor manicol, which x-ray crystallography has demonstrated functions by chelating divalent metal at the p66 RNase H active site. Finally, reduced DNA strand-transfer activity together with increased vinylogous urea sensitivity of p66/p51 heterodimers containing short p51 C-terminal deletions suggests an additional role for the p51 C terminus in nucleic acid binding that is compromised by inhibitor binding.

The glutamine side chain at position 91 on the β5a-β5b loop of human immunodeficiency virus type 1 reverse transcriptase is required for stabilizing the dNTP binding pocket

Earlier, we postulated that Gln91 of human immunodeficiency virus type 1 reverse transcriptase (HIV-1 RT) stabilizes the side chain of Tyr183 via hydrogen bonding interaction between O(H) of Tyr183 and CO of Q91 [Harris, D., et al. (1998) Biochemistry 37, 9630-9640]. To test this hypothesis, we generated mutant derivatives of Gln91 and analyzed their biochemical properties. The efficiency of reverse transcription was severely impaired by nonconservative substitution of Gln with Ala, while conservative substitution of Gln with Asn resulted in an approximately 70% loss of activity, a value similar to that observed with the Y183F mutation. The loss of polymerase activity from both Q91A and Q91N was significantly improved by a Met to Val substitution at position 184. Curiously, the Q91N mutant exhibited stringency in discriminating between correct and incorrect nucleotides, suggesting its possible interaction with residues influencing the flexibility of the dNTP binding pocket. In contrast, both double mutants, Q91A/M184V and Q91N/M184V, are found to be as error prone as the wild-type enzyme. We propose a model that suggests that subtle structural changes in the region due to mutation at position 91 may influence the stability of the side chain of Tyr183 in the catalytic YMDD motif of the enzyme, thus altering the active site geometry that may interfere in substrate recognition.

Multi-Site λ-dynamics for simulated Structure-Activity Relationship studies

Multi-Site λ-dynamics (MSλD) is a new free energy simulation method that is based on λ-dynamics. It has been developed to enable multiple substituents at multiple sites on a common ligand core to be modeled simultaneously and their free energies assessed. The efficacy of MSλD for estimating relative hydration free energies and relative binding affinties is demonstrated using three test systems. Model compounds representing multiple identical benzene, dihydroxybenzene and dimethoxybenzene molecules show total combined MSλD trajectory lengths of ~1.5 ns are sufficient to reliably achieve relative hydration free energy estimates within 0.2 kcal/mol and are less sensitive to the number of trajectories that are used to generate these estimates for hybrid ligands that contain up to ten substituents modeled at a single site or five substituents modeled at each of two sites. Relative hydration free energies among six benzene derivatives calculated from MSλD simulations are in very good agreement with those from alchemical free energy simulations (with average unsigned differences of 0.23 kcal/mol and R(2)=0.991) and experiment (with average unsigned errors of 1.8 kcal/mol and R(2)=0.959). Estimates of the relative binding affinities among 14 inhibitors of HIV-1 reverse transcriptase obtained from MSλD simulations are in reasonable agreement with those from traditional free energy simulations and experiment (average unsigned errors of 0.9 kcal/mol and R(2)=0.402). For the same level of accuracy and precision MSλD simulations are achieved ~20-50 times faster than traditional free energy simulations and thus with reliable force field parameters can be used effectively to screen tens to hundreds of compounds in structure-based drug design applications.

Dose-response curve slope is a missing dimension in the analysis of HIV-1 drug resistance

HIV-1 drug resistance is a major clinical problem. Resistance is evaluated using in vitro assays measuring the fold change in IC(50) caused by resistance mutations. Antiretroviral drugs are used at concentrations above IC(50), however, and inhibition at clinical concentrations can only be predicted from IC(50) if the shape of the dose-response curve is also known. Curve shape is influenced by cooperative interactions and is described mathematically by the slope parameter or Hill coefficient (m). Implicit in current analysis of resistance is the assumption that mutations shift dose-response curves to the right without affecting the slope. We show here that m is altered by resistance mutations. For reverse transcriptase and fusion inhibitors, single resistance mutations affect both slope and IC(50). For protease inhibitors, single mutations primarily affect slope. For integrase inhibitors, only IC(50) is affected. Thus, there are fundamental pharmacodynamic differences in resistance to different drug classes. Instantaneous inhibitory potential (IIP), the log inhibition of single-round infectivity at clinical concentrations, takes into account both slope and IC(50), and thus provides a direct measure of the reduction in susceptibility produced by mutations and the residual activity of drugs against resistant viruses. The standard measure, fold change in IC(50), does not correlate well with changes in IIP when mutations alter slope. These results challenge a fundamental assumption underlying current analysis of HIV-1 drug resistance and suggest that a more complete understanding of how resistance mutations reduce antiviral activity requires consideration of a previously ignored parameter, the dose-response curve slope.