Eastern extension of azoles as non-nucleoside inhibitors of HIV-1 reverse transcriptase; cyano group alternatives

Discovery of diarylpyrimidine derivatives bearing piperazine sulfonyl as potent HIV-1 nonnucleoside reverse transcriptase inhibitors

HIV-1 reverse transcriptase is one of the most attractive targets for the treatment of AIDS. However, the rapid emergence of drug-resistant strains and unsatisfactory drug-like properties seriously limit the clinical application of HIV-1 non-nucleoside reverse transcriptase inhibitors (NNRTIs). Here we show that a series of piperazine sulfonyl-bearing diarylpyrimidine-based NNRTIs were designed to improve the potency against wild-type and NNRTI-resistant strains by enhancing backbone-binding interactions. Among them, compound 18b1 demonstrates single-digit nanomolar potency against the wild-type and five mutant HIV-1 strains, which is significantly better than the approved drug etravirine. The co-crystal structure analysis and molecular dynamics simulation studies were conducted to explain the broad-spectrum inhibitory activity of 18b1 against reverse transcriptase variants. Besides, compound 18b1 demonstrates improved water solubility, cytochrome P450 liability, and other pharmacokinetic properties compared to the currently approved diarylpyrimidine (DAPY) NNRTIs. Therefore, we consider compound 18b1 a potential lead compound worthy of further study.

Unbinding Dynamics of Non-Nucleoside Inhibitors from HIV-1 Reverse Transcriptase

Non-nucleoside inhibitors of HIV-1 reverse transcriptase (NNRTIs), which bind to an allosteric site 10-15 Å from the polymerase active site, play a central role in anti-HIV chemotherapy. Though NNRTIs have been known for 30 years, the pathways by which they bind and unbind from HIV-RT have not been characterized. In crystal structures for complexes, three channels are found to extend from the NNRTI binding site to the exterior of the protein, while added mystery comes from the fact that the binding site is collapsed in the unliganded protein. To address this issue, metadynamics simulations have been performed to elucidate the unbinding of four NNRTIs from HIV-RT. A general and transferable collective variable defined by the distance between the center-of-mass (COM) of the binding pocket and COM of the ligand is used to follow the dynamics while minimizing the bias. The metadynamics also allows computation of the barriers to unbinding, which are compared with the observed potencies of the compounds in an antiviral assay.

Computer-aided discovery of anti-HIV agents

A review is provided on efforts in our laboratory over the last decade to discover anti-HIV agents. The work has focused on computer-aided design and synthesis of non-nucleoside inhibitors of HIV-1 reverse transcriptase (NNRTIs) with collaborative efforts on biological assaying and protein crystallography. Numerous design issues were successfully addressed including the need for potency against a wide range of viral variants, good aqueous solubility, and avoidance of electrophilic substructures. Computational methods including docking, de novo design, and free-energy perturbation (FEP) calculations made essential contributions. The result is novel NNRTIs with picomolar and low-nanomolar activities against wild-type HIV-1 and key variants that also show much improved solubility and lower cytotoxicity than recently approved drugs in the class.

Computational drug design strategies applied to the modelling of human immunodeficiency virus-1 reverse transcriptase inhibitors

Reverse transcriptase (RT) is a multifunctional enzyme in the human immunodeficiency virus (HIV)-1 life cycle and represents a primary target for drug discovery efforts against HIV-1 infection. Two classes of RT inhibitors, the nucleoside RT inhibitors (NRTIs) and the nonnucleoside transcriptase inhibitors are prominently used in the highly active antiretroviral therapy in combination with other anti-HIV drugs. However, the rapid emergence of drug-resistant viral strains has limited the successful rate of the anti-HIV agents. Computational methods are a significant part of the drug design process and indispensable to study drug resistance. In this review, recent advances in computer-aided drug design for the rational design of new compounds against HIV-1 RT using methods such as molecular docking, molecular dynamics, free energy calculations, quantitative structure-activity relationships, pharmacophore modelling and absorption, distribution, metabolism, excretion and toxicity prediction are discussed. Successful applications of these methodologies are also highlighted.

Structural optimization of pyridine-type DAPY derivatives to exploit the tolerant regions of the NNRTI binding pocket

Based on the crystallographic studies of diarylpyrimidines (DAPYs), we embarked on incorporating the hydrophilic piperidyl or morpholinyl group into the known DAPY derivatives bearing the pyridine moiety as a core structure, with the double aim to exploit additional interactions with the HIV-1 NNRTI binding pocket (NNIBP), as well as to improve the compound solubility. The antiviral evaluation result show that the most potent compounds I-8b2, I-8b3, I-8b4 and I-8c3 exhibited anti-HIV-1 (IIIB) strain activity ranging from 7.4 nM to 9.4 nM (SI = 168-1283), superior to FDA-approved drugs of nevirapine (NVP), lamivudine (3TC) and delavirdine (DLV), and comparable to etravirine (ETV), zidovudine (AZT) and efavirenz (EFV). Additionally, compounds I-8c2 and I-8c3 showed moderate activity against NNRTI resistant strains baring mutations K103N and Y181C with EC50 values of 6.2 μM and 6.8 μM, respectively. Preliminary structure-activity relationships (SARs), reverse transcriptase inhibition efficacy and molecular modeling of selected compounds are also presented. These outcomes support our design hypothesis and demonstrate that the piperidyl group modified pyridine-typed DAPY derivatives are highly potent NNRTIs with improved water solubility.

Searching for novel scaffold of triazole non-nucleoside inhibitors of HIV-1 reverse transcriptase

Azoles are a promising class of the new generation of HIV-1 non-nucleoside reverse transcriptase inhibitors (NNRTIs). From thousands of reported compounds, many possess the same basic structure of an aryl substituted azole ring linked by a thioglycolamide chain with another aromatic ring. In order to find novel extensions for this basic scaffold, we explored the 5-position substitution pattern of triazole NNRTIs using molecular docking followed by the synthesis of selected compounds. We found that heterocyclic substituents in the 5-position of the triazole ring are detrimental to the inhibitory activity of compounds with four-membered thioglycolamide linker and this substitution seems to be viable only for compounds with shorter two-membered linker. Promising compound, N-(4-carboxy-2-chlorophenyl)-2-((4-benzyl-5-methyl-4H-1,2,4-triazol-3-yl)sulfanyl)acetamide, with potent inhibitory activity and acceptable aqueous solubility has been identified in this study that could serve as lead scaffold for the development of novel water-soluble salts of triazole NNRTIs.

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

Enhanced Monte Carlo Sampling through Replica Exchange with Solute Tempering

With a view to improving the consistency of free energy perturbation calculations in Monte Carlo simulations of protein–ligand complexes, we have implemented the replica exchange with solute tempering (REST) method in the MCPRO software. By augmenting the standard REST approach with regular attempted jumps in selected dihedral angles, our combined method facilitates sampling of ligand binding modes that are separated by high free energy barriers and ensures that computed free energy changes are considerably less dependent on the starting conditions and the chosen mutation pathway than those calculated with standard Monte Carlo sampling. We have applied the enhanced sampling method to the calculation of the activities of seven non-nucleoside inhibitors of HIV-1 reverse transcriptase, and its Tyr181Cys variant, and have shown that a range of binding orientations is possible depending on the nature of the ligand and the presence of mutations at the binding site.

Design, synthesis, and biological evaluation of 1-phenylpyrazolo[3,4-e]pyrrolo[3,4-g]indolizine-4,6(1H,5H)-diones as new glycogen synthase kinase-3β inhibitors

Compound 5 was selected from our in-house library as a suitable starting point for the rational design of new GSK-3β inhibitors. MC/FEP calculations of 5 led to the identification of a structural class of new GSK-3β inhibitors. Compound 18 inhibited GSK-3β with an IC50 of 0.24 μM and inhibited tau phosphorylation in a cell-based assay. It proved to be a selective inhibitor of GSK-3 against a panel of 17 kinases and showed >10-fold selectivity against CDK2. Calculated physicochemical properties and Volsurf predictions suggested that compound 18 has the potential to diffuse passively across the blood-brain barrier.

Extension into the entrance channel of HIV-1 reverse transcriptase--crystallography and enhanced solubility

Non-nucleoside inhibitors of HIV-1 reverse transcriptase (HIV-RT) are reported that feature extension into the entrance channel near Glu138. Complexes of the parent anilinylpyrimidine 1 and the morpholinoethoxy analog 2j with HIV-RT have received crystallographic characterization confirming the designs. Measurement of aqueous solubilities of 2j, 2k, the parent triazene 2a, and other NNRTIs demonstrate profound benefits for addition of the morpholinyl substituent.

Substituted oxadiazoles: a patent review (2010 - 2012)

INTRODUCTION

The oxadiazoles represent a class of five-membered heterocyclic compounds which are of considerable interest in different areas of medicinal chemistry and drug discovery. Oxadiazoles can exist in different regioisomeric forms and employ in various agents with a broad range of biological activities. This review covers the work reported on various biological activities of oxadiazole derivatives from 2010 to 2012.

AREAS COVERED

Oxadiazole derivatives attract great attention due to their different kinds of pharmaceutical activities including antiviral, antimicrobial, anticancer, anticonvulsant, antidiabetic and anti-inflammatory activity. This paper provides a general review of oxadiazole derivatives published in international journals and patented between 2010 and 2012.

EXPERT OPINION

Oxadiazoles have been used frequently in drug-like molecules as bioisosteres for ester and amide functionalities and displayed numerous prominent pharmacological effects. The broad pharmacological profile of oxadiazole derivatives has attracted the attention of many researchers to explore this scaffold to its multiple potential against several activities. Therefore, oxadiazole motif is likely to be present in other therapeutic molecules in the future.

Optimization of benzyloxazoles as non-nucleoside inhibitors of HIV-1 reverse transcriptase to enhance Y181C potency

Design of non-nucleoside inhibitors of HIV-1 reverse transcriptase with improved activity towards Tyr181Cys containing variants was pursued with the assistance of free energy perturbation (FEP) calculations. Optimization of the 4-R substituent in 1 led to ethyl and isopropyl analogs 1e and 1f with 1-7 nM potency towards both the wild-type virus and a Tyr181C variant.

Identification of HIV inhibitors guided by free energy perturbation calculations

Free energy perturbation (FEP) theory coupled to molecular dynamics (MD) or Monte Carlo (MC) statistical mechanics offers a theoretically precise method for determining the free energy differences of related biological inhibitors. Traditionally requiring extensive computational resources and expertise, it is only recently that its impact is being felt in drug discovery. A review of computer-aided anti-HIV efforts employing FEP calculations is provided here that describes early and recent successes in the design of human immunodeficiency virus type 1 (HIV-1) protease and non-nucleoside reverse transcriptase inhibitors. In addition, our ongoing work developing and optimizing leads for small molecule inhibitors of cyclophilin A (CypA) is highlighted as an update on the current capabilities of the field. CypA has been shown to aid HIV-1 replication by catalyzing the cis/trans isomerization of a conserved Gly-Pro motif in the Nterminal domain of HIV-1 capsid (CA) protein. In the absence of a functional CypA, e.g., by the addition of an inhibitor such as cyclosporine A (CsA), HIV-1 has reduced infectivity. Our simulations of acylurea-based and 1-indanylketone-based CypA inhibitors have determined that their nanomolar and micromolar binding affinities, respectively, are tied to their ability to stabilize Arg55 and Asn102. A structurally novel 1-(2,6-dichlorobenzamido) indole core was proposed to maximize these interactions. FEP-guided optimization, experimental synthesis, and biological testing of lead compounds for toxicity and inhibition of wild-type HIV-1 and CA mutants have demonstrated a dose-dependent inhibition of HIV-1 infection in two cell lines. While the inhibition is modest compared to CsA, the results are encouraging.

Methyl effects on protein-ligand binding

The effects of addition of a methyl group to a lead compound on biological activity are examined. A literature analysis of >2000 cases reveals that an activity boost of a factor of 10 or more is found with an 8% frequency, and a 100-fold boost is a 1 in 200 event. Four cases in the latter category are analyzed in depth to elucidate any unusual aspects of the protein-ligand binding, distribution of water molecules, and changes in conformational energetics. The analyses include Monte Carlo/free-energy perturbation (MC/FEP) calculations for methyl replacements in inhibitor series for p38α MAP kinase, ACK1, PTP1B, and thrombin. Methyl substitutions ortho to an aryl ring can be particularly effective at improving activity by inducing a propitious conformational change. The greatest improvements in activity arise from coupling the conformational gain with the burial of the methyl group in a hydrophobic region of the protein.

Discovery of dimeric inhibitors by extension into the entrance channel of HIV-1 reverse transcriptase

Design of non-nucleoside inhibitors of HIV-1 reverse transcriptase is being pursued with computational guidance. Extension of azine-containing inhibitors into the entrance channel between Lys103 and Glu138 has led to the discovery of potent and structurally novel derivatives including dimeric inhibitors in an NNRTI-linker-NNRTI motif.

Computationally-guided optimization of a docking hit to yield catechol diethers as potent anti-HIV agents

A 5-μM docking hit has been optimized to an extraordinarily potent (55 pM) non-nucleoside inhibitor of HIV reverse transcriptase. Use of free energy perturbation (FEP) calculations to predict relative free energies of binding aided the optimizations by identifying optimal substitution patterns for phenyl rings and a linker. The most potent resultant catechol diethers feature terminal uracil and cyanovinylphenyl groups. A halogen bond with Pro95 likely contributes to the extreme potency of compound 42. In addition, several examples are provided illustrating failures of attempted grafting of a substructure from a very active compound onto a seemingly related scaffold to improve its activity.

Efficient discovery of potent anti-HIV agents targeting the Tyr181Cys variant of HIV reverse transcriptase

Non-nucleoside reverse transcriptase inhibitors (NNRTIs) that interfere with the replication of human immunodeficiency virus (HIV) are being pursued with guidance from molecular modeling including free-energy perturbation (FEP) calculations for protein-inhibitor binding affinities. The previously reported pyrimidinylphenylamine 1 and its chloro analogue 2 are potent anti-HIV agents; they inhibit replication of wild-type HIV-1 in infected human T-cells with EC(50) values of 2 and 10 nM, respectively. However, they show no activity against viral strains containing the Tyr181Cys (Y181C) mutation in HIV-RT. Modeling indicates that the problem is likely associated with extensive interaction between the dimethylallyloxy substituent and Tyr181. As an alternative, a phenoxy group is computed to be oriented in a manner diminishing the contact with Tyr181. However, this replacement leads to a roughly 1000-fold loss of activity for 3 (2.5 μM). The present report details the efficient, computationally driven evolution of 3 to novel NNRTIs with sub-10 nM potency toward both wild-type HIV-1 and Y181C-containing variants. The critical contributors were FEP substituent scans for the phenoxy and pyrimidine rings and recognition of potential benefits of addition of a cyanovinyl group to the phenoxy ring.

Imidazo[1,2-a]pyridin-3-amines as potential HIV-1 non-nucleoside reverse transcriptase inhibitors

During random screening of a small in-house library of compounds, certain substituted imidazo[1,2-a]pyridines were found to be weak allosteric inhibitors of HIV-1 reverse transcriptase (RT). A library of these compounds was prepared using the Groebke reaction and a subset of compounds prepared from 2-chlorobenzaldehyde, cyclohexyl isocyanide and a 6-substituted 2-aminopyridine showed good inhibitory activity in enzymatic (RT) and HIV anti-infectivity MAGI whole cell assays. The compound showing the best anti-HIV-1 IIIB whole cell activity (MAGI IC(50)=0.18 μM, IC(90)=1.06 μM), along with a good selectivity index (>800), was 2-(2-chlorophenyl)-3-(cyclohexylamino)imidazo[1,2-a]pyridine-5-carbonitrile 38.

HIV-1 NNRTIs: structural diversity, pharmacophore similarity, and implications for drug design

Nonnucleoside reverse transcriptase inhibitors (NNRTIs) nowadays represent very potent and most promising anti-AIDS agents that specifically target the HIV-1 reverse transcriptase (RT). However, the effectiveness of NNRTI drugs can be hampered by rapid emergence of drug-resistant viruses and severe side effects upon long-term use. Therefore, there is an urgent need to develop novel, highly potent NNRTIs with broad spectrum antiviral activity and improved pharmacokinetic properties, and more efficient strategies that facilitate and shorten the drug discovery process would be extremely beneficial. Fortunately, the structural diversity of NNRTIs provided a wide space for novel lead discovery, and the pharmacophore similarity of NNRTIs gave valuable hints for lead discovery and optimization. More importantly, with the continued efforts in the development of computational tools and increased crystallographic information on RT/NNRTI complexes, structure-based approaches using a combination of traditional medicinal chemistry, structural biology, and computational chemistry are being used increasingly in the design of NNRTIs. First, this review covers two decades of research and development for various NNRTI families based on their chemical scaffolds, and then describes the structural similarity of NNRTIs. We have attempted to assemble a comprehensive overview of the general approaches in NNRTI lead discovery and optimization reported in the literature during the last decade. The successful applications of medicinal chemistry strategies, crystallography, and computational tools for designing novel NNRTIs are highlighted. Future directions for research are also outlined.

Alchemical free energy methods for drug discovery: progress and challenges

Improved rational drug design methods are needed to lower the cost and increase the success rate of drug discovery and development. Alchemical binding free energy calculations, one potential tool for rational design, have progressed rapidly over the past decade, but still fall short of providing robust tools for pharmaceutical engineering. Recent studies, especially on model receptor systems, have clarified many of the challenges that must be overcome for robust predictions of binding affinity to be useful in rational design. In this review, inspired by a recent joint academic/industry meeting organized by the authors, we discuss these challenges and suggest a number of promising approaches for overcoming them.

Effects of Water Placement on Predictions of Binding Affinities for p38α MAP Kinase Inhibitors

Monte Carlo free energy perturbation (MC/FEP) calculations have been applied to compute the relative binding affinities of 17 congeneric pyridazo-pyrimidinone inhibitors of the protein p38α MAP kinase. Overall correlation with experiment was found to be modest when the complexes were hydrated using a traditional procedure with a stored solvent box. Significant improvements in accuracy were obtained when the MC/FEP calculations were repeated using initial solvent distributions optimized by the water placement algorithm JAWS. The results underscore the importance of accurate placement of water molecules in a ligand binding site for the reliable prediction of relative free energies of binding.