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.

Structural details of an interaction between cardiolipin and an integral membrane protein

Anionic lipids play a variety of key roles in biomembrane function, including providing the immediate environment for the integral membrane proteins that catalyze photosynthetic and respiratory energy transduction. Little is known about the molecular basis of these lipid–protein interactions. In this study, x-ray crystallography has been used to examine the structural details of an interaction between cardiolipin and the photoreaction center, a key light-driven electron transfer protein complex found in the cytoplasmic membrane of photosynthetic bacteria. X-ray diffraction data collected over the resolution range 30.0–2.1 Å show that binding of the lipid to the protein involves a combination of ionic interactions between the protein and the lipid headgroup and van der Waals interactions between the lipid tails and the electroneutral intramembrane surface of the protein. In the headgroup region, ionic interactions involve polar groups of a number of residues, the protein backbone, and bound water molecules. The lipid tails sit along largely hydrophobic grooves in the irregular surface of the protein. In addition to providing new information on the immediate lipid environment of a key integral membrane protein, this study provides the first, to our knowledge, high-resolution x-ray crystal structure for cardiolipin. The possible significance of this interaction between an integral membrane protein and cardiolipin is considered.

Real-time visualization of intracellular hydrodynamics in single living cells

Intracellular water concentrations in single living cells were visualized by nonlinear coherent anti-Stokes Raman scattering (CARS) microscopy. In combination with isotopic exchange measurements, CARS microscopy allowed the real-time observation of transient intracellular hydrodynamics at a high spatial resolution. Studies of the hydrodynamics in the microorganism Dictyostelium discoideum indicated the presence of a microscopic region near the plasma membrane where the mobility of water molecules is severely restricted. Modeling the transient hydrodynamics eventuated in the determination of cell-specific cytosolic diffusion and plasma membrane permeability constants. Our experiments demonstrate that CARS microscopy offers an invaluable tool for probing single-cell water dynamics.

Detection of alkali metal ions in DNA crystals using state-of-the-art X-ray diffraction experiments

The observation of light metal ions in nucleic acids crystals is generally a fortuitous event. Sodium ions in particular are notoriously difficult to detect because their X-ray scattering contributions are virtually identical to those of water and Na+…O distances are only slightly shorter than strong hydrogen bonds between well-ordered water molecules. We demonstrate here that replacement of Na+ by K+, Rb+ or Cs+ and precise measurements of anomalous differences in intensities provide a particularly sensitive method for detecting alkali metal ion-binding sites in nucleic acid crystals. Not only can alkali metal ions be readily located in such structures, but the presence of Rb+ or Cs+ also allows structure determination by the single wavelength anomalous diffraction technique. Besides allowing identification of high occupancy binding sites, the combination of high resolution and anomalous diffraction data established here can also pinpoint binding sites that feature only partial occupancy. Conversely, high resolution of the data alone does not necessarily allow differentiation between water and partially ordered metal ions, as demonstrated with the crystal structure of a DNA duplex determined to a resolution of 0.6 Å.

Surface geochemistry of the clay minerals

Clay minerals are layer type aluminosilicates that figure in terrestrial biogeochemical cycles, in the buffering capacity of the oceans, and in the containment of toxic waste materials. They are also used as lubricants in petroleum extraction and as industrial catalysts for the synthesis of many organic compounds. These applications derive fundamentally from the colloidal size and permanent structural charge of clay mineral particles, which endow them with significant surface reactivity. Unraveling the surface geochemistry of hydrated clay minerals is an abiding, if difficult, topic in earth sciences research. Recent experimental and computational studies that take advantage of new methodologies and basic insights derived from the study of concentrated ionic solutions have begun to clarify the structure of electrical double layers formed on hydrated clay mineral surfaces, particularly those in the interlayer region of swelling 2:1 layer type clay minerals. One emerging trend is that the coordination of interlayer cations with water molecules and clay mineral surface oxygens is governed largely by cation size and charge, similarly to a concentrated ionic solution, but the location of structural charge within a clay layer and the existence of hydrophobic patches on its surface provide important modulations. The larger the interlayer cation, the greater the influence of clay mineral structure and hydrophobicity on the configurations of adsorbed water molecules. This picture extends readily to hydrophobic molecules adsorbed within an interlayer region, with important implications for clay–hydrocarbon interactions and the design of catalysts for organic synthesis.

The Three-Dimensional Structure of Aspergillus niger Pectin Lyase B at 1.7-Å Resolution1

The three-dimensional structure of Aspergillus niger pectin lyase B (PLB) has been determined by crystallographic techniques at a resolution of 1.7 Å. The model, with all 359 amino acids and 339 water molecules, refines to a final crystallographic R factor of 16.5%. The polypeptide backbone folds into a large right-handed cylinder, termed a parallel β helix. Loops of various sizes and conformations protrude from the central helix and probably confer function. The largest loop of 53 residues folds into a small domain consisting of three antiparallel β strands, one turn of an α helix, and one turn of a 310 helix. By comparison with the structure of Erwinia chrysanthemi pectate lyase C (PelC), the primary sequence alignment between the pectate and pectin lyase subfamilies has been corrected and the active site region for the pectin lyases deduced. The substrate-binding site in PLB is considerably less hydrophilic than the comparable PelC region and consists of an extensive network of highly conserved Trp and His residues. The PLB structure provides an atomic explanation for the lack of a catalytic requirement for Ca2+ in the pectin lyase family, in contrast to that found in the pectate lyase enzymes. Surprisingly, however, the PLB site analogous to the Ca2+ site in PelC is filled with a positive charge provided by a conserved Arg in the pectin lyases. The significance of the finding with regard to the enzymatic mechanism is discussed.

Glass transition in DNA from molecular dynamics simulations.

Molecular dynamics simulations of the oligonucleotide duplex d(CGCGCG)2 in aqueous solution are used to investigate the glass transition phenomenon. The simulations were performed at temperatures in the 20 K to 340 K range. The mean square atomic fluctuations showed that the behavior of the oligonucleotide duplex was harmonic at low temperatures. A glass transition temperature at 223 K to 234 K was inferred for the oligonucleotide duplex, which is in agreement with experimental observations. The largest number of hydrogen bounds between the polar atoms of the oligonucleotide duplex and the water molecules was obtained at the glass transition temperature. With increasing temperature we observed a decrease in the average lifetime of the hydrogen bonds to water molecules.

A structure-based catalytic mechanism for the xanthine oxidase family of molybdenum enzymes.

The crystal structure of the xanthine oxidase-related molybdenum-iron protein aldehyde oxido-reductase from the sulfate reducing anaerobic Gram-negative bacterium Desulfovibrio gigas (Mop) was analyzed in its desulfo-, sulfo-, oxidized, reduced, and alcohol-bound forms at 1.8-A resolution. In the sulfo-form the molybdenum molybdopterin cytosine dinucleotide cofactor has a dithiolene-bound fac-[Mo, = O, = S, ---(OH2)] substructure. Bound inhibitory isopropanol in the inner compartment of the substrate binding tunnel is a model for the Michaelis complex of the reaction with aldehydes (H-C = O,-R). The reaction is proposed to proceed by transfer of the molybdenum-bound water molecule as OH- after proton transfer to Glu-869 to the carbonyl carbon of the substrate in concert with hydride transfer to the sulfido group to generate [MoIV, = O, -SH, ---(O-C = O, -R)). Dissociation of the carboxylic acid product may be facilitated by transient binding of Glu-869 to the molybdenum. The metal-bound water is replenished from a chain of internal water molecules. A second alcohol binding site in the spacious outer compartment may cause the strong substrate inhibition observed. This compartment is the putative binding site of large inhibitors of xanthine oxidase.

Potentiation of proton transfer function by electrostatic interactions in photosynthetic reaction centers from Rhodobacter sphaeroides: First results from site-directed mutation of the H subunit.

The x-ray crystallographic structure of the photosynthetic reaction center (RC) has proven critical in understanding biological electron transfer processes. By contrast, understanding of intraprotein proton transfer is easily lost in the immense richness of the details. In the RC of Rhodobacter (Rb.) sphaeroides, the secondary quinone (QB) is surrounded by amino acid residues of the L subunit and some buried water molecules, with M- and H-subunit residues also close by. The effects of site-directed mutagenesis upon RC turnover and quinone function have implicated several L-subunit residues in proton delivery to QB, although some species differences exist. In wild-type Rb. sphaeroides, Glu L212 and Asp L213 represent an inner shell of residues of particular importance in proton transfer to QB. Asp L213 is crucial for delivery of the first proton, coupled to transfer of the second electron, while Glu L212, possibly together with Asp L213, is necessary for delivery of the second proton, after the second electron transfer. We report here the first study, by site-directed mutagenesis, of the role of the H subunit in QB function. Glu H173, one of a cluster of strongly interacting residues near QB, including Asp L213, was altered to Gln. In isolated mutant RCs, the kinetics of the first electron transfer, leading to formation of the semiquinone, QB-, and the proton-linked second electron transfer, leading to the formation of fully reduced quinol, were both greatly retarded, as observed previously in the Asp L213 --> Asn mutant. However, the first electron transfer equilibrium, QA-QB <==> QAQB-, was decreased, which is opposite to the effect of the Asp L213 --> Asn mutation. These major disruptions of events coupled to proton delivery to QB were largely reversed by the addition of azide (N3-). The results support a major role for electrostatic interactions between charged groups in determining the protonation state of certain entities, thereby controlling the rate of the second electron transfer. It is suggested that the essential electrostatic effect may be to "potentiate" proton transfer activity by raising the pK of functional entities that actually transfer protons in a coupled fashion with the second electron transfer. Candidates include buried water (H3O+) and Ser L223 (serine-OH2+), which is very close to the O5 carbonyl of the quinone.

Interactions of protein antigens with antibodies.

There are now several crystal structures of antibody Fab fragments complexed to their protein antigens. These include Fab complexes with lysozyme, two Fab complexes with influenza virus neuraminidase, and three Fab complexes with their anti-idiotype Fabs. The pattern of binding that emerges is similar to that found with other protein-protein interactions, with good shape complementarity between the interacting surfaces and reasonable juxtapositions of polar residues so as to permit hydrogen-bond formation. Water molecules have been observed in cavities within the interface and on the periphery, where they often form bridging hydrogen bonds between antibody and antigen. For the most part the antigen is bound in the middle of the antibody combining site with most of the six complementarity-determining residues involved in binding. For the most studied antigen, lysozyme, the epitopes for four antibodies occupy approximately 45% of the accessible surface area. Some conformational changes have been observed to accompany binding in both the antibody and the antigen, although most of the information on conformational change in the latter comes from studies of complexes with small antigens.

Hydrated myoglobin's anharmonic fluctuations are not primarily due to dihedral transitions.

To characterize the functionally important anharmonic motions of proteins, simulations of carboxymyoglobin (MbCO) dynamics have been performed during which dihedral transitions were prohibited. Comparison of torsionally restrained and unrestrained protein dynamics simulated at three levels of hydration and at temperatures ranging from 100 to 400 K suggests that hydration "catalyzes" protein mobility by facilitating collective anharmonic motions that do not require dihedral transitions. When dihedral transitions were prohibited, dehydrated MbCO, to a good approximation, exhibited only harmonic fluctuations, whereas hydrated MbCO exhibited both harmonic and anharmonic motions. The fluctuation of helix centers of mass also remained highly anharmonic in the torsionally restrained hydrated system. Atomic mean-square fluctuation at 300 K was reduced upon prohibition of dihedral transitions by only 28% and 10% for MbCO hydrated by 350 and 3830 water molecules, respectively.

Features of MotA proton channel structure revealed by tryptophan-scanning mutagenesis.

The MotA protein of Escherichia coli is a component of the flagellar motors that functions in transmembrane proton conduction. Here, we report several features of MotA structure revealed by use of a mutagenesis-based approach. Single tryptophan residues were introduced at many positions within the four hydrophobic segments of MotA, and the effects on function were measured. Function was disrupted according to a periodic pattern that implies that the membrane-spanning segments are alpha-helices and that identifies the lipid-facing parts of each helix. The results support a hypothesis for MotA structure and mechanism in which water molecules form most of the proton-conducting pathway. The success of this approach in studying MotA suggests that it could be useful in structure-function studies of other integral membrane proteins.

Sticky ions in biological systems.

Aqueous gel sieving chromatography on Sephadex G-10 of the Group IA cations (Li+, Na+, K+, Rb+, Cs+) plus NH4+ as the Cl- salts, in combination with previous results for the halide anions (F-, Cl-, Br-, I-) as the Na+ salts [Washabaugh, M.W. & Collins, K.D. (1986) J. Biol. Chem. 261, 12477-12485], leads to the following conclusions. (i) The small monovalent ions (Li+, Na+, F-) flow through the gel with water molecules attached, whereas the large monovalent ions (K+, Rb+, Cs+, Cl-, Br-, I-) adsorb to the nonpolar surface of the gel, a process requiring partial dehydration of the ion and implying that these ions bind the immediately adjacent water molecules weakly. (ii) The transition from strong to weak hydration occurs at a radius of about 1.78 A for the monovalent anions, compared with a radius of about 1.06 A for the monovalent cations (using ionic radii), indicating that the anions are more strongly hydrated than the cations for a given charge density. (iii) The anions show larger deviations from ideal behavior (an elution position corresponding to the anhydrous molecular weight) than do the cations and dominate the chromatographic behavior of the neutral salts. These results are interpreted to mean that weakly hydrated ions (chaotropes) are "pushed" onto weakly hydrated surfaces by strong water-water interactions and that the transition from strong ionic hydration to weak ionic hydration occurs where the strength of ion-water interactions approximately equals the strength of water-water interactions in bulk solution.