Strongly hydrogen-bonded water molecules in the Schiff base region of rhodopsins.

In many rhodopsins, a positively charged retinal chromophore is stabilized by a negatively charged carboxylate, and the presence of bound water molecules has been found in the Schiff base region by X-ray crystallography of various rhodopsins. Low-temperature Fourier-transform infrared (FTIR) spectroscopy can directly monitor hydrogen-bonding alterations of internal water molecules of rhodopsins. In particular, we found that a bridged water molecule between the Schiff base and Asp 85 in bacteriorhodopsin (BR), a light-driven proton-pump protein, forms an extremely strong hydrogen bond. It is likely that a hydration switch of the water from Asp 85 to Asp 212 plays an important role in the proton transfer in the Schiff base region of BR. Comprehensive studies of archaeal and visual rhodopsins have revealed that strongly hydrogen-bonded water molecules are only found in the proteins exhibiting proton-pump activities. Strongly hydrogen-bonded water molecules and its transient weakening may be essential for the proton-pump function of rhodopsins.     

An evaluation of the hydration of lysozyme by an NMR titration method

In this study a new titration method is proposed to study the motional properties of water molecules in conjunction with globular proteins using proton NMR relaxation measurements. The method was applied to the study of the interaction of water with lysozyme and allowed identification of four water fractions-superbound water, polar-bound water, structured water and bulk water - in exchanged equilibrium. The titration demonstrated that 193 water molecules are hydrogen bonded directly to the lysozyme molecule. The combination of structured and bound water extends to 1.4 g H2O per g lysozyme and approx. two to three layers from the surface of the macromolecule. It is proposed that this structured water is related to non-isotropic water rotation in conjunction with hydrophobic patches and directly related to 'hydrophobic bonding' changes. Water amounts greater than 1.4 g H2O per g lysozyme are sufficiently distant from the macromolecule for motion to revert to that typical of water in bulk. The typical correlation times for water motion in the four fraction are: over 10(-6) s (superbound); 10(-9) s (polar bound); 10(-11) s (structured) and 10(-12) s (bulk). These results correlate well with results from other measurement techniques found in the literature.     

Localization of bound water in the solution structure of the immunoglobulin binding domain of streptococcal protein G. Evidence for solvent-induced helical distortion in solution.

The presence of bound water in the solution structure of the IgG binding domain of streptococcal protein G has been investigated by nuclear magnetic resonance using three-dimensional 1H rotating frame Overhauser 1H-15N multiple quantum coherence spectroscopy. The backbone amide protons of three residues, Ala20, Gln32 and Tyr33, are found to be in close proximity to bound water. Examination of the three-dimensional structure of the IgG binding domain indicates that in the vicinity of these three residues there are no backbone groups that do not already participate in hydrogen bonding and there are no suitably placed side-chain groups available for hydrogen bonding with water. As the lifetime of the bound water detected in this nuclear magnetic resonance experiment is greater than about one nanosecond, it is likely that the two bound water molecules participate in a bifurcating hydrogen bonding network comprising a CO-NH hydrogen bonded pair, such that the water molecule accepts a hydrogen bond from the NH proton and donates one to the carbonyl oxygen with the result that the amide proton is involved in a three center hydrogen bond. On the basis of the structure, one water molecule participates in such an interaction with the Ala20(NH)-Met1(CO) hydrogen bonded pair at the beginning of an anti-parallel beta-sheet, and the other with the Tyr33(NH)-Val29(CO) hydrogen bonded pair in the single alpha-helix. The latter, which is external and solvent accessible, is associated with a distortion in the alpha-helix centered around Tyr33 which consists of a significant increase in the CO(i-4)-N(i) and CO(i-4)-NH(i) distances relative to those in the rest of the helix, as well as a significant departure in the phi, psi angles of Tyr33 relative to regular helical geometry. Such solvent induced distortions in alpha-helices have been previously noticed in crystal structures and were postulated as possible folding intermediates for helical structures. The present observation of this phenomenon in solution indicates, however, that these water molecules are tightly bound and represent an integral part of the protein framework.     

Electron structure of plastoquinone and coupling of electron and proton transport in thylakoids of the higher plants

The energy dependence on hydrogen position for a system, consisting of plastoquinone (in different redox states) and histidine molecules was studied. The distance between the atoms forming the hydrogen bond, an oxygen of the quinone molecule and a nitrogen of histidine, was supposed to be fixed. It was shown that for neutral quinone the total energy is minimal when the hydrogen is bound to histidine; for reduced quinone, the probability of hydrogen binding to quinone and histidine is approximately equal (so that a hydrogen bond is formed) and on secondary reduction of plastoquinone, the hydrogen binds to it.     

Hydrogen-bonding alterations of the protonated Schiff base and water molecule in the chloride pump of Natronobacterium pharaonis

Halorhodopsin is a light-driven chloride ion pump. Chloride ion is bound in the Schiff base region of the retinal chromophore, and unidirectional chloride transport is probably enforced by the specific hydrogen-bonding interaction with the protonated Schiff base and internal water molecules. In this article, we study hydrogen-bonding alterations of the Schiff base and water molecules in halorhodopsin of Natronobacterium pharaonis (pHR) by assigning their N-D and O-D stretching vibrations in D(2)O, respectively. Highly accurate low-temperature Fourier transform infrared spectroscopy revealed that hydrogen bonds of the Schiff base and water molecules are weak in the unphotolyzed state, whereas they are strengthened upon retinal photoisomerization. Halide dependence of the stretching vibrations enabled us to conclude that the Schiff base forms a direct hydrogen bond with Cl(-) only in the K intermediate. Hydrogen bond of the Schiff base is further strengthened in the L(1) intermediate, whereas the halide dependence revealed that the acceptor is not Cl(-), but presumably a water molecule. Thus, it is concluded that the hydrogen-bonding interaction between the Schiff base and Cl(-) is not a driving force of the motion of Cl(-). Rather, the removal of its hydrogen bonds with the Schiff base and water(s) makes the environment around Cl(-) less polar in the L(1) intermediate, which presumably drives the motion of Cl(-) from its binding site to the cytoplasmic domain.     

Structure of the yeast cytochrome bc1 complex with a hydroxyquinone anion Qo site inhibitor bound

Bifurcated electron transfer during ubiquinol oxidation is the key reaction of cytochrome bc1 complex catalysis. Binding of the competitive inhibitor 5-n-heptyl-6-hydroxy-4,7-dioxobenzothiazole to the Qo site of the cytochrome bc1 complex from Saccharomyces cerevisiae was analyzed by x-ray crystallography. This alkylhydroxydioxobenzothiazole is bound in its ionized form as evident from the crystal structure and confirmed by spectroscopic analysis, consistent with a measured pKa = 6.1 of the hydroxy group in detergent micelles. Stabilizing forces for the hydroxyquinone anion inhibitor include a polarized hydrogen bond to the iron-sulfur cluster ligand His181 and on-edge interactions via weak hydrogen bonds with cytochrome b residue Tyr279. The hydroxy group of the latter contributes to stabilization of the Rieske protein in the b-position by donating a hydrogen bond. The reported pH dependence of inhibition with lower efficacy at alkaline pH is attributed to the protonation state of His181 with a pKa of 7.5. Glu272, a proposed primary ligand and proton acceptor of ubiquinol, is not bound to the carbonyl group of the hydroxydioxobenzothiazole ring but is rotated out of the binding pocket toward the heme bL propionate A, to which it is hydrogen-bonded via a single water molecule. The observed hydrogen bonding pattern provides experimental evidence for the previously proposed proton exit pathway involving the heme propionate and a chain of water molecules. Binding of the alkyl-6-hydroxy-4,7-dioxobenzothiazole is discussed as resembling an intermediate step of ubiquinol oxidation, supporting a single occupancy model at the Qo site.     

Molecular orbital study on the effect of the binding of water to base pair

The effects of binding water to base pairs was studied by means of the CNDO/2 molecular orbital method. The solvation energy is largest when water is bound as a proton donor and is smallest when it is stacked parallel to the plane of the base pair. The effects of two water molecules are nearly additive. The binding of one water molecule to the adenineuracil pair makes one of the two hydrogen bonds stronger and the other weaker. The change in the hydrogen bonding force is explained in terms of electrostatic and charge transfer energies. By the comparison with the adenine-cytosine pair, it is revealed that the binding of water to adenine serves to yield larger solvation energy for the complementary A-U pair than for the non-complementary A-C pair. It was also observed that the solvation energy due to the binding of water to pyrimidine was larger for A-C than for A-U.     

Refined structure of human carbonic anhydrase II at 2.0 Å resolution

The structure of human erythrocytic carbonic anhydrase II has been refined by constrained and restrained structure–factor least-squares refinement at 2.0 Å resolution. The conventional crystallographic R value is 17.3%. Of 167 solvent molecules associated with the protein, four are buried and stabilize secondary structure elements. The zinc ion is ligated to three histidyl residues and one water molecule in a nearly tetrahedral geometry. In addition to the zinc-bound water, seven more water molecules are identified in the active site. Assuming that Glu-106 is deprotonated at pH 8.5, some of the hydrogen bond donor–acceptor relations in the active site can be assigned and are described here in detail. The Oγ1 atom of Thr-199 donates its proton to the Oε1 atom of Glu-106 and can function as a hydrogen bond acceptor only in additional hydrogen bonds.     

A view of the three dimensional structure of globular proteins

A view of the three dimensional structure of globular proteins based on continuous networks of hydrogen bonds is proposed. Active sites of enzymes and ion sites are prominent and, within the networks, there are islands of hydrophobic regions giving an overall piebald effect to the appearance of the molecule. This point of view was originally suggested by the results of quantum mechanical computations on the coupling between hydrogen bonds. A formalism for the total energy of a globular protein in water is also suggested.The study of five lines of experimental evidence supports this suggestion. The analysis of the experimental X-ray data for ten globular proteins, using the NETWORK program, revealed the existence of these hydrogen bond networks; X-ray data showed that water molecules tend to occupy fixed positions relative to the protein molecule; a survey has shown that water molecules tend to occupy specific positions relative to the hydrogen bonding side chains; experimental evidence on the bulk properties of lysozyme showed that there exist tightly bound water molecules; graphics studies of the ribonucleaseA molecule demonstrated the networks and the piebald effect. This point of view is pictorially simple and, to illustrate the use of such networks, we discuss the simple ion pairs which occur as substructures within the networks.     

Inhibition of Ice Nucleus Growth in Water by Alanine Dipeptide

A molecular dynamics simulation has been carried out for the mixture of an ice nucleus, supercooled water and a molecule of alanine dipeptide (AD). The dipeptide molecule has been allocated near the nucleus surface which corresponds to the prism plane of ice crystal. The molecule is found to approach the ice surface so that the two hydrophilic sites on one side of the molecule (Oc2 and Hn1) are closest to the surface. The hydrogen bond between Hn1 site and the oxygen atom on the prism plane of the ice nucleus is expected. The perturbations of two hydrophilic sites (Oc1 and Hn2), which are surrounded by hydrophobic sites and are pointing away from the surface, attenuate the approach of water molecules to these sites. Thus, these water molecules diffuse. The hydrogen bond between the oxygen atoms on the prism plane and the hydrogen atoms of water molecules is attenuated by the diffusion.     

Structure of bacteriophage T4 lysozyme refined at 1.7 A resolution

The structure of the lysozyme from bacteriophage T4 has been refined at 1.7 A resolution to a crystallographic residual of 19.3%. The final model has bond lengths and bond angles that differ from "ideal" values by 0.019 A and 2.7 degrees, respectively. The crystals are grown from electron-dense phosphate solutions and the use of an appropriate solvent continuum substantially improved the agreement between the observed and calculated structure factors at low resolution. Apart from changes in the conformations of some side-chains, the refinement confirms the structure of the molecule as initially derived from a 2.4 A resolution electron density map. There are 118 well-ordered solvent molecules that are associated with the T4 lysozyme molecule in the crystal. Four of these are more-or-less buried. There is a clustering of water molecules within the active site cleft but, other than this, the solvent molecules are dispersed around the surface of the molecule and do not aggregate into ice-like structures or pentagonal or hexagonal clusters. The apparent motion of T4 lysozyme in the crystal can be interpreted in terms of significant interdomain motion corresponding to an opening and closing of the active site cleft. For the amino-terminal domain the motion can be described equally well (correlation coefficients approx. 0.87) as quasi-rigid-body motion either about a point or about an axis of rotation. The motion in the crystals of the carboxy-terminal domain is best described as rotation about an axis (correlation coefficient 0.80) although in this case the apparent motion seems to be influenced in part by crystal contacts and may be of questionable relevance to dynamics in solution.     

Apparent hydrogen bonding by strongly immobilized spin-labels

The hyperfine separations of nitroxide spin-labels which are tightly bound within hemoglobin exhibit a substantial temperature dependence even when the hemoglobin is immobilized by freezing or precipitation. It is shown that NO.--HX hydrogen bond formation by the spin-label within its binding site is a good explanation for the observed temperature dependence. Comparative studies using different hemoglobin derivatives and two different spin-labels suggest that the HX group may be some element of the protein matrix and that this hydrogen bond may be a factor in the stabilization of the label within its binding site. The hyperfine separation of a fatty acid spin probe incorporated into aqueous bilayer dispersons of dipalmitoylphosphatidylcholine also exhibits a temperature dependence at low temperature which is qualitatively similar to that of the spin-labeled hemoglobin systems. Saturation transfer electron paramagnetic resonance measurements indicate that label motion is not the source of this temperature dependence. A hydrogen-bond equilibrium between water molecules and the nitroxide NO. group appears to be a plausible source of the temperature-dependent hyperfine separation in the lipid bilayer system. Small amplitude torsional oscillation or librational motion by the nitroxide may also produce additional changes in the hyperfine separation which are difficult to distinguish from hydrogen-bonding effects under some circumstances. The apparent hydrogen-bond equilibrium exhibits a strong thermal and environmental dependence which may be of importance in a number of biophysical spin-label measurements.     

Hydrogen bond integrity between MHC class II molecules and bound peptide determines the intracellular fate of MHC class II molecules.

MHC class II molecules associate with peptides through pocket interactions and the formation of hydrogen bonds. The current paradigm suggests that the interaction of side chains of the peptide with pockets in the class II molecule is responsible for the formation of stable class II-peptide complexes. However, recent evidence has shown that the formation of hydrogen bonds between genetically conserved residues of the class II molecule and the main chain of the peptide contributes profoundly to peptide stability. In this study, we have used I-A(k), a class II molecule known to form strong pocket interactions with bound peptides, to probe the general importance of hydrogen bond integrity in peptide acquisition. Our studies have revealed that abolishing hydrogen bonds contributed by positions 81 or 82 in the beta-chain of I-A(k) results in class II molecules that are internally degraded when trafficked through proteolytic endosomal compartments. The presence of high-affinity peptides derived from either endogenous or exogenous sources protects the hydrogen bond-deficient variant from intracellular degradation. Together, these data indicate that disruption of the potential to form a complete hydrogen bond network between MHC class II molecules and bound peptides greatly diminishes the ability of class II molecules to bind peptides. The subsequent failure to stably acquire peptides leads to protease sensitivity of empty class II molecules, and thus to proteolytic degradation before export to the surface of APCs.     

Conserved Water-mediated Hydrogen Bond Network between TM-I, -II, -VI,and -VII in 7TM Receptor Activation

Five highly conserved polar residues connected by a number of structural water molecules together with two rotamer micro-switches, TrpVI:13 and TyrVII:20, constitute an extended hydrogen bond network between the intracellular segments of TM-I, -II, -VI, and -VII of 7TM receptors. Molecular dynamics simulations showed that, although the fewer water molecules in rhodopsin were relatively movable, the hydrogen bond network of the β2-adrenergic receptor was fully loaded with water molecules that were surprisingly immobilized between the two rotamer switches, both apparently being in their closed conformation. Manipulations of the rotamer state of TyrVII:20 and TrpVI:13 demonstrated that these residues served as gates for the water molecules at the intracellular and extracellular ends of the hydrogen bond network, respectively. TrpVI:13 at the bottom of the main ligand-binding pocket was shown to apparently function as a catching trap for water molecules. Mutational analysis of the β2-adrenergic receptor demonstrated that the highly conserved polar residues of the hydrogen bond network were all important for receptor signaling but served different functions, some dampening constitutive activity (AsnI:18, AspII:10, and AsnVII:13), whereas others (AsnVII:12 and AsnVII:16) located one helical turn apart and sharing a water molecule were shown to be essential for agonist-induced signaling. It is concluded that the conserved water hydrogen bond network of 7TM receptors constitutes an extended allosteric interface between the transmembrane segments being of crucial importance for receptor signaling and that part of the function of the rotamer micro-switches, TyrVII:20 and TrpVI:13, is to gate or trap the water molecules.     

The hydration of the neurotransmitter acetylcholine in aqueous solution

Neutron diffraction augmented with hydrogen isotope substitution has been used to examine the water structure around the acetylcholine molecular ion in aqueous solution. It is shown that the nearest-neighbor water molecules in the region around the trimethylammonium headgroup are located either in a ring around the central nitrogen atom or between the carbon atoms, forming a sheath around the onium group. Moreover the water molecules in this cavity do not bond to the onium group but rather form hydrogen bonds with water molecules in the surrounding aqueous environment. Given that in the bound state the onium headgroup must be completely desolvated, the absence of bonding between the onium headgroup and the surrounding water solvent may be selectively favorable to acetylcholine-binding in the receptor site. Away from the headgroup, pronounced hydrogen-bonding of water to the carbonyl oxygen is observed, but not to the ether oxygen in the acetylcholine chain.     

Proton conduction in phosphatidylethanolamine

The dc conductivity of polycrystalline phosphatidylethanolamine (PE) was measured in the temperature range 60–120°C. Since no conclusive evidence had so far been obtained for the presence of proton conduction in this phospholipid, hydrogen gas was shown in the present experiment to evolve during the electrolysis in its premelted state between 91 and 124°C. In this temperature range molecules assume rotation around the molecular axes and proton conduction of the Grotthus type takes place possibly along two chains of intermolecular hydrogen bonds running in parallel. Zwitter-ions behave cooperatively as proton donors and acceptors in transferring proton from molecule to molecule via the hydrogen bond networks. This efficient push-pull way of proton transferring seems to account for the fact that no polarization was observed in the dc conduction experiments. The amount of evolved gas appears to be not exactly in accordance with Faraday's law and discussions are made on possible causes for this slight deviation.     

The role of hydrogen bonding in the enzymatic reaction catalyzed by HIV-1 protease

The hydrogen-bond network in various stages of the enzymatic reaction catalyzed by HIV-1 protease was studied through quantum-classical molecular dynamics simulations. The approximate valence bond method was applied to the active site atoms participating directly in the rearrangement of chemical bonds. The rest of the protein with explicit solvent was treated with a classical molecular mechanics model. Two possible mechanisms were studied, general-acid/general-base (GA/GB) with Asp 25 protonated at the inner oxygen, and a direct nucleophilic attack by Asp 25. Strong hydrogen bonds leading to spontaneous proton transfers were observed in both reaction paths. A single-well hydrogen bond was formed between the peptide nitrogen and outer oxygen of Asp 125. The proton was diffusely distributed with an average central position and transferred back and forth on a picosecond scale. In both mechanisms, this interaction helped change the peptide-bond hybridization, increased the partial charge on peptidyl carbon, and in the GA/GB mechanism, helped deprotonate the water molecule. The inner oxygens of the aspartic dyad formed a low-barrier, but asymmetric hydrogen bond; the proton was not positioned midway and made a slightly elongated covalent bond, transferring from one to the other aspartate. In the GA/GB mechanism both aspartates may help deprotonate the water molecule. We observed the breakage of the peptide bond and found that the protonation of the peptidyl amine group was essential for the peptide-bond cleavage. In studies of the direct nucleophilic mechanism, the peptide carbon of the substrate and oxygen of Asp 25 approached as close as 2.3 A.     

Functional plasticity in the substrate binding site of beta-secretase

The aspartic protease beta-secretase (BACE) cleaves the amyloid precursor protein into a 42 residue beta-peptide, which is the principal biochemical marker of Alzheimer's disease. Multiple explicit-water molecular dynamics simulations of the apo and inhibitor bound structures of BACE indicate that both open- and closed-flap conformations are accessible at room temperature and should be taken into account for inhibitor design. Correlated motion is observed within each of the two lobes of BACE, as well as for the interfacial region. A self-inhibited conformation with the side chain of Tyr71 occupying the S(1) pocket is present in some of the unbound simulations. The reversible loss of the side chain hydrogen bond between the catalytic Asp32 and Ser35, due to the concomitant reorientation of the Ser35 hydroxyl group and a water molecule conserved in pepsin-like enzymes, provides further evidence for the suggestion that Ser35 assists in proton acceptance and release by Asp32 during catalysis.     

Refined crystal structure of cytoplasmic malate dehydrogenase at 2.5-A resolution

The molecular structure of cytoplasmic malate dehydrogenase from pig heart has been refined by alternating rounds of restrained least-squares methods and model readjustment on an interactive graphics system. The resulting structure contains 333 amino acids in each of the two subunits, 2 NAD molecules, 471 solvent molecules, and 2 large noncovalently bound molecules that are assumed to be sulfate ions. The crystallographic study was done on one entire dimer without symmetry restraints. Analysis of the relative position of the two subunits shows that the dimer does not obey exact 2-fold rotational symmetry; instead, the subunits are related by a 173 degrees rotation. The structure results in a R factor of 16.7% for diffraction data between 6.0 and 2.5 A, and the rms deviations from ideal bond lengths and angles are 0.017 A and 2.57 degrees, respectively. The bound coenzyme in addition to hydrophobic interactions makes numerous hydrogen bonds that either are directly between NAD and the enzyme or are with solvent molecules, some of which in turn are hydrogen bonded to the enzyme. The carboxamide group of NAD is hydrogen bonded to the side chain of Asn-130 and via a water molecule to the backbone nitrogens of Leu-157 and Asp-158 and to the carbonyl oxygen of Leu-154. Asn-130 is one of the corner residues in a beta-turn that contains the lone cis peptide bond in cytoplasmic malate dehydrogenase, situated between Asn-130 and Pro-131. The active site histidine, His-186, is hydrogen bonded from nitrogen ND1 to the carboxylate of Asp-158 and from its nitrogen NE2 to the sulfate ion bound in the putative substrate binding site. In addition to interacting with the active site histidine, this sulfate ion is also hydrogen bonded to the guanidinium group of Arg-161, to the carboxamide group of Asn-140, and to the hydroxyl group of Ser-241. It is speculated that the substrate, malate or oxaloacetate, is bound in the sulfate binding site with the substrate 1-carboxyl hydrogen bonded to the guanidinium group of Arg-161.     

Crystal structure of cyclodextrin glucanotransferase from alkalophilic Bacillus sp. 1011 complexed with 1-deoxynojirimycin at 2.0 A resolution

1-Deoxynojirimycin, a pseudo-monosaccharide, is a strong inhibitor of glucoamylase but a relatively weak inhibitor of cyclodextrin glucanotransferase (CGTase). To elucidate this difference, the crystal structure of the CGTase from alkalophilic Bacillus sp. 1011 complexed with 1-deoxynojirimycin was determined at 2.0 A resolution with the crystallographic R value of 0.154 (R(free) = 0.214). The asymmetric unit of the crystal contains two CGTase molecules and each molecule binds two 1-deoxynojirimycins. One 1-deoxynojirimycin molecule is bound to the active center by hydrogen bonds with catalytic residues and water molecules, but its binding mode differs from that expected in the substrate binding. Another 1-deoxynojirimycin found at the maltose-binding site 1 is bound to Asn-667 with a hydrogen bond and by stacking interaction with the indole moiety of Trp-662 of molecule 1 or Trp-616 of molecule 2. Comparison of this structure with that of the acarbose-CGTase complex suggested that the lack of stacking interaction with the aromatic side chain of Tyr-100 is responsible for the weak inhibition by 1-deoxynojirimycin of the enzymatic action of CGTase.