Ordered water molecules as key allosteric mediators in a cooperative dimeric hemoglobin

One of the most remarkable structural aspects of Scapharca dimeric hemoglobin is the disruption of a very well-ordered water cluster at the subunit interface upon ligand binding. We have explored the role of these crystallographically observed water molecules by site-directed mutagenesis and osmotic stress techniques. The isosteric mutation of Thr-72 → Val in the interface increases oxygen affinity more than 40-fold with a surprising enhancement of cooperativity. The only significant structural effect of this mutation is to destabilize two ordered water molecules in the deoxy interface. Wild-type Scapharca hemoglobin is strongly sensitive to osmotic conditions. Upon addition of glycerol, striking changes in Raman spectrum of the deoxy form are observed that indicate a transition toward the liganded form. Increased osmotic pressure, which lowers the oxygen affinity in human hemoglobin, raises the oxygen affinity of Scapharca hemoglobin regardless of whether the solute is glycerol, glucose, or sucrose. Analysis of these results provides an estimate of six water molecules lost upon oxygen binding to the dimer, in good agreement with eight predicted from crystal structures. These experiments suggest that the observed cluster of interfacial water molecules plays a crucial role in communication between subunits.

Cotransport of water by the Na+/glucose cotransporter

Water is transported across epithelial membranes in the absence of any hydrostatic or osmotic gradients. A prime example is the small intestine, where 10 liters of water are absorbed each day. Although water absorption is secondary to active solute transport, the coupling mechanism between solute and water flow is not understood. We have tested the hypothesis that water transport is directly linked to solute transport by cotransport proteins such as the brush border Na+/glucose cotransporter. The Na+/glucose cotransporter was expressed in Xenopus oocytes, and the changes in cell volume were measured under sugar-transporting and nontransporting conditions. We demonstrate that 260 water molecules are directly coupled to each sugar molecule transported and estimate that in the human intestine this accounts for 5 liters of water absorption per day. Other animal and plant cotransporters such as the Na+/Cl−/γ-aminobutyric acid, Na+/iodide and H+/amino acid transporters are also able to transport water and this suggests that cotransporters play an important role in water homeostasis.

Local osmotic gradients drive the water flux associated with Na+/glucose cotransport

It recently was proposed [Loo, D. D. F., Zeuthen, T., Chandy, G. & Wright, E. M. (1996) Proc. Natl. Acad. Sci. USA 93, 13367–13370] that SGLT1, the high affinity intestinal and renal sodium/glucose cotransporter carries water molecules along with the cosubstrates with a strict stoichiometry of two Na+, one glucose, and ≈220 water molecules per transport cycle. Using electrophysiology together with sensitive volumetric measurements, we investigated the nature of the driving force behind the cotransporter-mediated water flux. The osmotic water permeability of oocytes expressing human SGLT1 (Lp ± SE) averaged 3.8 ± 0.3 × 10−4 cm⋅s−1 (n = 15) and addition of 100 μM phlorizin (a specific SGLT1 inhibitor) reduced the permeability to 2.2 ± 0.2 × 10−4 cm⋅s−1 (n = 15), confirming the presence of a significant water permeability closely associated with the cotransporter. Addition of 5 mM α-methyl-glucose (αMG) induced an average inward current of 800 ± 10 nA at −50 mV and a water influx reaching 120 ± 20 pL cm−2 ⋅s−1 within 5–8 min. After rapidly inhibiting the Na+/glucose cotransport with phlorizin, the water flux remained significantly elevated, clearly indicating the presence of a local osmotic gradient (Δπ) estimated at 16 ± 2 mOsm. In short-term experiments, a rapid depolarization from −100 to 0 mV in the presence of αMG decreased the cotransport current by 94% but failed to produce a comparable reduction in the swelling rate. A mathematical model depicting the intracellular accumulation of transported osmolytes can accurately account for these observations. It is concluded that, in SGLT1-expressing oocytes, αMG-dependent water influx is induced by a local osmotic gradient by using both endogenous and SGLT1-dependent water permeability.

Large vortex-like structure of dipole field in computer models of liquid water and dipole-bridge between biomolecules

We propose a framework to describe the cooperative orientational motions of water molecules in liquid water and around solute molecules in water solutions. From molecular dynamics (MD) simulation a new quantity “site-dipole field” is defined as the averaged orientation of water molecules that pass through each spatial position. In the site-dipole field of bulk water we found large vortex-like structures of more than 10 Å in size. Such coherent patterns persist more than 300 ps although the orientational memory of individual molecules is quickly lost. A 1-ns MD simulation of systems consisting of two amino acids shows that the fluctuations of site-dipole field of solvent are pinned around the amino acids, resulting in a stable dipole-bridge between side-chains of amino acids. The dipole-bridge is significantly formed even for the side-chain separation of 14 Å, which corresponds to five layers of water. The way that dipole-bridge forms sensitively depends on the side-chain orientations and thereby explains the specificity in the solvent-mediated interactions between biomolecules.

Packing at the protein-water interface.

We have determined the packing efficiency at the protein-water interface by calculating the volumes of atoms on the protein surface and nearby water molecules in 22 crystal structures. We find that an atom on the protein surface occupies, on average, a volume approximately 7% larger than an atom of equivalent chemical type in the protein core. In these calculations, larger volumes result from voids between atoms and thus imply a looser or less efficient packing. We further find that the volumes of individual atoms are not related to their chemical type but rather to their structural location. More exposed atoms have larger volumes. Moreover, the packing around atoms in locally concave, grooved regions of protein surfaces is looser than that around atoms in locally convex, ridge regions. This as a direct manifestation of surface curvature-dependent hydration. The net volume increase for atoms on the protein surface is compensated by volume decreases in water molecules near the surface. These waters occupy volumes smaller than those in the bulk solvent by up to 20%; the precise amount of this decrease is directly related to the extent of contact with the protein.

Significance of bound water to local chain conformations in protein crystals.

We examine how the polypeptide chain in protein crystal structures exploits the multivalent hydrogen-bonding potential of bound water molecules. This shows that multiple interactions with a single water molecule tend to occur locally along the chain. A distinctive internal-coordinate representation of the local water-binding segments reveals several consensus conformations. The fractional water occupancy of each was found by comparison of the total number of conformations in the database regardless of the presence or absence of bound water. The water molecule appears particularly frequently in type II beta-turn geometries and an N-terminal helix feature. This work constitutes a first step into assessing not only the generality but also the significance of specific water binding in globular proteins.

Detection of one slowly exchanging substrate water molecule in the S3 state of photosystem II.

The exchangeability of the substrate water molecules at the catalytic site of water oxidation in photosystem II has been probed by isotope-exchange measurements using mass spectrometric detection of flash-induced oxygen evolution. A stirred sample chamber was constructed to reduce the lag time between injection of H2(18)O and the detecting flash by a factor of more than 1000 compared to the original experiments by R. Radmer and O. Ollinger [(1986) FEBS Lett. 195, 285-289]. Our data show that there is a slow (t1/2 approximately 500 ms, 10 degrees C) and a fast (t1/2 <25 ms, 10 degrees C) exchanging substrate water molecule in the S3 state of photosystem II. The slow exchange is coupled with an activation energy of about 75 kJ/mol and is discussed in terms of a terminal manganese oxo ligand, while the faster exchanging substrate molecule may represent a water molecule not directly bound to the manganese center.

Mapping strain exerted on blood vessel walls using deuterium double-quantum-filtered MRI

A technique is described for displaying distinct tissue layers of large blood vessel walls as well as measuring their mechanical strain. The technique is based on deuterium double-quantum-filtered (DQF) spectroscopic imaging. The effectiveness of the double-quantum filtration in suppressing the signal of bulk water is demonstrated on a phantom consisting of rat tail tendon fibers. Only intrafibrillar water is displayed, excluding all other signals of water molecules that reorient isotropically. One- and two-dimensional spectroscopic imaging of bovine aorta and coronary arteries show the characteristic DQF spectrum of each of the tissue layers. This property is used to obtain separate images of the outer layer, the tunica adventitia, or the intermediate layer, the tunica media, or both. To visualize the effect of elongation, the average residual quadrupole splitting <Δνq> is calculated for each pixel. Two-dimensional deuterium quadrupolar splitting images are obtained for a fully relaxed and a 55% elongated sample of bovine coronary artery. These images indicate that the strong effect of strain is associated with water molecules in the tunica adventitia whereas the DQF NMR signal of water in the tunica media is apparently strain-insensitive. After appropriate calibration, these average quadrupolar splitting images can be interpreted as strain maps.

Structural basis for chemical inhibition of human blood coagulation factor Xa

Factor Xa, the converting enzyme of prothrombin to thrombin, has emerged as an alternative (to thrombin) target for drug discovery for thromboembolic diseases. An inhibitor has been synthesized and the crystal structure of the complex between Des[1–44] factor Xa and the inhibitor has been determined by crystallographic methods in two different crystal forms to 2.3- and 2.4-Å resolution. The racemic mixture of inhibitor FX-2212, (2RS)-(3′-amidino-3-biphenylyl)-5-(4-pyridylamino)pentanoic acid, inhibits factor Xa activity by 50% at 272 nM in vitro. The S-isomer of FX-2212 (FX-2212a) was found to bind to the active site of factor Xa in both crystal forms. The biphenylamidine of FX-2212a occupies the S1-pocket, and the pyridine ring makes hydrophobic interactions with the factor Xa aryl-binding site. Several water molecules meditate inhibitor binding to residues in the active site. In contrast to the earlier crystal structures of factor Xa, such as those of apo-Des[1–45] factor Xa and Des[1–44] factor Xa in complex with a naphthyl inhibitor DX-9065a, two epidermal growth factor-like domains of factor Xa are well ordered in both our crystal forms as well as the region between the two domains, which recently was found to be the binding site of the effector cell protease receptor-1. This structure provides a basis for designing next generation inhibitors of factor Xa.

In disperse solution, “osmotic stress” is a restricted case of preferential interactions

In the practice of “osmotic stress,” the effect of excluded cosolvents on a biochemical equilibrium is interpreted as the number of water molecules participating in the reaction. This action is attributed to lowering of solvent water activity by the cosolvent. This concept of osmotic stress in disperse solution is erroneous: (i) A cosolvent cannot be both excluded and inert, i.e., noninteracting, because exclusion requires a positive free energy change; (ii) a decrease in water activity alone by addition of solute cannot affect an equilibrium when the reacting surface is in contact with the solvent; and (iii) osmotic stress in disperse solution is a restricted case of preferential interactions; the reaction is driven by the free energy of cosolvent exclusion, and the derived number of water molecules is solely a measure of the mutual perturbations of the chemical potentials of the cosolvent and the protein.

Crystal structure of 2,5-diketo-d-gluconic acid reductase A complexed with NADPH at 2.1-Å resolution

The three-dimensional structure of Corynebacterium 2,5-diketo-d-gluconic acid reductase A (2,5-DKGR A; EC 1.1.1.-), in complex with cofactor NADPH, has been solved by using x-ray crystallographic data to 2.1-Å resolution. This enzyme catalyzes stereospecific reduction of 2,5-diketo-d-gluconate (2,5-DKG) to 2-keto-l-gulonate. Thus the three-dimensional structure has now been solved for a prokaryotic example of the aldo–keto reductase superfamily. The details of the binding of the NADPH cofactor help to explain why 2,5-DKGR exhibits lower binding affinity for cofactor than the related human aldose reductase does. Furthermore, changes in the local loop structure near the cofactor suggest that 2,5-DKGR will not exhibit the biphasic cofactor binding characteristics observed in aldose reductase. Although the crystal structure does not include substrate, the two ordered water molecules present within the substrate-binding pocket are postulated to provide positional landmarks for the substrate 5-keto and 4-hydroxyl groups. The structural basis for several previously described active-site mutants of 2,5-DKGR A is also proposed. Recent research efforts have described a novel approach to the synthesis of l-ascorbate (vitamin C) by using a genetically engineered microorganism that is capable of synthesizing 2,5-DKG from glucose and subsequently is transformed with the gene for 2,5-DKGR. These modifications create a microorganism capable of direct production of 2-keto-l-gulonate from d-glucose, and the gulonate can subsequently be converted into vitamin C. In economic terms, vitamin C is the single most important specialty chemical manufactured in the world. Understanding the structural determinants of specificity, catalysis, and stability for 2,5-DKGR A is of substantial commercial interest.

Structure at 2.7 Å resolution of the Paracoccus denitrificans two-subunit cytochrome c oxidase complexed with an antibody FV fragment

The aa3 type cytochrome c oxidase consisting of the core subunits I and II only was isolated from the soil bacterium Paracoccus denitrificans and crystallized as complex with a monoclonal antibody Fv fragment. Crystals could be grown in the presence of a number of different nonionic detergents. However, only undecyl-β-d-maltoside and cyclohexyl-hexyl-β-d-maltoside yielded well-ordered crystals suitable for high resolution x-ray crystallographic studies. The crystals belong to space group P212121 and diffract x-rays to at least 2.5 Å (1 Å = 0.1 nm) resolution using synchrotron radiation. The structure was determined to a resolution of 2.7 Å using molecular replacement and refined to a crystallographic R-factor of 20.5% (Rfree = 25.9%). The refined model includes subunits I and II and the 2 chains of the Fv fragment, 2 heme A molecules, 3 copper atoms, and 1 Mg/Mn atom, a new metal (Ca) binding site, 52 tentatively identified water molecules, and 9 detergent molecules. Only four of the water molecules are located in the cytoplasmic half of cytochrome c oxidase. Most of them are near the interface of subunits I and II. Several waters form a hydrogen-bonded cluster, including the heme propionates and the Mg/Mn binding site. The Fv fragment binds to the periplasmic polar domain of subunit II and is critically involved in the formation of the crystal lattice. The crystallization procedure is well reproducible and will allow for the analysis of the structures of mechanistically interesting mutant cytochrome c oxidases.

Molecular switch in signal transduction: Reaction paths of the conformational changes in ras p21

Conformational changes in ras p21 triggered by the hydrolysis of GTP play an essential role in the signal transduction pathway. The path for the conformational change is determined by molecular dynamics simulation with a holonomic constraint directing the system from the known GTP-bound structure (with the γ-phosphate removed) to the GDP-bound structure. The simulation is done with a shell of water molecules surrounding the protein. In the switch I region, the side chain of Tyr-32, which undergoes a large displacement, moves through the space between loop 2 and the rest of the protein, rather than on the outside of the protein. As a result, the charged residues Glu-31 and Asp-33, which interact with Raf in the homologous RafRBD–Raps complex, remain exposed during the transition. In the switch II region, the conformational changes of α2 and loop 4 are strongly coupled. A transient hydrogen bonding complex between Arg-68 and Tyr-71 in the switch II region and Glu-37 in switch I region stabilizes the intermediate conformation of α2 and facilitates the unwinding of a helical turn of α2 (residues 66–69), which in turn permits the larger scale motion of loop 4. Hydrogen bond exchange between the protein and solvent molecules is found to be important in the transition. Possible functional implications of the results are discussed.

Streaming potential measurements in αβγ-rat epithelial Na+ channel in planar lipid bilayers

Streaming potentials across cloned epithelial Na+ channels (ENaC) incorporated into planar lipid bilayers were measured. We found that the establishment of an osmotic pressure gradient (Δπ) across a channel-containing membrane mimicked the activation effects of a hydrostatic pressure differential (ΔP) on αβγ-rENaC, although with a quantitative difference in the magnitude of the driving forces. Moreover, the imposition of a Δπ negates channel activation by ΔP when the Δπ was directed against ΔP. A streaming potential of 2.0 ± 0.7 mV was measured across αβγ-rat ENaC (rENaC)-containing bilayers at 100 mM symmetrical [Na+] in the presence of a 2 Osmol/kg sucrose gradient. Assuming single file movement of ions and water within the conduction pathway, we conclude that between two and three water molecules are translocated together with a single Na+ ion. A minimal effective pore diameter of 3 Å that could accommodate two water molecules even in single file is in contrast with the 2-Å diameter predicted from the selectivity properties of αβγ-rENaC. The fact that activation of αβγ-rENaC by ΔP can be reproduced by the imposition of Δπ suggests that water movement through the channel is also an important determinant of channel activity.

Proton uptake by bacterial reaction centers: The protein complex responds in a similar manner to the reduction of either quinone acceptor

In bacterial photosynthetic reaction centers, the protonation events associated with the different reduction states of the two quinone molecules constitute intrinsic probes of both the electrostatic interactions and the different kinetic events occurring within the protein in response to the light-generated introduction of a charge. The kinetics and stoichiometries of proton uptake on formation of the primary semiquinone QA− and the secondary acceptor QB− after the first and second flashes have been measured, at pH 7.5, in reaction centers from genetically modified strains and from the wild type. The modified strains are mutated at the L212Glu and/or at the L213Asp sites near QB; some of them carry additional mutations distant from the quinone sites (M231Arg → Leu, M43Asn → Asp, M5Asn → Asp) that compensate for the loss of L213Asp. Our data show that the mutations perturb the response of the protein system to the formation of a semiquinone, how distant compensatory mutations can restore the normal response, and the activity of a tyrosine residue (M247Ala → Tyr) in increasing and accelerating proton uptake. The data demonstrate a direct correlation between the kinetic events of proton uptake that are observed with the formation of either QA− or QB−, suggesting that the same residues respond to the generation of either semiquinone species. Therefore, the efficiency of transferring the first proton to QB is evident from examination of the pattern of H+/QA− proton uptake. This delocalized response of the protein complex to the introduction of a charge is coordinated by an interactive network that links the Q− species, polarizable residues, and numerous water molecules that are located in this region of the reaction center structure. This could be a general property of transmembrane redox proteins that couple electron transfer to proton uptake/release reactions.