Hydration in drug design. 2. Influence of local site surface shape on water binding

Summary If water molecules are strongly bound at a protein-ligand interface, they are unlikely to be displaced during ligand binding. Such water molecules can change the shape of the ligand binding site and thus affect strategies for drug design. To understand the nature of water binding, and factors influencing it, water molecules at the ligand binding sites of 26 high-resolution protein-ligand complexes have been examined here. Water molecules bound in deep grooves and cavities between the protein and the ligand are located in the indentations on the protein-site surface, but not in the indentations on the ligand surface. The majority of the water molecules bound in deep indentations on the protein-site surface make multiple polar contacts with the protein surface. This may indicate a strong binding of water molecules in deep indentations on protein-site surfaces. The local shape of the site surface may influence the binding of water molecules that mediate protein-ligand interactions.     

Hydration in drug design. 3. Conserved water molecules at the ligand-binding sites of homologous proteins

Summary Water molecules are known to play an important rôle in mediating protein-ligand interactions. If water molecules are conserved at the ligand-binding sites of homologous proteins, such a finding may suggest the structural importance of water molecules in ligand binding. Structurally conserved water molecules change the conventional definition of binding sites by changing the shape and complementarity of these sites. Such conserved water molecules can be important for site-directed ligand/drug design. Therefore, five different sets of homologous protein/protein-ligand complexes have been examined to identify the conserved water molecules at the ligand-binding sites. Our analysis reveals that there are as many as 16 conserved water molecules at the FAD binding site of glutathione reductase between the crystal structures obtained from human and E. coli. In the remaining four sets of high-resolution crystal structures, 2–4 water molecules have been found to be conserved at the ligand-binding sites. The majority of these conserved water molecules are either bound in deep grooves at the protein-ligand interface or completely buried in cavities between the protein and the ligand. All these water molecules, conserved between the protein/protein-ligand complexes from different species, have identical or similar apolar and polar interactions in a given set. The site residues interacting with the conserved water molecules at the ligand-binding sites have been found to be highly conserved among proteins from different species; they are more conserved compared to the other site residues interacting with the ligand. These water molecules, in general, make multiple polar contacts with protein-site residues.     

Role of bridging water molecules in GSK3β-inhibitor complexes: insights from QM/MM, MD, and molecular docking studies

The role of water molecules is increasingly gaining interest in drug design, and several studies have highlighted their paramount contributions to the specificity and the affinity of ligand binding. In this study, we employ the two-layer ONIOM-based quantum mechanics/molecular mechanics (QM/MM) calculations, molecular dynamics (MD) simulations, and molecular docking studies to investigate the effect of bridging water molecules at the GSK3β-inhibitors interfaces. The results obtained from the ONIOM geometry optimization and AIM analysis corroborated the presence of bridging water molecules that form hydrogen bonds with protein side chain of Thr138 and/or backbone of Gln185, and mediate interactions with inhibitors in the 10 selected GSK3β-inhibitor complexes. Subsequently, MD simulations carried out on a representative system of 1R0E demonstrated that the bridging water molecule is stable at the GSK3β-inhibitor interface and appears to contribute to the stability of the protein-inhibitor interactions. Furthermore, molecular docking studies of GSK3β-inhibitor complexes indicated that the inhibitors can increase binding affinities and the better docked conformation of inhibitors can be obtained by inclusion of the bridging water molecules, especially for the flexible inhibitors, in docking experiments into individual protein conformations. Our results elucidate the importance of bridging water molecules at the GSK3β-inhibitor interfaces and suggest that they might prove useful in rational drug design.     

Analysis of ligand-bound water molecules in high-resolution crystal structures of protein-ligand complexes

We have performed a comprehensive analysis of water molecules at the protein-ligand interfaces observed in 392 high-resolution crystal structures. There are a total of 1829 ligand-bound water molecules in these 392 complexes; 18% are surface water molecules, and 72% are interfacial water molecules. The number of ligand-bound water molecules in each complex structure ranges from 0 to 21 and has an average of 4.6. Of these interfacial water molecules, 76% are considered to be bridging water molecules, characterized by having polar interactions with both ligand and protein atoms. Among a number of factors that may influence the number of ligand-bound water molecules, the polar van der Waals (vdw) surface area of ligands has the highest Pearson linear correlation coefficient of 0.63. Our regression analysis predicted that one more ligand-bound water molecule is expected for every additional 24 A2 in the polar vdw surface area of the ligand. In contrast to the observation that the resolution is the primary factor influencing the number of water molecules in crystallographic models of proteins, we found that there is only a weak relationship between the number of ligand-bound water molecules and the resolution of the crystal structures. An analysis of the isotropic B factors of buried ligand-bound water molecules suggested that, when water molecules have fewer than two polar interactions with the protein-ligand complex, they are more mobile than protein atoms in the crystal structures; when they have more than three polar interactions, they are significantly less mobile than protein atoms.     

A new method for estimating the importance of hydrogen-bonding groups in the binding site of a protein

We introduce a new method to estimate the importance of hydrogen-bonding sitepoints in the binding site of a protein as part of a structure-based design strategy. Our method identifies hydrogen-bonding sitepoints within a binding pocket and ranks them according to both the accessibility of their hydrogen bonding regions to incoming ligands and their hydrogen-bonding strength. The combination of these components produces a prioritised list of sitepoints that are more likely to be involved in hydrogen bonding with an incoming ligand. A dataset of known protein-ligand interactions was used to compare the prioritisation of sitepoints identified by our method with those observed to be engaged in hydrogen bonding in their crystal structures. Our method was able to remove those sitepoints unable to bind the ligand due to a low accessibility or an unfavourable orientation and to award significantly higher hydrogen-bonding ranking values to those sitepoints observed to form hydrogen bonds. Our method can thus be used to identify hydrogen-bonding sitepoints that should be targeted preferentially in a drug design strategy.     

A method for fast energy estimation and visualization of protein-ligand interaction

Summary A new computational and graphical method for facilitating ligand-protein docking studies is developed on a three-dimensional computer graphics display. Various physical and chemical properties inside the ligand binding pocket of a receptor protein, whose structure is elucidated by X-ray crystal analysis, are calculated on three-dimensional grid points and are stored in advance. By utilizing those tabulated data, it is possible to estimate the non-bonded and electrostatic interaction energy and the number of possible hydrogen bonds between protein and ligand molecules in real time during an interactive docking operation. The method also provides a comprehensive visualization of the local environment inside the binding pocket.With this method, it becomes easier to find a roughly stable geometry of ligand molecules, and one can therefore make a rapid survey of the binding capability of many drug candidates. The method will be useful for drug design as well as for the examination of protein-ligand interactions.     

GREEN: A program package for docking studies in rational drug design

Summary A program package, GREEN, has been developed that enables docking studies between ligand molecules and a protein molecule. Based on the structure of the protein molecule, the physical and chemical environment of the ligand-binding site is expressed as three-dimensional grid-point data. The grid-point data are used for the real-time evaluation of the protein-ligand interaction energy, as well as for the graphical representation of the binding-site environment. The interactive docking operation is facilitated by various built-in functions, such as energy minimization, energy contribution analysis and logging of the manipulation trajectory. Interactive modeling functions are incorporated for designing new ligand molecules while considering the binding-site environment and the protein-ligand interaction. As an example of the application of GREEN, a docking study is presented on the complex between trypsin and a synthetic trypsin inhibitor. The program package will be useful for rational drug design, based on the 3D structure of the target protein.     

A water-swap reaction coordinate for the calculation of absolute protein-ligand binding free energies

The accurate prediction of absolute protein-ligand binding free energies is one of the grand challenge problems of computational science. Binding free energy measures the strength of binding between a ligand and a protein, and an algorithm that would allow its accurate prediction would be a powerful tool for rational drug design. Here we present the development of a new method that allows for the absolute binding free energy of a protein-ligand complex to be calculated from first principles, using a single simulation. Our method involves the use of a novel reaction coordinate that swaps a ligand bound to a protein with an equivalent volume of bulk water. This water-swap reaction coordinate is built using an identity constraint, which identifies a cluster of water molecules from bulk water that occupies the same volume as the ligand in the protein active site. A dual topology algorithm is then used to swap the ligand from the active site with the identified water cluster from bulk water. The free energy is then calculated using replica exchange thermodynamic integration. This returns the free energy change of simultaneously transferring the ligand to bulk water, as an equivalent volume of bulk water is transferred back to the protein active site. This, directly, is the absolute binding free energy. It should be noted that while this reaction coordinate models the binding process directly, an accurate force field and sufficient sampling are still required to allow for the binding free energy to be predicted correctly. In this paper we present the details and development of this method, and demonstrate how the potential of mean force along the water-swap coordinate can be improved by calibrating the soft-core Coulomb and Lennard-Jones parameters used for the dual topology calculation. The optimal parameters were applied to calculations of protein-ligand binding free energies of a neuraminidase inhibitor (oseltamivir), with these results compared to experiment. These results demonstrate that the water-swap coordinate provides a viable and potentially powerful new route for the prediction of protein-ligand binding free energies.     

Molecular docking with ligand attached water molecules

A novel approach to incorporate water molecules in protein-ligand docking is proposed. In this method, the water molecules display the same flexibility during the docking simulation as the ligand. The method solvates the ligand with the maximum number of water molecules, and these are then retained or displaced depending on energy contributions during the docking simulation. Instead of being a static part of the receptor, each water molecule is a flexible on/off part of the ligand and is treated with the same flexibility as the ligand itself. To favor exclusion of the water molecules, a constant entropy penalty is added for each included water molecule. The method was evaluated using 12 structurally diverse protein-ligand complexes from the PDB, where several water molecules bridge the ligand and the protein. A considerable improvement in successful docking simulations was found when including flexible water molecules solvating hydrogen bonding groups of the ligand. The method has been implemented in the docking program Molegro Virtual Docker (MVD).     

Involvement of water in carbohydrate-protein binding.

The interactions of trimannosides 1 and 2 with Con A were studied to reveal the effects of displacement of well-ordered water molecules on the thermodynamic parameters of protein-ligand complexation. Trisaccharide 2 is a derivative of 1, in which the hydroxyl at C-2 of the central mannose unit is replaced by a hydroxyethyl moiety. Upon binding, this moiety displaces a conserved water molecule present in the Con A binding site. Structural studies by NMR spectroscopy and MD simulations showed that the two compounds have very similar solution conformational properties. MD simulations of the complexes of Con A with 1 and 2 demonstrated that the hydroxyethyl side chain of 2 can establish the same hydrogen bonds in a low energy conformation with the protein binding site as those mediated by the water molecule in the complex of 1 with Con A. Isothermal titration microcalorimetry (ITC) measurements showed that 2 has a more favorable entropy of binding compared to 1. This term, which was expected, arises from the return of the highly ordered water molecule to bulk solution. The favorable entropy term was, however, offset by a relatively large unfavorable enthalpy term. This observation was rationalized by comparing the extent of hydrogen bond and solvation changes during binding. It is proposed that an indirect interaction through a water molecule will provide a larger number of hydrogen bonds in the complex that have higher occupancies than in bulk solution, thereby stabilizing the complex.     

人体端粒中(3+1)混合结构G-四链体稳定性的分子动力学模拟

利用分子动力学模拟方法, 考察了人体端粒中(3+1)混合结构G-四链体的结构及稳定性问题. 讨论了配位K+离子、药物分子(端粒抑素)和溶剂水分子对G-四链体的Hoogsteen氢键结构、π-π堆积作用的影响. 研究表明, K+离子与鸟嘌呤碱基上O6原子的配位作用减弱了对角鸟嘌呤间O6-O6的静电排斥作用, 使得相邻的四个鸟嘌呤能够以Hoogsteen氢键结合的方式形成具有近平面结构的稳定G-四平面. 另一方面, G-四平面间、G-四平面与药物分子间的π-π堆积作用降低了G-四链体复合物的总能, 有利于其稳定存在. 此外, 溶剂水分子主要分布在G-四链体的TTA环、骨架和糖环的周围, 使其位移涨落增大; 然而, 在3 ns动力学模拟中, 由于水分子没有进入到G-四链体的空腔中, 溶剂水对G-四平面的结构影响不明显.     

Prediction of binding constants of protein ligands: A fast method for the prioritization of hits obtained from de novo design or 3D database search programs

A dataset of 82 protein–ligand complexes of known 3D structure and binding constant K_i was analysed to elucidate the important factors that determine the strength of protein–ligand interactions. The following parameters were investigated: the number and geometry of hydrogen bonds and ionic interactions between the protein and the ligand, the size of the lipophilic contact surface, the flexibility of the ligand, the electrostatic potential in the binding site, water molecules in the binding site, cavities along the protein–ligand interface and specific interactions between aromatic rings. Based on these parameters, a new empirical scoring function is presented that estimates the free energy of binding for a protein–ligand complex of known 3D structure. The function distinguishes between buried and solvent accessible hydrogen bonds. It tolerates deviations in the hydrogen bond geometry of up to 0.25 Å in the length and up to 30 °Cs in the hydrogen bond angle without penalizing the score. The new energy function reproduces the binding constants (ranging from 3.7 × 10^-2 M to 1 × 10^-14 M, corresponding to binding energies between -8 and -80 kJ/mol) of the dataset with a standard deviation of 7.3 kJ/mol corresponding to 1.3 orders of magnitude in binding affinity. The function can be evaluated very fast and is therefore also suitable for the application in a 3D database search or de novo ligand design program such as LUDI. The physical significance of the individual contributions is discussed.     

Ligand-specific scoring functions: improved ranking of docking solutions

Molecular docking is a powerful computational method that has been widely used in many biomolecular studies to predict geometry of a protein-ligand complex. However, while its conformational search algorithms are usually able to generate correct conformation of a ligand in the binding site, the scoring methods often fail to discriminate it among many false variants. We propose to treat this problem by applying more precise ligand-specific scoring filters to re-rank docking solutions. In this way specific features of interactions between protein and different types of compounds can be implicitly taken into account. New scoring functions were constructed including hydrogen bonds, hydrophobic and hydrophilic complementarity terms. These scoring functions also discriminate ligands by the size of the molecule, the total hydrophobicity, and the number of peptide bonds for peptide ligands. Weighting coefficients of the scoring functions were adjusted using a training set of 60 protein-ligand complexes. The proposed method was then tested on the results of docking obtained for an additional 70 complexes. In both cases the success rate was 5-8% better compared to the standard functions implemented in popular docking software.     

Conserved core substructures in the overlay of protein-ligand complexes

The method of conserved core substructure matching (CSM) for the overlay of protein-ligand complexes is described. The method relies upon distance geometry to align structurally similar substructures without regard to sequence similarity onto substructures from a reference protein empirically selected to include key determinants of binding site location and geometry. The error in ligand position is reduced in reoriented ensembles generated with CSM when compared to other overlay methods. Since CSM can only succeed when the selected core substructure is geometrically conserved, misalignments only rarely occur. The method may be applied to reliably overlay large numbers of protein-ligand complexes in a way that optimizes ligand position at a specific binding site or subsite or to align structures from large and diverse protein families where the conserved binding site is localized to only a small portion of either protein. Core substructures may be complex and must be chosen with care. We have created a database of empirically selected core substructures to demonstrate the utility of CSM alignment of ligand binding sites in important drug targets. A Web-based interface can be used to apply CSM to align large collections of protein-ligand complexes for use in drug design using these substructures or to evaluate the use of alternative core substructures that may then be shared with the larger user community. Examples show the benefit of CSM in the practice of structure-based drug design.     

海藻糖和氨基酸之间相互作用的分子动力学模拟

虽然海藻糖已经广泛用于蛋白质稳定性研究,但海藻糖稳定蛋白质的作用机理尚不清晰.本文利用全原子分子动力学模拟研究了20种常见氨基酸和海藻糖之间的分子机理.结果表明,所有氨基酸,尤其是极性和带电氨基酸,均优先与水分子结合.相反,仅有疏水性氨基酸与海藻糖发生相互作用,尤其是芳香族和疏水性氨基酸的侧链更易于和海藻糖接触.所有氨基酸的主链与水分子接触的趋势一致.虽然氨基酸和海藻糖与水之间均形成氢键,但氨基酸和海藻糖之间的氢键相互作用要弱于氨基酸和水之间的氢键相互作用.上述分子模拟的结果对于海藻糖稳定蛋白质作用机理的解析及高效蛋白质稳定剂的理性设计具有非常重要的理论指导意义.     

海藻糖和氨基酸之间相互作用的分子动力学模拟

虽然海藻糖已经广泛用于蛋白质稳定性研究,但海藻糖稳定蛋白质的作用机理尚不清晰. 本文利用全原子分子动力学模拟研究了20种常见氨基酸和海藻糖之间的分子机理. 结果表明,所有氨基酸,尤其是极性和带电氨基酸,均优先与水分子结合. 相反,仅有疏水性氨基酸与海藻糖发生相互作用,尤其是芳香族和疏水性氨基酸的侧链更易于和海藻糖接触. 所有氨基酸的主链与水分子接触的趋势一致. 虽然氨基酸和海藻糖与水之间均形成氢键,但氨基酸和海藻糖之间的氢键相互作用要弱于氨基酸和水之间的氢键相互作用. 上述分子模拟的结果对于海藻糖稳定蛋白质作用机理的解析及高效蛋白质稳定剂的理性设计具有非常重要的理论指导意义.     

AcquaAlta: a directional approach to the solvation of ligand-protein complexes

Water molecules mediating polar interactions in ligand-protein complexes can substantially contribute to binding affinity and specificity. To account for such water molecules in computer-aided drug design, we performed an extensive search in the Cambridge Structural Database (CSD) to identify the geometrical criteria defining interactions of water molecules with ligand and protein. In addition, with ab initio calculations the propensity of ligand hydration was evaluated. Based on this information, we developed an algorithm (AcquaAlta) to reproduce water molecules bridging polar interactions between ligand and protein moieties. This approach was validated with 20 crystal structures and yielded a match of 76% between experimental and calculated water positions. When water molecules establishing only weak interactions with the protein were neglected, the match could be improved to 88%. Supported by a pharmacophore-based alignment tool, the solvation algorithm was then applied to the docking of oligopeptides to the periplasmic oligopeptide binding protein A (OppA). Calculated waters based on the crystal poses matched an average of 66% of the experimental waters. With water molecules calculated based on the docked ligands, the average match with the experimental waters dropped to 53%.     

Ligand-specific scoring functions: improved ranking of docking solutions1

Molecular docking is a powerful computational method that has been widely used in many biomolecular studies to predict geometry of a protein-ligand complex. However, while its conformational search algorithms are usually able to generate correct conformation of a ligand in the binding site, the scoring methods often fail to discriminate it among many false variants. We propose to treat this problem by applying more precise ligand-specific scoring filters to re-rank docking solutions. In this way specific features of interactions between protein and different types of compounds can be implicitly taken into account. New scoring functions were constructed including hydrogen bonds, hydrophobic and hydrophilic complementarity terms. These scoring functions also discriminate ligands by the size of the molecule, the total hydrophobicity, and the number of peptide bonds for peptide ligands. Weighting coefficients of the scoring functions were adjusted using a training set of 60 protein–ligand complexes. The proposed method was then tested on the results of docking obtained for an additional 70 complexes. In both cases the success rate was 5–8% better compared to the standard functions implemented in popular docking software.     

Dynamic properties of water around a protein-DNA complex from molecular dynamics simulations

Formation of protein-DNA complex is an important step in regulation of genes in living organisms. One important issue in this problem is the role played by water in mediating the protein-DNA interactions. In this work, we have carried out atomistic molecular dynamics simulations to explore the heterogeneous dynamics of water molecules present in different regions around a complex formed between the DNA binding domain of human TRF1 protein and a telomeric DNA. It is demonstrated that such heterogeneous water motions around the complex are correlated with the relaxation time scales of hydrogen bonds formed by those water molecules with the protein and DNA. The calculations reveal the existence of a fraction of extraordinarily restricted water molecules forming a highly rigid thin layer in between the binding motifs of the protein and DNA. It is further proved that higher rigidity of water layers around the complex originates from more frequent reformations of broken water-water hydrogen bonds. Importantly, it is found that the formation of the complex affects the transverse and longitudinal degrees of freedom of surrounding water molecules in a nonuniform manner.     

2D depiction of protein-ligand complexes

A method is presented for the automated preparation of schematic diagrams for protein-ligand complexes, in which the ligand is displayed in conventional 2D form, and the interactions to and between the residues in its vicinity are summarized in a concise and information-rich manner. The structural entities are arranged to maximize aesthetic ideals and to properly convey important distance relationships. The diagram is annotated with calculated hydrogen bonds, a substitution contour, solvent exposure, chelated metals, covalently bound linkages, pi-pi and pi-cation interactions, and, for series of complexes, conserved residues and interactions. Residues, cofactors, ions, and solvent components are drawn in cartoon form as adjuncts to the ligand. The method can be applied to aligned sets which contain multiple ligands, or multiple members of a protein family, in which case the ligand orientations and protein residue placement will show consistent trends throughout the series.