Large-Scale Molecular Dynamics Simulations of Cellular Compartments.

Molecular dynamics or MD simulation is gradually maturing into a tool for constructing in vivo models of living cells in atomistic details. The feasibility of such models is bolstered by integrating the simulations with data from microscopic, tomographic and spectroscopic experiments on exascale supercomputers, facilitated by the use of deep learning technologies. Over time, MD simulation has evolved from tens of thousands of atoms to over 100 million atoms comprising an entire cell organelle, a photosynthetic chromatophore vesicle from a purple bacterium. In this chapter, we present a step-by-step outline for preparing, executing and analyzing such large-scale MD simulations of biological systems that are essential to life processes. All scripts are provided via GitHub.

Phospholipid distributions in purple phototrophic bacteria and LH1-RC core complexes.

In contrast to plants, algae and cyanobacteria that contain glycolipids as the major lipid components in their photosynthetic membranes, phospholipids are the dominant lipids in the membranes of anoxygenic purple phototrophic bacteria. Although the phospholipid compositions in whole cells or membranes are known for a limited number of the purple bacteria, little is known about the phospholipids associated with individual photosynthetic complexes. In this study, we investigated the phospholipid distributions in both membranes and the light-harvesting 1-reaction center (LH1-RC) complexes purified from several purple sulfur and nonsulfur bacteria. 31P NMR was used for determining the phospholipid compositions and inductively coupled plasma atomic emission spectroscopy was used for measuring the total phosphorous contents. Combining these two techniques, we could determine the numbers of specific phospholipids in the purified LH1-RC complexes. A total of approximate 20-30 phospholipids per LH1-RC were detected as the tightly bound lipids in all species. The results revealed that while cardiolipin (CL) exists as a minor component in the membranes, it became the most abundant phospholipid in the purified core complexes and the sum of CL and phosphatidylglycerol accounted for more than two thirds of the total phospholipids for most species. Preferential association of these anionic phospholipids with the LH1-RC is discussed in the context of the recent high-resolution structure of this complex from Thermochromatium (Tch.) tepidum. The detergent lauryldimethylamine N-oxide was demonstrated to selectively remove phosphatidylethanolamine from the membrane of Tch. tepidum.

Probing the local lipid environment of the Rhodobacter sphaeroides cytochrome bc1 and Synechocystis sp. PCC 6803 cytochrome b6f complexes with styrene maleic acid.

Intracytoplasmic vesicles (chromatophores) in the photosynthetic bacterium Rhodobacter sphaeroides represent a minimal structural and functional unit for absorbing photons and utilising their energy for the generation of ATP. The cytochrome bc1 complex (cytbc1) is one of the four major components of the chromatophore alongside the reaction centre-light harvesting 1-PufX core complex (RC-LH1-PufX), the light-harvesting 2 complex (LH2), and ATP synthase. Although the membrane organisation of these complexes is known, their local lipid environments have not been investigated. Here we utilise poly(styrene-alt-maleic acid) (SMA) co-polymers as a tool to simultaneously determine the local lipid environments of the RC-LH1-PufX, LH2 and cytbc1 complexes. SMA has previously been reported to effectively solubilise complexes in lipid-rich membrane regions whilst leaving lipid-poor ordered protein arrays intact. Here we show that SMA solubilises cytbc1 complexes with an efficiency of nearly 70%, whereas solubilisation of RC-LH1-PufX and LH2 was only 10% and 22% respectively. This high susceptibility of cytbc1 to SMA solubilisation is consistent with this complex residing in a locally lipid-rich region. SMA solubilised cytbc1 complexes retain their native dimeric structure and co-purify with 56±6 phospholipids from the chromatophore membrane. We extended this approach to the model cyanobacterium Synechocystis sp. PCC 6803, and show that the cytochrome b6f complex (cytb6f) and Photosystem II (PSII) complexes are susceptible to SMA solubilisation, suggesting they also reside in lipid-rich environments. Thus, lipid-rich membrane regions could be a general requirement for cytbc1/cytb6f complexes, providing a favourable local solvent to promote rapid quinol/quinone binding and release at the Q0 and Qi sites.

Determination of Cell Doubling Times from the Return-on-Investment Time of Photosynthetic Vesicles Based on Atomic Detail Structural Models.

Cell doubling times of the purple bacterium Rhodobacter sphaeroides during photosynthetic growth are determined experimentally and computationally as a function of illumination. For this purpose, energy conversion processes in an intracytoplasmic membrane vesicle, the chromatophore, are described based on an atomic detail structural model. The cell doubling time and its illumination dependence are computed in terms of the return-on-investment (ROI) time of the chromatophore, determined computationally from the ATP production rate, and the mass ratio of chromatophores in the cell, determined experimentally from whole cell absorbance spectra. The ROI time is defined as the time it takes to produce enough ATP to pay for the construction of another chromatophore. The ROI time of the low light-growth chromatophore is 4.5-2.6 h for a typical illumination range of 10-100 µmol photons m-2 s-1, respectively, with corresponding cell doubling times of 8.2-3.9 h. When energy expenditure is considered as a currency, the benefit-to-cost ratio computed for the chromatophore as an energy harvesting device is 2-8 times greater than for photovoltaic and fossil fuel-based energy solutions and the corresponding ROI times are approximately 3-4 orders of magnitude shorter for the chromatophore than for synthetic systems.

Computational Methodologies for Real-Space Structural Refinement of Large Macromolecular Complexes.

The rise of the computer as a powerful tool for model building and refinement has revolutionized the field of structure determination for large biomolecular systems. Despite the wide availability of robust experimental methods capable of resolving structural details across a range of spatiotemporal resolutions, computational hybrid methods have the unique ability to integrate the diverse data from multimodal techniques such as X-ray crystallography and electron microscopy into consistent, fully atomistic structures. Here, commonly employed strategies for computational real-space structural refinement are reviewed, and their specific applications are illustrated for several large macromolecular complexes: ribosome, virus capsids, chemosensory array, and photosynthetic chromatophore. The increasingly important role of computational methods in large-scale structural refinement, along with current and future challenges, is discussed.

Light harvesting by lamellar chromatophores in Rhodospirillum photometricum.

Purple photosynthetic bacteria harvest light using pigment-protein complexes which are often arranged in pseudo-organelles called chromatophores. A model of a chromatophore from Rhodospirillum photometricum was constructed based on atomic force microscopy data. Molecular-dynamics simulations and quantum-dynamics calculations were performed to characterize the intercomplex excitation transfer network and explore the interplay between close-packing and light-harvesting efficiency.

Functional interfacing of Rhodospirillum rubrum chromatophores to a conducting support for capture and conversion of solar energy.

Owing to the considerable current interest in replacing fossil fuels with solar radiation as a clean, renewable, and secure energy source, light-driven electron transport in natural photosynthetic systems offers a valuable blueprint for conversion of sunlight to useful energy forms. In particular, intracytoplasmic membrane vesicles (chromatophores) from the purple bacterium Rhodospirillum rubrum provide a fully functional and robust photosynthetic apparatus, ideal for biophysical investigations of energy transduction and incorporation into biohybrid photoelectrochemical devices. These vesicular organelles, which arise by invagination of the cytoplasmic membrane, are the sites of the photochemical reaction centers and the light harvesting 1 (LH1) complex. The LH1 protein is responsible for collecting visible and near-IR radiant energy and funneling these excitations to the reaction center for conversion into a transmembrane charge separation. Here, we have investigated the morphology, fluorescence kinetics and photocurrent generation of chromatophores from Rsp. rubrum deposited directly onto gold surfaces in the absence of chemical surface modifications. Atomic force microscopy showed a significant coverage of the gold electrode surface by Rsp. rubrum chromatophores. By in situ fluorescence induction/relaxation measurements, a high retention of the quantum yield of photochemistry was demonstrated in the photoactive films. Chronoamperometric measurements showed that the assembled bioelectrodes were capable of generating sustained photocurrent under white light illumination at 220 mW/cm(2) with a maximum current of 1.5 µA/cm(2), which slowly declines in about 1 week. This study demonstrates the possibility of photoelectrochemical control of robust chromatophore preparations from Rsp. rubrum that paves the way for future incorporation into functional solar cells.

Near-IR absorbing solar cell sensitized with bacterial photosynthetic membranes.

Current interest in natural photosynthesis as a blueprint for solar energy conversion has led to the development of a biohybrid photovoltaic cell in which bacterial photosynthetic membrane vesicles (chromatophores) have been adsorbed to a gold electrode surface in conjunction with biological electrolytes (quinone [Q] and cytochrome c; Magis et al. [2010] Biochim. Biophys. Acta 1798, 637-645). Since light-driven current generation was dependent on an open circuit potential, we have tested whether this external potential could be replaced in an appropriately designed dye-sensitized solar cell (DSSC). Herein, we show that a DSSC system in which the organic light-harvesting dye is replaced by robust chromatophores from Rhodospirillum rubrum, together with Q and cytochrome c as electrolytes, provides band energies between consecutive interfaces that facilitate a unidirectional flow of electrons. Solar I-V testing revealed a relatively high I(sc) (short-circuit current) of 25 µA cm(-2) and the cell was capable of generating a current utilizing abundant near-IR photons (maximum at ca 880 nm) with greater than eight-fold higher energy conversion efficiency than white light. These studies represent a powerful demonstration of the photoexcitation properties of a biological system in a closed solid-state device and its successful implementation in a functioning solar cell.

Thermostability of Rhodopseudomonas viridis and Rhodospirillum rubrum chromatophores reflecting physiological conditions.

Relationships between growth conditions and thermostability were examined for photosynthetic inner membranes (chromatophores) from Rhodopseudomonas viridis and Rhodospirillum rubrum of which morphology, lipid composition, and protein/lipid rate are rather mutually different. Signals observed by differential scanning calorimetry of the chromatophores were correlated with thermal state transitions of the membrane components by reference to temperature dependencies of circular dichroism and absorption spectra of the purified supramolecule comprising a photoreaction center and surrounding light-harvesting pigment-protein complexes that are the prominent proteins in both membranes. The differential scanning calorimetry curves of those chromatophores exhibited different dependencies on growth stages and environmental temperatures. The obtained result appeared to reflect the differences in the protein/lipid rate and protein-lipid specificity between the two chromatophores.

Identification of chromatophore membrane protein complexes formed under different nitrogen availability conditions in Rhodospirillum rubrum.

The chromatophore membrane of the photosynthetic diazotroph Rhodospirillum rubrum is of vital importance for a number of central processes, including nitrogen fixation. Using a novel amphiphile, we have identified protein complexes present under different nitrogen availability conditions by the use of two-dimensional Blue Native/SDS-PAGE and NSI-LC-LTQ-Orbitrap mass spectrometry. We have identified several membrane protein complexes, including components of the ATP synthase, reaction center, light harvesting, and NADH dehydrogenase complexes. Additionally, we have identified differentially expressed proteins, such as subunits of the succinate dehydrogenase complex and other TCA cycle enzymes that are usually found in the cytosol, thus hinting at a possible association to the membrane in response to nitrogen deficiency. We propose a redox sensing mechanism that can influence the membrane subproteome in response to nitrogen availability.

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.

Study of effect of molecular mobility in chromatophore membranes of the bacterium E. shaposhnikovii on processes of photoinduced electron transport using the NMR-spin-echo method with isotope substitution and dehydration.

The effect of dehydration and (2)H2O/H2O isotope substitution on electron transport reactions and relaxation of proton-containing groups was studied in chromatophore membranes of Ectothiorhodospira shaposhnikovii. During dehydration (including isotope substitution of hydrate water) of preliminarily dehydrated isolated photosynthetic membranes there was a partial correlation between hydration intervals within which activation of electron transport from high-potential cytochrome c to photoactive bacteriochlorophyll dimer P890 of photosynthetic reaction center and variation of spin-lattice and spin-spin proton relaxation time was observed. Partial correlation between hydration intervals can be considered as evidence of correlation between mobility of non-water proton-containing groups with proton relaxation frequency approximately 10(8) sec(-1) with efficiency of electron transfer at the donor side of the chain.

Intrinsic curvature properties of photosynthetic proteins in chromatophores.

In purple bacteria, photosynthesis is carried out on large indentations of the bacterial plasma membrane termed chromatophores. Acting as primitive organelles, chromatophores are densely packed with the membrane proteins necessary for photosynthesis, including light harvesting complexes LH1 and LH2, reaction center (RC), and cytochrome bc(1). The shape of chromatophores is primarily dependent on species, and is typically spherical or flat. How these shapes arise from the protein-protein and protein-membrane interactions is still unknown. Now, using molecular dynamics simulations, we have observed the dynamic curvature of membranes caused by proteins in the chromatophore. A membrane-embedded array of LH2s was found to relax to a curved state, both for LH2 from Rps. acidophila and a homology-modeled LH2 from Rb. sphaeroides. A modeled LH1-RC-PufX dimer was found to develop a bend at the dimerizing interface resulting in a curved shape as well. In contrast, the bc(1) complex, which has not been imaged yet in native chromatophores, did not induce a preferred membrane curvature in simulation. Based on these results, a model for how the different photosynthetic proteins influence chromatophore shape is presented.

Heterogeneity of photosynthetic membranes from Rhodobacter capsulatus: size dispersion and ATP synthase distribution.

The density distribution of photosynthetic membrane vesicles (chromatophores) from Rhodobacter capsulatus has been studied by isopicnic centrifugation. The average vesicle diameters, examined by electron microscopy, varied between 61 and 72 nm in different density fractions (70 nm in unfractionated chromatophores). The ATP synthase catalytic activities showed maxima displaced toward the higher density fractions relative to bacteriochlorophyll, resulting in higher specific activities in those fractions (about threefold). The amount of ATP synthase, measured by quantitative Western blotting, paralleled the catalytic activities. The average number of ATP synthases per chromatophore, evaluated on the basis of the Western blotting data and of vesicle density analysis, ranged between 8 and 13 (10 in unfractionated chromatophores). Poisson distribution analysis indicated that the probability of chromatophores devoid of ATP synthase was negligible. The effects of ATP synthase inhibition by efrapeptin on the time course of the transmembrane electric potential (evaluated as carotenoid electrochromic response) and on ATP synthesis were studied comparatively. The ATP produced after a flash and the total charge associated with the proton flow coupled to ATP synthesis were more resistant to efrapeptin than the initial value of the phosphorylating currents, indicating that several ATP synthases are fed by protons from the same vesicle.

Atomic-level structural and functional model of a bacterial photosynthetic membrane vesicle.

The photosynthetic unit (PSU) of purple photosynthetic bacteria consists of a network of bacteriochlorophyll-protein complexes that absorb solar energy for eventual conversion to ATP. Because of its remarkable simplicity, the PSU can serve as a prototype for studies of cellular organelles. In the purple bacterium Rhodobacter sphaeroides the PSU forms spherical invaginations of the inner membrane, approximately 70 nm in diameter, composed mostly of light-harvesting complexes, LH1 and LH2, and reaction centers (RCs). Atomic force microscopy studies of the intracytoplasmic membrane have revealed the overall spatial organization of the PSU. In the present study these atomic force microscopy data were used to construct three-dimensional models of an entire membrane vesicle at the atomic level by using the known structure of the LH2 complex and a structural model of the dimeric RC-LH1 complex. Two models depict vesicles consisting of 9 or 18 dimeric RC-LH1 complexes and 144 or 101 LH2 complexes, representing a total of 3,879 or 4,464 bacteriochlorophylls, respectively. The in silico reconstructions permit a detailed description of light absorption and electronic excitation migration, including computation of a 50-ps excitation lifetime and a 95% quantum efficiency for one of the model membranes, and demonstration of excitation sharing within the closely packed RC-LH1 dimer arrays.

Isotopic replacement of pigments and a lipid in chlorosomes from Chlorobium limicola: characterization of the resultant chlorosomes.

Pigments including bacteriochlorophyll (BChl) c, carotenoids, and a trace of BChl a together with a lipid, monogalactosyl diglyceride (MGDG), were extracted with chloroform/methanol (1:1 v/v) from an aqueous suspension (50 mM Tris-HCl, pH 8.0) of chlorosomes from Chlorobium limicola; other lipids and proteins were left behind in the aqueous layer by funnel separation. The chloroform layer was dried by purging N2 gas, dissolved in methanol, and rapidly injected into the aqueous layer to reassemble chlorosomes. This technique has been developed to replace one-half of the inherent 12C-BChl c by 13C-BChl c to identify the intermolecular 13C...13C magnetic dipole correlation peaks (that are supposed to reduce their intensities to one-fourth by reducing the 13C-BChl c concentration into one-half) and to determine the structure of BChl c aggregates in the rod elements by means of solid-state NMR spectroscopy. The isotopically replaced chlorosomes were characterized (1) by sucrose density gradient centrifugation, zeta potential measurement, electron microscopy, and dynamic light scattering measurement to determine the morphology of chlorosomes, (2) by 13C NMR spectroscopy, electronic absorption and circular dichroism spectroscopies, and low-angle X-ray diffraction to determine the pigment assembly in the rod elements, and (3) by subpicosecond time-resolved absorption spectroscopy to determine the excited-state dynamics in the pigment assembly. The results characterized the reassembled chlorosomes to have (1) similar but longer morphological structures, (2) almost the same pigment assembly in the rod elements, and (3) basically the same excited-state dynamics in the pigment assembly.

The 17-propionate function of (bacterio)chlorophylls: biological implication of their long esterifying chains in photosynthetic systems.

Molecular structures of (bacterio)chlorophylls [= (B)Chls] in photosynthetic apparatus are surveyed, and a diversity of the ester groups of the 17-propionate substituent is particularly focused on in this review. In oxygenic photosynthetic species including green plants and algae, the ester of Chl molecules is limited to a phytyl group. Geranylgeranyl and farnesyl groups in addition to phytyl are observed in (B)Chl molecules inside photosynthetic proteins of anoxygenic bacteria. In main light-harvesting antennas of green bacteria (chlorosomes), a greater variety of ester groups including long straight chains are used in the composite BChl molecules. This diversity is ascribable to the fact that chlorosomal BChls self-aggregate to form a core part of chlorosomes without any specific interaction of oligopeptides. Biological significance of the long chains is discussed in photosynthetic apparatus, especially in chlorosomes.

Transmission electron microscopic study on supramolecular nanostructures of bacteriochlorophyll self-aggregates in chlorosomes of green photosynthetic bacteria.

Supramolecular nanostructures of bacteriochlorophyll (BChl) self-aggregates in major light-harvesting complexes (chlorosomes) of green photosynthetic bacteria were successfully observed by freeze-fracture transmission electron microscope. Rod-shaped nanostructures with approximately 10 nm in diameter could be visualized in three types of green sulfur bacteria (Chlorobium). Diameter of the rod-shaped nanostructures in Chlorobium chlorosomes was independent of the molecular structures of their light-harvesting pigments, namely BChl-c or d. In contrast, chlorosomes of the green filamentous bacterium Chloroflexus aurantiacus had rod-shaped nanostructures with approximately 5 nm in diameter. The present results support that BChl self-aggregates in chlorosomes form rod-shaped nanostructures called rod-elements with approximately 10- and 5-nm diameters for Chlorobium and Chloroflexus, respectively.

AFM studies of the supramolecular assembly of bacterial photosynthetic core-complexes.

The atomic force microscope (AFM) allows visualization of the assembly and molecular interactions of single proteins. Most recently, AFM images of bacterial membranes have revealed details of the supramolecular architecture of bacterial photosynthetic apparatus in different species. The near-native experimental conditions used in AFM imaging reduce artefacts and make AFM ideal for studying native conformations. High-resolution AFM of native membranes has revealed variation in core-complex architectures amongst species.

Structure, function, and wavelength selection in blue-absorbing proteorhodopsin.

The absorption maximum of blue proteorhodopsin (BPR) is the most blue-shifted of all retinal proteins found in archaea or bacteria, with the exception of sensory rhodopsin II (SRII). The absorption spectrum also exhibits a pH dependence larger than any other retinal protein. We examine the structural origins of these two properties of BPR by using optical spectroscopy, homology modeling, and molecular orbital theory. Bacteriorhodopsin (BR) and SRII are used as homology parents for comparative purposes. We find that the tertiary structure of BPR based on SRII is more realistic with respect to free energy, dynamic stability, and spectroscopic properties. Molecular orbital calculations including full single- and double-configuration interaction within the chromophore pi-electron system provide perspectives on the wavelength regulation in this protein and indicate that Arg-95, Gln-106, Glu-143, and Asp-229 play important, and in some cases pH-dependent roles. A possible model for the 0.22 eV red shift of BPR at low pH is examined, in which Glu-143 becomes protonated and releases Arg-95 to rotate up into the binding site, altering the electrostatic environment of the chromophore. At high pH, BPR has spectroscopic properties similar to SRII, but at low pH, BPR has spectroscopic properties more similar to BR. Nevertheless, SRII is a significantly better homology model for BPR and opens up the question of whether this protein serves as a proton pump, as commonly believed, or is a light sensor with structure-function properties more comparable to those of SRII. The function of BPR in the native organism is discussed with reference to the results of the homology model.