Analysis of solvent structure in proteins using neutron D2O-H2O solvent maps: pattern of primary and secondary hydration of trypsin.

A method of determining the water structure in protein crystals is described using neutron solvent difference maps. These maps are obtained by comparing the changes in diffracted intensities between two data sets, one in which H2O is the major solvent constituent, and a second in which D2O is the solvent medium. To a good first approximation, the protein atom contributions to the scattering intensities in both data sets are equal and cancel, but since H2O and D2O have very different neutron-scattering properties, their differences are accentuated to reveal an accurate representation of the solvent structure. The method also employs a series of density modification steps that impose known physical constraints on the density distribution function in the unit cell by making real space modifications directly to the density maps. Important attributes of the method are that (1) it is less subjective in the assignment of water positions than X-ray analysis; (2) there is threefold improvement in the signal-to-noise ratio for the solvent density; and (3) the iterative density modification produces a low-biased representation of the solvent density. Tests showed that water molecules with as low as 10% occupancy could be confidently assigned. About 300 water sites were assigned for trypsin from the refined solvent density; 140 of these sites were defined in the maps as discrete peaks, while the remaining were found within less-ordered channels of density. There is a very good correspondence between the sites in the primary hydration layer and waters found in the X-ray structure. Most water sites are clustered into H-bonding networks, many of which are found along intermolecular contact zones. The bound water is equally distributed between contacting apolar and polar atoms at the protein interface. A common occurrence at hydrophobic surfaces is that apolar atoms are circumvented by one or more waters that are part of a larger water network. When the effects on surface accessibility by neighboring molecules in the crystal lattice are taken into consideration, only about 29% of the surface does not interface ordered water. About 25% of the ordered water is found in the second hydration sphere. In many instances these waters bridge larger clusters of primary layer waters. It is apparent that, in certain regions of the crystal, the organization of ordered water reflects the characteristics of the crystal environment more than those of trypsin's surface alone.     

Solvent structure in vitamin B12 coenzyme crystals

Both the ordered and disordered solvent networks of vitamin B₁₂ coenzyme crystal hydrate have been generated by Monte Carlo simulation techniques. Several different potential functions have been use to model both water-water and water-solute (i.e., water-coenzyme) interactions. The results have been analysed in terms of the structural properties of the water networks, such as mean water oxygen and hydrogen positions, coordination of each water molecule, and maxima of probability density maps in all four asymmetric units of this crystal.The following results were found: (I) Within each asymmetric unit only one hydrogen bonding network was predicted although there were several hydrogen atom positions for any one solvent molecule (defined as maxima in probability density). (II) Reasonable agreement was obtained between predicted and experimental positions in the ordered solvent region, independent of the potential function used. (III) The positions of the calculated probability density maxima for the disordered channel region were different in different asymmetric units; this led to different simulated hydrogen bond networks which were not always consistent with the experimentally determined alternative (lower occupancy) sites.The results suggest that it is advisable to simulate more than one asymmetric unit if one wishes to look at disorder in the solvent regions. Probability density maps were qualitatively very useful for picturing these disordered regions. However, there were no significant differences between quantitative results predicted using either average atomic positions or maxima of the probability density distributions.Problems in quantifying agreement between experimental and predicted disordered solvent networks are discussed. The potential which included hydrogen atoms explicitly (EMPWI) seemed to give the best overall agreement, mainly because it was successful in predicting the unusually short hydrogen bonds which are found in this crystal.     

Computer simulation of hydration networks around amino acids

The networks of solvent hydrogen bonds around polar and apolar amino acids have been studied by computer simulation techniques using a non-pair additive model for the water molecules interactions. Analysis of the simulated aqueous solutions has shown the presence of water molecules which (a) form a bridge around individual polar solute atoms (self-bridging loops) and (b) form chains between different polar solute atoms (polar bridging chains). Some of these networks associated with polar solute atoms from pentagons but 4, 6 and 7 sided polygons are also seen. The water molecule close to apolar solute atoms (<4.0 Å) also form irregular networks with polygons of 4, 5, 6 and 7 sides. These networks are compared with those found experimentally in ice, clathrates and crystal hydrates of macromolecules.     

Water structure in vitamin B12 coenzyme crystals. II. Structural characteristics of the solvent networks.

The geometrical details of the solvent structure in vitamin B12 coenzyme crystals with respect to hydrogen bonding and nonbonded contacts, are described. The individual H-bond geometries varied over wide ranges, similar to those observed in small molecule structures. Large deviations from tetrahedral coordination were found around a majority of the waters. The mutual positions and orientations of the water molecules could not be adequately explained in terms of the H-bonding relationships present in the structure. However, additional investigations, which focused on the short range nonbonded contacts around water positions in a variety of crystal hydrates, revealed several structural regularities (Savage, 1986b). These features relate to the nonbonded O...O, H...O, and H...H interactions, and give rise to a set of repulsive restrictions that are seen to be very much stronger stereochemical restraints than those associated with H-bonding. The short-range restrictions appear largely to govern the local orientational correlations and packing arrangements of the water structure within the coenzyme (and other hydrate) crystals. In more general terms, the inclusion of the nonbonding relationships as well as the attractive H-bonding interactions, leads to a significant increase in our understanding of water structure(s). The repulsive restrictions can be used as stereochemical restraints in the interpretation and refinement of solvent structures within larger hydrate systems, such as protein crystals. They may also be included in potential functions used to simulate solvent structures in aqueous solutions and hydrate systems.     

Solvent structure in crystals of trypsin determined by X-ray and neutron diffraction.

The solvent structure in orthorhombic crystals of bovine trypsin has been independently determined by X-ray diffraction to 1.35 A resolution and by neutron diffraction to 2.1 A resolution. A consensus model of the water molecule positions was obtained using oxygen positions identified in the electron density map determined by X-ray diffraction, which were verified by comparison to D2O-H2O difference neutron scattering density. Six of 184 water molecules in the X-ray structure, all with B-factors greater than 50 A2, were found to be spurious after comparison with neutron results. Roughly two-thirds of the water of hydration expected from thermodynamic data for proteins was localized by neutron diffraction; approximately one-half of the water of hydration was located by X-ray diffraction. Polar regions of the protein are well hydrated, and significant D2O-H2O difference density is seen for a small number of water molecules in a second shell of hydration. Hydrogen bond lengths and angles calculated from unconstrained refinement of water positions are distributed about values typically seen in small molecule structures. Solvent models found in seven other bovine trypsin and trypsinogen and rat trypsin structures determined by X-ray diffraction were compared. Internal water molecules are well conserved in all trypsin structures including anionic rat trypsin, which is 65% homologous to bovine trypsin. Of the 22 conserved waters in trypsin, 19 were also found in trypsinogen, suggesting that they are located in regions of the apoprotein that are structurally conserved in the transition to the mature protein. Seven waters were displaced upon activation of trypsinogen. Water structure at crystal contacts is not generally conserved in different crystal forms. Three groups of integral structural water molecules are highly conserved in all solvent structures, including a spline of water molecules inserted between two beta-strands, which may resemble an intermediate in the formation of beta sheets during the folding of a protein.     

Fluid bilayer structure determination by the combined use of x-ray and neutron diffraction. II. "Composition-space" refinement method.

This is the second of two papers describing a method for the joint refinement of the structure of fluid bilayers using x-ray and neutron diffraction data. We showed in the first paper (Wiener, M. C., and S. H. White. 1990. Biophys. J. 59:162-173) that fluid bilayers generally consist of a nearly perfect lattice of thermally disordered unit cells and that the canonical resolution d/hmax is a measure of the widths of quasimolecular components represented by simple Gaussian functions. The thermal disorder makes possible a "composition space" representation in which the quasimolecular Gaussian distributions describe the number or probability of occupancy per unit length across the width of the bilayer of each component. This representation permits the joint refinement of neutron and x-ray lamellar diffraction data by means of a single quasimolecular structure that is fit simultaneously to both diffraction data sets. Scaling of each component by the appropriate neutron or x-ray scattering length maps the composition space profile to the appropriate scattering length space for comparison to experimental data. Other extensive properties, such as mass, can also be obtained by an appropriate scaling of the refined composition space structure. Based upon simple bilayer models involving crystal and liquid crystal structural information, we estimate that a fluid bilayer with hmax observed diffraction orders will be accurately represented by a structure with approximately hmax quasimolecular components. Strategies for assignment of quasimolecular components are demonstrated through detailed parsing of a phospholipid molecule based upon the one-dimensional projection of the crystal structure of dimyristoylphosphatidylcholine. Finally, we discuss in detail the number of experimental variables required for the composition space joint refinement. We find fluid bilayer structures to be marginally determined by the experimental data. The analysis of errors, which takes on particular importance under these circumstances, is also discussed.     

Optimization of linear disorder predictors yields tight association between crystallographic disorder and hydrophobicity

X-ray crystallographic protein structures often contain disordered regions that are observed as missing electron density. Diffraction data may give little or no direct evidence as to the specific nature of disordered regions. We have developed a weighted window-based disorder predictor optimized using crystallographic data. Performance of a predictor is strongly influenced by chain termini. Optimized score adjustment values for amino- and carboxy-terminal positions demonstrate a simple, monotonic relationship between disorder and residue distance from termini. This optimized disorder predictor performs similarly to DISOPRED2 on crystallographically disordered regions. Data-optimized residue disorder propensities show strong linear correlation with experimentally determined amino acid transfer energies between water and hydrogen-bonding organic solvents, which primarily reflect residue hydrophobicity (exemplified by the Nozaki-Tanford hydrophobicity scale). Disorder propensities do not correlate as well with transfer energies between water and apolar solvents, which primarily reflect a different hydropathic property: residue hydrophilicity (also reflected by the Kyte-Doolittle hydropathy scale). Our results suggest that while hydrophobic side-chain interactions are primarily involved in determining stability of the folded conformation, hydrogen bonding, and similar polar interactions are primarily involved in conformational and interaction specificity.     

Residence times of water molecules in the hydration sites of myoglobin

Hydration sites are high-density regions in the three-dimensional time-averaged solvent structure in molecular dynamics simulations and diffraction experiments. In a simulation of sperm whale myoglobin, we found 294 such high-density regions. Their positions appear to agree reasonably well with the distributions of waters of hydration found in 38 x-ray and 1 neutron high-resolution structures of this protein. The hydration sites are characterized by an average occupancy and a combination of residence time parameters designed to approximate a distribution of residence times. It appears that although the occupancy and residence times of the majority of sites are rather bulk-like, the residence time distribution is shifted toward the longer components, relative to bulk. The sites with particularly long residence times are located only in the cavities and clefts of the protein. This indicates that other factors, such as hydrogen bonds and hydrophobicity of underlying protein residues, play a lesser role in determining the residence times of the longest-lived sites.     

Specific assay for radiolabelled digoxin and its known apolar metabolites in biological fluids. I.

A specific assay method for radiolabelled digoxin and its known apolar metabolites in plasma, urine and saliva was developed. The assay permits the delineation of the pharmacokinetics of digoxin and its metabolites after single-dose administration of the drug to humans. Column chromatographic and solvent extraction procedures were used for the separation of apolar and polar compounds. Thin-layer chromatography was applied for the individual and specific assessment of digoxin and its apolar metabolites. Apolar and polar standards were used for quantitative assessments of all the procedures used. Accuracy and precision of the assay developed were evaluated in plasma, urine and saliva using biological samples spiked with known amounts of standards and by measuring replicates of biological samples obtained from pharmacokinetic studies with digoxin administration to humans.     

A review about nothing: Are apolar cavities in proteins really empty?

Cavities within proteins that are strictly apolar typically appear to be empty. It has been suggested, however, that water molecules may be present within such cavities but are too disordered to be seen in conventional crystallographic analyses. In contrast, it is argued here that solvent mobility will be limited by the size of the cavity and for this reason high‐occupancy solvent in cavities of typical volume should be readily detectable using X‐ray crystallography. Recent experimental studies of cavity hydration are reviewed. Such studies are consistent with theoretical predictions that it is energetically unfavorable to have a single water molecule in an apolar cavity. As apolar cavities become larger, a point is reached where it is favorable to have the cavity occupied by a cluster of mutually H‐bonded water molecules. The exact size of such a cavity in a protein is yet to be verified.     

Structural and thermodynamic analysis of the binding of solvent at internal sites in T4 lysozyme

To investigate the structural and thermodynamic basis of the binding of solvent at internal sites within proteins a number of mutations were constructed in T4 lysozyme. Some of these were designed to introduce new solvent-binding sites. Others were intended to displace solvent from preexisting sites. In one case Val-149 was replaced with alanine, serine, cysteine, threonine, isoleucine, and glycine. Crystallographic analysis shows that, with the exception of isoleucine, each of these substitutions results in the binding of solvent at a polar site that is sterically blocked in the wild-type enzyme. Mutations designed to perturb or displace a solvent molecule present in the native enzyme included the replacement of Thr-152 with alanine, serine, cysteine, valine, and isoleucine. Although the solvent molecule was moved in some cases by up to 1.7 A, in no case was it completely removed from the folded protein. The results suggest that hydrogen bonds from the protein to bound solvent are energy neutral. The binding of solvent to internal sites within proteins also appears to be energy neutral except insofar as the bound solvent may prevent a loss of energy due to potential hydrogen bonding groups that would otherwise be unsatisfied. The introduction of a solvent-binding site appears to require not only a cavity to accommodate the water molecule but also the presence of polar groups to help satisfy its hydrogen-bonding potential. It may be easier to design a site to accommodate two or more water molecules rather than one as the solvent molecules can then hydrogen-bond to each other. For similar reasons it is often difficult to design a point mutation that will displace a single solvent molecule from the core of a protein.     

Conformational characteristics of peptides and unanticipated results from crystal structure analyses

Preferred conformation and types of molecular folding are some of the topics that can be addressed by structure analysis using x-ray diffraction of single crystals. The conformations of small linear peptide molecules with 2-6 residues are affected by polarity of solvent, presence of water molecules, hydrogen bonding with neighboring molecules, and other packing forces. Larger peptides, both cyclic and linear, have many intramolecular hydrogen bonds, the effect of which outweighs any intermolecular attractions. Numerous polymorphs of decapeptides grown from a variety of solvents, with different cocrystallized solvents, show a constant conformation for each peptide. Large conformational changes occur, however, upon complexation with metal ions. A new form of free valinomycin grown from DMSO exhibits near three-fold symmetry with only three intramolecular hydrogen bonds. The peptide is in the form of a shallow bowl with a hydrophobic exterior. Near the bottom of the interior of the bowl are three carbonyl oxygens, spaced and directed so that they are in position to form three ligands to a K+, e.g., complexation can be completed by the three lobes containing the beta-bends closing over and encapsulating the K+ ion. In another example, free antamanide and the biologically inactive perhydro analogue, in which four phenyl groups become cyclic hexyl groups, have essentially the same folding of backbone and side chains. The conformation changes drastically upon complexation with Li+ or Na+. However, the metal ion complex of natural antamanide has a hydrophobic globlar form whereas the metal ion complex of the inactive perhydro analogue has a polar band around the middle. The structure results indicate that the antamanide molecule is in a complexed form during its biological activity. Single crystal x-ray diffraction structure analyses have identified the manner in which water molecules are essential to creating minipolar areas on apolar helices. Completely apolar peptides, such as membrane-active peptides, can acquire amphiphilic character by insertion of a water molecule into the helical backbone of Boc-Aib-Ala-Leu-Aib-Ala-Leu-Aib-Ala-Leu-Aib-OMe, for example. The C-terminal half assumes an alpha-helix conformation, whereas the N-terminal half is distorted by an insertion of a water molecule W(1) between N(Ala5) and O(Ala2), forming hydrogen bonds N(5)H...W(1) and W(1)...O(2). The distortion of the helix exposes C = O(Aib1) and C = O(Aib4) to the outside environment with the consequence of attracting additional water molecules. The leucyl side chains are on the other side of the molecule. Thus a helix with an apolar sequence can mimic an amphiphilic helix.     

Explicit solvent models in protein pKa calculations.

Continuum methods for calculation of protein electrostatics treat buried and ordered water molecules by one of two approximations; either the dielectric constant of regions containing ordered water molecules is equal to the bulk solvent dielectric constant, or it is equal to the protein dielectric constant though no fixed atoms are used to represent water molecules. A method for calculating the titration behavior of individual residues in proteins has been tested on models of hen egg white lysozyme containing various numbers of explicit water molecules. Water molecules were included based on hydrogen bonding, solvent accessibility, and/or proximity to titrating groups in the protein. Inclusion of water molecules significantly alters the calculated titration behavior of individual titrating sites, shifting calculated pKa values by up to 0.5 pH unit. Our results suggest that approximately one water molecule within hydrogen-bonding distance of each charged group should be included in protein electrostatics calculations.     

Refinement of the structure of recombinant rat intestinal fatty acid-binding apoprotein at 1.2-A resolution.

The three-dimensional structure of the 131-residue rat intestinal fatty acid-binding protein, without bound ligand (apoI-FABP), has been refined with x-ray diffraction data to a nominal resolution of 1.19 A. The final model has a conventional crystallographic R-factor of 16.9% for 34,290 unique reflections [a root mean square (r.m.s.) deviation for bond length of 0.012 A and a r.m.s. deviation of 2.368 degrees for bond angles]. Ninety-two residues are present as components of the protein's 10 anti-parallel beta-strands while 14 residues are part of its two short alpha-helices. The beta-strands and alpha-helices are organized into two nearly orthogonal beta-sheets. Particular attention has been placed in defining solvent structure and the structures of discretely disordered groups in this protein. Two hundred thirty-seven solvent molecules have been identified; 24 are located within apoI-FABP. The refined model includes alternate conformers for 228 protein atoms (109 main-chain, 119 side-chain) and 63 solvent molecules. We have found several aromatic side-chains with multiple conformations located near, or in, the protein's ligand binding site. This observation, along with the fact that these side-chains have a temperature factor that is relatively higher than that of other aromatic residues, suggests that they may be involved in the process of noncovalent binding of fatty acid. The absence of a true hydrophobic core in I-FABP suggests that its structural integrity may be maintained primarily by a hydrogen bonding network involving protein and solvent atoms.     

Structure of a ternary complex of an allosteric lactate dehydrogenase from Bacillus stearothermophilus at 2.5 A resolution.

We report the refined structure of a ternary complex of an allosterically activated lactate dehydrogenase, including the important active site loop. Eightfold non-crystallographic symmetry averaging was utilized to improve the density maps. Interactions between the protein and bound coenzyme and oxamate are described in relation to other studies using site-specific mutagenesis. Fructose 1,6-bisphosphate (FruP2) is bound to the enzyme across one of the 2-fold axes of the tetramer, with the two phosphate moieties interacting with two anion binding sites, one on each of two subunits, across this interface. However, because FruP2 binds at this special site, yet does not possess an internal 2-fold symmetry axis, the ligand is statistically disordered and binds to each site in two different orientations. Binding of FruP2 to the tetramer is signalled to the active site principally through two interactions with His188 and Arg173. His188 is connected to His195 (which binds the carbonyl group of the substrate) and Arg173 is connected to Arg171 (the residue that binds the carboxylate group of the substrate).     

On the dependence of molecular conformation on the type of solvent environment: a molecular dynamics study of cyclosporin A

The dependence of the conformation of cyclosporin A (CPA), a cyclic undecapeptide with potent immunosuppressive activity, on the type of solvent environment is examined using the computer simulation method of molecular dynamics (MD). Conformational and dynamic properties of CPA in aqueous solution are obtained from MD simulations of a CPA molecule dissolved in a box with water molecules. Corresponding properties of CPA in apolar solution are obtained from MD simulations of CPA in a box with carbontetrachloride. The results of these simulations in H2O and in CCl4 are compared to each other and to those of previous simulations of crystalline CPA and of an isolated CPA molecule. The conformation of the backbone of the cyclic polypeptide is basically independent of the type of solvent. In aqueous solution the beta-pleated sheet is slightly weaker and the gamma-turn is a bit less pronounced than in apolar solution. Side chains may adopt different conformations in different solvents. In apolar solution the hydrophobic side chain of the MeBmt residue is in an extended conformation with its hydroxyl group hydrogen bonded to the backbone carbonyl group. In aqueous solution this hydrophobic side chain folds over the core of the molecule and the mentioned hydrogen bond is broken in favor of hydrogen bonding to water molecules. The conformation obtained from the MD simulation in CCl4 nicely agrees with experimental atom-atom distance data as obtained from nmr experiments in chloroform. In aqueous solution the relaxation of atomic motion tends to be slower than in apolar solution.     

Mapping the polarity profiles of general anesthetic target sites using n-alkane-(alpha, omega)-diols

The effects of the homologous series of n-alkane-(alpha, omega)-diols have been studied on the inhibition of the purified firefly luciferase enzyme from Photinus pyralis, the inhibition of the purified bacterial luciferase enzyme from Vibrio harveyi, and the induction of general anesthesia in Xenopus laevis tadpoles. All but one of the diols tested were found to be reversible general anesthetics. The diols inhibited firefly luciferase by competing with its normal substrate firefly luciferin, and they inhibited bacterial luciferase by competing with the substrate n-decanal. For all but the smallest agent (1,4-butanediol), only a single diol molecule was found to be involved in the inhibition of the enzymes. Inhibition constants Ki were determined for the enzymes, and general anesthetic EC50 concentrations were determined for tadpoles. These data were then used in conjunction with previously determined n-alkane and n-alcohol data to calculate, as a function of chain length, the incremental standard Gibbs free energies delta (delta G0) for adding apolar -CH2- groups and for converting apolar terminal -CH3 groups to polar -CH2OH groups. The resulting plots of delta (delta G0) versus chain length gave a consistent mapping of the polarity profiles of the anesthetic-binding pockets. They clearly reveal the existence of two substantial and distinct polar regions in the anesthetic-binding pocket of firefly luciferase but only one such region for bacterial luciferase and for the unknown target sites underlying general anesthesia. The polarities and geometric properties of these different binding sites for straight-chain anesthetics are discussed in terms of simple models.     

Crystal structure of low humidity tetragonal lysozyme at 2.1-A resolution. Variability in hydration shell and its structural consequences

Tetragonal crystals of hen egg white lysozyme undergo a reversible transformation, accompanied by loss of water, when the relative humidity of the environment is reduced to about 90%. The structure of the low humidity form has been analyzed, using x-ray data collected at 88% relative humidity, in order to explore the variability in protein hydration caused by a change in the amount of water surrounding the protein molecule and the consequent conformational perturbations in the molecule. The structure has been refined by the restrained least-squares method to an R value of 0.162 for 6269 observed reflections in the 10-2.1-A resolution shell. The refined structure provides interesting examples for the variability in helical parameters, the role of interactions involving side chains and water in the stabilization of secondary structural features, and favorable specific hydration sites. The protein molecule as a whole moves slightly in the low humidity form from its position in the native crystals. The hydration shell tends to move along with the protein. Significant changes, however, occur in the hydration shell. These changes cause structural perturbations in the enzyme molecule, which are most pronounced in regions involved in substrate binding.     

Refined structure of c-phycocyanin from the cyanobacterium Synechococcus vulcanus at 1.6 A: insights into the role of solvent molecules in thermal stability and co-factor structure

The crystal structure of the light-harvesting phycobiliprotein, c-phycocyanin from the thermophilic cyanobacterium Synechococcus vulcanus has been refined to 1.6 A resolution based on the previously determined lower resolution structure (PDB entry 1I7Y). The improved data was collected using synchrotron radiation at 100 K. The significantly improved crystallographic data has lead to improved calculated electron density maps, allowing the unambiguous positioning of all protein and co-factor atoms and the positioning of 377 solvent molecules. The positions of solvent molecules at specific sites important for stabilization of different levels of self-assembly of the phycobilisome structure were identified and the bonding network is described. The presence of solvent molecules in the vicinity of the co-factors and in intermolecular spaces is identified and their possible roles are suggested. All three of the phycocyanobilin co-factors bind water molecules at specific sites between the propionic acid side chains. Molecular dynamic (MD) simulations support that these special waters have a role in stabilization of this conformation. On the basis of the crystal packing reported here and in comparison to other phycobiliprotein crystal forms, we have analyzed the roles of specific sites on the formation of the phycobilisome complex.     

Studies in lipid histochemistry. XIII. The OPA (osmiumtetroxide-periodic acid-alpha-naphthylamine) method for the detection of apolar lipids.

A new procedure for the detection of apolar lipids is described. It is a modification of the OTAN method (Adams, 1959) using periodic acid which oxidatively removes lower osmium derivatives from polar sites only, leaving those in apolar lipids intact and demonstrable with alpha-naphthylamine. Control steps for the exclusion of the possible interference of some less polar complex lipids and of lipopigments are described. The described technic is superior to the conventionally used sudan dyes due partly to the fact that only aqueous solutions are employed thus excluding any extraction of lipids, partly to the more distinct coloration.