Role of internal water molecules in bacteriorhodopsin

Internal water molecules are considered to play a crucial role in the functional processes of proton pump proteins. They may participate in hydrogen-bonding networks inside proteins that constitute proton pathways. In addition, they could participate in the switch reaction by mediating an essential proton transfer at the active site. Nevertheless, little has been known about the structure and function of internal water molecules in such proteins. Recent progress in infrared spectroscopy and X-ray crystallography provided new information on water molecules inside bacteriorhodopsin, the light-driven proton pump. The accumulated knowledge on bacteriorhodopsin in the last decade of the 20th century will lead to a realistic picture of internal water molecules at work in the 21st century. In this review, I describe how the role of water molecules has been studied in bacteriorhodopsin, and what should be known about the role of water molecules in the future.     

Molecular reaction mechanisms of proteins monitored by time-resolved FTIR-spectroscopy.

Time-resolved FTIR difference spectroscopy can provide a valuable insight into the molecular reaction mechanisms of proteins, especially membrane proteins. Isotopic labeling and site-directed mutagenesis allows an unequivocal assignment of IR absorption bands. Studies are presented which give insight into the proton pump mechanisms of proteins, especially bacteriorhodopsin. H-bonded network proton transfer via internal water molecules seems to be a general feature in proteins, also found in cytochrome c oxidase. Using caged GTP the intrinsic and GAP catalyzed GTPase activity of H-ras p21 is studied. Furthermore, protein folding reactions can be recorded with ns time-resolution.     

Thermodynamic stability of water molecules in the bacteriorhodopsin proton channel: a molecular dynamics free energy perturbation study.

The proton transfer activity of the light-driven proton pump, bacteriorhodopsin (bR) in the photochemical cycle might imply internal water molecules. The free energy of inserting water molecules in specific sites along the bR transmembrane channel has been calculated using molecular dynamics simulations based on a microscopic model. The existence of internal hydration is related to the free energy change on transfer of a water molecule from bulk solvent into a specific binding site. Thermodynamic integration and perturbation methods were used to calculate free energies of hydration for each hydrated model from molecular dynamics simulations of the creation of water molecules into specific protein-binding sites. A rigorous statistical mechanical formulation allowing the calculation of the free energy of transfer of water molecules from the bulk to a protein cavity is used to estimate the probabilities of occupancy in the putative bR proton channel. The channel contains a region lined primarily by nonpolar side-chains. Nevertheless, the results indicate that the transfer of four water molecules from bulk water to this apparently hydrophobic region is thermodynamically permitted. The column forms a continuous hydrogen-bonded chain over 12 A between a proton donor, Asp 96, and the retinal Schiff base acceptor. The presence of two water molecules in direct hydrogen-bonding association with the Schiff base is found to be strongly favorable thermodynamically. The implications of these results for the mechanism of proton transfer in bR are discussed.     

Water dynamics simulation as a tool for probing proton transfer pathways in a heptahelical membrane protein

The proton transfer pathway in a heptahelical membrane protein, the light-driven proton pump bacteriorhodopsin (BR), is probed by a combined approach of structural analysis of recent X-ray models and molecular dynamics (MD) simulations that provide the diffusion pathways of internal and external water molecules. Analyzing the hydrogen-bond contact frequencies of the water molecules with protein groups, the complete proton pathway through the protein is probed. Beside the well-known proton binding sites in the protein interior-the protonated Schiff base, Asp85 and Asp96, and the H(5)O(2) (+) complex stabilized by Glu204 and Glu194-the proton release and uptake pathways to the protein surfaces are described in great detail. Further residues were identified, by mutation of which the proposed pathways can be verified. In addition the diffusion pathway of water 502 from Lys216 to Asp96 is shown to cover the positions of the intruding waters 503 and 504 in the N-intermediate. The transiently established water chain in the N-state provides a proton pathway from Asp96 to the Schiff base in the M- to N-transition in a Grotthus-like mechanism, as concluded earlier from time-resolved Fourier transform infrared experiments [le Coutre et al., Proc Nat Acad Sci USA 1995;92:4962-4966].     

Conformational changes in the archaerhodopsin-3 proton pump: detection of conserved strongly hydrogen bonded water networks

Archaerhodopsin-3 (AR3) is a light-driven proton pump from Halorubrum sodomense, but little is known about its photocycle. Recent interest has focused on AR3 because of its ability to serve both as a high-performance, genetically-targetable optical silencer of neuronal activity and as a membrane voltage sensor. We examined light-activated structural changes of the protein, retinal chromophore, and internal water molecules during the photocycle of AR3. Low-temperature and rapid-scan time-resolved FTIR-difference spectroscopy revealed that conformational changes during formation of the K, M, and N photocycle intermediates are similar, although not identical, to bacteriorhodopsin (BR). Positive/negative bands in the region above 3,600 cm^− 1, which have previously been assigned to structural changes of weakly hydrogen bonded internal water molecules, were substantially different between AR3 and BR. This included the absence of positive bands recently associated with a chain of proton transporting water molecules in the cytoplasmic channel and a weakly hydrogen bonded water (W401), which is part of a hydrogen-bonded pentagonal cluster located near the retinal Schiff base. However, many of the broad IR continuum absorption changes below 3,000 cm^− 1 assigned to networks of water molecules involved in proton transport through cytoplasmic and extracellular portions in BR were very similar in AR3. This work and subsequent studies comparing BR and AR3 structural changes will help identify conserved elements in BR-like proton pumps as well as bioengineer AR3 to optimize neural silencing and voltage sensing.     

Calculation of proton transfers in Bacteriorhodopsin bR and M intermediates

Residue ionization states were calculated in nine crystal structures of bacteriorhodopsin trapped in bR, early M, and late M states by multiconformation continuum electrostatics. This combines continuum electrostatics and molecular mechanics, deriving equilibrium distributions of ionization states and polar residue and water positions. The three central cluster groups [retinal Schiff base (SB), Asp 85 and Asp 212] are ionized in bR structures while a proton has transferred from SB(+) to Asp 85(-) in late M structures matching experimental results. The proton shift in M is due to weaker SB(+)-ionized acid and more favorable SB(0)-ionized acid interactions following retinal isomerization. The proton release cluster (Glu 194 and Glu 204) binds one proton in bR, which is lost to water by pH 8 in late M. In bR the half-ionized state is stabilized by charge-dipole interactions while full ionization is disallowed by charge-charge repulsion between the closely spaced acids. In M the acids move apart, permitting full ionization. Arg 82 movement connects the proton shifts in the central and proton release clusters. Changes in total charge of the two clusters are coupled by direct long-range interactions. Separate calculations consider continuum or explicit water in internal cavities. The explicit waters and nearby polar residues can reorient to stabilize different charge distributions. Proton release to the low-pH, extracellular side of the protein occurs in these calculations where residue ionization remains at equilibrium with the medium. Thus, the key changes distinguishing the intermediates are indeed trapped in the structures.     

Interplay between local protein interactions and water bridging of a proton antenna carboxylate cluster

Proteins that bind protons at cell membrane interfaces often expose to the bulk clusters of carboxylate and histidine sidechains that capture protons transiently and, in proton transporters, deliver protons to an internal site. The protonation-coupled dynamics of bulk-exposed carboxylate clusters, also known as proton antennas, is poorly described. An essential open question is how water-mediated bridges between sidechains of the cluster respond to protonation change and facilitate transient proton storage. To address this question, here I studied the protonation-coupled dynamics at the proton-binding antenna of PsbO, a small extrinsinc subunit of the photosystem II complex, with atomistic molecular dynamics simulations and systematic graph-based analyses of dynamic protein and protein-water hydrogen-bond networks. The protonation of specific carboxylate groups is found to impact the dynamics of their local protein-water hydrogen-bond clusters. Regardless of the protonation state considered for PsbO, carboxylate pairs that can sample direct hydrogen bonding, or bridge via short hydrogen-bonded water chains, anchor to nearby basic or polar protein sidechains. As a result, carboxylic sidechains of the hypothesized antenna cluster are part of dynamic hydrogen bond networks that may rearrange rapidly when the protonation changes.     

Extended protein/water H-bond networks in photosynthetic water oxidation

Oxidation of water molecules in the photosystem II (PSII) protein complex proceeds at the manganese-calcium complex, which is buried deeply in the lumenal part of PSII. Understanding the PSII function requires knowledge of the intricate coupling between the water-oxidation chemistry and the dynamic proton management by the PSII protein matrix. Here we assess the structural basis for long-distance proton transfer in the interior of PSII and for proton management at its surface. Using the recent high-resolution crystal structure of PSII, we investigate prominent hydrogen-bonded networks of the lumenal side of PSII. This analysis leads to the identification of clusters of polar groups and hydrogen-bonded networks consisting of amino acid residues and water molecules. We suggest that long-distance proton transfer and conformational coupling is facilitated by hydrogen-bonded networks that often involve more than one protein subunit. Proton-storing Asp/Glu dyads, such as the D1-E65/D2-E312 dyad connected to a complex water-wire network, may be particularly important for coupling protonation states to the protein conformation. Clusters of carboxylic amino acids could participate in proton management at the lumenal surface of PSII. We propose that rather than having a classical hydrophobic protein interior, the lumenal side of PSII resembles a complex polyelectrolyte with evolutionary optimized hydrogen-bonding networks. This article is part of a Special Issue entitled: Photosynthesis Research for Sustainability: from Natural to Artificial.     

The role of protein-bound water molecules in microbial rhodopsins

Protein-bound internal water molecules are essential features of the structure and function of microbial rhodopsins. Besides structural stabilization, they act as proton conductors and even proton storage sites. Currently, the most understood model system exhibiting such features is bacteriorhodopsin (bR). During the last 20 years, the importance of water molecules for proton transport has been revealed through this protein. It has been shown that water molecules are as essential as amino acids for proton transport and biological function. In this review, we present an overview of the historical development of this research on bR. We furthermore summarize the recently discovered protein-bound water features associated with proton transport. Specifically, we discuss a pentameric water/amino acid arrangement close to the protonated Schiff base as central proton-binding site, a protonated water cluster as proton storage site at the proton-release site, and a transient linear water chain at the proton uptake site. We highlight how protein conformational changes reposition or reorient internal water molecules, thereby guiding proton transport. Last, we compare the water positions in bR with those in other microbial rhodopsins to elucidate how protein-bound water molecules guide the function of microbial rhodopsins. This article is part of a Special Issue entitled: Retinal Proteins — You can teach an old dog new tricks.     

A cluster of carboxylic groups in PsbO protein is involved in proton transfer from the water oxidizing complex of Photosystem II

The hypothesis presented here for proton transfer away from the water oxidation complex of Photosystem II (PSII) is supported by biochemical experiments on the isolated PsbO protein in solution, theoretical analyses of better understood proton transfer systems like bacteriorhodopsin and cytochrome oxidase, and the recently published 3D structure of PS II (Pdb entry 1S5L). We propose that a cluster of conserved glutamic and aspartic acid residues in the PsbO protein acts as a buffering network providing efficient acceptors of protons derived from substrate water molecules. The charge delocalization of the cluster ensures readiness to promptly accept the protons liberated from substrate water. Therefore protons generated at the catalytic centre of PSII need not be released into the thylakoid lumen as generally thought. The cluster is the beginning of a localized, fast proton transfer conduit on the lumenal side of the thylakoid membrane. Proton-dependent conformational changes of PsbO may play a role in the regulation of both supply of substrate water to the water oxidizing complex and the resultant proton transfer.     

A cluster of carboxylic groups in PsbO protein is involved in proton transfer from the water oxidizing complex of Photosystem II

The hypothesis presented here for proton transfer away from the water oxidation complex of Photosystem II (PSII) is supported by biochemical experiments on the isolated PsbO protein in solution, theoretical analyses of better understood proton transfer systems like bacteriorhodopsin and cytochrome oxidase, and the recently published 3D structure of PS II (Pdb entry 1S5L). We propose that a cluster of conserved glutamic and aspartic acid residues in the PsbO protein acts as a buffering network providing efficient acceptors of protons derived from substrate water molecules. The charge delocalization of the cluster ensures readiness to promptly accept the protons liberated from substrate water. Therefore protons generated at the catalytic centre of PSII need not be released into the thylakoid lumen as generally thought. The cluster is the beginning of a localized, fast proton transfer conduit on the lumenal side of the thylakoid membrane. Proton-dependent conformational changes of PsbO may play a role in the regulation of both supply of substrate water to the water oxidizing complex and the resultant proton transfer.     

Water molecules and hydrogen-bonded networks in bacteriorhodopsin--molecular dynamics simulations of the ground state and the M-intermediate

Protein crystallography provides the structure of a protein, averaged over all elementary cells during data collection time. Thus, it has only a limited access to diffusive processes. This article demonstrates how molecular dynamics simulations can elucidate structure-function relationships in bacteriorhodopsin (bR) involving water molecules. The spatial distribution of water molecules and their corresponding hydrogen-bonded networks inside bR in its ground state (G) and late M intermediate conformations were investigated by molecular dynamics simulations. The simulations reveal a much higher average number of internal water molecules per monomer (28 in the G and 36 in the M) than observed in crystal structures (18 and 22, respectively). We found nine water molecules trapped and 19 diffusive inside the G-monomer, and 13 trapped and 23 diffusive inside the M-monomer. The exchange of a set of diffusive internal water molecules follows an exponential decay with a 1/e time in the order of 340 ps for the G state and 460 ps for the M state. The average residence time of a diffusive water molecule inside the protein is approximately 95 ps for the G state and 110 ps for the M state. We have used the Grotthuss model to describe the possible proton transport through the hydrogen-bonded networks inside the protein, which is built up in the picosecond-to-nanosecond time domains. Comparing the water distribution and hydrogen-bonded networks of the two different states, we suggest possible pathways for proton hopping and water movement inside bR.     

Water and bacteriorhodopsin: structure, dynamics, and function

A wealth of information has been gathered during the past decades that water molecules do play an important role in the structure, dynamics, and function of bacteriorhodopsin (bR) and purple membrane. Light-induced structural alterations in bR as detected by X-ray and neutron diffraction at low and high resolution are discussed in relationship to the mechanism of proton pumping. The analysis of high resolution intermediate structures revealed photon-induced rearrangements of water molecules and hydrogen bonds concomitant with conformational changes in the chromophore and the protein. These observations led to an understanding of key features of the pumping mechanism, especially the vectoriality and the different modes of proton translocation in the proton release and uptake domain of bR. In addition, water molecules influence the function of bR via equilibrium fluctuations, which must occur with adequate amplitude so that energy barriers between conformational states can be overcome.     

Strongly hydrogen-bonded water molecules in the Schiff base region of rhodopsins.

In many rhodopsins, a positively charged retinal chromophore is stabilized by a negatively charged carboxylate, and the presence of bound water molecules has been found in the Schiff base region by X-ray crystallography of various rhodopsins. Low-temperature Fourier-transform infrared (FTIR) spectroscopy can directly monitor hydrogen-bonding alterations of internal water molecules of rhodopsins. In particular, we found that a bridged water molecule between the Schiff base and Asp 85 in bacteriorhodopsin (BR), a light-driven proton-pump protein, forms an extremely strong hydrogen bond. It is likely that a hydration switch of the water from Asp 85 to Asp 212 plays an important role in the proton transfer in the Schiff base region of BR. Comprehensive studies of archaeal and visual rhodopsins have revealed that strongly hydrogen-bonded water molecules are only found in the proteins exhibiting proton-pump activities. Strongly hydrogen-bonded water molecules and its transient weakening may be essential for the proton-pump function of rhodopsins.     

Atomic resolution structures of bacteriorhodopsin photocycle intermediates: the role of discrete water molecules in the function of this light-driven ion pump

High-resolution X-ray crystallographic studies of bacteriorhodopsin have tremendously advanced our understanding of this light-driven ion pump during the last 2 years, and emphasized the crucial role of discrete internal water molecules in the pump cycle. In the extracellular region an extensive three-dimensional hydrogen-bonded network of protein residues and seven water molecules leads from the buried retinal Schiff base via water 402 and the initial proton acceptor Asp85 to the membrane surface. Near Lys216 where the retinal binds, transmembrane helix G contains a pi-bulge that causes a non-proline kink. The bulge is stabilized by hydrogen bonding of the main chain carbonyl groups of Ala215 and Lys216 with two buried water molecules located in the otherwise very hydrophobic region between the Schiff base and the proton donor Asp96 in the cytoplasmic region. The M intermediate trapped in the D96N mutant corresponds to a late M state in the transport cycle, after protonation of Asp85 and release of a proton to the extracellular membrane surface, but before reprotonation of the deprotonated retinal Schiff base. The M intermediate from the E204Q mutant corresponds to an earlier M, as in this mutant the Schiff base deprotonates without proton release. The structures of these two M states reveal progressive displacements of the retinal, main chain and side chains induced by photoisomerization of the retinal to 13-cis,15-anti, and an extensive rearrangement of the three-dimensional network of hydrogen-bonded residues and bound water that accounts for the changed pK(a)s of the Schiff base, Asp85, the proton release group and Asp96. The structure for the M state from E204Q suggests, moreover, that relaxation of the steric conflicts of the distorted 13-cis,15-anti retinal plays a critical role in the reprotonation of the Schiff base by Asp96. Two additional waters now connect Asp96 to the carbonyl of residue 216, in what appears to be the beginning of a hydrogen-bonded chain that would later extend to the retinal Schiff base. Based on the ground state and M intermediate structures, models of the molecular events in the early part of the photocycle are presented, including a novel model which proposes that bacteriorhodopsin pumps hydroxide (OH(-)) ions from the extracellular to the cytoplasmic side.     

Dynamic water networks in cytochrome C oxidase from Paracoccus denitrificans investigated by molecular dynamics simulations

We present a molecular dynamics study of cytochrome c oxidase from Paracoccus denitrificans in the fully oxidized state, embedded in a fully hydrated dimyristoylphosphatidylcholine lipid bilayer membrane. Parallel simulations with different levels of protein hydration, 1.125 ns each in length, were carried out under conditions of constant temperature and pressure using three-dimensional periodic boundary conditions and full electrostatics to investigate the distribution and dynamics of water molecules and their corresponding hydrogen-bonded networks inside cytochrome c oxidase. The majority of the water molecules had residence times shorter than 100 ps, but a few water molecules are fixed inside the protein for up to 1.125 ns. The hydrogen-bonded network in cytochrome c oxidase is not uniformly distributed, and the degree of water arrangement is variable. The average number of solvent sites in the proton-conducting K- and D-pathways was determined. In contrast to single water files in narrow geometries we observe significant diffusion of individual water molecules along these pathways. The highly fluctuating hydrogen-bonded networks, combined with the significant diffusion of individual water molecules, provide a basis for the transfer of protons in cytochrome c oxidase, therefore leading to a better understanding of the mechanism of proton pumping.     

pK(a) Calculations suggest storage of an excess proton in a hydrogen-bonded water network in bacteriorhodopsin.

Calculations of protonation states and pK(a) values for the ionizable groups in the resting state of bacteriorhodopsin have been carried out using the recently available 1.55 A resolution X-ray crystallographic structure. The calculations are in reasonable agreement with the available experimental data for groups on or near the ion transport chain (the retinal Schiff base; Asp85, 96, 115, 212, and Arg82). In contrast to earlier studies using lower-resolution structural data, this agreement is achieved without manipulations of the crystallographically determined heavy-atom positions or ad hoc adjustments of the intrinsic pK(a) of the Schiff base. Thus, the theoretical methods used provide increased reliability as the input structural data are improved. Only minor effects on the agreement with experiment are found with respect to methodological variations, such as single versus multi-conformational treatment of hydrogen atom placements, or retaining the crystallographically determined internal water molecules versus treating them as high-dielectric cavities. The long-standing question of the identity of the group that releases a proton to the extracellular side of the membrane during the L-to-M transition of the photocycle is addressed by including as pH-titratable sites not only Glu204 and Glu194, residues near the extracellular side that have been proposed as the release group, but also an H(5)O(2)(+) molecule in a nearby cavity. The latter represents the recently proposed storage of the release proton in an hydrogen-bonded water network. In all calculations where this possibility is included, the proton is stored in the H(5)O(2)(+) rather than on either of the glutamic acids, thus establishing the plausibility on theoretical grounds of the storage of the release proton in bacteriorhodopsin in a hydrogen-bonded water network. The methods used here may also be applicable to other proteins that may store a proton in this way, such as the photosynthetic reaction center and cytochrome c oxidase.     

How environment supports a state: molecular dynamics simulations of two states in bacteriorhodopsin suggest lipid and water compensation

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.