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

TAR RNA selective targeting ruthenium(II) complexes as HIV-1 reverse transcriptase inhibitors: On exploring structure-activity relationships of multiple positions

HIV-1 reverse transcriptase (RT) inhibitors play a crucial role in the treatment of HIV by preventing the activity of the enzyme responsible for the replication of the virus. The HIV-1 Tat protein binds to transactivation response (TAR) RNA and recruits host factors to stimulate HIV-1 transcription. We have created a small library consisting of 4 × 6 polypyridyl Ru(II) complexes that selectively bind to TAR RNA, with targeting groups specific to HIV-1 TAR RNA. The molecule design was conducted by introducing hydroxyl or methoxy groups into an established potent TAR binder. The potential TAR binding ability was analysis from nature charge population and electrostatic potential by quantum chemistry calculations. Key modifications were found to be R1 and R3 groups. The most potent and selective TAR RNA binder was a3 with R1 = OH, R2 = H and R3 = Me. Through molecular recognition of hydrogen bonds and electrostatic attraction, they were able to firmly and selectively bind HIV-1 TAR RNA. Furthermore, they efficiently obstructed the contact between TAR RNA and Tat protein, and inhibited the reverse transcription activity of HIV-1 RT. The polypyridyl Ru(II) complexes were chemical and photo-stable, and sensitive and selective spectroscopic responses to TAR RNA. They exhibited little toxicity towards normal cells. Hence, this study might offer significant drug design approaches for researching AIDS and other illnesses associated with RT, including HCV, EBOV, and SARS-CoV-2. Moreover, it could contribute to fundamental research on the interactions of inorganic transition metal complexes with biomolecules.

In silico λ-dynamics predicts protein binding specificities to modified RNAs

RNA modifications shape gene expression through a smorgasbord of chemical changes to canonical RNA bases. Although numbering in the hundreds, only a few RNA modifications are well characterized, in part due to the absence of methods to identify modification sites. Antibodies remain a common tool to identify modified RNA and infer modification sites through straightforward applications. However, specificity issues can result in off-target binding and confound conclusions. This work utilizes in silico λ-dynamics to efficiently estimate binding free energy differences of modification-targeting antibodies between a variety of naturally occurring RNA modifications. Crystal structures of inosine and N6-methyladenosine (m6A) targeting antibodies bound to their modified ribonucleosides were determined and served as structural starting points. λ-Dynamics was utilized to predict RNA modifications that permit or inhibit binding to these antibodies. In vitro RNA-antibody binding assays supported the accuracy of these in silico results. High agreement between experimental and computed binding propensities demonstrated that λ-dynamics can serve as a predictive screen for antibody specificity against libraries of RNA modifications. More importantly, this strategy is an innovative way to elucidate how hundreds of known RNA modifications interact with biological molecules without the limitations imposed by in vitro or in vivo methodologies.

Covalent and noncovalent strategies for targeting Lys102 in HIV-1 reverse transcriptase

Reverse transcriptase (RT) is one of three key proteins responsible for the replication cycle of HIV-1 in the host. Several classes of inhibitors have been developed to target the enzyme, with non-nucleoside reverse transcriptase inhibitors forming first-line treatment. Previously, covalent RT inhibitors have been identified and found to bind irreversibly to commonly mutated residues such as Y181C. In this work we aim to circumvent the issue of NNRTI resistance through targeting K102, which has not yet been identified to confer drug resistance. As reported here, 34 compounds were synthesized and characterized biochemically and structurally with wild-type (WT) HIV-1 RT. Two of these inhibitors demonstrate covalent inhibition as evidenced by protein crystallography, enzyme kinetics, mass spectrometry, and antiviral potency in HIV-1 infected human T-cell assays.

Exploring novel HIV-1 reverse transcriptase inhibitors with drug-resistant mutants: A double mutant surprise

HIV-1 reverse transcriptase (RT) remains a key target for HIV drug development. As successful management of the disease requires lifelong treatment, the emergence of resistance mutations is inevitable, making development of new RT inhibitors, which remain effective against resistant variants crucial. To this end, previous computationally guided drug design efforts have resulted in catechol diether compounds, which inhibit wildtype RT with picomolar affinities and appear to be promising preclinical candidates. To confirm that these compounds remain potent against Y181C, a widespread mutation conferring resistance to first generation inhibitors, they were screened against the HIV-1 N119 clinical isolate, reported as a Y181C single mutant. In comparison to a molecular clone with the same mutation, N119 appears less susceptible to inhibition by our preclinical candidate compounds. A more detailed sequencing effort determined that N119 was misidentified and carries V106A in combination with Y181C. While both indolizine and naphthalene substituted catechol diethers are potent against the classical Y181C single mutant, the addition of V106A confers more resistance against the indolizine derivatives than the naphthalene derivatives. Crystal structures presented in this study highlight key features of the naphthyl group, which allow these compounds to remain potent in the double mutant, including stronger interactions with F227 and less reliance on V106 for stabilization of the ethoxy-uracil ring, which makes critical hydrogen bonds with other residues in the binding pocket.

Aromatic interactions with heterocycles in water

Conformationally well-defined supramolecular complexes that can be studied in different solvents provide a platform for separating and quantifying free energy contributions due to functional group interactions and desolvation. Here 1:1 complexes formed between four different calix[4]pyrrole receptors and eleven different pyridine N-oxide guests have been used to dissect the factors that govern aromatic interactions with heterocycles in water and in chloroform solution. 1H NMR spectroscopy shows that the three-dimensional structures of the complexes are fixed by four H-bonding interactions between the pyrrole donors at the bottom of the receptor and the N-oxide acceptor on the guest, locking the geometrical arrangement of interacting functional groups in the binding pocket at the other end of the receptor. An aromatic heterocycle on the guest makes two stacking interactions and two edge-to-face interactions with the side walls of the receptor. Chemical double mutant cycles were used to measure the free energy contribution of these four aromatic interactions to the overall stability of the complex. In chloroform, the aromatic interactions measured with pyridine, pyrimidine, furan, thiophene and thiazole are similar to the interactions with a phenyl group, but the effect of introducing a heteroatom depends on where it sits with respect to the aromatic side-walls of the cavity. A nitrogen lone pair directed into a π-face of the side-walls of the binding site leads to repulsive interactions of up to 8 kJ mol-1. In water, the heterocycle aromatic interactions are all significantly more favourable (by up to 12 kJ mol-1). For the non-polar heterocycles, furan and thiophene, the increase in interaction energy correlates directly with hydrophobicity, as measured by the free energy of transfer of the heterocycle from n-hexadecane into water (ΔG°(water-hex)). For the heterocycles with polar nitrogen H-bond acceptors, water can access cracks in the walls of the receptor binding site to solvate the edges of the heterocycles without significantly affecting the geometry of the aromatic interactions, and these nitrogen-water H-bonds stabilise the complexes by about 15 kJ mol-1. The results highlight the complexity of the solvation processes that govern molecular recognition in water.

An open-source molecular builder and free energy preparation workflow

Automated free energy calculations for the prediction of binding free energies of congeneric series of ligands to a protein target are growing in popularity, but building reliable initial binding poses for the ligands is challenging. Here, we introduce the open-source FEgrow workflow for building user-defined congeneric series of ligands in protein binding pockets for input to free energy calculations. For a given ligand core and receptor structure, FEgrow enumerates and optimises the bioactive conformations of the grown functional group(s), making use of hybrid machine learning/molecular mechanics potential energy functions where possible. Low energy structures are optionally scored using the gnina convolutional neural network scoring function, and output for more rigorous protein-ligand binding free energy predictions. We illustrate use of the workflow by building and scoring binding poses for ten congeneric series of ligands bound to targets from a standard, high quality dataset of protein-ligand complexes. Furthermore, we build a set of 13 inhibitors of the SARS-CoV-2 main protease from the literature, and use free energy calculations to retrospectively compute their relative binding free energies. FEgrow is freely available at https://github.com/cole-group/FEgrow, along with a tutorial.

Scorpion-Shaped Zinc Porphyrins as Tetrafunctional TAR RNA Predators and HIV-1 Reverse Transcriptase Inhibitors

HIV-1 reverse transcriptase (RT) inhibitors are fundamental to the discovery and development of anti-HIV drugs. Their main target is RT, and only a tiny number of them can bind to viral RNA. In this paper, five new Zn(II) porphyrin compounds were developed with different characters. ZnTPP4 has both the appearance and the functions of a scorpion with a rigid tail and stinger to selectively hunt HIV-1 TAR RNA based on the molecular recognition of hydrogen bonds, a fierce chelicera to bite RNA by metal coordination, mighty pedipalps to grasp the bound RNA by supramolecular inclusion, and a broad body maintaining the configuration of each functional area so that they can cooperate with each other and providing accommodation space for the bound RNA. This tetrafunctional Zn(II) porphyrin is relatively nontoxic to normal cells and can produce sensitive responses for RNA. Moreover, this work offers practical construction methodologies for medication of AIDS and other diseases closely related to RT like EBOV and SARS-CoV-2.

Addressing Intersite Coupling Unlocks Large Combinatorial Chemical Spaces for Alchemical Free Energy Methods

Alchemical free energy methods are playing a growing role in molecular design, both for computer-aided drug design of small molecules and for computational protein design. Multisite λ dynamics (MSλD) is a uniquely scalable alchemical free energy method that enables more efficient exploration of combinatorial alchemical spaces encountered in molecular design, but simulations have typically been limited to a few hundred ligands or sequences. Here, we focus on coupling between sites to enable scaling to larger alchemical spaces. We first discuss updates to the biasing potentials that facilitate MSλD sampling to include coupling terms and show that this can provide more thorough sampling of alchemical states. We then harness coupling between sites by developing a new free energy estimator based on the Potts models underlying direct coupling analysis, a method for predicting contacts from sequence coevolution, and find it yields more accurate free energies than previous estimators. The sampling requirements of the Potts model estimator scale with the square of the number of sites, a substantial improvement over the exponential scaling of the standard estimator. This opens up exploration of much larger alchemical spaces with MSλD for molecular design.

Contemporary Medicinal Chemistry Strategies for the Discovery and Development of Novel HIV-1 Non-nucleoside Reverse Transcriptase Inhibitors

Currently, HIV-1 non-nucleoside reverse transcriptase inhibitors (NNRTIs) are a major component of the highly active anti-retroviral therapy (HAART) regimen. However, the occurrence of drug-resistant strains and adverse reactions after long-term usage have inevitably compromised the clinical application of NNRTIs. Therefore, the development of novel inhibitors with distinct anti-resistance profiles and better pharmacological properties is still an enormous challenge. Herein, we summarize state-of-the-art medicinal chemistry strategies for the discovery of potent NNRTIs, such as structure-based design strategies, contemporary computer-aided drug design, covalent-binding strategies, and the application of multi-target-directed ligands. The strategies described here will facilitate the identification of promising HIV-1 NNRTIs.

Structural Studies and Structure Activity Relationships for Novel Computationally Designed Non-nucleoside Inhibitors and Their Interactions With HIV-1 Reverse Transcriptase

Reverse transcriptase (RT) from the human immunodeficiency virus continues to be an attractive drug target for antiretroviral therapy. June 2022 will commemorate the 30th anniversary of the first Human Immunodeficiency Virus (HIV) RT crystal structure complex that was solved with non-nucleoside reverse transcriptase inhibitor nevirapine. The release of this structure opened opportunities for designing many families of non-nucleoside reverse transcriptase inhibitors (NNRTIs). In paying tribute to the first RT-nevirapine structure, we have developed several compound classes targeting the non-nucleoside inhibitor binding pocket of HIV RT. Extensive analysis of crystal structures of RT in complex with the compounds informed iterations of structure-based drug design. Structures of seven additional complexes were determined and analyzed to summarize key interactions with residues in the non-nucleoside inhibitor binding pocket (NNIBP) of RT. Additional insights comparing structures with antiviral data and results from molecular dynamics simulations elucidate key interactions and dynamics between the nucleotide and non-nucleoside binding sites.

Long-acting and extended-release implant and nanoformulations with a synergistic antiretroviral two-drug combination controls HIV-1 infection in a humanized mouse model

The HIV pandemic has affected over 38 million people worldwide with close to 26 million currently accessing antiretroviral therapy (ART). A major challenge in the long-term treatment of HIV-1 infection is nonadherence to ART. Long-acting antiretroviral (LA-ARV) formulations, that reduce dosing frequency to less than once a day, are an urgent need that could tackle the adherence issue. Here, we have developed two LA-ART interventions, one an injectable nanoformulation, and the other, a removable implant, for the delivery of a synergistic two-drug ARV combination comprising a pre-clinical nonnucleoside reverse transcriptase inhibitor (NNRTI), Compound I, and the nucleoside reverse transcriptase inhibitor (NRTI), 4'-ethynyl-2-fluoro-2'-deoxyadenosine. The nanoformulation is poly(lactide-co-glycolide)-based and the implant is a copolymer of ω-pentadecalactone and p-dioxanone, poly(PDL-co-DO), a novel class of biocompatible, biodegradable materials. Both the interventions, packaged independently with each ARV, released sustained levels of the drugs, maintaining plasma therapeutic indices for over a month, and suppressed viremia in HIV-1-infected humanized mice for up to 42 days with maintenance of CD4+ T cells. These data suggest promise in the use of these new drugs as LA-ART formulations in subdermal implant and injectable mode.

Synthesis and biological evaluation of indole-based peptidomimetics as antibacterial agents against Gram-positive bacteria

The emergence of bacterial multidrug resistance and the lack of new antimicrobial agents urgently demand the discovery and development of novel antibacterials that avoid bacterial resistance. Antimicrobial peptidomimetics represent a promising approach for overcoming antibiotic resistance. Herein we report the synthesis and evaluation of indole-based amphiphilic antimicrobial peptidomimetics, bearing hydrophobic side chains and hydrophilic cationic moieties. Among these derivatives, compound 28 demonstrated potent antimicrobial activity against Gram-positive bacteria, low hemolytic activity and low toxicity towards mammalian cells, as well as good stability in salt conditions. Moreover, compound 28 showed the rapid killing of bacteria via membrane-targeting action without developing bacterial resistance. More importantly, compound 28 displayed high antimicrobial potency against Gram-positive bacteria in a murine model of bacterial keratitis, and was found to be more efficient than vancomycin. Thus, compound 28 had great potential as a promising lead compound for the treatment of Gram-positive bacterial infection.

Design, synthesis, and antiviral evaluation of novel piperidine-substituted arylpyrimidines as HIV-1 NNRTIs by exploring the hydrophobic channel of NNIBP

Herein, alkenylpiperidine and alkynylpiperidine moieties were introduced into the left wing of DAPYs (diarylpyrimidines) to explore the new site of the NNIBP (non-nucleoside inhibitor binding pocket) protein-solvent interface region via the structure-based drug design strategy. All the synthesized compounds displayed nanomolar to submicromolar activity against WT (wild-type) HIV-1. Among all, compound FT1 (EC₅₀ = 19 nM) was found to be the most active molecule, which is better than NVP (EC₅₀ = 0.10 μM). In addition, most of the compounds displayed micromolar activity against K103N and E138K mutant strains, while FT1 (EC_50(K103N) = 50 nM, EC_50(E138K) = 0.19 µM) still has the most effective activity. The molecular dynamics simulation studies revealed that the presence of pyridine moiety of FT1 was essential and played a significant role in its binding with RT (reverse transcriptase).

Generalizing the Discrete Gibbs Sampler-Based λ-Dynamics Approach for Multisite Sampling of Many Ligands

In this work, the discrete λ variant of the Gibbs sampler-based λ-dynamics (d-GSλD) method is developed to enable multiple functional group perturbations to be investigated at one or more sites of substitution off a common ligand core. The theoretical framework and special considerations for constructing discrete λ states for multisite d-GSλD are presented. The precision and accuracy of the d-GSλD method is evaluated with three test cases of increasing complexity. Specifically, methyl → methyl symmetric perturbations in water, 1,4-benzene hydration free energies and protein-ligand binding affinities for an example HIV-1 reverse transcriptase inhibitor series are computed with d-GSλD. Complementary MSλD calculations were also performed to compare with d-GSλD's performance. Excellent agreement between d-GSλD and MSλD is observed, with mean unsigned errors of 0.12 and 0.22 kcal/mol for computed hydration and binding free energy test cases, respectively. Good agreement with experiment is also observed, with errors of 0.5-0.7 kcal/mol. These findings support the applicability of the d-GSλD free energy method for a variety of molecular design problems, including structure-based drug design. Finally, a discussion of d-GSλD versus MSλD approaches is presented to compare and contrast features of both methods.

Exploiting the hydrophobic channel of the NNIBP: Discovery of novel diarylpyrimidines as HIV-1 NNRTIs against wild-type and K103N mutant viruses

To further explore the chemical space surrounding the "hydrophobic channel" of the NNRTI binding pocket (NNIBP), a new series of diarylpyrimidines (DAPYs) were designed and synthesized as potent HIV-1 non-nucleoside RT inhibitors (NNRTIs). The target compounds were evaluated for anti-HIV potency in MT-4 cells. Most of the synthesized DAPYs exhibited moderate to excellent activity against the HIV-1 wild-type (WT) strain with EC₅₀ values ranging from 16 nM to 0.722 µM. Interestingly, few compounds displayed remarkable activity in inhibiting K103N mutant virus with EC₅₀ values ranging from 39 nM to 1.708 µM. Notably, FS2 (EC_50(IIIB) = 16 nM, EC_50(K103N) = 39 nM, SI = 294) was identified as the most significant compound, which was considerably more potent than nevirapine, lamivudine, and comparable to zidovudine. Additionally, the HIV-1 RT inhibition assay confirmed their binding target. Preliminary structure-activity relationships (SARs) and molecular modeling studies were also performed, providing significant suggestions for further optimization.

Implementation of the QUBE Force Field in SOMD for High-Throughput Alchemical Free-Energy Calculations

The quantum mechanical bespoke (QUBE) force-field approach has been developed to facilitate the automated derivation of potential energy function parameters for modeling protein-ligand binding. To date, the approach has been validated in the context of Monte Carlo simulations of protein-ligand complexes. We describe here the implementation of the QUBE force field in the alchemical free-energy calculation molecular dynamics simulation package SOMD. The implementation is validated by demonstrating the reproducibility of absolute hydration free energies computed with the QUBE force field across the SOMD and GROMACS software packages. We further demonstrate, by way of a case study involving two series of non-nucleoside inhibitors of HIV-1 reverse transcriptase, that the availability of QUBE in a modern simulation package that makes efficient use of graphics processing unit acceleration will facilitate high-throughput alchemical free-energy calculations.

Drug Design Strategies to Avoid Resistance in Direct-Acting Antivirals and Beyond

Drug resistance is prevalent across many diseases, rendering therapies ineffective with severe financial and health consequences. Rather than accepting resistance after the fact, proactive strategies need to be incorporated into the drug design and development process to minimize the impact of drug resistance. These strategies can be derived from our experience with viral disease targets where multiple generations of drugs had to be developed to combat resistance and avoid antiviral failure. Significant efforts including experimental and computational structural biology, medicinal chemistry, and machine learning have focused on understanding the mechanisms and structural basis of resistance against direct-acting antiviral (DAA) drugs. Integrated methods show promise for being predictive of resistance and potency. In this review, we give an overview of this research for human immunodeficiency virus type 1, hepatitis C virus, and influenza virus and the lessons learned from resistance mechanisms of DAAs. These lessons translate into rational strategies to avoid resistance in drug design, which can be generalized and applied beyond viral targets. While resistance may not be completely avoidable, rational drug design can and should incorporate strategies at the outset of drug development to decrease the prevalence of drug resistance.

Covalent Inhibition of Wild-Type HIV-1 Reverse Transcriptase Using a Fluorosulfate Warhead

Covalent inhibitors of wild-type HIV-1 reverse transcriptase (CRTIs) are reported. Three compounds derived from catechol diether non-nucleoside inhibitors (NNRTIs) with addition of a fluorosulfate warhead are demonstrated to covalently modify Tyr181 of HIV-RT. X-ray crystal structures for complexes of the CRTIs with the enzyme are provided, which fully demonstrate the covalent attachment, and confirmation is provided by appropriate mass shifts in ESI-TOF mass spectra. The three CRTIs and six noncovalent analogues are found to be potent inhibitors with both IC₅₀ values for in vitro inhibition of WT RT and EC₅₀ values for cytopathic protection of HIV-1-infected human T-cells in the 5-320 nM range.

Assembly of functionalized π-extended indolizine polycycles through dearomative [3+2] cycloaddition/oxidative decarbonylation

Reported herein is an unexpected construction of functionalized π-extended indolizine polycycles through a one-pot two-step cascade process comprising the base-promoted dearomative [3+2] cycloaddition of quinilinium salts and 3-alkenyl oxindoles, followed by a DDQ-mediated oxidative decarbonylation. Moreover, we could achieve the substrate-controlled diverse synthesis of structurally strained cyclopropane spirooxindole by using pyridinium salts as starting materials.

Highly stereoselective dearomative [3 + 2] cycloadditon of cyclic pyridinium ylides to access spiro-indolizidine scaffolds

A novel type of pyridinium salt bearing a EWG on the pyridine was developed as an efficient pyridinium ylide precursor in the [3 + 2] cycloaddition with nitroolefins to construct various spiro-indolizidine scaffolds via a dearomative pathway.

Palladium catalyzed synthesis of indolizines via the carbonylative coupling of bromopyridines, imines and alkynes

We report herein the development of a palladium-catalyzed, multicomponent synthesis of indolizines. The reaction proceeds via the carbonylative formation of a high energy, mesoionic pyridine-based 1,3-dipole, which can undergo spontaneous cycloaddition with alkynes. Overall, this provides a route to prepare indolizines in a modular fashion from combinations of commercially available or easily generated reagents: 2-bromopyridines, imines and alkynes.

A palladium catalyzed, multicomponent synthesis of indolizines is described via the carbon monoxide driven generation of reactive, pyridine-based 1,3-dipoles.

Structure activity relationship towards design of cryptosporidium specific thymidylate synthase inhibitors

Cryptosporidiosis is a human gastrointestinal disease caused by protozoans of the genus Cryptosporidium, which can be fatal in immunocompromised individuals. The essential enzyme, thymidylate synthase (TS), is responsible for de novo synthesis of deoxythymidine monophosphate. The TS active site is relatively conserved between Cryptosporidium and human enzymes. In previous work, we identified compound 1, (2-amino-4-oxo-4,7-dihydro-pyrrolo[2,3-d]pyrimidin-methyl-phenyl-l-glutamic acid), as a promising selective Cryptosporidium hominis TS (ChTS) inhibitor. In the present study, we explore the structure-activity relationship around 1 glutamate moiety by synthesizing and biochemically evaluating the inhibitory activity of analogues against ChTS and human TS (hTS). X-Ray crystal structures were obtained for compounds bound to both ChTS and hTS. We establish the importance of the 2-phenylacetic acid moiety methylene linker in optimally positioning compounds 23, 24, and 25 within the active site. Moreover, through the comparison of structural data for 5, 14, 15, and 23 bound in both ChTS and hTS identified that active site rigidity is a driving force in determining inhibitor selectivity.

Molecular and cellular studies evaluating a potent 2-cyanoindolizine catechol diether NNRTI targeting wildtype and Y181C mutant HIV-1 reverse transcriptase

The development of efficacious NNRTIs for HIV/AIDS therapy is commonly met with the emergence of drug resistant strains, including the Y181C variant. Using a computationally-guided approach, we synthesized the catechol diether series of NNRTIs, which display sub-nanomolar potency in cellular assays. Among the most potent were a series of 2-cyanoindolizine substituted catechol diethers, including Compound 1. We present here a thorough evaluation of this compound, including biochemical, cellular, and structural studies. The compound demonstrates low nanomolar potency against both WT and Y181C HIV-1 RT in in vitro and cellular assays. Our crystal structures of both the wildtype and mutant forms of RT in complex with Compound 1 allow the interrogation of this compound's features that allow it to maintain strong efficacy against the drug resistant mutant. Among these are compensatory shifts in the NNRTI binding pocket, persistence of multiple hydrogen bonds, and van der Waals contacts throughout the binding site. Further, the fluorine at the C6 position of the indolizine moiety makes multiple favorable interactions with both RT forms. The present study highlights the indolizine-substituted catechol diether class of NNRTIs as promising therapeutic candidates possessing optimal pharmacological properties and significant potency against multiple RT variants.

Current structure-based methods for designing non-nucleoside reverse transcriptase inhibitors

In 2019, structure-based methods continue to guide the design of novel antiretroviral therapies targeting HIV reverse transcriptase. This Review summarizes key findings from reverse transcriptase–non-nucleoside reverse transcriptase inhibitor analog crystal structure complexes reported from 2015 to 2019. Results from the literature and structure analysis have informed new ideas for structure-guided non-nucleoside reverse transcriptase inhibitor drug design.

Converting SMILES to Stacking Interaction Energies

Predicting the strength of stacking interactions involving heterocycles is vital for several fields, including structure-based drug design. While quantum chemical computations can provide accurate stacking interaction energies, these come at a steep computational cost. To address this challenge, we recently developed quantitative predictive models of stacking interactions between druglike heterocycles and the aromatic amino acids Phe, Tyr, and Trp (DOI: 10.1021/jacs.9b00936 ). These models depend on heterocycle descriptors derived from electrostatic potentials (ESPs) computed using density functional theory and provide accurate stacking interactions without the need for expensive computations on stacked dimers. Herein, we show that these ESP-based descriptors can be reliably evaluated directly from the atom connectivity of the heterocycle, providing a means of predicting both the descriptors and the potential for a given heterocycle to engage in stacking interactions without resorting to any quantum chemical computations. This enables the rapid conversion of simple molecular representations (e.g., SMILES) directly into accurate stacking interaction energies using a freely available online tool, thereby providing a way to rank the stacking abilities of large sets of heterocycles.

Structural and pharmacological evaluation of a novel non-nucleoside reverse transcriptase inhibitor as a promising long acting nanoformulation for treating HIV

Combination antiretroviral therapy (cART) has been proven effective in inhibiting human immunodeficiency virus type 1 (HIV-1) infection and has significantly improved the health outcomes in acquired immune deficiency syndrome (AIDS) patients. The therapeutic benefits of cART have been challenged because of the toxicity and emergence of drug-resistant HIV-1 strains along with lifelong patient compliance resulting in non-adherence. These issues also hinder the clinical benefits of non-nucleoside reverse transcriptase inhibitors (NNRTIs), which are one of the vital components of cART for the treatment of HIV-1 infection. In this study, using a computational and structural based drug design approach, we have discovered an effective HIV -1 NNRTI, compound I (Cmpd I) that is very potent in biochemical assays and which targets key residues in the allosteric binding pocket of wild-type (WT)-RT as revealed by structural studies. Furthermore, Cmpd I exhibited very potent antiviral activity in HIV-1 infected T cells, lacked cytotoxicity (therapeutic index >100,000), and no significant off-target effects were noted in pharmacological assays. To address the issue of non-adherence, we developed a long-acting nanoformulation of Cmpd I (Cmpd I-NP) using poly (lactide-coglycolide) (PLGA) particles. The pharmacokinetic studies of free and nanoformulated Cmpd I were carried out in BALB/c mice. Intraperitoneal administration of Cmpd I and Cmpd I-NP in BALB/c mice revealed prolonged serum residence time of 48 h and 30 days, respectively. The observed serum concentrations of Cmpd I in both cases were sufficient to provide >97% inhibition in HIV-1 infected T-cells. The significant antiviral activity along with favorable pharmacological and pharmacokinetic profile of Cmpd I, provide compelling and critical support for its further development as an anti-HIV therapeutic agent.

Computation of protein-ligand binding free energies using quantum mechanical bespoke force fields

A quantum mechanical bespoke molecular mechanics force field is derived for the L99A mutant of T4 lysozyme and used to compute absolute binding free energies of six benzene analogs to the protein. Promising agreement between theory and experiment highlights the potential for future use of system-specific force fields in computer-aided drug design.

Identification of Dihydrofuro[3,4- d]pyrimidine Derivatives as Novel HIV-1 Non-Nucleoside Reverse Transcriptase Inhibitors with Promising Antiviral Activities and Desirable Physicochemical Properties

To address drug resistance to HIV-1 non-nucleoside reverse transcriptase inhibitors (NNRTIs), a series of novel diarylpyrimidine (DAPY) derivatives targeting "tolerant region I" and "tolerant region II" of the NNRTIs binding pocket (NNIBP) were designed utilizing a structure-guided scaffold-hopping strategy. The dihydrofuro[3,4- d]pyrimidine derivatives 13c2 and 13c4 proved to be exceptionally potent against a wide range of HIV-1 strains carrying single NNRTI-resistant mutations (EC₅₀ = 0.9-8.4 nM), which were remarkably superior to that of etravirine (ETV). Meanwhile, both compounds exhibited comparable activities with ETV toward the virus with double mutations F227L+V106A and K103N+Y181C. Furthermore, the most active compound 13c2 showed favorable pharmacokinetic properties with an oral bioavailability of 30.96% and a half-life of 11.1 h, which suggested that 13c2 is worth further investigation as a novel NNRTI to circumvent drug resistance.

A computationally designed binding mode flip leads to a novel class of potent tri-vector cyclophilin inhibitors

Cyclophilins (Cyps) are a major family of drug targets that are challenging to prosecute with small molecules because the shallow nature and high degree of conservation of the active site across human isoforms offers limited opportunities for potent and selective inhibition. Herein a computational approach based on molecular dynamics simulations and free energy calculations was combined with biophysical assays and X-ray crystallography to explore a flip in the binding mode of a reported urea-based Cyp inhibitor. This approach enabled access to a distal pocket that is poorly conserved among key Cyp isoforms, and led to the discovery of a new family of sub-micromolar cell-active inhibitors that offer unprecedented opportunities for the development of next-generation drug therapies based on Cyp inhibition. The computational approach is applicable to a broad range of organic functional groups and could prove widely enabling in molecular design.

Targeting the hydrophobic channel of NNIBP: discovery of novel 1,2,3-triazole-derived diarylpyrimidines as novel HIV-1 NNRTIs with high potency against wild-type and K103N mutant virus

Novel 1,2,3-triazole-derived diarylpyrimidines were discovered as potent HIV-1 NNRTIs.

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.

Equilibrium Model of Drug-Modulated GagPol-Embedded HIV-1 Reverse Transcriptase Dimerization to Enhance Premature Protease Activation

Lack of effective strategies for killing cells latently infected with HIV-1 limits the eradication of AIDS. Unfortunately, current antiretroviral inhibitors are designed to target virus production but not latent infection. Interestingly, some non-nucleoside reverse transcriptase inhibitors (NNRTIs) have shown off-design effects, specifically, premature activation of HIV-1 protease (PR) within virus-infected cells that induces apoptosis. Here, we analyze an equilibrium model of HIV-1 reverse transcriptase (RT) binding to NNRTIs to understand the optimal binding characteristics that enhance RT dimerization within embedded GagPol dimers. This would allow NNRTIs to act as PR autoactivation enhancers (PAEs). We compute that ∼700-fold enhancement is theoretically possible by PAEs. Both a strong drug-dimer binding affinity (K_D12 < 100 nM) and relatively weaker drug-monomer affinity (K_D2/K_D12 > 10) are required for significant enhancement (∼50-fold or more) relative to the drug-free dimer concentration within a drug concentration limit of 10 μM. Our approach rationalizes the observed effects of efavirenz on premature activation of PR and may be useful to guide the design of suitable drug candidates and their optimal dosage regimens for this therapy class.

Predicting Binding Free Energies in a Large Combinatorial Chemical Space Using Multisite λ Dynamics

In this study, we demonstrate the extensive scalability of the biasing potential replica exchange multisite λ dynamics (BP-REX MSλD) free energy method by calculating binding affinities for 512 inhibitors to HIV Reverse Transcriptase (HIV-RT). This is the largest exploration of chemical space using free energy methods known to date, requires only a few simulations, and identifies 55 new inhibitor designs against HIV-RT predicted to be at least as potent as a tight binding reference compound (i.e., as potent as 56 nM). We highlight that BP-REX MSλD requires an order of magnitude less computational resources than conventional free energy methods while maintaining a similar level of precision, overcomes the inherent poor scalability of conventional free energy methods, and enables the exploration of combinatorially large chemical spaces in the context of in silico drug discovery.

Discovery of novel diarylpyrimidines as potent HIV-1 NNRTIs by investigating the chemical space of a less explored “hydrophobic channel”

We described the identification of novel HIV-1 NNRTIs via exploration of the chemical space of a seldom explored “hydrophobic channel”.

From in silico hit to long-acting late-stage preclinical candidate to combat HIV-1 infection

The HIV-1 pandemic affecting over 37 million people worldwide continues, with nearly one-half of the infected population on highly active antiretroviral therapy (HAART). Major therapeutic challenges remain because of the emergence of drug-resistant HIV-1 strains, limitations because of safety and toxicity with current HIV-1 drugs, and patient compliance for lifelong, daily treatment regimens. Nonnucleoside reverse transcriptase inhibitors (NNRTIs) that target the viral polymerase have been a key component of the current HIV-1 combination drug regimens; however, these issues hamper them. Thus, the development of novel more effective NNRTIs as anti-HIV-1 agents with fewer long-term liabilities, efficacy on new drug-resistant HIV-1 strains, and less frequent dosing is crucial. Using a computational and structure-based design strategy to guide lead optimization, a 5 µM virtual screening hit was transformed to a series of very potent nanomolar to picomolar catechol diethers. One representative, compound I, was shown to have nanomolar activity in HIV-1-infected T cells, potency on clinically relevant HIV-1 drug-resistant strains, lack of cytotoxicity and off-target effects, and excellent in vivo pharmacokinetic behavior. In this report, we show the feasibility of compound I as a late-stage preclinical candidate by establishing synergistic antiviral activity with existing HIV-1 drugs and clinical candidates and efficacy in HIV-1-infected humanized [human peripheral blood lymphocyte (Hu-PBL)] mice by completely suppressing viral loads and preventing human CD4⁺ T-cell loss. Moreover, a long-acting nanoformulation of compound I [compound I nanoparticle (compound I-NP)] in poly(lactide-coglycolide) (PLGA) was developed that shows sustained maintenance of plasma drug concentrations and drug efficacy for almost 3 weeks after a single dose.

Relative Binding Free Energy Calculations in Drug Discovery: Recent Advances and Practical Considerations

Accurate in silico prediction of protein-ligand binding affinities has been a primary objective of structure-based drug design for decades due to the putative value it would bring to the drug discovery process. However, computational methods have historically failed to deliver value in real-world drug discovery applications due to a variety of scientific, technical, and practical challenges. Recently, a family of approaches commonly referred to as relative binding free energy (RBFE) calculations, which rely on physics-based molecular simulations and statistical mechanics, have shown promise in reliably generating accurate predictions in the context of drug discovery projects. This advance arises from accumulating developments in the underlying scientific methods (decades of research on force fields and sampling algorithms) coupled with vast increases in computational resources (graphics processing units and cloud infrastructures). Mounting evidence from retrospective validation studies, blind challenge predictions, and prospective applications suggests that RBFE simulations can now predict the affinity differences for congeneric ligands with sufficient accuracy and throughput to deliver considerable value in hit-to-lead and lead optimization efforts. Here, we present an overview of current RBFE implementations, highlighting recent advances and remaining challenges, along with examples that emphasize practical considerations for obtaining reliable RBFE results. We focus specifically on relative binding free energies because the calculations are less computationally intensive than absolute binding free energy (ABFE) calculations and map directly onto the hit-to-lead and lead optimization processes, where the prediction of relative binding energies between a reference molecule and new ideas (virtual molecules) can be used to prioritize molecules for synthesis. We describe the critical aspects of running RBFE calculations, from both theoretical and applied perspectives, using a combination of retrospective literature examples and prospective studies from drug discovery projects. This work is intended to provide a contemporary overview of the scientific, technical, and practical issues associated with running relative binding free energy simulations, with a focus on real-world drug discovery applications. We offer guidelines for improving the accuracy of RBFE simulations, especially for challenging cases, and emphasize unresolved issues that could be improved by further research in the field.