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

Semiaromatic Polyamide-Based Membrane in Forward Osmosis: Molecular Insights

Despite the increased interest in forward osmosis (FO) in recent years, the technology's advancement in commercial and industrial applications has been hampered by the absence of suitable FO membranes and ideal draw solutes, which demands the exploration of new membranes and novel draw solutes targeted for some specific applications. In this context, we considered a semiaromatic polyamide (SAPA) for an application where monovalent salt can be permeated but has high selectivity toward divalent salt and excellent water permeability. In this regard, we constructed an atomistic model for the membrane via a heuristic approach using an equilibrated mixture of hydrolyzed trimesoyl chloride and piperazine monomers and performed nonequilibrium molecular dynamics simulations on the SAPA membrane in the FO process to understand the structural properties and performance of the membrane at the atomistic level. We used pure water as the feed and Na₂SO₄ as the draw solution. It is observed that the SAPA membrane shows excellent water permeability and no reverse draw solute flux. To further test the dynamics of salt ions inside the membranes, we performed two distinct equilibrium simulations on systems consisting of either monovalent salt, such as NaCl, or divalent salt, such as Na₂SO₄. The atomistic details of the interactions between the functional groups of the membrane and salt ions provided in this work can inspire further experiments on SAPA membranes in the context of separation of monovalent and divalent salts, which have applications in the treatment of textile industry wastewater.

Archaeal Lipids Regulating the Trimeric Structure Dynamics of Bacteriorhodopsin for Efficient Proton Release and Uptake

S-TGA-1 and PGP-Me are native archaeal lipids associated with the bacteriorhodopsin (bR) trimer and contribute to protein stabilization and native dynamics for proton transfer. However, little is known about the underlying molecular mechanism of how these lipids regulate bR trimerization and efficient photocycling. Here, we explored the specific binding of S-TGA-1 and PGP-Me with the bR trimer and elucidated how specific interactions modulate the bR trimeric structure and proton release and uptake using long-term atomistic molecular dynamic simulations. Our results showed that S-TGA-1 and PGP-Me are essential for stabilizing the bR trimer and maintaining the coherent conformational dynamics necessary for proton transfer. The specific binding of S-TGA-1 with W80 and K129 regulates proton release on the extracellular surface by forming a “Glu-shared” model. The interaction of PGP-Me with K40 ensures proton uptake by accommodating the conformation of the helices to recruit enough water molecules on the cytoplasmic side. The present study results could fill in the theoretical gaps of studies on the functional role of archaeal lipids and could provide a reference for other membrane proteins containing similar archaeal lipids.

Long-distance perturbation on Schiff base-counterion interactions by His30 and the extracellular Na+-binding site in Krokinobacter rhodopsin 2

Krokinobacter rhodopsin 2 (KR2), a light-driven Na+ pump, is a dual-functional protein, pumping protons in the absence of Na+ when K+ or larger alkali metal ions are present. A specific mutation in helix A near the extracellular Na+ binding site, H30A, eliminates its proton pumping ability. We induced structural changes in H30A by altering the alkali metal ion bound at the extracellular binding site, and observed a strong electrostatic interaction between the Schiff base and counterion and torsion around the Schiff base as revealed by solid-state nuclear magnetic resonance (NMR) and Fourier transform infrared (FTIR) spectroscopies. The strong interaction when His30 was absent and no ion bound at the extracellular binding site disabled retinal reisomerization, as was shown with flash-photolysis, forming a small amount of only a K-like intermediate. This revealed why H30A lacks the proton pumping function. Long-distance perturbation of the binding site and Schiff base revealed that a non-transported ion binding at the extracellular site is essential for pumping.

Microbial Rhodopsins

Microbial rhodopsins (MRs) are a large family of photoactive membrane proteins, found in microorganisms belonging to all kingdoms of life, with new members being constantly discovered. Among the MRs are light-driven proton, cation and anion pumps, light-gated cation and anion channels, and various photoreceptors. Due to their abundance and amenability to studies, MRs served as model systems for a great variety of biophysical techniques, and recently found a great application as optogenetic tools. While the basic aspects of microbial rhodopsins functioning have been known for some time, there is still a plenty of unanswered questions. This chapter presents and summarizes the available knowledge, focusing on the functional and structural studies.

How Mg2+ ion and water network affect the stability and structure of non-Watson-Crick base pairs in E. coli loop E of 5S rRNA: a molecular dynamics and reference interaction site model (RISM) study

The non-Watson-Crick (non-WC) base pairs of Escherichia coli loop E of 5S rRNA are stabilized by Mg²⁺ ions through water-mediated interaction. It is important to know the synergic role of Mg²⁺ and the water network surrounding Mg²⁺ in stabilizing the non-WC base pairs of RNA. For this purpose, free energy change of the system is calculated using molecular dynamics (MD) simulation as Mg²⁺ is pulled from RNA, which causes disturbance of the water network. It was found that Mg²⁺ remains hexahydrated unless it is close to or far from RNA. In the pentahydrated form, Mg²⁺ interacts directly with RNA. Water network has been identified by two complimentary methods; MD followed by a density-based clustering algorithm and three-dimensional-reference interaction site model. These two methods gave similar results. Identification of water network around Mg²⁺ and non-WC base pairs gives a clue to the strong effect of water network on the stability of this RNA. Based on sequence analysis of all Eubacteria 5s rRNA, we propose that hexahydrated Mg²⁺ is an integral part of this RNA and geometry of base pairs surrounding it adjust to accommodate the [Formula: see text]. Overall the findings from this work can help in understanding the basis of the complex structure and stability of RNA with non-WC base pairs.

Proton-coupled electron transfer and the role of water molecules in proton pumping by cytochrome c oxidase

Molecular oxygen acts as the terminal electron sink in the respiratory chains of aerobic organisms. Cytochrome c oxidase in the inner membrane of mitochondria and the plasma membrane of bacteria catalyzes the reduction of oxygen to water, and couples the free energy of the reaction to proton pumping across the membrane. The proton-pumping activity contributes to the proton electrochemical gradient, which drives the synthesis of ATP. Based on kinetic experiments on the O-O bond splitting transition of the catalytic cycle (A → P(R)), it has been proposed that the electron transfer to the binuclear iron-copper center of O2 reduction initiates the proton pump mechanism. This key electron transfer event is coupled to an internal proton transfer from a conserved glutamic acid to the proton-loading site of the pump. However, the proton may instead be transferred to the binuclear center to complete the oxygen reduction chemistry, which would constitute a short-circuit. Based on atomistic molecular dynamics simulations of cytochrome c oxidase in an explicit membrane-solvent environment, complemented by related free-energy calculations, we propose that this short-circuit is effectively prevented by a redox-state-dependent organization of water molecules within the protein structure that gates the proton transfer pathway.

Nanoscale electron transport and photodynamics enhancement in lipid-depleted bacteriorhodopsin monomers

Potential future use of bacteriorhodopsin (bR) as a solid-state electron transport (ETp) material requires the highest possible active protein concentration. To that end we prepared stable monolayers of protein-enriched bR on a conducting HOPG substrate by lipid depletion of the native bR. The ETp properties of this construct were then investigated using conducting probe atomic force microscopy at low bias, both in the ground dark state and in the M-like intermediate configuration, formed upon excitation by green light. Photoconductance modulation was observed upon green and blue light excitation, demonstrating the potential of these monolayers as optoelectronic building blocks. To correlate protein structural changes with the observed behavior, measurements were made as a function of pressure under the AFM tip, as well as humidity. The junction conductance is reversible under pressure changes up to ∼300 MPa, but above this pressure the conductance drops irreversibly. ETp efficiency is enhanced significantly at >60% relative humidity, without changing the relative photoactivity significantly. These observations are ascribed to changes in protein conformation and flexibility and suggest that improved electron transport pathways can be generated through formation of a hydrogen-bonding network.

Exploring cavity dynamics in biomolecular systems

Background

The internal cavities of proteins are dynamic structures and their dynamics may be associated with conformational changes which are required for the functioning of the protein. In order to study the dynamics of these internal protein cavities, appropriate tools are required that allow rapid identification of the cavities as well as assessment of their time-dependent structures.

Results

In this paper, we present such a tool and give results that illustrate the applicability for the analysis of molecular dynamics trajectories. Our algorithm consists of a pre-processing step where the structure of the cavity is computed from the Voronoi diagram of the van der Waals spheres based on coordinate sets from the molecular dynamics trajectory. The pre-processing step is followed by an interactive stage, where the user can compute, select and visualize the dynamic cavities. Importantly, the tool we discuss here allows the user to analyze the time-dependent changes of the components of the cavity structure. An overview of the cavity dynamics is derived by rendering the dynamic cavities in a single image that gives the cavity surface colored according to its time-dependent dynamics.

Conclusion

The Voronoi-based approach used here enables the user to perform accurate computations of the geometry of the internal cavities in biomolecules. For the first time, it is possible to compute dynamic molecular paths that have a user-defined minimum constriction size. To illustrate the usefulness of the tool for understanding protein dynamics, we probe the dynamic structure of internal cavities in the bacteriorhodopsin proton pump.

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.

Direct observation of self-assembled chain-like water structures in a nanoscopic water meniscus

Sawtooth-like oscillatory forces generated by water molecules confined between two oxidized silicon surfaces were observed using a cantilever-based optical interfacial force microscope when the two surfaces approached each other in ambient environments. The humidity-dependent oscillatory amplitude and periodicity were 3-12 nN and 3-4 water diameters, respectively. Half of each period was matched with a freely jointed chain model, possibly suggesting that the confined water behaved like a bundle of water chains. The analysis also indicated that water molecules self-assembled to form chain-like structures in a nanoscopic meniscus between two hydrophilic surfaces in air. From the friction force data measured simultaneously, the viscosity of the chain-like water was estimated to be between 10(8) and 10(10) times greater than that of bulk water. The suggested chain-like structure resolves many unexplained properties of confined water at the nanometer scale, thus dramatically improving the understanding of a variety of water systems in nature.

Deprotonation of D96 in bacteriorhodopsin opens the proton uptake pathway

Despite extensive investigation, the precise mechanism controlling the opening of the cytoplasmic proton uptake pathway in bacteriorhodopsin (bR) has remained a mystery. From an analysis of the X-ray structure of the D96G/F171C/F219L triple mutant of bR and 60 independent molecular dynamics simulations of bR photointermediates, we report that the deprotonation of D96, a key residue in proton transfer reactions, serves two roles that occur sequentially. First, D96 donates a proton to the Schiff base. Subsequently, the deprotonation of D96 serves to "unlatch" the cytoplasmic side. The latching function of D96 appears to be remarkably robust, functioning to open hydration channels in all photointermediate structures. These results suggest that the protonation state of D96 may be the critical biophysical cue controlling the opening and closing of the cytoplasmic half-channel in bR. We suspect that this protonation-switch mechanism could also be utilized in other proton pumps to minimize backflow and reinforce directionality.

Molecular Dynamics Computer Simulations of Multidrug RND Efflux Pumps

Over-expression of multidrug efflux pumps of the Resistance Nodulation Division (RND) protein super family counts among the main causes for microbial resistance against pharmaceuticals. Understanding the molecular basis of this process is one of the major challenges of modern biomedical research, involving a broad range of experimental and computational techniques. Here we review the current state of RND transporter investigation employing molecular dynamics simulations providing conformational samples of transporter components to obtain insights into the functional mechanism underlying efflux pump-mediated antibiotics resistance in Escherichia coli and Pseudomonas aeruginosa.

Molecular dynamics simulations reveal proton transfer pathways in cytochrome C-dependent nitric oxide reductase

Nitric oxide reductases (NORs) are membrane proteins that catalyze the reduction of nitric oxide (NO) to nitrous oxide (N₂O), which is a critical step of the nitrate respiration process in denitrifying bacteria. Using the recently determined first crystal structure of the cytochrome c-dependent NOR (cNOR) [Hino T, Matsumoto Y, Nagano S, Sugimoto H, Fukumori Y, et al. (2010) Structural basis of biological N2O generation by bacterial nitric oxide reductase. Science 330: 1666–70.], we performed extensive all-atom molecular dynamics (MD) simulations of cNOR within an explicit membrane/solvent environment to fully characterize water distribution and dynamics as well as hydrogen-bonded networks inside the protein, yielding the atomic details of functionally important proton channels. Simulations reveal two possible proton transfer pathways leading from the periplasm to the active site, while no pathways from the cytoplasmic side were found, consistently with the experimental observations that cNOR is not a proton pump. One of the pathways, which was newly identified in the MD simulation, is blocked in the crystal structure and requires small structural rearrangements to allow for water channel formation. That pathway is equivalent to the functional periplasmic cavity postulated in cbb₃ oxidase, which illustrates that the two enzymes share some elements of the proton transfer mechanisms and confirms a close evolutionary relation between NORs and C-type oxidases. Several mechanisms of the critical proton transfer steps near the catalytic center are proposed.

Denitrification is an anaerobic process performed by several bacteria as an alternative to aerobic respiration. A key intermediate step is catalyzed by the nitric oxide reductase (NOR) enzyme, which is situated in the cytoplasmic membrane. Proton delivery to the catalytic site inside NOR is an important part of its functioning. In this work we use molecular dynamics simulations to describe water distribution and to identify proton transfer pathways in cNOR. Our results reveal two channels from the periplasmic side of the membrane and none from the cytoplasmic side, indicating that cNOR is not a proton pump. It is our hope that these results will provide a basis for further experimental and computational studies aimed to understand details of the NOR mechanism. Furthermore, this work sheds light on the molecular evolution of respiratory enzymes.

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.

Forced unbinding of GPR17 ligands from wild type and R255I mutant receptor models through a computational approach

Background

GPR17 is a hybrid G-protein-coupled receptor (GPCR) activated by two unrelated ligand families, extracellular nucleotides and cysteinyl-leukotrienes (cysteinyl-LTs), and involved in brain damage and repair. Its exploitment as a target for novel neuro-reparative strategies depends on the elucidation of the molecular determinants driving binding of purinergic and leukotrienic ligands. Here, we applied docking and molecular dynamics simulations (MD) to analyse the binding and the forced unbinding of two GPR17 ligands (the endogenous purinergic agonist UDP and the leukotriene receptor antagonist pranlukast from both the wild-type (WT) receptor and a mutant model, where a basic residue hypothesized to be crucial for nucleotide binding had been mutated (R255I) to Ile.

Results

MD suggested that GPR17 nucleotide binding pocket is enclosed between the helical bundle and extracellular loop (EL) 2. The driving interaction involves R255 and the UDP phosphate moiety. To support this hypothesis, steered MD experiments showed that the energy required to unbind UDP is higher for the WT receptor than for R255I. Three potential binding sites for pranlukast where instead found and analysed. In one of its preferential docking conformations, pranlukast tetrazole group is close to R255 and phenyl rings are placed into a subpocket highly conserved among GPCRs. Pulling forces developed to break polar and aromatic interactions of pranlukast were comparable. No differences between the WT receptor and the R255I receptor were found for the unbinding of pranlukast.

Conclusions

These data thus suggest that, in contrast to which has been hypothesized for nucleotides, the lack of the R255 residue doesn't affect the binding of pranlukast a crucial role for R255 in binding of nucleotides to GPR17. Aromatic interactions are instead likely to play a predominant role in the recognition of pranlukast, suggesting that two different binding subsites are present on GPR17.

Dynamics of the internal water molecules in squid rhodopsin

Understanding the mechanism of G-protein coupled receptors action is of major interest for drug design. The visual rhodopsin is the prototype structure for the family A of G-protein coupled receptors. Upon photoisomerization of the covalently bound retinal chromophore, visual rhodopsins undergo a large-scale conformational change that prepares the receptor for a productive interaction with the G-protein. The mechanism by which the local perturbation of the retinal cis-trans isomerization is transmitted throughout the protein is not well understood. The crystal structure of the visual rhodopsin from squid solved recently suggests that a chain of water molecules extending from the retinal toward the cytoplasmic side of the protein may play a role in the signal transduction from the all-trans retinal geometry to the activated receptor. As a first step toward understanding the role of water in rhodopsin function, we performed a molecular dynamics simulation of squid rhodopsin embedded in a hydrated bilayer of polyunsaturated lipid molecules. The simulation indicates that the water molecules present in the crystal structure participate in favorable interactions with side chains in the interhelical region and form a persistent hydrogen-bond network in connecting Y315 to W274 via D80.

Calcium binding to the purple membrane: A molecular dynamics study

The purple membrane (PM) is a specialized membrane patch found in halophilic archaea, containing the photoreceptor bacteriorhodopsin (bR). It is long known that calcium ions bind to the PM, but their position and role remain elusive to date. Molecular dynamics simulations in conjunction with a highly detailed model of the PM have been used to investigate the stability of calcium ions placed at three proposed cation binding sites within bR, one near the Schiff base, one in the region of the proton release group, and one near Glu9. The simulations suggest that, of the sites investigated, the binding of calcium ions was most likely at the proton release group. Binding in the region of the Schiff base, while possible, was associated with significant changes in local geometry. Calcium ions placed near Glu9 in the interior of bR (simultaneously to a Ca(2+) near the Schiff base and another one near the Glu194-Glu204 site) were not stable. The results obtained are discussed in relation to recent experimental observations and theoretical considerations.

Potential proton-release channels in bacteriorhodopsin

The protein bacteriorhodopsin pumps protons across a bacterial membrane; its pumping cycle is triggered by the photoisomerization of a retinal cofactor and involves multiple proton-transfer reactions between intermittent protonation sites. These transfers are either direct or mediated by hydrogen-bonded networks, which may include internal water molecules. The terminal step of the proton-transfer sequence is the proton release from a pocket near Glu194 and Glu204 to the extracellular bulk during the transition from the L to the M photointermediate states. The polar and charged side chains connecting these two regions in the crystal structures show no structural changes between the initial bR state and the L/M states, and no intermittent protonation changes have been detected so far in this region. Based on biomolecular simulations, we propose two potential proton-release channels, which connect the release pocket to the extracellular medium. In simulations of the L photointermediate we observe bulk water entering these channels and forming transient hydrogen-bonded networks, which could serve as fast deprotonation pathways from the release pocket to the bulk via a Grotthuss mechanism. For the first channel, we find that the triple Arg7, Glu9, and Tyr79 acts as a valve, thereby gating water uptake and release. The second channel has two release paths, which split at the position Asn76/Pro77 underneath the release group. Here, water molecules either exchange directly with the bulk or diffuse within the protein towards Arg 134/Lys129, where the exchange with the bulk occurs.

Dynamics of voltage profile in enzymatic ion transporters, demonstrated in electrokinetics of proton pumping rhodopsin

H(+)-pumping rhodopsins mediate a primordial conversion of light to metabolic energy. Bacteriorhodopsin from Halobacterium salinarium is the first identified and (biochemically) best-studied H(+)-pumping rhodopsin. The electrical properties of H(+)-pumping rhodopsins, however, are known in more detail for the homolog Acetabularia rhodopsin, isolated from the eukaryotic green alga Acetabularia acetabulum. Based on data from Acetabularia rhodopsin we present a general reaction kinetic model of H(+)-pumping rhodopsins with only seven independent parameters, which fits the kinetic properties of photocurrents as functions of light, transmembrane voltage, internal and external pH, and time. The model describes fast photoisomerization of retinal with simultaneous H(+) transfer to an H(+) acceptor, reprotonation of retinal from the intracellular face via an H(+) donor, and proton release to the extracellular space via an H(+) release complex. The voltage sensitivities of the individual reaction steps and their temporal changes are treated here by a novel approach, whereby--as in an Ohmic voltage divider--the effective portions of the total transmembrane voltage decrease with the relative velocities of the individual reaction steps. This analysis quantitatively infers dynamic changes of the voltage profile and of the pK values of the H(+)-binding sites involved.

Dynamical heterogeneity of specific amino acids in bacteriorhodopsin

Components of biological macromolecules, complexes and membranes are animated by motions occurring over a wide range of time and length scales, the synergy of which is at the basis of biological activity. Understanding biological function thus requires a detailed analysis of the underlying dynamical heterogeneity. Neutron scattering, using specific isotope labeling, and molecular dynamics simulations were combined in order to study the dynamics of specific amino acid types in bacteriorhodopsin within the purple membrane (PM) of Halobacterium salinarum. Motions of leucine, isoleucine and tyrosine residues on the pico- to nanosecond time scale were examined separately as a function of temperature from 20 to 300 K. The dynamics of the three residue types displayed different temperature dependence: isoleucine residues have larger displacements compared to the global PM above 120 K; leucine residues have displacements similar to that of PM in the entire temperature range studied; and tyrosine residues have displacements smaller than that of the average membrane in an intermediate temperature range. Experimental features were mostly well reproduced by molecular dynamics simulations performed at five temperatures, which allowed the dynamical characterisation of the amino acids under study as a function of local environment. The resulting dynamical map of bacteriorhodopsin revealed that movements of a specific residue are determined by both its environment and its residue type.

Microbial rhodopsins: scaffolds for ion pumps, channels, and sensors

Microbial rhodopsins have been intensively researched for the last three decades. Since the discovery of bacteriorhodopsin, the scope of microbial rhodopsins has been considerably extended, not only in view of the large number of family members, but also their functional properties as pumps, sensors, and channels. In this review, we give a short overview of old and newly discovered microbial rhodopsins, the mechanism of signal transfer and ion transfer, and we discuss structural and mechanistic aspects of phototaxis.

Structures and spectral signatures of protonated water networks in bacteriorhodopsin

Networks of internal water molecules are thought to provide proton transfer pathways in many enzymatic and photosynthetic reactions. Extremely broad absorption continua observed in recent IR spectroscopic measurements on the photodriven proton pump bacteriorhodopsin (BR) suggest such networks may also serve as proton storage and release sites for these reactions. By combining electronic structure calculations with molecular mechanical force fields, we examine the dynamics and the resulting IR spectra of two protonated water networks, H+.(H2O)3 and H+.(H2O)4, in the release pocket of the initial state of BR, which possibly serve as proton donors to the extracellular surface. For both network sizes, topologically similar structures are found, which are anchored at residues E194 and E204 and stabilized by additional hydrogen bonds from neighboring protein side chains. These protonated water networks assume neither the classic Zundel nor Eigen motives but prefer wire-like topologies. Upon gauging calculated IR spectra of finite clusters with experimental gas-phase data, it is possible to link spectral features computed for these chain-like structures in the initial state of the BR photocycle to the measured absorption continua, in particular for the larger H+.(H2O)4 network. Furthermore, the free energy of proton dislocation along these chains is found to be within the range that is easily accessible at room temperature because of fluctuations.

Suppression of the back proton-transfer from Asp85 to the retinal Schiff base in bacteriorhodopsin: A theoretical analysis of structural elements

The transfer of a proton from the retinal Schiff base to the nearby Asp85 protein group is an essential step in the directional proton-pumping by bacteriorhodopsin. To avoid the wasteful back reprotonation of the Schiff base from Asp85, the protein must ensure that, following Schiff base deprotonation, the energy barrier for back proton-transfer from Asp85 to the Schiff base is larger than that for proton-transfer from the Schiff base to Asp85. Here, three structural elements that may contribute to suppressing the back proton-transfer from Asp85 to the Schiff base are investigated: (i) retinal twisting; (ii) hydrogen-bonding distances in the active site; and (iii) the number and location of internal water molecules. The impact of the pattern of bond twisting on the retinal deprotonation energy is dissected by performing an extensive set of quantum-mechanical calculations. Structural rearrangements in the active site, such as changes of the Thr89:Asp85 distance and relocation of water molecules hydrogen-bonding to the Asp85 acceptor group, may participate in the mechanism which ensures that following the transfer of the Schiff base proton to Asp85 the protein proceeds with the subsequent photocycle steps, and not with back proton transfer from Asp85 to the Schiff base.

External electric control of the proton pumping in bacteriorhodopsin

Comparative analysis of the photoelectric response of dried films of purple membranes (PM) depending on their degree of orientation is presented. Time dependence of the photo-induced protein electric response signal (PERS) of oriented and non-oriented films to a single laser pulse in the presence of the external electric field (EEF) was experimentally determined. The signal does not appear in the non-oriented films when the EEF is absent, whereas the PERS of the oriented PM films demonstrates the variable polarity on the microsecond time scale. In the presence of the EEF the PERS of the non-oriented film rises exponentially preserving the same polarization. The polarization of the PERS changes by changing the polarity of the EEF with no influence on the time constant of the PERS kinetics. The EEF effect on the PERS of the oriented films is more complicated. By subtracting the PERS when EEF not equal 0 from the PERS when EEF = 0 the resulting signal is comparable to that of the non-oriented films. Generalizing the experimental data we conclude that the EEF influence is of the same origin for the films of any orientation. To explain the experimental results the two-state model is suggested. It assumes that the EEF directionally changes the pK(a) values of the Schiff base (SB) and of the proton acceptor aspartic acid D85 in bacteriorhodopsin. Because of that the SB-->D85 proton transfer might be blocked and consequently the L-->M intermediate transition should vanish. Thus, on the characteristic time scale tau( L --> M ) approximately 30 micros; both intermediates, the M intermediate, appearing under normal conditions, and the L intermediate as persisting under the blocked conditions when D85 is protonated, should coexist in the film. The total PERS is a result of the potentials corresponding to the electrogenic products of intermediates L and M that are of the opposite polarity. It is concluded that the ratio of bacteriorhodopsin concentrations corresponding to the L and M intermediates is driven by the EEF and, consequently, it should define the PERS of the non-oriented films. According to this model the orientation degree of the film could be evaluated by describing the PERS.

Model for Proton Transport Coupled to Protein Conformational Change: Application to Proton Pumping in the Bacteriorhodopsin Photocycle

A modeling method is presented for protein systems in which proton transport is coupled to conformational change, as in proton pumps and in motors driven by the proton-motive force. Previously developed methods for calculating pKa values in proteins using a macroscopic dielectric model are extended beyond the equilibrium case to a master-equation model for the time evolution of the system through states defined by ionization microstate and a discrete set of conformers. The macroscopic dielectric model supplies free energy changes for changes of protonation microstate, while the method for obtaining the energetics of conformational change and the relaxation rates, the other ingredients needed for the master equation, are system dependent. The method is applied to the photoactivated proton pump, bacteriorhodopsin, using conformational free energy differences from experiment and treating relaxation rates through three adjustable parameters. The model is found to pump protons with an efficiency relatively insensitive to parameter choice over a wide range of parameter values, and most of the main features of the known photocycle from very early M to the return to the resting state are reproduced. The boundaries of these parameter ranges are such that short-range proton transfers are faster than longer-range ones, which in turn are faster than conformational changes. No relaxation rates depend on conformation. The results suggest that an "accessibility switch", while not ruled out, is not required and that vectorial proton transport can be achieved through the coupling of the energetics of ionization and conformational states.

pH dependence of light-driven proton pumping by an archaerhodopsin from Tibet: comparison with bacteriorhodopsin

The pH-dependence of photocycle of archaerhodopsin 4 (AR4) was examined, and the underlying proton pumping mechanism investigated. AR4 is a retinal-containing membrane protein isolated from a strain of halobacteria from a Tibetan salt lake. It acts as a light-driven proton pump like bacteriorhodopsin (BR). However, AR4 exhibits an "abnormal" feature--the time sequence of proton release and uptake is reversed at neutral pH. We show here that the temporal sequence of AR4 reversed to "normal"--proton release preceding proton uptake--when the pH is increased above 8.6. We estimated the pK(a) of the proton release complex (PRC) in the M-intermediate to be approximately 8.4, much higher than 5.7 of wide-type BR. The pH-dependence of the rate constant of M-formation shows that the pK(a) of PRC in the initial state of AR4 is approximately 10.4, whereas it is 9.7 in BR. Thus in AR4, the chromophore photoisomerization and subsequent proton transport from the Schiff base to Asp-85 is coupled to a decrease in the pK(a) of PRC from 10.4 to 8.4, which is 2 pK units less than in BR (4 units). This weakened coupling accounts for the lack of early proton release at neutral pH and the reversed time sequence of proton release and uptake in AR4. Nevertheless the PRC in AR4 effectively facilitates deprotonation of primary proton acceptor and recovery of initial state at neutral pH. We found also that all pK(a)s of the key amino acid residues in AR4 were elevated compared to those of BR.

The Cytomatrix as a Cooperative System of Macromolecular and Water Networks

Water was called by Szent-Gyorgi "life's mater and matrix, mother and medium." This chapter considers both aspects of his statement. Many astrobiologists argue that some, if not all, of Earth's water arrived during cometary bombardments. Amorphous water ices of comets possibly facilitated organization of complex organic molecules, kick-starting prebiotic evolution. In Gaian theory, Earth retains its water as a consequence of biological activity. The cell cytomatrix is a proteinaceous matrix/lattice incorporating the cytoskeleton, a pervasive, holistic superstructural network that integrates metabolic pathways. Enzymes of metabolic pathways are ordered in supramolecular clusters (metabolons) associated with cytoskeleton and/or membranes. Metabolic intermediates are microchanneled through metabolons without entering a bulk aqueous phase. Rather than being free in solution, even major signaling ions are probably clustered in association with the cytomatrix. Chloroplasts and mitochondria, like bacteria and archaea, also contain a cytoskeletal lattice, metabolons, and channel metabolites. Eukaryotic metabolism is mathematically a scale-free or small-world network. Enzyme clusters of bacterial origin are incorporated at a pathway level that is architecturally archaean. The eucaryotic cell may be a product of serial endosymbiosis, a chimera. Cell cytoplasm is approximately 80% water. Water is indisputably a conserved structural element of proteins, essential to their folding, specificity, ligand binding, and to enzyme catalysis. The vast literature of organized cell water has long argued that the cytomatrix and cell water are an entire system, a continuum, or gestalt. Alternatives are offered to mainstream explanations of cell electric potentials, ion channel, enzyme, and motor protein function, in terms of high-order cooperative systems of ions, water, and macromolecules. This chapter describes some prominent concepts of organized cell water, including vicinal water network theory, the association-induction hypothesis, wave-cluster theory, phase-gel transition theories, and theories of low- and high-density water polymorphs.

Molecular dynamics simulations of proteins in lipid bilayers

With recent advances in X-ray crystallography of membrane proteins promising many new high-resolution structures, molecular dynamics simulations will become increasingly valuable for understanding membrane protein function, as they can reveal the dynamic behavior concealed in the static structures. Dramatic increases in computational power, in synergy with more efficient computational methodologies, now allow us to carry out molecular dynamics simulations of any structurally known membrane protein in its native environment, covering timescales of up to 0.1 micros. At the frontiers of membrane protein simulations are ion channels, aquaporins, passive and active transporters, and bioenergetic proteins.

Application of classical molecular dynamics for evaluation of proton transfer mechanism on a protein

Proton transfer reactions on surfaces are prevalent in biology, chemistry and physics. In the present study, we employed classical Molecular Dynamics simulations to search for the presence of transient configurations that enable proton transfer, or proton sharing, between adjacent carboxylate groups on the protein surface. The results demonstrate that, during random fluctuations of the residues on the surface, there are repeated situations in which nearby carboxylates either share a common proton through a hydrogen bond, or are connected by a few water molecules that form conducting networks. These networks do not extend out of the common Coulomb cage of the participating residues and the lifetimes of the bridged structures are sufficiently long to allow passage of a proton between the carboxylates. The detection of domains capable of supporting a rapid proton transfer on a protein supports the notion that clusters of carboxylates are the operative elements of proton collecting antennae, as in bacteriorhodopsin, cytochrome c oxidase or the photosynthetic reaction center.

Functional waters in intraprotein proton transfer monitored by FTIR difference spectroscopy

Much progress has been made in our understanding of water molecule reactions on surfaces, proton solvation in gas-phase water clusters and proton transfer through liquids. Compared with our advanced understanding of these physico-chemical systems, much less is known about individual water molecules and their cooperative behaviour in heterogeneous proteins during enzymatic reactions. Here we use time-resolved Fourier transform infrared spectroscopy (trFTIR) and in situ H2(18)O/H2(16)O exchange FTIR to determine how the membrane protein bacteriorhodopsin uses the interplay among strongly hydrogen-bonded water molecules, a water molecule with a dangling hydroxyl group and a protonated water cluster to transfer protons. The precise arrangement of water molecules in the protein matrix results in a controlled Grotthuss proton transfer, in contrast to the random proton migration that occurs in liquid water. Our findings support the emerging paradigm that intraprotein water molecules are as essential for biological functions as amino acids.