Structural clues to the mechanism of ion pumping in bacteriorhodopsin

Probing Structure and Reaction Dynamics of Proteins Using Time-Resolved Resonance Raman Spectroscopy

The mechanistic understanding of protein functions requires insight into the structural and reaction dynamics. To elucidate these processes, a variety of experimental approaches are employed. Among them, time-resolved (TR) resonance Raman (RR) is a particularly versatile tool to probe processes of proteins harboring cofactors with electronic transitions in the visible range, such as retinal or heme proteins. TR RR spectroscopy offers the advantage of simultaneously providing molecular structure and kinetic information. The various TR RR spectroscopic methods can cover a wide dynamic range down to the femtosecond time regime and have been employed in monitoring photoinduced reaction cascades, ligand binding and dissociation, electron transfer, enzymatic reactions, and protein un- and refolding. In this account, we review the achievements of TR RR spectroscopy of nearly 50 years of research in this field, which also illustrates how the role of TR RR spectroscopy in molecular life science has changed from the beginning until now. We outline the various methodological approaches and developments and point out current limitations and potential perspectives.

Hydronium Ions Accompanying Buried Acidic Residues Lead to High Apparent Dielectric Constants in the Interior of Proteins

Internal ionizable groups are known to play important roles in protein functions. A mystery that has attracted decades of extensive experimental and theoretical studies is the apparent dielectric constants experienced by buried ionizable groups, which are much higher than values expected for protein interiors. Many interpretations have been proposed, such as water penetration, conformational relaxation, local unfolding, protein intrinsic backbone fluctuations, etc. However, these interpretations conflict with many experimental observations. The virtual mixture of multiple states (VMMS) simulation method developed in our lab provides a direct approach for studying the equilibrium of multiple chemical states and can monitor p K_a values along simulation trajectories. Through VMMS simulations of staphylococcal nuclease (SNase) variants with internal Asp or Glu residues, we discovered that cations were attracted to buried deprotonated acidic groups and the presence of the nearby cations were essential to reproduce experimentally measured p K_a values. This finding, combined with structural analysis and validation simulations, suggests that the proton released from a deprotonation process stays near the deprotonated group inside proteins, possibly in the form of a hydronium ion. The existence of a proton near a buried charge has many implications in our understanding of protein functions.

A Coupled Ionization-Conformational Equilibrium Is Required To Understand the Properties of Ionizable Residues in the Hydrophobic Interior of Staphylococcal Nuclease

Ionizable residues in the interior of proteins play essential roles, especially in biological energy transduction, but are relatively rare and seem incompatible with the complex and polar environment. We perform a comprehensive study of the internal ionizable residues on 21 variants of staphylococcal nuclease with internal Lys, Glu, or Asp residues. Using pH replica exchange molecular dynamics simulations, we find that, in most cases, the pK_a values of these internal ionizable residues are shifted significantly from their values in solution. Our calculated results are in excellent agreement with the experimental observations of the Garcia-Moreno group. We show that the interpretation of the experimental pK_a values requires the study of not only protonation changes but also conformational changes. The coupling between the protonation and conformational equilibria suggests a mechanism for efficient pH-sensing and regulation in proteins. This study provides new physical insights into how internal ionizable residues behave in the hydrophobic interior of proteins.

Probing Small Protonated Water Clusters in Acetonitrile Solutions by 1H NMR

In a previous publication by Kalish et al. (J. Phys. Chem. A 115 (2011) 4063) the existence of well defined small protonated water clusters in acetonitrile has been established by IR spectroscopy. Here we report on a 1H NMR study of triflic acid, CF3SO3H, in acetonitrile-water solutions. Using NMR we are able to corroborate the general solvation scheme we have proposed for the hydrated proton in acetonitrile as a function of the molar ratio between the strong mineral acid and water, n = [H2O]/[acid]. According to this scheme, backed now by both IR absorption spectroscopy and NMR measurements, the very strong triflic acid completely dissociates in acetonitrile/water solutions to yield the aqueous proton and the triflate anion when n > 1. Furthermore, increasing n results in the proton solvated in increasingly larger water clusters formed within the acetonitrile solution.

Clearly distinguishable by NMR are the smallest protonated water clusters, the protonated water monomer, H3 +O, and the protonated water dimer, H5 +O2, which dominate the solution for n = 1,2,3. For larger n the NMR study indicates the gradual increase of the average protonated water cluster size as a function of n while the proton inner solvation core more closely retaining the characteristics of a deformed protonated water dimer, (H2O-H+⋯OH2) s than that of the protonated water monomer (H3 +O) s .

The high dielectric constant of staphylococcal nuclease is encoded in its structural architecture

The pK(a) values of Lys-66, Glu-66, and Asp-66 buried in the interior of the staphylococcal nuclease Δ+PHS variant were reported to be shifted by as much as 5 pK(a) units from their normal values. Reproducing the pK(a) of these buried ionizable residues using continuum electrostatic calculations required the use of a high protein dielectric constant of 10 or higher. The apparent high dielectric constant has been rationalized as a consequence of a local structural reorganization or increased fluctuations in the microenvironment of the mutation site (Chimenti, M. S., et al. J. Mol. Biol. 2011, 405, 361-377). We have calculated the dielectric constant of Δ+PHS and the Lys-66, Asp-66, and Glu-66 mutants from first principles using the Kirkwood-Fröhlich equation and discovered that staphylococcal nuclease has a naturally high dielectric constant ranging from 20 to 30. This high dielectric constant does not change significantly with the mutation of residue 66 or with the ionization of the mutated residues. Calculation of the spatial dependence of the dielectric constant for the microenvironment of residue-66 produces a value of about 10, which matches well with the apparent dielectric constant needed to reproduce pK(a) values from continuum electrostatic calculations. Our results suggest an alternative explanation that the high dielectric constant of staphylococcal nuclease is a property resulting from the intrinsic backbone fluctuations originating from its structural architecture.

Apparent stoichiometry of water in proton hydration and proton dehydration reactions in CH3CN/H2O solutions

Gradual solvation of protons by water is observed in liquids by mixing strong mineral acids with various amounts of water in acetonitrile solutions, a process which promotes rapid dissociation of the acids in these solutions. The stoichiometry of the reaction XH(+) + n(H(2)O) = X + (H(2)O)(n)H(+) was studied for strong mineral acids (negatively charged X, X = ClO(4)¯, Cl¯, Br¯, I¯, CF(3)SO(3)¯) and for strong cationic acids (uncharged X, X = R*NH(2), H(2)O). We have found by direct quantitative analysis preference of n = 2 over n = 1 for both groups of proton transfer reactions at relatively low water concentrations in acetonitrile. At high water concentrations, we have found that larger water solvates must also be involved in the solvation of the proton while the spectral features already observed for n = 2, H(+)(H(2)O)(2), remain almost unchanged at large n values up to at least 10 M of water.

The effects of pK(a) tuning on the thermodynamics and kinetics of folding: design of a solvent-shielded carboxylate pair at the a-position of a coiled-coil

The tuning of the pK(a) of ionizable residues plays a critical role in various protein functions, such as ligand-binding, catalysis, and allostery. Proteins harness the free energy of folding to position ionizable groups in highly specific environments that strongly affect their pK(a) values. To investigate the interplay among protein folding kinetics, thermodynamics, and pK(a) modulation, we introduced a pair of Asp residues at neighboring interior positions of a coiled-coil. A single Asp residue was replaced for an Asn side chain at the a-position of the coiled-coil from GCN4, which was also crosslinked at the C-terminus via a flexible disulfide bond. The thermodynamic and kinetic stability of the system was measured by circular dichroism and stopped-flow fluorescence as a function of pH and concentration of guanidine HCl. Both sets of data are consistent with a two-state equilibrium between fully folded and unfolded forms. Distinct pK(a) values of 6.3 and 5.35 are assigned to the first and second protonation of the Asp pair; together they represent an energetic difference of 5 kcal/mol relative to the protonation of two Asp residues with unperturbed pK(a) values. Analysis of the rate data as a function of pH and denaturant concentration allowed calculation of the kinetic constants for the conformational transitions of the peptide with the Asp residues in the doubly protonated, singly protonated, and unprotonated forms. The doubly and singly protonated forms fold rapidly, and a ϕ-value analysis shows that their contribution to folding occurs subsequent to the transition state ensemble for folding. By contrast, the doubly charged state shows a reduced rate of folding and a ϕ-value near 0.5 indicative of a repulsive interaction, and possibly also heterogeneity in the transition state ensemble.

Conformational consequences of ionization of Lys, Asp, and Glu buried at position 66 in staphylococcal nuclease

The pK(a) values measured previously for the internal Lys-66, Asp-66, and Glu-66 in variants of a highly stable form of staphylococcal nuclease are shifted by as many as 5 pK(a) units relative to normal pK(a) values in water. These shifts cannot be reproduced with continuum electrostatics calculations with static structures unless the protein is treated with high dielectric constants near 10. These high apparent dielectric constants are inconsistent with the highly hydrophobic microenvironments of the ionizable moieties in crystal structures. To examine the origins of these high apparent dielectric constants, we showed that the pK(a) values of these internal residues are sensitive to the global stability of the protein; the shifts tend to be smaller in less stable forms of nuclease. This implies that the apparent dielectric constants reported by these internal ionizable groups are high because they reflect conformational reorganization coupled to their ionization. To detect this directly, acid-base titrations monitored with Trp fluorescence and near-UV and far-UV CD spectroscopy were performed on variants with Lys-66, Glu-66, or Asp-66 in background proteins with different stabilities. Conformational reorganization coupled to the ionization of the internal groups was spectroscopically detectable, especially in the less stable background proteins. The data show that to improve the accuracy of structure-based pK(a) calculations of internal groups the calculations will have to treat explicitly all structural reorganization coupled to ionization. The data also suggest a novel approach to mapping the folding free energy landscape of proteins by using internal ionizable groups to stabilize partially unfolded states.

Microscopic pKa analysis of Glu286 in cytochrome c oxidase (Rhodobacter sphaeroides): toward a calibrated molecular model

As stringent tests for the molecular model and computational protocol, microscopic pK(a) calculations are performed for the key residue, Glu286, in cytochrome c oxidase (CcO) using a combined quantum mechanical/molecular mechanical (QM/MM) potential and a thermodynamic integration protocol. The impact of the number of water molecules in the hydrophobic cavity and protonation state of several key residues (e.g., His334, Cu(B)-bound water, and PRD(a3)) on the computed microscopic pK(a) values of Glu286 has been systematically examined. To help evaluate the systematic errors in the QM/MM-based protocol, microscopic pK(a) calculations have also been carried out for sites in a soluble protein (Asp70 in T4 lysozyme) and a better-characterized membrane protein (Asp85 in bacteriorhodopsin). Overall, the results show a significant degree of internal consistency and reproducibility that support the effectiveness of the computational framework. Although the number of water molecules in the hydrophobic cavity does not greatly influence the computed pK(a) of Glu286, the protonation states of several residues, some of which are rather far away, have more significant impacts. Adopting the standard protonation state for all titratable residues leaves a large net charge on the system and a significantly elevated pK(a) for Glu286, highlighting that any attempt to address the energetics of proton transfers in CcO at a microscopic level should carefully select the protonation state of residues, even those not in the immediate neighborhood of the active site. The calculations indirectly argue against the deprotonation of His334 for the proton pumping process, although further studies that explicitly compute its pK(a) are required for a more conclusive statement. Finally, the deprotonated Glu286 is found to be in a stable water-mediated connection with PRD(a3) for at least several nanoseconds when this presumed pumping site is protonated. This does not support the proposed role of Glu286 as a robust gating valve that prevents proton leakage, although a conclusive statement awaits a more elaborate characterization of the Glu286-PRD(a3) connectivity with free energy simulations and a protonated PRD(a3). The large sets of microscopic simulations performed here have provided useful guidance to the establishment of a meaningful molecular model and effective computational protocol for explicitly analyzing the proton transfer kinetics in CcO, which is required for answering key questions regarding the pumping function of this fascinating and complex system.

X-ray diffraction of bacteriorhodopsin photocycle intermediates (Review)

Recent advances in the crystallography of bacteriorhodopsin, the light-driven proton pump, have yielded structural models for all intermediates of the photochemical cycle. For seven of the species, X-ray diffraction data were collected from trapped photostationary states in crystals, and for the two remaining ones the structures of selected mutants are available. The changes of the retinal chromophore, protein and bound water describe, at an atomic level, how accommodation of the twisted photoisomerized retinal to its binding site causes de-protonation of the retinal Schiff base and initiates cascades of gradual conformational rearrangements of the protein. One cascade propagates in the extracellular direction and results in proton release, and the other in the cytoplasmic direction and results in side-chain and main-chain rearrangements, formation of a chain of hydrogen-bonded water, and proton uptake from the bulk. Such local-global conformational coupling, with gradual spreading of a local perturbation over the rest of the protein, might be the uniting principle of transporters and receptors.

Retinal conformation and dynamics in activation of rhodopsin illuminated by solid-state H NMR spectroscopy

Solid-state NMR spectroscopy gives a powerful avenue for investigating G protein-coupled receptors and other integral membrane proteins in a native-like environment. This article reviews the use of solid-state (2)H NMR to study the retinal cofactor of rhodopsin in the dark state as well as the meta I and meta II photointermediates. Site-specific (2)H NMR labels have been introduced into three regions (methyl groups) of retinal that are crucially important for the photochemical function of rhodopsin. Despite its phenomenal stability (2)H NMR spectroscopy indicates retinal undergoes rapid fluctuations within the protein binding cavity. The spectral lineshapes reveal the methyl groups spin rapidly about their three-fold (C(3)) axes with an order parameter for the off-axial motion of SC(3) approximately 0.9. For the dark state, the (2)H NMR structure of 11-cis-retinal manifests torsional twisting of both the polyene chain and the beta-ionone ring due to steric interactions of the ligand and the protein. Retinal is accommodated within the rhodopsin binding pocket with a negative pretwist about the C11=C12 double bond. Conformational distortion explains its rapid photochemistry and reveals the trajectory of the 11-cis to trans isomerization. In addition, (2)H NMR has been applied to study the retinylidene dynamics in the dark and light-activated states. Upon isomerization there are drastic changes in the mobility of all three methyl groups. The relaxation data support an activation mechanism whereby the beta-ionone ring of retinal stays in nearly the same environment, without a large displacement of the ligand. Interactions of the beta-ionone ring and the retinylidene Schiff base with the protein transmit the force of the retinal isomerization. Solid-state (2)H NMR thus provides information about the flow of energy that triggers changes in hydrogen-bonding networks and helix movements in the activation mechanism of the photoreceptor.

High tolerance for ionizable residues in the hydrophobic interior of proteins

Internal ionizable groups are quite rare in water-soluble globular proteins. Presumably, this reflects the incompatibility between charges and the hydrophobic environment in the protein interior. Here we show that proteins can have an inherently high tolerance for internal ionizable groups. The 25 internal positions in staphylococcal nuclease were substituted one at a time with Lys, Glu, or Asp without abolishing enzymatic activity and without detectable changes in the conformation of the protein. Similar results with substitutions of 6 randomly chosen internal positions in ribonuclease H with Lys and Glu suggest that the ability of proteins to tolerate internal ionizable groups might be a property common to many proteins. Eighty-six of the 87 substitutions made were destabilizing, but in all but one case the proteins remained in the native state at neutral pH. By comparing the stability of each variant protein at two different pH values it was established that the pK(a) values of most of the internal ionizable groups are shifted; many of the internal ionizable groups are probably neutral at physiological pH values. These studies demonstrate that special structural adaptations are not needed for ionizable groups to exist stably in the hydrophobic interior of proteins. The studies suggest that enzymes and other proteins that use internal ionizable groups for functional purposes could have evolved through the random accumulation of mutations that introduced ionizable groups at internal positions, followed by evolutionary adaptation and optimization to modulate stability, dynamics, and other factors necessary for function.

Crystallographic study of hydration of an internal cavity in engineered proteins with buried polar or ionizable groups

Although internal water molecules are essential for the structure and function of many proteins, the structural and physical factors that govern internal hydration are poorly understood. We have examined the molecular determinants of internal hydration systematically, by solving the crystal structures of variants of staphylococcal nuclease with Gln-66, Asn-66, and Tyr-66 at cryo (100 K) and room (298 K) temperatures, and comparing them with existing cryo and room temperature structures of variants with Glu-66, Asp-66, Lys-66, Glu-92 or Lys-92 obtained under conditions of pH where the internal ionizable groups are in the neutral state. At cryogenic temperatures the polar moieties of all these internal side chains are hydrated except in the cases of Lys-66 and Lys-92. At room temperature the internal water molecules were observed only in variants with Glu-66 and Tyr-66; water molecules in the other variants are probably present but they are disordered and therefore undetectable crystallographically. Each internal water molecule establishes between 3 and 5 hydrogen bonds with the protein or with other internal water molecules. The strength of interactions between internal polar side chains and water molecules seems to decrease from carboxylic acids to amides to amines. Low temperature, low cavity volume, and the presence of oxygen atoms in the cavity increase the positional stability of internal water molecules. This set of structures and the physical insight they contribute into internal hydration will be useful for the development and benchmarking of computational methods for artificial hydration of pockets, cavities, and active sites in proteins.

Characterisation of Schiff base and chromophore in green proteorhodopsin by solid-state NMR

The proteorhodopsin family consists of hundreds of homologous retinal containing membrane proteins found in bacteria in the photic zone of the oceans. They are colour tuned to their environment and act as light-driven proton pumps with a potential energetic and regulatory function. Precise structural details are still unknown. Here, the green proteorhodopsin variant has been selected for a chemical shift analysis of retinal and Schiff base by solid-state NMR. Our data show that the chromophore exists in mainly all-trans configuration in the proteorhodopsin ground state. The optical absorption maximum together with retinal and Schiff base chemical shifts indicate a strong interaction network between chromophore and opsin.

Solid-state 2H NMR spectroscopy of retinal proteins in aligned membranes

Solid-state 2H NMR spectroscopy gives a powerful avenue to investigating the structures of ligands and cofactors bound to integral membrane proteins. For bacteriorhodopsin (bR) and rhodopsin, retinal was site-specifically labeled by deuteration of the methyl groups followed by regeneration of the apoprotein. 2H NMR studies of aligned membrane samples were conducted under conditions where rotational and translational diffusion of the protein were absent on the NMR time scale. The theoretical lineshape treatment involved a static axial distribution of rotating C-C2H3 groups about the local membrane frame, together with the static axial distribution of the local normal relative to the average normal. Simulation of solid-state 2H NMR lineshapes gave both the methyl group orientations and the alignment disorder (mosaic spread) of the membrane stack. The methyl bond orientations provided the angular restraints for structural analysis. In the case of bR the retinal chromophore is nearly planar in the dark- and all-trans light-adapted states, as well upon isomerization to 13-cis in the M state. The C13-methyl group at the "business end" of the chromophore changes its orientation to the membrane upon photon absorption, moving towards W182 and thus driving the proton pump in energy conservation. Moreover, rhodopsin was studied as a prototype for G protein-coupled receptors (GPCRs) implicated in many biological responses in humans. In contrast to bR, the retinal chromophore of rhodopsin has an 11-cis conformation and is highly twisted in the dark state. Three sites of interaction affect the torsional deformation of retinal, viz. the protonated Schiff base with its carboxylate counterion; the C9-methyl group of the polyene; and the beta-ionone ring within its hydrophobic pocket. For rhodopsin, the strain energy and dynamics of retinal as established by 2H NMR are implicated in substituent control of activation. Retinal is locked in a conformation that is twisted in the direction of the photoisomerization, which explains the dark stability of rhodopsin and allows for ultra-fast isomerization upon absorption of a photon. Torsional strain is relaxed in the meta I state that precedes subsequent receptor activation. Comparison of the two retinal proteins using solid-state 2H NMR is thus illuminating in terms of their different biological functions.

Proton transfer 200 years after von Grotthuss: insights from ab initio simulations

In the last decade, ab initio simulations and especially Car-Parrinello molecular dynamics have significantly contributed to the improvement of our understanding of both the physical and chemical properties of water, ice, and hydrogen-bonded systems in general. At the heart of this family of in silico techniques lies the crucial idea of computing the many-body interactions by solving the electronic structure problem "on the fly" as the simulation proceeds, which circumvents the need for pre-parameterized potential models. In particular, the field of proton transfer in hydrogen-bonded networks greatly benefits from these technical advances. Here, several systems of seemingly quite different nature and of increasing complexity, such as Grotthuss diffusion in water, excited-state proton-transfer in solution, phase transitions in ice, and protonated water networks in the membrane protein bacteriorhodopsin, are discussed in the realms of a unifying viewpoint.

Dynamic structure of retinylidene ligand of rhodopsin probed by molecular simulations

Rhodopsin is currently the only available atomic-resolution template for understanding biological functions of the G protein-coupled receptor (GPCR) family. The structural basis for the phenomenal dark state stability of 11-cis-retinal bound to rhodopsin and its ultrafast photoreaction are active topics of research. In particular, the beta-ionone ring of the retinylidene inverse agonist is crucial for the activation mechanism. We analyzed a total of 23 independent, 100 ns all-atom molecular dynamics simulations of rhodopsin embedded in a lipid bilayer in the microcanonical (N,V,E) ensemble. Analysis of intramolecular fluctuations predicts hydrogen-out-of-plane (HOOP) wagging modes of retinal consistent with those found in Raman vibrational spectroscopy. We show that sampling and ergodicity of the ensemble of simulations are crucial for determining the distribution of conformers of retinal bound to rhodopsin. The polyene chain is rigidly locked into a single, twisted conformation, consistent with the function of retinal as an inverse agonist in the dark state. Most surprisingly, the beta-ionone ring is mobile within its binding pocket; interactions are non-specific and the cavity is sufficiently large to enable structural heterogeneity. We find that retinal occupies two distinct conformations in the dark state, contrary to most previous assumptions. The beta-ionone ring can rotate relative to the polyene chain, thereby populating both positively and negatively twisted 6-s-cis enantiomers. This result, while unexpected, strongly agrees with experimental solid-state (2)H NMR spectra. Correlation analysis identifies the residues most critical to controlling mobility of retinal; we find that Trp265 moves away from the ionone ring prior to any conformational transition. Our findings reinforce how molecular dynamics simulations can challenge conventional assumptions for interpreting experimental data, especially where existing models neglect conformational fluctuations.

Studies of the Bacteriorhodopsin Photocycle without the Use of Light: Clues to Proton Transfer Coupled Reactions

In the photochemical cycle of bacteriorhodopsin, the light-driven proton pump of halobacteria, only the first step, the isomerization of the all-trans retinal to 13-cis, is dependent on illumination. Because the steps that accomplish the translocation of a proton during the ensuing reaction sequence of intermediate states are thermal reactions, they have direct analogies with such steps in other ion pumps. In a surprisingly large number of cases, the reactions of the photocycle could be studied without using light. This review recounts experiments of this kind, and what they contribute to understanding the transport mechanism of this pump, and perhaps indirectly other ion pumps as well.

Electrostatic and steric interactions determine bacteriorhodopsin single-molecule biomechanics

Bacteriorhodopsin (bR) is a haloarchaeal membrane protein that converts the energy of single photons into large structural changes to directionally pump protons across purple membrane. This is achieved by a complex combination of local dynamic interactions controlling bR biomechanics at the submolecular level, producing efficient amplification of the retinal photoisomerization. Using single molecule force spectroscopy at different salt concentrations, we show that tryptophan (Trp) residues use steric specific interactions to create a rigid scaffold in bR extracellular region and are responsible for the main unfolding barriers. This scaffold, which encloses the retinal, controls bR local mechanical properties and anchors the protein into the membrane. Furthermore, the stable Trp-based network allows ion binding to two specific sites on the extracellular loops (BC and FG), which are involved in proton release and lateral transport. In contrast, the cytoplasmic side of bR is mainly governed by relatively weak nonspecific electrostatic interactions that provide the flexibility necessary for large cytoplasmic structural rearrangements during the photocycle. The presence of an extracellular Trp-based network tightly enclosing the retinal seems common to most haloarchaeal rhodopsins, and could be relevant to their exceptional efficiency.

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.

"Proton holes" in long-range proton transfer reactions in solution and enzymes: A theoretical analysis

Proton transfers are fundamental to chemical processes in solution and biological systems. Often, the well-known Grotthuss mechanism is assumed where a series of sequential "proton hops" initiates from the donor and combines to produce the net transfer of a positive charge over a long distance. Although direct experimental evidence for the sequential proton hopping has been obtained recently, alternative mechanisms may be possible in complex molecular systems. To understand these events, all accessible protonation states of the mediating groups should be considered. This is exemplified by transfers through water where the individual water molecules can exist in three protonation states (water, hydronium, and hydroxide); as a result, an alternative to the Grotthuss mechanism for a proton transfer through water is to generate a hydroxide by first protonating the acceptor and then transfer the hydroxide toward the donor through water. The latter mechanism can be most generally described as the transfer of a "proton hole" from the acceptor to the donor where the "hole" characterizes the deprotonated state of any mediating molecule. This pathway is distinct and is rarely considered in the discussion of proton-transfer processes. Using a calibrated quantum mechanical/molecular mechanical (QM/MM) model and an effective sampling technique, we study proton transfers in two solution systems and in Carbonic Anhydrase II. Although the relative weight of the "proton hole" and Grotthuss mechanisms in a specific system is difficult to determine precisely using any computational approach, the current study establishes an energetics motivated framework that hinges on the donor/acceptor pKa values and electrostatics due to the environment to argue that the "proton hole" transfer is likely as important as the classical Grotthuss mechanism for proton transport in many complex molecular systems.

Laser-induced transient grating analysis of dynamics of interaction between sensory rhodopsin II D75N and the HtrII transducer

The interaction between sensory rhodopsin II (SRII) and its transducer HtrII was studied by the time-resolved laser-induced transient grating method using the D75N mutant of SRII, which exhibits minimal visible light absorption changes during its photocycle, but mediates normal phototaxis responses. Flash-induced transient absorption spectra of transducer-free D75N and D75N joined to 120 amino-acid residues of the N-terminal part of the SRII transducer protein HtrII (DeltaHtrII) showed only one spectrally distinct K-like intermediate in their photocycles, but the transient grating method resolved four intermediates (K(1)-K(4)) distinct in their volumes. D75N bound to HtrII exhibited one additional slower kinetic species, which persists after complete recovery of the initial state as assessed by absorption changes in the UV-visible region. The kinetics indicate a conformationally changed form of the transducer portion (designated Tr*), which persists after the photoreceptor returns to the unphotolyzed state. The largest conformational change in the DeltaHtrII portion was found to cause a DeltaHtrII-dependent increase in volume rising in 8 micros in the K(4) state and a drastic decrease in the diffusion coefficient (D) of K(4) relatively to those of the unphotolyzed state and Tr*. The magnitude of the decrease in D indicates a large structural change, presumably in the solvent-exposed HAMP domain of DeltaHtrII, where rearrangement of interacting molecules in the solvent would substantially change friction between the protein and the solvent.

Three strategically placed hydrogen-bonding residues convert a proton pump into a sensory receptor

In haloarchaea, light-driven ion transporters have been modified by evolution to produce sensory receptors that relay light signals to transducer proteins controlling motility behavior. The proton pump bacteriorhodopsin and the phototaxis receptor sensory rhodopsin II (SRII) differ by 74% of their residues, with nearly all conserved residues within the photoreactive retinal-binding pocket in the membrane-embedded center of the proteins. Here, we show that three residues in bacteriorhodopsin replaced by the corresponding residues in SRII enable bacteriorhodopsin to efficiently relay the retinal photoisomerization signal to the SRII integral membrane transducer (HtrII) and induce robust phototaxis responses. A single replacement (Ala-215-Thr), bridging the retinal and the membrane-embedded surface, confers weak phototaxis signaling activity, and the additional two (surface substitutions Pro-200-Thr and Val-210-Tyr), expected to align bacteriorhodopsin and HtrII in similar juxtaposition as SRII and HtrII, greatly enhance the signaling. In SRII, the three residues form a chain of hydrogen bonds from the retinal's photoisomerized C(13)=C(14) double bond to residues in the membrane-embedded alpha-helices of HtrII. The results suggest a chemical mechanism for signaling that entails initial storage of energy of photoisomerization in SRII's hydrogen bond between Tyr-174, which is in contact with the retinal, and Thr-204, which borders residues on the SRII surface in contact with HtrII, followed by transfer of this chemical energy to drive structural transitions in the transducer helices. The results demonstrate that evolution accomplished an elegant but simple conversion: The essential differences between transport and signaling proteins in the rhodopsin family are far less than previously imagined.

Role of the cytoplasmic domain in Anabaena sensory rhodopsin photocycling: vectoriality of Schiff base deprotonation

Light-induced electric signals in intact E. coli cells generated by heterologously expressed full-length and C-terminally truncated versions of Anabaena sensory rhodopsin (ASR) demonstrate that the charge movements within the membrane-embedded part of the molecule are stringently controlled by the cytoplasmic domain. In particular, truncation inverts the direction of proton movement during Schiff base deprotonation from outward to cytoplasmic. Truncation also alters faster charge movements that occur before Schiff base deprotonation. Asp(217) as previously shown by FTIR serves as a proton acceptor in the truncated ASR but not in the full-length version, and its mutation to Asn restores the natural outward direction of proton movement. Introduction of a potential negative charge (Ser(86) to Asp) on the cytoplasmic side favors a cytoplasmic direction of proton release from the Schiff base. In contrast, mutation of the counterion Asp(75) to Glu reverses the photocurrent to the outward direction in the truncated pigment, and in both truncated and full-length versions accelerates Schiff base deprotonation more than 10-fold. The communication between the cytoplasmic domain and the membrane-embedded photoactive site of ASR demonstrated here is likely to derive from the receptor's use of a cytoplasmic protein for signal transduction, as has been suggested previously from binding studies.

G-protein-coupled receptor domain overexpression in Halobacterium salinarum: long-range transmembrane interactions in heptahelical membrane proteins

The aminergic alpha(2b)-adrenergic receptor (alpha(2b)-AR) third intracellular loop (alpha(2b)-AR 3i) mediates receptor subcellular compartmentalization and signal transduction processes via ligand-dependent interaction with G(i)- and G(o)- proteins. To understand the structural origins of these processes we engineered several lengths of alpha(2b)-AR 3i into the third intracellular loop of the proton pump bacteriorhodopsin (bR) and produced the fusion proteins in quantities suitable for physical studies. The fusion proteins were expressed in the Archaeon Halobacterium salinarum and purified. A highly expressed fusion protein was crystallized from bicelles and diffracted to low resolution on an in-house diffractometer. The bR-alpha(2b)-AR 3i(203-292) protein possessed a photocycle slightly perturbed from that of the wild-type bR. The first half of the fusion protein photocycle, correlated with proton release, is accelerated by a factor of 3, whereas the second half, correlated with proton uptake, is slightly slower than wild-type bR. In addition, there is a large decrease in the pK(a), (from 9.6 to 8.3) of the terminal proton release group in the unphotolyzed state of bR-alpha(2b)-AR 3i as deduced from the pH-dependence of the M-formation. Perturbation of a cytoplasmic loop has thus resulted in the perturbation of proton release at the extracellular surface. The current work indicates that long-range and highly coupled intramolecular interactions exist that are capable of "transducing" structural perturbations (e.g., signals) across the cellular membrane. This gene fusion approach may have general applicability for physical studies of G-protein-coupled receptor domains in the context of the bR structural scaffold.

Structural and energetic determinants of primary proton transfer in bacteriorhodopsin

In the light-driven bacteriorhodopsin proton pump, the first proton transfer step is from the retinal Schiff base to a nearby carboxylate group. The mechanism of this transfer step is highly controversial, in particular whether a direct proton jump is allowed. Here, we review the structural and energetic determinants of the direct proton transfer path computed by using a combined quantum mechanical/molecular mechanical approach. Both protein flexibility and electrostatic interactions play an important role in shaping the proton transfer energy profile. Detailed analysis of the energetics of putative transitions in the first half of the photocycle focuses on two elements that determine the likelihood that a given configuration of the active site is populated during the proton-pumping cycle. First, the rate-limiting barrier for proton transfer must be consistent with the kinetics of the photocycle. Second, the active-site configuration must be compatible with a productive overall pumping cycle.

Excitation of the L intermediate of bacteriorhodopsin: electric responses to test x-ray structures

The L intermediate of bacteriorhodopsin was excited, and its electrical response was measured. Two positive components were found in it with respect to the direction of proton pumping: an unresolved fast component, and a slower one (tau=7 micros) of small amplitude. The fast component was assigned to a charge motion corresponding to reisomerization of the retinal moiety, whereas the slow one was attributed to charge rearrangements reestablishing the ground state. Because three x-ray crystallographic structures have recently been reported for the L intermediate, it seemed important to calculate the intramolecular dipole moment changes associated to bR-->L for all three structures, so as to compare them with similar quantities determined from the electrical signals. The results are discussed in terms of amino acid side chains possibly contributing to the observed effect. We propose to use electrical signals as a verification tool for intermediate structures of the photocycle, and thus for molecular models of proton pumping.

Structural evolution of the chromophore in the primary stages of trans/cis isomerization in photoactive yellow protein

We have studied the structural changes induced by optical excitation of the chromophore in wild-type photoactive yellow protein (PYP) in liquid solution with a combined approach of polarization-sensitive ultrafast infrared spectroscopy and density functional theory calculations. We identify the nuC8-C9 marker modes for solution phase PYP in the P and I0 states, from which we derive that the first intermediate state I0 that appears with a 3 ps time constant can be characterized to have a cis geometry. This is the first unequivocal demonstration that the formation of I0 correlates with the conversion from the trans to the cis state. For the P and I0 states we compare the experimentally measured vibrational band patterns and anisotropies with calculations and find that for both trans and cis configurations the planarity of the chromophore has a strong influence. The C7=C8-(C9=O)-S moiety of the chromophore in the dark P state has a trans geometry with the C=O group slightly tilted out-of-plane, in accordance with the earlier reported structure obtained in an X-ray diffraction study of PYP crystals. In the case of I0, experiment and theory are only in agreement when the C7=C8-(C9=O)-S moiety has a planar configuration. We find that the carboxylic side group of Glu46 that is hydrogen-bonded to the chromophore phenolate oxygen does not alter its orientation on going from the electronic ground P state, via the electronic excited P state to the intermediate I0 state, providing conclusive experimental evidence that the primary stages of PYP photoisomerization involve flipping of the enone thioester linkage without significant relocation of the phenolate moiety.

ULTRAFAST CHEMISTRY: Using Time-Resolved Vibrational Spectroscopy for Interrogation of Structural Dynamics

Time-resolved infrared (IR) and Raman spectroscopy elucidates molecular structure evolution during ultrafast chemical reactions. Following vibrational marker modes in real time provides direct insight into the structural dynamics, as is evidenced in studies on intramolecular hydrogen transfer, bimolecular proton transfer, electron transfer, hydrogen bonding during solvation dynamics, bond fission in organometallic compounds and heme proteins, cis-trans isomerization in retinal proteins, and transformations in photochromic switch pairs. Femtosecond IR spectroscopy monitors the site-specific interactions in hydrogen bonds. Conversion between excited electronic states can be followed for intramolecular electron transfer by inspection of the fingerprint IR- or Raman-active vibrations in conjunction with quantum chemical calculations. Excess internal vibrational energy, generated either by optical excitation or by internal conversion from the electronic excited state to the ground state, is observable through transient frequency shifts of IR-active vibrations and through nonequilibrium populations as deduced by Raman resonances.

Regeneration and inhibition of proton pumping activity of bacteriorhodopsin blue membrane by cationic amine anesthetics

Bacteriorhodopsin (bR) is the prototype of an integral membrane protein with seven membrane-spanning alpha-helices and serves as a model of the G-protein-coupled drug receptors. This study is aimed at reaching a greater understanding of the role of amine local anesthetic cations on the proton transport in the bR protein, and furthermore, the functional role of "the cation" in the proton pumping mechanism. The effect of the amine anesthetic cations on the proton pump in the bR blue membrane was compared with those by divalent (Ca2+, Mg2+ and Mn2+) and monovalent metal cations (Li+, Na+, K+ and Cs+), which are essential for the correct functioning of the proton pumping of the bR protein. The results suggest that the interacting site of the divalent cation to the bR membrane may differ from that of the monovalent metal cation. The electric current profile of the bR blue membrane in the presence of the amine anesthetic cations was biphasic, involving the generation and inhibition of the proton pumping activity in a concentration-dependent manner. The extent of the regeneration of the proton pump by the additives increased in the order of monovalent metal cation<monovalent amine anesthetic cation<divalent metal cation. We found that organic cations such as the amine anesthetics can also regenerate the proton pump in the bR protein. The inhibition of proton transport in the bR protein by the anesthetic cations was elucidated using the wild type, the E204Q and the D96N mutated bRs. The hydrophobic interaction of the amine anesthetics with the bR protein plays an important part in inhibiting the bR proton pump.

Rational design of lipid for membrane protein crystallization

The lipidic cubic phase has been used to grow crystals of membrane proteins for high-resolution structure determination. However, the original, so-called, in meso method does not work reliably at low temperatures, where proteins are generally more stable, because the hosting lipid turns solid. The need existed therefore for a lipid that forms the cubic phase and that supports crystal growth at low temperatures. We created a database of phase diagrams and used it to design such a lipid. X-ray diffraction showed that the new lipid exhibits designed phase behavior. Further, it produces diffraction quality membrane protein crystals by the in meso method at 6 degrees C. This demonstrates that lipidic materials, like their protein counterparts are amenable to rational design. The same approach as used in this study should find application in extending the range of membrane proteins crystallizable by the in meso method and in tailoring transport of cubic phases for controlled delivery and uptake.

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.

Protons trap NR1/NR2B NMDA receptors in a nonconducting state

NMDA receptors are highly expressed in the CNS and are involved in excitatory synaptic transmission, as well as synaptic plasticity. Given that overstimulation of NMDA receptors can cause cell death, it is not surprising that these channels are under tight control by a series of inhibitory extracellular ions, including zinc, magnesium, and H+. We studied the inhibition by extracellular protons of recombinant NMDA receptor NR1/NR2B single-channel and macroscopic responses in transiently transfected human embryonic kidney HEK 293 cells using patch-clamp techniques. We report that proton inhibition proceeds identically in the absence or presence of agonist, which rules out the possibility that protonation inhibits receptors by altering coagonist binding. The response of macroscopic currents in excised patches to rapid jumps in pH was used to estimate the microscopic association and dissociation rates for protons, which were 1.4 x 10(9) m(-1) sec(-1) and 110-196 sec(-1), respectively (K(d) corresponds to pH 7.2). Protons reduce the open probability without altering the time course of desensitization or deactivation. Protons appear to slow at least one time constant describing the intra-activation shut-time histogram and modestly reduce channel open time, which we interpret to reflect a reduction in the overall channel activation rate and possible proton-induced termination of openings. This is consistent with a modest proton-dependent slowing of the macroscopic response rise time. From these data, we propose a physical model of proton inhibition that can describe macroscopic and single-channel properties of NMDA receptor function over a range of pH values.

Correlation of the O-intermediate rate with the pKa of Asp-75 in the dark, the counterion of the Schiff base of Pharaonis phoborhodopsin (sensory rhodopsin II)

Pharaonis phoborhodopsin (ppR), also called pharaonis sensory rhodopsin II, NpSRII, is a photoreceptor of negative phototaxis in Natronomonas (Natronobacterium) pharaonis. The photocycle rate of ppR is slow compared to that of bacteriorhodopsin, despite the similarity in their x-ray structures. The decreased rate of the photocycle of ppR is a result of the longer lifetime of later photo-intermediates such as M- (ppR(M)) and O-intermediates (ppR(O)). In this study, mutants were prepared in which mutated residues were located on the extracellular surface (P182, P183, and V194) and near the Schiff base (T204) including single, triple (P182S/P183E/V194T), and quadruple mutants. The decay of ppR(O) of the triple mutant was accelerated approximately 20-times from 690 ms for the wild-type to 36 ms. Additional mutation resulting in a triple mutant at the 204th position such as T204C or T204S further decreased the decay half-time to 6.6 or 8 ms, almost equal to that of bacteriorhodopsin. The decay half-times of the ppR(O) of mutants (11 species) and those of the wild-type were well-correlated with the pK(a) value of Asp-75 in the dark for the respective mutants as spectroscopically estimated, although there are some exceptions. The implications of these observations are discussed in detail.

Five Residues in the HtrI Transducer Membrane-proximal Domain Close the Cytoplasmic Proton-conducting Channel of Sensory Rhodopsin I

Transducer-free sensory rhodopsins carry out light-driven proton transport in Halobacterium salinarum membranes. Transducer binding converts the proton pumps to signal-relay devices in which the transport is inhibited. In sensory rhodopsin I (SRI) binding of its cognate transducer HtrI inhibits transport by closing a cytoplasmic proton-conducting channel necessary for proton uptake during the SRI photochemical reaction cycle. To investigate the channel closure, a series of HtrI mutants truncated in the membrane-proximal cytoplasmic portion of an SRI-HtrI fusion were constructed and expressed in H. salinarum membranes. We found that binding of the membrane-embedded portion of HtrI is insufficient for channel closure, whereas cytoplasmic extension of the second HtrI transmembrane helix by 13 residues blocks proton conduction through the channel as well as full-length HtrI. Specifically the closure activity is localized in this 13-residue membrane-proximal cytoplasmic domain to the 5 final residues, each of which incrementally contributes to reduction of proton conductivity. Moreover, these same residues in the dark incrementally and proportionally increase the pKa of the Asp-76 counterion to the protonated Schiff base chromophore in the membrane-embedded photoactive site. We conclude that this critical region of HtrI alters the dark conformation of SRI as well as light-induced channel opening. The 5 residues in HtrI correspond in position to 5 residues demonstrated on the homologous NpHtrII to interact with the E-F loop of its cognate receptor NpSRII in the accompanying article (Yang, C.-S., Sineshchekov, O., Spudich, E. N., and Spudich, J. L. (2004) J. Biol. Chem. 279, 42970-42976). These results strongly suggest that the membrane-proximal region of Htr proteins interact with their cognate sensory rhodopsin cytoplasmic domains as part of the signal-relay coupling between the proteins.

Bacteriorhodopsin

Fourier transform infrared and Raman spectroscopy, solid-state NMR, and X-ray crystallography have contributed detailed information about the structural changes in the proton transport cycle of the light-driven pump, bacteriorhodopsin. The results over the past few years add up to a step-by-step description of the configurational changes of the photoisomerized retinal, how these changes result in internal proton transfers and the release of a proton to the extracellular surface and uptake on the other side, as well as the conservation and transformation of excess free energy during the cycle.

X-ray crystallographic analysis of lipid-protein interactions in the bacteriorhodopsin purple membrane

The past decade has witnessed increasingly detailed insights into the structural mechanism of the bacteriorhodopsin photocycle. Concurrently, there has been much progress within our knowledge pertaining to the lipids of the purple membrane, including the discovery of new lipids and the overall effort to localize and identify each lipid within the purple membrane. Therefore, there is a need to classify this information to generalize the findings. We discuss the properties and roles of haloarchaeal lipids and present the structural data as individual case studies. Lipid-protein interactions are discussed in the context of structure-function relationships. A brief discussion of the possibility that bacteriorhodopsin functions as a light-driven inward hydroxide pump rather than an outward proton pump is also presented.