Prediction of the water content in protein binding sites

HydraProt: A New Deep Learning Tool for Fast and Accurate Prediction of Water Molecule Positions for Protein Structures

Water molecules are integral to the structural stability of proteins and vital for facilitating molecular interactions. However, accurately predicting their precise position around protein structures remains a significant challenge, making it a vibrant research area. In this paper, we introduce HydraProt (deep Hydration of Proteins), a novel methodology for predicting precise positions of water molecule oxygen atoms around protein structures, leveraging two interconnected deep learning architectures: a 3D U-net and a Multi-Layer Perceptron (MLP). Our approach starts by introducing a coarse voxel-based representation of the protein, which allows for rapid sampling of candidate water positions via the 3D U-net. These water positions are then assessed by embedding the water-protein relationship in the Euclidean space by means of an MLP. Finally, a postprocessing step is applied to further refine the MLP predictions. HydraProt surpasses existing state-of-the-art approaches in terms of precision and recall and has been validated on large data sets of protein structures. Notably, our method offers rapid inference runtime and should constitute the method of choice for protein structure studies and drug discovery applications. Our pretrained models, data, and the source code required to reproduce these results are accessible at https://github.com/azamanos/HydraProt.

Evidences for reaction mechanism of 9DB1 DNA catalyst

The first reported reaction mechanism of a DNAzyme, i.e. 9DB1, by using molecular dynamics (MD) simulations includes some ambiguities. We try to overcome some of these ambiguous aspects such as the role of mono and divalent metal ions and observed metal rescue effects by surveying the role of functional groups of original 9DB1 and a variety of its rate conserving and rate decreasing mutations via MD simulations. Conformational differences of these two distinct groups are responsible for their opposite rate trends. Blocking of the OH3' of acceptor nucleotide from effective attack by its hydrogen bond to O4' of donor nucleotide is observed in rate decreasing mutations. Our simulations manifest the role of Na+ and Mg2+ ions in bringing close to each other the ligated atoms. These findings along with observed conformational changes explain carefully the reported metal rescue effects for some phosphate groups.

Molecular docking, simulations of animal peptides against the envelope protein of Dengue virus

Peptides are biologically active, small molecules with high specificity in its mode of action that can be effective at nanomolar concentrations. Peptide-based antiviral medicines have already been licensed and used to treat human immunodeficiency virus (HIV), influenza virus and hepatitis C virus. So far, no peptide drug has been approved for antiviral treatment against Dengue virus, and many are under clinical trials. Therefore, developing a reasonable peptide against the Dengue virus Envelope protein structure will be a successful strategy for treating Dengue. Hence, we investigated protein-protein docking interactions between 215 peptides retrieved from the AVP database against the envelope protein using Cluspro and HADDOCK followed by the evaluation of their allegenicity, toxicity and physicochemical characteristics investigation. Further validation of the protein-peptide complexes was performed with Molecular dynamics simulations and Molecular Mechanics Poisson-Boltzmann surface area (MMPBSA) analysis on the hit inhibitors. This study revealed that Indolicidin (-75.026 ± 1.54 KJ/mol) and Human Neutrophil peptide-1 (-71.6551 ± 2.112 KJ/mol) shows higher negative ΔG binding implicating their relative stabilization in the protein-peptide interactions during 100 ns of dynamic simulations. Also, both the peptides exhibited desirable physicochemical properties and were nonallergenic. Hence, we further aim to test these peptides by in vitro and in vivo studies to confirm their efficacy against Dengue virus.Communicated by Ramaswamy H. Sarma.

HydraMap v.2: Prediction of Hydration Sites and Desolvation Energy with Refined Statistical Potentials

The complex network of water molecules within the binding pocket of a target protein undergoes alterations upon ligand binding, presenting a significant challenge for conventional molecular modeling methods to accurately characterize and compute the associated energy changes. We have previously developed an empirical method, HydraMap (J. Chem. Inf. Model. 2020, 60, 4359-4375), which employs statistical potentials to predict hydration sites and compute desolvation energy, achieving a reasonable balance between accuracy and speed. In this work, we present its improved version, namely, HydraMap v.2. We updated the statistical potentials for protein-water interactions through an analysis of 17 042 crystal protein structures. We also introduced a new feature to evaluate ligand-water interactions by incorporating statistical potentials derived from the solvated structures of 9878 small organic molecules produced by molecular dynamics simulations. By combining these potentials, HydraMap v.2 can predict and compare the hydration sites in a binding pocket before and after ligand binding, identifying key water molecules involved in the binding process, such as those forming bridging hydrogen bonds and unstable ones that can be replaced. We demonstrated the application of HydraMap v.2 in explaining the structure-activity relationship of a panel of MCL-1 inhibitors. The desolvation energies calculated by summing the energy change of each hydration site before and after ligand binding showed good correlation with known ligand binding affinities on six target proteins. In conclusion, HydraMap v.2 offers a cost-effective solution for estimating the desolvation energy during protein-ligand binding and also is practical in guiding lead optimization in structure-based drug discovery.

Recent Developments in Ultralarge and Structure-Based Virtual Screening Approaches

Drug development is a wide scientific field that faces many challenges these days. Among them are extremely high development costs, long development times, and a small number of new drugs that are approved each year. New and innovative technologies are needed to solve these problems that make the drug discovery process of small molecules more time and cost efficient, and that allow previously undruggable receptor classes to be targeted, such as protein-protein interactions. Structure-based virtual screenings (SBVSs) have become a leading contender in this context. In this review, we give an introduction to the foundations of SBVSs and survey their progress in the past few years with a focus on ultralarge virtual screenings (ULVSs). We outline key principles of SBVSs, recent success stories, new screening techniques, available deep learning-based docking methods, and promising future research directions. ULVSs have an enormous potential for the development of new small-molecule drugs and are already starting to transform early-stage drug discovery.

WaterKit: Thermodynamic Profiling of Protein Hydration Sites

Water desolvation is one of the key components of the free energy of binding of small molecules to their receptors. Thus, understanding the energetic balance of solvation and desolvation resulting from individual water molecules can be crucial when estimating ligand binding, especially when evaluating different molecules and poses as done in High-Throughput Virtual Screening (HTVS). Over the most recent decades, several methods were developed to tackle this problem, ranging from fast approximate methods (usually empirical functions using either discrete atom-atom pairwise interactions or continuum solvent models) to more computationally expensive and accurate ones, mostly based on Molecular Dynamics (MD) simulations, such as Grid Inhomogeneous Solvation Theory (GIST) or Double Decoupling. On one hand, MD-based methods are prohibitive to use in HTVS to estimate the role of waters on the fly for each ligand. On the other hand, fast and approximate methods show an unsatisfactory level of accuracy, with low agreement with results obtained with the more expensive methods. Here we introduce WaterKit, a new grid-based sampling method with explicit water molecules to calculate thermodynamic properties using the GIST method. Our results show that the discrete placement of water molecules is successful in reproducing the position of crystallographic waters with very high accuracy, as well as providing thermodynamic estimates with accuracy comparable to more expensive MD simulations. Unlike these methods, WaterKit can be used to analyze specific regions on the protein surface, (such as the binding site of a receptor), without having to hydrate and simulate the whole receptor structure. The results show the feasibility of a general and fast method to compute thermodynamic properties of water molecules, making it well-suited to be integrated in high-throughput pipelines such as molecular docking.

Uncovering water effects in protein-ligand recognition: importance in the second hydration shell and binding kinetics

Developing a ligand with high affinity for a specific protein target is essential for drug design, and water molecules are well known to play a key role in protein-drug recognition. However, predicting the role of particularly ordered water molecules in drug binding remains challenging. Furthermore, hydration free energy contributed from the water network, including the second shell of water molecules, is far from being well studied. In this research we focused on these aspects to accurately and efficiently evaluate water effects in protein-ligand binding affinity. We developed a new strategy using a free-energy calculation method, VM2. We successfully predicted the stable ordered water molecules in a number of protein systems: PDE 10a, HSP90, tryptophan synthase (TRPS), CDK2 and Factor Xa. In some of these, the second shell of water molecules appeared to be critical in protein-ligand binding. We also applied the strategy to largely improve binding free energy calculation using the MM/PBSA method. When applying MM/PBSA alone for two systems, CDK2 and Factor Xa, the computed binding free energy resulted in poor to moderate R² values with experimental data. However, including water free energy correction greatly improved the free energy calculation. Furthermore, our work helped to explain how xk263 is a 1000 times faster binder to HIVp than ritonavir, a potentially useful tool for investigating binding kinetics. Our studies reveal the importance of fully considering water effects in therapeutic developments in pharmaceutical and biotechnology industries and for fundamental research in protein-ligand recognition.

Water Networks in Complexes between Proteins and FDA-Approved Drugs

Water molecules at protein-ligand interfaces are often of significant pharmaceutical interest, owing in part to the entropy which can be released upon the displacement of an ordered water by a therapeutic compound. Protein structures may not, however, completely resolve all critical bound water molecules, or there may be no experimental data available. As such, predicting the location of water molecules in the absence of a crystal structure is important in the context of rational drug design. Grand canonical Monte Carlo (GCMC) is a computational technique that is gaining popularity for the simulation of buried water sites. In this work, we assess the ability of GCMC to accurately predict water binding locations, using a dataset that we have curated, containing 108 unique structures of complexes between proteins and Food and Drug Administration (FDA)-approved small-molecule drugs. We show that GCMC correctly predicts 81.4% of nonbulk crystallographic water sites to within 1.4 Å. However, our analysis demonstrates that the reported performance of water prediction methods is highly sensitive to the way in which the performance is measured. We also find that crystallographic water sites with more protein/ligand hydrogen bonds and stronger electron density are more reliably predicted by GCMC. An analysis of water networks revealed that more than half of the structures contain at least one ligand-contacting water network. In these cases, displacement of a water site by a ligand modification might yield unexpected results if the larger network is destabilized. Cooperative effects between waters should therefore be explicitly considered in structure-based drug design.

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.

A Guide to In Silico Drug Design

The drug discovery process is a rocky path that is full of challenges, with the result that very few candidates progress from hit compound to a commercially available product, often due to factors, such as poor binding affinity, off-target effects, or physicochemical properties, such as solubility or stability. This process is further complicated by high research and development costs and time requirements. It is thus important to optimise every step of the process in order to maximise the chances of success. As a result of the recent advancements in computer power and technology, computer-aided drug design (CADD) has become an integral part of modern drug discovery to guide and accelerate the process. In this review, we present an overview of the important CADD methods and applications, such as in silico structure prediction, refinement, modelling and target validation, that are commonly used in this area.

Absolute Binding Free Energy Calculations for Buried Water Molecules

Water often plays a key role in mediating protein-ligand interactions. Understanding contributions from active-site water molecules to binding thermodynamics of a ligand is important in predicting binding free energies for ligand optimization. In this work, we tested a non-equilibrium switching method for absolute binding free energy calculations on water molecules in binding sites of 13 systems. We discuss the lessons we learned about identified issues that affected our calculations and ways to address them. This work fits with our larger focus on how to do accurate ligand binding free energy calculations when water rearrangements are very slow, such as rearrangements due to ligand modification (as in relative free energy calculations) or ligand binding (as in absolute free energy calculations). The method studied in this work can potentially be used to account for limited water sampling via providing endpoint corrections to free energy calculations using our calculated binding free energy of water.

Predicting Conserved Water Molecules in Binding Sites of Proteins Using Machine Learning Methods and Combining Features

Water molecules play an important role in many biological processes in terms of stabilizing protein structures, assisting protein folding, and improving binding affinity. It is well known that, due to the impacts of various environmental factors, it is difficult to identify the conserved water molecules (CWMs) from free water molecules (FWMs) directly as CWMs are normally deeply embedded in proteins and form strong hydrogen bonds with surrounding polar groups. To circumvent this difficulty, in this work, the abundance of spatial structure information and physicochemical properties of water molecules in proteins inspires us to adopt machine learning methods for identifying the CWMs. Therefore, in this study, a machine learning framework to identify the CWMs in the binding sites of the proteins was presented. First, by analyzing water molecules' physicochemical properties and spatial structure information, six features (i.e., atom density, hydrophilicity, hydrophobicity, solvent-accessible surface area, temperature B-factors, and mobility) were extracted. Those features were further analyzed and combined to reach a higher CWM identification rate. As a result, an optimal feature combination was determined. Based on this optimal combination, seven different machine learning models (including support vector machine (SVM), K-nearest neighbor (KNN), decision tree (DT), logistic regression (LR), discriminant analysis (DA), naïve Bayes (NB), and ensemble learning (EL)) were evaluated for their abilities in identifying two categories of water molecules, i.e., CWMs and FWMs. It showed that the EL model was the desired prediction model due to its comprehensive advantages. Furthermore, the presented methodology was validated through a case study of crystal 3skh and extensively compared with Dowser++. The prediction performance showed that the optimal feature combination and the desired EL model in our method could achieve satisfactory prediction accuracy in identifying CWMs from FWMs in the proteins' binding sites.

Predicting Protein-Ligand Docking Structure with Graph Neural Network

Modern day drug discovery is extremely expensive and time consuming. Although computational approaches help accelerate and decrease the cost of drug discovery, existing computational software packages for docking-based drug discovery suffer from both low accuracy and high latency. A few recent machine learning-based approaches have been proposed for virtual screening by improving the ability to evaluate protein-ligand binding affinity, but such methods rely heavily on conventional docking software to sample docking poses, which results in excessive execution latencies. Here, we propose and evaluate a novel graph neural network (GNN)-based framework, MedusaGraph, which includes both pose-prediction (sampling) and pose-selection (scoring) models. Unlike the previous machine learning-centric studies, MedusaGraph generates the docking poses directly and achieves from 10 to 100 times speedup compared to state-of-the-art approaches, while having a slightly better docking accuracy.

Enhancing Sampling of Water Rehydration on Ligand Binding: A Comparison of Techniques

Water often plays a key role in protein structure, molecular recognition, and mediating protein-ligand interactions. Thus, free energy calculations must adequately sample water motions, which often proves challenging in typical MD simulation time scales. Thus, the accuracy of methods relying on MD simulations ends up limited by slow water sampling. Particularly, as a ligand is removed or modified, bulk water may not have time to fill or rearrange in the binding site. In this work, we focus on several molecular dynamics (MD) simulation-based methods attempting to help rehydrate buried water sites: BLUES, using nonequilibrium candidate Monte Carlo (NCMC); grand, using grand canonical Monte Carlo (GCMC); and normal MD. We assess the accuracy and efficiency of these methods in rehydrating target water sites. We selected a range of systems with varying numbers of waters in the binding site, as well as those where water occupancy is coupled to the identity or binding mode of the ligand. We analyzed the rehydration of buried water sites in binding pockets using both clustering of trajectories and direct analysis of electron density maps. Our results suggest both BLUES and grand enhance water sampling relative to normal MD and grand is more robust than BLUES, but also that water sampling remains a major challenge for all of the methods tested. The lessons we learned for these methods and systems are discussed.

The oxytocin signaling complex reveals a molecular switch for cation dependence

Oxytocin (OT) and vasopressin (AVP) are conserved peptide signaling hormones that are critical for diverse processes including osmotic homeostasis, reproduction, lactation and social interaction. OT acts through the oxytocin receptor (OTR), a magnesium-dependent G protein-coupled receptor that is a therapeutic target for treatment of postpartum hemorrhage, dysfunctional labor and autism. However, the molecular mechanisms that underlie OTR activation by OT and the dependence on magnesium remain unknown. Here we present the wild-type active-state structure of human OTR bound to OT and miniG_q/i determined by cryo-EM. The structure reveals a unique activation mechanism adopted by OTR involving both the formation of a Mg²⁺ coordination complex between OT and the receptor, and disruption of transmembrane helix 7 (TM7) by OT. Our functional assays demonstrate the role of TM7 disruption and provide the mechanism of full agonism by OT and partial agonism by OT analogs. Furthermore, we find that the identity of a single cation-coordinating residue across vasopressin family receptors determines whether the receptor is cation-dependent. Collectively, these results demonstrate how the Mg²⁺-dependent OTR is activated by OT, provide essential information for structure-based drug discovery efforts and shed light on the molecular determinants of cation dependence of vasopressin family receptors throughout the animal kingdom.

Plasticity in ligand recognition at somatostatin receptors

Somatostatin is a signaling peptide that plays a pivotal role in physiologic processes relating to metabolism and growth through its actions at somatostatin receptors (SSTRs). Members of the SSTR subfamily, particularly SSTR2, are key drug targets for neuroendocrine neoplasms, with synthetic peptide agonists currently in clinical use. Here, we show the cryogenic-electron microscopy structures of active-state SSTR2 in complex with heterotrimeric Gi3 and either the endogenous ligand SST14 or the FDA-approved drug octreotide. Complemented by biochemical assays and molecular dynamics simulations, these structures reveal key details of ligand recognition and receptor activation at SSTRs. We find that SSTR ligand recognition is highly diverse, as demonstrated by ligand-induced conformational changes in ECL2 and substantial sequence divergence across subtypes in extracellular regions. Despite this complexity, we rationalize several known sources of SSTR subtype selectivity and identify an additional interaction for specific binding. These results provide valuable insights for structure-based drug discovery at SSTRs.

Fumarato and Phthalato Bridged Dinuclear Metal-Organic Cu(II) and Mn(II) Compounds involving Infinite Fumarate-water Assemblies and Unusual Structure-guiding H-bonded Synthons: Antiproliferative Evaluation and Theoretical Studies

Two new dinuclear coordination compounds viz. [Cu2(µ-fum)(phen)2(H2O)6](fum)•6H2O (1) and [Mn2(µ-phth)2(phen)4]•2H2O (2) (phen = 1,10-phenanthroline, fum = fumarate and phth = phthalate) have been synthesized and characterized by elemental analysis, single...

Accurate Prediction of Hydration Sites of Proteins Using Energy Model With Atom Embedding

We propose a method based on neural networks to accurately predict hydration sites in proteins. In our approach, high-quality data of protein structures are used to parametrize our neural network model, which is a differentiable score function that can evaluate an arbitrary position in 3D structures on proteins and predict the nearest water molecule that is not present. The score function is further integrated into our water placement algorithm to generate explicit hydration sites. In experiments on the OppA protein dataset used in previous studies and our selection of protein structures, our method achieves the highest model quality in terms of F1 score, compared to several previous studies.

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.

G-protein activation by a metabotropic glutamate receptor

Family C G-protein-coupled receptors (GPCRs) operate as obligate dimers with extracellular domains that recognize small ligands, leading to G-protein activation on the transmembrane (TM) domains of these receptors by an unknown mechanism¹. Here we show structures of homodimers of the family C metabotropic glutamate receptor 2 (mGlu2) in distinct functional states and in complex with heterotrimeric G_i. Upon activation of the extracellular domain, the two transmembrane domains undergo extensive rearrangement in relative orientation to establish an asymmetric TM6-TM6 interface that promotes conformational changes in the cytoplasmic domain of one protomer. Nucleotide-bound G_i can be observed pre-coupled to inactive mGlu2, but its transition to the nucleotide-free form seems to depend on establishing the active-state TM6-TM6 interface. In contrast to family A and B GPCRs, G-protein coupling does not involve the cytoplasmic opening of TM6 but is facilitated through the coordination of intracellular loops 2 and 3, as well as a critical contribution from the C terminus of the receptor. The findings highlight the synergy of global and local conformational transitions to facilitate a new mode of G-protein activation.

Mapping Water Thermodynamics on Drug Candidates via Molecular Building Blocks: a Strategy to Improve Ligand Design and Rationalize SAR

The consideration of interactions involving water molecules in protein-ligand binding is widely appreciated in drug discovery nowadays. However, it is not ultimately clear how insights about these interactions translate into molecular design concepts. In this work, we introduce a computational strategy that, trained with high-precision experimental data, allows for the decomposition of water-related thermodynamic properties into chemically relevant building blocks (BBs) of a given ligand scaffold. For each of these BBs, a score based on solvation energy and entropy is computed, thus enabling the analysis of solvent-related affinity contributions for individual BBs. We find the nonvariable BB in a congeneric ligand pair to have a larger impact on the binding affinity than the variable part thus suggesting strong cooperative effects. Furthermore, we find enhanced solute-solvent interactions for a BB due to the presence of a C-F bond. Our investigation may be used to design drug molecules with tailored solvent thermodynamic properties.

Outliers in SAR and QSAR: 3. Importance of considering the role of water molecules in protein-ligand interactions and quantitative structure-activity relationship studies

It is frequently mentioned that QSARs have not generally lived up to expectations, especially in cases where high predictability is expected yet failed to deliver satisfactory results. Even though outliers can provide an increased opportunity in drug discovery research, outliers in SAR and QSAR can contort predictions and affect the accuracy if proper attention is not given. The percentages of outliers in QSARs have not changed appreciably over the last decade. In our previous studies, we suggested two possible sources of outliers in SAR and QSAR. In this paper, we suggest an additional possible source of outliers in QSAR. We presented several literature examples that show one or more water molecules that play a critical role in protein-ligand binding interactions as observed in their crystal structures. These examples illustrate that failing to account for the effects of water molecules in protein-ligand interactions could mislead interpretation and possibly yield outliers in SAR and QSAR. Examples include cases where QSAR, considering the role of water molecules in protein-ligand crystal structures, provided deeper insight into the understanding and interpretation of the developed QSAR.

Enhancing water sampling of buried binding sites using nonequilibrium candidate Monte Carlo

Water molecules can be found interacting with the surface and within cavities in proteins. However, water exchange between bulk and buried hydration sites can be slow compared to simulation timescales, thus leading to the inefficient sampling of the locations of water. This can pose problems for free energy calculations for computer-aided drug design. Here, we apply a hybrid method that combines nonequilibrium candidate Monte Carlo (NCMC) simulations and molecular dynamics (MD) to enhance sampling of water in specific areas of a system, such as the binding site of a protein. Our approach uses NCMC to gradually remove interactions between a selected water molecule and its environment, then translates the water to a new region, before turning the interactions back on. This approach of gradual removal of interactions, followed by a move and then reintroduction of interactions, allows the environment to relax in response to the proposed water translation, improving acceptance of moves and thereby accelerating water exchange and sampling. We validate this approach on several test systems including the ligand-bound MUP-1 and HSP90 proteins with buried crystallographic waters removed. We show that our BLUES (NCMC/MD) method enhances water sampling relative to normal MD when applied to these systems. Thus, this approach provides a strategy to improve water sampling in molecular simulations which may be useful in practical applications in drug discovery and biomolecular design.

Biologically relevant and energetically significant cooperative ternary (π–π)2/(π–π)1/(π–π)2 assemblies and fascinating discrete (H2O)21 clusters in isostructural 2,5-pyridine dicarboxylato Co(ii) and Zn(ii) phenanthroline compounds: antiproliferative evaluation and theoretical studies

Cytotoxicity in cancer cells with structure activity relationship has been explored in isostructural Co(ii) and Zn(ii) compounds involving energetically significant cooperative (π–π)₂/(π–π)₁/(π–π)₂ assemblies and fascinating (H₂O)₂₁ clusters.

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.

Thermodynamics and Kinetics of Drug-Target Binding by Molecular Simulation

Computational studies play an increasingly important role in chemistry and biophysics, mainly thanks to improvements in hardware and algorithms. In drug discovery and development, computational studies can reduce the costs and risks of bringing a new medicine to market. Computational simulations are mainly used to optimize promising new compounds by estimating their binding affinity to proteins. This is challenging due to the complexity of the simulated system. To assess the present and future value of simulation for drug discovery, we review key applications of advanced methods for sampling complex free-energy landscapes at near nonergodicity conditions and for estimating the rate coefficients of very slow processes of pharmacological interest. We outline the statistical mechanics and computational background behind this research, including methods such as steered molecular dynamics and metadynamics. We review recent applications to pharmacology and drug discovery and discuss possible guidelines for the practitioner. Recent trends in machine learning are also briefly discussed. Thanks to the rapid development of methods for characterizing and quantifying rare events, simulation's role in drug discovery is likely to expand, making it a valuable complement to experimental and clinical approaches.

aquaPELE: A Monte Carlo-Based Algorithm to Sample the Effects of Buried Water Molecules in Proteins

Water is frequently found inside proteins, carrying out important roles in catalytic reactions or molecular recognition tasks. Therefore, computational models that aim to study protein-ligand interactions usually have to include water effects through explicit or implicit approaches to obtain reliable results. While full explicit models might be too computationally daunting for some applications, implicit models are normally faster but omit some of the most important contributions of water. This is the case of our in-house software, called protein energy landscape exploration (PELE), which uses implicit models to speed up conformational explorations as much as possible; the lack of explicit water sampling, however, limits its model. In this work, we confront this problem with the development of aquaPELE. It is a new algorithm that extends the exploration capabilities while keeping efficiency as it employs a mixed implicit/explicit approach to also take into account the effects of buried water molecules. With an additional Monte Carlo (MC) routine, a set of explicit water molecules is perturbed inside protein cavities and their effects are dynamically adjusted to the current state of the system. As a result, this implementation can be used to predict the principal hydration sites or the rearrangement and displacement of conserved water molecules upon the binding of a ligand. We benchmarked this new tool focusing on estimating ligand binding modes and hydration sites in cavities with important interfacial water molecules, according to crystallographic structures. Results suggest that aquaPELE sets a fast and reliable alternative for molecular recognition studies in systems with a strong water-dependency.

Alchemical Binding Free Energy Calculations in AMBER20: Advances and Best Practices for Drug Discovery

Predicting protein-ligand binding affinities and the associated thermodynamics of biomolecular recognition is a primary objective of structure-based drug design. Alchemical free energy simulations offer a highly accurate and computationally efficient route to achieving this goal. While the AMBER molecular dynamics package has successfully been used for alchemical free energy simulations in academic research groups for decades, widespread impact in industrial drug discovery settings has been minimal because of the previous limitations within the AMBER alchemical code, coupled with challenges in system setup and postprocessing workflows. Through a close academia-industry collaboration we have addressed many of the previous limitations with an aim to improve accuracy, efficiency, and robustness of alchemical binding free energy simulations in industrial drug discovery applications. Here, we highlight some of the recent advances in AMBER20 with a focus on alchemical binding free energy (BFE) calculations, which are less computationally intensive than alternative binding free energy methods where full binding/unbinding paths are explored. In addition to scientific and technical advances in AMBER20, we also describe the essential practical aspects associated with running relative alchemical BFE calculations, along with recommendations for best practices, highlighting the importance not only of the alchemical simulation code but also the auxiliary functionalities and expertise required to obtain accurate and reliable results. This work is intended to provide a contemporary overview of the scientific, technical, and practical issues associated with running relative BFE simulations in AMBER20, with a focus on real-world drug discovery applications.

Guiding Conventional Protein-Ligand Docking Software with Convolutional Neural Networks

The high-performance computational techniques have brought significant benefits for drug discovery efforts in recent decades. One of the most challenging problems in drug discovery is the protein-ligand binding pose prediction. To predict the most stable structure of the complex, the performance of conventional structure-based molecular docking methods heavily depends on the accuracy of scoring or energy functions (as an approximation of affinity) for each pose of the protein-ligand docking complex to effectively guide the search in an exponentially large solution space. However, due to the heterogeneity of molecular structures, the existing scoring calculation methods are either tailored to a particular data set or fail to exhibit high accuracy. In this paper, we propose a convolutional neural network (CNN)-based model that learns to predict the stability factor of the protein-ligand complex and exhibits the ability of CNNs to improve the existing docking software. Evaluated results on PDBbind data set indicate that our approach reduces the execution time of the traditional docking-based method while improving the accuracy. Our code, experiment scripts, and pretrained models are available at https://github.com/j9650/MedusaNet.

Comparing Fragment Binding Poses Prediction Using HSP90 as a Key Study: When Bound Water Makes the Difference

Fragment-Based Drug Discovery (FBDD) approaches have gained popularitynot only in industry but also in academic research institutes. However, the computational prediction of the binding mode adopted by fragment-like molecules within a protein binding site is still a very challenging task. One of the most crucial aspects of fragment binding is related to the large amounts of bound waters in the targeted binding pocket. The binding affinity of fragmentsmay not be sufficientto displace the bound water molecules. In the present work, we confirmed the importance of the bound water molecules in the correct prediction of the fragment binding mode. Moreover, we investigate whether the use of methods based on explicit solvent molecular dynamics simulations can improve the accuracy of fragment posing. The protein chosen for this study is HSP-90.

Enhancing Water Sampling in Free Energy Calculations with Grand Canonical Monte Carlo

The prediction of protein-ligand binding affinities using free energy perturbation (FEP) is becoming increasingly routine in structure-based drug discovery. Most FEP packages use molecular dynamics (MD) to sample the configurations of proteins and ligands, as MD is well-suited to capturing coupled motion. However, MD can be prohibitively inefficient at sampling water molecules that are buried within binding sites, which has severely limited the domain of applicability of FEP and its prospective usage in drug discovery. In this paper, we present an advancement of FEP that augments MD with grand canonical Monte Carlo (GCMC), an enhanced sampling method, to overcome the problem of sampling water. We accomplished this without degrading computational performance. On both old and newly assembled data sets of protein-ligand complexes, we show that the use of GCMC in FEP is essential for accurate and robust predictions for ligand perturbations that disrupt buried water.

Metadynamics as a Postprocessing Method for Virtual Screening with Application to the Pseudokinase Domain of JAK2

With standard scoring methods, top-ranked compounds from virtual screening by docking often turn out to be inactive. For this reason, metadynamics, a method used to sample rare events, was studied to further evaluate docking poses with the aim of reducing false positives. Specifically, virtual screening was performed with Glide SP to seek potential molecules to bind to the ATP site in the pseudokinase domain of JAK2 kinase, and promising compounds were selected from the top-ranked 1000 based on visualization. Rescoring with Glide XP, GOLD, and MM/GBSA was unable to differentiate well between active and inactive compounds. Metadynamics was then used to gauge the relative binding affinity from the required time or the potential of mean force needed to dissociate the ligand from the bound complex. With consideration of previously known binders of varying affinities, metadynamics was able to differentiate between the most active compounds and inactive or weakly active ones, and it could identify correctly most of the selected virtual screening compounds as false positives. Thus, metadynamics has the potential to be a viable postprocessing method for virtual screening, minimizing the expense of buying or synthesizing inactive compounds.

An assessment of water placement algorithms in quantum mechanics/molecular mechanics modeling: the case of rhodopsins' first spectral absorption band maxima

Quantum mechanics/molecular mechanics (QM/MM) models are a widely used tool to obtain detailed insight into the properties and functioning of proteins. The outcome of QM/MM studies heavily depends on the quality of the applied QM/MM model. Prediction and right placement of internal water molecules in protein cavities is one of the critical parts of any QM/MM model construction. Herein, we performed a systematic study of four protein hydration algorithms. We tested these algorithms for their ability to predict X-ray-resolved water molecules for a set of membrane photosensitive rhodopsin proteins, as well as the influence of the applied water placement algorithms on the QM/MM calculated absorption maxima (λ_max) of these proteins. We used 49 rhodopsins and their intermediates with available X-ray structures as the test set. We found that a proper choice of hydration algorithms and setups is needed to predict functionally important water molecules in the chromophore-binding cavity of rhodopsins, such as the water cluster in the N-H region of bacteriorhodopsin or two water molecules in the binding pocket of bovine visual rhodopsin. The QM/MM calculated λ_max of rhodopsins is also quite sensitive to the applied protein hydration protocols. The best methodology allows obtaining an 18.0 nm average value for the absolute deviation of the calculated λ_max from the experimental λ_max. Although the major effect of water molecules on λ_max originates from the water molecules located in the binding pocket, the water molecules outside the binding pocket also affect the calculated λ_max mainly by causing a reorganization of the protein structure. The results reported in this study can be used for the evaluation and further development of hydration methodologies, in general, and rhodopsin QM/MM models, in particular.

GemSpot: A Pipeline for Robust Modeling of Ligands into Cryo-EM Maps

Producing an accurate atomic model of biomolecule-ligand interactions from maps generated by cryoelectron microscopy (cryo-EM) often presents challenges inherent to the methodology and the dynamic nature of ligand binding. Here, we present GemSpot, an automated pipeline of computational chemistry methods that take into account EM map potentials, quantum mechanics energy calculations, and water molecule site prediction to generate candidate poses and provide a measure of the degree of confidence. The pipeline is validated through several published cryo-EM structures of complexes in different resolution ranges and various types of ligands. In all cases, at least one identified pose produced both excellent interactions with the target and agreement with the map. GemSpot will be valuable for the robust identification of ligand poses and drug discovery efforts through cryo-EM.

Applications of water molecules for analysis of macromolecule properties

Water molecules maintain proteins' structures, functions, stabilities and dynamics. They can occupy certain positions or pass quickly via a protein's interior. Regardless of their behaviour, water molecules can be used for the analysis of proteins' structural features and biochemical properties. Here, we present a list of several software programs that use the information provided by water molecules to: i) analyse protein structures and provide the optimal positions of water molecules for protein hydration, ii) identify high-occupancy water sites in order to analyse ligand binding modes, and iii) detect and describe tunnels and cavities. The analysis of water molecules' distribution and trajectories sheds a light on proteins' interactions with small molecules, on the dynamics of tunnels and cavities, on protein composition and also on the functionality, transportation network and location of functionally relevant residues. Finally, the correct placement of water molecules in protein crystal structures can significantly improve the reliability of molecular dynamics simulations.

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.

A machine learning based intramolecular potential for a flexible organic molecule

Here, we employ the kernel regression machine learning technique to construct an analytical potential that reproduces the quantum mechanical potential energy surface of a small, flexible, drug-like molecule, 3-(benzyloxy)pyridin-2-amine.

Modelling flexible protein–ligand binding in p38α MAP kinase using the QUBE force field

The accuracy of quantum mechanical bespoke (QUBE) force fields for protein–ligand binding free energy calculations are benchmarked against experiment.

Adipato bridged novel hexanuclear Cu(ii) and polymeric Co(ii) coordination compounds involving cooperative supramolecular assemblies and encapsulated guest water clusters in a square grid host: antiproliferative evaluation and theoretical studies

Novel adipato bridged Cu(ii) and Co(ii) complexes synthesized by considering cytotoxicity, apoptosis, ROS generation, molecular docking and pharmacophore features.

The In Silico Drug Discovery Toolbox: Applications in Lead Discovery and Optimization

BACKGROUND

Discovery and development of a new drug is a long lasting and expensive journey that takes around 20 years from starting idea to approval and marketing of new medication. Despite R&D expenditures have been constantly increasing in the last few years, the number of new drugs introduced into market has been steadily declining. This is mainly due to preclinical and clinical safety issues, which still represent about 40% of drug discontinuation. To cope with this issue, a number of in silico techniques are currently being used for an early stage evaluation/prediction of potential safety issues, allowing to increase the drug-discovery success rate and reduce costs associated with the development of a new drug.

METHODS

In the present review, we will analyse the early steps of the drug-discovery pipeline, describing the sequence of steps from disease selection to lead optimization and focusing on the most common in silico tools used to assess attrition risks and build a mitigation plan.

RESULTS

A comprehensive list of widely used in silico tools, databases, and public initiatives that can be effectively implemented and used in the drug discovery pipeline has been provided. A few examples of how these tools can be problem-solving and how they may increase the success rate of a drug discovery and development program have been also provided. Finally, selected examples where the application of in silico tools had effectively contributed to the development of marketed drugs or clinical candidates will be given.

CONCLUSION

The in silico toolbox finds great application in every step of early drug discovery: (i) target identification and validation; (ii) hit identification; (iii) hit-to-lead; and (iv) lead optimization. Each of these steps has been described in details, providing a useful overview on the role played by in silico tools in the decision-making process to speed-up the discovery of new drugs.

Inclusion of enclosed hydration effects in the binding free energy estimation of dopamine D3 receptor complexes

Confined hydration and conformational flexibility are some of the challenges encountered for the rational design of selective antagonists of G-protein coupled receptors. We present a set of C3-substituted (-)-stepholidine derivatives as potent binders of the dopamine D3 receptor. The compounds are characterized biochemically, as well as by computer modeling using a novel molecular dynamics-based alchemical binding free energy approach which incorporates the effect of the displacement of enclosed water molecules from the binding site. The free energy of displacement of specific hydration sites is obtained using the Hydration Site Analysis method with explicit solvation. This work underscores the critical role of confined hydration and conformational reorganization in the molecular recognition mechanism of dopamine receptors and illustrates the potential of binding free energy models to represent these key phenomena.

Discovery of Orexant and Anorexant Agents with Indazole Scaffold Endowed with Peripheral Antiedema Activity

The endocannabinoid system represents an integrated neuronal network involved in the control of several organisms' functions, such as feeding behavior. A series of hybrids of 5-(4-chlorophenyl)-1-(2,4-dichloro-phenyl)-4-methyl-N-(piperidin-1-yl)-1H-pyrazole-3-carboxamide (mimonabant), a well-known inverse agonist of the type-1 cannabinoid receptor (CB1), once used as an antiobesity drug, and the N-(2S)-substitutes of 1-[(4-fluorophenyl)methyl]indazole-3-carboxamide with 1-amino-3-methyl-1-oxobutane (AB-Fubinaca), 1-amino-3,3-dimethyl-1-oxobutane (ADB-Fubinaca), and 3-methylbutanoate (AMB-Fubinaca), endowed with potent agonistic activity towards cannabinoid receptors CB1 and CB2 were in solution as C-terminal amides, acids, methyl esters and N-methyl amides. These compounds have been studied by binding assays to cannabinoid receptors and by functional receptor assays, using rat brain membranes in vitro. The most active among them as an agonist, (S)-1-(2,4-dichlorobenzyl)-N-(3,3-dimethyl-1-(methylamino)-1-oxobutan-2-yl)-1H-indazole-3-carboxamide (LONI11), and an antagonist, (S)-2-(1-(2,4-dichlorobenzyl)-1H-indazole-3-carboxamido)-3-methylbutanoic acid (LONI4), were tested in vivo in mic, to evaluate their ability to stimulate or suppress feeding behavior after intraperitoneal (i.p.) administration. For a LONI11 formalin test and a tail flick test after an administration by the subcutaneous (s.c.) and intracerebroventricular (i.c.v.) routes, respectively, were also carried out in vivo in mice to investigate the antinociceptive property at the central and peripheral levesl. We observed a significant orexant effect for LONI11 and an intense anorexant effect for (S)-methyl 2-(1-(2,4-dichlorobenzyl)-1H-indazole-3-carboxamido)-3,3-dimethylbutanoate (LONI2) and LONI4. In zymosan-induced edema and hyperalgesia, LONI11 reduced the percent of paw volume increase and paw latency after s.c. administration, also suggesting a possible peripheral anti-inflammatory activity.

Time-Domain Analysis of Molecular Dynamics Trajectories Using Deep Neural Networks: Application to Activity Ranking of Tankyrase Inhibitors

Molecular dynamics simulations provide valuable insights into the behavior of molecular systems. Extending the recent trend of using machine learning techniques to predict physicochemical properties from molecular dynamics data, we propose to consider the trajectories as multidimensional time series represented by 2D tensors containing the ligand-protein interaction descriptor values for each time step. Similar in structure to the time series encountered in modern approaches for signal, speech, and natural language processing, these time series can be directly analyzed using long short-term memory (LSTM) recurrent neural networks or convolutional neural networks (CNNs). The predictive regression models for the ligand-protein affinity were built for a subset of the PDBbind v.2017 database and applied to inhibitors of tankyrase, an enzyme of the poly(ADP-ribose)-polymerase (PARP) family that can be used in the treatment of colorectal cancer. As an additional test set, a subset of the Community Structure-Activity Resource (CSAR) data set was used. For comparison, the random forest and simple neural network models based on the crystal pose or the trajectory-averaged descriptors were used, as well as the commonly employed docking and molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) scores. Convolutional neural networks based on the 2D tensors of ligand-protein interaction descriptors for short (2 ns) trajectories provide the best accuracy and predictive power, reaching the Spearman rank correlation coefficient of 0.73 and Pearson correlation coefficient of 0.70 for the tankyrase test set. Taking into account the recent increase in computational power of modern GPUs and relatively low computational complexity of the proposed approach, it can be used as an advanced virtual screening filter for compound prioritization.