Conserved water molecules in a large family of microbial ribonucleases.

We systematically analyzed the crystallographically determined water molecules of all known structures of RNase T1 and compared them to the ordered solvent in a large number of related microbial nucleases. To assess the crystallographers' impact on the interpretation of the solvent structure, we independently refined five validation structures from diffraction data derived from five isomorphous crystals of RNase T1. We also compared the positions of water molecules found in 11 published isomorphous RNase T1 inhibitor complexes. These data suggest that the positions of most of the waters located on the surface of a protein and that are well-determined in the experimental electron density maps are determined primarily by crystal packing forces. Water molecules with less well-defined electron density are in general unique to one or a small number of crystal structures. Only a small number of the well-defined waters are found to be independent of the crystal environment. These waters have a low accessible surface area and B-factor, and tend to be conserved in the crystal structures of a number of evolutionary related ribonucleases as well. A single water molecule is found conserved in all known microbial ribonucleases.     

Modeling the Hydration Layer around Proteins: HyPred

Protein hydration plays an integral role in determining protein function and stability. We develop a simple method with atomic level precision for predicting the solvent density near the surface of a protein. A set of proximal radial distribution functions are defined and calculated for a series of different atom types in proteins using all-atom, explicit solvent molecular dynamic simulations for three globular proteins. A major improvement in predicting the hydration layer is found when the protein is held immobile during the simulations. The distribution functions are used to develop a model for predicting the hydration layer with sub-1-Ångstrom resolution without the need for additional simulations. The model and the distribution functions for a given protein are tested in their ability to reproduce the hydration layer from the simulations for that protein, as well as those for other proteins and for simulations in which the protein atoms are mobile. Predictions for the density of water in the hydration shells are then compared with high occupancy sites observed in crystal structures. The accuracy of both tests demonstrates that the solvation model provides a basis for quantitatively understanding protein solvation and thereby predicting the hydration layer without additional simulations.     

DRoP: Automated detection of conserved solvent-binding sites on proteins

Water and ligand binding play critical roles in the structure and function of proteins, yet their binding sites and significance are difficult to predict a priori. Multiple solvent crystal structures (MSCS) is a method where several X-ray crystal structures are solved, each in a unique solvent environment, with organic molecules that serve as probes of the protein surface for sites evolved to bind ligands, while the first hydration shell is essentially maintained. When superimposed, these structures contain a vast amount of information regarding hot spots of protein-protein or protein-ligand interactions, as well as conserved water-binding sites retained with the change in solvent properties. Optimized mining of this information requires reliable structural data and a consistent, objective analysis tool. Detection of related solvent positions (DRoP) was developed to automatically organize and rank the water or small organic molecule binding sites within a given set of structures. It is a flexible tool that can also be used in conserved water analysis given multiple structures of any protein independent of the MSCS method. The DRoP output is an HTML format list of the solvent sites ordered by conservation rank in its population within the set of structures, along with renumbered and recolored PDB files for visualization and facile analysis. Here, we present a previously unpublished set of MSCS structures of bovine pancreatic ribonuclease A (RNase A) and use it together with published structures to illustrate the capabilities of DRoP.     

Role of invariant water molecules in retaining the active site geometry of β-lactamase: a molecular dynamics simulation study

Water molecules play a critical role in stabilising the three-dimensional architecture, dynamics and function of biological macromolecules. Comparative analysis of structurally similar proteins has shown that there are water molecules conserved in the same relative positions and make similar hydrogen bonds with proteins in all crystal structures. These invariant water molecules are essential for the maintenance of the native structure of proteins. The present study explores the role of invariant water molecules to maintain the active site geometry of β-lactamase enzyme. Thirteen crystal structures of class-A β-lactamase from Staphylococcus aureus have been used in this study. Molecular dynamics simulations of the protein structures were performed in hydrated as well as in dehydrated conditions. The analysis showed that significant changes occur in the active site geometry due to dehydration. These changes can be attributed to the removal of water molecules at the active site.     

Lining the pockets of kinases and phosphatases

The regulation of the activity of kinases and phosphatases is an essential aspect of intracellular signal transduction. Recently determined structures of AGC protein kinases, including isoforms of PKB, PKC, GRK and ROCK, indicate that occupancy of a hydrophobic pocket in the kinase N-lobe by a segment of the protein immediately C terminal to the kinase domain provides a mechanism for regulating kinase activity. In addition, crystal structures of Aurora-A and Aurora-B, which are closely related to AGC family kinases, in complex with their activators, TPX2 and INCENP, respectively, show how allosteric kinase activation is achieved by the binding of the activator protein to an equivalent hydrophobic pocket. Hence, regulation of kinase activity by analogous interactions is a shared regulatory mechanism of these kinases. Two crystal structures have explained the molecular basis of PKA anchoring through its regulatory subunits by members of the AKAP family of scaffold proteins. AKAPs can also interact directly with protein kinase and phosphatase catalytic domains. The crystal structure of the PP1 catalytic subunit in complex with the targeting subunit MYPT1 indicates that there is also scope for intimate phosphatase regulation by scaffold proteins.     

Conservation of solvent-binding sites in 10 crystal forms of T4 lysozyme.

Solvent-binding sites were compared in 10 different crystal forms of phage T4 lysozyme that were refined using data from 2.6 A to 1.7 A resolution. The sample included 18 crystallographically independent lysozyme molecules. Despite different crystallization conditions, variable crystal contacts, changes due to mutation, and varying attention to solvent during crystallographic refinement, 62% of the 20 most frequently occupied sites were conserved. Allowing for potential steric interference from neighboring molecules in the crystal lattice, this fraction increased to 79% of the sites. There was, however, no solvent-binding site that was occupied in all 18 lysozyme molecules. A buried double site was occupied in 17 instances and 2 other internal sites were occupied 15 times. Apart from these buried sites, the most frequently occupied sites were often at the amino-termini of alpha-helices. Solvent molecules at the most conserved sites tended to have crystallographic thermal factors lower than average, but atoms with low B-factors were not restricted to these sites. Although superficial inspection may suggest that only 50-60% (or less) of solvent-binding sites are conserved in different crystal forms of a protein, it appears that many sites appear to be empty either because of steric interference or because the apparent occupancy of a given site can vary from crystal to crystal. The X-ray method of identifying sites is somewhat subjective and tends to result in specification only of those solvent molecules that are well ordered and bound with high occupancy, even though there is clear evidence for solvent bound at many additional sites.(ABSTRACT TRUNCATED AT 250 WORDS)     

Structure of cubic insulin crystals in glucose solutions.

X-ray structures of cubic insulin crystals in high concentrations of glucose at different pH levels and temperatures have been refined to high resolution. We have identified one glucose-binding site near the N-terminus of the A-chain whose occupancy is pH dependent. The effects of reduced water activity on the ordered protein and solvent structures have been examined. Our analysis showed no notable conformational changes in the ordered protein structures or ordered solvent molecules near the protein surface, but the presence of glucose does have a significant effect on the overall density distribution of the bulk solvent in the solvent-accessible volume. We compared the structure of cubic insulin at room temperature and liquid-nitrogen temperature, under identical solvent conditions, using glucose as a cryoprotectant. In this case, we found that the average temperature factor of the protein is reduced and more water molecules can be identified, but there are no significant changes in the protein conformation.     

Role of flexibility and polarity as determinants of the hydration of internal cavities and pockets in proteins

Molecular dynamics simulations of Staphylococcal nuclease and of 10 variants with internal polar or ionizable groups were performed to investigate systematically the molecular determinants of hydration of internal cavities and pockets in proteins. In contrast to apolar cavities in rigid carbon structures, such as nanotubes or buckeyballs, internal cavities in proteins that are large enough to house a few water molecules will most likely be dehydrated unless they contain a source of polarity. The water content in the protein interior can be modulated by the flexibility of protein elements that interact with water, which can impart positional disorder to water molecules, or bias the pattern of internal hydration that is stabilized. This might explain differences in the patterns of hydration observed in crystal structures obtained at cryogenic and room temperature conditions. The ability of molecular dynamics simulations to determine the most likely sites of water binding in internal pockets and cavities depends on its efficiency in sampling the hydration of internal sites and alternative protein and water conformations. This can be enhanced significantly by performing multiple molecular dynamics simulations as well as simulations started from different initial hydration states.     

Deducing hydration sites of a protein from molecular dynamics simulations

Invariant water molecules that are of structural or functional importance to proteins are detected from their presence in the same location in different crystal structures of the same protein or closely related proteins. In this study we have investigated the location of invariant water molecules from MD simulations of ribonuclease A, HIV1-protease and Hen egg white lysozyme. Snapshots of MD trajectories represent the structure of a dynamic protein molecule in a solvated environment as opposed to the static picture provided by crystallography. The MD results are compared to an analysis on crystal structures. A good correlation is observed between the two methods with more than half the hydration sites identified as invariant from crystal structures featuring as invariant in the MD simulations which include most of the functionally or structurally important residues. It is also seen that the propensities of occupying the various hydration sites on a protein for structures obtained from MD and crystallographic studies are different. In general MD simulations can be used to predict invariant hydration sites when there is a paucity of crystallographic data or to complement crystallographic results.     

Water polygons in high‐resolution protein crystal structures

We have analyzed the interstitial water (ISW) structures in 1500 protein crystal structures deposited in the Protein Data Bank that have greater than 1.5 Å resolution with less than 90% sequence similarity with each other. We observed varieties of polygonal water structures composed of three to eight water molecules. These polygons may represent the time‐ and space‐averaged structures of “stable” water oligomers present in liquid water, and their presence as well as relative population may be relevant in understanding physical properties of liquid water at a given temperature. On an average, 13% of ISWs are localized enough to be visible by X‐ray diffraction. Of those, averages of 78% are water molecules in the first water layer on the protein surface. Of the localized ISWs beyond the first layer, almost half of them form water polygons such as trigons, tetragons, as well as expected pentagons, hexagons, higher polygons, partial dodecahedrons, and disordered networks. Most of the octagons and nanogons are formed by fusion of smaller polygons. The trigons are most commonly observed. We suggest that our observation provides an experimental basis for including these water polygon structures in correlating and predicting various water properties in liquid state.     

Structural consequences of effector binding to the T state of aspartate carbamoyltransferase: crystal structures of the unligated and ATP- and CTP-complexed enzymes at 2.6-A resolution

The crystal structure of Escherichia coli aspartate carbamoyltransferase complexed with adenosine 5'-triphosphate (ATP) has been solved by molecular replacement and has been refined to a crystallographic residual of 0.17 at 2.6-A resolution by using the computer program X-PLOR. The unit cell dimensions of this crystal form are a = b = 122.2 A and c = 143.3 A and the space group is P321. Although the c-axis unit cell dimension is approximately 1 A longer than the corresponding dimension of the CTP-ligated P321 crystal form (c = 142.2 A), the ATP-ligated enzyme adopts a T-like quaternary structure. The base moiety of ATP interacts with residues Glu10, Ile12, and Lys60 while the ribose is near Asp19 and Lys60; the triphosphate entity is bound to Lys94, although His20 and Arg96 are nearby. We observe a higher occupancy for ATP in the allosteric site of the R1 regulatory chain in comparison to the occupancy of the R6 allosteric site. These crystallographically independent sites are related by a molecular 2-fold axis. There are other violations of the noncrystallographic symmetry that are similar to those observed in the refined CTP-ligated aspartate carbamoyltransferase structure. These infringements on the molecular symmetry might be the result of intermolecular interactions in the crystal. To ensure the most meaningful comparison with the ATP-ligated structure, we refined the previously reported CTP-bound and unligated structures to crystallographic residuals between 0.17 and 0.18 using X-PLOR. These X-PLOR refined structures are not significantly different from the initial structures that had been crystallographically refined by a restrained least-squares method. After making all possible comparisons between the CTP- and ATP-ligated and the unligated T-state structures, we find that the most significant differences are located at the allosteric sites and in small changes in the quaternary structures. At the allosteric site, the binding of CTP and ATP successively enlarges the nucleotide binding cavity, particularly in the vicinity of the base. The changes in the quaternary structure can be characterized by an increase in the separation of the catalytic trimers by approximately 0.5 A as ATP binds to the unligated T structure. On the basis of these structural studies, we discuss the relationships between the conformational differences in the allosteric site and the small changes in the quaternary structure within the T form to the possible mechanisms for CTP inhibition and ATP activation.     

Mechanism of proton transport in bacteriorhodopsin from crystallographic structures of the K,L, M1, M2, and M2' intermediates of the photocycle

We produced the L intermediate of the photocycle in a bacteriorhodopsin crystal in photo-stationary state at 170 K with red laser illumination at 60% occupancy, and determined its structure to 1.62 A resolution. With this model, high-resolution structural information is available for the initial bacteriorhodopsin, as well as the first five states in the transport cycle. These states involve photo-isomerization of the retinal and its initial configurational changes, deprotonation of the retinal Schiff base and the coupled release of a proton to the extracellular membrane surface, and the switch event that allows reprotonation of the Schiff base from the cytoplasmic side. The six structural models describe the transformations of the retinal and its interaction with water 402, Asp85, and Asp212 in atomic detail, as well as the displacements of functional residues farther from the Schiff base. The changes provide rationales for how relaxation of the distorted retinal causes movements of water and protein atoms that result in vectorial proton transfers to and from the Schiff base.     

Interaction between bound water molecules and local protein structures: A statistical analysis of the hydrogen bond structures around bound water molecules

Water molecules play an important role in protein folding and protein interactions through their structural association with proteins. Examples of such structural association can be found in protein crystal structures, and can often explain protein functionality in the context of structure. We herein report the systematic analysis of the local structures of proteins interacting with water molecules, and the characterization of their geometric features. We first examined the interaction of water molecules with a large local interaction environment by comparing the preference of water molecules in three regions, namely, the protein–protein interaction (PPI) interfaces, the crystal contact (CC) interfaces, and the non‐interfacial regions. High preference of water molecules to the PPI and CC interfaces was found. In addition, the bound water on the PPI interface was more favorably associated with the complex interaction structure, implying that such water‐mediated structures may participate in the shaping of the PPI interface. The pairwise water‐mediated interaction was then investigated, and the water‐mediated residue–residue interaction potential was derived. Subsequently, the types of polar atoms surrounding the water molecules were analyzed, and the preference of the hydrogen bond acceptor was observed. Furthermore, the geometries of the structures interacting with water were analyzed, and it was found that the major structure on the protein surface exhibited planar geometry rather than tetrahedral geometry. Several previously undiscovered characteristics of water–protein interactions were unfolded in this study, and are expected to lead to a better understanding of protein structure and function. Proteins 2016; 84:43–51. © 2015 Wiley Periodicals, Inc.     

The Dynameomics rotamer library: amino acid side chain conformations and dynamics from comprehensive molecular dynamics simulations in water

We have recently completed systematic molecular dynamics simulations of 807 different proteins representing 95% of the known autonomous protein folds in an effort we refer to as Dynameomics. Here we focus on the analysis of side chain conformations and dynamics to create a dynamic rotamer library. Overall this library is derived from 31,000 occurrences of each of 86,217 different residues, or 2.7 × 10(9) rotamers. This dynamic library has 74% overlap of rotamer distributions with rotamer libraries derived from static high-resolution crystal structures. Seventy-five percent of the residues had an assignable primary conformation, and 68% of the residues had at least one significant alternate conformation. The average correlation time for switching between rotamers ranged from 22 ps for Met to over 8 ns for Cys; this time decreased 20-fold on the surface of the protein and modestly for dihedral angles further from the main chain. Side chain S(2) axis order parameters were calculated and they correlated well with those derived from NMR relaxation experiments (R = 0.9). Relationships relating the S(2) axis order parameters to rotamer occupancy were derived. Overall the Dynameomics rotamer library offers a comprehensive depiction of side chain rotamer preferences and dynamics in solution, and more realistic distributions for dynamic proteins in solution at ambient temperature than libraries derived from crystal structures, in particular charged surface residues are better represented. Details of the rotamer library are presented here and the library itself can be downloaded at http://www.dynameomics.org.     

Crystal structures of engrailed homeodomain mutants: implications for stability and dynamics

We report the crystal structures and biophysical characterization of two stabilized mutants of the Drosophila Engrailed homeodomain that have been engineered to minimize electrostatic repulsion. Four independent copies of each mutant occupy the crystal lattice, and comparison of these structures illustrates variation that can be partly ascribed to networks of correlated conformational adjustments. Central to one network is leucine 26 (Leu26), which occupies alternatively two side chain rotameric conformations (-gauche and trans) and different positions within the hydrophobic core. Similar sets of conformational substates are observed in other Engrailed structures and in another homeodomain. The pattern of structural adjustments can account for NMR relaxation data and sequence co-variation networks in the wider homeodomain family. It may also explain the dysfunction associated with a P26L mutation in the human ARX homeodomain protein. Finally, we observe a novel dipolar interaction between a conserved tryptophan and a water molecule positioned along the normal to the indole ring. This interaction may explain the distinctive fluorescent properties of the homeodomain family.     

The role of conserved waters in conformational transitions of Q61H K-ras

To investigate the stability and functional role of long-residence water molecules in the Q61H variant of the signaling protein K-ras, we analyzed all available Ras crystal structures and conformers derived from a series of independent explicit solvent molecular dynamics (MD) simulations totaling 1.76 μs. We show that the protein samples a different region of phase space in the presence and absence of several crystallographically conserved and buried water molecules. The dynamics of these waters is coupled with the local as well as the global motions of the protein, in contrast to less buried waters whose exchange with bulk is only loosely coupled with the motion of loops in their vicinity. Aided by two novel reaction coordinates involving the distance (d) between the C(α) atoms of G60 at switch 2 and G10 at the P-loop and the N-C(α)-C-O dihedral (ξ) of G60, we further show that three water molecules located in lobe1, at the interface between the lobes and at lobe2, are involved in the relative motion of residues at the two lobes of Q61H K-ras. Moreover, a d/ξ plot classifies the available Ras x-ray structures and MD-derived K-ras conformers into active GTP-, intermediate GTP-, inactive GDP-bound, and nucleotide-free conformational states. The population of these states and the transition between them is modulated by water-mediated correlated motions involving the functionally critical switch 2, P-loop and helix 3. These results suggest that water molecules act as allosteric ligands to induce a population shift among distinct switch 2 conformations that differ in effector recognition.     

第一类鱼抗冻蛋白对冰晶生长界面层结构的影响

根据冰晶在水溶液中生长的基本热力学性质,应用多层界面模型,分别得到了冰晶在纯水及抗冻蛋白溶液中生长界面层的吉布斯自由能.由冰晶生长界面层的吉布斯自由能,分析了冰晶在三种不同第一类鱼抗冻蛋白分子溶液中,热平衡状态下生长界面层的微观平衡结构,发现冰晶在抗冻蛋白溶液中生长与其在纯水中生长相比,界面层结构有明显变化,结合抗冻蛋...     

Protein--protein recognition: juxtaposition of domain and interface cores in immunoglobulins and other sandwich-like proteins

Structural analysis of a non-redundant data set of 47 immunoglobulin (Ig) proteins was carried out using a combination of criteria: atom--atom contact compatibility, position occupancy rate, conservation of residue type and positional conservation in 3D space. Our analysis shows that roughly half of the interface positions between the light and heavy chains are specific to individual structures while the other half are conserved across the database. The tendency for conservation of a primary subset of positions holds true for the intra-domain faces as well. These subsets, with an average of 12 conserved positions and a contact surface of 630 A(2), delineate the inter- and intra-domain core, a refined instrument with a reduced target for analysis of sheet--sheet interactions in sandwich-like proteins. Employing this instrument, we find that a majority of Ig interface core positions are adjoined in sequence to domain core positions. This was derived independent of geometric considerations, however beta-sheet side-chain geometry clearly dictates it. The geometric wedding of the domain and interface cores supports the concept of a rigid-like substructure on the protein surface involved in complex formation and indicates a close relationship between surface determinants and those involved in protein folding of Ig domains. The definitions developed for the Ig interface and domain cores proved satisfactory to extract first-approximation cores for a group of 24 non-Ig sandwich-like proteins, treated as individual structures due to their diverse strand topologies. We show that the same rule of positional connectivity between the rigid domain core and interface core extends generally to sandwich-like proteins interacting in a sheet--sheet fashion. The non-Ig structures were used as templates to analyze sandwich-like interfaces of unresolved homologous proteins using a database merging structure and sequence conservation.     

Crystal structures of E. coli laccase CueO at different copper concentrations

CueO protein is a hypothetical bacterial laccase and a good laccase candidate for large scale industrial application. Four CueO crystal structures were determined at different copper concentrations. Low copper occupancy in apo-CueO and slow copper reconstitution process in CueO with exogenous copper were demonstrated. These observations well explain the copper dependence of CueO oxidase activity. Structural comparison between CueO and other three fungal laccase proteins indicates that Glu106 in CueO constitutes the primary counter-work for reconstitution of the trinuclear copper site. Mutation of Glu106 to a Phe enhanced CueO oxidation activity and supported this hypothesis. In addition, an extra alpha-helix from Leu351 to Gly378 covers substrate biding pocket of CueO and might compromises the electron transfer from substrate to type I copper.