An analysis of the relationship between hydration and protein-DNA interactions.

Eleven protein-DNA crystal structures were analyzed to test the hypothesis that hydration sites predicted in the first hydration shell of DNA mark the positions where protein residues hydrogen-bond to DNA. For nine of those structures, protein atoms, which form hydrogen bonds to DNA bases, were found within 1.5 A of the predicted hydration positions in 86% of the interactions. The correspondence of the predicted hydration sites with the hydrogen-bonded protein side chains was significantly higher for bases inside the conserved DNA recognition sequences than outside those regions. In two CAP-DNA complexes, predicted base hydration sites correctly marked 71% (within 1.5 A) of protein atoms, which form hydrogen bonds to DNA bases. Phosphate hydration was compared to actual protein binding sites in one CAP-DNA complex with 78% marked contacts within 2.0 A. These data suggest that hydration sites mark the binding sites at protein-DNA interfaces.     

Phosphate group binding "cup" of PLP-dependent and non-PLP-dependent enzymes: leitmotif and variations

Pyridoxal-5'-phosphate (PLP) is widely used by many enzymes in reactions where amino acids are interconverted. Whereas the role of the pyridoxal ring in catalysis is well understood, the functional role of the single phosphate group in PLP has been less studied. Here we construct unambiguous connection diagrams that describe the interactions among the three non-ester phosphate oxygen atoms of PLP and surrounding atoms from the protein binding site and from water molecules, the so-called phosphate group binding "cup". These diagrams provide a simple means to identify common recognition motifs for the phosphate group in both similar and different protein folds. Diagrams were constructed and compared in the cases of five newly determined structures of PLP-dependent transferases (fold type I enzymes) and, additionally, two non-PLP protein complexes (indole-3-glycerol phosphate synthase (IGPS) with bound indole-3-glycerol phosphate (IGP) and old yellow enzyme (OYE) with bound flavin mononucleotide (FMN)). A detailed comparison of the diagrams shows that three positions out of ten in the structure of the phosphate group binding "cup" contain invariant atoms, while seven others are occupied by conserved atom types. This level of similarity was also observed in the fold type III (TIM beta/alpha-barrel) enzymes that bind three different ligands: PLP, IGP and FMN.     

Identification of the conserved spatial position of key active-site atoms in glycoside hydrolase 13 family members

A computational study on the glycoside hydrolase 13 (GH13) family of the CAZy database has been carried out at the atomic level in order to identify the conserved positions that may be responsible for recognition of the substrate. Analysis with substrate analog-, inhibitor-, or product-bound 3D structures was carried out to find the atomic spatial arrangement of the amino acids that make −2, −1, +1, and +2 subsites and water oxygen atoms around the ligand. The identified conserved positions of subsites were independent from the nature of the amino acid. The −1 and +1 subsites have more conserved positions than the −2 and +2 subsites. Some of the clusters of the −1 and +1 subsites have atoms of the same chemical nature. A spatially conserved position for water, which is stabilized by a hydrogen bond with the carboxyl group of a proton donor (Glu) and Asp of the catalytic triad, was found in the −1 subsite of 75% of the enzymes subjected to analysis. This position could be the region of hydrolytic water.     

Conserved water molecules in X-ray structures highlight the role of water in intramolecular and intermolecular interactions

Water molecules immobilized on a protein or DNA surface are known to play an important role in intramolecular and intermolecular interactions. Comparative analysis of related three-dimensional (3D) structures allows to predict the locations of such water molecules on the protein surface. We have developed and implemented the algorithm WLAKE detecting "conserved" water molecules, i.e. those located in almost the same positions in a set of superimposed structures of related proteins or macromolecular complexes. The problem is reduced to finding maximal cliques in a certain graph. Despite exponential algorithm complexity, the program works appropriately fast for dozens of superimposed structures. WLAKE was used to predict functionally significant water molecules in enzyme active sites (transketolases) as well as in intermolecular (ETS-DNA complexes) and intramolecular (thiol-disulfide interchange protein) interactions. The program is available online at http://monkey.belozersky.msu.ru/~evgeniy/wLake/wLake.html.     

Continuum electrostatic analysis of preferred solvation sites around proteins in solution

To understand water-protein interactions in solution, the electrostatic field is calculated by solving the Poisson-Boltzmann equation, and the free energy surface of water is mapped by translating and rotating an explicit water molecule around the protein. The calculation is applied to T4 lysozyme with data available on the conservation of solvent binding sites in 18 crystallographically independent molecules. The free energy maps around the ordered water sites provide information on the relationship between water positions in crystal structure and in solution. Results show that almost all conserved sites and the majority of nonconserved sites are within 1.3 A of local free energy minima. This finding is in sharp contrast to the behavior of randomly placed water molecules in the boundary layer, which, on the average, must travel more than 3 A to the nearest free energy minimum. Thus, the solvation sites are at least partially determined by protein-water interactions rather than by crystal packing alone. The characteristic water residence times, obtained from the free energies at the local minima, are in good agreement with nuclear magnetic resonance experiments. Only about half of the potential sites show up as ordered water in the 1.7 A resolution X-ray structure. Crystal packing interactions can stabilize weak or mobile potential sites (in fact, some ordered water positions are not close to free energy minima) or can prevent water from occupying certain sites. Apart from a few buried water molecules that are strong binders, the free energies are not very different for conserved and nonconserved sites. We show that conservation of a water site between two crystals occurs if the positions of protein atoms, primarily contributing to the free energy at the local minimum, do not substantially change from one structure to the other. This requirement can be correlated with the nature of the side chain contacting the water molecule in the site.     

A quantitative analysis of interfacial amino acid conservation in protein-protein hetero complexes

A long-standing question in molecular biology is whether interfaces of protein-protein complexes are more conserved than the rest of the protein surfaces. Although it has been reported that conservation can be used as an indicator for predicting interaction sites on proteins, there are recent reports stating that the interface regions are only slightly more conserved than the rest of the protein surfaces, with conservation signals not being statistically significant enough for predicting protein-protein binding sites. In order to properly address these controversial reports we have studied a set of 28 well resolved hetero complex structures of proteins that consists of transient and non-transient complexes. The surface positions were classified into four conservation classes and the conservation index of the surface positions was quantitatively analyzed. The results indicate that the surface density of highly conserved positions is significantly higher in the protein-protein interface regions compared with the other regions of the protein surface. However, the average conservation index of the patches in the interface region is not significantly higher compared with other surface regions of the protein structures. This finding demonstrates that the number of conserved residue positions is a more appropriate indicator for predicting protein-protein binding sites than the average conservation index in the interacting region. We have further validated our findings on a set of 59 benchmark complex structures. Furthermore, an analysis of 19 complexes of antigen-antibody interactions shows that there is no conservation of amino acid positions in the interacting regions of these complexes, as expected, with the variable region of the immunoglobulins interacting mostly with the antigens. Interestingly, antigen interacting regions also have a higher number of non-conserved residue positions in the interacting region than the rest of the protein surface.     

Atomic and residue hydrophilicity in the context of folded protein structures

Water-protein interactions drive protein folding, stabilize the folded structure, and influence molecular recognition and catalysis. We analyzed the closest protein contacts of 10,837 water molecules in crystallographic structures to define a specific hydrophilicity scale reflecting specific rather than bulk solvent interactions. The tendencies of different atom and residue types to be the nearest protein neighbors of bound water molecules correlated with other hydrophobicity scales, verified the relevance of crystallographically determined water positions, and provided a direct experimental measure of water affinity in the context of the folded protein. This specific hydrophilicity was highly correlated with hydrogen-bonding capacity, and correlated better with experimental than computationally derived measures of partitioning between aqueous and organic phases. Atoms with related chemistry clustered with respect to the number of bound water molecules. Neutral and negatively charged oxygen atoms were the most hydrophilic, followed by positively-charged then neutral nitrogen atoms, followed by carbon and sulfur atoms. Agreement between observed side-chain specific hydrophilicity values and values derived from the atomic hydrophilicity scale showed that hydrophilicity values can be synthesized for different functional groups, such as unusual side or main chains, discontinuous epitopes, and drug molecules. Two methods of atomic hydrophilicity analysis provided a measure of complementarity in the interfaces of trypsin:pancreatic trypsin inhibitor and HIV protease:U-75875 inhibitor complexes. © 1995 Wiley-Liss, Inc.     

A fragment library based on Gaussian mixtures predicting favorable molecular interactions

Here, a protein atom-ligand fragment interaction library is described. The library is based on experimentally solved structures of protein-ligand and protein-protein complexes deposited in the Protein Data Bank (PDB) and it is able to characterize binding sites given a ligand structure suitable for a protein. A set of 30 ligand fragment types were defined to include three or more atoms in order to unambiguously define a frame of reference for interactions of ligand atoms with their receptor proteins. Interactions between ligand fragments and 24 classes of protein target atoms plus a water oxygen atom were collected and segregated according to type. The spatial distributions of individual fragment - target atom pairs were visually inspected in order to obtain rough-grained constraints on the interaction volumes. Data fulfilling these constraints were given as input to an iterative expectation-maximization algorithm that produces as output maximum likelihood estimates of the parameters of the finite Gaussian mixture models. Concepts of statistical pattern recognition and the resulting mixture model densities are used (i) to predict the detailed interactions between Chlorella virus DNA ligase and the adenine ring of its ligand and (ii) to evaluate the "error" in prediction for both the training and validation sets of protein-ligand interaction found in the PDB. These analyses demonstrate that this approach can successfully narrow down the possibilities for both the interacting protein atom type and its location relative to a ligand fragment.     

Analysis of structural water and CH···π interactions in HIV-1 protease and PTP1B complexes using a hydrogen bond prediction tool,HBPredicT

A hydrogen bond prediction tool HBPredicT is developed for detecting structural water molecules and CH···π interactions in PDB files of protein-ligand complexes. The program adds the missing hydrogen atoms to the protein, ligands, and oxygen atoms of water molecules and subsequently all the hydrogen bonds in the complex are located using specific geometrical criteria. Hydrogen bonds are classified into various types based on (i) donor and acceptor atoms, and interactions such as (ii) protein-protein, (iii) protein-ligand, (iv) protein-water, (v) ligand-water, (vi) water-water, and (vii) protein-water-ligand. Using the information in category (vii), the water molecules which form hydrogen bonds with the ligand and the protein simultaneously–the structural water–is identified and retrieved along with the associated ligand and protein residues. For CH···π interactions, the relevant portions of the corresponding structures are also extracted in the output. The application potential of this program is tested using 19 HIV-1 protease and 11 PTP1B inhibitor complexes. All the systems showed presence of structural water molecules and in several cases, the CH···π interaction between ligand and protein are detected. A rare occurrence of CH···π interactions emanating from both faces of a phenyl ring of the inhibitor is identified in HIV-1 protease 1D4L.     

Identification of RNA-recognizing modules on the surface of ribosomal proteins

Properties of specific interaction between ribosomal proteins and ribosomal RNAs were analyzed and a method for determination of "recognizing modules" on the protein surface was proposed. The method is based on the search of protein atoms making conserved H-bonds with RNA and forming an invariant spatial structure in homologous rRNA-protein complexes and in the isolated protein. A potential of the method is demonstrated on the determination of the recognizing modules on the surfaces of ribosomal proteins S8, S15 and L5.     

Functional attributes of the phosphate group binding cup of pyridoxal phosphate-dependent enzymes

Twenty-four structures of pyridoxal-5'-phosphate (PLP)-dependent enzymes that represent five different folds are shown to share a common recognition pattern for the phosphate group of their PLP-ligands. All atoms that interact with the phosphate group of PLP in these proteins are organized within a two-layer structure so that the first interacting layer contains from five to seven atoms and parallel with this is a second layer containing from three to seven interacting atoms. In order to identify features of the phosphate-binding site common to PLP-dependent enzymes, a simple procedure is described that assigns relative positions to all interacting atoms unambiguously, such that the networks of interactions for different proteins can be compared. On the basis of these diagrams for 24 enzyme-cofactor complexes, a detailed comparison of the two-layer structures of PLP-dependent enzymes, with both similar and different folds, was made. A majority of the structurally defined PLP-dependent proteins use the same atom types in analogous "key" positions to bind their PLP-ligands. In some instances, proteins use water molecules when a key position is unoccupied. A similar two-layer recognition pattern extends to protein recognition of at least one other, non-PLP ligand, glucosamine 6-phosphate. We refer to this three-dimensional recognition pattern as the phosphate-binding cup. In general, the phosphate-binding cup provides a very stable anchoring point for PLP. When numerous water molecules occur within the cup, however, then the phosphate group of PLP participates directly in the enzymatic reactions with inorganic phosphate replacing the water molecules of the cup. With glucosamine-6-phosphate synthase, the water molecules of the phosphate-binding cup facilitate the entry of substrate and the exit of product.     

Neutron crystallographic studies of T4 lysozyme at cryogenic temperature

Bacteriophage T4 lysozyme (T4L) has been used as a paradigm for seminal biophysical studies on protein structure, dynamics, and stability. Approximately 700 mutants of this protein and their respective complexes have been characterized by X‐ray crystallography; however, despite the high resolution diffraction limits attained in several studies, no hydrogen atoms were reported being visualized in the electron density maps. To address this, a 2.2 Å‐resolution neutron data set was collected at 80 K from a crystal of perdeuterated T4L pseudo‐wild type. We describe a near complete atomic structure of T4L, which includes the positions of 1737 hydrogen atoms determined by neutron crystallography. The cryogenic neutron model reveals explicit detail of the hydrogen bonding interactions in the protein, in addition to the protonation states of several important residues.     

Geometric analysis and comparison of protein-DNA interfaces: why is there no simple code for recognition?

Structural studies of protein-DNA complexes have shown that there are many distinct families of DNA-binding proteins, and have shown that there is no simple "code" describing side-chain/base interactions. However, systematic analysis and comparison of protein-DNA complexes has been complicated by the diversity of observed contacts, the sheer number of complexes currently available and the absence of any consistent method of comparison that retains detailed structural information about the protein-DNA interface. To address these problems, we have developed geometric methods for characterizing the local structural environment in which particular side-chain/base interactions are observed. In particular, we develop methods for analyzing and comparing spatial relationships at the protein-DNA interface. Our method involves attaching local coordinate systems to the DNA bases and to the C(alpha) atoms of the peptide backbone (these are relatively rigid structural units). We use these tools to consider how the position and orientation of the polypeptide backbone (with respect to the DNA) helps to determine what contacts are possible at any given position in a protein-DNA complex. Here, we focus on base contacts that are made in the major groove, and we use spatial relationships in analyzing: (i) the observed patterns of side-chain/base interactions; (ii) observed helix docking orientations; (iii) family/subfamily relationships among DNA-binding proteins; and (iv) broader questions about evolution, altered specificity mutants and the limits for the design of new DNA-binding proteins. Our analysis, which highlights differences in spatial relationships in different complexes and at different positions in a complex, helps explain why there is no simple, general code for protein-DNA recognition.     

Knowledge based modelling of homologous proteins, Part II: Rules for the conformations of substituted sidechains

This paper describes a rapid, automated procedure which can be used for model building sidechains using (i) spatial information from sidechains in topologically equivalent positions as far as such a correlation is observed, and then (ii) most probable conformations of the sidechains in the respective secondary structure type. Analysis of topologically equivalent residues in the structurally conserved regions of a family of proteins implies that the spatial positions of the atoms in the sidechains rather than conformations should be considered when model building. Rules for the modelling of all 20 side-chains from each other in alpha-helical, beta-sheet and loop regions--a total of 1200--are established. Cluster analysis is used on positional data from the sidechain atoms of structurally equivalent residues in an homologous family to guide modelling. The most probable conformation for the sidechain is used for modelling atoms where no useful guidance is obtainable from equivalent sidechains of the homologous proteins. In order to test the procedure we have modelled the sidechains of the residues in the structurally conserved regions of myoglobin from four other globins. The automated procedure described here has been incorporated into the program COMPOSER.     

The carbohydrate-binding site in galectin-3 is preorganized to recognize a sugarlike framework of oxygens: ultra-high-resolution structures and water dynamics

The recognition of carbohydrates by proteins is a fundamental aspect of communication within and between living cells. Understanding the molecular basis of carbohydrate-protein interactions is a prerequisite for the rational design of synthetic ligands. Here we report the high- to ultra-high-resolution crystal structures of the carbohydrate recognition domain of galectin-3 (Gal3C) in the ligand-free state (1.08 ? at 100 K, 1.25 ? at 298 K) and in complex with lactose (0.86 ?) or glycerol (0.9 ?). These structures reveal striking similarities in the positions of water and carbohydrate oxygen atoms in all three states, indicating that the binding site of Gal3C is preorganized to coordinate oxygen atoms in an arrangement that is nearly optimal for the recognition of β-galactosides. Deuterium nuclear magnetic resonance (NMR) relaxation dispersion experiments and molecular dynamics simulations demonstrate that all water molecules in the lactose-binding site exchange with bulk water on a time scale of nanoseconds or shorter. Nevertheless, molecular dynamics simulations identify transient water binding at sites that agree well with those observed by crystallography, indicating that the energy landscape of the binding site is maintained in solution. All heavy atoms of glycerol are positioned like the corresponding atoms of lactose in the Gal3C complexes. However, binding of glycerol to Gal3C is insignificant in solution at room temperature, as monitored by NMR spectroscopy or isothermal titration calorimetry under conditions where lactose binding is readily detected. These observations make a case for protein cryo-crystallography as a valuable screening method in fragment-based drug discovery and further suggest that identification of water sites might inform inhibitor design.     

Synthesis and Th1-immunostimulatory activity of α-galactosylceramide analogues bearing a halogen-containing or selenium-containing acyl chain

A novel series of CD1d ligand α-galactosylceramides (α-GalCers) were synthesized by incorporation of the heavy atoms Br and Se in the acyl chain backbone of α-galactosyl-N-cerotoylphytosphingosine. The synthetic analogues are potent CD1d ligands and stimulate mouse invariant natural killer T (iNKT) cells to selectively enhance Th1 cytokine production. These synthetic analogues would be efficient X-ray crystallographic probes to disclose precise atomic positions of alkyl carbons and lipid–protein interactions in KRN7000/CD1d complexes.     

Conserved water molecules in MHC class-I molecules and their putative structural and functional roles

A set of conserved water positions making direct contacts with the alpha1 and alpha2 domains of the MHC class-I protein was identified by a cluster analysis in 12 high-resolution crystal structures of proteins from different allele types and different species, comprising human, mouse and rat. The analysis revealed a total of 63 clusters, corresponding to water molecules, whose positions are conserved in half or more of the analyzed structures. Analysis of these clusters shows that the most conserved water positions-those appearing in the largest fraction of the structures-were also the most accurately defined, as measured by their normalized crystallographic B-factor. Not too surprisingly, these positions displayed better overlap and formed more H-bonds with the protein. In a second part of this work, a detailed analysis is presented of three of the most conserved water positions and their putative structural and functional roles are discussed. The most highly conserved of the three appears to play an important role in stabilizing the conformation of a twisted beta-turn between residues 118 and 122 (numbering of HLA-B3501, PDB code 1A1N). An equivalent water molecule was found to be associated with a similar beta-turn in 43 unrelated structures surveyed in the PDB, leading to the suggestion that this water molecule plays an important structural role in this type of turn. The second water molecule makes hydrogen bonds with residues lining pocket B in the peptide-binding groove and is suggested to play a role in modulating peptide recognition. The third highly conserved water molecule is located at the first kink of the alpha2 helix, possibly playing a role in determining the position of the N-terminal segment of that helix, which also carries side chains in contact with the bound peptide. This information on conserved water positions in MHC class-I molecules should be helpful in modeling interactions with bound peptide antigens and in designing new peptides with tailor-made affinities.     

High resolution crystal structures of HIV-1 protease with a potent non-peptide inhibitor (UIC-94017) active against multi-drug-resistant clinical strains

The compound UIC-94017 (TMC-114) is a second-generation HIV protease inhibitor with improved pharmacokinetics that is chemically related to the clinical inhibitor amprenavir. UIC-94017 is a broad-spectrum potent inhibitor active against HIV-1 clinical isolates with minimal cytotoxicity. We have determined the high-resolution crystal structures of UIC-94017 in complexes with wild-type HIV-1 protease (PR) and mutant proteases PR(V82A) and PR(I84V) that are common in drug-resistant HIV. The structures were refined at resolutions of 1.10-1.53A. The crystal structures of PR and PR(I84V) with UIC-94017 ternary complexes show that the inhibitor binds to the protease in two overlapping positions, while the PR(V82A) complex had one ordered inhibitor. In all three structures, UIC-94017 forms hydrogen bonds with the conserved main-chain atoms of Asp29 and Asp30 of the protease. These interactions are proposed to be critical for the potency of this compound against HIV isolates that are resistant to multiple protease inhibitors. Other small differences were observed in the interactions of the mutants with UIC-94017 as compared to PR. PR(V82A) showed differences in the position of the main-chain atoms of residue 82 compared to PR structure that better accommodated the inhibitor. Finally, the 1.10A resolution structure of PR(V82A) with UIC-94017 showed an unusual distribution of electron density for the catalytic aspartate residues, which is discussed in relation to the reaction mechanism.     

Synthesis and characterization of novel Cu(II)-bipyridine complexes having functional groups and their application toward molecular recognition

Two kinds of Cu(II) complexes having 2,2′-bipyridine derivatives with two 1-naphthoylamide groups or two ethyl dimethylmalonylamide moieties at 6 and 6′ positions as ligands were prepared and characterized by X-ray crystallography and spectroscopic methods. Those ligands bound to the Cu(II) centers in a tetradentate fashion including two amide oxygen atoms in the equatorial planes. Those complexes were found to recognize carboxylic acids as guest molecules by coordination and additional non-covalent interactions, including intramolecular π-π interactions or hydrogen bonding.