Multiparameter optimization in CNS drug discovery: design of pyrimido[4,5-d]azepines as potent 5-hydroxytryptamine 2C (5-HT2C) receptor agonists with exquisite functional selectivity over 5-HT2A and 5-HT2B receptors.

A series of 4-substituted pyrimido[4,5-d]azepines that are potent, selective 5-HT2C receptor partial agonists is described. A rational medicinal chemistry design strategy to deliver CNS penetration coupled with SAR-based optimization of selectivity and agonist potency provided compounds with the desired balance of preclinical properties. Lead compounds 17 (PF-4479745) and 18 (PF-4522654) displayed robust pharmacology in a preclinical canine model of stress urinary incontinence (SUI) and no measurable functional agonism at the key selectivity targets 5-HT2A and 5-HT2B in relevant tissue-based assay systems. Utilizing recent advances in the structural biology of GPCRs, homology modeling has been carried out to rationalize binding and agonist efficacy of these compounds.

Variational Implicit-Solvent Modeling of Host-Guest Binding: A Case Study on Cucurbit[7]uril|

The synthetic host cucurbit[7]uril (CB[7]) binds aromatic guests or metal complexes with ultrahigh affinity compared with that typically displayed in protein-ligand binding. Due to its small size, CB[7] serves as an ideal receptor-ligand system for developing computational methods for molecular recognition. Here, we apply the recently developed variational implicit-solvent model (VISM), numerically evaluated by the level-set method, to study hydration effects in the high-affinity binding of the B2 bicyclo[2.2.2]octane derivative to CB[7]. For the unbound host, we find that the host cavity favors the hydrated state over the dry state due to electrostatic effects. For the guest binding, we find reasonable agreement to experimental binding affinities. Dissection of the individual VISM free-energy contributions shows that the major driving forces are water-mediated hydrophobic interactions and the intrinsic (vacuum) host-guest van der Waals interactions. These findings are in line with recent experiments and molecular dynamics simulations with explicit solvent. It is expected that the level-set VISM, with further refinement on the electrostatic descriptions, can efficiently predict molecular binding and recognition in a wide range of future applications.

Inactivating mutation in histone deacetylase 3 stabilizes its active conformation.

Histone deacetylases (HDACs), together with histone acetyltransferases (HATs), regulate gene expression by modulating the acetylation level of chromatin. HDAC3 is implicated in many important cellular processes, particularly in cancer cell proliferation and metastasis, making inhibition of HDAC3 a promising epigenetic treatment for certain cancers. HDAC3 is activated upon complex formation with both inositol tetraphosphate (IP4) and the deacetylase-activating domain (DAD) of multi-protein nuclear receptor corepressor complexes. In previous studies, we have shown that binding of DAD and IP4 to HDAC3 significantly restricts its conformational space towards its stable ternary complex conformation, and suggest this to be the active conformation. Here, we report a single mutation of HDAC3 that is capable of mimicking the stabilizing effects of DAD and IP4, without the presence of either. This mutation, however, results in a total loss of deacetylase activity, prompting a closer evaluation of our understanding of the activation of HDAC3.

Substrate-dependent dynamics of UDP-galactopyranose mutase: Implications for drug design.

Trypanosoma cruzi is the causative agent of Chagas disease, a neglected tropical disease that represents one of the major health challenges of the Latin American countries. Successful efforts were made during the last few decades to control the transmission of this disease, but there is still no treatment for the 10 million adults in the chronic phase of the disease. In T. cruzi, as well as in other pathogens, the flavoenzyme UDP-galactopyranose mutase (UGM) catalyzes the conversion of UDP-galactopyranose to UDP-galactofuranose, a precursor of the cell surface ß-galactofuranose that is involved in the virulence of the pathogen. The fact that UGM is not present in humans makes inhibition of this enzyme a good approach in the design of new Chagas therapeutics. By performing a series of computer simulations of T. cruzi UGM in the presence or absence of an active site ligand, we address the molecular details of the mechanism that controls the uptake and retention of the substrate. The simulations suggest a modular mechanism in which each moiety of the substrate controls the flexibility of a different protein loop. Furthermore, the calculations indicate that interactions with the substrate diphosphate moiety are especially important for stabilizing the closed active site. This hypothesis is supported with kinetics measurements of site-directed mutants of T. cruzi UGM. Our results extend our knowledge of UGM dynamics and offer new alternatives for the prospective design of drugs.

Population based reweighting of scaled molecular dynamics.

Molecular dynamics simulation using enhanced sampling methods is one of the powerful computational tools used to explore protein conformations and free energy landscapes. Enhanced sampling methods often employ either an increase in temperature or a flattening of the potential energy surface to rapidly sample phase space, and a corresponding reweighting algorithm is used to recover the Boltzmann statistics. However, potential energies of complex biomolecules usually involve large fluctuations on a magnitude of hundreds of kcal/mol despite minimal structural changes during simulation. This leads to noisy reweighting statistics and complicates the obtainment of accurate final results. To overcome this common issue in enhanced conformational sampling, we propose a scaled molecular dynamics method, which modifies the biomolecular potential energy surface and employs a reweighting scheme based on configurational populations. Statistical mechanical theory is applied to derive the reweighting formula, and the canonical ensemble of simulated structures is recovered accordingly. Test simulations on alanine dipeptide and the fast folding polypeptide Chignolin exhibit sufficiently enhanced conformational sampling and accurate recovery of free energy surfaces and thermodynamic properties. The results are comparable to long conventional molecular dynamics simulations and exhibit better recovery of canonical statistics over methods which employ a potential energy term in reweighting.

w-REXAMD: A Hamiltonian Replica Exchange Approach to Improve Free Energy Calculations for Systems with Kinetically Trapped Conformations.

Free energy governs the equilibrium extent of many biological processes. High barriers separating free energy minima often limit the sampling in molecular dynamics (MD) simulations, leading to inaccurate free energies. Here, we demonstrate enhanced sampling and improved free energy calculations, relative to conventional MD, using windowed accelerated MD within a Hamiltonian replica exchange framework (w-REXAMD). We show that for a case in which multiple conformations are separated by large free energy barriers, w-REXAMD is a useful enhanced sampling technique, without any necessary reweighting.

On the Role of Dewetting Transitions in Host-Guest Binding Free Energy Calculations.

We use thermodynamic integration (TI) and explicit solvent molecular dynamics (MD) simulation to estimate the absolute free energy of host-guest binding. In the unbound state, water molecules visit all of the internally accessible volume of the host, which is fully hydrated on all sides. Upon binding of an apolar guest, the toroidal host cavity is fully dehydrated; thus, during the intermediate λ stages along the integration, the hydration of the host fluctuates between hydrated and dehydrated states. Estimating free energies by TI can be especially challenging when there is a considerable difference in hydration between the two states of interest. We investigate these aspects using the popular TIP3P and TIP4P water models. TI free energy estimates through MD largely depend on water-related interactions, and water dynamics significantly affect the convergence of binding free energy calculations. Our results indicate that wetting/dewetting transitions play a major role in slowing the convergence of free energy estimation. We employ two alternative approaches-one analytical and the other empirically based on actual MD sampling-to correct for the standard state free energy. This correction is sizable (up to 4 kcal/mol), and the two approaches provide corrections that differ by about 1 kcal/mol. For the system considered here, the TIP4P water model combined with an analytical correction for the standard state free energy provides higher overall accuracy. This observation might be transferable to other systems in which water-related contributions dominate the binding process.

Structural insight into the separate roles of inositol tetraphosphate and deacetylase-activating domain in activation of histone deacetylase 3.

Histone deacetylases (HDACs) repress transcription by deacetylating acetyllysines on specific histone tails. HDAC3 is implicated in neurodegenerative diseases, certain leukemias, and even in disrupting HIV-1 latency. A recent crystal structure of HDAC3 in complex with the deacetylase-activating domain (DAD) of its corepressor complex revealed an inositol tetraphosphate (IP4) molecule at the protein-protein interface. IP4 was shown to play an important, yet mechanistically ambiguous, role in the activity of HDAC3. The goal of this article is to explore the conformational ensemble of HDAC3 in its inactive apo state and in the presence of each or both of DAD and IP4. Using triplicate, 100 ns molecular dynamic simulations, we study the apo, ternary, and intermediate DAD-bound or IP4-bound HDAC3 states. We find that a population-shift effect is induced by the presence of each corepressor, and is most notable in the presence of both. Our results offer new insights into the change in dynamics necessary for the activation of HDAC3 and highlight the roles of IP4 and DAD in this process.

Antibacterial drug leads targeting isoprenoid biosynthesis.

With the rise in resistance to antibiotics such as methicillin, there is a need for new drugs. We report here the discovery and X-ray crystallographic structures of 10 chemically diverse compounds (benzoic, diketo, and phosphonic acids, as well as a bisamidine and a bisamine) that inhibit bacterial undecaprenyl diphosphate synthase, an essential enzyme involved in cell wall biosynthesis. The inhibitors bind to one or more of the four undecaprenyl diphosphate synthase inhibitor binding sites identified previously, with the most active leads binding to site 4, outside the catalytic center. The most potent leads are active against Staphylococcus aureus [minimal inhibitory concentration (MIC)(90) ∼0.25 µg/mL], and one potently synergizes with methicillin (fractional inhibitory concentration index = 0.25) and is protective in a mouse infection model. These results provide numerous leads for antibacterial development and open up the possibility of restoring sensitivity to drugs such as methicillin, using combination therapies.

Drug hypersensitivity caused by alteration of the MHC-presented self-peptide repertoire.

Idiosyncratic adverse drug reactions are unpredictable, dose-independent and potentially life threatening; this makes them a major factor contributing to the cost and uncertainty of drug development. Clinical data suggest that many such reactions involve immune mechanisms, and genetic association studies have identified strong linkages between drug hypersensitivity reactions to several drugs and specific HLA alleles. One of the strongest such genetic associations found has been for the antiviral drug abacavir, which causes severe adverse reactions exclusively in patients expressing the HLA molecular variant B*57:01. Abacavir adverse reactions were recently shown to be driven by drug-specific activation of cytokine-producing, cytotoxic CD8(+) T cells that required HLA-B*57:01 molecules for their function; however, the mechanism by which abacavir induces this pathologic T-cell response remains unclear. Here we show that abacavir can bind within the F pocket of the peptide-binding groove of HLA-B*57:01, thereby altering its specificity. This provides an explanation for HLA-linked idiosyncratic adverse drug reactions, namely that drugs can alter the repertoire of self-peptides presented to T cells, thus causing the equivalent of an alloreactive T-cell response. Indeed, we identified specific self-peptides that are presented only in the presence of abacavir and that were recognized by T cells of hypersensitive patients. The assays that we have established can be applied to test additional compounds with suspected HLA-linked hypersensitivities in vitro. Where successful, these assays could speed up the discovery and mechanistic understanding of HLA-linked hypersensitivities, and guide the development of safer drugs.

Protecting High Energy Barriers: A New Equation to Regulate Boost Energy in Accelerated Molecular Dynamics Simulations.

Molecular dynamics (MD) is one of the most common tools in computational chemistry. Recently, our group has employed accelerated molecular dynamics (aMD) to improve the conformational sampling over conventional molecular dynamics techniques. In the original aMD implementation, sampling is greatly improved by raising energy wells below a predefined energy level. Recently, our group presented an alternative aMD implementation where simulations are accelerated by lowering energy barriers of the potential energy surface. When coupled with thermodynamic integration simulations, this implementation showed very promising results. However, when applied to large systems, such as proteins, the simulation tends to be biased to high energy regions of the potential landscape. The reason for this behavior lies in the boost equation used since the highest energy barriers are dramatically more affected than the lower ones. To address this issue, in this work, we present a new boost equation that prevents oversampling of unfavorable high energy conformational states. The new boost potential provides not only better recovery of statistics throughout the simulation but also enhanced sampling of statistically relevant regions in explicit solvent MD simulations.

Large-scale conformational changes of Trypanosoma cruzi proline racemase predicted by accelerated molecular dynamics simulation.

Chagas' disease, caused by the protozoan parasite Trypanosoma cruzi (T. cruzi), is a life-threatening illness affecting 11-18 million people. Currently available treatments are limited, with unacceptable efficacy and safety profiles. Recent studies have revealed an essential T. cruzi proline racemase enzyme (TcPR) as an attractive candidate for improved chemotherapeutic intervention. Conformational changes associated with substrate binding to TcPR are believed to expose critical residues that elicit a host mitogenic B-cell response, a process contributing to parasite persistence and immune system evasion. Characterization of the conformational states of TcPR requires access to long-time-scale motions that are currently inaccessible by standard molecular dynamics simulations. Here we describe advanced accelerated molecular dynamics that extend the effective simulation time and capture large-scale motions of functional relevance. Conservation and fragment mapping analyses identified potential conformational epitopes located in the vicinity of newly identified transient binding pockets. The newly identified open TcPR conformations revealed by this study along with knowledge of the closed to open interconversion mechanism advances our understanding of TcPR function. The results and the strategy adopted in this work constitute an important step toward the rationalization of the molecular basis behind the mitogenic B-cell response of TcPR and provide new insights for future structure-based drug discovery.

Enhanced Lipid Diffusion and Mixing in Accelerated Molecular Dynamics.

Accelerated molecular dynamics (aMD) is an enhanced sampling technique that expedites conformational space sampling by reducing the barriers separating various low-energy states of a system. Here, we present the first application of the aMD method on lipid membranes. Altogether, ∼1.5 µs simulations were performed on three systems: a pure POPC bilayer, a pure DMPC bilayer, and a mixed POPC:DMPC bilayer. Overall, the aMD simulations are found to produce significant speedup in trans-gauche isomerization and lipid lateral diffusion versus those in conventional MD (cMD) simulations. Further comparison of a 70-ns aMD run and a 300-ns cMD run of the mixed POPC:DMPC bilayer shows that the two simulations yield similar lipid mixing behaviors, with aMD generating a 2-3-fold speedup compared to cMD. Our results demonstrate that the aMD method is an efficient approach for the study of bilayer structural and dynamic properties. On the basis of simulations of the three bilayer systems, we also discuss the impact of aMD parameters on various lipid properties, which can be used as a guideline for future aMD simulations of membrane systems.

From Zn to Mn: the study of novel manganese-binding groups in the search for new drugs against tuberculosis.

In most eubacteria, apicomplexans, and most plants, including the causal agents for diseases such as malaria, leprosy, and tuberculosis, the methylerythritol phosphate pathway is the route for the biosynthesis of the C(5) precursors to the essential isoprenoid class of compounds. Owing to their absence in humans, the enzymes of the methylerythritol phosphate pathway have become attractive targets for drug discovery. This work investigates a new class of inhibitors against the second enzyme of the pathway, 1-deoxy-D-xylulose 5-phosphate reductoisomerase. Inhibition of this enzyme may involve the chelation of a crucial active site Mn ion, and the metal-chelating moieties studied here have previously been shown to be successful in application to the zinc-dependent metalloproteinases. Quantum mechanics and docking calculations presented in this work suggest the transferability of these metal-chelating compounds to Mn-containing 1-deoxy-D-xylulose 5-phosphate reductoisomerase enzyme, as a promising starting point to the development of potent inhibitors.

POVME: an algorithm for measuring binding-pocket volumes.

Researchers engaged in computer-aided drug design often wish to measure the volume of a ligand-binding pocket in order to predict pharmacology. We have recently developed a simple algorithm, called POVME (POcket Volume MEasurer), for this purpose. POVME is Python implemented, fast, and freely available. To demonstrate its utility, we use the new algorithm to study three members of the matrix-metalloproteinase family of proteins. Despite the structural similarity of these proteins, differences in binding-pocket dynamics are easily identified.

From sensors to silencers: quinoline- and benzimidazole-sulfonamides as inhibitors for zinc proteases.

Derived from the extensive work in the area of small molecule zinc(II) ion sensors, chelating fragment libraries of quinoline- and benzimidazole-sulfonamides have been prepared and screened against several different zinc(II)-dependent matrix metalloproteinases (MMPs). The fragments show impressive inhibition of these metalloenzymes and preferences for different MMPs based on the nature of the chelating group. The findings show that focused chelator libraries are a powerful strategy for the discovery of lead fragments for metalloprotein inhibition.

Coupling Constant pH Molecular Dynamics with Accelerated Molecular Dynamics.

An extension of the constant pH method originally implemented by Mongan et al. (J. Comput. Chem.2004, 25, 2038-2048) is proposed in this study. This adapted version of the method couples the constant pH methodology with the enhanced sampling technique of accelerated molecular dynamics, in an attempt to overcome the sampling issues encountered with current standard constant pH molecular dynamics methods. Although good results were reported by Mongan et al. on application of the standard method to the hen egg-white lysozyme (HEWL) system, residues which possess strong interactions with neighboring groups tend to converge slowly, resulting in the reported inconsistencies for predicted pK(a) values, as highlighted by the authors. The application of the coupled method described in this study to the HEWL system displays improvements over the standard version of the method, with the improved sampling leading to faster convergence and producing pK(a) values in closer agreement to those obtained experimentally for the more slowly converging residues.

Including receptor flexibility and induced fit effects into the design of MMP-2 inhibitors.

Matrix metalloproteinases (MMPs) comprise a class of flexible proteins required for normal tissue remodeling. Overexpression of MMPs is associated with a wide range of pathophysiological processes, including vascular disease, multiple sclerosis, Alzheimer's disease, and cancer. Nearly all MMP inhibitors have failed in clinical trials, in part due to lack of specificity. Due to the highly dynamic molecular motions of the MMP-2 binding pockets, the rational drug design of MMP inhibitors has been very challenging. To address these challenges, in the current study we combine computer docking with molecular dynamics (MD) simulations in order to incorporate receptor-flexibility and induced-fit effects into the drug-design process. Our strategy identifies molecular fragments predicted to target multiple MMP-2 binding pockets.

Thioamide hydroxypyrothiones supersede amide hydroxypyrothiones in potency against anthrax lethal factor.

Anthrax lethal factor (LF) is a critical virulence factor in the pathogenesis of anthrax. A structure-activity relationship (SAR) of potential lethal factor inhibitors (LFi) is presented in which the zinc-binding group (ZBG), linker, and backbone moieties for a series of hydroxypyrone-based compounds were systematically varied. It was found that hydroxypyrothione ZBGs generate more potent inhibitors than hydroxypyrone ZBGs. Furthermore, coupling the hydroxypyrothione to a backbone group via a thioamide bond improves potency when compared to an amide linker. QM/MM studies show that the thioamide bond in these inhibitors allows for the formation of two additional hydrogen bonds with the protein active site. In both types of hydroxypyrothione compounds, ligand efficiencies of 0.29-0.54 kcal mol(-1) per heavy atom were achieved. The results highlight the need for a better understanding to optimize the interplay between the ZBG, linker, and backbone to get improved LFi.

Inhibition of cathepsin B by Au(I) complexes: a kinetic and computational study.

Gold(I) compounds have been used in the treatment of rheumatoid arthritis for over 80 years, but the biological targets and the structure-activity relationships of these drugs are not well understood. Of particular interest is the molecular mechanism behind the antiarthritic activity of the orally available drug triethylphosphine(2,3,4,6-tetra-O-acetyl-beta-l-D-thiopyranosato-S) gold(I) (auranofin, Ridaura). The cathepsin family of lysosomal, cysteine-dependent enzymes is an attractive biological target of Au(I) and is inhibited by auranofin and auranofin analogs with reasonable potency. Here we employ a combination of experimental and computational investigations into the effect of changes in the phosphine ligand of auranofin on its in vitro inhibition of cathepsin B. Sequential replacement of the ethyl substituents of triethylphosphine by phenyl groups leads to increasing potency in the resultant Au(I) complexes, due in large part to favorable interactions of the more sterically bulky Au(I)-PR3 fragments with the enzyme active site.