Long Dissociation of Bictegravir from HIV-1 Integrase-DNA Complexes.

The HIV integrase (IN) strand transfer inhibitor (INSTI) bictegravir (BIC) has a long dissociation half-life (t1/2) from wild-type IN-DNA complexes: BIC 163 hr > dolutegravir (DTG) 96 hr > raltegravir (RAL) 10 hr > elvitegravir (EVG) 3.3 hr. In cells, BIC had more durable antiviral activity against wild-type HIV after drug washout than RAL or EVG. BIC also had a longer t1/2 and maintained longer antiviral activity after drug washout than DTG with the clinically relevant resistance IN mutant G140S+Q148H. Structural analyses indicate that BIC makes more contacts with the IN-DNA complex than DTG mainly via its bicyclic ring system which may contribute to more prolonged residence time and resilience against many resistance mutations.

Evaluation of Free Energy Calculations for the Prioritization of Macrocycle Synthesis.

A tremendous research and development effort was exerted toward combating chronic hepatitis C, ultimately leading to curative oral treatments, all of which are targeting viral proteins. Despite the advantage of numerous targets allowing for broad hepatitis C virus (HCV) genotype coverage, the only host target inhibitors that advanced into clinical development were Cyclosporin A based cyclophilin inhibitors. While cyclosporin-based molecules typically require a fermentation process, Gilead successfully pursued a fully synthetic, oral program based on Sanglifehrin A. The drug discovery process, though greatly helped by facile crystallography, was still hampered by the limitations in the accuracy of predictive computational methods for prioritizing compound ideas. Recent advances in accuracy and speed of free energy perturbation (FEP) methods, however, are attractive for prioritizing and derisking synthetically challenging molecules and potentially could have had a significant impact on the speed of the development of this program. Here in our simulated prospective study, the binding free energies of 26 macrocyclic cyclophilin inhibitors were blindly predicted using FEP+ to test this hypothesis. The predictions had a low mean unsigned error (MUE) (1.1 kcal/mol) and accurately reproduced many design decisions from the program, suggesting that FEP+ has the potential to drive synthetic chemistry efforts by more accurately ranking compounds with nonintuitive structure-activity relationships (SARs).

Targeting prostate cancer with compounds possessing dual activity as androgen receptor antagonists and HDAC6 inhibitors.

While enzalutamide and abiraterone are approved for treatment of metastatic castration-resistant prostate cancer (mCRPC), approximately 20-40% of patients have no response to these agents. It has been stipulated that the lack of response and the development of secondary resistance to these drugs may be due to the presence of AR splice variants. HDAC6 has a role in regulating the androgen receptor (AR) by modulating heat shock protein 90 (Hsp90) acetylation, which controls the nuclear localization and activation of the AR in androgen-dependent and independent scenarios. With dual-acting AR-HDAC6 inhibitors it should be possible to target patients who don't respond to enzalutamide. Herein, we describe the design, synthesis and biological evaluation of dual-acting compounds which target AR and are also specific towards HDAC6. Our efforts led to compound 10 which was found to have potent dual activity (HDAC6 IC50=0.0356µM and AR binding IC50=<0.03µM). Compound 10 was further evaluated for antagonist and other cell-based activities, in vitro stability and pharmacokinetics.

Computer-aided drug design of falcipain inhibitors: virtual screening, structure-activity relationships, hydration site thermodynamics, and reactivity analysis.

Falcipains (FPs) are hemoglobinases of Plasmodium falciparum that are validated targets for the development of antimalarial chemotherapy. A combined ligand- and structure-based virtual screening of commercial databases was performed to identify structural analogs of virtual screening hits previously discovered in our laboratory. A total of 28 low micromolar inhibitors of FP-2 and FP-3 were identified and the structure-activity relationship (SAR) in each series was elaborated. The SAR of the compounds was unusually steep in some cases and could not be explained by a traditional analysis of the ligand-protein interactions (van der Waals, electrostatics, and hydrogen bonds). To gain further insights, a statistical thermodynamic analysis of explicit solvent in the ligand binding domains of FP-2 and FP-3 was carried out to understand the roles played by water molecules in binding of these inhibitors. Indeed, the energetics associated with the displacement of water molecules upon ligand binding explained some of the complex trends in the SAR. Furthermore, low potency of a subset of FP-2 inhibitors that could not be understood by the water energetics was explained in the context of poor chemical reactivity of the reactive centers of these compounds. The present study highlights the importance of considering energetic contributors to binding beyond traditional ligand-protein interactions.

Improving the Prediction of Absolute Solvation Free Energies Using the Next Generation OPLS Force Field.

Explicit solvent molecular dynamics free energy perturbation simulations were performed to predict absolute solvation free energies of 239 diverse small molecules. We use OPLS2.0, the next generation OPLS force field, and compare the results with popular small molecule force fields-OPLS_2005, GAFF, and CHARMm-MSI. OPLS2.0 produces the best correlation with experimental data (R(2) = 0.95, slope = 0.96) and the lowest average unsigned errors (0.7 kcal/mol). Important classes of compounds that performed suboptimally with OPLS_2005 show significant improvements.

Prediction of Absolute Solvation Free Energies using Molecular Dynamics Free Energy Perturbation and the OPLS Force Field.

The accurate prediction of protein-ligand binding free energies is a primary objective in computer-aided drug design. The solvation free energy of a small molecule provides a surrogate to the desolvation of the ligand in the thermodynamic process of protein-ligand binding. Here, we use explicit solvent molecular dynamics free energy perturbation to predict the absolute solvation free energies of a set of 239 small molecules, spanning diverse chemical functional groups commonly found in drugs and drug-like molecules. We also compare the performance of absolute solvation free energies obtained using the OPLS_2005 force field with two other commonly used small molecule force fields-general AMBER force field (GAFF) with AM1-BCC charges and CHARMm-MSI with CHelpG charges. Using the OPLS_2005 force field, we obtain high correlation with experimental solvation free energies (R(2) = 0.94) and low average unsigned errors for a majority of the functional groups compared to AM1-BCC/GAFF or CHelpG/CHARMm-MSI. However, OPLS_2005 has errors of over 1.3 kcal/mol for certain classes of polar compounds. We show that predictions on these compound classes can be improved by using a semiempirical charge assignment method with an implicit bond charge correction.

Computations of Absolute Solvation Free Energies of Small Molecules Using Explicit and Implicit Solvent Model.

Accurate determination of absolute solvation free energy plays a critical role in numerous areas of biomolecular modeling and drug discovery. A quantitative representation of ligand and receptor desolvation, in particular, is an essential component of current docking and scoring methods. Furthermore, the partitioning of a drug between aqueous and nonpolar solvents is one of the important factors considered in pharmacokinetics. In this study, the absolute hydration free energy for a set of 239 neutral ligands spanning diverse chemical functional groups commonly found in drugs and drug-like candidates is calculated using the molecular dynamics free energy perturbation method (FEP/MD) with explicit water molecules, and compared to experimental data as well as its counterparts obtained using implicit solvent models. The hydration free energies are calculated from explicit solvent simulations using a staged FEP procedure permitting a separation of the total free energy into polar and nonpolar contributions. The nonpolar component is further decomposed into attractive (dispersive) and repulsive (cavity) components using the Weeks-Chandler-Anderson (WCA) separation scheme. To increase the computational efficiency, all of the FEP/MD simulations are generated using a mixed explicit/implicit solvent scheme with a relatively small number of explicit TIP3P water molecules, in which the influence of the remaining bulk is incorporated via the spherical solvent boundary potential (SSBP). The performances of two fixed-charge force fields designed for small organic molecules, the General Amber force field (GAFF), and the all-atom CHARMm-MSI, are compared. Because of the crucial role of electrostatics in solvation free energy, the results from various commonly used charge generation models based on the semiempirical (AM1-BCC) and QM calculations [charge fitting using ChelpG and RESP] are compared. In addition, the solvation free energies of the test set are also calculated using Poisson-Boltzmann (PB) and Generalized Born model of solvation (GB), which are two widely used continuum electrostatic implicit solvent models. The protocol for running the absolute solvation free energy calculations used throughout is automated as much as possible, with minimum user intervention, so that it can be used in large-scale analysis and force field optimization.

Rescoring docking hit lists for model cavity sites: predictions and experimental testing.

Molecular docking computationally screens thousands to millions of organic molecules against protein structures, looking for those with complementary fits. Many approximations are made, often resulting in low "hit rates." A strategy to overcome these approximations is to rescore top-ranked docked molecules using a better but slower method. One such is afforded by molecular mechanics-generalized Born surface area (MM-GBSA) techniques. These more physically realistic methods have improved models for solvation and electrostatic interactions and conformational change compared to most docking programs. To investigate MM-GBSA rescoring, we re-ranked docking hit lists in three small buried sites: a hydrophobic cavity that binds apolar ligands, a slightly polar cavity that binds aryl and hydrogen-bonding ligands, and an anionic cavity that binds cationic ligands. These sites are simple; consequently, incorrect predictions can be attributed to particular errors in the method, and many likely ligands may actually be tested. In retrospective calculations, MM-GBSA techniques with binding-site minimization better distinguished the known ligands for each cavity from the known decoys compared to the docking calculation alone. This encouraged us to test rescoring prospectively on molecules that ranked poorly by docking but that ranked well when rescored by MM-GBSA. A total of 33 molecules highly ranked by MM-GBSA for the three cavities were tested experimentally. Of these, 23 were observed to bind--these are docking false negatives rescued by rescoring. The 10 remaining molecules are true negatives by docking and false positives by MM-GBSA. X-ray crystal structures were determined for 21 of these 23 molecules. In many cases, the geometry prediction by MM-GBSA improved the initial docking pose and more closely resembled the crystallographic result; yet in several cases, the rescored geometry failed to capture large conformational changes in the protein. Intriguingly, rescoring not only rescued docking false positives, but also introduced several new false positives into the top-ranking molecules. We consider the origins of the successes and failures in MM-GBSA rescoring in these model cavity sites and the prospects for rescoring in biologically relevant targets.

Computational study of IAG-nucleoside hydrolase: determination of the preferred ground state conformation and the role of active site residues.

The mechanism of action of inosine-adenosine-guanosine nucleoside hydrolase (IAG-NH) has been investigated by long-term molecular dynamics (MD) simulation in TIP3P water using stochastic boundary conditions. Special attention has been given to the role of leaving group pocket residues and conformation of the bound substrate at the active site of IAG-NH. We also describe the positioning of the residues of an important flexible loop at the active site, which was previously unobservable by X-ray crystallography due to high B-factors. Five MD simulations have been performed with the Enzyme x Substrate complexes: Enzyme x anti-Adenosine with Asp40-COOH [E(40H) x Ade(a)], Enzyme x anti-Adenosine with Asp40-COO- [E(40) x Ade(a)], Enzyme x syn-Adenosine with Asp40-COOH [E(40H) x Ade(s)], Enzyme x syn-Adenosine with Asp40-COO- [E(40) x Ade(s)], and Enzyme x anti-Inosine with Asp40-COO- [E(40) x Ino(a)]. Overall, the structures generated from the MD simulation of E(40H) x Ade(s) preserve the catalytically important hydrogen bonds as well as electrostatic and hydrophobic interactions to provide a plausible catalytic structure. When deprotonated Asp40 (Asp4-COO-) is present, the active site is open to water solvent which interferes with the base stacking between Trp83 and nucleobase. A calculation using Poisson-Boltzmann equation module supports that Asp40 indeed has an elevated pK(app). Solvent accessible surface area (SASA) calculations on all the five MD structures shows that systems with protonated Asp40, namely, E(40H) x Ade(a) and E(40H) x Ade(s), have zero SASA. It is found that a water molecule is hydrogen-bonded to the N7 of the nucleobase and is probably the essential general acid to protonate the departing nucleobase anion. The N7-bonded water is in turn hydrogen-bonded to waters in a channel, held in place by the residues of the flexible loop, Tyr257, His247, and Cys245. Using normal-mode analysis with elastic network model, we find that the flexible loop explores a conformational space much larger than in the MD trajectory, leading to a "gating"-like motion with respect to the active site.

Molecular dynamics studies of ground state and intermediate of the hyperthermophilic indole-3-glycerol phosphate synthase.

Indole-3-glycerol phosphate synthase catalyzes the terminal ring closure step in tryptophan biosynthesis. In this paper, we compare the results from molecular dynamics (MD) simulations of enzyme-bound substrate at 298, 333, 363, and 385 K and the enzyme-bound intermediate at 385 K, solvated in TIP3P water box with a CHARMM force field. Results from MD simulations agree with experimental studies supporting the observation that Lys-110 is the general acid. Based on its location in the active site during the MD simulations, Glu-210 warrants classification as the general base instead of the previously proposed Glu-159. We find that the relative population of the reactive enzyme-substrate Michaelis conformers [near attack conformers (NACs)] with temperature correlates well (correlation coefficient of 0.96) with the relative activity of this thermophilic enzyme. At higher temperature, the enzyme-substrate electrostatic interaction favors the binding of the substrate in NAC conformation, whereas, at lower temperature, the substrate is distorted and bound in a nonreactive conformation. This change is reflected in the approximately 1,100-fold increase in population of NACs at 385 K relative to 298 K. The easily determined population of NACs at given temperature tells much about the thermophilic property of the enzyme. Thus, the hyperthermophilic enzyme has evolved to have optimum activity at high temperatures, and, with lowering of the temperature, the electrostatic interaction at the active site is enhanced and the structure is deformed. This model can be regarded as a general explanation for the activity of hyperthermophilic enzymes.

Computational study of the ground state of thermophilic indole glycerol phosphate synthase: structural alterations at the active site with temperature.

Hyperthermophlic indole-3-glycerol phosphate synthase (IGPS) catalyzes the terminal ring-closure step in tryptophan biosynthesis. In this paper, we compare the results from the molecular dynamics (MD) simulation of enzyme-bound substrate at 298 K (E.S298) and 385 K (E.S385) solvated in the TIP3P water box using the CHARMM force field to address the question of the structural change of the Enzyme. Substrate complex with temperature. The population of the reactive Enzyme. Substrate conformers (near attack conformers or NACs) increases by approximately 1100-fold in going from room temperature (E.S298) to high temperature (E.S385). This increased population of NAC conformers in the Michaelis complex correlates well with the increase in rate in going from 298 to 385 K. The positioning of the two active site residues Lys53 and Lys110 controls binding of the substrate in the favorable orientation for general acid-catalyzed intramolecular ring formation reaction. It can be concluded that the NAC formation allowing general acid catalysis has much to do with the temperature dependence of the free energy of reaction.

Probing the Ser-Ser-Lys catalytic triad mechanism of peptide amidase: computational studies of the ground state, transition state, and intermediate.

Peptide amidase (Pam), a hydrolytic enzyme that belongs to the amidase signature (AS) family, selectively catalyzes the hydrolysis of the C-terminal amide bond (CO-NH(2)) of peptides. The recent availability of the X-ray structures of Pam, fatty acid amide hydrolase, and malonamidase E2 has led to the proposal of a novel Ser-Ser-Lys catalytic triad mechanism for the amide hydrolysis by the AS enzymes. The molecular dynamics (MD) simulations using the CHARMM force field were performed to explore the catalytic mechanism of Pam. The 1.8 A X-ray crystal structure of Pam in complex with the amide analogue of chymostatin was chosen for the initial coordinates for the MD simulations. The five systems that were investigated are as follows: (i) enzyme.substrate with Lys123-NH(2), (ii) enzyme.substrate with Lys123-NH(3)(+), (iii) enzyme.substrate with Lys123-NH(3)(+) and Ser226-O(-), (iv) enzyme.transition state, and (v) enzyme.tetrahedral intermediate. Our data support the presence of the hydrogen bonding network among the catalytic triad residues, Ser226, Ser202, and Lys123, where Ser226 acts as the nucleophile and Ser202 bridges Ser226 and Lys123. The MD simulation supports the catalytic role of the crystallographic waters, Wat1 and Wat2. In all the systems that have been studied, the backbone amide nitrogens of Asp224 and Thr223 create an oxyanion hole by hydrogen bonding to the terminal amide oxygen of the substrate, and stabilize the oxyanion tetrahedral intermediate. The results from both our computational investigation and previously published experimental pH profile support two mechanisms. In a mechanism that is relevant at lower pH, the Lys123-NH(3)(+)-Ser202 dyad provides structural support to the catalytic residue Ser226, which in turn carries out a nucleophilic attack at the substrate amide carbonyl in concert with Wat1-mediated deprotonation and stabilization of the tetrahedral transition state by the oxyanion hole. In the mechanism operating at higher pH, the Lys123-NH(2)-Ser202 catalytic dyad acts as a general base to assist addition of Ser226 to the substrate amide carbonyl. The results from the MD simulation of the tetrahedral intermediate state show that both Ser202 and Lys123 are possible candidates for protonation of the leaving group, NH(2), to form the acyl-enzyme intermediate.