Bacteriorhodopsin: the functional details of a molecular machine are being resolved

A mid-infrared laser microscope for the time-resolved study of light-induced protein conformational changes

We have developed a confocal laser microscope operating in the mid-infrared range for the study of light-sensitive proteins, such as rhodopsins. The microscope features a co-aligned infrared and visible illumination path for the selective excitation and probing of proteins located in the IR focus only. An external-cavity tunable quantum cascade laser provides a wavelength tuning range (5.80-6.35 µm or 1570-1724 cm-1) suitable for studying protein conformational changes as a function of time delay after visible light excitation with a pulsed LED. Using cryogen-free detectors, the relative changes in the infrared absorption of rhodopsin thin films around 10-4 have been observed with a time resolution down to 30 ms. The measured full-width at half maximum of the Airy disk at λ = 6.08 µm in transmission mode with a confocal arrangement of apertures is 6.6 µm or 1.1λ. Dark-adapted sample replacement at the beginning of each photocycle is then enabled by exchanging the illuminated thin-film location with the microscope mapping stage synchronized to data acquisition and LED excitation and by averaging hundreds of time traces acquired in different nearby locations within a homogeneous film area. We demonstrate that this instrument provides crucial advantages for time-resolved IR studies of rhodopsin thin films with a slow photocycle. Time-resolved studies of inhomogeneous samples may also be possible with the presented instrument.

Common Mechanism of Activated Catalysis in P-loop Fold Nucleoside Triphosphatases-United in Diversity

To clarify the obscure hydrolysis mechanism of ubiquitous P-loop-fold nucleoside triphosphatases (Walker NTPases), we analysed the structures of 3136 catalytic sites with bound Mg-NTP complexes or their analogues. Our results are presented in two articles; here, in the second of them, we elucidated whether the Walker A and Walker B sequence motifs-common to all P-loop NTPases-could be directly involved in catalysis. We found that the hydrogen bonds (H-bonds) between the strictly conserved, Mg-coordinating Ser/Thr of the Walker A motif ([Ser/Thr]WA) and aspartate of the Walker B motif (AspWB) are particularly short (even as short as 2.4 ångströms) in the structures with bound transition state (TS) analogues. Given that a short H-bond implies parity in the pKa values of the H-bond partners, we suggest that, in response to the interactions of a P-loop NTPase with its cognate activating partner, a proton relocates from [Ser/Thr]WA to AspWB. The resulting anionic [Ser/Thr]WA alkoxide withdraws a proton from the catalytic water molecule, and the nascent hydroxyl attacks the gamma phosphate of NTP. When the gamma-phosphate breaks away, the trapped proton at AspWB passes by the Grotthuss relay via [Ser/Thr]WA to beta-phosphate and compensates for its developing negative charge that is thought to be responsible for the activation barrier of hydrolysis.

Structure-based insights into evolution of rhodopsins

Rhodopsins, most of which are proton pumps generating transmembrane electrochemical proton gradients, span all three domains of life, are abundant in the biosphere, and could play a crucial role in the early evolution of life on earth. Whereas archaeal and bacterial proton pumps are among the best structurally characterized proteins, rhodopsins from unicellular eukaryotes have not been well characterized. To fill this gap in the current understanding of the proton pumps and to gain insight into the evolution of rhodopsins using a structure-based approach, we performed a structural and functional analysis of the light-driven proton pump LR (Mac) from the pathogenic fungus Leptosphaeria maculans. The first high-resolution structure of fungi rhodopsin and its functional properties reveal the striking similarity of its membrane part to archaeal but not to bacterial rhodopsins. We show that an unusually long N-terminal region stabilizes the protein through direct interaction with its extracellular loop (ECL2). We compare to our knowledge all available structures and sequences of outward light-driven proton pumps and show that eukaryotic and archaeal proton pumps, most likely, share a common ancestor.

Zabelskii et al. present a structural and functional analysis of the lightdriven proton pump LR (Mac) from the fungus Leptosphaeria maculans. Their findings indicate that the archaeal ancestry of eukaryotic type 1 rhodopsins, and that the archaeal host of the proto-mitochondrial endosymbiont was capable of light-driven proton pumping.

Expansion of the "Sodium World" through Evolutionary Time and Taxonomic Space

In 1986, Vladimir Skulachev and his colleagues coined the term "Sodium World" for the group of diverse organisms with sodium (Na)-based bioenergetics. Albeit only few such organisms had been discovered by that time, the authors insightfully noted that "the great taxonomic variety of organisms employing the Na-cycle points to the ubiquitous distribution of this novel type of membrane-linked energy transductions". Here we used tools of bioinformatics to follow expansion of the Sodium World through the evolutionary time and taxonomic space. We searched for those membrane protein families in prokaryotic genomes that correlate with the use of the Na-potential for ATP synthesis by different organisms. In addition to the known Na-translocators, we found a plethora of uncharacterized protein families; most of them show no homology with studied proteins. In addition, we traced the presence of Na-based energetics in many novel archaeal and bacterial clades, which were recently identified by metagenomic techniques. The data obtained support the view that the Na-based energetics preceded the proton-dependent energetics in evolution and prevailed during the first two billion years of the Earth history before the oxygenation of atmosphere. Hence, the full capacity of Na-based energetics in prokaryotes remains largely unexplored. The Sodium World expanded owing to the acquisition of new functions by Na-translocating systems. Specifically, most classes of G-protein-coupled receptors (GPCRs), which are targeted by almost half of the known drugs, appear to evolve from the Na-translocating microbial rhodopsins. Thereby the GPCRs of class A, with 700 representatives in human genome, retained the Na-binding site in the center of the transmembrane heptahelical bundle together with the capacity of Na-translocation. Mathematical modeling showed that the class A GPCRs could use the energy of transmembrane Na-potential for increasing both their sensitivity and selectivity. Thus, GPCRs, the largest protein family coded by human genome, stem from the Sodium World, which encourages exploration of other Na-dependent enzymes of eukaryotes.

The Dynamics of the Neuropeptide Y Receptor Type 1 Investigated by Solid-State NMR and Molecular Dynamics Simulation

We report data on the structural dynamics of the neuropeptide Y (NPY) G-protein-coupled receptor (GPCR) type 1 (Y1R), a typical representative of class A peptide ligand GPCRs, using a combination of solid-state NMR and molecular dynamics (MD) simulation. First, the equilibrium dynamics of Y1R were studied using 15N-NMR and quantitative determination of 1H-13C order parameters through the measurement of dipolar couplings in separated-local-field NMR experiments. Order parameters reporting the amplitudes of the molecular motions of the C-H bond vectors of Y1R in DMPC membranes are 0.57 for the Cα sites and lower in the side chains (0.37 for the CH2 and 0.18 for the CH3 groups). Different NMR excitation schemes identify relatively rigid and also dynamic segments of the molecule. In monounsaturated membranes composed of longer lipid chains, Y1R is more rigid, attributed to a higher hydrophobic thickness of the lipid membrane. The presence of an antagonist or NPY has little influence on the amplitude of motions, whereas the addition of agonist and arrestin led to a pronounced rigidization. To investigate Y1R dynamics with site resolution, we conducted extensive all-atom MD simulations of the apo and antagonist-bound state. In each state, three replicas with a length of 20 μs (with one exception, where the trajectory length was 10 μs) were conducted. In these simulations, order parameters of each residue were determined and showed high values in the transmembrane helices, whereas the loops and termini exhibit much lower order. The extracellular helix segments undergo larger amplitude motions than their intracellular counterparts, whereas the opposite is observed for the loops, Helix 8, and termini. Only minor differences in order were observed between the apo and antagonist-bound state, whereas the time scale of the motions is shorter for the apo state. Although these relatively fast motions occurring with correlation times of ns up to a few µs have no direct relevance for receptor activation, it is believed that they represent the prerequisite for larger conformational transitions in proteins.

Dynamics of Bacteriorhodopsin in the Dark-Adapted State from Solution Nuclear Magnetic Resonance Spectroscopy

To achieve efficient proton pumping in the light-driven proton pump bacteriorhodopsin (bR), the protein must be tightly coupled to the retinal to rapidly convert retinal isomerization into protein structural rearrangements. Methyl group dynamics of bR embedded in lipid nanodiscs were determined in the dark-adapted state, and were found to be mostly well ordered at the cytosolic side. Methyl groups in the M145A mutant of bR, which displays only 10 % residual proton pumping activity, are less well ordered, suggesting a link between side-chain dynamics on the cytosolic side of the bR cavity and proton pumping activity. In addition, slow conformational exchange, attributed to low frequency motions of aromatic rings, was indirectly observed for residues on the extracellular side of the bR cavity. This may be related to reorganization of the water network. These observations provide a detailed picture of previously undescribed equilibrium dynamics on different time scales for ground-state bR.

Infrared Difference Spectroscopy of Proteins: From Bands to Bonds

Infrared difference spectroscopy probes vibrational changes of proteins upon their perturbation. Compared with other spectroscopic methods, it stands out by its sensitivity to the protonation state, H-bonding, and the conformation of different groups in proteins, including the peptide backbone, amino acid side chains, internal water molecules, or cofactors. In particular, the detection of protonation and H-bonding changes in a time-resolved manner, not easily obtained by other techniques, is one of the most successful applications of IR difference spectroscopy. The present review deals with the use of perturbations designed to specifically change the protein between two (or more) functionally relevant states, a strategy often referred to as reaction-induced IR difference spectroscopy. In the first half of this contribution, I review the technique of reaction-induced IR difference spectroscopy of proteins, with special emphasis given to the preparation of suitable samples and their characterization, strategies for the perturbation of proteins, and methodologies for time-resolved measurements (from nanoseconds to minutes). The second half of this contribution focuses on the spectral interpretation. It starts by reviewing how changes in H-bonding, medium polarity, and vibrational coupling affect vibrational frequencies, intensities, and bandwidths. It is followed by band assignments, a crucial aspect mostly performed with the help of isotopic labeling and site-directed mutagenesis, and complemented by integration and interpretation of the results in the context of the studied protein, an aspect increasingly supported by spectral calculations. Selected examples from the literature, predominately but not exclusively from retinal proteins, are used to illustrate the topics covered in this review.

Design of a light-gated proton channel based on the crystal structure of Coccomyxa rhodopsin

The light-driven proton pump Coccomyxa subellipsoidea rhodopsin (CsR) provides-because of its high expression in heterologous host cells-an opportunity to study active proton transport under controlled electrochemical conditions. In this study, solving crystal structure of CsR at 2.0-Å resolution enabled us to identify distinct features of the membrane protein that determine ion transport directivity and voltage sensitivity. A specific hydrogen bond between the highly conserved Arg⁸³ and the nearby nonconserved tyrosine (Tyr¹⁴) guided our structure-based transformation of CsR into an operational light-gated proton channel (CySeR) that could potentially be used in optogenetic assays. Time-resolved electrophysiological and spectroscopic measurements distinguished pump currents from channel currents in a single protein and emphasized the necessity of Arg⁸³ mobility in CsR as a dynamic extracellular barrier to prevent passive conductance. Our findings reveal that molecular constraints that distinguish pump from channel currents are structurally more confined than was generally expected. This knowledge might enable the structure-based design of novel optogenetic tools, which derive from microbial pumps and are therefore ion specific.

Protein dynamics observed by tunable mid-IR quantum cascade lasers across the time range from 10ns to 1s

We have developed a spectrometer based on tunable quantum cascade lasers (QCLs) for recording time-resolved absorption spectra of proteins in the mid-infrared range. We illustrate its performance by recording time-resolved difference spectra of bacteriorhodopsin in the carboxylic range (1800-1700cm^-1) and on the CO rebinding reaction of myoglobin (1960-1840cm^-1), at a spectral resolution of 1cm^-1. The spectrometric setup covers the time range from 4ns to nearly a second with a response time of 10-15ns. Absorption changes as low as 1×10^-4 are detected in single-shot experiments at t>1μs, and of 5×10^-6 in kinetics obtained after averaging 100 shots. While previous time-resolved IR experiments have mostly been conducted on hydrated films of proteins, we demonstrate here that the brilliance of tunable quantum cascade lasers is superior to perform ns time-resolved experiments even in aqueous solution (H₂O).

A natural light-driven inward proton pump

Light-driven outward H+ pumps are widely distributed in nature, converting sunlight energy into proton motive force. Here we report the characterization of an oppositely directed H+ pump with a similar architecture to outward pumps. A deep-ocean marine bacterium, Parvularcula oceani, contains three rhodopsins, one of which functions as a light-driven inward H+ pump when expressed in Escherichia coli and mouse neural cells. Detailed mechanistic analyses of the purified proteins reveal that small differences in the interactions established at the active centre determine the direction of primary H+ transfer. Outward H+ pumps establish strong electrostatic interactions between the primary H+ donor and the extracellular acceptor. In the inward H+ pump these electrostatic interactions are weaker, inducing a more relaxed chromophore structure that leads to the long-distance transfer of H+ to the cytoplasmic side. These results demonstrate an elaborate molecular design to control the direction of H+ transfers in proteins.

New ultrarapid-scanning interferometer for FT-IR spectroscopy with microsecond time-resolution

A novel Fourier-transform infrared (FT-IR) rapid-scan spectrometer has been developed (patent pending EP14194520.4) which yields 1000 times higher time resolution as compared to conventional rapid-scanning spectrometers. The central element to achieve faster scanning rates is based on a sonotrode whose front face represents the movable mirror of the interferometer. A prototype spectrometer with a time resolution of 13 μs was realized, capable of fully automated long-term measurements with a flow cell for liquid samples, here a photosynthetic membrane protein in solution. The performance of this novel spectrometer is demonstrated by recording the photoreaction of bacteriorhodopsin initiated by a short laser pulse that is synchronized to the data recording. The resulting data are critically compared to those obtained by step-scan spectroscopy and demonstrate the relevance of performing experiments on proteins in solution. The spectrometer allows for future investigations of fast, non-repetitive processes, whose investigation is challenging to step-scan FT-IR spectroscopy.

Conversion of a light-driven proton pump into a light-gated ion channel

Interest in microbial rhodopsins with ion pumping activity has been revitalized in the context of optogenetics, where light-driven ion pumps are used for cell hyperpolarization and voltage sensing. We identified an opsin-encoding gene (CsR) in the genome of the arctic alga Coccomyxa subellipsoidea C-169 that can produce large photocurrents in Xenopus oocytes. We used this property to analyze the function of individual residues in proton pumping. Modification of the highly conserved proton shuttling residue R83 or its interaction partner Y57 strongly reduced pumping power. Moreover, this mutation converted CsR at moderate electrochemical load into an operational proton channel with inward or outward rectification depending on the amino acid substitution. Together with molecular dynamics simulations, these data demonstrate that CsR-R83 and its interacting partner Y57 in conjunction with water molecules forms a proton shuttle that blocks passive proton flux during the dark-state but promotes proton movement uphill upon illumination.

Eukaryotic G protein-coupled receptors as descendants of prokaryotic sodium-translocating rhodopsins

Abstract

Microbial rhodopsins and G-protein coupled receptors (GPCRs, which include animal rhodopsins) are two distinct (super) families of heptahelical (7TM) membrane proteins that share obvious structural similarities but no significant sequence similarity. Comparison of the recently solved high-resolution structures of the sodium-translocating bacterial rhodopsin and various Na⁺-binding GPCRs revealed striking similarity of their sodium-binding sites. This similarity allowed us to construct a structure-guided sequence alignment for the two (super)families, which highlighted their evolutionary relatedness. Our analysis supports a common underlying molecular mechanism for both families that involves a highly conserved aromatic residue playing a pivotal role in rotation of the 6th transmembrane helix.

Reviewers

This article was reviewed by Oded Beja, G. P. S. Raghava and L. Aravind.

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-015-0091-4) contains supplementary material, which is available to authorized users.

Proton transfer and protein conformation dynamics in photosensitive proteins by time-resolved step-scan Fourier-transform infrared spectroscopy

Monitoring the dynamics of protonation and protein backbone conformation changes during the function of a protein is an essential step towards understanding its mechanism. Protonation and conformational changes affect the vibration pattern of amino acid side chains and of the peptide bond, respectively, both of which can be probed by infrared (IR) difference spectroscopy. For proteins whose function can be repetitively and reproducibly triggered by light, it is possible to obtain infrared difference spectra with (sub)microsecond resolution over a broad spectral range using the step-scan Fourier transform infrared technique. With -10(2)-10(3) repetitions of the photoreaction, the minimum number to complete a scan at reasonable spectral resolution and bandwidth, the noise level in the absorption difference spectra can be as low as -10(-) (4), sufficient to follow the kinetics of protonation changes from a single amino acid. Lower noise levels can be accomplished by more data averaging and/or mathematical processing. The amount of protein required for optimal results is between 5-100 µg, depending on the sampling technique used. Regarding additional requirements, the protein needs to be first concentrated in a low ionic strength buffer and then dried to form a film. The protein film is hydrated prior to the experiment, either with little droplets of water or under controlled atmospheric humidity. The attained hydration level (g of water / g of protein) is gauged from an IR absorption spectrum. To showcase the technique, we studied the photocycle of the light-driven proton-pump bacteriorhodopsin in its native purple membrane environment, and of the light-gated ion channel channelrhodopsin-2 solubilized in detergent.

The G-protein-coupled neuropeptide Y receptor type 2 is highly dynamic in lipid membranes as revealed by solid-state NMR spectroscopy

In spite of the recent success in crystallizing several G-protein-coupled receptors (GPCRs), a comprehensive biophysical characterization of these molecules under physiological conditions also requires the study of the molecular dynamics of these proteins. The molecular mobility of the human neuropeptide Y receptor type 2 reconstituted into dimyristoylphosphatidylcholine (DMPC) membranes was investigated by means of solid-state NMR spectroscopy. Static (15) N NMR spectra show that the receptor performs axially symmetric motions in the membrane, and several residues undergo large amplitude fluctuations. This was confirmed by quantitative measurements of the motional (1) H,(13) C order parameter of the CH, CH2 , and CH3 groups. In directly polarized (13) C NMR experiments, these order parameters showed astonishingly low values of SCH =0.55, S CH 2=0.33, and S CH 3=0.17, which corresponds to segmental amplitudes of approximately 50° in the backbone and approximately 50-60° in the side chain. At physiological temperature, (2) H NMR spectra of the deuterated receptor showed a narrow component that is indicative of molecular order parameters of S≤0.3 superimposed with a very broad spectrum that could stem from the transmembrane α-helices. These results suggest that the crystal structures of GPCRs only represent a static snapshot of these highly mobile molecules, which undergo significant structural fluctuations with relatively large amplitudes in a liquid-crystalline membrane at physiological temperature.

Creating electrochemical gradients by light: from bio-inspired concepts to photoelectric conversion

Artificial light harvesting systems can be used to convert light into electrochemical gradients and photocurrents.

Light-induced Difference Terahertz Spectroscopy

Visible/near-infraredlaser-induced difference spectroscopy basedon a time-domain terahertz system has beendeveloped, and used to study copperpathancyonine. We find that the absorptionpeak of this molecule at 1.08 THz changessignificantly under 790 nm laserexcitation, suggesting that we haveobserved the first evidence of vibrationalmode changes in the THz range induced byvisible/near-infrared light.

Conformational dynamics of a membrane transport protein probed by H/D exchange and covalent labeling: the glycerol facilitator

Glycerol facilitator (GF) is a tetrameric membrane protein responsible for the selective permeation of glycerol and water. Each of the four GF subunits forms a transmembrane channel. Every subunit consists of six helices that completely span the lipid bilayer, as well as two half-helices (TM7 and TM3). X-ray crystallography has revealed that the selectivity of GF is due to its unique amphipathic channel interior. To explore the structural dynamics of GF, we employ hydrogen/deuterium exchange (HDX) and oxidative labeling with mass spectrometry (MS). HDX-MS reveals that transmembrane helices are generally more protected than extramembrane segments, consistent with data previously obtained for other membrane proteins. Interestingly, TM7 does not follow this trend. Instead, this half-helix undergoes rapid deuteration, indicative of a highly dynamic local structure. The oxidative labeling behavior of most GF residues is consistent with the static crystal structure. However, the side chains of C134 and M237 undergo labeling although they should be inaccessible according to the X-ray structure. In agreement with our HDX-MS data, this observation attests to the fact that TM7 is only marginally stable. We propose that the highly mobile nature of TM7 aids in the efficient diffusion of guest molecules through the channel ("molecular lubrication"). In the absence of such dynamics, host-guest molecular recognition would favor semipermanent binding of molecules inside the channel, thereby impeding transport. The current work highlights the complementary nature of HDX, covalent labeling, and X-ray crystallography for the characterization of membrane proteins.

Hydrogen exchange mass spectrometry of bacteriorhodopsin reveals light-induced changes in the structural dynamics of a biomolecular machine

Many proteins act as molecular machines that are fuelled by a nonthermal energy source. Examples include transmembrane pumps and stator-rotor complexes. These systems undergo cyclic motions (CMs) that are being driven along a well-defined conformational trajectory. Superimposed on these CMs are thermal fluctuations (TFs) that are coupled to stochastic motions of the solvent. Here we explore whether the TFs of a molecular machine are affected by the occurrence of CMs. Bacteriorhodopsin (BR) is a light-driven proton pump that serves as a model system in this study. The function of BR is based on a photocycle that involves trans/cis isomerization of a retinal chromophore, as well as motions of transmembrane helices. Hydrogen/deuterium exchange (HDX) mass spectrometry was used to monitor the TFs of BR, focusing on the monomeric form of the protein. Comparative HDX studies were conducted under illumination and in the dark. The HDX kinetics of BR are dramatically accelerated in the presence of light. The isotope exchange rates and the number of backbone amides involved in EX2 opening transitions increase roughly 2-fold upon illumination. In contrast, light/dark control experiments on retinal-free protein produced no discernible differences. It can be concluded that the extent of TFs in BR strongly depends on photon-driven CMs. The light-induced differences in HDX behavior are ascribed to protein destabilization. Specifically, the thermodynamic stability of the dark-adapted protein is estimated to be 5.5 kJ mol(-1) under the conditions of our work. This value represents the free energy difference between the folded state F and a significantly unfolded conformer U. Illumination reduces the stability of F by 2.2 kJ mol(-1). Mechanical agitation caused by isomerization of the chromophore is transferred to the surrounding protein scaffold, and subsequently, the energy dissipates into the solvent. Light-induced retinal motions therefore act analogously to an internal heat source that promotes the occurrence of TFs. Overall, our data highlight the potential of HDX methods for probing the structural dynamics of molecular machines under "engine on" and "engine off" conditions.

Time-resolved FT-IR Spectroscopy of Membrane Proteins

Time-resolved Fourier transform infrared spectroscopy (FT-IR) offers distinct advantages concerning restrictions pertinent to biomolecules. In particular, it is possible to monitor the temporal evolution of the reaction mechanism of complex machineries as membrane proteins, where other techniques encounter significant experimental difficulties. Here, we present the classical principles and experimental realizations of time-resolved FT-IR spectroscopy together with recent developments employed in our laboratory. Examples from applications to retinal proteins are reviewed that underline the impact of time-resolved FT-IR spectroscopy on the understanding of protein reactions on the level of single bonds.

Site-directed mutagenesis combined with oxidative methionine labeling for probing structural transitions of a membrane protein by mass spectrometry

Exposure of the membrane protein bacteriorhodopsin (BR) to SDS induces partial breakdown of the native conformation. The exact structural properties of this SDS state remain a matter of debate, despite its widespread use in BR folding experiments. The current work employs hydroxyl radical (·OH) labeling in conjunction with mass spectrometry (MS)-based peptide mapping for probing the solvent accessibility of individual BR segments in the presence of SDS. Previous work revealed methionine sulfoxide formation to be the dominant oxidative pathway. Those data suggested extensive unfolding of helices A and D in SDS. Unfortunately, the lack of Met residues in helices C and F implies that no direct information on the behavior of the latter two elements could be obtained. Here, we address this problem by employing two variants with additional Met residues, L93M (helix C) and V179M (helix F). The oxidation behavior of the resulting 11 methionines can be grouped into three categories: (1) extensively labeled both in native BR and in SDS (loop residues M32, M68, and M163); (2) protected in the native state but not in SDS (M20, M118); (3) always protected (M56, M60, M93, M145, M179, M209). These data show that a solvent-inaccessible core is retained in SDS. This core consists of partially intact helices B, C, E, F, and G. The termini of these helices are highly dynamic and/or unraveled, particularly on the cytoplasmic side. Overall, this work demonstrates how the use of engineered ·OH labeling sites can provide insights into structural properties of membrane proteins.

Dynamics of reassembled thioredoxin studied by magic angle spinning NMR: snapshots from different time scales

Solid-state NMR spectroscopy can be used to probe internal protein dynamics in the absence of the overall molecular tumbling. In this study, we report (15)N backbone dynamics in differentially enriched 1-73(U-(13)C,(15)N)/74-108(U-(15)N) reassembled thioredoxin on multiple time scales using a series of 2D and 3D MAS NMR experiments probing the backbone amide (15)N longitudinal relaxation, (1)H-(15)N dipolar order parameters, (15)N chemical shift anisotropy (CSA), and signal intensities in the temperature-dependent and (1)H T(2)'-filtered NCA experiments. The spin-lattice relaxation rates R(1) (R(1) = 1/T(1)) were observed in the range from 0.012 to 0.64 s(-1), indicating large site-to-site variations in dynamics on pico- to nanosecond time scales. The (1)H-(15)N dipolar order parameters, <S>, and (15)N CSA anisotropies, delta(sigma), reveal the backbone mobilities in reassembled thioredoxin, as reflected in the average <S> = 0.89 +/- 0.06 and delta(sigma) = 92.3 +/- 5.2 ppm, respectively. From the aggregate of experimental data from different dynamics methods, some degree of correlation between the motions on the different time scales has been suggested. Analysis of the dynamics parameters derived from these solid-state NMR experiments indicates higher mobilities for the residues constituting irregular secondary structure elements than for those located in the alpha-helices and beta-sheets, with no apparent systematic differences in dynamics between the alpha-helical and beta-sheet residues. Remarkably, the dipolar order parameters derived from the solid-state NMR measurements and the corresponding solution NMR generalized order parameters display similar qualitative trends as a function of the residue number. The comparison of the solid-state dynamics parameters to the crystallographic B-factors has identified the contribution of static disorder to the B-factors. The combination of longitudinal relaxation, dipolar order parameter, and CSA line shape analyses employed in this study provides snapshots of dynamics and a new insight on the correlation of these motions on multiple time scales.

Side chain resonances in static oriented proton-decoupled 15N solid-state NMR spectra of membrane proteins

Proton-decoupled (15)N solid-state NMR spectra are used to analyze the structure, dynamics, and membrane topology of proteins uniformly labeled with (15)N. Preparation of the proteins by bacterial overexpression results in the labeling not only of the backbone amides but also of nitrogens localized within the side chains of arginine, glutamine, tryptophan, asparagines, lysines, and histidines. Most of these side chain resonances appear in the spectral region of the anisotropic backbone amides, and residual intensities have been observed also in cross-polarization spectra. In the past this issue has received little attention although it can cause ambiguities during assignment. Here we show that by combining cross-polarization and Hahn echo solid-state NMR experiments, it is possible to differentiate between side chain and backbone resonances. This is demonstrated using experimental and simulated (15)N spectra of oriented purple membranes, diphtheria toxin T domain and Bcl-x(L).

Monitoring light-induced structural changes of Channelrhodopsin-2 by UV-visible and Fourier transform infrared spectroscopy

Channelrhodopsin-2 (ChR2) is a microbial type rhodopsin and a light-gated cation channel that controls phototaxis in Chlamydomonas. We expressed ChR2 in COS-cells, purified it, and subsequently investigated this unusual photoreceptor by flash photolysis and UV-visible and Fourier transform infrared difference spectroscopy. Several transient photoproducts of the wild type ChR2 were identified, and their kinetics and molecular properties were compared with those of the ChR2 mutant E90Q. Based on the spectroscopic data we developed a model of the photocycle comprising six distinguishable intermediates. This photocycle shows similarities to the photocycle of the ChR2-related Channelrhodopsin of Volvox but also displays significant differences. We show that molecular changes include retinal isomerization, changes in hydrogen bonding of carboxylic acids, and large alterations of the protein backbone structure. These alterations are stronger than those observed in the photocycle of other microbial rhodopsins like bacteriorhodopsin and are related to those occurring in animal rhodopsins. UV-visible and Fourier transform infrared difference spectroscopy revealed two late intermediates with different time constants of tau = 6 and 40 s that exist during the recovery of the dark state. The carboxylic side chain of Glu(90) is involved in the slow transition. The molecular changes during the ChR2 photocycle are discussed with respect to other members of the rhodopsin family.

Hydration dynamics in a partially denatured ensemble of the globular protein human alpha-lactalbumin investigated with molecular dynamics simulations

Atomistic molecular dynamics simulations are used to probe changes in the nature and subnanosecond dynamical behavior of solvation waters that accompany partial denaturation of the globular protein, human alpha-lactalbumin. A simulated ensemble of subcompact conformers, similar to the molten globule state of human alpha-lactalbumin, demonstrates a marginal increase in the amount of surface solvation relative to the native state. This increase is accompanied by subtle but distinct enhancement in surface water dynamics, less favorable protein-water interactions, and a marginal decrease in the anomalous behavior of solvation water dynamics. The extent of solvent influx is not proportional to the increased surface area, and the partially denatured conformers are less uniformly solvated compared to their native counterpart. The observed solvation in partially denatured conformers is lesser in extent compared to earlier experimental estimates in molten globule states, and is consistent with more recent descriptions based on nuclear magnetic relaxation dispersion studies.

An extended dead-end elimination algorithm to determine gap-free lists of low energy states

Proteins are flexible systems and commonly populate several functionally important states. To understand protein function, these states and their energies have to be identified. We introduce an algorithm that allows the determination of a gap-free list of the low energy states. This algorithm is based on the dead-end elimination (DEE) theorem and is termed X-DEE (extended DEE). X-DEE is applicable to discrete systems whose state energy can be formulated as pairwise interaction between sites and their intrinsic energies. In this article, the computational performance of X-DEE is analyzed and discussed. X-DEE is implemented to determine the lowest energy protonation states of proteins, a problem to which DEE has not been applied so far. We use X-DEE to calculate a list of low energy protonation states for two bacteriorhodopsin structures that represent the first proton transfer step of the bacteriorhodopsin photocycle.

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.

Collagen dynamics in articular cartilage under osmotic pressure

Cartilage is a complex biological tissue consisting of collagen, proteoglycans and water. The structure and molecular mobility of the collagen component of cartilage were studied by (13)C solid-state NMR spectroscopy as a function of hydration. The hydration level of cartilage was adjusted between fully hydrated ( approximately 80 wt% H(2)O) and highly dehydrated ( approximately 30 wt% H(2)O) using the osmotic stress technique. Thus, the conditions of mechanical load could be simulated and the response of the tissue macromolecules to mechanical stress is reported. From the NMR measurements, the following results were obtained. (i) Measurements of motionally averaged dipolar (1)H-(13)C couplings were carried out to study the segmental mobility in cartilage collagen at full hydration. Backbone segments undergo fast motions with amplitudes of approximately 35 degrees whereas the collagen side-chains are somewhat more mobile with amplitudes between 40 and 50 degrees . In spite of the high water content of cartilage, collagen remains essentially rigid. (ii) No chemical shift changes were observed in (13)C cross-polarization magic angle spinning spectra of cartilage tissue at varying hydration indicating that the collagen structure was not altered by application of high osmotic stress. (iii) The (1)H-(13)C dipolar coupling values detected for collagen signals respond to dehydration. The dipolar coupling values gradually increase upon cartilage dehydration, reaching rigid limit values at approximately 30 wt% H(2)O. This indicates that collagen is essentially dehydrated in cartilage tissue under very high mechanical load, which provides insights into the elastic properties of cartilage collagen, although the mechanical pressures applied here exceed the physiological limit.

Nonequilibrium membrane fluctuations driven by active proteins

We extend a model for nonthermal membrane undulations driven by active (adenosine triphosphate-dependent or light-harvesting) membrane proteins [N. Gov, Phys. Rev. Lett. 93, 268104 (2004)]. The present model accounts for the fact that proteins can diffuse laterally across the membrane surface and that individual proteins are expected to exert forces preferentially in one normal direction over the other (due to their orientation within the bilayer). The addition of these effects alters the scaling of fluctuation amplitudes with system size. Additionally, theoretical arguments and dynamic simulations both suggest that, in certain regimes, the probability distribution of fluctuation amplitudes is expected to be non-Gaussian (in contrast to thermal systems).

Bacteriorhodopsin-like proteins of eubacteria and fungi: the extent of conservation of the haloarchaeal proton-pumping mechanism

A stereotypical image of a retinal-binding proton pump derived from extensive studies of halobacterial ion-transporting and sensory rhodopsins is a fast-cycling protein which possesses two strategically placed carboxylic acids serving as proton donor and acceptor for the retinal Schiff base. We review recent biophysical and bioinformatic data on the novel eubacterial and eucaryotic rhodopsins to analyze the extent of conservation of the haloarchaeal mechanism of transmembrane proton transport. We show that only the most essential elements of the haloarchaeal proton-pumping machinery are conserved universally, and that a mere presence of these elements in primary structures does not guarantee the proton-pumping ability.

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

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

Time-Resolved Step-Scan FT-IR Spectroscopy: Focus on Multivariate Curve Resolution

The present paper describes the application of step-scan FT-IR spectroscopy in combination with chemometric analysis of the spectral data for the study of the photocycle of bacteriorhodopsin. The focus is on the performance of this instrumentation for time-resolved experiments. Three-dimensional data-spectra recorded over time-are studied using various factor analysis techniques, e.g., singular values decomposition, evolving factor analysis, and multivariate curve resolution based on alternating least squares. Transient intermediates formed in the time domain ranging from 1 micros to 6.6 ms are clearly detected through reliable pure time evolving profiles. At the same time, pure difference absorbance spectra are provided. As a result, valuable information about transitions and dynamics of the protein can be extracted. We conclude first that step-scan FT-IR spectroscopy is a useful technique for the direct study of difficult photochemical systems. Second, and this is the essential motivation of this paper, chemometrics provide a step forward in the description of the photointermediates.

Activity, stability and conformational flexibility of seed coat soybean peroxidase

Seed coat soybean peroxidase (SBP) belongs to class III of the plant peroxidase superfamily that includes the classical peroxidase, namely horseradish peroxidase (HRP). We have measured the catalytic activity (k(cat)) and catalytic efficiency (k(cat)/K(M)) of SBP and that of HRP-C for the oxidation of ABTS [2,2'-azino-bis-(3-ethylbenzthiazoline-6-sulphonate)] by hydrogen peroxide at 25 degrees C. We observed that the k(cat) and k(cat)/K(M) values for SBP are much higher than those for HRP-C at all pH values, rendering SBP a more potent peroxidase. This is attributed to the relatively more solvent exposed delta-meso heme edge in SBP. We observed that the maximum catalytic activity and conformational stability of SBP is at pH approximately 5.5. A pH maximum of 5.0 for the catalytic activity of SBP has recently been reported. Estimation of secondary structural elements at various pH values indicated that there is a maximal reduction of beta-strands and beta-turns at pH 5.5 causing the heme to be further exposed to the solvent and increasing the overall conformational flexibility of the protein.

Vibrational spectroscopy of an algal Phot-LOV1 domain probes the molecular changes associated with blue-light reception

The LOV1 domain of the blue light Phot1-receptor (phototropin homolog) from Chlamydomonas reinhardtii has been studied by vibrational spectroscopy. The FMN modes of the dark state of LOV1 were identified by preresonance Raman spectroscopy and assigned to molecular vibrations. By comparing the blue-light-induced FTIR difference spectrum with the preresonance Raman spectrum, most of the differences are due to FMN modes. Thus, we exclude large backbone changes of the protein that might occur during the phototransformation of the dark state LOV1-447 into the putative signaling state LOV1-390. Still, the presence of smaller amide difference bands cannot be excluded but may be masked by overlapping FMN modes. The band at 2567 cm(-1) is assigned to the S-H stretching vibration of C57, the residue that forms the transient thio-adduct with the chromophore FMN. The occurrence of this band is evidence that C57 is protonated in the dark state of LOV1. This result challenges conclusions from the homologous LOV2 domain from oat that the thiolate of the corresponding cysteine is the reactive species.

Thermal unfolding of soybean peroxidase. Appropriate high denaturant concentrations induce cooperativity allowing the correct measurement of thermodynamic parameters

We have earlier reported that both guanidine hydrochloride (GdnHCl)-induced and heat-induced unfolding of seed coat soybean peroxidase (SBP), monitored by far UV CD, show single step transition. However, although GdnHCl-induced unfolding follows a two-state pathway, the heat-induced denaturation proceeds through intermediates as indicated by the very low cooperativity of transition. In the former case, analysis of the data based on the two-state model gives true thermodynamic parameters, whereas underestimated values are obtained in the latter case. Available complex equations also cannot be applied for the analysis of the thermal unfolding of SBP due to the absence of separate transitions for the intermediates. In the present study, we report a method to obtain true thermodynamic parameters from thermal transition curves of SBP using the two-state model. When SBP is subjected to thermal unfolding at high GdnHCl concentrations (5.8-6.9 M), cooperative behavior is observed, which allowed the analysis by the two-state model to determine their thermodynamic parameters. We then obtained the thermodynamic parameters in the absence of GdnHCl by extrapolating the graph of linear dependence of DeltaH(m) on T(m) to the T(m) corresponding to 0 m GdnHCl. Another key point for checking the validity of our method was the fact that the unfolded state of SBP generated by either heat or GdnHCl is the same by which we could cross-check our results with that obtained from GdnHCl unfolding. Having obtained the true thermodynamic parameters, we report a detailed thermodynamic study of SBP. Further we address the effect of heme in the thermal unfolding mechanism of SBP.

Proton transfer dynamics on the surface of the late M state of bacteriorhodopsin

The cytoplasmic surface of the BR (initial) state of bacteriorhodopsin is characterized by a cluster of three carboxylates that function as a proton-collecting antenna. Systematic replacement of most of the surface carboxylates indicated that the cluster is made of D104, E161, and E234 (Checover, S., Y. Marantz, E. Nachliel, M. Gutman, M. Pfeiffer, J. Tittor, D. Oesterhelt, and N. Dencher. 2001. Biochemistry. 40:4281-4292), yet the BR state is a resting configuration; thus, its proton-collecting antenna can only indicate the presence of its role in the photo-intermediates where the protein is re-protonated by protons coming from the cytoplasmic matrix. In the present study we used the D96N and the triple (D96G/F171C/F219L) mutant for monitoring the proton-collecting properties of the protein in its late M state. The protein was maintained in a steady M state by continuous illumination and subjected to reversible pulse protonation caused by repeated excitation of pyranine present in the reaction mixture. The re-protonation dynamics of the pyranine anion was subjected to kinetic analysis, and the rate constants of the reaction of free protons with the surface groups and the proton exchange reactions between them were calculated. The reconstruction of the experimental signal indicated that the late M state of bacteriorhodopsin exhibits an efficient mechanism of proton delivery to the unoccupied-most basic-residue on its cytoplasmic surface (D38), which exceeds that of the BR configuration of the protein. The kinetic analysis was carried out in conjunction with the published structure of the M state (Sass, H., G. Büldt, R. Gessenich, D. Hehn, D. Neff, R. Schlesinger, J. Berendzen, and P. Ormos. 2000. Nature. 406:649-653), the model that resolves most of the cytoplasmic surface. The combination of the kinetic analysis and the structural information led to identification of two proton-conducting tracks on the protein's surface that are funneling protons to D38. One track is made of the carboxylate moieties of residues D36 and E237, while the other is made of D102 and E232. In the late M state the carboxylates of both tracks are closer to D38 than in the BR (initial) state, accounting for a more efficient proton equilibration between the bulk and the protein's proton entrance channel. The triple mutant resembles in the kinetic properties of its proton conducting surface more the BR-M state than the initial state confirming structural similarities with the BR-M state and differences to the BR initial state.

Proton translocation by bacteriorhodopsin in the absence of substantial conformational changes

Unlike wild-type bacteriorhodopsin (BR), the BR triple mutant D96G/F171C/F219L has been shown to undergo only minor structural rearrangements during its photocycle. Nonetheless, the mutant is capable of transporting protons at a rate of 125(+/-40) H+/BR per minute under light-saturating conditions. Light adaptation of the triple mutant's retinal proceeds in a pH-dependent manner up to a maximum of 63% all-trans. These two findings imply that the transport activity of the triple mutant comprises 66% of the wild-type activity. Time-resolved spectroscopy reveals that the identity and sequence of intermediates in the photocycle of the triple mutant in the all-trans configuration correspond to that of wild-type BR. The only differences relate to a slower rise and decay of the M and O intermediates, and a significant spectral contribution from a 13-cis component. No indication for accumulation of the N intermediate is found under a variety of conditions that normally favor the formation of this species in wild-type BR. The Fourier transform infrared (FTIR) spectrum of the M intermediate in the triple mutant resembles that of wild type. Minor changes in the amide I region during the photocycle suggest that only small movements of the protein backbone occur. Electron microscopy reveals large differences in conformation between the unilluminated state of the mutant protein and wild-type but no light-induced changes in time-resolved measurements. Evidently, proton transport by the triple mutant does not require the major conformational rearrangements that occur on the same time-scale with wild-type. Thus, we conclude that large conformational changes observed in the photocycle of the wild-type and many BR mutants are not a prerequisite for the change in accessibility of the Schiff base nitrogen atom that must occur during vectorial catalysis to allow proton transport.