Classification of water molecules in protein binding sites

How to replace the residual solvation shell of polar active site residues to achieve nanomolar inhibition of tRNA-guanine transglycosylase

In a computational and structural study, we investigated a series of 4-substituted lin-benzoguanines that are potent inhibitors of tRNA-guanine transglycosylase (TGT), a putative target for the treatment of shigellosis. At first glance, it appears self-evident that the placement of a positively charged ligand functional group between the carboxylate groups of two adjacent aspartate residues in the glycosylase catalytic center leads to enhanced ligand binding. The concomitant displacement of water molecules that partially solvate the aspartates prior to ligand binding appears to result as a consequence of this. However, the case study presented herein shows that this premise is much too superficial. Placement of a likely positively charged amino group at such a pivotal position, interfering with the residual water solvation shell, is at best cost-neutral compared with the unsubstituted parent ligand not conflicting with the residual water shell. A ligand that orients a hydroxy group in this position shows even decreased binding. Based on the cost-neutral placement of the amino functionality, hydrophobic side chains can now be further attached to fill, with increasing potency, a small hydrophobic pocket remote to the aspartates. Any attempts to cross the pivotal position between both aspartates with nonpolar scaffolds reveals only decreased binding, even though the waters of the residual solvation shell are successfully repelled. This surprising observation fostered a detailed analysis of the role of water molecules involved in the residual solvation of polar active site residues. Their geometry and putative replacement in the binding pocket of TGT has been studied by a comparative database analysis, computational active site mapping, and a series of crystal structure analyses. Furthermore, conformational preferences of attached hydrophobic moieties explain their contribution to a gradual increase in binding affinity.

Energetics of displacing water molecules from protein binding sites: consequences for ligand optimization

A strategy in drug design is to consider enhancing the affinity of lead molecules with structural modifications that displace water molecules from a protein binding site. Because success of the approach is uncertain, clarification of the associated energetics was sought in cases where similar structural modifications yield qualitatively different outcomes. Specifically, free-energy perturbation calculations were carried out in the context of Monte Carlo statistical mechanics simulations to investigate ligand series that feature displacement of ordered water molecules in the binding sites of scytalone dehydratase, p38-alphaMAP kinase, and EGFR kinase. The change in affinity for a ligand modification is found to correlate with the ease of displacement of the ordered water molecule. However, as in the EGFR example, the binding affinity may diminish if the free-energy increase due to the removal of the bound water molecule is not more than compensated by the additional interactions of the water-displacing moiety. For accurate computation of the effects of ligand modifications, a complete thermodynamic analysis is shown to be needed. It requires identification of the location of water molecules in the protein-ligand interface and evaluation of the free-energy changes associated with their removal and with the introduction of the ligand modification. Direct modification of the ligand in free-energy calculations is likely to trap the ordered molecule and provide misleading guidance for lead optimization.

Prediction of the water content in protein binding sites

An efficient molecular simulation methodology has been developed to determine the positioning of water molecules in the binding site of a protein or protein-ligand complex. Occupancies and absolute binding free energies of water molecules are computed using a statistical thermodynamics approach. The methodology, referred to as Just Add Water Molecules (JAWS), features "theta-water" molecules that can appear and disappear on a binding-site grid. Key approximations render the technique far more efficient than conventional free energy simulations. Testing of JAWS on five diverse examples (neuraminidase, scytalone dehydratase, major urinary protein 1, beta-lactoglobulin, and COX-2) demonstrates its accuracy in locating hydration sites in comparison to results from high-resolution crystal structures. Possible applications include aid in refinement of protein crystal structures, drug lead optimization, setup of docking calculations, and simulations of protein-ligand complexes.

Light chain somatic mutations change thermodynamics of binding and water coordination in the HyHEL-10 family of antibodies

Thermodynamic and structural studies addressed the increased affinity due to L-chain somatic mutations in the HyHEL-10 family of affinity matured IgG antibodies, using ITC, SPR with van't Hoff analysis, and X-ray crystallography. When compared to the parental antibody H26L26, the H26L10 and H26L8 chimeras binding to lysozyme showed an increase in favorable DeltaG(o) of -1.2+/-0.1 kcal mol(-1) and -1.3+/-0.1 kcal mol(-1), respectively. Increase in affinity of the H26L10 chimera was due to a net increase in favorable enthalpy change with little difference in change in entropy compared to H26L26. The H26L8 chimera exhibited the greatest increase in favorable enthalpy but also showed an increase in unfavorable entropy change, with the result being that the affinities of both chimeras were essentially equivalent. Site-directed L-chain mutants identified the shared somatic mutation S30G as the dominant contributor to increasing affinity to lysozyme. This mutation was not influenced by H-chain somatic mutations. Residue 30L is at the periphery of the binding interface and S30G effects an increase in hydrophobicity and decrease in H-bonding ability and size, but does not make any new energetically important antigen contacts. A new 1.2-A structure of the H10L10-HEL complex showed changes in the pattern of both inter- and intra-molecular water bridging with no other significant structural alterations near the binding interface compared to the H26L26-HEL complex. These results highlight the necessity for investigating both the structure and the thermodynamics associated with introduced mutations, in order to better assess and understand their impact on binding. Furthermore, it provides an important example of how backbone flexibility and water-bridging may favorably influence the thermodynamics of an antibody-antigen interaction.

A proline-based neuraminidase inhibitor: DFT studies on the zwitterion conformation, stability and formation

The designs of potent neuraminidase (NA) inhibitors are an efficient way to deal with the recent "2009 H1N1" influenza epidemic. In this work, density functional calculations were employed to study the conformation, stability and formation of the zwitterions of 5-[(1R,2S)-1-(acetylamino)-2-methoxy-2-methylpentyl]-4-[(1Z)-1-propenyl]-(4S,5R)-d-proline (BL), a proline-based NA inhibitor. Compared to proline, the zwitterion stability of BL is enhanced by 1.76 kcal mol(-1) due to the introduction of functional groups. However, the zwitterion of BL will not represent a local minimum on the potential energy surface until the number of water molecules increases up to two (n = 2). With the addition of two and three water molecules, the energy differences between the zwitterions and corresponding canonical isomers were calculated at 3.13 and -1.54 kcal mol(-1), respectively. The zwitterions of BL are mainly stabilized by the H-bonds with the water molecules, especially in the case of three water molecules where the carboxyl-O atoms are largely coordination-saturated by three H-bonds of medium strengths, causing the zwitterion stability even superior to the canonical isomer. With the presence of two and three water molecules, the energy barriers for the conversion processes from the canonical isomers to the zwitterions are equal to 4.96 and 3.13 kcal mol(-1), respectively. It indicated that the zwitterion formation is facile to take place with addition of two molecules and further facilitated by more water molecules. Besides, the zwitterion formation of BL is finished in a single step, different from other NA inhibitors. Owing to the above advantages, BL is a good NA inhibitor candidate and more attention should be paid to explorations of BL-based drugs.

Entropic cost of protein-ligand binding and its dependence on the entropy in solution

Two theoretical formulations are proposed and compared for the loss of translational and rotational entropy upon protein-ligand binding in water. The two theories share the same approach to evaluate the translational and rotational entropy of the ligand when bound. The potential of the bound ligand is modeled by six harmonic oscillators that are parametrized from the force and torque magnitudes measured in a molecular dynamics simulation, yielding vibrational and librational entropies. In the aqueous phase, the theories differ because there is no unique way to assign the total entropy to molecules in solution. In one approach, the ligand is allowed unrestricted access to the full solution volume at the standard concentration and is assigned the same translational and rotational entropy as if it were an ideal gas. We term this a "molecule-frame" (MF) theory because it considers configurational space in the reference frame of the molecule of interest. The entropy of the solvent is penalized because it is excluded from the molecule's volume. In the second theory, all molecules including the solvent are confined by their neighbors in mean-field configurational volumes. This we term a "system-frame" (SF) theory because the configurational space available to all molecules is considered in the reference frame of the whole system. Molecules have vibrational and librational entropy in the same way as they do when bound. In addition, the discrete size of the solvent molecules quantizes the configurational space into an effective number of minima according to the solute molecule's standard concentration and the mean volume of a solvent molecule. This leads to the cratic entropy expressed in terms of the solute molecule's mole fraction. The equivalent number of minima in rotational space depends on both the solute molecule's volume and the solvent molecule's volume. This leads to an equation for the orientational entropy based on the proposed concept of "angle fraction". The MF and SF theories are applied to calculate the translational and rotational entropy losses involved in the formation of six different protein-ligand complexes, in two of which the ligand is water. The MF entropy losses range from -80 to -142 J K(-1) mol(-1) for ligands at the 1 M standard-state concentration and from -52 to -63 J K(-1) mol(-1) for water at the 55.6 M standard-state concentration. They depend logarithmically on both the number and strength of interactions between the ligand and protein through the forces and torques. This is observed to lead to moderate dependencies on the ligands' moments of inertia and masses. The SF entropy losses are smaller and range from -50 to -75 J K(-1) mol(-1) for ligands at the 1 M standard-state concentration and from 0 to -12 J K(-1) mol(-1) for water. They depend logarithmically on the ligand solvent's molecular volume and weakly on the relative strengths of the ligand's interactions with the protein and water. The cratic entropy loss in water at the standard concentration is constant and is also demonstrated to be implicit in MF theories. Entropy losses from the two approaches are also compared with those from other computational approaches and with experiment. The use of the force and torque magnitudes leads to smaller bound volumes than are obtained from ligand-displacement approaches. The general agreement of the SF entropy losses with those from experiment suggests that the SF theory is more consistent with the assumptions made in experimental measurements than the MF solvation theories, which would require a compensating entropy gain in the solvent in order to agree.

X-ray crystal structures of Phanerochaete chrysosporium Laminarinase 16A in complex with products from lichenin and laminarin hydrolysis

The 1,3(4)-beta-D-glucanases of glycoside hydrolase family 16 provide useful examples of versatile yet specific protein-carbohydrate interactions. In the present study, we report the X-ray structures of the 1,3(4)-beta-D-glucanase Phanerochaete chrysosporium Laminarinase 16A in complex with beta-glucan products from laminarin (1.6 A) and lichenin (1.1 A) hydrolysis. The G6G3G3G glucan, in complex with the enzyme, showed a beta-1,6 branch in the acceptor site. The G4G3G ligand-protein complex showed that there was no room for a beta-1,6 branch in the -1 or -2 subsites; furthermore, the distorted residue in the -1 subsite and the glucose in the -2 subsite required a beta-1,3 bond between them. These are the first X-ray crystal structures of any 1,3(4)-beta-D-glucanase in complex with glucan products. They provide details of both substrate and product binding in support of earlier enzymatic evidence.

Binding of small-molecule ligands to proteins: "what you see" is not always "what you get"

We review insights from computational studies of affinities of ligands binding to proteins. The power of structural biology is in translating knowledge of protein structures into insights about their forces, binding, and mechanisms. However, the complementary power of computer modeling is in showing "the rest of the story" (i.e., how motions and ensembles and alternative conformers and the entropies and forces that cannot be seen in single molecular structures also contribute to binding affinities). Upon binding to a protein, a ligand can bind in multiple orientations; the protein or ligand can be deformed by the binding event; waters, ions, or cofactors can have unexpected involvement; and conformational or solvation entropies can sometimes play large and otherwise unpredictable roles. Computer modeling is helping to elucidate these factors.

A molecular dynamics approach to study the importance of solvent in protein interactions

Water constitutes the cellular environment for biomolecules to interact. Solvent is important for protein folding and stability, and it is also known to actively participate in many catalytic processes in the cell. However, solvent is often ignored in molecular recognition and not taken into account in protein-protein interaction studies and rational design. Previously we developed SCOWLP, a database and its web application (http://www.scowlp.org), to perform studies on the contribution of solvent to protein interface definition in all protein complexes of the PDB. We introduced the concept of wet spots, interfacial residues interacting only through one water molecule, which were shown to considerably enrich protein interface descriptions. Analysis of interfacial solvent in a nonredundant dataset of protein complexes suggested the importance of including interfacial water molecules in protein interaction studies. In this work we use a molecular dynamics approach to gain deeper insights into solvent contribution to protein interfaces. We characterize the dynamic and energetic properties of water-mediated protein interactions by comparing different interfacial interaction types (direct, dual and wet spot) at residue and solvent level. For this purpose, we perform an analysis of 17 representative complexes from two protein families of different interface nature. Energetically wet spots are quantitatively comparable to other residues in interfaces, and their mobility is shown to be lower than protein surface residues. The residence time of water molecules in wet spots sites is higher than of those on the surface of the protein. In terms of free energy, though wet-spots-forming water molecules are very heterogeneous, their contribution to the free energy of complex formation is considerable. We find that water molecules can play an important role in interaction conservation in protein interfaces by allowing sequence variability in the corresponding binding partner, and we discuss the important implications of our observations related to the use of the correlated mutations concept in protein interactions studies. The results obtained in this work help to deepen our understanding of the physico-chemical nature underlying protein-protein interactions and strengthen the idea of using the wet spots concept to qualitatively improve the accuracy of folding, docking and rational design algorithms.

Computations of standard binding free energies with molecular dynamics simulations

An increasing number of studies have reported computations of the standard (absolute) binding free energy of small ligands to proteins using molecular dynamics (MD) simulations and explicit solvent molecules that are in good agreement with experiments. This encouraging progress suggests that physics-based approaches hold the promise of making important contributions to the process of drug discovery and optimization in the near future. Two types of approaches are principally used to compute binding free energies with MD simulations. The most widely known is the alchemical double decoupling method, in which the interaction of the ligand with its surroundings are progressively switched off. It is also possible to use a potential of mean force (PMF) method, in which the ligand is physically separated from the protein receptor. For both of these computational approaches, restraining potentials may be activated and released during the simulation for sampling efficiently the changes in translational, rotational, and conformational freedom of the ligand and protein upon binding. Because such restraining potentials add bias to the simulations, it is important that their effects be rigorously removed to yield a binding free energy that is properly unbiased with respect to the standard state. A review of recent results is presented, and differences in computational methods are discussed. Examples of computations with T4-lysozyme mutants, FKBP12, SH2 domain, and cytochrome P450 are discussed and compared. Remaining difficulties and challenges are highlighted.

Computational evaluation of protein-small molecule binding

Determining protein-small molecule binding affinity is a key component of present-day rational drug discovery. To circumvent the time, labor, and materials costs associated with experimental protein-small molecule binding assays, a variety of structure-based computational methods have been developed for determining protein-small molecule binding affinities. These methods can be placed in one of two classes: accurate but slow (Class 1), and fast but approximate (Class 2). Class 1 methods, which explicitly take into account protein flexibility and include an atomic-level description of solvation, are capable of quantitatively reproducing experimental protein-small molecule absolute binding free energies. However, Class 1 computational requirements make screening thousands to millions of small molecules against a protein, as required for rational drug design, infeasible for the foreseeable future. Class 2 methods, on the contrary, are sufficiently fast to perform such inhibitor screening, yet they suffer from limited descriptions of protein flexibility and solvation, which in turn limit their ability to select and rank-order small molecules by computed binding affinities. This review presents an overview of Class 1 and Class 2 methods, and avenues of research in Class 2 methods aimed at bringing them closer to Class 1 accuracy.

Exploiting ordered waters in molecular docking

A current weakness in docking is the treatment of water-mediated protein-ligand interactions. We explore switching ordered water molecules "on" and "off" during docking screens of a large library. The method assumes additivity and scales linearly with the number of waters sampled despite the exponential growth in configurations. It is tested for ligand enrichment against 24 targets, exploring up to 256 water configurations. Water inclusion increased enrichment substantially for 12 targets, while most others were largely unaffected.

Computation of binding free energy with molecular dynamics and grand canonical Monte Carlo simulations

The binding of a ligand to a receptor is often associated with the displacement of a number of bound water molecules. When the binding site is exposed to the bulk region, this process may be sampled adequately by standard unbiased molecular dynamics trajectories. However, when the binding site is deeply buried and the exchange of water molecules with the bulk region may be difficult to sample, the convergence and accuracy in free energy perturbation (FEP) calculations can be severely compromised. These problems are further compounded when a reduced system including only the region surrounding the binding site is simulated. To address these issues, we couple molecular dynamics (MD) with grand canonical Monte Carlo (GCMC) simulations to allow the number of water to fluctuate during an alchemical FEP calculation. The atoms in a spherical inner region around the binding pocket are treated explicitly while the influence of the outer region is approximated using the generalized solvent boundary potential (GSBP). At each step during thermodynamic integration, the number of water in the inner region is equilibrated with GCMC and energy data generated with MD is collected. Free energy calculations on camphor binding to a deeply buried pocket in cytochrome P450cam, which causes about seven water molecules to be expelled, are used to test the method. It concluded that solvation free energy calculations with the GCMC/MD method can greatly improve the accuracy of the computed binding free energy compared to simulations with fixed number of water.

Towards the development of universal, fast and highly accurate docking/scoring methods: a long way to go

Accelerating the drug discovery process requires predictive computational protocols capable of reducing or simplifying the synthetic and/or combinatorial challenge. Docking-based virtual screening methods have been developed and successfully applied to a number of pharmaceutical targets. In this review, we first present the current status of docking and scoring methods, with exhaustive lists of these. We next discuss reported comparative studies, outlining criteria for their interpretation. In the final section, we describe some of the remaining developments that would potentially lead to a universally applicable docking/scoring method.

Ligand-protein docking with water molecules

The presence of water molecules plays an important role in the accuracy of ligand-protein docking predictions. Comprehensive docking simulations have been performed on a large set of ligand-protein complexes whose crystal structures contain water molecules in their binding sites. Only those water molecules found in the immediate vicinity of both the ligand and the protein were considered. We have investigated whether prior optimization of the orientation of water molecules in either the presence or absence of the bound ligand has any effect on the accuracy of docking predictions. We have observed a statistically significant overall increase in accuracy when water molecules are included during docking simulations and have found this to be independent of the method of optimization of the orientation of water molecules. These results confirm the importance of including water molecules whenever possible in a ligand-protein docking simulation. Our findings also reveal that prior optimization of the orientation of water molecules, in the absence of any bound ligand, does not have a detrimental effect on the improved accuracy of ligand-protein docking. This is important, given the use of docking simulations to predict the binding modes of new ligands or drug molecules.

Synthesis of urea picket porphyrins and their use in the elucidation of the role buried solvent plays in the selectivity and stoichiometry of anion binding receptors

The synthesis of alpha,alpha-5,10-diurea and alpha,alpha,alpha-5,10,15-triurea picket porphyrins are detailed in this report. In previous reports, these porphyrins, along with alpha,alpha,alpha,alpha-5,10,15,20-tetraurea picket porphyrin, were used to demonstrate the important role one buried solvent molecule plays in the selectivity and stoichiometry of binding inorganic anions. Building on prior work, this report discusses the results of acetate anion binding studies between tetra- and diurea picket porphyrins (the latter does not contain a buried solvent molecule in the anion-receptor complex), compares differences in thermodynamic data obtained from van't Hoff plots of a porphyrin anion receptor able to utilize buried solvent in its binding motif with one that does not, and compares the crystal structure of a tetraurea porphyrin 1-chloride anion complex that contains buried solvent with new X-ray crystal structures of tetraurea porphyrin 1-dichloride or bisdihydrogenphosphate anion complexes that contain no buried solvent. Data from our previous work, and the work described herein, demonstrates that one buried solvent molecule provides stability to the receptor-anion complex that is similar in energy to a moderately strong hydrogen bond.