Analysis of biomolecular solvation sites by 3D-RISM theory

Transformation of peptides to small molecules in medicinal chemistry: Challenges and opportunities

Peptides are native binders involved in numerous physiological life procedures, such as cellular signaling, and serve as ready-made regulators of biochemical processes. Meanwhile, small molecules compose many drugs owing to their outstanding advantages of physiochemical properties and synthetic convenience. A novel field of research is converting peptides into small molecules, providing a convenient portable solution for drug design or peptidomic research. Endowing properties of peptides onto small molecules can evolutionarily combine the advantages of both moieties and improve the biological druggability of molecules. Herein, we present eight representative recent cases in this conversion and elaborate on the transformation process of each case. We discuss the innovative technological methods and research approaches involved, and analyze the applicability conditions of the approaches and methods in each case, guiding further modifications of peptides to small molecules. Finally, based on the aforementioned cases, we summarize a general procedure for peptide-to-small molecule modifications, listing the technological methods available for each transformation step and providing our insights on the applicable scenarios for these methods. This review aims to present the progress of peptide-to-small molecule modifications and propose our thoughts and perspectives for future research in this field.

RISMiCal: A software package to perform fast RISM/3D-RISM calculations

Solvent plays an essential role in a variety of chemical, physical, and biological processes that occur in the solution phase. The reference interaction site model (RISM) and its three-dimensional extension (3D-RISM) serve as powerful computational tools for modeling solvation effects in chemical reactions, biological functions, and structure formations. We present the RISM integrated calculator (RISMiCal) program package, which is based on RISM and 3D-RISM theories with fast GPU code. RISMiCal has been developed as an integrated RISM/3D-RISM program that has interfaces with external programs such as Gaussian16, GAMESS, and Tinker. Fast 3D-RISM programs for single- and multi-GPU codes written in CUDA would enhance the availability of these hybrid methods because they require the performance of many computationally expensive 3D-RISM calculations. We expect that our package can be widely applied for chemical and biological processes in solvent. The RISMiCal package is available at https://rismical-dev.github.io.

Free Energy Density of a Fluid and Its Role in Solvation and Binding

The concept that a fluid has a position-dependent free energy density appears in the literature but has not been fully developed or accepted. We set this concept on an unambiguous theoretical footing via the following strategy. First, we set forth four desiderata that should be satisfied by any definition of the position-dependent free energy density, f(R), in a system comprising only a fluid and a rigid solute: its volume integral, plus the fixed internal energy of the solute, should be the system free energy; it deviates from its bulk value, fbulk, near a solute but should asymptotically approach fbulk with increasing distance from the solute; it should go to zero where the solvent density goes to zero; and it should be well-defined in the most general case of a fluid made up of flexible molecules with an arbitrary interaction potential. Second, we use statistical thermodynamics to formulate a definition of the free energy density that satisfies these desiderata. Third, we show how any free energy density satisfying the desiderata may be used to analyze molecular processes in solution. In particular, because the spatial integral of f(R) equals the free energy of the system, it can be used to compute free energy changes that result from the rearrangement of solutes as well as the forces exerted on the solutes by the solvent. This enables the use of a thermodynamic analysis of water in protein binding sites to inform ligand design. Finally, we discuss related literature and address published concerns regarding the thermodynamic plausibility of a position-dependent free energy density. The theory presented here has applications in theoretical and computational chemistry and may be further generalizable beyond fluids, such as to solids and macromolecules.

Molecular Insights into the Impact of Mutations on the Binding Affinity of Targeted Covalent Inhibitors of BTK

Targeted covalent inhibitors (TCIs) have witnessed a significant resurgence in recent years, particularly in the kinase drug discovery field for treating diverse clinical indications. The inhibition of Bruton's tyrosine kinase (BTK) for treating B-cell cancers is a classic example where TCIs such as ibrutinib have had breakthroughs in targeted therapy. However, selectivity remains challenging, and the emergence of resistance mutations is a critical concern for clinical efficacy. Computational methods that can accurately predict the impact of mutations on inhibitor binding affinity could prove helpful in informing targeted approaches─providing insights into drug resistance mechanisms. In addition, such systems could help guide the systematic evaluation and impact of mutations in disease models for optimal experimental design. Here, we have employed in silico physics-based methods to understand the effects of mutations on the binding affinity and conformational dynamics of select TCIs of BTK. The TCIs studied include ibrutinib, acalabrutinib, and zanubrutinib─all of which are FDA-approved drugs for treating multiple forms of leukemia and lymphoma. Our results offer useful molecular insights into the structural determinants, thermodynamics, and conformational energies that impact ligand binding for this biological target of clinical relevance.

The Three-Dimensional Reference Interaction Site Model Approach as a Promising Tool for Studying Hydrated Viruses and Their Complexes with Ligands

Viruses are the most numerous biological form living in any ecosystem. Viral diseases affect not only people but also representatives of fauna and flora. The latest pandemic has shown how important it is for the scientific community to respond quickly to the challenge, including critically assessing the viral threat and developing appropriate measures to counter this threat. Scientists around the world are making enormous efforts to solve these problems. In silico methods, which allow quite rapid obtention of, in many cases, accurate information in this field, are effective tools for the description of various aspects of virus activity, including virus-host cell interactions, and, thus, can provide a molecular insight into the mechanism of virus functioning. The three-dimensional reference interaction site model (3D-RISM) seems to be one of the most effective and inexpensive methods to compute hydrated viruses, since the method allows us to provide efficient calculations of hydrated viruses, remaining all molecular details of the liquid environment and virus structure. The pandemic challenge has resulted in a fast increase in the number of 3D-RISM calculations devoted to hydrated viruses. To provide readers with a summary of this literature, we present a systematic overview of the 3D-RISM calculations, covering the period since 2010. We discuss various biophysical aspects of the 3D-RISM results and demonstrate capabilities, limitations, achievements, and prospects of the method using examples of viruses such as influenza, hepatitis, and SARS-CoV-2 viruses.

Identification of potential agonist-like molecules for α2-adrenergic receptor by multi-layer virtual screening to combat sinusitis

Sinusitis is one of the most common respiratory inflammatory conditions and a significant health issue that affects millions of people worldwide with a global prevalence of 10-15%. The side effects of available drug regimens of sinus infection demand the urgent development of new drug candidates to combat sinusitis. With the aim of identifying new drug-like candidates to control sinus, we have conducted multifold comprehensive screening of drug-like molecules targeting α2-adrenergic receptor (α2-AR), which serve as the primary drug target in sinusitis. By structure-based virtual screening of in-house compound's database, ten molecules (CP1-CP10) with agonistic effects for α2-AR were selected, and their binding mechanism with critical residues of α2-AR and their physicochemical properties were studied. Moreover, the process of receptor activation by these compounds and the conformational changes in α2-AR caused by these molecules, were further explored by molecular dynamic simulation. The MM-PBSA estimated free energies of compounds are higher than that of reference agonist (ΔGTOTAL = -39.0 kcal/mol). Among all, CP2-CP3, CP7-CP8 and CP6 have the highest binding free energies of -78.9 kcal/mol, -77.3 kcal/mol, -75.60 kcal/mol, -64.8 kcal/mol, and -61.6 kcal/mol, respectively. While CP4 (-55.0 kcal/mol), CP5 (-49.2 kcal/mol), CP9 (-54.8 ± 0.07 kcal/mol), CP10 (-56.7 ± 0.10 kcal/mol) and CP1 (-46.0 ± 0.08 kcal/mol) also exhibited significant binding free energies. These energetically favorable binding energies indicate strong binding affinity of our compounds for α2-AR as compared to known partial agonist. Therefore, these molecules can serve as excellent drug-like candidates for sinusitis.

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.

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.

Computational Strategies for Entropy Modeling in Chemical Processes

Computational simulations of entropy are important in understanding the thermodynamic forces that drive chemical reactions on a molecular scale. In recent years, various algorithms have been developed and applied in conjunction with molecular modeling techniques to evaluate the change of entropy in solvation, hydrophobic interactions, and chemical reactions. The aim of this review is to highlight four specific computational entropy calculation methods: normal mode analysis, free volume theory, two-phase thermodynamics, and configurational entropy modeling. The technical aspects, applications, and limitations of each method will be discussed in detail.

Protein 3D Hydration: A Case of Bovine Pancreatic Trypsin Inhibitor

Characterization of the hydrated state of a protein is crucial for understanding its structural stability and function. In the present study, we have investigated the 3D hydration structure of the protein BPTI (bovine pancreatic trypsin inhibitor) by molecular dynamics (MD) and the integral equation method in the three-dimensional reference interaction site model (3D-RISM) approach. Both methods have found a well-defined hydration layer around the protein and revealed the localization of BPTI buried water molecules corresponding to the X-ray crystallography data. Moreover, under 3D-RISM calculations, the obtained positions of waters bound firmly to the BPTI sites are in reasonable agreement with the experimental results mentioned above for the BPTI crystal form. The analysis of the 3D hydration structure (thickness of hydration shell and hydration numbers) was performed for the entire protein and its polar and non-polar parts using various cut-off distances taken from the literature as well as by a straightforward procedure proposed here for determining the thickness of the hydration layer. Using the thickness of the hydration shell from this procedure allows for calculating the total hydration number of biomolecules properly under both methods. Following this approach, we have obtained the thickness of the BPTI hydration layer of 3.6 Å with 369 water molecules in the case of MD simulation and 3.9 Å with 333 water molecules in the case of the 3D-RISM approach. The above procedure was also applied for a more detailed description of the BPTI hydration structure near the polar charged and uncharged radicals as well as non-polar radicals. The results presented for the BPTI as an example bring new knowledge to the understanding of protein hydration.

Structure of the MRAS-SHOC2-PP1C phosphatase complex

RAS-MAPK signalling is fundamental for cell proliferation and is altered in most human cancers1-3. However, our mechanistic understanding of how RAS signals through RAF is still incomplete. Although studies revealed snapshots for autoinhibited and active RAF-MEK1-14-3-3 complexes4, the intermediate steps that lead to RAF activation remain unclear. The MRAS-SHOC2-PP1C holophosphatase dephosphorylates RAF at serine 259, resulting in the partial displacement of 14-3-3 and RAF-RAS association3,5,6. MRAS, SHOC2 and PP1C are mutated in rasopathies-developmental syndromes caused by aberrant MAPK pathway activation6-14-and SHOC2 itself has emerged as potential target in receptor tyrosine kinase (RTK)-RAS-driven tumours15-18. Despite its importance, structural understanding of the SHOC2 holophosphatase is lacking. Here we determine, using X-ray crystallography, the structure of the MRAS-SHOC2-PP1C complex. SHOC2 bridges PP1C and MRAS through its concave surface and enables reciprocal interactions between all three subunits. Biophysical characterization indicates a cooperative assembly driven by the MRAS GTP-bound active state, an observation that is extendible to other RAS isoforms. Our findings support the concept of a RAS-driven and multi-molecular model for RAF activation in which individual RAS-GTP molecules recruit RAF-14-3-3 and SHOC2-PP1C to produce downstream pathway activation. Importantly, we find that rasopathy and cancer mutations reside at protein-protein interfaces within the holophosphatase, resulting in enhanced affinities and function. Collectively, our findings shed light on a fundamental mechanism of RAS biology and on mechanisms of clinically observed enhanced RAS-MAPK signalling, therefore providing the structural basis for therapeutic interventions.

Relative Binding Free-Energy Calculations at Lipid-Exposed Sites: Deciphering Hot Spots

Relative binding free-energy (RBFE) calculations are experiencing resurgence in the computer-aided drug design of novel small molecules due to performance gains allowed by cutting-edge molecular mechanic force fields and computer hardware. Application of RBFE to soluble proteins is becoming a routine, while recent studies outline necessary steps to successfully apply RBFE at the orthosteric site of membrane-embedded G-protein-coupled receptors (GPCRs). In this work, we apply RBFE to a congeneric series of antagonists that bind to a lipid-exposed, extra-helical site of the P2Y1 receptor. We find promising performance of RBFE, such that it may be applied in a predictive manner on drug discovery programs targeting lipid-exposed sites. Further, by the application of the microkinetic model, binding at a lipid-exposed site can be split into (1) membrane partitioning of the drug molecule followed by (2) binding at the extra-helical site. We find that RBFE can be applied to calculate the free energy of each step, allowing the uncoupling of observed binding free energy from the influence of membrane affinity. This protocol may be used to identify binding hot spots at extra-helical sites and guide drug discovery programs toward optimizing intrinsic activity at the target.

Accelerating Molecular Dynamics Enrichments of High-Affinity Ligands for Proteins

Molecular docking algorithms are used to seek the most active compounds from a pool of ligands. In principle, molecular dynamics (MD) simulations with accurate physical potentials and sampling could yield better enrichments, but they are computationally expensive. Here, we describe a method called MELD-Bracket that utilizes biased replica exchange ladders in MD in order to compete different ligands against each other within a fast bracket style "binding tournament". MELD-Bracket finds best-binders rapidly when ligands are well separated in their binding affinities.

Benchmarking Free Energy Calculations in Liquid Aliphatic Ketone Solvents Using the 3D-RISM-KH Molecular Solvation Theory

The three-dimensional reference interaction site model of the molecular solvation theory with the Kovalenko–Hirata closure is used to calculate the free energy of solvation of organic solutes in liquid aliphatic ketones. The ketone solvent sites were modeled using a modified united-atom force field. The successful application of these solvation models in calculating ketone–water partition coefficients of a large number of solutes supports the validation and benchmarking reported here.

Accurate prediction of hydration free energies and solvation structures using molecular density functional theory with a simple bridge functional

This paper assesses the ability of molecular density functional theory to predict efficiently and accurately the hydration free energies of molecular solutes and the surrounding microscopic water structure. A wide range of solutes were investigated, including hydrophobes, water as a solute, and the FreeSolv database containing 642 drug-like molecules having a variety of shapes and sizes. The usual second-order approximation of the theory is corrected by a third-order, angular-independent bridge functional. The overall functional is parameter-free in the sense that the only inputs are bulk water properties, independent of the solutes considered. These inputs are the direct correlation function, compressibility, liquid-gas surface tension, and excess chemical potential of the solvent. Compared to molecular simulations with the same force field and the same fixed solute geometries, the present theory is shown to describe accurately the solvation free energy and structure of both hydrophobic and hydrophilic solutes. Overall, the method yields a precision of order 0.5 k_BT for the hydration free energies of the FreeSolv database, with a computer speedup of 3 orders of magnitude. The theory remains to be improved for a better description of the H-bonding structure and the hydration free energy of charged solutes.

Field-Template, QSAR, Ensemble Molecular Docking, and 3D-RISM Solvation Studies Expose Potential of FDA-Approved Marine Drugs as SARS-CoVID-2 Main Protease Inhibitors

Currently, SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) has infected people among all countries and is a pandemic as declared by the World Health Organization (WHO). SARS-CoVID-2 main protease is one of the therapeutic drug targets that has been shown to reduce virus replication, and its high-resolution 3D structures in complex with inhibitors have been solved. Previously, we had demonstrated the potential of natural compounds such as serine protease inhibitors eventually leading us to hypothesize that FDA-approved marine drugs have the potential to inhibit the biological activity of SARS-CoV-2 main protease. Initially, field-template and structure–activity atlas models were constructed to understand and explain the molecular features responsible for SARS-CoVID-2 main protease inhibitors, which revealed that Eribulin Mesylate, Plitidepsin, and Trabectedin possess similar characteristics related to SARS-CoVID-2 main protease inhibitors. Later, protein–ligand interactions are studied using ensemble molecular-docking simulations that revealed that marine drugs bind at the active site of the main protease. The three-dimensional reference interaction site model (3D-RISM) studies show that marine drugs displace water molecules at the active site, and interactions observed are favorable. These computational studies eventually paved an interest in further in vitro studies. Finally, these findings are new and indeed provide insights into the role of FDA-approved marine drugs, which are already in clinical use for cancer treatment as a potential alternative to prevent and treat infected people with SARS-CoV-2.

Tuning Local Hydration Enables a Deeper Understanding of Protein-Ligand Binding: The PP1-Src Kinase Case

Water plays a key role in biomolecular recognition and binding. Despite the development of several computational and experimental approaches, it is still challenging to comprehensively characterize water-mediated effects on the binding process. Here, we investigate how water affects the binding of Src kinase to one of its inhibitors, PP1. Src kinase is a target for treating several diseases, including cancer. We use biased molecular dynamics simulations, where the hydration of predetermined regions is tuned at will. This computational technique efficiently accelerates the SRC-PP1 binding simulation and allows us to identify several key and yet unexplored aspects of the solvent's role. This study provides a further perspective on the binding phenomenon, which may advance the current drug design approaches for the development of new kinase inhibitors.

Molecular Recognition and Self-Organization in Life Phenomena Studied by a Statistical Mechanics of Molecular Liquids, the RISM/3D-RISM Theory

There are two molecular processes that are essential for living bodies to maintain their life: the molecular recognition, and the self-organization or self-assembly. Binding of a substrate by an enzyme is an example of the molecular recognition, while the protein folding is a good example of the self-organization process. The two processes are further governed by the other two physicochemical processes: solvation and the structural fluctuation. In the present article, the studies concerning the two molecular processes carried out by Hirata and his coworkers, based on the statistical mechanics of molecular liquids or the RISM/3D-RISM theory, are reviewed.

Instantaneous generation of protein hydration properties from static structures

Complex molecular simulation methods are typically required to calculate the thermodynamic properties of biochemical systems. One example thereof is the thermodynamic profiling of (de)solvation of proteins, which is an essential driving force for protein-ligand and protein-protein binding. The thermodynamic state of water molecules depends on its enthalpic and entropic components; the latter is governed by dynamic properties of the molecule. Here, we developed, to the best of our knowledge, two novel machine learning methods based on deep neural networks that are able to generate the converged thermodynamic state of dynamic water molecules in the heterogeneous protein environment based solely on the information of the static protein structure. The applicability of our machine learning methods to predict the hydration information is demonstrated in two different studies, the qualitative analysis and quantitative prediction of structure-activity relationships, and the prediction of protein-ligand binding modes.

Protein Matrix Control of Reaction Center Excitation in Photosystem II

Photosystem II (PSII) is a multisubunit pigment-protein complex that uses light-induced charge separation to power oxygenic photosynthesis. Its reaction center chromophores, where the charge transfer cascade is initiated, are arranged symmetrically along the D1 and D2 core polypeptides and comprise four chlorophyll (PD1, PD2, ChlD1, ChlD2) and two pheophytin molecules (PheoD1 and PheoD2). Evolution favored productive electron transfer only via the D1 branch, with the precise nature of primary excitation and the factors that control asymmetric charge transfer remaining under investigation. Here we present a detailed atomistic description for both. We combine large-scale simulations of membrane-embedded PSII with high-level quantum-mechanics/molecular-mechanics (QM/MM) calculations of individual and coupled reaction center chromophores to describe reaction center excited states. We employ both range-separated time-dependent density functional theory and the recently developed domain based local pair natural orbital (DLPNO) implementation of the similarity transformed equation of motion coupled cluster theory with single and double excitations (STEOM-CCSD), the first coupled cluster QM/MM calculations of the reaction center. We find that the protein matrix is exclusively responsible for both transverse (chlorophylls versus pheophytins) and lateral (D1 versus D2 branch) excitation asymmetry, making ChlD1 the chromophore with the lowest site energy. Multipigment calculations show that the protein matrix renders the ChlD1 → PheoD1 charge-transfer the lowest energy excitation globally within the reaction center, lower than any pigment-centered local excitation. Remarkably, no low-energy charge transfer states are located within the "special pair" PD1-PD2, which is therefore excluded as the site of initial charge separation in PSII. Finally, molecular dynamics simulations suggest that modulation of the electrostatic environment due to protein conformational flexibility enables direct excitation of low-lying charge transfer states by far-red light.

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.

Water Thermodynamics of Peptide Toxin Binding Sites on Ion Channels

Peptide toxins isolated from venomous creatures, long prized as research tools due to their innate potency for ion channels, are emerging as drugs as well. However, it remains challenging to understand why peptide toxins bind with high potency to ion channels, to identify residues that are key for activity, and to improve their affinities via mutagenesis. We use WaterMap, a molecular dynamics simulation-based method, to gain computational insight into these three questions by calculating the locations and thermodynamic properties of water molecules in the peptide toxin binding sites of five ion channels. These include an acid-sensing ion channel, voltage-gated potassium channel, sodium channel in activated and deactivated states, transient-receptor potential channel, and a nicotinic receptor whose structures were recently determined by crystallography and cryo-electron microscopy (cryo-EM). All channels had water sites in the peptide toxin binding site, and an average of 75% of these sites were stable (low-energy), and 25% were unstable (medium or high energy). For the sodium channel, more unstable water sites were present in the deactivated state structure than the activated. Additionally, for each channel, unstable water sites coincided with the positions of peptide toxin residues that previous mutagenesis experiments had shown were important for activity. Finally, for the sodium channel in the deactivated state, unstable water sites were present in the peptide toxin binding pocket but did not overlap with the peptide toxin, suggesting that future experimental efforts could focus on targeting these sites to optimize potency.

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.

Elucidating the multiple roles of hydration for accurate protein-ligand binding prediction via deep learning

Accurate and efficient prediction of protein-ligand interactions has been a long-lasting dream of practitioners in drug discovery. The insufficient treatment of hydration is widely recognized to be a major limitation for accurate protein-ligand scoring. Using an integration of molecular dynamics simulations on thousands of protein structures with novel big-data analytics based on convolutional neural networks and deep Taylor decomposition, we consistently identify here three different patterns of hydration to be essential for protein-ligand interactions. In addition to desolvation and water-mediated interactions, the formation of enthalpically favorable networks of first-shell water molecules around solvent-exposed ligand moieties is identified to be essential for protein-ligand binding. Despite being currently neglected in drug discovery, this hydration phenomenon could lead to new avenues in optimizing the free energy of ligand binding. Application of deep neural networks incorporating hydration to docking provides 89% accuracy in binding pose ranking, an essential step for rational structure-based drug design.

Utilizing Grand Canonical Monte Carlo Methods in Drug Discovery

The concepts behind targeting waters for potency and selectivity gains have been well documented and explored, although maximizing such potential gains can prove to be challenging. This problem is exacerbated in cases where there are multiple interacting waters, wherein perturbation of one water can affect the free energy landscape of the remaining waters. Knowing the right modification a priori is challenging, and computational approaches are ideally suited to help answer the key question of which substitution is best to try. Here, we use Grand Canonical Monte Carlo and the recent Grand Canonical Alchemical Perturbation methods to both understand and predict the effect of ligand-mediated water displacement when more than one water molecule is involved, as well as to understand how exploiting water networks can help govern selectivity.

Role of Displacing Confined Solvent in the Conformational Equilibrium of β-Cyclodextrin

This study investigates the role of hydration and its relationship to the conformational equilibrium of the host molecule β-cyclodextrin. Molecular dynamics simulations indicate that the unbound β-cyclodextrin exhibits two state behavior in explicit solvent due to the opening and closing of its cavity. In implicit solvent, these transitions are not observed, and there is one dominant conformation of β-cyclodextrin with an open cavity. Based on these observations, we investigate the hypothesis that the expulsion of thermodynamically unfavorable water molecules into the bulk plays an important role in controlling the accessibility of the closed macrostate at room temperature. We compare the results of the molecular mechanics analytical generalized Born plus nonpolar solvation approach to those obtained through grid inhomogeneous solvation theory analysis with explicit solvation to elucidate the thermodynamic forces at play. The work illustrates the use of continuum solvent models to tease out solvation effects related to the inhomogeneity and the molecular nature of water and demonstrates the key role of the thermodynamics of enclosed hydration in driving the conformational equilibrium of molecules in solution.

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.

A molecular reconstruction approach to site-based 3D-RISM and comparison to GIST hydration thermodynamic maps in an enzyme active site

Computed, high-resolution, spatial distributions of solvation energy and entropy can provide detailed information about the role of water in molecular recognition. While grid inhomogeneous solvation theory (GIST) provides rigorous, detailed thermodynamic information from explicit solvent molecular dynamics simulations, recent developments in the 3D reference interaction site model (3D-RISM) theory allow many of the same quantities to be calculated in a fraction of the time. However, 3D-RISM produces atomic-site, rather than molecular, density distributions, which are difficult to extract physical meaning from. To overcome this difficulty, we introduce a method to reconstruct molecular density distributions from atomic-site density distributions. Furthermore, we assess the quality of the resulting solvation thermodynamics density distributions by analyzing the binding site of coagulation Factor Xa with both GIST and 3D-RISM. We find good qualitative agreement between the methods for oxygen and hydrogen densities as well as direct solute-solvent energetic interactions. However, 3D-RISM predicts lower energetic and entropic penalties for moving water from the bulk to the binding site.

Characterizing hydration sites in protein-ligand complexes towards the design of novel ligands

Water is an essential part of protein binding sites and mediates interactions to ligands. Its displacement by ligand parts affects the free binding energy of resulting protein-ligand complexes. Therefore the characterization of solvation properties is important for design. Of particular interest is the propensity of localized water to be favorably displaced by a ligand. This review discusses two popular computational approaches addressing these questions, namely WaterMap based on statistical mechanics analysis of MD simulations and 3D RISM based on integral equation theory of liquids. The theoretical background and recent applications in structure-based design will be presented.

Discovery of Small-Molecule Antagonists of the H3K9me3 Binding to UHRF1 Tandem Tudor Domain

Ubiquitin-like with PHD and RING finger domains 1 (UHRF1) is a multidomain protein that plays a critical role in maintaining DNA methylation patterns through concurrent recognition of hemimethylated DNA and histone marks by various domains, and recruitment of DNA methyltransferase 1 (DNMT1). UHRF1 is overexpressed in various cancers, including breast cancer. The tandem tudor domain (TTD) of UHRF1 specifically and tightly binds to histone H3 di- or trimethylated at lysine 9 (H3K9me2 or H3K9me3, respectively), and this binding is essential for UHRF1 function. We developed an H3K9me3 peptide displacement assay, which was used to screen a library of 44,000 compounds for small molecules that disrupt the UHRF1-H3K9me3 interaction. This screen resulted in the identification of NV01, which bound to UHRF1-TTD with a K_d value of 5 μM. The structure of UHRF1-TTD in complex with NV01 confirmed binding to the H3K9me3-binding pocket. Limited structure-based optimization of NV01 led to the discovery of NV03 (K_d of 2.4 μM). These well-characterized small-molecule antagonists of the UHRF1-H3K9me2/3 interaction could be valuable starting chemical matter for developing more potent and cell-active probes toward further characterizing UHRF1 function, with possible applications as anticancer therapeutics.

A Systematic Analysis of the Binding Affinity between the Pim-1 Kinase and Its Inhibitors Based on the MM/3D-RISM/KH Method

A systematic study of the binding affinities of 16 lead compounds targeting the Pim-1 kinase based on the 3D-RISM/KH theory and MD simulations is reported. The results show a correlation coefficient R = 0.69 between the theoretical and experimental values of the binding free energy. This demonstrates that the method is applicable to the problem of compound screening and lead optimization, for which relative values of the free energy among the compounds have significance. We elucidate the contribution of the ligand fragments to the binding free energy. Our results indicate that the interactions between the residues and the triazolo[4,3-b]pyridazine scaffold as well as the phenyl ring of the ligand molecule make significant contributions to stabilization of the complex. Using the 3D-RISM/KH theory, we further analyze the probability distribution of a ligand fragment around the protein-ligand complex in which the substituent around the phenyl ring is removed from the ligand. The results demonstrate that the 3D-RISM/KH theory is capable of predicting the position of substitution on a ligand that has a higher affinity to a target protein.

Testing inhomogeneous solvation theory in structure-based ligand discovery

Binding-site water is often displaced upon ligand recognition, but is commonly neglected in structure-based ligand discovery. Inhomogeneous solvation theory (IST) has become popular for treating this effect, but it has not been tested in controlled experiments at atomic resolution. To do so, we turned to a grid-based version of this method, GIST, readily implemented in molecular docking. Whereas the term only improves docking modestly in retrospective ligand enrichment, it could be added without disrupting performance. We thus turned to prospective docking of large libraries to investigate GIST's impact on ligand discovery, geometry, and water structure in a model cavity site well-suited to exploring these terms. Although top-ranked docked molecules with and without the GIST term often overlapped, many ligands were meaningfully prioritized or deprioritized; some of these were selected for testing. Experimentally, 13/14 molecules prioritized by GIST did bind, whereas none of the molecules that it deprioritized were observed to bind. Nine crystal complexes were determined. In six, the ligand geometry corresponded to that predicted by GIST, for one of these the pose without the GIST term was wrong, and three crystallographic poses differed from both predictions. Notably, in one structure, an ordered water molecule with a high GIST displacement penalty was observed to stay in place. Inclusion of this water-displacement term can substantially improve the hit rates and ligand geometries from docking screens, although the magnitude of its effects can be small and its impact in drug binding sites merits further controlled studies.

Binding Pose Flip Explained via Enthalpic and Entropic Contributions

The anomalous binding modes of five highly similar fragments of TIE2 inhibitors, showing three distinct binding poses, are investigated. We report a quantitative rationalization for the changes in binding pose based on molecular dynamics simulations. We investigated five fragments in complex with the transforming growth factor β receptor type 1 kinase domain. Analyses of these simulations using Grid Inhomogeneous Solvation Theory (GIST), pK_A calculations, and a tool to investigate enthalpic differences upon binding unraveled the various thermodynamic contributions to the different binding modes. While one binding mode flip can be rationalized by steric repulsion, the second binding pose flip revealed a different protonation state for one of the ligands, leading to different enthalpic and entropic contributions to the binding free energy. One binding pose is stabilized by the displacement of entropically unfavored water molecules (binding pose determined by solvation entropy), ligands in the other binding pose are stabilized by strong enthalpic interactions, overcompensating the unfavorable water entropy in this pose (binding pose determined by enthalpic interactions). This analysis elucidates unprecedented details determining the flipping of the binding modes, which can elegantly explain the experimental findings for this system.

A combined treatment of hydration and dynamical effects for the modeling of host-guest binding thermodynamics: the SAMPL5 blinded challenge

As part of the SAMPL5 blinded experiment, we computed the absolute binding free energies of 22 host-guest complexes employing a novel approach based on the BEDAM single-decoupling alchemical free energy protocol with parallel replica exchange conformational sampling and the AGBNP2 implicit solvation model specifically customized to treat the effect of water displacement as modeled by the Hydration Site Analysis method with explicit solvation. Initial predictions were affected by the lack of treatment of ionic charge screening, which is very significant for these highly charged hosts, and resulted in poor relative ranking of negatively versus positively charged guests. Binding free energies obtained with Debye-Hückel treatment of salt effects were in good agreement with experimental measurements. Water displacement effects contributed favorably and very significantly to the observed binding affinities; without it, the modeling predictions would have grossly underestimated binding. The work validates the implicit/explicit solvation approach employed here and it shows that comprehensive physical models can be effective at predicting binding affinities of molecular complexes requiring accurate treatment of conformational dynamics and hydration.

How Mg2+ ion and water network affect the stability and structure of non-Watson-Crick base pairs in E. coli loop E of 5S rRNA: a molecular dynamics and reference interaction site model (RISM) study

The non-Watson-Crick (non-WC) base pairs of Escherichia coli loop E of 5S rRNA are stabilized by Mg²⁺ ions through water-mediated interaction. It is important to know the synergic role of Mg²⁺ and the water network surrounding Mg²⁺ in stabilizing the non-WC base pairs of RNA. For this purpose, free energy change of the system is calculated using molecular dynamics (MD) simulation as Mg²⁺ is pulled from RNA, which causes disturbance of the water network. It was found that Mg²⁺ remains hexahydrated unless it is close to or far from RNA. In the pentahydrated form, Mg²⁺ interacts directly with RNA. Water network has been identified by two complimentary methods; MD followed by a density-based clustering algorithm and three-dimensional-reference interaction site model. These two methods gave similar results. Identification of water network around Mg²⁺ and non-WC base pairs gives a clue to the strong effect of water network on the stability of this RNA. Based on sequence analysis of all Eubacteria 5s rRNA, we propose that hexahydrated Mg²⁺ is an integral part of this RNA and geometry of base pairs surrounding it adjust to accommodate the [Formula: see text]. Overall the findings from this work can help in understanding the basis of the complex structure and stability of RNA with non-WC base pairs.