Is water structure around hydrophobic groups clathrate-like?

The term "clathrate structure" is quantified for solvation of nonpolar groups by enumerating hydrogen-bonded ring sizes both in the solvation shell and through the shell-bulk interface and comparing it to a bulk control using the ST4 water model. For clathrate-like structure to be evident, the distributions along the hydrophobic surface are expected to be dominated by pentagons, with significant depletion of hexagons and larger polygons. While the distribution in this region is indeed distinguished by a large number of pentagons, there are significant contributions from hexagons and larger rings as well. Calculated polygon distributions through the shell-bulk interface indicate that when water structure is highly cooperative along the hydrophobic surface, hydrogen-bonded pathways leading back into bulk are then reduced. These results are qualitatively consistent with the observation that hydrophobicity is proportional to the nonpolar solute surface area.

Direct evidence for modified solvent structure within the hydration shell of a hydrophobic amino acid.

Neutron scattering experiments are used to determine scattering profiles for aqueous solutions of hydrophobic and hydrophilic amino acid analogs. Solutions of hydrophobic solutes show a shift in the main diffraction peak to smaller angle as compared with pure water, whereas solutions of hydrophilic solutes do not. The same difference for solutions of hydrophobic and hydrophilic side chains is also predicted by molecular dynamics simulations. The neutron scattering curves of aqueous solutions of hydrophobic amino acids at room temperature are qualitatively similar to differences between the liquid molecular structure functions measured for ambient and supercooled water. The nonpolar solute-induced expansion of water structure reported here is also complementary to recent neutron experiments where compression of aqueous solvent structure has been observed at high salt concentration.

Raman spectral evidence for hydration forces between collagen triple helices

Hydration forces are thought to result from the energetic cost of water rearrangement near macromolecular surfaces. Raman spectra, collected on the same collagen samples on which these forces were measured, reveal a continuous change in water hydrogen-bonding structure as a function of separation between collagen triple helices. The varying spectral parameters track the force-distance curve. The energetic cost of water “restructuring,” estimated from the spectra, is consistent with the measured energy of intermolecular interaction. These correlations support the idea that the change in water structure underlies the exponentially varying forces seen in this system at least over the 13–18-Å range of interaxial separations.

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.

Computer determination of peptide conformations in water: different roads to structure.

Fragments of proteins (short peptides) that "fold" suggest a mechanism of how complete conformational search in protein folding is avoided. We used a computational method to determine structures of two foldable peptides in explicit water: RVEW and CSVTC. The optimization starts from random structures and no experimental constraints are used. In agreement with NMR data, the simulations find a hydrophobic pair (Val/Trp) in REVW. The structure of CSVTC is induced by a surface water that bridges two amide hydrogens, a drive to structure hypothesized by Ben-Naim [Ben-Naim, A. (1990) J. Chem. Phys. 93, 8196-8210] that is largely ignored in studies of folding. Tendency to structure in short peptide chains suggests a mechanism for the formation of short-range nucleation sites in protein folding.

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.

Crystal structure of Bacillus subtilis YabJ, a purine regulatory protein and member of the highly conserved YjgF family

The yabJ gene in Bacillus subtilis is required for adenine-mediated repression of purine biosynthetic genes in vivo and codes for an acid-soluble, 14-kDa protein. The molecular mechanism of YabJ is unknown. YabJ is a member of a large, widely distributed family of proteins of unknown biochemical function. The 1.7-Å crystal structure of YabJ reveals a trimeric organization with extensive buried hydrophobic surface and an internal water-filled cavity. The most important finding in the structure is a deep, narrow cleft between subunits lined with nine side chains that are invariant among the 25 most similar homologs. This conserved site is proposed to be a binding or catalytic site for a ligand or substrate that is common to YabJ and other members of the YER057c/YjgF/UK114 family of proteins.

The crystal structure of the Rev binding element of HIV-1 reveals novel base pairing and conformational variability

The crystal and molecular structure of an RNA duplex corresponding to the high affinity Rev protein binding element (RBE) has been determined at 2.1-Å resolution. Four unique duplexes are present in the crystal, comprising two structural variants. In each duplex, the RNA double helix consists of an annealed 12-mer and 14-mer that form an asymmetric internal loop consisting of G-G and G-A noncanonical base pairs and a flipped-out uridine. The 12-mer strand has an A-form conformation, whereas the 14-mer strand is distorted to accommodate the bulges and noncanonical base pairing. In contrast to the NMR model of the unbound RBE, an asymmetric G-G pair with N2-N7 and N1-O6 hydrogen bonding, is formed in each helix. The G-A base pairing agrees with the NMR structure in one structural variant, but forms a novel water-mediated pair in the other. A backbone flip and reorientation of the G-G base pair is required to assume the RBE conformation present in the NMR model of the complex between the RBE and the Rev peptide.

Crystal structure of an anti-carbohydrate antibody directed against Vibrio cholerae O1 in complex with antigen: Molecular basis for serotype specificity

The crystal structure of the murine Fab S-20-4 from a protective anti-cholera Ab specific for the lipopolysaccharide Ag of the Ogawa serotype has been determined in its unliganded form and in complex with synthetic fragments of the Ogawa O-specific polysaccharide (O-SP). The upstream terminal O-SP monosaccharide is shown to be the primary antigenic determinant. Additional perosamine residues protrude outwards from the Ab surface and contribute only marginally to the binding affinity and specificity. A complementary water-excluding hydrophobic interface and five Ab–Ag hydrogen bonds are crucial for carbohydrate recognition. The structure reported here explains the serotype specificity of anti-Ogawa Abs and provides a rational basis toward the development of a synthetic carbohydrate-based anti-cholera vaccine.

Molecular phylogenetic analysis of evolutionary trends in stonefly wing structure and locomotor behavior

Insects in the order Plecoptera (stoneflies) use a form of two-dimensional aerodynamic locomotion called surface skimming to move across water surfaces. Because their weight is supported by water, skimmers can achieve effective aerodynamic locomotion even with small wings and weak flight muscles. These mechanical features stimulated the hypothesis that surface skimming may have been an intermediate stage in the evolution of insect flight, which has perhaps been retained in certain modern stoneflies. Here we present a phylogeny of Plecoptera based on nucleotide sequence data from the small subunit rRNA (18S) gene. By mapping locomotor behavior and wing structural data onto the phylogeny, we distinguish between the competing hypotheses that skimming is a retained ancestral trait or, alternatively, a relatively recent loss of flight. Our results show that basal stoneflies are surface skimmers, and that various forms of surface skimming are distributed widely across the plecopteran phylogeny. Stonefly wings show evolutionary trends in the number of cross veins and the thickness of the cuticle of the longitudinal veins that are consistent with elaboration and diversification of flight-related traits. These data support the hypothesis that the first stoneflies were surface skimmers, and that wing structures important for aerial flight have become elaborated and more diverse during the radiation of modern stoneflies.

Photochemically induced nuclear spin polarization in reaction centers of photosystem II observed by 13C-solid-state NMR reveals a strongly asymmetric electronic structure of the P680.+ primary donor chlorophyll

We report 13C magic angle spinning NMR observation of photochemically induced dynamic nuclear spin polarization (photo- CIDNP) in the reaction center (RC) of photosystem II (PS2). The light-enhanced NMR signals of the natural abundance 13C provide information on the electronic structure of the primary electron donor P680 (chlorophyll a molecules absorbing around 680 nm) and on the pz spin density pattern in its oxidized form, P680⨥. Most centerband signals can be attributed to a single chlorophyll a (Chl a) cofactor that has little interaction with other pigments. The chemical shift anisotropy of the most intense signals is characteristic for aromatic carbon atoms. The data reveal a pronounced asymmetry of the electronic spin density distribution within the P680⨥. PS2 shows only a single broad and intense emissive signal, which is assigned to both the C-10 and C-15 methine carbon atoms. The spin density appears shifted toward ring III. This shift is remarkable, because, for monomeric Chl a radical cations in solution, the region of highest spin density is around ring II. It leads to a first hypothesis as to how the planet can provide itself with the chemical potential to split water and generate an oxygen atmosphere using the Chl a macroaromatic cycle. A local electrostatic field close to ring III can polarize the electronic charge and associated spin density and increase the redox potential of P680 by stabilizing the highest occupied molecular orbital, without a major change of color. This field could be produced, e.g., by protonation of the keto group of ring V. Finally, the radical cation electronic structure in PS2 is different from that in the bacterial RC, which shows at least four emissive centerbands, indicating a symmetric spin density distribution over the entire bacteriochlorophyll macrocycle.

Covalent intermediate trapped in 2-keto-3-deoxy-6- phosphogluconate (KDPG) aldolase structure at 1.95-Å resolution

2-Keto-3-deoxy-6-phosphogluconate (KDPG) aldolase catalyzes the reversible cleavage of KDPG to pyruvate and glyceraldehyde-3-phosphate. The enzyme is a class I aldolase whose reaction mechanism involves formation of Schiff base intermediates between Lys-133 and a keto substrate. A covalent adduct was trapped by flash freezing KDPG aldolase crystals soaked with 10 mM pyruvate in acidic conditions at pH 4.6. Structure determination to 1.95-Å resolution showed that pyruvate had undergone nucleophilic attack with Lys-133, forming a protonated carbinolamine intermediate, a functional Schiff base precursor, which was stabilized by hydrogen bonding with active site residues. Carbinolamine interaction with Glu-45 indicates general base catalysis of several rate steps. Stereospecific addition is ensured by aromatic interaction of Phe-135 with the pyruvate methyl group. In the native structure, Lys-133 donates all of its hydrogen bonds, indicating the presence of an ɛ-ammonium salt group. Nucleophilic activation is postulated to occur by proton transfer in the monoprotonated zwitterionic pair (Glu-45/Lys-133). Formation of the zwitterionic pair requires prior side chain rearrangement by protonated Lys-133 to displace a water molecule, hydrogen bonded to the zwitterionic residues.

The arginine finger of RasGAP helps Gln-61 align the nucleophilic water in GAP-stimulated hydrolysis of GTP

The Ras family of GTPases is a collection of molecular switches that link receptors on the plasma membrane to signaling pathways that regulate cell proliferation and differentiation. The accessory GTPase-activating proteins (GAPs) negatively regulate the cell signaling by increasing the slow intrinsic GTP to GDP hydrolysis rate of Ras. Mutants of Ras are found in 25–30% of human tumors. The most dramatic property of these mutants is their insensitivity to the negative regulatory action of GAPs. All known oncogenic mutants of Ras map to a small subset of amino acids. Gln-61 is particularly important because virtually all mutations of this residue eliminate sensitivity to GAPs. Despite its obvious importance for carcinogenesis, the role of Gln-61 in the GAP-stimulated GTPase activity of Ras has remained a mystery. Our molecular dynamics simulations of the p21ras–p120GAP–GTP complex suggest that the local structure around the catalytic region can be different from that revealed by the x-ray crystal structure. We find that the carbonyl oxygen on the backbone of the arginine finger supplied in trans by p120GAP (Arg-789) interacts with a water molecule in the active site that is forming a bridge between the NH2 group of the Gln-61 and the γ-phosphate of GTP. Thus, Arg-789 may play a dual role in generating the nucleophile as well as stabilizing the transition state for P—O bond cleavage.

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.