The role of water structure in conformational changes of nucleic acids in ambient and high-pressure conditions.
(25/15056)
This review describes and summarizes data on the structure and properties of water under normal conditions, at high salt concentration and under high pressure. We correlate the observed conformational changes in nucleic acids with changes in water structure and activity, and suggest a mechanism of conformational transitions of nucleic acids which accounts for changes in the water structure. From the biophysical, biochemical and crystallographic data we conclude that the Z-DNA form can be induced only at low water activity produced by high salt concentrations or high pressure, and accompanied by the stabilizing conjugative effect of the cytidine O4' electrons of the CG base pairs. (+info)
Hydrophobic hydration of amphipathic peptides.
(26/15056)
Biomolecular surfaces and interfaces are commonly found with apolar character. The hydrophobic effect thus plays a crucial role in processes involving association with biomolecular surfaces in the cellular environment. By computer simulation, we compared the hydrogen bonding structures and energetics of the proximal hydration shells of the monomer and dimer from a recent study of an extrinsic membrane peptide, melittin. The two peptides were studied in their amphipathic alpha-helical forms, which possess extended hydrophobic surfaces characterized by different topography. The topography of the peptide-water interface was found to be critical in determining the enthalpic nature of hydrophobic hydration. This topographical dependence has far-reaching implications in the regulation of bioactivities in the presence of amphipathicity. This result also engenders reconsideration of the validity of using free energy parameters that depend solely on the chemical nature of constituent moieties in characterizing hydrophobic hydration of proteins and biomolecules in general. (+info)
An alamethicin channel in a lipid bilayer: molecular dynamics simulations.
(27/15056)
We present the results of 2-ns molecular dynamics (MD) simulations of a hexameric bundle of Alm helices in a 1-palmitoyl-2-oleoylphosphatidylcholine bilayer. These simulations explore the dynamic properties of a model of a helix bundle channel in a complete phospholipid bilayer in an aqueous environment. We explore the stability and conformational dynamics of the bundle in a phospholipid bilayer. We also investigate the effect on bundle stability of the ionization state of the ring of Glu18 side chains. If all of the Glu18 side chains are ionised, the bundle is unstable; if none of the Glu18 side chains are ionized, the bundle is stable. pKA calculations suggest that either zero or one ionized Glu18 is present at neutral pH, correlating with the stable form of the helix bundle. The structural and dynamic properties of water in this model channel were examined. As in earlier in vacuo simulations (Breed et al., 1996 .Biophys. J. 70:1643-1661), the dipole moments of water molecules within the pore were aligned antiparallel to the helix dipoles. This contributes to the stability of the helix bundle. (+info)
Defining the transmembrane helix of M2 protein from influenza A by molecular dynamics simulations in a lipid bilayer.
(28/15056)
Integral membrane proteins containing at least one transmembrane (TM) alpha-helix are believed to account for between 20% and 30% of most genomes. There are several algorithms that accurately predict the number and position of TM helices within a membrane protein sequence. However, these methods tend to disagree over the beginning and end residues of TM helices, posing problems for subsequent modeling and simulation studies. Molecular dynamics (MD) simulations in an explicit lipid and water environment are used to help define the TM helix of the M2 protein from influenza A virus. Based on a comparison of the results of five different secondary structure prediction algorithms, three different helix lengths (an 18mer, a 26mer, and a 34mer) were simulated. Each simulation system contained 127 POPC molecules plus approximately 3500-4700 waters, giving a total of approximately 18,000-21,000 atoms. Two simulations, each of 2 ns duration, were run for the 18mer and 26mer, and five separate simulations were run for the 34mer, using different starting models generated by restrained in vacuo MD simulations. The total simulation time amounted to 11 ns. Analysis of the time-dependent secondary structure of the TM segments was used to define the regions that adopted a stable alpha-helical conformation throughout the simulation. This analysis indicates a core TM region of approximately 20 residues (from residue 22 to residue 43) that remained in an alpha-helical conformation. Analysis of atomic density profiles suggested that the 18mer helix revealed a local perturbation of the lipid bilayer. Polar side chains on either side of this region form relatively long-lived H-bonds to lipid headgroups and water molecules. (+info)
Simulation analysis of the retinal conformational equilibrium in dark-adapted bacteriorhodopsin.
(29/15056)
In dark-adapted bacteriorhodopsin (bR) the retinal moiety populates two conformers: all-trans and (13,15)cis. Here we examine factors influencing the thermodynamic equilibrium and conformational transition between the two forms, using molecular mechanics and dynamics calculations. Adiabatic potential energy mapping indicates that whereas the twofold intrinsic torsional potentials of the C13==C14 and C15==N16 double bonds favor a sequential torsional pathway, the protein environment favors a concerted, bicycle-pedal mechanism. Which of these two pathways will actually occur in bR depends on the as yet unknown relative weight of the intrinsic and environmental effects. The free energy difference between the conformers was computed for wild-type and modified bR, using molecular dynamics simulation. In the wild-type protein the free energy of the (13,15)cis retinal form is calculated to be 1.1 kcal/mol lower than the all-trans retinal form, a value within approximately kBT of experiment. In contrast, in isolated retinal the free energy of the all-trans state is calculated to be 2.1 kcal/mol lower than (13,15)cis. The free energy differences are similar to the adiabatic potential energy differences in the various systems examined, consistent with an essentially enthalpic origin. The stabilization of the (13,15)cis form in bR relative to the isolated retinal molecule is found to originate from improved protein-protein interactions. Removing internal water molecules near the Schiff base strongly stabilizes the (13,15)cis form, whereas a double mutation that removes negative charges in the retinal pocket (Asp85 to Ala; Asp212 to Ala) has the opposite effect. (+info)
Simulation study of a gramicidin/lipid bilayer system in excess water and lipid. I. Structure of the molecular complex.
(30/15056)
This paper reports on a simulation of a gramicidin channel inserted into a fluid phase DMPC bilayer with 100 lipid molecules. Two lipid molecules per leaflet were removed to insert the gramicidin, so the resulting preparation had 96 lipid molecules and 3209 water molecules. Constant surface tension boundary conditions were employed. Like previous simulations with a lower lipid/gramicidin ratio (Woolf, T. B., and B. Roux. 1996. Proteins: Struct., Funct., Genet. 24:92-114), it is found that tryptophan-water hydrogen bonds are more common than tryptophan-phospholipid hydrogen bonds. However, one of the tryptophan NH groups entered into an unusually long-lived hydrogen bonding pattern with two glycerol oxygens of one of the phospholipid molecules. Comparisons were made between the behavior of the lipids adjacent to the channel with those farther away. It was found that hydrocarbon chains of lipids adjacent to the channel had higher-order parameters than those farther away. The thickness of the lipid bilayer immediately adjacent to the channel was greater than it was farther away. In general, the lipids adjacent to the membrane had similar orientations to those seen by Woolf and Roux, while those farther away had similar orientations to those pertaining before the insertion of the gramicidin. A corollary to this observation is that the thickness of the hydrocarbon region adjacent to the gramicidin was much thicker than what other studies have identified as the "hydrophobic length" of the gramicidin channel. (+info)
Simulation study of a gramicidin/lipid bilayer system in excess water and lipid. II. Rates and mechanisms of water transport.
(31/15056)
A gramicidin channel in a fluid phase DMPC bilayer with excess lipid and water has been simulated. By use of the formal correspondence between diffusion and random walk, a permeability for water through the channel was calculated, and was found to agree closely with the experimental results of Rosenberg and Finkelstein (Rosenberg, P.A., and A. Finkelstein. 1978. J. Gen. Physiol. 72:327-340; 341-350) for permeation of water through gramicidin in a phospholipid membrane. By using fluctuation analysis, components of resistance to permeation were computed for movement through the channel interior, for the transition step at the channel mouth where the water molecule solvation environment changes, and for the process of diffusion up to the channel mouth. The majority of the resistance to permeation appears to occur in the transition step at the channel mouth. A significant amount is also due to structurally based free energy barriers within the channel. Only small amounts are due to local friction within the channel or to diffusive resistance for approaching the channel mouth. (+info)
Histone-hydrocarbon interaction. Partition of histones in aqueous two-phase systems containing poly(ethylene glycol)-bound hydrocarbons.
(32/15056)
The hydrophobic properties of histones have been examined with help of the two-phase partition technique using dextran-poly(ethylene glycol)-water systems. We have found that different fatty acid esters of poly(ethylene glycol) interact with total histones in a manner similar to proteins of the type beta-lactoglobulin and serum albumins. Thus the maximum interaction occurs when the fatty acid contains 16-18 carbon atoms. With less than eight carbon atoms in the polymer-bound fatty acid, no histone-hydrocarbon interaction is observed. The interaction of the five individual histone fractions with palmitate depends on the type of salt used and on its concentration. We suggest that the histones can be divided into three groups with decreasing hydrophobic properties: H3, H2a greater than H4, H2b greater than H1. (+info)