Influence of Water Content and Temperature on Molecular Mobility and Intracellular Glasses in Seeds and Pollen1

Although the occurrence of intracellular glasses in seeds and pollen has been established, physical properties such as rotational correlation times and viscosity have not been studied extensively. Using electron paramagnetic resonance spectroscopy, we examined changes in the molecular mobility of the hydrophilic nitroxide spin probe 3-carboxy-proxyl during melting of intracellular glasses in axes of pea (Pisum sativum L.) seeds and cattail (Typha latifolia L.) pollen. The rotational correlation time of the spin probe in intracellular glasses of both organisms was approximately 10−3 s. Using the distance between the outer extrema of the electron paramagnetic resonance spectrum (2Azz) as a measure of molecular mobility, we found a sharp increase in mobility at a definite temperature during heating. This temperature increased with decreasing water content of the samples. Differential scanning calorimetry data on these samples indicated that this sharp increase corresponded to melting of the glassy matrix. Molecular mobility was found to be inversely correlated with storage stability. With decreasing water content, the molecular mobility reached a minimum, and increased again at very low water content. Minimum mobility and maximum storage stability occurred at a similar water content. This correlation suggests that storage stability might be at least partially controlled by molecular mobility. At low temperatures, when storage longevity cannot be determined on a realistic time scale, 2Azz measurements can provide an estimate of the optimum storage conditions.

Molecular and Physiological Responses to Water Deficit in Drought-Tolerant and Drought-Sensitive Lines of Sunflower1 : Accumulation of Dehydrin Transcripts Correlates with Tolerance

To investigate correlations between phenotypic adaptation to water limitation and drought-induced gene expression, we have studied a model system consisting of a drought-tolerant line (R1) and a drought-sensitive line (S1) of sunflowers (Helianthus annuus L.) subjected to progressive drought. R1 tolerance is characterized by the maintenance of shoot cellular turgor. Drought-induced genes (HaElip1, HaDhn1, and HaDhn2) were previously identified in the tolerant line. The accumulation of the corresponding transcripts was compared as a function of soil and leaf water status in R1 and S1 plants during progressive drought. In leaves of R1 plants the accumulation of HaDhn1 and HaDhn2 transcripts, but not HaElip1 transcripts, was correlated with the drought-adaptive response. Drought-induced abscisic acid (ABA) concentration was not associated with the varietal difference in drought tolerance. Stomata of both lines displayed similar sensitivity to ABA. ABA-induced accumulation of HaDhn2 transcripts was higher in the tolerant than in the sensitive genotype. HaDhn1 transcripts were similarly accumulated in the tolerant and in the sensitive plants in response to ABA, suggesting that additional factors involved in drought regulation of HaDhn1 expression might exist in tolerant plants.

Role of a critical water in scytalone dehydratase-catalyzed reaction

Scytalone dehydratase (EC 4.2.1.94) catalyzes the dehydration of two important intermediates in the biosynthesis of melanin, and it functions without metal ions or any cofactors. Using molecular orbital theory, we have examined the role of a critical water molecule in the mechanism of scytalone dehydratase. The water, together with an internal hydrogen bonding, contributes significantly to the stabilization of the transition state (or the enolate intermediate). The role of two active site tyrosines (Tyr-50 and Tyr-30) is (i) to hold the critical water in place so that it may stabilize the transition state without much structural rearrangement during the catalytic reaction, and (ii) to polarize the water, making it a better general acid. The stereochemistry of the scytalone dehydratase-catalyzed dehydration is also discussed.

A molecular level picture of the stabilization of A-DNA in mixed ethanol–water solutions

Advances in computer power, methodology, and empirical force fields now allow routine “stable” nanosecond-length molecular dynamics simulations of DNA in water. The accurate representation of environmental influences on structure remains a major, unresolved issue. In contrast to simulations of A-DNA in water (where an A-DNA to B-DNA transition is observed) and in pure ethanol (where disruption of the structure is observed), A-DNA in ≈85% ethanol solution remains in a canonical A-DNA geometry as expected. The stabilization of A-DNA by ethanol is likely due to disruption of the spine of hydration in the minor groove and the presence of ion-mediated interhelical bonds and extensive hydration across the major groove.

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.

Chiral molecular self-assembly of phospholipid tubules: A circular dichroism study

We report on spectroscopic studies of the chiral structure in phospholipid tubules formed in mixtures of alcohol and water. Synthetic phospholipids containing diacetylenic moieties in the acyl chains self-assemble into hollow, cylindrical tubules in appropriate conditions. Circular dichroism provides a direct measure of chirality of the molecular structure. We find that the CD spectra of tubules formed in mixtures of alcohol and water depends strongly on the alcohol used and the lipid concentration. The relative spectral intensity of different circular dichroism bands correlates with the number of bilayers observed using microscopy. The results provide experimental evidence that tubule formation is based on chiral packing of the lipid molecules and that interbilayer interactions are important to the tubule structure.

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.

The importance of reactant positioning in enzyme catalysis: A hybrid quantum mechanics/molecular mechanics study of a haloalkane dehalogenase

Hybrid quantum mechanics/molecular mechanics calculations using Austin Model 1 system-specific parameters were performed to study the SN2 displacement reaction of chloride from 1,2-dichloroethane (DCE) by nucleophilic attack of the carboxylate of acetate in the gas phase and by Asp-124 in the active site of haloalkane dehalogenase from Xanthobacter autotrophicus GJ10. The activation barrier for nucleophilic attack of acetate on DCE depends greatly on the reactants having a geometry resembling that in the enzyme or an optimized gas-phase structure. It was found in the gas-phase calculations that the activation barrier is 9 kcal/mol lower when dihedral constraints are used to restrict the carboxylate nucleophile geometry to that in the enzyme relative to the geometries for the reactants without dihedral constraints. The calculated quantum mechanics/molecular mechanics activation barriers for the enzymatic reaction are 16.2 and 19.4 kcal/mol when the geometry of the reactants is in a near attack conformer from molecular dynamics and in a conformer similar to the crystal structure (DCE is gauche), respectively. This haloalkane dehalogenase lowers the activation barrier for dehalogenation of DCE by 2–4 kcal/mol relative to the single point energies of the enzyme's quantum mechanics atoms in the gas phase. SN2 displacements of this sort in water are infinitely slower than in the gas phase. The modest lowering of the activation barrier by the enzyme relative to the reaction in the gas phase is consistent with mutation experiments.

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.

Quantification of the hydrophobic interaction by simulations of the aggregation of small hydrophobic solutes in water

The hydrophobic interaction, the tendency for nonpolar molecules to aggregate in solution, is a major driving force in biology. In a direct approach to the physical basis of the hydrophobic effect, nanosecond molecular dynamics simulations were performed on increasing numbers of hydrocarbon solute molecules in water-filled boxes of different sizes. The intermittent formation of solute clusters gives a free energy that is proportional to the loss in exposed molecular surface area with a constant of proportionality of 45 ± 6 cal/mol⋅Å2. The molecular surface area is the envelope of the solute cluster that is impenetrable by solvent and is somewhat smaller than the more traditional solvent-accessible surface area, which is the area transcribed by the radius of a solvent molecule rolled over the surface of the cluster. When we apply a factor relating molecular surface area to solvent-accessible surface area, we obtain 24 cal/mol⋅Å2. Ours is the first direct calculation, to our knowledge, of the hydrophobic interaction from molecular dynamics simulations; the excellent qualitative and quantitative agreement with experiment proves that simple van der Waals interactions and atomic point-charge electrostatics account for the most important driving force in biology.

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.

Intrinsic compressibility and volume compression in solvated proteins by molecular dynamics simulation at high pressure.

Constant pressure and temperature molecular dynamics techniques have been employed to investigate the changes in structure and volumes of two globular proteins, superoxide dismutase and lysozyme, under pressure. Compression (the relative changes in the proteins' volumes), computed with the Voronoi technique, is closely related with the so-called protein intrinsic compressibility, estimated by sound velocity measurements. In particular, compression computed with Voronoi volumes predicts, in agreement with experimental estimates, a negative bound water contribution to the apparent protein compression. While the use of van der Waals and molecular volumes underestimates the intrinsic compressibilities of proteins, Voronoi volumes produce results closer to experimental estimates. Remarkably, for two globular proteins of very different secondary structures, we compute identical (within statistical error) protein intrinsic compressions, as predicted by recent experimental studies. Changes in the protein interatomic distances under compression are also investigated. It is found that, on average, short distances compress less than longer ones. This nonuniform contraction underlines the peculiar nature of the structural changes due to pressure in contrast with temperature effects, which instead produce spatially uniform changes in proteins. The structural effects observed in the simulations at high pressure can explain protein compressibility measurements carried out by fluorimetric and hole burning techniques. Finally, the calculation of the proteins static structure factor shows significant shifts in the peaks at short wavenumber as pressure changes. These effects might provide an alternative way to obtain information concerning compressibilities of selected protein regions.

The human AQP4 gene: definition of the locus encoding two water channel polypeptides in brain.

The aquaporin family of membrane water transport proteins are expressed in diverse tissues, and in brain the predominant water channel protein is AQP4. Here we report the isolation and characterization of the human AQP4 cDNAs and genomic DNA. Two cDNAs were isolated corresponding to the two initiating methionines (M1 in a 323-aa polypeptide and M23 in a 301-aa polypeptide) previously identified in rat [Jung, J.S., Bhat, R.V., Preston, G.M., Guggino, W.B. & Agre, P. (1994) Proc. Natl. Acad. Sci. USA 91, 13052-13056]. Similar to other aquaporins, the AQP4 gene is composed of four exons encoding 127, 55, 27, and 92 amino acids separated by introns of 0.8, 0.3, and 5.2 kb. Unlike other aquaporins, an alternative coding initiation sequence (designated exon 0) was located 2.7 kb upstream of exon 1. When spliced together, M1 and the subsequent 10 amino acids are encoded by exon 0; the next 11 amino acids and M23 are encoded by exon 1. Transcription initiation sites have been mapped in the proximal promoters of exons 0 and 1. RNase protection revealed distinct transcripts corresponding to M1 and M23 mRNAs, and AQP4 immunoblots of cerebellum demonstrated reactive polypeptides of 31 and 34 kDa. Using a P1 and a lambda EMBL subclone, the chromosomal site of the human AQP4 gene was mapped to chromosome 18 at the junction of q11.2 and q12.1 by fluorescence in situ hybridization. These studies may now permit molecular characterization of AQP4 during human development and in clinical disorders.

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.