Robust Classification of “Relevant” Water Molecules in Putative Protein Binding Sites

Enthalpic Classification of Water Molecules in Target-Ligand Binding

Water molecules play various roles in target-ligand binding. For example, they can be replaced by the ligand and leave the surface of the binding pocket or stay conserved in the interface and form bridges with the target. While experimental techniques supply target-ligand complex structures at an increasing rate, they often have limitations in the measurement of a detailed water structure. Moreover, measurements of binding thermodynamics cannot distinguish between the different roles of individual water molecules. However, such a distinction and classification of the role of individual water molecules would be key to their application in drug design at atomic resolution. In this study, we investigate a quantitative approach for the description of the role of water molecules during ligand binding. Starting from complete hydration structures of the free and ligand-bound target molecules, binding enthalpy scores are calculated for each water molecule using quantum mechanical calculations. A statistical evaluation showed that the scores can distinguish between conserved and displaced classes of water molecules. The classification system was calibrated and tested on more than 1000 individual water positions. The practical tests of the enthalpic classification included important cases of antiviral drug research on HIV-1 protease inhibitors and the Influenza A ion channel. The methodology of classification is based on open source program packages, Gromacs, Mopac, and MobyWat, freely available to the scientific community.

Merging cultures and disciplines to create a drug discovery ecosystem at Virginia commonwealth university: Medicinal chemistry, structural biology, molecular and behavioral pharmacology and computational chemistry

The Department of Medicinal Chemistry, together with the Institute for Structural Biology, Drug Discovery and Development, at Virginia Commonwealth University (VCU) has evolved, organically with quite a bit of bootstrapping, into a unique drug discovery ecosystem in response to the environment and culture of the university and the wider research enterprise. Each faculty member that joined the department and/or institute added a layer of expertise, technology and most importantly, innovation, that fertilized numerous collaborations within the University and with outside partners. Despite moderate institutional support with respect to a typical drug discovery enterprise, the VCU drug discovery ecosystem has built and maintained an impressive array of facilities and instrumentation for drug synthesis, drug characterization, biomolecular structural analysis and biophysical analysis, and pharmacological studies. Altogether, this ecosystem has had major impacts on numerous therapeutic areas, such as neurology, psychiatry, drugs of abuse, cancer, sickle cell disease, coagulopathy, inflammation, aging disorders and others. Novel tools and strategies for drug discovery, design and development have been developed at VCU in the last five decades; e.g., fundamental rational structure-activity relationship (SAR)-based drug design, structure-based drug design, orthosteric and allosteric drug design, design of multi-functional agents towards polypharmacy outcomes, principles on designing glycosaminoglycans as drugs, and computational tools and algorithms for quantitative SAR (QSAR) and understanding the roles of water and the hydrophobic effect.

HINT, a code for understanding the interaction between biomolecules: a tribute to Donald J. Abraham

A long-lasting goal of computational biochemists, medicinal chemists, and structural biologists has been the development of tools capable of deciphering the molecule-molecule interaction code that produces a rich variety of complex biomolecular assemblies comprised of the many different simple and biological molecules of life: water, small metabolites, cofactors, substrates, proteins, DNAs, and RNAs. Software applications that can mimic the interactions amongst all of these species, taking account of the laws of thermodynamics, would help gain information for understanding qualitatively and quantitatively key determinants contributing to the energetics of the bimolecular recognition process. This, in turn, would allow the design of novel compounds that might bind at the intermolecular interface by either preventing or reinforcing the recognition. HINT, hydropathic interaction, was a model and software code developed from a deceptively simple idea of Donald Abraham with the close collaboration with Glen Kellogg at Virginia Commonwealth University. HINT is based on a function that scores atom-atom interaction using LogP, the partition coefficient of any molecule between two phases; here, the solvents are water that mimics the cytoplasm milieu and octanol that mimics the protein internal hydropathic environment. This review summarizes the results of the extensive and successful collaboration between Abraham and Kellogg at VCU and the group at the University of Parma for testing HINT in a variety of different biomolecular interactions, from proteins with ligands to proteins with DNA.

Screening and Biological Evaluation of Soluble Epoxide Hydrolase Inhibitors: Assessing the Role of Hydrophobicity in the Pharmacophore-Guided Search of Novel Hits

The human soluble epoxide hydrolase (sEH) is a bifunctional enzyme that modulates the levels of regulatory epoxy lipids. The hydrolase activity is carried out by a catalytic triad located at the center of a wide L-shaped binding site, which contains two hydrophobic subpockets at both sides. On the basis of these structural features, it can be assumed that desolvation is a major factor in determining the maximal achievable affinity that can be attained for this pocket. Accordingly, hydrophobic descriptors may be better suited to the search of novel hits targeting this enzyme. This study examines the suitability of quantum mechanically derived hydrophobic descriptors in the discovery of novel sEH inhibitors. To this end, three-dimensional quantitative structure-activity relationship (3D-QSAR) pharmacophores were generated by combining electrostatic and steric or alternatively hydrophobic and hydrogen-bond parameters in conjunction with a tailored list of 76 known sEH inhibitors. The pharmacophore models were then validated by using two external sets chosen (i) to rank the potency of four distinct series of compounds and (ii) to discriminate actives from decoys, using in both cases datasets taken from the literature. Finally, a prospective study was performed including a virtual screening of two chemical libraries to identify new potential hits, which were subsequently experimentally tested for their inhibitory activity on human, rat, and mouse sEH. The use of hydrophobic-based descriptors led to the identification of six compounds as inhibitors of the human enzyme with IC50 < 20 nM, including two with IC50 values of 0.4 and 0.7 nM. The results support the use of hydrophobic descriptors as a valuable tool in the search of novel scaffolds that encode a proper hydrophilic/hydrophobic distribution complementary to the target's binding site.

Targeting human thymidylate synthase: Ensemble-based virtual screening for drug repositioning and the role of water

A drug repositioning computational approach was carried to search inhibitors for human thymidylate synthase. An ensemble-based virtual screening of FDA-approved drugs showed the drugs Imatinib, Lumacaftor and Naldemedine to be likely candidates for repurposing. The role of water in the drug-receptor interactions was revealed by the application of an extended AutoDock scoring function that included the water forcefield. The binding affinity scores when hydrated ligands were docked were improved in the drugs considered. Further binding free energy calculations based on the Molecular Mechanics Poisson-Boltzmann Surface Area method revealed that Imatinib, Lumacaftor and Naldemedine scored -130.7 ± 28.1, -210.6 ± 29.9 and -238.0 ± 25.4 kJ/mol, respectively, showing good binding affinity for the candidates considered. Overall, the analysis of the molecular dynamics trajectory of the receptor-drug complexes revealed stable structures for Imatinib, Lumacaftor and Naldemedine, for the entire simulation time.

Water: new aspect of hydrogen bonding in the solid state

All water-water contacts in the crystal structures from the Cambridge Structural Database with d OO ≤ 4.0 Å have been found. These contacts were analysed on the basis of their geometries and interaction energies from CCSD(T)/CBS calculations. The results show 6729 attractive water-water contacts, of which 4717 are classical hydrogen bonds (d OH ≤ 3.0 Å and α ≥ 120°) with most being stronger than -3.3 kcal mol-1. Beyond the region of these hydrogen bonds, there is a large number of attractive interactions (2062). The majority are antiparallel dipolar interactions, where the O-H bonds of two water molecules lying in parallel planes are oriented antiparallel to each other. Developing geometric criteria for these antiparallel dipoles (β1, β2 ≥ 160°, 80 ≤ α ≤ 140° and T HOHO > 40°) yielded 1282 attractive contacts. The interaction energies of these antiparallel oriented water molecules are up to -4.7 kcal mol-1, while most of the contacts have interaction energies in the range -0.9 to -2.1 kcal mol-1. This study suggests that the geometric criteria for defining attractive water-water interactions should be broader than the classical hydrogen-bonding criteria, a change that may reveal undiscovered and unappreciated interactions controlling molecular structure and chemistry.

Combination of pose and rank consensus in docking-based virtual screening: the best of both worlds

The use of high-throughput docking (HTD) in the drug discovery pipeline is today widely established. In spite of methodological improvements in docking accuracy (pose prediction), scoring power, ranking power, and screening power in HTD remain challenging. In fact, pose prediction is of critical importance in view of the pose-dependent scoring process, since incorrect poses will necessarily decrease the ranking power of scoring functions. The combination of results from different docking programs (consensus scoring) has been shown to improve the performance of HTD. Moreover, it has been also shown that a pose consensus approach might also result in database enrichment. We present a new methodology named Pose/Ranking Consensus (PRC) that combines both pose and ranking consensus approaches, to overcome the limitations of each stand-alone strategy. This approach has been developed using four docking programs (ICM, rDock, Auto Dock 4, and PLANTS; the first one is commercial, the other three are free). We undertook a thorough analysis for the best way of combining pose and rank strategies, and applied the PRC to a wide range of 34 targets sampling different protein families and binding site properties. Our approach exhibits an improved systematic performance in terms of enrichment factor and hit rate with respect to either pose consensus or consensus ranking alone strategies at a lower computational cost, while always ensuring the recovery of a suitable number of ligands. An analysis using four free docking programs (replacing ICM by Auto Dock Vina) displayed comparable results.

The new methodology named Pose/Ranking Consensus (PRC) combines both pose and ranking consensus strategies. It displays an enhanced performance in terms of enrichment factor and hit rate, ensuring the recovery of a suitable number of ligands.

Ligand binding: evaluating the contribution of the water molecules network using the Fragment Molecular Orbital method

Water molecules play a crucial role in protein-ligand binding, and many tools exist that aim to predict the position and relative energies of these important, but challenging participants of biomolecular recognition. The available tools are, in general, capable of predicting the location of water molecules. However, predicting the effects of their displacement is still very challenging. In this work, a linear-scaling quantum mechanics-based approach was used to assess water network energetics and the changes in network stability upon ligand structural modifications. This approach offers a valuable way to improve understanding of SAR data and help guide compound design.

Water molecules at protein-drug interfaces: computational prediction and analysis methods

The fundamental importance of water molecules at drug-protein interfaces is now widely recognised and a significant feature in structure-based drug design. Experimental methods for analysing the role of water in drug binding have many challenges, including the accurate location of bound water molecules in crystal structures, and problems in resolving specific water contributions to binding thermodynamics. Computational analyses of binding site water molecules provide an alternative, and in principle complete, structural and thermodynamic picture, and their use is now commonplace in the pharmaceutical industry. In this review, we describe the computational methodologies that are available and discuss their strengths and weaknesses. Additionally, we provide a critical analysis of the experimental data used to validate the methods, regarding the type and quality of experimental structural data. We also discuss some of the fundamental difficulties of each method and suggest directions for future study.

Structural understanding of 5-(4-hydroxy-phenyl)-N-(2-(5-methoxy-1H-indol-3-yl)-ethyl)-3-oxopentanamide as a neuroprotectant for Alzheimer's disease

In our continuing efforts to develop novel neuroprotectants for Alzheimer's disease (AD), a series of analogs based on a lead compound that was recently shown to target the mitochondrial complex I were designed, synthesized and biologically characterized to understand the structure features that are important for neuroprotective activities. The results from a cellular AD model highlighted the important roles of the 4-OH on the phenyl ring and the 5-OCH3 on the indole ring of the lead compound. The results also demonstrated that the β-keto moiety can be modified to retain or improve the neuroprotective activity. Docking studies of selected analogs to the FMN site of mitochondrial complex I also supported the observed neuroprotective activities. Collectively, the results provide further information to guide optimization and development of analogs based on this chemical scaffold as neuroprotectants with a novel mechanism of action for AD.

Nonsteroidal ecdysone receptor agonists use a water channel for binding to the ecdysone receptor complex EcR/USP

The ecdysone receptor (EcR) possesses the remarkable capacity to adapt structurally to different types of ligands. EcR binds ecdysteroids, including 20-hydroxyecdysone (20E), as well as nonsteroidal synthetic agonists such as insecticidal dibenzoylhydrazines (DBHs). Here, we report the crystal structures of the ligand-binding domains of Heliothis virescens EcR/USP bound to the DBH agonist BYI09181 and to the imidazole-type compound BYI08346. The region delineated by helices H7 and H10 opens up to tightly fit a phenyl ring of the ligands to an extent that depends on the bulkiness of ring substituent. In the structure of 20E-bound EcR, this part of the ligand-binding pocket (LBP) contains a channel filled by water molecules that form an intricate hydrogen bond network between 20E and LBP. The water channel present in the nuclear receptor bound to its natural hormone acts as a critical molecular adaptation spring used to accommodate synthetic agonists inside its binding cavity.

GridSolvate: A Web Server for the Prediction of Biomolecular Hydration Properties

We present a novel web server, named gridSolvate, dedicated to the prediction of biomolecular hydration properties. Given a solute in atomic representation, such as a protein or protein–ligand complex, the server determines positions and excess chemical potential of buried and first hydration shell water molecules. Calculations are based on our semiexplicit hydration model that provides computational efficiency close to implicit solvent approaches, yet captures a number of physical effects unique to explicit solvent representation. The model was introduced and validated before in the context of bulk hydration of drug-like solutes and determination of protein hydration sites. Current methodological developments merge those two avenues into a single, easily accessible tool. Here, we focus on the server’s ability to predict water distribution and affinity within protein–ligand interfaces. We demonstrate that with possibly minimal user intervention the server correctly predicts the locations of 77% of interface water molecules in an external set of test structures. The server is freely available at https://gsolvate.biomod.cent.uw.edu.pl.

Identification of 7,8-dihydroxy-3-phenylcoumarin as a reversible monoamine oxidase enzyme inhibitor

We herein report the biological evaluation of 3-arylcoumarin derivatives (3a-l) as potential human monoamine oxidase-A and -B (hMAO-A and hMAO-B) inhibitors. The result indicated that 7,8-dihydroxy-3-(4-nitrophenyl)coumarin (3j) was most effective against MAO-A (inhibition concentration [IC₅₀ ] = 6.46 ± 0.02 µM) and MAO-B (IC₅₀  = 3.8 ± 0.3 µM) enzymes than other synthesized compounds and reference compounds (pargyline and moclobemide). Furthermore, compound (3j) showed (a) nonselectivity against hMAO enzymes, (b) reversible hMAO enzymes inhibition, and (c) neuroprotection against H₂ O₂ -treated human neuroblastoma (N2a) cells. Finally, a molecular modeling study revealed that the hMAO enzymes inhibitory activity of the compound (3j) may be due to the orientation where the nitro (NO₂ ) group lies deep into the receptor and the phenyl ring directed toward flavin adenosine dinucleotide via hydrogen bond interaction, and possible π-π interaction with various important residues. Thus, the results of the present study demonstrate that compound (3j) can be considered as a promising scaffold for the development of hMAO-A and hMAO-B inhibitors.

Characterization and Discovery of a Selective Small-Molecule Modulator of Mitochondrial Complex I Targeting a Unique Binding Site

Mitochondrial dysfunction has been recognized as an essential contributor to many human diseases including neurodegenerative disorders. However, the exact pathological role of mitochondrial dysfunction, especially in mitochondrial reactive oxygen species-associated oxidative stress, remains elusive, partially due to the lack of chemical probes with well-defined mechanisms of action. Herein, we describe the characterization and discovery of a rationally designed small molecule ZCM-I-1 as a selective modulator of the production of reactive oxygen species from mitochondrial complex I that does not alter mitochondrial membrane potential and bioenergetics. Chemical biology studies employing photoaffinity probes derived from ZCM-I-1 demonstrated its novel mechanism of action of modulating complex I via interactions with the flavin mononucleotide site, proximal in the reaction pathway within complex I.

High-Throughput Docking Using Quantum Mechanical Scoring

Today high-throughput docking is one of the most commonly used computational tools in drug lead discovery. While there has been an impressive methodological improvement in docking accuracy, docking scoring still remains an open challenge. Most docking programs are rooted in classical molecular mechanics. However, to better characterize protein-ligand interactions, the use of a more accurate quantum mechanical (QM) description would be necessary. In this work, we introduce a QM-based docking scoring function for high-throughput docking and evaluate it on 10 protein systems belonging to diverse protein families, and with different binding site characteristics. Outstanding results were obtained, with our QM scoring function displaying much higher enrichment (screening power) than a traditional docking method. It is acknowledged that developments in quantum mechanics theory, algorithms and computer hardware throughout the upcoming years will allow semi-empirical (or low-cost) quantum mechanical methods to slowly replace force-field calculations. It is thus urgently needed to develop and validate novel quantum mechanical-based scoring functions for high-throughput docking toward more accurate methods for the identification and optimization of modulators of pharmaceutically relevant targets.

Advances in the Treatment of Explicit Water Molecules in Docking and Binding Free Energy Calculations

Background:

The inclusion of direct effects mediated by water during the ligandreceptor recognition is a hot-topic of modern computational chemistry applied to drug discovery and development. Docking or virtual screening with explicit hydration is still debatable, despite the successful cases that have been presented in the last years. Indeed, how to select the water molecules that will be included in the docking process or how the included waters should be treated remain open questions.

Objective:

In this review, we will discuss some of the most recent methods that can be used in computational drug discovery and drug development when the effect of a single water, or of a small network of interacting waters, needs to be explicitly considered.

Results:

Here, we analyse the software to aid the selection, or to predict the position, of water molecules that are going to be explicitly considered in later docking studies. We also present software and protocols able to efficiently treat flexible water molecules during docking, including examples of applications. Finally, we discuss methods based on molecular dynamics simulations that can be used to integrate docking studies or to reliably and efficiently compute binding energies of ligands in presence of interfacial or bridging water molecules.

Conclusions:

Software applications aiding the design of new drugs that exploit water molecules, either as displaceable residues or as bridges to the receptor, are constantly being developed. Although further validation is needed, workflows that explicitly consider water will probably become a standard for computational drug discovery soon.

Molecular Docking Using Quantum Mechanical-Based Methods

Computational methods are a powerful and consolidated tool in the early stage of the drug lead discovery process. Among these techniques, high-throughput molecular docking has proved to be extremely useful in identifying novel bioactive compounds within large chemical libraries. In the docking procedure, the predominant binding mode of each small molecule within a target binding site is assessed, and a docking score reflective of the likelihood of binding is assigned to them. These methods also shed light on how a given hit could be modified in order to improve protein-ligand interactions and are thus able to guide lead optimization. The possibility of reducing time and cost compared to experimental approaches made this technology highly appealing. Due to methodological developments and the increase of computational power, the application of quantum mechanical methods to study macromolecular systems has gained substantial attention in the last decade. A quantum mechanical description of the interactions involved in molecular association of biomolecules may lead to better accuracy compared to molecular mechanics, since there are many physical phenomena that cannot be correctly described within a classical framework, such as covalent bond formation, polarization effects, charge transfer, bond rearrangements, halogen bonding, and others, that require electrons to be explicitly accounted for. Considering the fact that quantum mechanics-based approaches in biomolecular simulation constitute an active and important field of research, we highlight in this work the recent developments of quantum mechanical-based molecular docking and high-throughput docking.

Prediction of models for ordered solvent in macromolecular structures by a classifier based upon resolution-independent projections of local feature data

Current software tools for the automated building of models for macromolecular X-ray crystal structures are capable of assembling high-quality models for ordered macromolecule and small-molecule scattering components with minimal or no user supervision. Many of these tools also incorporate robust functionality for modelling the ordered water molecules that are found in nearly all macromolecular crystal structures. However, no current tools focus on differentiating these ubiquitous water molecules from other frequently occurring multi-atom solvent species, such as sulfate, or the automated building of models for such species. PeakProbe has been developed specifically to address the need for such a tool. PeakProbe predicts likely solvent models for a given point (termed a `peak') in a structure based on analysis (`probing') of its local electron density and chemical environment. PeakProbe maps a total of 19 resolution-dependent features associated with electron density and two associated with the local chemical environment to a two-dimensional score space that is independent of resolution. Peaks are classified based on the relative frequencies with which four different classes of solvent (including water) are observed within a given region of this score space as determined by large-scale sampling of solvent models in the Protein Data Bank. Designed to classify peaks generated from difference density maxima, PeakProbe also incorporates functionality for identifying peaks associated with model errors or clusters of peaks likely to correspond to multi-atom solvent, and for the validation of existing solvent models using solvent-omit electron-density maps. When tasked with classifying peaks into one of four distinct solvent classes, PeakProbe achieves greater than 99% accuracy for both peaks derived directly from the atomic coordinates of existing solvent models and those based on difference density maxima. While the program is still under development, a fully functional version is publicly available. PeakProbe makes extensive use of cctbx libraries, and requires a PHENIX licence and an up-to-date phenix.python environment for execution.

AMMOS2: a web server for protein-ligand-water complexes refinement via molecular mechanics

AMMOS2 is an interactive web server for efficient computational refinement of protein-small organic molecule complexes. The AMMOS2 protocol employs atomic-level energy minimization of a large number of experimental or modeled protein-ligand complexes. The web server is based on the previously developed standalone software AMMOS (Automatic Molecular Mechanics Optimization for in silico Screening). AMMOS utilizes the physics-based force field AMMP sp4 and performs optimization of protein-ligand interactions at five levels of flexibility of the protein receptor. The new version 2 of AMMOS implemented in the AMMOS2 web server allows the users to include explicit water molecules and individual metal ions in the protein-ligand complexes during minimization. The web server provides comprehensive analysis of computed energies and interactive visualization of refined protein-ligand complexes. The ligands are ranked by the minimized binding energies allowing the users to perform additional analysis for drug discovery or chemical biology projects. The web server has been extensively tested on 21 diverse protein-ligand complexes. AMMOS2 minimization shows consistent improvement over the initial complex structures in terms of minimized protein-ligand binding energies and water positions optimization. The AMMOS2 web server is freely available without any registration requirement at the URL: http://drugmod.rpbs.univ-paris-diderot.fr/ammosHome.php.

Lipophilicity in drug design: an overview of lipophilicity descriptors in 3D-QSAR studies

The pharmacophore concept is a fundamental cornerstone in drug discovery, playing a critical role in determining the success of in silico techniques, such as virtual screening and 3D-QSAR studies. The reliability of these approaches is influenced by the quality of the physicochemical descriptors used to characterize the chemical entities. In this context, a pivotal role is exerted by lipophilicity, which is a major contribution to host-guest interaction and ligand binding affinity. Several approaches have been undertaken to account for the descriptive and predictive capabilities of lipophilicity in 3D-QSAR modeling. Recent efforts encode the use of quantum mechanical-based descriptors derived from continuum solvation models, which open novel avenues for gaining insight into structure-activity relationships studies.

Drug Discovery

Modern chemistry foundations were made in between the 18th and 19th centuries and have been extended in 20th century. R&D towards synthetic chemistry was introduced during the 1960s. Development of new molecular drugs from the herbal plants to synthetic chemistry is the fundamental scientific improvement. About 10-14 years are needed to develop a new molecule with an average cost of more than $800 million. Pharmaceutical industries spend the highest percentage of revenues, but the achievement of desired molecular entities into the market is not increasing proportionately. As a result, an approximate of 0.01% of new molecular entities are approved by the FDA. The highest failure rate is due to inadequate efficacy exhibited in Phase II of the drug discovery and development stage. Innovative technologies such as combinatorial chemistry, DNA sequencing, high-throughput screening, bioinformatics, computational drug design, and computer modeling are now utilized in the drug discovery. These technologies can accelerate the success rates in introducing new molecular entities into the market.

4-Amino-7,8-dihydro-1,6-naphthyridin-5(6 H)-ones as Inhaled Phosphodiesterase Type 4 (PDE4) Inhibitors: Structural Biology and Structure-Activity Relationships

Rational design of a novel template of naphthyridinones rapidly led to PDE4 inhibitors with subnanomolar enzymatic potencies. X-ray crystallography confirmed the binding mode of this novel template. We achieved compounds with double-digit picomolar enzymatic potencies through further structure-based design by targeting both the PDE4 enzyme metal-binding pocket and occupying the solvent-filled pocket. A strategy for lung retention and long duration of action based on low aqueous solubility was followed. In vivo efficacies were measured in a rat lung neutrophilia model by suspension microspray and dry powder administration. Suspension microspray of potent compounds showed in vivo efficacy with a clear dose-response. Despite sustained lung levels, dry powder administration performed much less well and without proper dose-response, highlighting clear differences between the two formulations. This indicates a deficiency in the low aqueous solubility strategy for long duration lung efficacy.

Prediction of Ordered Water Molecules in Protein Binding Sites from Molecular Dynamics Simulations: The Impact of Ligand Binding on Hydration Networks

Water plays a major role in ligand binding and is attracting increasing attention in structure-based drug design. Water molecules can make large contributions to binding affinity by bridging protein–ligand interactions or by being displaced upon complex formation, but these phenomena are challenging to model at the molecular level. Herein, networks of ordered water molecules in protein binding sites were analyzed by clustering of molecular dynamics (MD) simulation trajectories. Locations of ordered waters (hydration sites) were first identified from simulations of high resolution crystal structures of 13 protein–ligand complexes. The MD-derived hydration sites reproduced 73% of the binding site water molecules observed in the crystal structures. If the simulations were repeated without the cocrystallized ligands, a majority (58%) of the crystal waters in the binding sites were still predicted. In addition, comparison of the hydration sites obtained from simulations carried out in the absence of ligands to those identified for the complexes revealed that the networks of ordered water molecules were preserved to a large extent, suggesting that the locations of waters in a protein–ligand interface are mainly dictated by the protein. Analysis of >1000 crystal structures showed that hydration sites bridged protein–ligand interactions in complexes with different ligands, and those with high MD-derived occupancies were more likely to correspond to experimentally observed ordered water molecules. The results demonstrate that ordered water molecules relevant for modeling of protein–ligand complexes can be identified from MD simulations. Our findings could contribute to development of improved methods for structure-based virtual screening and lead optimization.

Effective lead optimization targeting the displacement of bridging receptor–ligand water molecules

Enhancing the binding affinities of ligands by means of lead modifications that displace bridging water molecules at protein–ligand interfaces is an important and widely studied lead optimization strategy.

The Roles of Water in the Protein Matrix: A Largely Untapped Resource for Drug Discovery

The value of thoroughly understanding the thermodynamics specific to a drug discovery/design study is well known. Over the past decade, the crucial roles of water molecules in protein structure, function, and dynamics have also become increasingly appreciated. This Perspective explores water in the biological environment by adopting its point of view in such phenomena. The prevailing thermodynamic models of the past, where water was seen largely in terms of an entropic gain after its displacement by a ligand, are now known to be much too simplistic. We adopt a set of terminology that describes water molecules as being "hot" and "cold", which we have defined as being easy and difficult to displace, respectively. The basis of these designations, which involve both enthalpic and entropic water contributions, are explored in several classes of biomolecules and structural motifs. The hallmarks for characterizing water molecules are examined, and computational tools for evaluating water-centric thermodynamics are reviewed. This Perspective's summary features guidelines for exploiting water molecules in drug discovery.

Waterdock 2.0: Water placement prediction for Holo-structures with a pymol plugin

Water is often found to mediate interactions between a ligand and a protein. It can play a significant role in orientating the ligand within a binding pocket and contribute to the free energy of binding. It would thus be extremely useful to be able to accurately predict the position and orientation of water molecules within a binding pocket. Recently, we developed the WaterDock protocol that was able to predict 97% of the water molecules in a test set. However, this approach generated false positives at a rate of over 20% in most cases and whilst this might be acceptable for some applications, in high throughput scenarios this is not desirable. Here we tackle this problem via the inclusion of knowledge regarding the solvation structure of ligand functional groups. We call this new protocol WaterDock2 and demonstrate that this protocol maintains a similar true positive rate to the original implementation but is capable of reducing the false-positive rate by over 50%. To improve the usability of the method, we have also developed a plugin for the popular graphics program PyMOL. The plugin also contains an implementation of the original WaterDock.

Diving into the Water: Inducible Binding Conformations for BRD4, TAF1(2), BRD9, and CECR2 Bromodomains

The biological role played by non-BET bromodomains remains poorly understood, and it is therefore imperative to identify potent and highly selective inhibitors to effectively explore the biology of individual bromodomain proteins. A ligand-efficient nonselective bromodomain inhibitor was identified from a 6-methyl pyrrolopyridone fragment. Small hydrophobic substituents replacing the N-methyl group were designed directing toward the conserved bromodomain water pocket, and two distinct binding conformations were then observed. The substituents either directly displaced and rearranged the conserved solvent network, as in BRD4(1) and TAF1(2), or induced a narrow hydrophobic channel adjacent to the lipophilic shelf, as in BRD9 and CECR2. The preference of distinct substituents for individual bromodomains provided selectivity handles useful for future lead optimization efforts for selective BRD9, CECR2, and TAF1(2) inhibitors.

A Real-World Perspective on Molecular Design

We present a series of small molecule drug discovery case studies where computational methods were prospectively employed to impact Roche research projects, with the aim of highlighting those methods that provide real added value. Our brief accounts encompass a broad range of methods and techniques applied to a variety of enzymes and receptors. Most of these are based on judicious application of knowledge about molecular conformations and interactions: filling of lipophilic pockets to gain affinity or selectivity, addition of polar substituents, scaffold hopping, transfer of SAR, conformation analysis, and molecular overlays. A case study of sequence-driven focused screening is presented to illustrate how appropriate preprocessing of information enables effective exploitation of prior knowledge. We conclude that qualitative statements enabling chemists to focus on promising regions of chemical space are often more impactful than quantitative prediction.

Water Sites, Networks, And Free Energies with Grand Canonical Monte Carlo

Water molecules play integral roles in the formation of many protein-ligand complexes, and recent computational efforts have been focused on predicting the thermodynamic properties of individual waters and how they may be exploited in rational drug design. However, when water molecules form highly coupled hydrogen-bonding networks, there is, as yet, no method that can rigorously calculate the free energy to bind the entire network or assess the degree of cooperativity between waters. In this work, we report theoretical and methodological developments to the grand canonical Monte Carlo simulation technique. Central to our results is a rigorous equation that can be used to calculate efficiently the binding free energies of water networks of arbitrary size and complexity. Using a single set of simulations, our methods can locate waters, estimate their binding affinities, capture the cooperativity of the water network, and evaluate the hydration free energy of entire protein binding sites. Our techniques have been applied to multiple test systems and compare favorably to thermodynamic integration simulations and experimental data. The implications of these methods in drug design are discussed.

Prediction of Water Binding to Protein Hydration Sites with a Discrete, Semiexplicit Solvent Model

Buried water molecules are ubiquitous in protein structures and are found at the interface of most protein-ligand complexes. Determining their distribution and thermodynamic effect is a challenging yet important task, of great of practical value for the modeling of biomolecular structures and their interactions. In this study, we present a novel method aimed at the prediction of buried water molecules in protein structures and estimation of their binding free energies. It is based on a semiexplicit, discrete solvation model, which we previously introduced in the context of small molecule hydration. The method is applicable to all macromolecular structures described by a standard all-atom force field, and predicts complete solvent distribution within a single run with modest computational cost. We demonstrate that it indicates positions of buried hydration sites, including those filled by more than one water molecule, and accurately differentiates them from sterically accessible to water but void regions. The obtained estimates of water binding free energies are in fair agreement with reference results determined with the double decoupling method.

Rational Design of Benzylidenehydrazinyl-Substituted Thiazole Derivatives as Potent Inhibitors of Human Dihydroorotate Dehydrogenase with in Vivo Anti-arthritic Activity

Human dihydroorotate dehydrogenase (hDHODH) is an attractive therapeutic target for the treatment of rheumatoid arthritis, transplant rejection and other autoimmune diseases. Based on the X-ray structure of hDHODH in complex with lead compound 7, a series of benzylidenehydrazinyl-substituted thiazole derivatives as potent inhibitors of hDHODH were designed and synthesized, of which 19 and 30 were the most potent with IC₅₀ values in the double-digit nanomolar range. Moreover, compound 19 displayed significant anti-arthritic effects and favorable pharmacokinetic profiles in vivo. Further X-ray structure and SAR analyses revealed that the potencies of the designed inhibitors were partly attributable to additional water-mediated hydrogen bond networks formed by an unexpected buried water between hDHODH and the 2-(2-methylenehydrazinyl)thiazole scaffold. This work not only elucidates promising scaffolds targeting hDHODH for the treatment of rheumatoid arthritis, but also demonstrates that the water-mediated hydrogen bond interaction is an important factor in molecular design and optimization.

Molecular docking and structure-based drug design strategies

Pharmaceutical research has successfully incorporated a wealth of molecular modeling methods, within a variety of drug discovery programs, to study complex biological and chemical systems. The integration of computational and experimental strategies has been of great value in the identification and development of novel promising compounds. Broadly used in modern drug design, molecular docking methods explore the ligand conformations adopted within the binding sites of macromolecular targets. This approach also estimates the ligand-receptor binding free energy by evaluating critical phenomena involved in the intermolecular recognition process. Today, as a variety of docking algorithms are available, an understanding of the advantages and limitations of each method is of fundamental importance in the development of effective strategies and the generation of relevant results. The purpose of this review is to examine current molecular docking strategies used in drug discovery and medicinal chemistry, exploring the advances in the field and the role played by the integration of structure- and ligand-based methods.

Malonate-based inhibitors of mammalian serine racemase: kinetic characterization and structure-based computational study

Overactivation of NMDA receptors has been implicated in various neuropathological conditions, including brain ischaemia, neurodegenerative disorders and epilepsy. Production of d-serine, an NMDA receptor co-agonist, from l-serine is catalyzed in vivo by the pyridoxal-5'-phosphate (PLP)-dependent enzyme serine racemase. Specific inhibition of this enzyme has been proposed as a promising strategy for treatment of neurological conditions caused by NMDA receptor dysfunction. Here we present the synthesis and activity analysis of a series of malonate-based inhibitors of mouse serine racemase (mSR). The compounds possessed IC50 values ranging from 40 ± 11 mM for 2,2-bis(hydroxymethyl)malonate down to 57 ± 1 μM for 2,2-dichloromalonate, the most effective competitive mSR inhibitor known to date. The structure-activity relationship of the whole series in the human orthologue (hSR) was interpreted using Glide docking, WaterMap analysis of hydration and quantum mechanical calculations based on the X-ray structure of the hSR/malonate complex. Docking into the hSR active site with three thermodynamically favourable water molecules was able to discern qualitatively between good and weak inhibitors. Further improvement in ranking was obtained using advanced PM6-D3H4X/COSMO semiempirical quantum mechanics-based scoring which distinguished between the compounds with IC50 better/worse than 2 mM. We have thus not only found a new potent hSR inhibitor but also worked out a computer-assisted protocol to rationalize the binding affinity which will thus aid in search for more effective SR inhibitors. Novel, potent hSR inhibitors may represent interesting research tools as well as drug candidates for treatment of diseases associated with NMDA receptor overactivation.

Comparative assessment of scoring functions on an updated benchmark: 1. Compilation of the test set

Scoring functions are often applied in combination with molecular docking methods to predict ligand binding poses and ligand binding affinities or to identify active compounds through virtual screening. An objective benchmark for assessing the performance of current scoring functions is expected to provide practical guidance for the users to make smart choices among available methods. It can also elucidate the common weakness in current methods for future improvements. The primary goal of our comparative assessment of scoring functions (CASF) project is to provide a high-standard, publicly accessible benchmark of this type. Our latest study, i.e., CASF-2013, evaluated 20 popular scoring functions on an updated set of protein-ligand complexes. This data set was selected out of 8302 protein-ligand complexes recorded in the PDBbind database (version 2013) through a fairly complicated process. Sample selection was made by considering the quality of complex structures as well as binding data. Finally, qualified complexes were clustered by 90% similarity in protein sequences. Three representative complexes were chosen from each cluster to control sample redundancy. The final outcome, namely, the PDBbind core set (version 2013), consists of 195 protein-ligand complexes in 65 clusters with binding constants spanning nearly 10 orders of magnitude. In this data set, 82% of the ligand molecules are "druglike" and 78% of the protein molecules are validated or potential drug targets. Correlation between binding constants and several key properties of ligands are discussed. Methods and results of the scoring function evaluation will be described in a companion work in this issue (doi: 10.1021/ci500081m ).

Targeting cystalysin, a virulence factor of treponema denticola-supported periodontitis

Cystalysin from Treponema denticola is a pyridoxal 5'-phosphate dependent lyase that catalyzes the formation of pyruvate, ammonia, and sulfide from cysteine. It is a virulence factor in adult periodontitis because its reaction contributes to hemolysis, which sustains the pathogen. Therefore, it was proposed as a potential antimicrobial target. To identify specific inhibitors by structure-based in silico methods, we first validated the crystal structure of cystalysin as a reliable starting point for the design of ligands. By using single-crystal absorption microspectrophotometry, we found that the enzyme in the crystalline state, with respect to that in solution, exhibits: 1) the same absorption spectra for the catalytic intermediates, 2) a close pKa value for the residue controlling the keto enamine ionization, and 3) similar reactivity with glycine, L-serine, L-methionine, and the nonspecific irreversible inhibitor aminoethoxyvinylglycine. Next, we screened in silico a library of 9357 compounds with the Fingerprints for Ligands and Proteins (FLAP) software, by using the three-dimensional structure of cystalysin as a template. From the library, 17 compounds were selected and experimentally evaluated by enzyme assays and spectroscopic methods. Two compounds were found to competitively inhibit recombinant T. denticola cystalysin, with inhibition constant (Ki ) values of 25 and 37 μM. One of them exhibited a minimum inhibitory concentration (MIC) value of 64 μg mL(-1) on Moraxella catarrhalis ATCC 23246, which proves its ability to cross bacterial membranes.

Intuitive, but not simple: including explicit water molecules in protein-protein docking simulations improves model quality

Characterizing the nature of interaction between proteins that have not been experimentally cocrystallized requires a computational docking approach that can successfully predict the spatial conformation adopted in the complex. In this work, the Hydropathic INTeractions (HINT) force field model was used for scoring docked models in a data set of 30 high-resolution crystallographically characterized "dry" protein-protein complexes and was shown to reliably identify native-like models. However, most current protein-protein docking algorithms fail to explicitly account for water molecules involved in bridging interactions that mediate and stabilize the association of the protein partners, so we used HINT to illuminate the physical and chemical properties of bridging waters and account for their energetic stabilizing contributions. The HINT water Relevance metric identified the "truly" bridging waters at the 30 protein-protein interfaces and we utilized them in "solvated" docking by manually inserting them into the input files for the rigid body ZDOCK program. By accounting for these interfacial waters, a statistically significant improvement of ∼24% in the average hit-count within the top-10 predictions the protein-protein dataset was seen, compared to standard "dry" docking. The results also show scoring improvement, with medium and high accuracy models ranking much better than incorrect ones. These improvements can be attributed to the physical presence of water molecules that alter surface properties and better represent native shape and hydropathic complementarity between interacting partners, with concomitantly more accurate native-like structure predictions.