Protein-induced membrane curvature investigated through molecular dynamics flexible fitting

Direct visualization of exciton reequilibration in the LH1 and LH2 complexes of Rhodobacter sphaeroides by multipulse spectroscopy

The dynamics of the excited states of the light-harvesting complexes LH1 and LH2 of Rhodobacter sphaeroides are governed, mainly, by the excitonic nature of these ring-systems. In a pump-dump-probe experiment, the first pulse promotes LH1 or LH2 to its excited state and the second pulse dumps a portion of the excited state. By selective dumping, we can disentangle the dynamics normally hidden in the excited-state manifold. We find that by using this multiple-excitation technique we can visualize a 400-fs reequilibration reflecting relaxation between the two lowest exciton states that cannot be directly explored by conventional pump-probe. An oscillatory feature is observed within the exciton reequilibration, which is attributed to a coherent motion of a vibrational wavepacket with a period of ∼150 fs. Our disordered exciton model allows a quantitative interpretation of the observed reequilibration processes occurring in these antennas.

Förster energy transfer theory as reflected in the structures of photosynthetic light-harvesting systems

Förster’s theory of resonant energy transfer underlies a fundamental process in nature, namely the harvesting of sunlight by photosynthetic life forms. The theoretical framework developed by Förster and others describes how electronic excitation migrates in the photosynthetic apparatus of plants, algae, and bacteria from light absorbing pigments to reaction centers where light energy is utilized for the eventual conversion into chemical energy. The demand for highest possible efficiency of light harvesting appears to have shaped the evolution of photosynthetic species from bacteria to plants which, despite a great variation in architecture, display common structural themes founded on the quantum physics of energy transfer as described first by Förster. Herein, Förster’s theory of excitation transfer is summarized, including recent extensions, and the relevance of the theory to photosynthetic systems as evolved in purple bacteria, cyanobacteria, and plants is demonstrated. Förster’s energy transfer formula, as used widely today in many fields of science, is also derived.

Monomeric RC-LH1 core complexes retard LH2 assembly and intracytoplasmic membrane formation in PufX-minus mutants of Rhodobacter sphaeroides

In the model photosynthetic bacterium Rhodobacter sphaeroides domains of light-harvesting 2 (LH2) complexes surround and interconnect dimeric reaction centre-light-harvesting 1-PufX (RC-LH1-PufX) 'core' complexes, forming extensive networks for energy transfer and trapping. These complexes are housed in spherical intracytoplasmic membranes (ICMs), which are assembled in a stepwise process where biosynthesis of core complexes tends to dominate the early stages of membrane invagination. The kinetics of LH2 assembly were measured in PufX mutants that assemble monomeric core complexes, as a consequence of either a twelve-residue N-terminal truncation of PufX (PufXΔ12) or the complete removal of PufX (PufX(-)). Lower rates of LH2 assembly and retarded maturation of membrane invagination were observed for the larger and less curved ICM from the PufX(-) mutant, consistent with the proposition that local membrane curvature, initiated by arrays of bent RC-LH1-PufX dimers, creates a favourable environment for stable assembly of LH2 complexes. Transmission electron microscopy and high-resolution atomic force microscopy were used to examine ICM morphology and membrane protein organisation in these mutants. Some partitioning of core and LH2 complexes was observed in PufX(-) membranes, resulting in locally ordered clusters of monomeric RC-LH1 complexes. The distribution of core and LH2 complexes in the three types of membrane examined is consistent with previous models of membrane curvature and domain formation (Frese et al., 2008), which demonstrated that a combination of crowding and asymmetries in sizes and shapes of membrane protein complexes drives membrane organisation.

Experimental evidence that the membrane-spanning helix of PufX adopts a bent conformation that facilitates dimerisation of the Rhodobacter sphaeroides RC-LH1 complex through N-terminal interactions

The PufX polypeptide is an integral component of some photosynthetic bacterial reaction center-light harvesting 1 (RC-LH1) core complexes. Many aspects of the structure of PufX are unresolved, including the conformation of its long membrane-spanning helix and whether C-terminal processing occurs. In the present report, NMR data recorded on the Rhodobacter sphaeroides PufX in a detergent micelle confirmed previous conclusions derived from equivalent data obtained in organic solvent, that the α-helix of PufX adopts a bent conformation that would allow the entire helix to reside in the membrane interior or at its surface. In support of this, it was found through the use of site-directed mutagenesis that increasing the size of a conserved glycine on the inside of the bend in the helix was not tolerated. Possible consequences of this bent helical structure were explored using a series of N-terminal deletions. The N-terminal sequence ADKTIFNDHLN on the cytoplasmic face of the membrane was found to be critical for the formation of dimers of the RC-LH1 complex. It was further shown that the C-terminus of PufX is processed at an early stage in the development of the photosynthetic membrane. A model in which two bent PufX polypeptides stabilise a dimeric RC-LH1 complex is presented, and it is proposed that the N-terminus of PufX from one half of the dimer engages in electrostatic interactions with charged residues on the cytoplasmic surface of the LH1α and β polypeptides on the other half of the dimer.

Applications of the molecular dynamics flexible fitting method

In recent years, cryo-electron microscopy (cryo-EM) has established itself as a key method in structural biology, permitting the structural characterization of large biomolecular complexes in various functional states. The data obtained through single-particle cryo-EM has recently seen a leap in resolution thanks to landmark advances in experimental and computational techniques, resulting in sub-nanometer resolution structures being obtained routinely. The remaining gap between these data and revealing the mechanisms of molecular function can be closed through hybrid modeling tools that incorporate known atomic structures into the cryo-EM data. One such tool, molecular dynamics flexible fitting (MDFF), uses molecular dynamics simulations to combine structures from X-ray crystallography with cryo-EM density maps to derive atomic models of large biomolecular complexes. The structures furnished by MDFF can be used subsequently in computational investigations aimed at revealing the dynamics of the complexes under study. In the present work, recent applications of MDFF are presented, including the interpretation of cryo-EM data of the ribosome at different stages of translation and the structure of a membrane-curvature-inducing photosynthetic complex.

Recognition of the regulatory nascent chain TnaC by the ribosome

Regulatory nascent chains interact with the ribosomal exit tunnel and modulate their own translation. To characterize nascent chain recognition by the ribosome at the atomic level, extensive molecular dynamics simulations of TnaC, the leader peptide of the tryptophanase operon, inside the exit tunnel were performed for an aggregate time of 2.1 mus. The simulations, complemented by quantum chemistry calculations, suggest that the critical TnaC residue W12 is recognized by the ribosome via a cation-pi interaction, whereas TnaC's D16 forms salt bridges with ribosomal proteins. The simulations also show that TnaC-mediated translational arrest does not involve a swinging of ribosomal protein L22, as previously proposed. Furthermore, bioinformatic analyses and simulations suggest nascent chain elements that may prevent translational arrest in various organisms. Altogether, the current study unveils atomic-detail interactions that explain the role of elements of TnaC and the ribosome essential for translational arrest.

Native architecture of the photosynthetic membrane from Rhodobacter veldkampii

The photosynthetic membrane in purple bacteria contains several pigment-protein complexes that assure light capture and establishment of the chemiosmotic gradient. The bioenergetic tasks of the photosynthetic membrane require the strong interaction between these various complexes. In the present work, we acquired the first images of the native outer membrane architecture and the supramolecular organization of the photosynthetic apparatus in vesicular chromatophores of Rhodobacter (Rb.) veldkampii. Mixed with LH2 (light-harvesting complex 2) rings, the PufX-containing LH1-RC (light-harvesting complex 1--reaction center) core complexes appear as C-shaped monomers, with random orientations in the photosynthetic membrane. Within the LH1 fence surrounding the RC, a remarkable gap that is probably occupied (or partially occupied) by PufX is visualized. Sequence alignment revealed that one specific region in PufX may be essential for PufX-induced core dimerization. In this region of ten amino acids in length all Rhodobacter species had five conserved amino acids, with the exception of Rb. veldkampii. Our findings provide direct evidence that the presence of PufX in Rb. veldkampii does not directly govern the dimerization of LH1-RC core complexes in the native membrane. It is indicated, furthermore, that the high membrane curvature of Rb. veldkampii chromatophores (Rb. veldkampii features equally small vesicular chromatophores alike Rb. sphaeroides) is not due to membrane bending induced by dimeric RC-LH1-PufX cores, as it has been proposed in Rb. sphaeroides.

Energy Transfer Dynamics in an RC-LH1-PufX Tubular Photosynthetic Membrane

Light absorption and the subsequent transfer of excitation energy are the first two steps of the photosynthetic process, carried out by protein-bound pigments, mainly bacteriochlorophylls (BChls), in photosynthetic bacteria. BChls are anchored in light-harvesting (LH) complexes, such as light-harvesting complex I (LH1), which directly associates with the reaction center (RC), forming the RC-LH1 core complex. In Rhodobacter sphaeroides, RC-LH1 core complexes contain an additional protein, PufX, and assemble into dimeric RC-LH1-PufX core complexes. In the absence of light-harvesting complexes II, the former complexes can aggregate into a helically ordered tubular photosynthetic membrane. We examined the excitation transfer dynamics in a single RC-LH1-PufX core complex dimer using the hierarchical equations of motion for dissipative quantum dynamics that accurately, yet computationally costly, treat the coupling between BChls and their protein environment. A widely employed description, generalized Förster theory, was also used to calculate the transfer rates of the same excitonic system in order to verify the accuracy of this computationally cheap method. Additionally, in light of the structural uncertainties in the Rhodobacter sphaeroides RC-LH1-PufX core complex, geometrical alterations were introduced in the BChl organization. It is shown that the energy transfer dynamics is not affected by the considered changes in the BChl organization, and that generalized Förster theory provides accurate transfer rates. An all-atom model for a tubular photosynthetic membrane is then constructed on the basis of electron microscopy data, and the overall energy transfer properties of this membrane are computed.

Photosynthetic vesicle architecture and constraints on efficient energy harvesting

Photosynthetic chromatophore vesicles found in some purple bacteria constitute one of the simplest light-harvesting systems in nature. The overall architecture of chromatophore vesicles and the structural integration of vesicle function remain poorly understood despite structural information being available on individual constituent proteins. An all-atom structural model for an entire chromatophore vesicle is presented, which improves upon earlier models by taking into account the stoichiometry of core and antenna complexes determined by the absorption spectrum of intact vesicles in Rhodobacter sphaeroides, as well as the well-established curvature-inducing properties of the dimeric core complex. The absorption spectrum of low-light-adapted vesicles is shown to correspond to a light-harvesting-complex 2 to reaction center ratio of 3:1. A structural model for a vesicle consistent with this stoichiometry is developed and used in the computation of excitonic properties. Considered also is the packing density of antenna and core complexes that is high enough for efficient energy transfer and low enough for quinone diffusion from reaction centers to cytochrome bc(1) complexes.

Self-assembly of photosynthetic membranes

Bacterial photosynthetic membranes, also known as chromatophores, are tightly packed with integral membrane proteins that work together to carry out photosynthesis. Chromatophores display a wide range of cellular morphologies; spherical, tubular, and lamellar chromatophores have all been observed in different bacterial species, or with different protein constituents. Through recent computational modeling and simulation, it has been demonstrated that the light-harvesting complexes abundant in chromatophores induce local membrane curvature via multiple mechanisms. These protein complexes assemble to generate a global curvature and sculpt the chromatophores into various cellular-scale architectures.

Membrane invagination in Rhodobacter sphaeroides is initiated at curved regions of the cytoplasmic membrane, then forms both budded and fully detached spherical vesicles

The purple phototrophic bacteria synthesize an extensive system of intracytoplasmic membranes (ICM) in order to increase the surface area for absorbing and utilizing solar energy. Rhodobacter sphaeroides cells contain curved membrane invaginations. In order to study the biogenesis of ICM in this bacterium mature (ICM) and precursor (upper pigmented band - UPB) membranes were purified and compared at the single membrane level using electron, atomic force and fluorescence microscopy, revealing fundamental differences in their morphology, protein organization and function. Cryo-electron tomography demonstrates the complexity of the ICM of Rba. sphaeroides. Some ICM vesicles have no connection with other structures, others are found nearer to the cytoplasmic membrane (CM), often forming interconnected structures that retain a connection to the CM, and possibly having access to the periplasmic space. Near-spherical single invaginations are also observed, still attached to the CM by a 'neck'. Small indents of the CM are also seen, which are proposed to give rise to the UPB precursor membranes upon cell disruption. 'Free-living' ICM vesicles, which possess all the machinery for converting light energy into ATP, can be regarded as bacterial membrane organelles.

Regulation of the protein-conducting channel by a bound ribosome

During protein synthesis, it is often necessary for the ribosome to form a complex with a membrane-bound channel, the SecY/Sec61 complex, in order to translocate nascent proteins across a cellular membrane. Structural data on the ribosome-channel complex are currently limited to low-resolution cryo-electron microscopy maps, including one showing a bacterial ribosome bound to a monomeric SecY complex. Using that map along with available atomic-level models of the ribosome and SecY, we have determined, through molecular dynamics flexible fitting (MDFF), an atomic-resolution model of the ribosome-channel complex. We characterized computationally the sites of ribosome-SecY interaction within the complex and determined the effect of ribosome binding on the SecY channel. We also constructed a model of a ribosome in complex with a SecY dimer by adding a second copy of SecY to the MDFF-derived model. The study involved 2.7-million-atom simulations over altogether nearly 50 ns.

Membrane curvature induced by aggregates of LH2s and monomeric LH1s

The photosynthetic apparatus of purple bacteria is contained within organelles called chromatophores, which form as extensions of the cytoplasmic membrane. The shape of these chromatophores can be spherical (as in Rhodobacter sphaeroides), lamellar (as in Rhodopseudomonas acidophila and Phaeospirillum molischianum), or tubular (as in certain Rb. sphaeroides mutants). Chromatophore shape is thought to be influenced by the integral membrane proteins Light Harvesting Complexes I and II (LH1 and LH2), which pack tightly together in the chromatophore. It has been suggested that the shape of LH2, together with its close packing in the membrane, induces membrane curvature. The mechanism of LH2-induced curvature is explored via molecular dynamics simulations of multiple LH2 complexes in a membrane patch. LH2s from three species-Rb. sphaeroides, Rps. acidophila, and Phsp. molischianum-were simulated in different packing arrangements. In each case, the LH2s pack together and tilt with respect to neighboring LH2s in a way that produces an overall curvature. This curvature appears to be driven by a combination of LH2's shape and electrostatic forces that are modulated by the presence of well-conserved cytoplasmic charged residues, the removal of which inhibits LH2 curvature. The interaction of LH2s and an LH1 monomer is also explored, and it suggests that curvature is diminished by the presence of LH1 monomers. The implications of our results for chromatophore shape are discussed.

Structural model and excitonic properties of the dimeric RC-LH1-PufX complex from Rhodobacter sphaeroides

The light-harvesting apparatus of the purple bacterial photosynthetic unit consists of a pool of peripheral light-harvesting complexes that transfer excitation energy to a reaction center (RC) via the surrounding pigment-protein complex LH1. Recent electron microscopy and atomic force microscopy studies have revealed that RC-LH1 units of Rhodobacter sphaeroides form membrane-bending dimeric complexes together with the polypeptide PufX. We present a structural model for these RC-LH1-PufX dimeric complexes constructed using the molecular dynamics flexible fitting method based on an EM density map. The arrangement of the LH1 BChls displays a distortion near the proposed location of the PufX polypeptide. The resulting atomic model for BChl arrays is used to compute the excitonic properties of the dimeric RC-LH1 complex. A comparison is presented between the structural and excitonic features of the S-shaped dimeric BChl array of Rhodobacter sphaeroides and the circular BChl arrangement found in other purple bacteria.