Proton permeation of lipid bilayers

Proton permeation of the lipid bilayer barrier has two unique features. First, permeability coefficients measured at neutral pH ranges are six to seven orders of magnitude greater than expected from knowledge of other monovalent cations. Second, proton conductance across planar lipid bilayers varies at most by a factor of 10 when pH is varied from near 1 to near 11. Two mechanisms have been proposed to account for this anomalous behavior: proton conductance related to contaminants of lipid bilayers, and proton translocation along transient hydrogen-bonded chains (tHBC) of associated water molecules in the membrane. The weight of evidence suggests that trace contaminants may contribute to proton conductance across planar lipid membranes at certain pH ranges, but cannot account for the anomalous proton flux in liposome systems.Two new results will be reported here which were designed to test the tHBC model. These include measurements of relative proton/potassium permeability in the gramicidin channel, and plots of proton flux against the magnitude of pH gradients. (1) The relative permeabilities of protons and potassium through the gramicidin channel, which contains a single strand of hydrogenbonded water molecules, were found to differ by at least four orders of magnitude when measured at neutral pH ranges. This result demonstrates that a hydrogen-bonded chain of water molecules can provide substantial discrimination between protons and other cations. It was also possible to calculate that if approximately 7% of bilayer water was present in a transient configuration similar to that of the gramicidin channel, it could account for the measured proton flux. (2) The plot of proton conductance against pH gradient across liposome membranes was superlinear, a result that is consistent with one of three alternative tHBC models for proton conductance described by Nagle elsewhere in this volume.     

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

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.     

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.     

Attenuation of proton currents by methanol in a dioxolane-linked gramicidin A channel in different lipid bilayers.

The mobility of protons in a dioxolane-linked gramicidin A channel (D1) is comparable to the mobility of protons in aqueous solutions (Cukierman, S., E. P. Quigley, and D. S. Crumrine. 1997. Biophys. J. 73:2489-2502). Aliphatic alcohols decrease the mobility of H+ in aqueous solutions. In this study, the effects of methanol on proton conduction through D1 channels were investigated in different lipid bilayers and at different HCl concentrations. Methanol attenuated H+ currents in a voltage-independent manner. Attenuation of proton currents was also independent of H+ concentrations in solution. In phospholipid bilayers, methanol decreased the single channel conductance to protons without affecting the binding affinity of protons to bilayers. In glycerylmonooleate membranes, the attenuation of single channel proton conductances qualitatively resembled the decrease of conductivities of HCl solutions by methanol. However, in both types of lipid bilayers, single channel proton conductances through D1 channels were considerably more attenuated than the conductivities of different HCl solutions. This suggests that methanol modulates single proton currents through D1 channels. It is proposed that, on average, one methanol molecule binds to a D1 channel, and attenuates H+ conductance. The Gibbs free energy of this process (DeltaG0) is approximately 1.2 kcal/mol, which is comparable to the free energy of decrease of HCl conductivity in methanol solutions (1.6 kcal/mol). Apolar substances like urea and glucose that do not transport protons in HCl solutions and do not permeate D1 channels decreased solution conductivity and single channel conductance by a considerably larger proportion than methanol. Cs+ currents through D1 channels were considerably less (fivefold) attenuated by methanol than proton currents. It is proposed that methanol partitions inside the pore of gramicidin channels and delays the transfer of protons between water and methanol molecules, causing a significant attenuation of the single channel proton conductance. Gramicidin channels offer an interesting experimental model to study proton hopping along a single chain of water molecules interrupted by a single methanol molecule.     

Three ways in, one way out: water dynamics in the trans-membrane domains of the inner membrane translocase AcrB

Powered by proton-motive force, the inner membrane translocase AcrB is the engine of the AcrAB-TolC efflux pump in Escherichia coli. As proton conduction in proteins occurs along hydrogen-bonded networks of polar residues and water molecules, knowledge of the protein-internal water distribution and water-interacting residues allows drawing conclusions to possible pathways of proton conduction. Here, we report a series of 6× 50 ns independent molecular dynamics simulations of asymmetric AcrB embedded in a phospholipid/water environment. Simulating each monomer in its proposed protonation state, we calculated for each trans-membrane domain the average water distribution, identified residues interacting with these waters and quantified each residue's frequency of water hydrogen bond contact. Combining this information we find three possible routes of proton transfer connecting a continuously hydrated region of known key residues in the TMD interior to bulk water by one cytoplasmic and up to three periplasm water channels in monomer B and A. We find that water access of the trans-membrane domains is regulated by four groups of residues in a combination of side chain re-orientations and shifts of trans-membrane helices. Our findings support a proton release event via Arg971 during the C intermediate or in the transition to A, and proton uptake occurring in the A or B state or during a so far unknown intermediate in between B and C where cytoplasmic water access is still possible. Our simulations suggest experimentally testable hypotheses, which have not been investigated so far.     

Protein, lipid and water organization in bacteriorhodopsin crystals: a molecular view of the purple membrane at 1.9 A resolution.

BACKGROUND: Bacteriorhodopsin (bR) from Halobacterium salinarum is a proton pump that converts the energy of light into a proton gradient that drives ATP synthesis. The protein comprises seven transmembrane helices and in vivo is organized into purple patches, in which bR and lipids form a crystalline two-dimensional array. Upon absorption of a photon, retinal, which is covalently bound to Lys216 via a Schiff base, is isomerized to a 13-cis,15-anti configuration. This initiates a sequence of events - the photocycle - during which a proton is transferred from the Schiff base to Asp85, followed by proton release into the extracellular medium and reprotonation from the cytoplasmic side. RESULTS: The structure of bR in the ground state was solved to 1.9 A resolution from non-twinned crystals grown in a lipidic cubic phase. The structure reveals eight well-ordered water molecules in the extracellular half of the putative proton translocation pathway. The water molecules form a continuous hydrogen-bond network from the Schiff-base nitrogen (Lys216) to Glu194 and Glu204 and includes residues Asp85, Asp212 and Arg82. This network is involved both in proton translocation occurring during the photocycle, as well as in stabilizing the structure of the ground state. Nine lipid phytanyl moieties could be modeled into the electron-density maps. Matrix-assisted laser desorption/ionization mass spectrometry (MALDI-MS) analysis of single crystals demonstrated the presence of four different charged lipid species. CONCLUSIONS: The structure of protein, lipid and water molecules in the crystals represents the functional entity of bR in the purple membrane of the bacteria at atomic resolution. Proton translocation from the Schiff base to the extracellular medium is mediated by a hydrogen-bond network that involves charged residues and water molecules.     

Molecular dynamics simulation of lipid packing and mobility in bilayer membranes

A model of lipid bilayer membrane in water has been developed. Parameters have been selected that allow molecular dynamics simulation of lipid bilayers in the all-atom approximation. The calculated indices of packing and mobility of lipid molecules for the liquid crystalline state of the bilayer agree well with the experimental data. Based on the model of the liquid crystalline state of the membrane, a system in the gel-like state has been constructed. The gel-state model reproduces well the packing of lipids in real bilayers, whereas the mobility of molecules proves to be overestimated.     

Proton conduction in lecithins. II. Effect of hydration

The effect of hydration on the electrical properties of lecithins was studied using dielectric spectroscopy in the low and audio frequency range. The samples are characterized by two impedances which are related to electrode and bulk processes. Adsorption of water increases the bulk conduction by serval orders of magnitude, while its influence on electrode conduction is less and the electrode capacitance is practically unchanged. The dependence of the sample impedance on the electric field was studied. The results confirm the existence of electrode events, charge transport at the electrode is accelerated by the electric field. The temperature dependence of the bulk conduction was measured. The activation energy of conduction depends on the phase state and is in connection with the mobility of the lipid molecules. An ionic (proton) mechanism of conduction is suggested which involves head group rotation.     

Computer simulation of ion channel gating: the M(2) channel of influenza A virus in a lipid bilayer

The transmembrane fragment of the influenza virus M(2) protein forms a homotetrameric channel that transports protons. In this paper, we use molecular dynamics simulations to help elucidate the mechanism of channel gating by four histidines that occlude the channel lumen in the closed state. We test two competing hypotheses. In the "shuttle" mechanism, the delta nitrogen atom on the extracellular side of one histidine is protonated by the incoming proton, and, subsequently, the proton on the epsilon nitrogen atom is released on the opposite side. In the "water-wire" mechanism, the gate opens because of electrostatic repulsion between four simultaneously biprotonated histidines. This allows for proton transport along the water wire that penetrates the gate. For each system, composed of the channel embedded in a hydrated phospholipid bilayer, a 1.3-ns trajectory was obtained. It is found that the states involved in the shuttle mechanism, which contain either single-protonated histidines or a mixture of single-protonated histidines plus one biprotonated residue, are stable during the simulations. Furthermore, the orientations and dynamics of water molecules near the gate are conducive to proton transfer. In contrast, the fully biprotonated state is not stable. Additional simulations show that if only two histidines are biprotonated, the channel deforms but the gate remains closed. These results support the shuttle mechanism but not the gate-opening mechanism of proton gating in M(2).     

Structural changes in lipid membranes and collagen irradiated with UV light and the protective effect of plant extracts

The results of experimental studies on the effect of UV irradiation on collagen, artificial lipid membranes, and rat skin, as well as the protective effect of plant extracts from UV radiation are presented. The irradiation of collagen and lipid membranes with solar and artificial UV light leads to structural changes in these objects. In particular, collagen molecules denature and transfer into a new conformational state. The effect of UV light on lipid membranes and liposomes leads to a disturbance of membrane structure, which is connected with a decrease in the number of lipid molecules involved in the cooperative transition from gel into a liquid crystal state. The components of plant extracts (mainly flavonoids) absorb UV radiation in the erythem-forming spectral area and block the destructive processes occurring in collagen and lipids.     

Pore Hydration States of KcsA Potassium Channels in Membranes

Water-filled hydrophobic cavities in channel proteins serve as gateways for transfer of ions across membranes, but their properties are largely unknown. We determined water distributions along the conduction pores in two tetrameric channels embedded in lipid bilayers using neutron diffraction: potassium channel KcsA and the transmembrane domain of M2 protein of influenza A virus. For the KcsA channel in the closed state, the distribution of water is peaked in the middle of the membrane, showing water in the central cavity adjacent to the selectivity filter. This water is displaced by the channel blocker tetrabutyl-ammonium. The amount of water associated with the channel was quantified, using neutron diffraction and solid state NMR. In contrast, the M2 proton channel shows a V-shaped water profile across the membrane, with a narrow constriction at the center, like the hourglass shape of its internal surface. These two types of water distribution are therefore very different in their connectivity to the bulk water. The water and protein profiles determined here provide important evidence concerning conformation and hydration of channels in membranes and the potential role of pore hydration in channel gating.     

Structure and dynamics of a proton wire: a theoretical study of H+ translocation along the single-file water chain in the gramicidin A channel.

The rapid translocation of H+ along a chain of hydrogen-bonded water molecules, or proton wire, is thought to be an important mechanism for proton permeation through transmembrane channels. Computer simulations are used to study the properties of the proton wire formed by the single-file waters in the gramicidin A channel. The model includes the polypeptidic dimer, with 22 water molecules and one excess proton. The dissociation of the water molecules is taken into account by the "polarization model" of Stillinger and co-workers. The importance of quantum effects due to the light mass of the hydrogen nuclei is examined with the use of discretized Feynman path integral molecular dynamics simulations. Results show that the presence of an excess proton in the pore orients the single-file water molecules and affects the geometry of water-water hydrogen bonding interactions. Rather than a well-defined hydronium ion OH3+ in the single-file region, the protonated species is characterized by a strong hydrogen bond resembling that of O2H5+. The quantum dispersion of protons has a small but significant effect on the equilibrium structure of the hydrogen-bonded water chain. During classical trajectories, proton transfer between consecutive water molecules is a very fast spontaneous process that takes place in the subpicosecond time scale. The translocation along extended regions of the chain takes place neither via a totally concerted mechanism in which the donor-acceptor pattern would flip over the entire chain in a single step, nor via a succession of incoherent hops between well-defined intermediates. Rather, proton transfer in the wire is a semicollective process that results from the subtle interplay of rapid hydrogen-bond length fluctuations along the water chain. These rapid structural fluctuations of the protonated single file of waters around an average position and the slow movements of the average position of the excess proton along the channel axis occur on two very different time scales. Ultimately, it is the slow reorganization of hydrogen bonds between single-file water molecules and channel backbone carbonyl groups that, by affecting the connectivity and the dynamics of the single-file water chain, also limits the translocation of the proton across the pore.     

Water wires in atomistic models of the Hv1 proton channel

The voltage-gated proton channel (Hv1) is homologous to the voltage-sensing domain (VSD) of voltage-gated potassium (Kv) channels but lacks a separate pore domain. The Hv1 monomer has dual functions: it gates the proton current and also serves as the proton conduction pathway. To gain insight into the structure and dynamics of the yet unresolved proton permeation pathway, we performed all-atom molecular dynamics simulations of two different Hv1 homology models in a lipid bilayer in excess water. The structure of the Kv1.2-Kv2.1 paddle-chimera VSD was used as template to generate both models, but they differ in the sequence alignment of the S4 segment. In both models, we observe a water wire that extends through the membrane, whereas the corresponding region is dry in simulations of the Kv1.2-Kv2.1 paddle-chimera. We find that the kinetic stability of the water wire is dependent upon the identity and location of the residues lining the permeation pathway, in particular, the S4 arginines. A measurement of water transport kinetics indicates that the water wire is a relatively static feature of the permeation pathway. Taken together, our results suggest that proton conduction in Hv1 may occur via Grotthuss hopping along a robust water wire, with exchange of water molecules between inner and outer ends of the permeation pathway minimized by specific water-protein interactions. This article is part of a Special Issue entitled: Membrane protein structure and function.     

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.     

Molecular mechanism of H+ conduction in the single-file water chain of the gramicidin channel

The conduction of protons in the hydrogen-bonded chain of water molecules (or "proton wire") embedded in the lumen of gramicidin A is studied with molecular dynamics free energy simulations. The process may be described as a "hop-and-turn" or Grotthuss mechanism involving the chemical exchange (hop) of hydrogen nuclei between hydrogen-bonded water molecules arranged in single file in the lumen of the pore, and the subsequent reorganization (turn) of the hydrogen-bonded network. Accordingly, the conduction cycle is modeled by two complementary steps corresponding respectively to the translocation 1) of an ionic defect (H+) and 2) of a bonding defect along the hydrogen-bonded chain of water molecules in the pore interior. The molecular mechanism and the potential of mean force are analyzed for each of these two translocation steps. It is found that the mobility of protons in gramicidin A is essentially determined by the fine structure and the dynamic fluctuations of the hydrogen-bonded network. The translocation of H+ is mediated by spontaneous (thermal) fluctuations in the relative positions of oxygen atoms in the wire. In this diffusive mechanism, a shallow free-energy well slightly favors the presence of the excess proton near the middle of the channel. In the absence of H+, the water chain adopts either one of two polarized configurations, each of which corresponds to an oriented donor-acceptor hydrogen-bond pattern along the channel axis. Interconversion between these two conformations is an activated process that occurs through the sequential and directional reorientation of water molecules of the wire. The effect of hydrogen-bonding interactions between channel and water on proton translocation is analyzed from a comparison to the results obtained previously in a study of model nonpolar channels, in which such interactions were missing. Hydrogen-bond donation from water to the backbone carbonyl oxygen atoms lining the pore interior has a dual effect: it provides a coordination of water molecules well suited both to proton hydration and to high proton mobility, and it facilitates the slower reorientation or turn step of the Grotthuss mechanism by stabilizing intermediate configurations of the hydrogen-bonded network in which water molecules are in the process of flipping between their two preferred, polarized states. This mechanism offers a detailed molecular model for the rapid transport of protons in channels, in energy-transducing membrane proteins, and in enzymes.     

Proton conduction in gramicidin A and in its dioxolane-linked dimer in different lipid bilayers.

Gramicidin A (gA) molecules were covalently linked with a dioxolane ring. Dioxolane-linked gA dimers formed ion channels, selective for monovalent cations, in planar lipid bilayers. The main goal of this study was to compare the functional single ion channel properties of natural gA and its covalently linked dimer in two different lipid bilayers and HCl concentrations (10-8000 mM). Two ion channels with different gating and conductance properties were identified in bilayers from the product of dimerization reaction. The most commonly observed and most stable gramicidin A dimer is the main object of this study. This gramicidin dimer remained in the open state most of the time, with brief closing flickers (tau(closed) approximately 30 micros). The frequency of closing flickers increased with transmembrane potential, making the mean open time moderately voltage dependent (tau(open) changed approximately 1.43-fold/100 mV). Such gating behavior is markedly different from what is seen in natural gA channels. In PEPC (phosphatidylethanolamine-phosphatidylcholine) bilayers, single-channel current-voltage relationships had an ohmic behavior at low voltages, and a marked sublinearity at relatively higher voltages. This behavior contrasts with what was previously described in GMO (glycerylmonooleate) bilayers. In PEPC bilayers, the linear conductance of single-channel proton currents at different proton concentrations was essentially the same for both natural and gA dimers. g(max) and K(D), obtained from fitting experimental points to a Langmuir adsorption isotherm, were approximately 1500 pS and 300 mM, respectively, for both the natural gA and its dimer. In GMO bilayers, however, proton affinities of gA and the dioxolane-dimer were significantly lower (K(D) of approximately 1 and 1.5 M, respectively), and the g(max) higher (approximately 1750 and 2150 pS, respectively) than in PEPC bilayers. Furthermore, the relationship between single-channel conductance and proton concentration was linear at low bulk concentrations of H+ (0.01-2 M) and saturated at concentrations of more than 3 M. It is concluded that 1) The mobility of protons in gramicidin A channels in different lipid bilayers is remarkably similar to proton mobilities in aqueous solutions. In particular, at high concentrations of HCl, proton mobilities in gramicidin A channel and in solution differ by only 25%. 2) Differences between proton conductances in gramicidin A channels in GMO and PEPC cannot be explained by surface charge effects on PEPC membranes. It is proposed that protonated phospholipids adjacent to the mouth of the pore act as an additional source of protons for conduction through gA channels in relation to GMO bilayers. 3) Some experimental results cannot be reconciled with simple alterations in access resistance to proton flow in gA channels. Said differences could be explained if the structure and/or dynamics of water molecules inside gramicidin A channels is modulated by the lipid environment and by modifications in the structure of gA channels. 4) The dioxolane ring is probably responsible for the closing flickers seen in the dimer channel. However, other factors can also influence closing flickers.     

Structural insights into the lipid transfer mechanism of a non‐specific lipid transfer protein

The non‐specific lipid transfer proteins (nsLTPs) are multifunctional seed proteins engaged in several different physiological processes. The nsLTPs are stabilized by four disulfide bonds and exhibit a characteristic hydrophobic cavity, which is the primary lipid binding site. While these proteins are known to transfer lipids between membranes, the mechanism of lipid transfer has remained elusive. Four crystal structures of nsLTP from Solanum melongena, one in the apo‐state and three myristic acid bound states were determined. Among the three lipid bound states, two lipid molecules were bound on the nsLTP surface at different positions and one was inside the cavity. The lipid‐dependent conformational changes leading to opening of the cavity were revealed based on structural and spectroscopic data. The surface‐bound lipid represented a transient intermediate state and the lipid ultimately moved inside the cavity through the cavity gate as revealed by molecular dynamics simulations. Two critical residues in the loop regions played possible ‘gating’ role in the opening and closing of the cavity. Antifungal activity and membrane permeabilization effect of nsLTP against Fusarium oxysporum suggested that it could possibly involve in bleaching out the lipids. Collectively, these studies support a model of lipid transfer mechanism by nsLTP via intermediate states.     

Structural and dynamic mechanisms for the function and inhibition of the M2 proton channel from influenza A virus

The M2 proton channel from influenza A virus, a prototype for a class of viral ion channels known as viroporins, conducts protons along a chain of water molecules and ionizable sidechains, including His37. Recent studies highlight a delicate interplay between protein folding, proton binding, and proton conduction through the channel. Drugs inhibit proton conduction by binding to an aqueous cavity adjacent to M2's proton-selective filter, thereby blocking access of proton to the filter, and altering the energetic landscape of the channel and the energetics of proton-binding to His37.     

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.