Functional waters in intraprotein proton transfer monitored by FTIR difference spectroscopy

Modeling the membrane potential generation of bacteriorhodopsin

The archaeon Halobacterium salinarum can grow phototrophically with only light as its energy source. It uses the retinal containing and light-driven proton pump bacteriorhodopsin to enhance the membrane potential which drives the ATP synthase. Therefore, a model of the membrane potential generation of bacteriorhodopsin is of central importance to the development of a mathematical model of the bioenergetics of H. salinarum. To measure the current produced by bacteriorhodopsin at different light intensities and clamped voltages, we expressed the gene in Xenopus laevis oocytes. We present current-voltage measurements and a mathematical model of the current-voltage relationship of bacteriorhodopsin and its generation of the membrane potential. The model consists of three intermediate states, the BR, L, and M states, and comparisons between model predictions and experimental data show that the L to M reaction must be inhibited by the membrane potential. The model is not able to fit the current-voltage measurements when only the M to BR phase is membrane potential dependent, while it is able to do so when either only the L to M reaction or both reactions (L to M and M to BR) are membrane potential dependent. We also show that a decay term is necessary for modeling the rate of change of the membrane potential.

Controlled ionic condensation at the surface of a native extremophile membrane

At the nanoscale level biological membranes present a complex interface with the solvent. The functional dynamics and relative flexibility of membrane components together with the presence of specific ionic effects can combine to create exciting new phenomena that challenge traditional theories such as the Derjaguin-Landau-Verwey-Overbeek (DLVO) theory or models interpreting the role of ions in terms of their ability to structure water (structure making/breaking). Here we investigate ionic effects at the surface of a highly charged extremophile membrane composed of a proton pump (bacteriorhodopsin) and archaeal lipids naturally assembled into a 2D crystal. Using amplitude-modulation atomic force microscopy (AM-AFM) in solution, we obtained sub-molecular resolution images of ion-induced surface restructuring of the membrane. We demonstrate the presence of a stiff cationic layer condensed at its extracellular surface. This layer cannot be explained by traditional continuum theories. Dynamic force spectroscopy experiments suggest that it is produced by electrostatic correlation mediated by a Manning-type condensation of ions. In contrast, the cytoplasmic surface is dominated by short-range repulsive hydration forces. These findings are relevant to archaeal bioenergetics and halophilic adaptation. Importantly, they present experimental evidence of a natural system that locally controls its interactions with the surrounding medium and challenges our current understanding of biological interfaces.

A residue outside the active site CXXC motif regulates the catalytic efficiency of Glutaredoxin 3

The glutaredoxin (Grx) family of oxidoreductases has a conserved residue at position 8 that varies between Arginine in Grx1 and Lysine in Grx3. It has been proposed that this Arg/Lys change is the main cause for the 35 mV difference in redox potential between the two enzymes. To gain insights into the catalytic machinery of Grx3 and directly evaluate the role of residue 8 in the catalysis of thiol-disulfide exchange by this enzyme, we synthesized the "wild type" enzyme (sGrx3), and four analogues substituting the lysine at position 8 with arginine, ornithine (Orn), citrulline (Cit) and norvaline (Nva). The redox potential and equilibration kinetics with thioredoxin (Trx1) were determined for each enzyme by fluorescence intensity. While minor effects on redox potential were observed, we found that residue 8 had a more marked effect on the catalytic efficiency of this enzyme. Surprisingly, truncation of the functional group resulted in a more efficient enzyme, Lys8Nva, exhibiting rate constants that are an order of magnitude higher than sGrx3 for both forward and reverse reactions. These observations pose the question why would a residue that reduces the rate of enzyme turnover be evolutionarily conserved? The significant changes in the kinetic parameters suggest that this position plays an important role in the thiol-disulfide exchange reaction by affecting the nucleophilic thiolate through electrostatic or hydrogen bonding interactions. Since the reduced Grx has an exposed thiol that could easily be alkylated, either Arg or Lys could act as a gatekeeper that deters unwanted electrophiles from attacking the active site thiolate.

Protein hydration dynamics and molecular mechanism of coupled water-protein fluctuations

Protein surface hydration is fundamental to its structural stability and flexibility, and water-protein fluctuations are essential to biological function. Here, we report a systematic global mapping of water motions in the hydration layer around a model protein of apomyoglobin in both native and molten globule states. With site-directed mutagenesis, we use intrinsic tryptophan as a local optical probe to scan the protein surface one at a time with single-site specificity. With femtosecond resolution, we examined 16 mutants in two states and observed two types of water-network relaxation with distinct energy and time distributions. The first water motion results from the local collective hydrogen-bond network relaxation and occurs in a few picoseconds. The initial hindered motions, observed in bulk water in femtoseconds, are highly suppressed and drastically slow down due to structured water-network collectivity in the layer. The second water-network relaxation unambiguously results from the lateral cooperative rearrangements in the inner hydration shell and occurs in tens to hundreds of picoseconds. Significantly, this longtime dynamics is the coupled interfacial water-protein motions and is the direct measurement of such cooperative fluctuations. These local protein motions, although highly constrained, are necessary to assist the longtime water-network relaxation. A series of correlations of hydrating water dynamics and coupled fluctuations with local protein's chemical and structural properties were observed. These results are significant and reveal various water behaviors in the hydration layer with wide heterogeneity. We defined a solvation speed and an angular speed to quantify the water-network rigidity and local protein flexibility, respectively. We also observed that the dynamic hydration layer extends to more than 10 A. Finally, from native to molten globule states, the hydration water networks loosen up, and the protein locally becomes more flexible with larger global plasticity and partial unfolding.

A G protein-coupled receptor at work: the rhodopsin model

G protein-coupled receptors (GPCRs) are ubiquitous signal transducers in cell membranes, as well as important drug targets. Interaction with extracellular agonists turns the seven transmembrane helix (7TM) scaffold of a GPCR into a catalyst for GDP and GTP exchange in heterotrimeric Galphabetagamma proteins. Activation of the model GPCR, rhodopsin, is triggered by photoisomerization of its retinal ligand. From the augmentation of biochemical and biophysical studies by recent high-resolution 3D structures, its activation intermediates can now be interpreted as the stepwise engagement of protein domains. Rearrangement of TM5-TM6 opens a crevice at the cytoplasmic side of the receptor into which the C terminus of the Galpha subunit can bind. The Galpha C-terminal helix is used as a transmission rod to the nucleotide binding site. The mechanism relies on dynamic interactions between conserved residues and could therefore be common to other GPCRs.

Probability distributions of hydration water molecules around polar protein atoms obtained by a database analysis

Hydration structures on protein surfaces are visualized by high-resolution cryogenic X-ray crystallography. We calculated the probability distributions of 4,831,570 hydration water molecules found around the 4,214,227 polar atoms in main chains and hydrophilic side chains from the 17,984 crystal structures in the Protein Data Bank. The structures are refined using the diffraction data collected below 150 K and at resolutions of better than 2.2 A. The calculated distributions were nonrandom but condensed into a few clusters. The clusters were decomposed into the distance and angular distributions by viewing from the polar coordinate system. The major peaks in the clusters were almost located along the directions of the N-H and O-H bonds or the lone pairs of oxygen atoms. The Gaussian fitting method was applied for the distribution profiles to evaluate quantitatively the peak positions and the widths. The parameters characterizing the distributions apparently depended on the hydrogen-bond partners of water molecules and on the modes whether the water molecules acted as donors or acceptors of protons. This led to propose the different roles of NH(n) (n = 1, 3), OH, and CO groups in protein hydration and possible in protein-ligand and protein-protein interaction: While C horizontal lineO groups appear to control the H-bond distances, NH(n) groups likely limit the angular range of H-bonds. The OH groups have both characteristics. In addition, it was also demonstrated that polar protein atoms were arranged to satisfy the tetrahedral hydrogen-bond geometry of water molecules, suggesting essential roles of water molecules in the folding process and in the stabilization of protein structures. These probability distributions are probably one of fundamental data to better understand the roles of hydration water molecules in the folding process and the stability of proteins in solution.

Conformational heterogeneity and propagation of structural changes in the LOV2/Jalpha domain from Avena sativa phototropin 1 as recorded by temperature-dependent FTIR spectroscopy

Phototropins control phototropism, chloroplast movement, stomatal opening, and leaf expansion in plants. Phototropin 1 (phot1) is composed of a kinase domain linked to two blue light-sensing domains, LOV2 and LOV1, which bind flavin mononucleotide. Disruption of the interaction between the LOV2 domain and a helical segment named Jalpha, joining LOV to the kinase domain, induces the subsequent kinase activity of phototropin 1 and further-downstream signal transduction. Here we study the effects of temperature and hydration on the light-triggered signal propagation in the phot1 LOV2 domain of Avena sativa (AsLOV2/Jalpha), using Fourier transform infrared spectroscopy to unravel part of the molecular mechanism of phototropin 1. We report that AsLOV2/Jalpha shows an intense signal in the amide I and II regions, arising mainly from beta-sheet changes and the unbinding of the Jalpha helix from the Per-ARNT-Sim core and its subsequent partial unfolding. Importantly, these structural changes only occur under conditions of full hydration and at temperatures above 280 K. We characterized a newly isolated low-hydration intermediate that shows a downshift of high-frequency amide I signals and that possibly corresponds to loop tightening, without large beta-sheet or Jalpha structural changes. In addition, we report a heterogeneity in AsLOV2/Jalpha involving two different C(4)=O conformer populations, coexisting in the dark state and characterized by C(4)=O carbonyl frequencies at 1712 cm(-1) and 1694 cm(-1) that are attributable to a single H-bond and two H-bonds at this site, respectively. Such conformers display slightly shifted absorption spectra and cause a splitting of the 475-nm band in the ultraviolet/visible spectra of LOV domains at low temperature.

Residues in the H+ translocation site define the pKa for sugar binding to LacY

A remarkably high pKa of approximately 10.5 has been determined for sugar-binding affinity to the lactose permease of Escherichia coli (LacY), indicating that, under physiological conditions, substrate binds to fully protonated LacY. We have now systematically tested site-directed replacements for the residues involved in sugar binding, as well as H+ translocation and coupling, in order to determine which residues may be responsible for this alkaline pKa. Mutations in the sugar-binding site (Glu126, Trp151, Glu269) markedly decrease affinity for sugar but do not alter the pKa for binding. In contrast, replacements for residues involved in H+ translocation (Arg302, Tyr236, His322, Asp240, Glu325, Lys319) exhibit pKa values for sugar binding that are either shifted toward neutral pH or independent of pH. Values for the apparent dissociation constant for sugar binding (K(d)(app)) increase greatly for all mutants except neutral replacements for Glu325 or Lys319, which are characterized by remarkably high affinity sugar binding (i.e., low K(d)(app)) from pH 5.5 to pH 11. The pH dependence of the on- and off-rate constants for sugar binding measured directly by stopped-flow fluorometry implicates k(off) as a major factor for the affinity change at alkaline pH and confirms the effects of pH on K(d)(app) inferred from steady-state fluorometry. These results indicate that the high pKa for sugar binding by wild-type LacY cannot be ascribed to any single amino acid residue but appears to reside within a complex of residues involved in H+ translocation. There is structural evidence for water bound in this complex, and the water could be the site of protonation responsible for the pH dependence of sugar binding.

Understanding dynamics of myoglobin in heterogeneous aqueous environments using coupled water fractions

This work presents an analysis of near environment of myoglobin (Mb) in different aqueous solutions (in the presence of NaCl, sucrose, trehalose, urea, and glycerol) using the coupled water fractions measured using a quartz crystal microbalance (QCM). The secondary structural features of the protein from circular dichroic (CD) spectroscopy and the coupled water fractions give important clues to the overall dynamics of the protein. Using time resolved fluorescence, these leads have been applied to understand the observed lifetime relaxations of Mb. Though the time scales of observation of coupled water and the lifetimes are very different, our study suggests that the trends in coupled water fraction seem to be good indicators for regulation of the relaxation dynamics of the protein. The relaxations generally show a triphasic distribution of time scales. The initial relaxation in the picoseconds time scale represents the local motions of coupled water followed by a slightly slower decay in hundreds of picoseconds attributable to coupled water-'quasi free' water interactions. The third nanosecond lifetime is due to changes in transitions in isomers of hydrated protein. The dynamics of coupled water in Mb with NaCl is the fastest (around 21 ps) and is slowest in glycerol (250 ps). The results strongly indicate that it is the resident times of water molecules that play a dominant role in the overall stability of protein in a particular hydrated isomer and not just always the number of such water molecules in the hydrated protein.

The cytochrome ba3 oxygen reductase from Thermus thermophilus uses a single input channel for proton delivery to the active site and for proton pumping

The heme-copper oxygen reductases are redox-driven proton pumps that generate a proton motive force in both prokaryotes and mitochondria. These enzymes have been divided into 3 evolutionarily related groups: the A-, B- and C-families. Most experimental work on proton-pumping mechanisms has been performed with members of the A-family. These enzymes require 2 proton input pathways (D- and K-channels) to transfer protons used for oxygen reduction chemistry and for proton pumping, with the D-channel transporting all pumped protons. In this work we use site-directed mutagenesis to demonstrate that the ba(3) oxygen reductase from Thermus thermophilus, a representative of the B-family, does not contain a D-channel. Rather, it utilizes only 1 proton input channel, analogous to that of the A-family K-channel, and it delivers protons to the active site for both O2 chemistry and proton pumping. Comparison of available subunit I sequences reveals that the only structural elements conserved within the oxygen reductase families that could perform these functions are active-site components, namely the covalently linked histidine-tyrosine, the Cu(B) and its ligands, and the active-site heme and its ligands. Therefore, our data suggest that all oxygen reductases perform the same chemical reactions for oxygen reduction and comprise the essential elements of the proton-pumping mechanism (e.g., the proton-loading and kinetic-gating sites). These sites, however, cannot be located within the D-channel. These results along with structural considerations point to the A-propionate region of the active-site heme and surrounding water molecules as the proton-loading site.

Structural waters define a functional channel mediating activation of the GPCR, rhodopsin

Structural water molecules may act as prosthetic groups indispensable for proper protein function. In the case of allosteric activation of G protein-coupled receptors (GPCRs), water likely imparts structural plasticity required for agonist-induced signal transmission. Inspection of structures of GPCR superfamily members reveals the presence of conserved embedded water molecules likely important to GPCR function. Coupling radiolytic hydroxyl radical labeling with rapid H(2)O(18) solvent mixing, we observed no exchange of these structural waters with bulk solvent in either ground state or for the Meta II or opsin states. However, the radiolysis approach permitted labeling of selected side chain residues within the transmembrane helices and revealed activation-induced changes in local structural constraints likely mediated by dynamics of both water and protein. These results suggest both a possible general mechanism for water-dependent communication in family A GPCRs based on structural conservation, and a strategy for probing membrane protein structure.

Site-specific difference 2D-IR spectroscopy of bacteriorhodopsin

We demonstrate the extension of the principle of difference Fourier transform infrared (FTIR) spectroscopy to difference 2D-IR spectroscopy. To this end, we measure difference 2D-IR spectra of the protein bacteriorhodopsin in its early J- and K-intermediates. By comparing with the static 2D-IR spectrum of the protonated Schiff base of all-trans retinal, we demonstrate that the 2D-IR spectrum of the all-trans retinal chromophore in bacteriorhodopsin can be measured with the background from the remainder of the protein completely suppressed. We discuss several models to interpret the detailed line shape of the difference 2D-IR spectrum.

Efforts toward the direct experimental characterization of enzyme microenvironments: tyrosine100 in dihydrofolate reductase

State secrets: Site-specific deuteration and FTIR studies reveal that Tyr100 in dihydrofolate reductase plays an important role in catalysis, with a strong electrostatic coupling occurring between Tyr100 and the charge that develops in the hydride-transfer transition state (see picture, NADP(+) purple, Tyr100 green). However, relaying correlated motions that facilitate catalysis from distal sites of the protein to the hydride donor may also be involved.

Long-distance proton transfer with a break in the bacteriorhodopsin active site

Bacteriorhodopsin is a proton-pumping membrane protein found in the plasma membrane of the archaeon Halobacterium salinarium. Light-induced isomerization of the retinal chromophore from all-trans to 13-cis leads to a sequence of five conformation-coupled proton transfer steps and the net transport of one proton from the cytoplasmic to the extracellular side of the membrane. The mechanism of the long-distance proton transfer from the primary acceptor Asp85 to the extracellular proton release group during the O --> bR is poorly understood. Experiments suggest that this long-distance transfer could involve a transient state [O] in which the proton resides on the intermediate carrier Asp212. To assess whether the transient protonation of Asp212 participates in the deprotonation of Asp85, we performed hybrid Quantum Mechanics/Molecular Mechanics proton transfer calculations using different protein structures and with different retinal geometries and active site water molecules. The structural models were assessed by computing UV-vis excitation energies and C=O vibrational frequencies. The results indicate that a transient [O] conformer with protonated Asp212 could indeed be sampled during the long-distance proton transfer to the proton release group. Our calculations suggest that, in the starting proton transfer state O, the retinal is strongly twisted and at least three water molecules are present in the active site.

Conformations of NhaA, the Na+/H+ exchanger from Escherichia coli, in the pH-activated and ion-translocating states

NhaA, the main sodium-proton exchanger in the inner membrane of Escherichia coli, regulates the cytosolic concentrations of H+ and Na+. It is inactive at acidic pH, becomes active between pH 6 and pH 7, and reaches maximum activity at pH 8. By cryo-electron microscopy of two-dimensional crystals grown at pH 4 and incubated at higher pH, we identified two sequential conformational changes in the protein in response to pH or substrate ions. The first change is induced by a rise in pH from 6 to 7 and marks the transition from the inactive state to the pH-activated state. pH activation, which precedes the ion-induced conformational change, is accompanied by an overall expansion of the NhaA monomer and a local ordering of the N-terminus. The second conformational change is induced by the substrate ions Na+ and Li+ at pH above 7 and involves a 7-A displacement of helix IVp. This movement would cause a charge imbalance at the ion-binding site that may trigger the release of the substrate ion and open a periplasmic exit channel.

Active water in protein-protein communication within the membrane: the case of SRII-HtrII signal relay

We detect internal water molecules in a membrane-embedded receptor-transducer complex and demonstrate water structure changes during formation of the signaling state. Time-resolved FTIR spectroscopy reveals stimulus-induced repositioning of one or more structurally active water molecules to a significantly more hydrophobic environment in the signaling state of the sensory rhodopsin II (SRII)-transducer (HtrII) complex. These waters, distinct from bound water molecules within the SRII receptor, appear to be in the middle of the transmembrane interface region near the Tyr199(SRII)-Asn74(HtrII) hydrogen bond. We conclude that water potentially plays an important role in the SRII --> HtrII signal transfer mechanism in the membrane's hydrophobic core.

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.

How does a membrane protein achieve a vectorial proton transfer via water molecules?

We present a detailed mechanism for the proton transfer from a protein-bound protonated water cluster to the bulk water directed by protein side chains in the membrane protein bacteriorhodopsin. We use a combined approach of time-resolved Fourier transform infrared spectroscopy, molecular dynamics simulations, and X-ray structure analysis to elucidate the functional role of a hydrogen bond between Ser193 and Glu204. These two residues seal the internal protonated water cluster from the bulk water and the protein surface. During the photocycle of bacteriorhodopsin, a transient protonation of Glu204 leads to a breaking of this hydrogen bond. This breaking opens the gate to the extracellular bulk water, leading to a subsequent proton release from the protonated water cluster. We show in detail how the protein achieves vectorial proton transfer via protonated water clusters in contrast to random proton transfer in liquid water.

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.

Water-mediated spectral shifts in rhodopsin and bathorhodopsin

The role of water molecules in spectral tuning of proteins has been left largely unexplored. This topic is important because changing hydrogen bond patterns during the activation process may lead to spectral shifts which can be of diagnostic value for the underlying structures. Arguments put forward in this article are based on spectral shift calculations of the rhodopsin and bathorhodopsin chromophore due to wat2a and 2b in the presence and absence of the counterion and of the amino acids lining the rhodopsin binding pocket. They show, among others, that a single water molecule can shift the absorbance by up to 0.1 eV or 34 nm depending on the environment of the chromophore.