Dataset Information

Dynamics of water molecules in the bacteriorhodopsin trimer in explicit lipid/water environment.

ABSTRACT: Protonated networks of internal water molecules appear to be involved in proton transfer in various integral membrane proteins. High-resolution x-ray studies of protein crystals at low temperature deliver mean positions of most internal waters, but only limited information about fluctuations within such H-bonded networks formed by water and residues. The question arises as to how water molecules behave inside and on the surface of a fluctuating membrane protein under more physiological conditions. Therefore, as an example, long-time molecular dynamics simulations of bacteriorhodopsin were performed with explicit membrane/water environment. Based on a recent x-ray model the bacteriorhodopsin trimer was inserted in a fully solvated 16 x 16 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC)-bilayer patch, resulting in a system of approximately 84,000 atoms. Unrestrained molecular dynamics calculations of 5 ns were performed using the GROMACS package and force field. Mean water densities were computed to describe the anisotropic distribution of internal water molecules. In the whole protein two larger areas of higher water density are identified. They are located between the central proton binding site, the Schiff base, and the extracellular proton release site. Separated by Arg-82 these water clusters could provide a proton release pathway in a Grotthus-like mechanism as indicated by a continuum absorbance change observed during the photocycle by time-resolved Fourier transform infrared spectroscopy. Residues are identified which are H-bonded to the water clusters and are therefore facilitating proton conduction. Their influence on proton transfer via the H-bonded network as indicated by the continuum absorbance change is predicted. This may explain why several site-directed mutations alter the proton release kinetics without a direct involvement in proton transfer.

SUBMITTER: Kandt C

PROVIDER: S-EPMC1303921 | biostudies-literature | 2004 Feb

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Dynamics of water molecules in the bacteriorhodopsin trimer in explicit lipid/water environment.

Kandt Christian C Schlitter Jürgen J Gerwert Klaus K

Biophysical journal 20040201 2

Protonated networks of internal water molecules appear to be involved in proton transfer in various integral membrane proteins. High-resolution x-ray studies of protein crystals at low temperature deliver mean positions of most internal waters, but only limited information about fluctuations within such H-bonded networks formed by water and residues. The question arises as to how water molecules behave inside and on the surface of a fluctuating membrane protein under more physiological condition ...[more]

PMID: 14747309

Similar Datasets

Project description:The light-driven proton pump bacteriorhodopsin (bR) is a transmembrane protein that uses large conformational changes for proton transfer from the cytoplasmic to the extracellular regions. Crystal structures, due to their solvent conditions, do not resolve the effect of lipid molecules on these protein conformational changes. To begin to understand the molecular details behind such large conformational changes, we simulated two conformations of wild-type bacteriorhodopsin, one of the dark-adapted state and the second of an intermediate (M(O)) state, each within an explicit dimyristoyl-phosphatidylcholine (DMPC) lipid bilayer. The simulations included all-hydrogen and all-atom representations of protein, lipid, and water and were performed for 20 ns. We investigate the equilibrium properties and the dynamic motions of the two conformations in the lipid setting. We note that the conformational state of the M(O) intermediate bR remains markedly different from the dark-adapted bR state in that the M(O) intermediate shows rearrangement of the cytoplasmic portions of helices C, F, and G, and nearby loops. This difference in the states remained throughout the simulations, and the results are stable on the molecular dynamics timescale and provide an illustration of the changes in both lipid and water that help to stabilize a particular state. Our analysis focuses on how the environment adjusts to these two states and on how the dynamics of the helices, loops, and water molecules can be related to the pump mechanism of bacteriorhodopsin. For example, water generally behaves in the same manner on the extracellular sides of both simulations but is decreased in the cytoplasmic region of the M(O) intermediate. We suspect that the different water behavior is closely related to the fluctuations of microcavities volume in the protein interior, which is strongly coupled to the collective motion of the protein. Our simulation result suggests that experimental observation can be useful to verify a decreased number of waters in the cytoplasmic regions of the late-intermediate stages by measuring the rate of water exchange with the interior of the protein.

Project description:Lipid membranes are more than just barriers between cell compartments; they provide molecular environments with a finely tuned balance between hydrophilic and hydrophobic interactions that enable proteins to dynamically fold and self-assemble to regulate biological function. Characterizing dynamics at the lipid-water interface is essential to understanding molecular complexities from the thermodynamics of liquid-liquid phase separation down to picosecond-scale reorganization of interfacial hydrogen-bond networks.Ultrafast vibrational spectroscopy, including two-dimensional infrared (2D IR) and vibrational sum-frequency generation (VSFG) spectroscopies, is a powerful tool to examine picosecond interfacial dynamics. Two-dimensional IR spectroscopy provides a bond-centered view of dynamics with subpicosecond time resolutions, as vibrational frequencies are highly sensitive to the local environment. Recently, 2D IR spectroscopy has been applied to carbonyl and phosphate vibrations intrinsically located at the lipid-water interface. Interface-specific VSFG spectroscopy probes the water vibrational modes directly, accessing H-bond strength and water organization at lipid headgroup positions. Signals in VSFG arise from the interfacial dipole contributions, directly probing headgroup ordering and water orientation to provide a structural view of the interface.In this Account we discuss novel applications of ultrafast spectroscopy to lipid membranes, a field that has experienced significant growth over the past decade. In particular, ultrafast experiments now offer a molecular perspective on increasingly complex membranes. The powerful combination of ultrafast, interface-selective spectroscopy and simulations opens up new routes to understanding multicomponent membranes and their function. This Account highlights key prevailing views that have emerged from recent experiments: (1) Water dynamics at the lipid-water interface are slow compared to those of bulk water as a result of disrupted H-bond networks near the headgroups. (2) Peptides, ions, osmolytes, and cosolvents perturb interfacial dynamics, indicating that dynamics at the interface are affected by bulk solvent dynamics and vice versa. (3) The interfacial environment is generally dictated by the headgroup structure and orientation, but hydrophobic interactions within the acyl chains also modulate interfacial dynamics. Ultrafast spectroscopy has been essential to characterizing the biophysical chemistry of the lipid-water interface; however, challenges remain in interpreting congested spectra as well as designing appropriate model systems to capture the complexity of a membrane environment.

Project description:The mechanism of acid-catalyzed keto-enol tautomerization of ?-cyclopentanedione (CPD) in solution is studied computationally. Reaction profiles are first calculated for a limited solvation environment using ab initio and density functional methods. Barrier heights for systems including up to one hydration shell of explicit water molecules depend strongly on the number of waters involved in proton transfer and to a lesser but significant extent on the number of waters forming hydrogen bonds with waters in the proton-transfer chain (each such water reduces the barrier by 4.4 kcal/mol on average). Barriers of 8-13 kcal/mol were obtained when a full or nearly full hydration shell was present, consistent with calculations for nonacid-catalyzed keto-enol tautomerization of related molecules. The presence of HOCl reduced the barrier by 4.5 kcal/mol in relation to the gas phase, consistent with the well-known principle that keto-enol tautomerization can be acid- or base-catalyzed. The reaction was also modeled beginning with snapshots of reactant conformations taken from a 300 K molecular dynamics simulation of CPD, HOCl, and 324 explicit waters. Reaction profiles were calculated at a QM/MM level with waters in the first hydration shell either fixed or energy-minimized at each step along the reaction coordinate. A substantial variation in barrier height was observed in both cases, depending primarily on electrostatic interactions (hydrogen bonding) with first-hydration-shell waters and, to a lesser extent, on electrostatic interactions with more distant waters and geometric distortion effects. For the lowest barriers, the extent of barrier reduction by waters involved in proton transfer is consistent with the limited solvation results, but further barrier reduction due to hydrogen bonding to waters involved in proton transfer is not observed. It is postulated that this is because highly flexible structures such as extensive hydrogen bonding networks optimal for reaction are entropically disfavored and so may not contribute significantly to the observed reaction rate.

Dataset Information

Dynamics of water molecules in the bacteriorhodopsin trimer in explicit lipid/water environment.

Publications

Dynamics of water molecules in the bacteriorhodopsin trimer in explicit lipid/water environment.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets