High resolution 13C-solid state NMR of bacteriorhodopsin: assignment of specific aspartic acids and structural implications of single site mutations

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.

Reconstruction, modeling & analysis of Halobacterium salinarum R-1 metabolism

We present a genome-scale metabolic reconstruction for the extreme halophile Halobacterium salinarum. The reconstruction represents a summary of the knowledge regarding the organism's metabolism, and has already led to new research directions and improved the existing annotation. We used the network for computational analysis and studied the aerobic growth of the organism using dynamic simulations in media with 15 available carbon and energy sources. Simulations resulted in predictions for the internal fluxes, which describe at the molecular level how the organism lives and grows. We found numerous indications that cells maximized energy production even at the cost of longer term concerns such as growth prospects. Simulations showed a very low carbon incorporation rate of only approximately 15%. All of the supplied nutrients were simultaneously degraded, unexpectedly including five which are essential. These initially surprising behaviors are likely adaptations of the organism to its natural environment where growth occurs in blooms. In addition, we also examined specific aspects of metabolism, including how each of the supplied carbon and energy sources is utilized. Finally, we investigated the consequences of the model assumptions and the network structure on the quality of the flux predictions.

Magnetic resonance studies of the bacteriorhodopsin pump cycle

Active transport requires the alternation of substrate uptake and release with a switch in the access of the substrate binding site to the two sides of the membrane. Both the transfer and switch aspects of the photocycle have been subjects of magnetic resonance studies in bacteriorhodopsin. The results for ion transfer indicate that the Schiff base of the chromophore is hydrogen bonded before, during, and after its deprotonation. This suggests that the initial complex counterion of the Schiff base decomposes in such a way that the Schiff base carries its immediate hydrogen-bonding partner with it as it rotates during the first half of the photocycle. If so, bacteriorhodopsin acts as an inward-directed hydroxide pump rather than as an outward-directed proton pump. The studies of the access switch explore both protein-based and chromophore-based mechanisms. Combined with evidence from functional studies of mutants and other forms of spectroscopy, the results suggest that maintaining access to the extracellular side of the protein after photoisomerization involves twisting of the chromophore and that the decisive switch in access to the cytoplasmic side results from relaxation of the chromophore when the constraints on the Schiff base are released by decomposition of the complex counterion.

Approaches to determining membrane protein structures to high resolution: do selections of subpopulations occur?

Three different methods are currently used for the study of high-resolution structures of membrane proteins: X-ray crystallography, electron crystallography, and nuclear magnetic resonance (NMR) spectroscopy. Thus far, all methods combined have yielded a rather modest number of crystal structures that have been solved at the atomic level. It is hypothesized here that different methods may select different populations of proteins on the basis of various properties. Thus, protein stability may be a significant factor in the formation of three-dimensional (3D) crystals from detergent solutions, since exposure of hydrophobic protein zones to water may cause structural perturbation or denaturation in conformationally labile proteins. This is different in the formation of two-dimensional (2D) crystals where a protein remains protected in its native membrane environment. A biological selection mechanism may therefore be operative in that highly ordered lattices may form only if strong protein-protein interactions are relevant in vivo, thereby limiting the number of proteins that are amenable to electron crystallography. Keeping a protein in a bilayer environment throughout 3D crystallization maintains the lateral pressure existing in native membranes. This can be accomplished by using lipidic cubic phases. Alternatively, the hydrophobic interface of a membrane protein may be spared from contact with water by crystallization from organic solvents where the polar caps are protected in reverse micelles by using appropriate detergents. Some of the criteria that are useful in optimizing the various approaches are given. While the usefulness of complementary methods seems obvious, the results presented may be particularly critical in recognizing key problems in other structural approaches.

Bioenergetics of the Archaea

In the late 1970s, on the basis of rRNA phylogeny, Archaea (archaebacteria) was identified as a distinct domain of life besides Bacteria (eubacteria) and Eucarya. Though forming a separate domain, Archaea display an enormous diversity of lifestyles and metabolic capabilities. Many archaeal species are adapted to extreme environments with respect to salinity, temperatures around the boiling point of water, and/or extremely alkaline or acidic pH. This has posed the challenge of studying the molecular and mechanistic bases on which these organisms can cope with such adverse conditions. This review considers our cumulative knowledge on archaeal mechanisms of primary energy conservation, in relationship to those of bacteria and eucarya. Although the universal principle of chemiosmotic energy conservation also holds for Archaea, distinct features have been discovered with respect to novel ion-transducing, membrane-residing protein complexes and the use of novel cofactors in bioenergetics of methanogenesis. From aerobically respiring Archaea, unusual electron-transporting supercomplexes could be isolated and functionally resolved, and a proposal on the organization of archaeal electron transport chains has been presented. The unique functions of archaeal rhodopsins as sensory systems and as proton or chloride pumps have been elucidated on the basis of recent structural information on the atomic scale. Whereas components of methanogenesis and of phototrophic energy transduction in halobacteria appear to be unique to Archaea, respiratory complexes and the ATP synthase exhibit some chimeric features with respect to their evolutionary origin. Nevertheless, archaeal ATP synthases are to be considered distinct members of this family of secondary energy transducers. A major challenge to future investigations is the development of archaeal genetic transformation systems, in order to gain access to the regulation of bioenergetic systems and to overproducers of archaeal membrane proteins as a prerequisite for their crystallization.

Magic angle spinning NMR spectroscopy of membrane proteins

The passage of molecules and information across cell membranes is mediated largely by membrane-spanning proteins acting as channels, pumps, receptors and enzymes. These proteins perform many tasks: they control electrochemical gradients across the membrane, receive signals from the environment or from other cells, convert light energy into chemical signals, transport small molecules into and out of cells, and harness proton gradients to generate the energy consumed in metabolism. Indeed, of the estimated 50000–100000 genes in the human genome, fully 20–40 % are thought to encode integral membrane proteins. If one also includes membrane-associated proteins, which are attached to the membrane surface through fatty acyl chains or electrostatic interactions, this percentage is likely to be much higher.

Reversal of the surface charge asymmetry in purple membrane due to single amino acid substitutions

Twenty-seven mutant bacteriorhodopsin's were screened to determine the PKa for reversal of the permanent electric dipole moment. The photoelectric response of an aqueous purple-membrane suspension was used to determine the direction of the purple-membrane dipole moment as a function of pH. The pK(a) for the dipole reversal of wild-type bacteriorhodopsin is 4.5. Six of the 27 mutant bacteriorhodopsin's were found to have a pK(a) for dipole reversal larger than that of wild-type bacteriorhodopsin. Two of these mutants, L93T and L93W, involve a neutral amino acid substitution in the interior of the protein. The direction of the purple-membrane permanent electric dipole moment is determined by the purple-membrane surface charge asymmetry. We conclude that these two substitutions, which do not involve charge replacement, alter the pK(a) for the reversal of the purple-membrane surface charge asymmetry. We suggest that these changes to the pK(a) are due to altered protein folding at the surface of the purple-membrane induced by single-site substitutions in the protein interior.

Molecular dynamics study of early picosecond events in the bacteriorhodopsin photocycle: dielectric response, vibrational cooling and the J, K intermediates

Molecular dynamics simulations have been carried out to study the J625 and K590 intermediates of bacteriorhodopsin's (bRs) photocycle starting from a refined structure of bR568. The coupling between the electronic states of retinal and the protein matrix is characterized by the energy difference delta E(t) between the excited state and the ground state to which the protein contributes through the Coulomb interaction. Our simulations indicate that the J625 intermediate is related to a polarization of the protein matrix due to the brief (200 fs) change of retinal's charge distribution in going to the excited state and back to the ground state, and that the rise time of the K590 intermediate is determined by vibrational cooling of retinal.

Molecular dynamics study of the M412 intermediate of bacteriorhodopsin

Molecular dynamics simulations have been carried out to study the M412 intermediate of bacteriorhodopsin's (bR) photocycle. The simulations start from two simulated structures for the L550 intermediate of the photocycle, one involving a 13-cis retinal with strong torsions, the other a 13,14-dicis retinal, from which the M412 intermediate is initiated through proton transfer to Asp-85. The simulations are based on a refined structure of bR568 obtained through all-atom molecular dynamics simulations and placement of 16 waters inside the protein. The structures of the L550 intermediates were obtained through simulated photoisomerization and subsequent molecular dynamics, and simulated annealing. Our simulations reveal that the M412 intermediate actually comprises a series of conformations involving 1) a motion of retinal; 2) protein conformational changes; and 3) diffusion and reconfiguration of water in the space between the retinal Schiff base nitrogen and the Asp-96 side group. (1) turns the retinal Schiff base nitrogen from an early orientation toward Asp-85 to a late orientation toward Asp-96; (2) disconnects the hydrogen bond network between retinal and Asp-85 and tilts the helix F of bR, enlarging bR's cytoplasmic channel; (3) adds two water molecules to the three water molecules existing in the cytoplasmic channel at the bR568 stage and forms a proton conduction pathway. The conformational change (2) of the protein involves a 60 degrees bent of the cytoplasmic side of helix F and is induced through a break of a hydrogen bond between Tyr-185 and a water-side group complex in the counterion region.

Membrane protein structure: the contribution and potential of novel solid state NMR approaches

Alternative methods for describing molecular detail for large integral membrane proteins are required in the absence of routine crystallographic approaches. Novel solid state NMR methods, devised for the study of large molecular assemblies, are now finding applications in biological systems, including integral membrane proteins. Wild-type and genetically engineered proteins can be investigated and detailed information about side chains, prosthetic groups, ligands (e.g. drugs) and binding sites can be deduced. The molecular structure and dynamics of selected parts of the proteins are accessible by a range of different solid state NMR approaches. Inter- and intra-atomic distances can be determined rather accurately (within ångströms) and the orientation of molecular bonds (within 2 degrees) can be measured in ideal cases. Here, a brief description of the methods is given and then some specific examples described with an indication of the future potential for the approaches in studying membrane proteins. It is anticipated that this emerging NMR methodology will be more widely used in the future, not only for resolving local structure, but also for more expansive descriptions of membrane protein structure at atomic resolution.

Molecular dynamics study of the 13-cis form (bR548) of bacteriorhodopsin and its photocycle

The structure and the photocycle of bacteriorhodopsin (bR) containing 13-cis,15-syn retinal, so-called bR548, has been studied by means of molecular dynamics simulations performed on the complete protein. The simulated structure of bR548 was obtained through isomerization of in situ retinal around both its C13-C14 and its C15-N bond starting from the simulated structure of bR568 described previously, containing all-trans,15-anti retinal. After a 50-ps equilibration, the resulting structure of bR548 was examined by replacing retinal by analogues with modified beta-ionone rings and comparing with respective observations. The photocycle of bR548 was simulated by inducing a rapid 13-cis,15-anti-->all-trans,15-syn isomerization through a 1-ps application of a potential that destabilizes the 13-cis isomer. The simulation resulted in structures consistent with the J, K, and L intermediates observed in the photocycle of bR548. The results offer an explanation of why an unprotonated retinal Schiff base intermediate, i.e., an M state, is not formed in the bR548 photocycle. The Schiff base nitrogen after photoisomerization of bR548 points to the intracellular rather than to the extracellular site. The simulations suggest also that leakage from the bR548 to the bR568 cycle arises due to an initial 13-cis,15-anti-->all-trans,15-anti photoisomerization.

Two-dimensional crystallization of Escherichia coli-expressed bacteriorhodopsin and its D96N variant: high resolution structural studies in projection

Highly ordered two-dimensional (2-D) crystals of Escherichia coli-expressed bacteriorhodopsin analog (e-bR) and its D96N variant (e-D96N) reconstituted in Halobacterium halobium lipids have been obtained by starting with the opsin protein purified in the denaturing detergent sodium dodecyl sulfate. These crystals embedded in glucose show electron diffraction in projection to better than 3.0 A at room temperature. This is the first instance that expressed bR or a variant has been crystallized in 2-D arrays showing such high order. The crystal lattice is homologous to that in wild-type bR (w-bR) in purple membranes (PM) and permit high resolution analyses of the structure of the functionally impaired D96N variant. The e-bR crystal is isomorphous to that in PM with an overall averaged fractional change of 12.7% (26-3.6-A resolution) in the projection structure factors. The projection difference Fourier map e-bR-PM at 3.6-A resolution indicates small conformational changes equivalent to movement of approximately < 7 C-atoms distributed within and in the neighborhood of the protein envelope. This result shows that relative to w-bR there are no global structural rearrangements in e-bR at this 3.6 A resolution level. The e-D96N crystal is isomorphous to the e-bR crystal with a smaller (9.2%) overall averaged fractional change in the structure factors. The significant structural differences between e-D96N and e-bR are concentrated at high resolution (5-3.6 A); however, these changes are small as quantified from the 3.6 A resolution e-D96N-e-bR Fourier difference map. The difference map showed no statistically significant peaks or valleys within 5 A in projection from the site of D96 substitution on helix C. Elsewhere within the protein envelope the integrated measure of peaks or valleys was < approximately 3 C-atom equivalents. Thus, our results show that for the isosteric substitution of Asp96 by Asn, the molecular conformation of bR in its ground state is essentially unaltered. Therefore, the known effect of D96N on the slowed M412 decay is not due to ground-state structural perturbations.

Molecular dynamics study of the proton pump cycle of bacteriorhodopsin

Retinal isomerization reactions, which are functionally important in the proton pump cycle of bacteriorhodopsin, were studied by molecular dynamics simulations performed on the complete protein. Retinal isomerizations were simulated in situ to account for the effects of the retinal-protein interactions. The protein structure employed was that described in Nonella et al. [Nonella, M., Windemuth, A., & Schulten, K. (1991) Photochem. Photobiol. 54, 937-948]. We investigated two mechanisms suggested previously for the proton pump cycle, the 13-cis isomerization model (C-T model) and the 13,14-dicis isomerization model. According to these models, retinal undergoes an all-trans-->13-cis or an all-trans-->13,14-dicis photoisomerization as the primary step of the pump cycle. From the simulations emerged a consistent picture of isomerization reactions and their control through the retinal-protein interactions which favors the 13,14-dicis isomerization model. Electrostatic interactions between the protonated Schiff base and its counterion are found to direct the stereochemistry of retinal in the photocycle: this and other interactions steer retinal toward the 13,14-dicis geometry in the primary photoreaction, toward the 13-cis geometry after its deprotonation, and to the all-trans isomeric form after its reprotonation. We also propose a catalytic mechanism involving hydrogen bonding of the Schiff base to main chain oxygen atoms of Val-49 and Thr-89 for the 13-cis-->all-trans thermal reisomerization of retinal. The all-trans-->13-cis primary photoreaction required by the "C-T" model was found to be inhibited by the Schiff base-counterion interaction, but the possibility of such a reaction can not be excluded. In order to investigate the "C-T" model, we enforced an all-trans-->13-cis photoisomerization in a simulation and monitored the subsequent protein conformational changes. The effects of internal water molecules on retinal isomerization reactions were studied by placing 16 water molecules in the proton conduction channel. The results indicate that water affects the nature of the Schiff base counterion and the nature of the primary photoreaction. Water chains, formed between positively and negatively charged protein groups in the proton conduction channel, are suggested to be involved in the reprotonation and deprotonation of retinal.

Charge transport of ion pumps on lipid bilayer membranes

Ion pumps create ion gradients across cell membranes while consuming light energy or chemical energy. The ion gradients are used by the corresponding cell types for passive-ion transport via ion channels or carriers or for accumulation of nutrients like sugar or amino acids via cotransport systems or antiporters.

Asp85 is the only internal aspartic acid that gets protonated in the M intermediate and the purple-to-blue transition of bacteriorhodopsin. A solid-state 13C CP-MAS NMR investigation

High-resolution solid-state 13C NMR spectra of the ground state and M intermediate of the bacteriorhodopsin mutant D96N with the isotope label at [4-13C]Asp and [11-13C]Trp were recorded. The NMR spectra show that Asp85 is protonated in the M intermediate. The environment of Asp85 is quite hydrophobic. On the other hand, Asp212 remains deprotonated and a slight shift to lower field indicates a more hydrophilic environment. Asp85 also protonates in the purple-to-blue transition of bacteriorhodopsin in the deionized membrane, where it experiences a similar environment to M. The shift of Trp resonances in M reflect a conformational change of the protein in forming the M intermediate.

NMR studies of retinal proteins

A review is given of the use of nuclear magnetic resonance (NMR) spectroscopy to study bacteriorhodopsin and bovine rhodopsin. Solution and solid-state approaches are included. The studies of the bacterial proton pump examine the chromophore, the peptide backbone, and the protein side chains. The studies of the bovine visual pigment are limited to the chromophore. Various forms of each pigment are considered. Both structural and dynamic features are addressed.

High-resolution solid state 13C NMR of bacteriorhodopsin: characterization of [4-13C]Asp resonances

Solid state 13C nuclear magnetic resonance measurements of bacteriorhodopsin labeled with [4-13C]Asp show that resonances of single amino acids can be resolved. In order to assign and characterize the resonances of specific Asp residues, three different approaches were used. (1) Determination of the chemical shift anisotropy from side-band intensities provides information about the protonation state of Asp residues. (2) Relaxation studies and T1 filtering allow one to discriminate between resonances with different mobility. (3) A comparison of the spectra of light- and dark-adapted bacteriorhodopsin provides evidence for resonances from aspartic acid residues in close neighborhood of the chromophore. In agreement with other investigations, four resonances are assigned to internal residues. Two of them are protonated in the ground state up to pH 10 (Asp96 and Asp115). All other detected resonances, including Asp85 and Asp212, are due to deprotonated aspartic acid. Two lines due to the two internal deprotonated groups change upon dark and light adaptation, whereas the protonated Asp residues are unaffected.