Standard free energy of releasing a localized water molecule from the binding pockets of proteins: double-decoupling method

Localized water molecules in the binding pockets of proteins play an important role in noncovalent association of proteins and small drug compounds. At times, the dominant contribution to the binding free energy comes from the release of localized water molecules in the binding pockets of biomolecules. Therefore, to quantify the energetic importance of these water molecules for drug design purposes, we have used the double-decoupling approach to calculate the standard free energy of tying up a water molecule in the binding pockets of two protein complexes. The double-decoupling approach is based on the underlying principle of statistical thermodynamics. We have calculated the standard free energies of tying up the water molecule in the binding pockets of these complexes to be favorable. These water molecules stabilize the protein-drug complexes by interacting with the ligands and binding pockets. Our results offer ideas that could be used in optimizing protein-drug interactions, by designing ligands that are capable of targeting localized water molecules in protein binding sites. The resulting free energy of ligand binding could benefit from the potential free energy gain accompanying the release of these water molecules. Furthermore, we have examined the theoretical background of the double-decoupling method and its connection to the molecular dynamics thermodynamic integration techniques.     

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.     

Nonionic side chains modulate the affinity and specificity of binding between functionalized polyamines and structured RNA

This Communication introduces side-chain-bearing polyamines as molecules for selective recognition of folded RNA structures. The complex folded structures associated with RNA create binding pockets for proteins, and also binding sites for small molecules. Developing organic molecules that can bind RNA with high affinity and specificity is a challenge that must be overcome for RNA to be considered a viable drug target. In this work, six polyamines with different side chains were synthesized to test for effects on binding affinity and specificity to TAR RNA and RRE RNA of HIV. Binding interactions between polyamines and RNAs were examined using two footprinting assays, based on terbium-induced cleavage and magnesium-catalyzed cleavage at higher pH. The binding constants and the binding specificity were highly dependent on the side chains of the polyamines, demonstrating that this class of molecules is a very promising starting point for development of highly selective RNA-binding ligands.     

Structural analysis of charge discrimination in the binding of inhibitors to human carbonic anhydrases I and II

Despite the similarity in the active site pockets of carbonic anhydrase (CA) isozymes I and II, the binding affinities of benzenesulfonamide inhibitors are invariably higher with CA II as compared to CA I. To explore the structural basis of this molecular recognition phenomenon, we have designed and synthesized simple benzenesulfonamide inhibitors substituted at the para position with positively charged, negatively charged, and neutral functional groups, and we have determined the affinities and X-ray crystal structures of their enzyme complexes. The para-substituents are designed to bind in the midsection of the 15 A deep active site cleft, where interactions with enzyme residues and solvent molecules are possible. We find that a para-substituted positively charged amino group is more poorly tolerated in the active site of CA I compared with CA II. In contrast, a para-substituted negatively charged carboxylate substituent is tolerated equally well in the active sites of both CA isozymes. Notably, enzyme-inhibitor affinity increases upon neutralization of inhibitor charged groups by amidation or esterification. These results inform the design of short molecular linkers connecting the benzenesulfonamide group and a para-substituted tail group in "two-prong" CA inhibitors: an optimal linker segment will be electronically neutral, yet capable of engaging in at least some hydrogen bond interactions with protein residues and/or solvent. Microcalorimetric data reveal that inhibitor binding to CA I is enthalpically less favorable and entropically more favorable than inhibitor binding to CA II. This contrasting behavior may arise in part from differences in active site desolvation and the conformational entropy of inhibitor binding to each isozyme active site.     

<Emphasis Type="Italic">In-silico</Emphasis> Screening using Flexible Ligand Binding Pockets: A Molecular Dynamics-based Approach

Summary In-silico screening of flexible ligands against flexible ligand binding pockets (LBP) is an emerging approach in structure-based drug discovery. Here, we describe a molecular dynamics (MD) based docking approach to investigate the influence on the high-throughput in-silico screening of small molecules against flexible ligand binding pockets. In our approach, an ensemble of 51 energetically favorable structures of the LBP of human estrogen receptor α (hERα) were collected from 3 ns MD simulations. In-silico screening of 3500 endocrine disrupting compounds against these flexible ligand binding pockets resulted in thousands of ER–ligand complexes of which 582 compounds were unique. Detailed analysis of MD generated structures showed that only 17 of the LBP residues significantly contribute to the overall binding pocket flexibility. Using the flexible LBP conformations generated, we have identified 32 compounds that bind better to the flexible ligand-binding pockets compared to the crystal structure. These compounds, though chemically divergent, are structurally similar to the natural hormone. Our MD-based approach in conjunction with grid–based distributed computing could be applied routinely for in-silico screening of large databases against any given target.     

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.     

MCSS-based predictions of RNA binding sites

The diversity of RNA tertiary structures provides the basis for specific recognition by proteins or small molecules. To investigate the structural basis and the energetics which control RNA-ligand interactions, favorable RNA binding sites are identified using the MCSS method, which has been employed previously only for protein receptors. Two different RNAs for which the structures have been determined by NMR spectroscopy were examined: two structures of the TAR RNA which contains an arginine binding site, and the structure of the 16S rRNA which contains an aminoglycoside binding site (paromomycin). In accord with the MCSS methodology, the functional groups representing the entire ligand or only part of it (one residue in the case of the aminoglycosides) are first replicated and distributed with random positions and orientations around the target and then energy minimized in the force field of the target RNA. The Coulombic term and the dielectric constant of the force field are adjusted to approximate the effects of solvent-screening and counterions. Optimal force field parameters are determined to reproduce the binding mode of arginine to the TAR RNA. The more favorable binding sites for each residue of the aminoglycoside ligands are then calculated and compared with the binding sites observed experimentally. The predictability of the method is evaluated and refinements are proposed to improve its accuracy. Received: 24 April 1998 / Accepted: 4 August 1998 / Published online: 7 December 1998     

Hydration in drug design. 1. Multiple hydrogen-bonding features of water molecules in mediating protein-ligand interactions

Summary Water is known to play an important rôle in the recognition and stabilization of the interaction between a ligand and its site. This has important implications for drug design. Analyses of 19 high-resolution crystal structures of protein-ligand complexes reveal the multiple hydrogen-bonding feature of water molecules mediating protein-ligand interactions. Most of the water molecules (nearly 80%) involved in bridging the protein and the ligand can make three or more hydrogen bonds when distance and bond angles are used as criteria to define hydrogen-bonding interactions. Isotropic B-factors have been used to take into account the mobility of water molecules. The water molecules at binding sites bridge the protein and ligand, and interact with other water molecules to form a complex network of interconnecting hydrogen bonds. Some water molecules at the site do not directly bridge between the protein and the ligand, but may contribute indirectly to the stability of the complex by holding bridging water molecules in the right position through a network of hydrogen bonds. These water networks are probably crucial for the stability of the protein-ligand complex and are important for any site-directed drug design strategies.     

Characterization of the galectin-1 carbohydrate recognition domain in terms of solvent occupancy

Human galectin-1, a galactosil-terminal sugar binding soluble protein, is a potent multifunctional effector that participates in specific protein-carbohydrate and protein-protein interactions. Recent studies revealed that it plays a key role as a modulator of cellular differentiation and immunological response. In this work, we have investigated the solvation properties of the carbohydrate recognition domain of Gal-1 by means of molecular dynamics simulations. Water sites (ws) were identified in terms of radial and angular distribution functions, and properties such as water residence times, interaction energies, and free-energy contributions were evaluated for those sites. Our results allowed us to correlate the thermodynamic properties of the ws and their binding pattern with the N-acetilgalactoside ligand. These results let us further infer that the water molecules located at the ws, which exhibit much more favorable binding, are the ones replaced by -OH groups of the sugar.     

Solvent models for protein–ligand binding: Comparison of implicit solvent poisson and surface generalized born models with explicit solvent simulations

Solvent effects play a crucial role in mediating the interactions between proteins and their ligands. Implicit solvent models offer some advantages for modeling these interactions, but they have not been parameterized on such complex problems, and therefore, it is not clear how reliable they are. We have studied the binding of an octapeptide ligand to the murine MHC class I protein using both explicit solvent and implicit solvent models. The solvation free energy calculations are more than 10³ faster using the Surface Generalized Born implicit solvent model compared to FEP simulations with explicit solvent. For some of the electrostatic calculations needed to estimate the binding free energy, there is near quantitative agreement between the explicit and implicit solvent model results; overall, the qualitative trends in the binding predicted by the explicit solvent FEP simulations are reproduced by the implicit solvent model. With an appropriate choice of reference system based on the binding of the discharged ligand, electrostatic interactions are found to enhance the binding affinity because the favorable Coulomb interaction energy between the ligand and protein more than compensates for the unfavorable free energy cost of partially desolvating the ligand upon binding. Some of the effects of protein flexibility and thermal motions on charging the peptide in the solvated complex are also considered. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 591–607, 2001     

What induces pocket openings on protein surface patches involved in protein–protein interactions?

We previously showed for the proteins BCL-X_L, IL-2, and MDM2 that transient pockets at their protein–protein binding interfaces can be identified by applying the PASS algorithm to molecular dynamics (MD) snapshots. We now investigated which aspects of the natural conformational dynamics of proteins induce the formation of such pockets. The pocket detection protocol was applied to three different conformational ensembles for the same proteins that were extracted from MD simulations of the inhibitor bound crystal conformation in water and the free crystal/NMR structure in water and in methanol. Additional MD simulations studied the impact of backbone mobility. The more efficient CONCOORD or normal mode analysis (NMA) techniques gave significantly smaller pockets than MD simulations, whereas tCONCOORD generated pockets comparable to those observed in MD simulations for two of the three systems. Our findings emphasize the influence of solvent polarity and backbone rearrangements on the formation of pockets on protein surfaces and should be helpful in future generation of transient pockets as putative ligand binding sites at protein–protein interfaces. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Thermodynamics of buried water clusters at a protein-ligand binding interface

The structure of the complex of cyclophilin A (CypA) with cyclosporin A (CsA, 1) shows a cluster of four water molecules buried at the binding interface, which is rearranged when CsA is replaced by (5-hydroxynorvaline)-2-cyclosporin (2). The thermodynamic contributions of each bound water molecule in the two complexes are explored with the inhomogeneous fluid solvation theory and molecular dynamics simulations. Water (WTR) 133 in complex 1 contributes little to the binding affinity, while WTR6 and 7 in complex 2 play an essential role in mediating protein-ligand binding with a hydrogen bond network. The calculations reveal that the rearrangement of the water molecules contributes favorably to the binding affinity, even though one of them is displaced going from ligand 1 to 2. Another favorable contribution comes from the larger protein-ligand interactions of ligand 2. However, these favorable contributions are not sufficient to overcome the unfavorable desolvation free energy change and the conformational entropy of the hydroxylpropyl group of ligand 2 in the complex, leading to a lower binding affinity of ligand 2. These physical insights may be useful in the development of improved scoring functions for binding affinity prediction.     

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%.     

Why not consider a spherical protein? Implications of backbone hydrogen bonding for protein structure and function

The intrinsic ability of protein structures to exhibit the geometric features required for molecular function in the absence of evolution is examined in the context of three systems: the reference set of real, single domain protein structures, a library of computationally generated, compact homopolypeptides, artificial structures with protein-like secondary structural elements, and quasi-spherical random proteins packed at the same density as proteins but lacking backbone secondary structure and hydrogen bonding. Without any evolutionary selection, the library of artificial structures has similar backbone hydrogen bonding, global shape, surface to volume ratio and statistically significant structural matches to real protein global structures. Moreover, these artificial structures have native like ligand binding cavities, and a tiny subset has interfacial geometries consistent with native-like protein-protein interactions and DNA binding. In contrast, the quasi-spherical random proteins, being devoid of secondary structure, have a lower surface to volume ratio and lack ligand binding pockets and intermolecular interaction interfaces. Surprisingly, these quasi-spherical random proteins exhibit protein like distributions of virtual bond angles and almost all have a statistically significant structural match to real protein structures. This implies that it is local chain stiffness, even without backbone hydrogen bonding, and compactness that give rise to the likely completeness of the library solved single domain protein structures. These studies also suggest that the packing of secondary structural elements generates the requisite geometry for intermolecular binding. Thus, backbone hydrogen bonding plays an important role not only in protein structure but also in protein function. Such ability to bind biological molecules is an inherent feature of protein structure; if combined with appropriate protein sequences, it could provide the non-zero background probability for low-level function that evolution requires for selection to occur.     

HIV-1整合酶与芳香二酮酸类抑制剂相互作用的分子模拟研究

用分子对接方法研究了一系列芳香二酮酸类抑制剂与HIV-1整合酶的识别及相互作用. 结果表明, 抑制剂结合到整合酶Asp64～Leu68, Thr115～Phe121, Gln148～Lys159和Mg2＋所构成的口袋区, 抑制机理与5CITEP相似. 采用分子动力学模拟和MM/PBSA方法计算了芳香二酮酸类抑制剂与整合酶之间的结合自由能, 计算结果与实验值相吻合, 平均绝对偏差为3.6 kJ/mol, 体系范德华相互作用和溶剂化效应的非极性项是利于形成复合物的主要因素. 相关性分析结果表明, 结合自由能值与疏水相互作用有较强的线性相关(R＝0.61), 基于此, 用多元线性回归方法给出了一个能较强预测芳香二酮酸类抑制剂与HIV-1整合酶的结合自由能预测模型, 为后续基于抑制剂结构的抗HIV-1药物分子设计提供指导.     

RNA and the Small Molecule World

A key role in essential cellular processes is played by RNA molecules, and these are attractive targets for drug design. The functional diversity of RNA can be attributed to the sophisticated three-dimensional structures it assumes. These intricate folds create potential binding pockets for ions, low molecular weight ligands, and proteins. Recent experiments have demonstrated that small molecules such as tobramycin ( 1 ) can regulate gene expression in living cells through specific interactions with a messenger RNA (mRNA).     

Aggregates of poly-functional amphiphilic molecules in water and oil phases

The solvation and aggregate formation of complex amphiphilic molecules such as tetra-acids in polar and nonpolar phases are studied via Molecular Dynamics simulations. The nonpolar core of tetra-acid molecules is found to be effectively impermeable for water molecules resulting in a low solubility in the polar solvent, while nonpolar solvent molecules sufficiently solvate the amphiphilic molecules considered, enabling an open conformation of their molecular structure. The rigidity of the core region of the tetra-acid molecules has been found to play a crucial role in their behavior in both polar and nonpolar phases. In the polar phase, simulations have shown that tetra-acids form micelle-like structures with a small aggregation number, confirming previous experimental work. The identification of a case of study in which micelle-like structures can form only with a small aggregation number enables the study via Molecular Dynamics of micelle-micelle interactions. Micelle stability and dispersion in the polar phase under different conditions can be therefore investigated. In the nonpolar phase, the preferential interactions between carboxyl groups, the affinity of the tetra-acids with the solvent molecules, and the structural characteristics of the central core moiety of the tetra-acids have been found to possibly induce a web like array, or network.     

WATsite: Hydration site prediction program with PyMOL interface

Water molecules that mediate protein–ligand interactions or are released from the binding site on ligand binding can contribute both enthalpically and entropically to the free energy of ligand binding. To elucidate the thermodynamic profile of individual water molecules and their potential contribution to ligand binding, a hydration site analysis program WATsite was developed together with an easy‐to‐use graphical user interface based on PyMOL. WATsite identifies hydration sites from a molecular dynamics simulation trajectory with explicit water molecules. The free energy profile of each hydration site is estimated by computing the enthalpy and entropy of the water molecule occupying a hydration site throughout the simulation. The results of the hydration site analysis can be displayed in PyMOL. A key feature of WATsite is that it is able to estimate the protein desolvation free energy for any user specified ligand. The WATsite program and its PyMOL plugin are available free of charge from http://people.pnhs.purdue.edu/~mlill/software . © 2014 Wiley Periodicals, Inc.     

Water molecules inside protein structure affect binding of monosaccharides with HIV‐1 antibody 2G12

Water molecules inside biomolecules constitute integral parts of their structure and participate in the functions of the proteins. Some of the X‐ray crystallographic data are insufficient for analyzing a series of ligand–protein complexes in the same condition. We theoretically investigated antibody binding abilities of saccharide ligands and the effects of the inner water molecules of ligand–antibody complexes. Classical molecular dynamics and quantum chemical simulations using a model with possible water molecules inside the protein were performed with saccharide ligands and Human Immunodeficiency Virus 1 neutralizing antibody 2G12 complexes to estimate how inner water molecules of the protein affect the dynamics of the complexes as well as the ligand–antibody interaction. Our results indicate the fact that d ‐fructose's strong affinity to the antibody was partly due to the good retentiveness of solvent water molecules of the ligand and its stability of the ligand's conformation and relative position in the active site. © 2016 Wiley Periodicals, Inc.