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.     

Pharmacophore anchor models of ATAT1 to discover potential inhibitors and lead optimization

Post-translation modification of microtubules is associated with many diseases like cancer. Alpha Tubulin Acetyltransferase 1 (ATAT1) is a major enzyme that acetylates ‘Lys-40’ in alpha-tubulin on the luminal side of microtubules and is a drug target that lacks inhibitors. Here, we developed pharmacophore anchor models of ATAT1 which were constructed statistically using thousands of docked compounds, for drug design and investigating binding mechanisms. Our models infer the compound moiety preferences with the physico-chemical properties for the ATAT1 binding site. The results from the pharmacophore anchor models show the three main sub-pockets, including S1 acetyl site, S2 adenine site, and S3 diphosphate site with anchors, where conserved moieties interact with respective sub-pocket residues in each site and help in guiding inhibitor discovery. We validated these key anchors by analyzing 162 homologous protein sequences (>99 species) and over 10 structures with various bound ligands and mutations. Our results were consistent with previous works also providing new interesting insights. Our models applied in virtual screening predicted several ATAT1 potential inhibitors. We believe that our model is useful for future inhibitor discovery and for guiding lead optimization.     

Pharmacophore fingerprint-based approach to binding site subpocket similarity and its application to bioisostere replacement

Bioisosteres have been defined as structurally different molecules or substructures that can form comparable intermolecular interactions, and therefore, fragments that bind to similar protein structures exhibit a degree of bioisosterism. We present KRIPO (Key Representation of Interaction in POckets): a new method for quantifying the similarities of binding site subpockets based on pharmacophore fingerprints. The binding site fingerprints have been optimized to improve their performance for both intra- and interprotein family comparisons. A range of attributes of the fingerprints was considered in the optimization, including the placement of pharmacophore features, whether or not the fingerprints are fuzzified, and the resolution and complexity of the pharmacophore fingerprints (2-, 3-, and 4-point fingerprints). Fuzzy 3-point pharmacophore fingerprints were found to represent the optimal balance between computational resource requirements and the identification of potential replacements. The complete PDB was converted into a database comprising almost 300?000 optimized fingerprints of local binding sites together with their associated ligand fragments. The value of the approach is demonstrated by application to two crystal structures from the Protein Data Bank: (1) a MAP kinase P38 structure in complex with a pyridinylimidazole inhibitor ( 1A9U ) and (2) a complex of thrombin with melagatran ( 1K22 ). Potentially valuable bioisosteric replacements for all subpockets of the two studied protein are identified.     

X-ray spectroscopy study of the Cu(II)GHK peptide complex

The Cu(II)–Gly–His–Lys (Glycyl–Histidyl–Lysine) complex is of interest as a model peptide to test the methodology for studying the structure of metal sites in proteins, in particular, the copper binding site in amyloid-β. X-ray absorption spectra of the Cu(II)GHK aqueous solution are measured. The stability of the complex under X-ray radiation is controlled by optical spectroscopy. The structural models with different copper site coordination constructed based on the crystallographic structure are considered. Two optimal models are selected from the analysis of the theoretical X-ray absorption spectra of the constructed structures. The structural parameters of the selected models are optimized. It is found that the spectrum of the five-coordinated model with water molecules in the equatorial and axial positions “down” (with Cu–O distances of 1.97 Å and 2.31 Å respectively) has the best agreement with the experiment.     

De novo ligand design with explicit water molecules: an application to bacterial neuraminidase

Most computer-aided drug design methods ignore the presence of crystallographically-determined water molecules in the binding site of a target protein. In this paper, our de novo ligand design methods are applied to the X-ray crystal structure of bacterial neuraminidase in the presence of some selected water molecules. We have found that, for this particular protein, the complete removal of all bound water molecules leads to difficulties in generating any potential ligands if the unsatisfied hydrogen-bonding sitepoints left by removing these water molecules are to be satisfied by a ligand. As more of the crystallographically determined water molecules are allowed in the binding site, it becomes much easier to generate ligands in larger numbers and with wider chemical diversity. This example shows that, in some cases, bound water molecules can be more accessible for hydrogen bonding to an incoming ligand than the actual protein binding sitepoints associated with them. From the point of view of de novo ligand design, water molecules can thus act as versatile amphiprotic hydrogen-bonding sitepoints and reduce the conformational constraints of a particular binding site.     

Including tightly-bound water molecules in de novo drug design. Exemplification through the in silico generation of poly(ADP-ribose)polymerase ligands

Different strategies for the in silico generation of ligand molecules in the binding site of poly(ADP-ribose)polymerase (PARP) were studied in order to observe the effect of the targeting and displacement of tightly bound water molecules. Several molecular scaffolds were identified as having better interactions in the binding site when targeting one or two tightly bound water molecules in the NAD binding site. Energy calculations were conducted in order to assess the ligand-protein and ligand-water-protein interactions of different functional groups of the generated ligands. These calculations were used to evaluate the energetic consequences of the presence of tightly bound water molecules and to identify those that contribute favorably to the binding of ligands.     

Using pharmacophore models to gain insight into structural binding and virtual screening: an application study with CDK2 and human DHFR

This study provides results from two case studies involving the application of the HypoGenRefine algorithm within Catalyst for the automated generation of excluded volume from ligand information alone. A limitation of pharmacophore feature hypothesis alone is that activity prediction is based purely on the presence and arrangement of pharmacophoric features; steric effects remained unaccounted. Recently reported studies have illustrated the usefulness of combining excluded volumes to the pharmacophore models. In general, these excluded volumes attempt to penalize molecules occupying steric regions that are not occupied by active molecules. The HypoGenRefine algorithm in Catalyst accounts for steric effects on activity, based on the targeted addition of excluded volume features to the pharmacophores. The automated inclusion of excluded volumes to pharmacophore models has been applied to two systems: CDK2 and human DHFR. These studies are used as examples to illustrate how ligands could bind in the protein active site with respect to allowed and disallowed binding regions. Additionally, automated refinement of the pharmacophore with these excluded volume features provides a more selective model to reduce false positives and a better enrichment rate in virtual screening.     

Steric and allosteric effects of fatty acids on the binding of warfarin to human serum albumin revealed by molecular dynamics and free energy calculations

Human serum albumin (HSA) binds with drugs and fatty acids (FAs). This study was initiated to elucidate the relationship between the warfarin binding affinity of HSA and the positions of bound FA molecules. Molecular dynamics simulations of 11 HSA-warfarin-myristate complexes were performed. HSA-warfarin binding free energy was then calculated for each of the complexes by the molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) method. The results indicated that the magnitude of the binding free energy was smaller in HSA-warfarin complexes that had 4 or more myristate molecules than in complexes with no myristate molecules. The unfavorable effect on the HSA-warfarin binding affinity was caused sterically by the binding of a myristate molecule to the FA binding site closest to the warfarin binding site. On the other hand, the magnitude of HSA-warfarin binding free energy was largest when 3 myristate molecules were bound to the high-affinity sites. The strongest HSA-warfarin binding was attributable to favorable entropic contribution related to larger atomic fluctuations of the amino acid residues at the warfarin binding site. In the binding of 2 myristate molecules to the sites with the highest and second-highest affinities, allosteric modulation that enhanced electrostatic interactions between warfarin and some of the amino acid residues around the warfarin binding site was observed. This study clarified the structural and energetic properties of steric/allosteric effects of FAs on the HSA-warfarin binding affinity and illustrated the approach to analyze protein-ligand interactions in situations such that multiple ligands bind to the other sites of the protein.     

Supramolecular photochirogenesis with biomolecules. Mechanistic studies on the enantiodifferentiation for the photocyclodimerization of 2-anthracenecarboxylate mediated by bovine serum albumin

Photophysics and photochemistry of 2-anthracenecarboxylate (AC) bound to bovine serum albumin (BSA) were investigated in detail for the first time by electronic absorption, circular dichroism (CD), steady-state and time-resolved fluorescence, fluorescence quenching, and product analysis studies. Through the spectroscopic investigations, it was revealed that the four independent binding pockets of BSA, which are known to accommodate 1, 3, 2, and 3 AC molecules in the order of decreasing affinity, are distinctly different in hydrophobicity, chiral environment, and accessibility. Interestingly, AC bound to site 1 gave highly structured fluorescence with dual lifetimes of 4.8 and 2.1 ns in an intensity ratio of 3:2, which may be assigned to the existence of two positional or orientational isomers within the very hydrophobic site 1. In contrast, the lifetime of AC in site 2 was much longer (13.3 ns), and ACs in sites 3 and 4 have broader fluorescence spectra with lifetimes that were practically indistinguishable from that in bulk water (15.8 ns). Although each of sites 2-4 simultaneously binds multiple AC molecules, no CD exciton coupling or static fluorescence quenching was detected, indicating that ACs bound to each site are not in close proximity to each other. Quenching studies with nitromethane further confirmed the significant difference in accessibility among the binding sites; thus, ACs bound to sites 1 and 2 are highly protected from the attack of the quencher, affording 32 and 10 times smaller rate constants than that for free AC in water. Product studies in the presence and absence of nitromethane more clearly revealed the photochirogenic performance of each binding site. Although the addition of nitromethane did not greatly alter the product distribution, the enantiomeric excesses (ee's) of chiral cycloadducts 2 and 3 were critically manipulated by selectively retarding the photoreaction occurring at the more accessible binding sites. Thus, the highest ee of 38% was obtained for 2 in the presence of 18 mM nitromethane, while the highest ee of 58% was attained for 3 in the absence of nitromethane, both at [AC]/[BSA]=3.6.     

Classification of water molecules in protein binding sites

Water molecules play a crucial role in mediating the interaction between a ligand and a macromolecular receptor. An understanding of the nature and role of each water molecule in the active site of a protein could greatly increase the efficiency of rational drug design approaches: if the propensity of a water molecule for displacement can be determined, then synthetic effort may be most profitably applied to the design of specific ligands with the displacement of this water molecule in mind. In this paper, a thermodynamic analysis of water molecules in the binding sites of six proteins, each complexed with a number of inhibitors, is presented. Two classes of water molecules were identified: those conserved and not displaced by any of the ligands, and those that are displaced by some ligands. The absolute binding free energies of 54 water molecules were calculated using the double decoupling method, with replica exchange thermodynamic integration in Monte Carlo simulations. It was found that conserved water molecules are on average more tightly bound than displaced water molecules. In addition, Bayesian statistics is used to calculate the probability that a particular water molecule may be displaced by an appropriately designed ligand, given the calculated binding free energy of the water molecule. This approach therefore allows the numerical assessment of whether or not a given water molecule should be targeted for displacement as part of a rational drug design strategy.     

Exploring experimental sources of multiple protein conformations in structure-based drug design

Receptor flexibility must be incorporated into structure-based drug design in order to portray a more accurate representation of a protein in solution. Our approach is to generate pharmacophore models based on multiple conformations of a protein and is very similar to solvent mapping of hot spots. Previously, we had success using computer-generated conformations of apo human immunodeficiency virus-1 protease (HIV-1p). Here, we examine the use of an NMR ensemble versus a collection of crystal structures, and we compare back to our previous study based on computer-generated conformations. To our knowledge, this is the first direct comparison of an NMR ensemble and a collection of crystal structures to incorporate protein flexibility in structure-based drug design. To provide an accurate comparison between the experimental sources, we used bound structures for our multiple protein structure (MPS) pharmacophore models. The models from an NMR ensemble and a collection of crystal structures were both able to discriminate known HIV-1p inhibitors from decoy molecules and displayed superior performance over models created from single conformations of the protein. Although the active-site conformations were already predefined by bound ligands, the use of MPS allows us to overcome the cross-docking problem and generate a model that does not simply reproduce the chemical characteristics of a specific ligand class. We show that there is more structural variation between 28 structures in an NMR ensemble than 90 crystal structures bound to a variety of ligands. MPS models from both sources performed well, but the model determined using the NMR ensemble appeared to be the most general yet accurate representation of the active site. This work encourages the use of NMR models in structure-based design.     

Fatty acid binding sites of human and bovine albumins: Differences observed by spin probe ESR

Bovine and human serum albumins and recombinant human albumin, all non-covalently complexed with 5- and 16-doxyl stearic acids, were investigated by ESR spectroscopy in solution over a range of pH values (5.5–8.0) and temperatures (25–50 °C), with respect to the allocation and mobility of fatty acid (FA) molecules bound to the proteins and conformation of the binding sites. In all proteins bound FA undergo a permanent intra-albumin migration between the binding sites and inter-domain residence. Nature identity of the recombinant human albumin to its serum-derived analog was observed. However, the binding sites of bovine albumin appeared shorter in length and wider in diameter than those of human albumin. Presumably, less tightly folded domains in bovine albumin allow better penetration of water molecules in the interior of the globule that resulted in higher activation energy of FA dissociation from the binding site. Thus, the sensitive technique based on ESR non-covalent spin labeling allowed quantitative analysis and reliable comparison of the fine features of binding proteins.     

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.     

Snooker: a structure-based pharmacophore generation tool applied to class A GPCRs

G-protein coupled receptors (GPCRs) are important drug targets for various diseases and of major interest to pharmaceutical companies. The function of individual members of this protein family can be modulated by the binding of small molecules at the extracellular side of the structurally conserved transmembrane (TM) domain. Here, we present Snooker, a structure-based approach to generate pharmacophore hypotheses for compounds binding to this extracellular side of the TM domain. Snooker does not require knowledge of ligands, is therefore suitable for apo-proteins, and can be applied to all receptors of the GPCR protein family. The method comprises the construction of a homology model of the TM domains and prioritization of residues on the probability of being ligand binding. Subsequently, protein properties are converted to ligand space, and pharmacophore features are generated at positions where protein ligand interactions are likely. Using this semiautomated knowledge-driven bioinformatics approach we have created pharmacophore hypotheses for 15 different GPCRs from several different subfamilies. For the beta-2-adrenergic receptor we show that ligand poses predicted by Snooker pharmacophore hypotheses reproduce literature supported binding modes for ～75% of compounds fulfilling pharmacophore constraints. All 15 pharmacophore hypotheses represent interactions with essential residues for ligand binding as observed in mutagenesis experiments and compound selections based on these hypotheses are shown to be target specific. For 8 out of 15 targets enrichment factors above 10-fold are observed in the top 0.5% ranked compounds in a virtual screen. Additionally, prospectively predicted ligand binding poses in the human dopamine D3 receptor based on Snooker pharmacophores were ranked among the best models in the community wide GPCR dock 2010.     

Mapping the binding site of a large set of quinazoline type EGF-R inhibitors using molecular field analyses and molecular docking studies

In the current work, three-dimensional QSAR studies for one large set of quinazoline type epidermal growth factor receptor (EGF-R) inhibitors were conducted using two types of molecular field analysis techniques: comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). These compounds belonging to six different structural classes were randomly divided into a training set of 122 compounds and a test set of 13 compounds. The statistical results showed that the 3D-QSAR models derived from CoMFA were superior to those generated from CoMSIA. The most optimal CoMFA model after region focusing bears significant cross-validated r(2)(cv) of 0.60 and conventional r(2) of 0.92. The predictive power of the best CoMFA model was further validated by the accurate estimation to these compounds in the external test set, and the mean agreement of experimental and predicted log(IC(50)) values of the inhibitors is 0.6 log unit. Separate CoMFA models were conducted to evaluate the influence of different partial charges (Gasteiger-Marsili, Gasteiger-Hückel, MMFF94, ESP-AM1, and MPA-AM1) on the statistical quality of the models. The resulting CoMFA field map provides information on the geometry of the binding site cavity and the relative weights of various properties in different site pockets for each of the substrates considered. Moreover, in the current work, we applied MD simulations combined with MM/PBSA (Molecular mechanics/Possion-Boltzmann Surface Area) to determine the correct binding mode of the best inhibitor for which no ligand-protein crystal structure was present. To proceed, we define the following procedure: three hundred picosecond molecular dynamics simulations were first performed for the four binding modes suggested by DOCK 4.0 and manual docking, and then MM/PBSA was carried out for the collected snapshots. The most favorable binding mode identified by MM/PBSA has a binding free energy about 10 kcal/mol more favorable than the second best one. The most favorable binding mode identified by MM/PBSA can give satisfactory explanation of the SAR data of the studied molecules and is in good agreement with the contour maps of CoMFA. The most favorable binding mode suggests that with the quinazoline-based inhibitor, the N3 atom is hydrogen-bonded to a water molecule which, in turn, interacts with Thr 766, not Thr 830 as proposed by Wissner et al. (J. Med. Chem. 2000, 43, 3244). The predicted complex structure of quinazoline type inhibitor with EGF-R as well as the pharmacophore mapping from CoMFA can interpret the structure activities of the inhibitors well and afford us important information for structure-based drug design.     

Interaction of bis(1-anilino-8-naphthalenesulfonate) with yeast hexokinase: a steady-state fluorescence study

Bis(1-analino-8-naphthalenesulfonate) (bis-ANS) is a useful probe for hydrophobic areas on protein molecules and it has been proposed that it has a general affinity for the nucleotide binding site(s). There appear to be two different classes of binding sites for bis-ANS on hexokinase and these can be tentatively assigned as primary and secondary binding sites. The rate of binding of bis-ANS at the primary binding site is fast, whereas binding at secondary site(s) is slow. The slow increase in the fluorescence intensity on binding with bis-ANS is not due to conformational change in the enzyme, which may lead to the increase in the quantum yield of the bound dye. Circular dichroism measurements indicate that there is no significant change in the secondary structure on binding with this probe. In the presence of saturating amounts of glucose, the increase in fluorescence intensity due to binding at the secondary binding site(s) is significantly lowered. This indicates that glucose-induced conformational change has been sensed by this probe. From kinetic studies, it has been observed that bis-ANS is an effective competitive inhibitor of yeast hexokinase with respect to ATP. The stoichiometry of binding of this fluorescent probe is about one per subunit at the primary site both in the presence and absence of glucose, and the dissociation constant of bis-ANS is unaffected by glucose. It is possible to decrease significantly the amount of fluorescence intensity at the primary site by nucleotides. These results indicate that bis-ANS interacts at the site where nucleotide interacts. Energy transfer experiments indicate the proximity of some tryptophan(s) and bound bis-ANS molecule(s).     

Heterogeneity and liability of Pb(II) complexation by humic substances: practical interpretation tools

The complexation of Pb(II) by natural organic matter (NOM) is better described by taking into account the dependence of the strength of binding on metal loading conditions. The utility of a linear differential equilibrium function for interpretation of metal ion binding data is demonstrated. This approach considers the binding intensity (log K*) as a function of metal ion loading (theta = bound metal/binding site concentration). Three methods for calculating this function are presented: -direct calculation from metal titration curves, -direct calculation from polarograms, and -compilation of data derived from interpretation of complexation in terms of one- or two- binding sites (e.g. Scatchard analysis), i.e. Cc (complexation capacity = effective site concentration)-K pairs. Heterogeneity also impacts on the apparent lability of complexes; complexes formed at the lowest metal loadings are the least labile.     

Identification of Chemokine Ligands by Biochemical Fragmentation and Simulated Peptide Evolution

Short linear peptides can overcome certain limitations of small molecules for targeting protein–protein interactions (PPIs). Herein, the interaction between the human chemokine CCL19 with chemokine receptor CCR7 was investigated to obtain receptor‐derived CCL19‐binding peptides. After identifying a linear binding site of CCR7, five hexapeptides binding to CCL19 in the low micromolar to nanomolar range were designed, guided by pharmacophore and lipophilicity screening of computationally generated peptide libraries. The results corroborate the applicability of the computational approach and the chosen selection criteria to obtain short linear peptides mimicking a protein–protein interaction site.