Protein-induced membrane curvature investigated through molecular dynamics flexible fitting

Sulfoquinovosyl diacylglycerol is required for dimerisation of the Rhodobacter sphaeroides reaction centre-light harvesting 1 core complex

The reaction centre-light harvesting 1 (RC-LH1) core complex is indispensable for anoxygenic photosynthesis. In the purple bacterium Rhodobacter (Rba.) sphaeroides RC-LH1 is produced both as a monomer, in which 14 LH1 subunits form a C-shaped antenna around 1 RC, and as a dimer, where 28 LH1 subunits form an S-shaped antenna surrounding 2 RCs. Alongside the five RC and LH1 subunits, an additional polypeptide known as PufX provides an interface for dimerisation and also prevents LH1 ring closure, introducing a channel for quinone exchange that is essential for photoheterotrophic growth. Structures of Rba. sphaeroides RC-LH1 complexes revealed several new components; protein-Y, which helps to form the quinone channel; protein-Z, of unknown function and seemingly unique to dimers; and a tightly bound sulfoquinovosyl diacylglycerol (SQDG) lipid that interacts with two PufX arginine residues. This lipid lies at the dimer interface alongside weak density for a second molecule, previously proposed to be an ornithine lipid. In this work we have generated strains of Rba. sphaeroides lacking protein-Y, protein-Z, SQDG or ornithine lipids to assess the roles of these previously unknown components in the assembly and activity of RC-LH1. We show that whilst the removal of either protein-Y, protein-Z or ornithine lipids has only subtle effects, SQDG is essential for the formation of RC-LH1 dimers but its absence has no functional effect on the monomeric complex.

Mechanism of Membrane Curvature Induced by SNX1: Insights from Molecular Dynamics Simulations

SNX proteins have been found to induce membrane remodeling to facilitate the generation of transport carriers in endosomal pathways. However, the molecular mechanism of membrane bending and the role of lipids in the bending process remain elusive. Here, we conducted coarse-grained molecular dynamics simulations to investigate the role of the three structural modules (PX, BAR, and AH) of SNX1 and the PI3P lipids in membrane deformation. We observed that the presence of all three domains is essential for SNX1 to achieve a stable membrane deformation. BAR is capable of remodeling the membrane through the charged residues on its concave surface, but it requires PX and AH to establish stable membrane binding. AH penetrates into the lipid membrane, thereby promoting the induction of membrane curvature; however, it is inadequate on its own to maintain membrane bending. PI3P lipids are also indispensable for membrane remodeling, as they play a dominant role in the interactions of lipids with the BAR domain. Our results enhance the comprehension of the molecular mechanism underlying SNX1-induced membrane curvature and help future studies of curvature-inducing proteins.

Relationship between Protein-Induced Membrane Curvature and Membrane Thermal Undulation

This work studied the membrane curvature generated by anchored proteins lacking amphipathic helices and intrinsic morphologies, including the Epsin N-terminal homology domain, intrinsically disordered C-terminal domain, and truncated C-terminal fragments, by using coarse-grained molecular dynamics simulations. We found that anchored proteins can stabilize the thermal undulation of membranes at a wavelength five times the protein's binding size. This proportional connection is governed by the membrane bending rigidity and protein density. Extended intrinsically disordered proteins with relatively high hydrophobicity favor colliding with the membrane, leading to a much larger binding size, and show superiority in generating membrane curvature at low density over folded proteins.

Molecular Modeling is an Enabling Approach to Complement and Enhance Channelopathy Research

Hundreds of human membrane proteins form channels that transport necessary ions and compounds, including drugs and metabolites, yet details of their normal function or how function is altered by genetic variants to cause diseases are often unknown. Without this knowledge, researchers are less equipped to develop approaches to diagnose and treat channelopathies. High-resolution computational approaches such as molecular modeling enable researchers to investigate channelopathy protein function, facilitate detailed hypothesis generation, and produce data that is difficult to gather experimentally. Molecular modeling can be tailored to each physiologic context that a protein may act within, some of which may currently be difficult or impossible to assay experimentally. Because many genomic variants are observed in channelopathy proteins from high-throughput sequencing studies, methods with mechanistic value are needed to interpret their effects. The eminent field of structural bioinformatics integrates techniques from multiple disciplines including molecular modeling, computational chemistry, biophysics, and biochemistry, to develop mechanistic hypotheses and enhance the information available for understanding function. Molecular modeling and simulation access 3D and time-dependent information, not currently predictable from sequence. Thus, molecular modeling is valuable for increasing the resolution with which the natural function of protein channels can be investigated, and for interpreting how genomic variants alter them to produce physiologic changes that manifest as channelopathies. © 2022 American Physiological Society. Compr Physiol 12:3141-3166, 2022.

Cryo-EM structure of the dimeric Rhodobacter sphaeroides RC-LH1 core complex at 2.9 Å: the structural basis for dimerisation

The dimeric reaction centre light-harvesting 1 (RC-LH1) core complex of Rhodobacter sphaeroides converts absorbed light energy to a charge separation, and then it reduces a quinone electron and proton acceptor to a quinol. The angle between the two monomers imposes a bent configuration on the dimer complex, which exerts a major influence on the curvature of the membrane vesicles, known as chromatophores, where the light-driven photosynthetic reactions take place. To investigate the dimerisation interface between two RC-LH1 monomers, we determined the cryogenic electron microscopy structure of the dimeric complex at 2.9 Å resolution. The structure shows that each monomer consists of a central RC partly enclosed by a 14-subunit LH1 ring held in an open state by PufX and protein-Y polypeptides, thus enabling quinones to enter and leave the complex. Two monomers are brought together through N-terminal interactions between PufX polypeptides on the cytoplasmic side of the complex, augmented by two novel transmembrane polypeptides, designated protein-Z, that bind to the outer faces of the two central LH1 β polypeptides. The precise fit at the dimer interface, enabled by PufX and protein-Z, by C-terminal interactions between opposing LH1 αβ subunits, and by a series of interactions with a bound sulfoquinovosyl diacylglycerol lipid, bring together each monomer creating an S-shaped array of 28 bacteriochlorophylls. The seamless join between the two sets of LH1 bacteriochlorophylls provides a path for excitation energy absorbed by one half of the complex to migrate across the dimer interface to the other half.

Amphiphilic nanoparticles generate curvature in lipid membranes and shape liposome-liposome interfaces

We show by molecular dynamics that amphiphilic Au nanoparticles (NP) with a diameter of 4 nm generate curvature in phosphatidylcholine lipid membranes. NPs generate negative curvature when they adsorb on the membrane surface but, as they get spontaneously and progressively embedded into the membrane core, the curvature turns positive. As membrane embedding is kinetically slow, both configurations can be observed by Cryo-EM. NP-induced curvature explains the peculiar structure of liposome-liposome interfaces in presence of NPs.

On the vibrational free energy of hydrated proteins

When the hydration shell of a protein is filled with at least 0.6 gram of water per gram of protein, a significant anti-correlation between the vibrational free energy and the potential energy of energy-minimized conformers is observed. This means that low potential energy, well-hydrated, protein conformers tend to be more rigid than high-energy ones. On the other hand, in the case of CASP target 624, when its hydration shell is filled, a significant energy gap is observed between the crystal structure and the best conformers proposed during the prediction experiment, strongly suggesting that including explicit water molecules may help identifying unlikely conformers among good-looking ones.

Hybrid methods for combined experimental and computational determination of protein structure

Knowledge of protein structure is paramount to the understanding of biological function, developing new therapeutics, and making detailed mechanistic hypotheses. Therefore, methods to accurately elucidate three-dimensional structures of proteins are in high demand. While there are a few experimental techniques that can routinely provide high-resolution structures, such as x-ray crystallography, nuclear magnetic resonance (NMR), and cryo-EM, which have been developed to determine the structures of proteins, these techniques each have shortcomings and thus cannot be used in all cases. However, additionally, a large number of experimental techniques that provide some structural information, but not enough to assign atomic positions with high certainty have been developed. These methods offer sparse experimental data, which can also be noisy and inaccurate in some instances. In cases where it is not possible to determine the structure of a protein experimentally, computational structure prediction methods can be used as an alternative. Although computational methods can be performed without any experimental data in a large number of studies, inclusion of sparse experimental data into these prediction methods has yielded significant improvement. In this Perspective, we cover many of the successes of integrative modeling, computational modeling with experimental data, specifically for protein folding, protein-protein docking, and molecular dynamics simulations. We describe methods that incorporate sparse data from cryo-EM, NMR, mass spectrometry, electron paramagnetic resonance, small-angle x-ray scattering, Förster resonance energy transfer, and genetic sequence covariation. Finally, we highlight some of the major challenges in the field as well as possible future directions.

Inducible intracellular membranes: molecular aspects and emerging applications

Membrane remodeling and phospholipid biosynthesis are normally tightly regulated to maintain the shape and function of cells. Indeed, different physiological mechanisms ensure a precise coordination between de novo phospholipid biosynthesis and modulation of membrane morphology. Interestingly, the overproduction of certain membrane proteins hijack these regulation networks, leading to the formation of impressive intracellular membrane structures in both prokaryotic and eukaryotic cells. The proteins triggering an abnormal accumulation of membrane structures inside the cells (or membrane proliferation) share two major common features: (1) they promote the formation of highly curved membrane domains and (2) they lead to an enrichment in anionic, cone-shaped phospholipids (cardiolipin or phosphatidic acid) in the newly formed membranes. Taking into account the available examples of membrane proliferation upon protein overproduction, together with the latest biochemical, biophysical and structural data, we explore the relationship between protein synthesis and membrane biogenesis. We propose a mechanism for the formation of these non-physiological intracellular membranes that shares similarities with natural inner membrane structures found in α-proteobacteria, mitochondria and some viruses-infected cells, pointing towards a conserved feature through evolution. We hope that the information discussed in this review will give a better grasp of the biophysical mechanisms behind physiological and induced intracellular membrane proliferation, and inspire new applications, either for academia (high-yield membrane protein production and nanovesicle production) or industry (biofuel production and vaccine preparation).

How a lipid bilayer membrane responds to an oscillating nanoparticle: Promoted membrane undulation and directional wave propagation

Mechanical forces acting on a plasma membrane are of essential importance to cellular functioning via inducing delicate change of the membrane shape with the underlying mechanism yet to be elucidated. Here, we introduce an oscillating nanoparticle (NP) interaction with a lipid bilayer membrane, using the coarse-grained simulation to investigate the dynamic membrane response to constrained mechanical stimulation, which is ubiquitous in biology. Our results demonstrate that, the membrane responds to an oscillating NP by generating nanoscale undulation waves, which immediately propagate through the membrane. In dynamics, propagation of the generated membrane undulation waves always starts from flattening of the region where the NP locates, thus producing a lateral force to propel the waves away from the point of stimulation. The speed of membrane undulation wave propagation is proportional to that of NP oscillation and accelerated by increasing the integral membrane surface tension, suggesting that both the membrane bending and stretching contribute to the energy driving the unique response of membrane undulation wave propagation.

Atoms to Phenotypes: Molecular Design Principles of Cellular Energy Metabolism

We report a 100-million atom-scale model of an entire cell organelle, a photosynthetic chromatophore vesicle from a purple bacterium, that reveals the cascade of energy conversion steps culminating in the generation of ATP from sunlight. Molecular dynamics simulations of this vesicle elucidate how the integral membrane complexes influence local curvature to tune photoexcitation of pigments. Brownian dynamics of small molecules within the chromatophore probe the mechanisms of directional charge transport under various pH and salinity conditions. Reproducing phenotypic properties from atomistic details, a kinetic model evinces that low-light adaptations of the bacterium emerge as a spontaneous outcome of optimizing the balance between the chromatophore's structural integrity and robust energy conversion. Parallels are drawn with the more universal mitochondrial bioenergetic machinery, from whence molecular-scale insights into the mechanism of cellular aging are inferred. Together, our integrative method and spectroscopic experiments pave the way to first-principles modeling of whole living cells.

Membrane destruction and phospholipid extraction by using two-dimensional MoS2 nanosheets

The interaction of two-dimensional (2D) nanomaterials and bacterial membranes has attracted tremendous attention in antibacterial applications. Various peculiarities of 2D nanomaterials may lead to multiple mechanisms of their interactions with membranes. Here, we investigated the interaction between molybdenum disulfide (MoS2) nanosheets and the bacterial membrane by using both theoretical and experimental approaches. Molecular dynamics simulation presented that MoS2 nanosheets can disrupt the structure of the lipid membrane by making dents on its surface and extracting phospholipid molecules to reduce the integrity of the membrane. This is attributed to the combination of the dispersion interaction of lipid tails with S atoms and the electrostatic interactions of lipid head groups with the Mo and S atoms in the lateral edges of the MoS2 nanosheet. Scanning electron microscopy and transmission electron microscopy confirmed the dents and the destruction of the cell membrane, which would lead to the loss of cytoplasm and the death of bacteria. It should be noted that the phenomenon where MoS2 induces a dent is different from the direct insertion of graphene-based nanomaterials, which might be due to the thicker and stiffer structure of MoS2. Therefore, we believe that the molecular interactions of 2D nanomaterials with bacterial membranes should be highly correlated with their structural characteristics. This newly discovered mechanism of MoS2 nanomaterials to disrupt the cell membrane may promote the application of transition metal dichalcogenide (TMD) nanomaterials in designing remarkable antibacterial materials in the near future.

Biophysics of membrane curvature remodeling at molecular and mesoscopic lengthscales

At the micron scale, where cell organelles display an amazing complexity in their shape and organization, the physical properties of a biological membrane can be better-understood using continuum models subject to thermal (stochastic) undulations. Yet, the chief orchestrators of these complex and intriguing shapes are a specialized class of membrane associating often peripheral proteins called curvature remodeling proteins (CRPs) that operate at the molecular level through specific protein-lipid interactions. We review multiscale methodologies to model these systems at the molecular as well as at the mesoscopic and cellular scales, and also present a free energy perspective of membrane remodeling through the organization and assembly of CRPs. We discuss the morphological space of nearly planar to highly curved membranes, methods to include thermal fluctuations, and review studies that model such proteins as curvature fields to describe the emergent curved morphologies. We also discuss several mesoscale models applied to a variety of cellular processes, where the phenomenological parameters (such as curvature field strength) are often mapped to models of real systems based on molecular simulations. Much insight can be gained from the calculation of free energies of membranes states with protein fields, which enable accurate mapping of the state and parameter values at which the membrane undergoes morphological transformations such as vesiculation or tubulation. By tuning the strength, anisotropy, and spatial organization of the curvature-field, one can generate a rich array of membrane morphologies that are highly relevant to shapes of several cellular organelles. We review applications of these models to budding of vesicles commonly seen in cellular signaling and trafficking processes such as clathrin mediated endocytosis, sorting by the ESCRT protein complexes, and cellular exocytosis regulated by the exocyst complex. We discuss future prospects where such models can be combined with other models for cytoskeletal assembly, and discuss their role in understanding the effects of cell membrane tension and the mechanics of the extracellular microenvironment on cellular processes.

Membrane vesiculation induced by proteins of the dengue virus envelope studied by molecular dynamics simulations

Biological membranes are continuously remodeled in the cell by specific membrane-shaping machineries to form, for example, tubes and vesicles. We examine fundamental mechanisms involved in the vesiculation processes induced by a cluster of envelope (E) and membrane (M) proteins of the dengue virus (DENV) using molecular dynamics simulations and a coarse-grained model. We show that an arrangement of three E-M heterotetramers (EM3) works as a bending unit and an ordered cluster of five such units generates a closed vesicle, reminiscent of the virus budding process. In silico mutagenesis of two charged residues of the anchor helices of the envelope proteins of DENV shows that Arg-471 and Arg-60 are fundamental to produce bending stress on the membrane. The fine-tuning between the size of the EM3 unit and its specific bending action suggests this protein unit is an important factor in determining the viral particle size.

Iterative Molecular Dynamics-Rosetta Membrane Protein Structure Refinement Guided by Cryo-EM Densities

Knowing atomistic details of proteins is essential not only for the understanding of protein function but also for the development of drugs. Experimental methods such as X-ray crystallography, NMR, and cryo-electron microscopy (cryo-EM) are the preferred forms of protein structure determination and have achieved great success over the most recent decades. Computational methods may be an alternative when experimental techniques fail. However, computational methods are severely limited when it comes to predicting larger macromolecule structures with little sequence similarity to known structures. The incorporation of experimental restraints in computational methods is becoming increasingly important to more reliably predict protein structure. One such experimental input used in structure prediction and refinement is cryo-EM densities. Recent advances in cryo-EM have arguably revolutionized the field of structural biology. Our previously developed cryo-EM-guided Rosetta-MD protocol has shown great promise in the refinement of soluble protein structures. In this study, we extended cryo-EM density-guided iterative Rosetta-MD to membrane proteins. We also improved the methodology in general by picking models based on a combination of their score and fit-to-density during the Rosetta model selection. By doing so, we have been able to pick models superior to those with the previous selection based on Rosetta score only and we have been able to further improve our previously refined models of soluble proteins. The method was tested with five membrane spanning protein structures. By applying density-guided Rosetta-MD iteratively we were able to refine the predicted structures of these membrane proteins to atomic resolutions. We also showed that the resolution of the density maps determines the improvement and quality of the refined models. By incorporating high-resolution density maps (∼4 Å), we were able to more significantly improve the quality of the models than when medium-resolution maps (6.9 Å) were used. Beginning from an average starting structure root mean square deviation (RMSD) to native of 4.66 Å, our protocol was able to refine the structures to bring the average refined structure RMSD to 1.66 Å when 4 Å density maps were used. The protocol also successfully refined the HIV-1 CTD guided by an experimental 5 Å density map.

Engineering of a calcium-ion binding site into the RC-LH1-PufX complex of Rhodobacter sphaeroides to enable ion-dependent spectral red-shifting

The reaction centre-light harvesting 1 (RC-LH1) complex of Thermochromatium (Tch.) tepidum has a unique calcium-ion binding site that enhances thermal stability and red-shifts the absorption of LH1 from 880nm to 915nm in the presence of calcium-ions. The LH1 antenna of mesophilic species of phototrophic bacteria such as Rhodobacter (Rba.) sphaeroides does not possess such properties. We have engineered calcium-ion binding into the LH1 antenna of Rba. sphaeroides by progressively modifying the native LH1 polypeptides with sequences from Tch. tepidum. We show that acquisition of the C-terminal domains from LH1 α and β of Tch. tepidum is sufficient to activate calcium-ion binding and the extent of red-shifting increases with the proportion of Tch. tepidum sequence incorporated. However, full exchange of the LH1 polypeptides with those of Tch. tepidum results in misassembled core complexes. Isolated α and β polypeptides from our most successful mutant were reconstituted in vitro with BChl a to form an LH1-type complex, which was stabilised 3-fold by calcium-ions. Additionally, carotenoid specificity was changed from spheroidene found in Rba. sphaeroides to spirilloxanthin found in Tch. tepidum, with the latter enhancing in vitro formation of LH1. These data show that the C-terminal LH1 α/β domains of Tch. tepidum behave autonomously, and are able to transmit calcium-ion induced conformational changes to BChls bound to the rest of a foreign antenna complex. Thus, elements of foreign antenna complexes, such as calcium-ion binding and blue/red switching of absorption, can be ported into Rhodobacter sphaeroides using careful design processes.

The C-terminus of PufX plays a key role in dimerisation and assembly of the reaction center light-harvesting 1 complex from Rhodobacter sphaeroides

In bacterial photosynthesis reaction center-light-harvesting 1 (RC-LH1) complexes trap absorbed solar energy by generating a charge separated state. Subsequent electron and proton transfers form a quinol, destined to diffuse to the cytochrome bc₁ complex. In bacteria such as Rhodobacter (Rba.) sphaeroides and Rba. capsulatus the PufX polypeptide creates a channel for quinone/quinol traffic across the LH1 complex that surrounds the RC, and it is therefore essential for photosynthetic growth. PufX also plays a key role in dimerization of the RC-LH1-PufX core complex, and the structure of the Rba. sphaeroides complex shows that the PufX C-terminus, particularly the region from X49-X53, likely mediates association of core monomers. To investigate this putative interaction we analysed mutations PufX R49L, PufX R53L, PufX R49/53L and PufX G52L by measuring photosynthetic growth, fractionation of detergent-solubilised membranes, formation of 2-D crystals and electron microscopy. We show that these mutations do not affect assembly of PufX within the core or photosynthetic growth but they do prevent dimerization, consistent with predictions from the RC-LH1-PufX structure. We obtained low resolution structures of monomeric core complexes with and without PufX, using electron microscopy of negatively stained single particles and 3D reconstruction; the monomeric complex with PufX corresponds to one half of the dimer structure whereas LH1 completely encloses the RC if the gene encoding PufX is deleted. On the basis of the insights gained from these mutagenesis and structural analyses we propose a sequence for assembly of the dimeric RC-LH1-PufX complex.

CHARMM-GUI MDFF/xMDFF Utilizer for Molecular Dynamics Flexible Fitting Simulations in Various Environments

X-ray crystallography and cryo-electron microscopy are two popular methods for the structure determination of biological molecules. Atomic structures are derived through the fitting and refinement of an initial model into electron density maps constructed by both experiments. Two computational approaches, MDFF and xMDFF, have been developed to facilitate this process by integrating the experimental data with molecular dynamics simulation. However, the setup of an MDFF/xMDFF simulation requires knowledge of both experimental and computational methods, which is not straightforward for nonexpert users. In addition, sometimes it is desirable to include realistic environments, such as explicit solvent and lipid bilayers during the simulation, which poses another challenge even for expert users. To alleviate these difficulties, we have developed MDFF/xMDFF Utilizer in CHARMM-GUI that helps users to set up an MDFF/xMDFF simulation. The capability of MDFF/xMDFF Utilizer is greatly enhanced by integration with other CHARMM-GUI modules, including protein structure manipulation, a diverse set of lipid types, and all-atom CHARMM and coarse-grained PACE force fields. With this integration, various simulation environments are available for MDFF Utilizer (vacuum, implicit/explicit solvent, and bilayers) and xMDFF Utilizer (vacuum and solution). In this work, three examples are shown to demonstrate the usage of MDFF/xMDFF Utilizer.

Overall energy conversion efficiency of a photosynthetic vesicle

The chromatophore of purple bacteria is an intracellular spherical vesicle that exists in numerous copies in the cell and that efficiently converts sunlight into ATP synthesis, operating typically under low light conditions. Building on an atomic-level structural model of a low-light-adapted chromatophore vesicle from Rhodobacter sphaeroides, we investigate the cooperation between more than a hundred protein complexes in the vesicle. The steady-state ATP production rate as a function of incident light intensity is determined after identifying quinol turnover at the cytochrome bc1 complex (cytb⁢c1) as rate limiting and assuming that the quinone/quinol pool of about 900 molecules acts in a quasi-stationary state. For an illumination condition equivalent to 1% of full sunlight, the vesicle exhibits an ATP production rate of 82 ATP molecules/s. The energy conversion efficiency of ATP synthesis at illuminations corresponding to 1%–5% of full sunlight is calculated to be 0.12–0.04, respectively. The vesicle stoichiometry, evolutionarily adapted to the low light intensities in the habitat of purple bacteria, is suboptimal for steady-state ATP turnover for the benefit of protection against over-illumination.

DOI:http://dx.doi.org/10.7554/eLife.09541.001

Photosynthesis, or the conversion of light energy into chemical energy, is a process that powers almost all life on Earth. Plants and certain bacteria share similar processes to perform photosynthesis, though the purple bacterium Rhodobacter sphaeroides uses a photosynthetic system that is much less complex than that in plants. Light harvesting inside the bacterium takes place in up to hundreds of compartments called chromatophores. Each chromatophore in turn contains hundreds of cooperating proteins that together absorb the energy of sunlight and convert and store it in molecules of ATP, the universal energy currency of all cells.

The chromatophore of primitive purple bacteria provides a model for more complex photosynthetic systems in plants. Though researchers had characterized its individual components over the years, less was known about the overall architecture of the chromatophore and how its many components work together to harvest light energy efficiently and robustly. This knowledge would provide insight into the evolutionary pressures that shaped the chromatophore and its ability to work efficiently at different light intensities.

Sener et al. now present a highly detailed structural model of the chromatophore of purple bacteria based on the findings of earlier studies. The model features the position of every atom of the constituent proteins and is used to examine how energy is transferred and converted. Sener et al. describe the sequence of energy conversion steps and calculate the overall energy conversion efficiency, namely how much of the light energy arriving at the microorganism is stored as ATP.

These calculations show that the chromatophore is optimized to produce chemical energy at low light levels typical of purple bacterial habitats, and dissipate excess energy to avoid being damaged under brighter light. The chromatophore’s architecture also displays robustness against perturbations of its components. In the future, the approach used by Sener et al. to describe light harvesting in this bacterial compartment can be applied to more complex systems, such as those in plants.

DOI:http://dx.doi.org/10.7554/eLife.09541.002

Curvature-undulation coupling as a basis for curvature sensing and generation in bilayer membranes

We present coarse-grained molecular dynamics simulations of the epsin N-terminal homology domain interacting with a lipid bilayer and demonstrate a rigorous theoretical formalism and analysis method for computing the induced curvature field in varying concentrations of the protein in the dilute limit. Our theory is based on the description of the height-height undulation spectrum in the presence of a curvature field. We formulated an objective function to compare the acquired undulation spectrum from the simulations to that of the theory. We recover the curvature field parameters by minimizing the objective function even in the limit where the protein-induced membrane curvature is of the same order as the amplitude due to thermal undulations. The coupling between curvature and undulations leads to significant predictions: (i) Under dilute conditions, the proteins can sense a site of spontaneous curvature at distances much larger than their size; (ii) as the density of proteins increases the coupling focuses and stabilizes the curvature field to the site of the proteins; and (iii) the mapping of the protein localization and the induction of a stable curvature is a cooperative process that can be described through a Hill function.

Molecular dynamics-based refinement and validation for sub-5 Å cryo-electron microscopy maps

Two structure determination methods, based on the molecular dynamics flexible fitting (MDFF) paradigm, are presented that resolve sub-5 Å cryo-electron microscopy (EM) maps with either single structures or ensembles of such structures. The methods, denoted cascade MDFF and resolution exchange MDFF, sequentially re-refine a search model against a series of maps of progressively higher resolutions, which ends with the original experimental resolution. Application of sequential re-refinement enables MDFF to achieve a radius of convergence of ~25 Å demonstrated with the accurate modeling of β-galactosidase and TRPV1 proteins at 3.2 Å and 3.4 Å resolution, respectively. The MDFF refinements uniquely offer map-model validation and B-factor determination criteria based on the inherent dynamics of the macromolecules studied, captured by means of local root mean square fluctuations. The MDFF tools described are available to researchers through an easy-to-use and cost-effective cloud computing resource on Amazon Web Services.

DOI:http://dx.doi.org/10.7554/eLife.16105.001

To understand the roles that proteins and other large molecules play inside cells, it is important to determine their structures. One of the techniques that researchers can use to do this is called cryo-electron microscopy (cryo-EM), which rapidly freezes molecules to fix them in position before imaging them in fine detail.

The cryo-EM images are like maps that show the approximate position of atoms. These images must then be processed in order to build a three-dimensional model of the protein that shows how its atoms are arranged relative to each other. One computational approach called Molecular Dynamics Flexible Fitting (MDFF) works by flexibly fitting possible atomic structures into cryo-EM maps. Although this approach works well with relatively undetailed (or ‘low resolution’) cryo-EM images, it struggles to handle the high-resolution cryo-EM maps now being generated.

Singharoy, Teo, McGreevy et al. have now developed two MDFF methods – called cascade MDFF and resolution exchange MDFF – that help to resolve atomic models of biological molecules from cryo-EM images. Each method can refine poorly guessed models into ones that are consistent with the high-resolution experimental images. The refinement is achieved by interpreting a range of images that starts with a ‘fuzzy’ image. The contrast of the image is then progressively improved until an image is produced that has a resolution that is good enough to almost distinguish individual atoms.

The method works because each cryo-EM image shows not just one, but a collection of atomic structures that the molecule can take on, with the fuzzier parts of the image representing the more flexible parts of the molecule. By taking into account this flexibility, the large-scale features of the protein structure can be determined first from the fuzzier images, and increasing the contrast of the images allows smaller-scale refinements to be made to the structure.

The MDFF tools have been designed to be easy to use and are available to researchers at low cost through cloud computing platforms. They can now be used to unravel the structure of many different proteins and protein complexes including those involved in photosynthesis, respiration and protein synthesis.

DOI:http://dx.doi.org/10.7554/eLife.16105.002

Genetic and Ultrastructural Analysis Reveals the Key Players and Initial Steps of Bacterial Magnetosome Membrane Biogenesis

Magnetosomes of magnetotactic bacteria contain well-ordered nanocrystals for magnetic navigation and have recently emerged as the most sophisticated model system to study the formation of membrane bounded organelles in prokaryotes. Magnetosome biosynthesis is thought to begin with the formation of a dedicated compartment, the magnetosome membrane (MM), in which the biosynthesis of a magnetic mineral is strictly controlled. While the biomineralization of magnetosomes and their subsequent assembly into linear chains recently have become increasingly well studied, the molecular mechanisms and early stages involved in MM formation remained poorly understood. In the Alphaproteobacterium Magnetospirillum gryphiswaldense, approximately 30 genes were found to control magnetosome biosynthesis. By cryo-electron tomography of several key mutant strains we identified the gene complement controlling MM formation in this model organism. Whereas the putative magnetosomal iron transporter MamB was most crucial for the process and caused the most severe MM phenotype upon elimination, MamM, MamQ and MamL were also required for the formation of wild-type-like MMs. A subset of seven genes (mamLQBIEMO) combined within a synthetic operon was sufficient to restore the formation of intracellular membranes in the absence of other genes from the key mamAB operon. Tracking of de novo magnetosome membrane formation by genetic induction revealed that magnetosomes originate from unspecific cytoplasmic membrane locations before alignment into coherent chains. Our results indicate that no single factor alone is essential for MM formation, which instead is orchestrated by the cumulative action of several magnetosome proteins.

One of the most intriguing examples for membrane-bounded prokaryotic organelles are magnetosomes which consist of well-ordered chains of perfectly shaped magnetic nanocrystals that in many aquatic bacteria serve as geomagnetic field sensors to direct their swimming towards microoxic zones at the bottom of natural waters. In the model bacterium Magnetospirillum gryphiswaldense and related magnetotactic microorganisms, magnetosomes are formed by a complex pathway that is orchestrated by more than 30 genes. However, the initial and most crucial step of magnetosome biosynthesis, formation and differentiation of a dedicated intracellular membrane compartment for controlled biomineralization of magnetite crystals, remained only poorly understood. By ultrastructural analysis of several mutants and genetic induction of de novo magnetosome synthesis, we identified the key determinants and early steps of magnetosome membrane biogenesis. Our results suggest that formation of intracellular membranes in bacteria is mediated by a cumulative action of several factors, but apparently is differently controlled than intracellular membrane remodeling in eukaryotic cells.

Advances in the molecular dynamics flexible fitting method for cryo-EM modeling

Molecular Dynamics Flexible Fitting (MDFF) is an established technique for fitting all-atom structures of molecules into corresponding cryo-electron microscopy (cryo-EM) densities. The practical application of MDFF is simple but requires a user to be aware of and take measures against a variety of possible challenges presented by each individual case. Some of these challenges arise from the complexity of a molecular structure or the limited quality of available structural models and densities to be interpreted, while others stem from the intricacies of MDFF itself. The current article serves as an overview of the strategies that have been developed since MDFF's inception to overcome common challenges and successfully perform MDFF simulations.

Simulations of cellulose translocation in the bacterial cellulose synthase suggest a regulatory mechanism for the dimeric structure of cellulose

In addition to suggesting a mechanism for regulating cellulose structure, molecular simulations indicate translocation is not rate-limiting for cellulose biosynthesis.

The processive cycle of the bacterial cellulose synthase (Bcs) includes the addition of a single glucose moiety to the end of a growing cellulose chain followed by the translocation of the nascent chain across the plasma membrane. The mechanism of this translocation and its precise location within the processive cycle are not well understood. In particular, the molecular details of how a polymer (cellulose) whose basic structural unit is a dimer (cellobiose) can be constructed by adding one monomer (glucose) at a time are yet to be elucidated. Here, we have utilized molecular dynamics simulations and free energy calculations to the shed light on these questions. We find that translocation forward by one glucose unit is quite favorable energetically, giving a free energy stabilization of greater than 10 kcal mol^–1. In addition, there is only a small barrier to translocation, implying that translocation is not rate limiting within the Bcs processive cycle (given experimental rates for cellulose synthesis in vitro). Perhaps most significantly, our results also indicate that steric constraints at the transmembrane tunnel entrance regulate the dimeric structure of cellulose. Namely, when a glucose molecule is added to the cellulose chain in the same orientation as the acceptor glucose, the terminal glucose freely rotates upon forward motion, thus suggesting a regulatory mechanism for the dimeric structure of cellulose. We characterize both the conserved and non-conserved enzyme–polysaccharide interactions that drive translocation, and find that 20 of the 25 residues that strongly interact with the translocating cellulose chain in the simulations are well conserved, mostly with polar or aromatic side chains. Our results also allow for a dynamical analysis of the role of the so-called ‘finger helix’ in cellulose translocation that has been observed structurally. Taken together, these findings aid in the elucidation of the translocation steps of the Bcs processive cycle and may be widely relevant to polysaccharide synthesizing or degrading enzymes that couple catalysis with chain translocation.

Development and dynamics of the photosynthetic apparatus in purple phototrophic bacteria

The purple bacterium Rhodobacter sphaeroides provides a useful model system for studies of the assembly and dynamics of bacterial photosynthetic membranes. For the nascent developing membrane, proteomic analyses showed an ~2-fold enrichment in general membrane assembly factors, compared to chromatophores. When the protonophore carbonyl-cyanide m-chlorophenyl-hydrazone (CCCP) was added to an ICM inducing culture, an ~2-fold elevation in spectral counts vs. the control was seen for the SecA translocation ATPase, the preprotein translocase SecY, SecD and SecF insertion components, and chaperonins DnaJ and DnaK, which act early in the assembly process. It is suggested that these factors accumulated with their nascent polypeptides, as putative assembly intermediates in a functionally arrested state. Since in Synechocystis PCC 6803, a link has been established between Chl delivery involving the high-light HilD protein and the SecY/YidC-requiring cotranslational insertion of nascent polypeptides, such a connection between BChl biosynthesis and insertion and folding of nascent Rba. sphaeroides BChl binding proteins is likely to also occur. AFM imaging studies of the formation of the reaction center (RC)-light harvesting 1 (LH1) complex suggested a cooperative assembly mechanism in which, following the association between the RC template and the initial LH1 unit, addition of successive LH1 units to the RC drives the assembly process to completion. Alterations in membrane dynamics as the developing membrane becomes filled with LH2-rings were assessed by fluorescence induction/relaxation kinetics, which showed a slowing in RC electron transfer rate thought to mainly reflect alterations in donor side electron transfer. This was attributed to an increased distance for electron flow in cytochrome c2 between the RC and cytochrome bc1 complexes, as suggested in the current structural models. This article is part of a Special Issue entitled Organization and dynamics of bioenergetic systems in bacteria, edited by Prof Conrad Mullineaux.

Structure based computational assessment of channel properties of assembled ORF-8a from SARS-CoV

ORF 8a is a short 39 amino acid bitopic membrane protein encoded by severe acute respiratory syndrome causing corona virus (SARS‐CoV). It has been identified to increase permeability of the lipid membrane for cations. Permeability is suggested to occur due to the assembly of helical bundles. Computational models of a pentameric assembly of 8a peptides are generated using the first 22 amino acids, which include the transmembrane domain. Low energy structures reveal a hydrophilic pore mantled by residues Thr‐8, and −18, Ser‐11, Cys‐13, and Arg‐22. Potential of mean force (PMF) profiles for mono (Na⁺, K⁺, Cl⁻) and divalent (Ca²⁺) ions along the pore are calculated. The data support experimental findings of a weak cation selectivity of the channel. Calculations on 8a are compared to data derived for a pentameric bundle consisting of the M2 helices of the bacterial pentameric ligand gated ion channel GLIC (3EHZ). PMF curves of both, bundles 8a and M2, show sigmoidal shaped profiles. In comparison to the data for the M2‐GLIC model, data of the 8a bundle show lower amplitude of the PMF values between maximum and minimum and less discrimination amongst ions. Proteins 2015; 83:300–308. © 2014 Wiley Periodicals, Inc.

Flexibility in the insulin receptor ectodomain enables docking of insulin in crystallographic conformation observed in a hormone-bound microreceptor

Insulin binding to the insulin receptor (IR) is the first key step in initiating downstream signaling cascades for glucose homeostasis in higher organisms. The molecular details of insulin recognition by IR are not yet completely understood, but a picture of hormone/receptor interactions at one of the epitopes (Site 1) is beginning to emerge from recent structural evidence. However, insulin-bound structures of truncated IR suggest that crystallographic conformation of insulin cannot be accommodated in the full IR ectodomain due to steric overlap of insulin with the first two type III fibronectin domains (F1 and F2), which are contributed to the insulin binding-pocket by the second subunit in the IR homodimer. A conformational change in the F1-F2 pair has thus been suggested. In this work, we present an all-atom structural model of complex of insulin and the IR ectodomain, where no structural overlap of insulin with the receptor domains (F1 and F2) is observed. This structural model was arrived at by flexibly fitting parts of our earlier insulin/IR all-atom model into the simulated density maps of crystallized constructs combined with conformational sampling from apo-IR solution conformations. Importantly, our experimentally-consistent model helps rationalize yet unresolved Site.

xMDFF: molecular dynamics flexible fitting of low-resolution X-ray structures

X-ray crystallography remains the most dominant method for solving atomic structures. However, for relatively large systems, the availability of only medium-to-low-resolution diffraction data often limits the determination of all-atom details. A new molecular dynamics flexible fitting (MDFF)-based approach, xMDFF, for determining structures from such low-resolution crystallographic data is reported. xMDFF employs a real-space refinement scheme that flexibly fits atomic models into an iteratively updating electron-density map. It addresses significant large-scale deformations of the initial model to fit the low-resolution density, as tested with synthetic low-resolution maps of D-ribose-binding protein. xMDFF has been successfully applied to re-refine six low-resolution protein structures of varying sizes that had already been submitted to the Protein Data Bank. Finally, via systematic refinement of a series of data from 3.6 to 7 Å resolution, xMDFF refinements together with electrophysiology experiments were used to validate the first all-atom structure of the voltage-sensing protein Ci-VSP.

Organization in photosynthetic membranes of purple bacteria in vivo: the role of carotenoids

Photosynthesis in purple bacteria is performed by pigment-protein complexes that are closely packed within specialized intracytoplasmic membranes. Here we report on the influence of carotenoid composition on the organization of RC-LH1 pigment-protein complexes in intact membranes and cells of Rhodobacter sphaeroides. Mostly dimeric RC-LH1 complexes could be isolated from strains expressing native brown carotenoids when grown under illuminated/anaerobic conditions, or from strains expressing green carotenoids when grown under either illuminated/anaerobic or dark/semiaerobic conditions. However, mostly monomeric RC-LH1 complexes were isolated from strains expressing the native photoprotective red carotenoid spheroidenone, which is synthesized during phototrophic growth in the presence of oxygen. Despite this marked difference, linear dichroism (LD) and light-minus-dark LD spectra of oriented intact intracytoplasmic membranes indicated that RC-LH1 complexes are always assembled in ordered arrays, irrespective of variations in the relative amounts of isolated dimeric and monomeric RC-LH1 complexes. We propose that part of the photoprotective response to the presence of oxygen mediated by synthesis of spheroidenone may be a switch of the structure of the RC-LH1 complex from dimers to monomers, but that these monomers are still organized into the photosynthetic membrane in ordered arrays. When levels of the dimeric RC-LH1 complex were very high, and in the absence of LH2, LD and ∆LD spectra from intact cells indicated an ordered arrangement of RC-LH1 complexes. Such a degree of ordering implies the presence of highly elongated, tubular membranes with dimensions requiring orientation along the length of the cell and in a proportion larger than previously observed.

GPU-accelerated analysis and visualization of large structures solved by molecular dynamics flexible fitting

Hybrid structure fitting methods combine data from cryo-electron microscopy and X-ray crystallography with molecular dynamics simulations for the determination of all-atom structures of large biomolecular complexes. Evaluating the quality-of-fit obtained from hybrid fitting is computationally demanding, particularly in the context of a multiplicity of structural conformations that must be evaluated. Existing tools for quality-of-fit analysis and visualization have previously targeted small structures and are too slow to be used interactively for large biomolecular complexes of particular interest today such as viruses or for long molecular dynamics trajectories as they arise in protein folding. We present new data-parallel and GPU-accelerated algorithms for rapid interactive computation of quality-of-fit metrics linking all-atom structures and molecular dynamics trajectories to experimentally-determined density maps obtained from cryo-electron microscopy or X-ray crystallography. We evaluate the performance and accuracy of the new quality-of-fit analysis algorithms vis-à-vis existing tools, examine algorithm performance on GPU-accelerated desktop workstations and supercomputers, and describe new visualization techniques for results of hybrid structure fitting methods.

The architecture of Rhodobacter sphaeroides chromatophores

The chromatophores of Rhodobacter (Rb.) sphaeroides represent a minimal bio-energetic system, which efficiently converts light energy into usable chemical energy. Despite extensive studies, several issues pertaining to the morphology and molecular architecture of this elemental energy conversion system remain controversial or unknown. To tackle these issues, we combined electron microscope tomography, immuno-electron microscopy and atomic force microscopy. We found that the intracellular Rb. sphaeroides chromatophores form a continuous reticulum rather than existing as discrete vesicles. We also found that the cytochrome bc1 complex localizes to fragile chromatophore regions, which most likely constitute the tubular structures that interconnect the vesicles in the reticulum. In contrast, the peripheral light-harvesting complex 2 (LH2) is preferentially hexagonally packed within the convex vesicular regions of the membrane network. Based on these observations, we propose that the bc1 complexes are in the inter-vesicular regions and surrounded by reaction center (RC) core complexes, which in turn are bounded by arrays of peripheral antenna complexes. This arrangement affords rapid cycling of electrons between the core and bc1 complexes while maintaining efficient excitation energy transfer from LH2 domains to the RCs.

Variation in supramolecular organisation of the photosynthetic membrane of Rhodobacter sphaeroides induced by alteration of PufX

In purple bacteria of the genus Rhodobacter (Rba.), an LH1 antenna complex surrounds the photochemical reaction centre (RC) with a PufX protein preventing the LH1 complex from completely encircling the RC. In membranes of Rba. sphaeroides, RC-LH1 complexes associate as dimers which in turn assemble into longer range ordered arrays. The present work uses linear dichroism (LD) and dark-minus-light difference LD (ΔLD) to probe the organisation of genetically altered RC-LH1 complexes in intact membranes. The data support previous proposals that Rba. capsulatus, and Rba. sphaeroides heterologously expressing the PufX protein from Rba. capsulatus, produce monomeric core complexes in membranes that lack long-range order. Similarly, Rba. sphaeroides with a point mutation in the Gly 51 residue of PufX, which is located on the membrane-periplasm interface, assembles mainly non-ordered RC-LH1 complexes that are most likely monomeric. All the Rba. sphaeroides membranes in their ΔLD spectra exhibited a spectral fingerprint of small degree of organisation implying the possibility of ordering influence of LH1, and leading to an important conclusion that PufX itself has no influence on ordering RC-LH1 complexes, as long-range order appears to be induced only through its role of configuring RC-LH1 complexes into dimers.

The role of shape complementarity in the protein-protein interactions

We use a dissipative particle dynamic simulation to investigate the effects of shape complementarity on the protein-protein interactions. By monitoring different kinds of protein shape-complementarity modes, we gave a clear mechanism to reveal the role of the shape complementarity in the protein-protein interactions, i.e., when the two proteins with shape complementarity approach each other, the conformation of lipid chains between two proteins would be restricted significantly. The lipid molecules tend to leave the gap formed by two proteins to maximize the configuration entropy, and therefore yield an effective entropy-induced protein-protein attraction, which enhances the protein aggregation. In short, this work provides an insight into understanding the importance of the shape complementarity in the protein-protein interactions especially for protein aggregation and antibody-antigen complexes. Definitely, the shape complementarity is the third key factor affecting protein aggregation and complex, besides the electrostatic-complementarity and hydrophobic complementarity.

Membrane undulation induced by NS4A of Dengue virus: a molecular dynamics simulation study

Nonstructural protein 4A (NS4A) of Dengue virus (DENV) is a membrane protein involved in rearrangements of the endoplasmic reticulum membrane that are required for formation of replication vesicles. NS4A is composed most likely of three membrane domains. The N- and C-terminal domains are supposed to traverse the lipid membrane whereas the central one is thought to reside on the membrane surface, thus forming a u-shaped protein. All three membrane domains are proposed to be helical by secondary structure prediction programs. After performing multi nanosecond molecular dynamics (MD) simulations at various temperatures (300, 310, and 315.15 K) with each of the individual domains, they are used in a docking approach to define putative association motifs of the transmembrane domains (TMDs). Two structures of the u-shaped protein are generated by separating two assembled TMDs linking them with the membrane-attached domain. Lipid undulation is monitored with the structures embedded in a fully hydrated lipid bilayer applying multiple 200 ns MD simulations at 310 K. An intact structure of the protein supports membrane undulation. The strong unwinding of the helices in the domain-linking section of one of the structures lowers its capability to induce membrane curvature. Unwinding of the link region is due to interactions of two tryptophan residues, Trp-96 and 104. These results provide first insights into the membrane-altering properties of DENV NS4A.

Structural implications of hydrogen-bond energetics in membrane proteins revealed by high-pressure spectroscopy

The light-harvesting 1 (LH1) integral membrane complex of Rhodobacter sphaeroides provides a convenient model system in which to examine the poorly understood role of hydrogen bonds (H-bonds) as stabilizing factors in membrane protein complexes. We used noncovalently bound arrays of bacteriochlorophyll chromophores within native and genetically modified variants of LH1 complexes to monitor local changes in the chromophore binding sites induced by externally applied hydrostatic pressure. Whereas membrane-bound complexes demonstrated very high resilience to pressures reaching 2.1 GPa, characteristic discontinuous shifts and broadenings of the absorption spectra were observed around 1 GPa for detergent-solubilized proteins, in similarity to those observed when specific (α or β) H-bonds between the chromophores and the surrounding protein were selectively removed by mutagenesis. These pressure effects, which were reversible upon decompression, allowed us to estimate the rupture energies of H-bonds to the chromophores in LH1 complexes. A quasi-independent, additive role of H-bonds in the α- and β-sublattices in reinforcing the wild-type LH1 complex was established. A comparison of a reaction-center-deficient LH1 complex with complexes containing reaction centers also demonstrated a stabilizing effect of the reaction center. This study thus provides important insights into the design principles of natural photosynthetic complexes.

Investigation of photosynthetic membrane structure using atomic force microscopy

Photosynthetic processes, including light capture, electron transfer, and energy conversion, are not only ensured by the activities of individual photosynthetic complexes but also substantially determined and regulated by the composition and assembly of the overall photosynthetic apparatus at the supramolecular level. In recent years, atomic force microscopy (AFM) has matured as a unique and powerful tool for directly assessing the supramolecular assembly of integral membrane protein complexes in their native membrane environment at submolecular resolution. This review highlights the major contributions and advances of AFM studies to our understanding of the structure of the bacterial photosynthetic machinery and its regulatory arrangement during chromatic adaptation. AFM topographs of other biological membrane systems and potential future applications of AFM are also discussed.

Enhanced sampling and overfitting analyses in structural refinement of nucleic acids into electron microscopy maps

Flexible fitting computational algorithms are often useful to interpret low-resolution maps of many macromolecular complexes generated by electron microscopy (EM) imaging. One such atomistic simulation technique is molecular dynamics flexible fitting (MDFF), which has been widely applied to generate structural models of large ribonucleoprotein assemblies such as the ribosome. We have previously shown that MDFF simulations of globular proteins are sensitive to the resolution of the target EM map and the strength of restraints used to preserve the secondary structure elements during fitting (Vashisth, H.; et al. Structure 2012, 20, 1453-1462). In this work, we aim to systematically examine the quality of structural models of various nucleic acids obtained via MDFF by varying the map resolution and the strength of structural restraints. We also demonstrate how an enhanced conformational sampling technique for proteins, temperature-accelerated molecular dynamics (TAMD), can be combined with MDFF for the structural refinement of nucleic acids in EM maps. Finally, we also demonstrate application of TAMD-assisted MDFF (TAMDFF) on a RNA/protein complex and suggest that TAMDFF is a viable strategy for enhanced conformational fitting in target maps of ribonucleoprotein complexes.

Consensus among multiple approaches as a reliability measure for flexible fitting into cryo-EM data

Cryo-electron microscopy (cryo-EM) can provide low-resolution density maps of large macromolecular assemblies. As the number of structures deposited in the Protein Data Bank by fitting a high-resolution structure into a low-resolution cryo-EM map is increasing, there is a need to revise the protocols and improve the measures for fitting. A recent study suggested using a combination of multiple automated flexible fitting approaches to improve the interpretation of cryo-EM data. The current work further explores the use of multiple approaches by validating this "consensus" fitting approach and deriving a local reliability measure. Here four different flexible fitting approaches are applied for fitting an initial structure into a simulated density map of known target structure from a dataset of proteins. It is found that the models produced from different approaches often have a consensus in conformation and are also near to the target structure, whereas cases not showing consensus are away from the target. A high correlation is also observed between the RMSF profiles calculated with respect to the average and the target structures, which indicates that the relation between consensus and accuracy can also be extended to a per-residue level. Therefore, the RMSF among the fitted models is proposed as a local reliability measure, which can be used to assess the reliability of the fit at specific regions. Hence, we encourage the community to use consensus flexible fitting with different methods to report on local reliability of the resulting models and improve the interpretation of cryo-EM data.

Using enhanced sampling and structural restraints to refine atomic structures into low-resolution electron microscopy maps

For a variety of problems in structural biology, low-resolution maps generated by electron microscopy imaging are often interpreted with the help of various flexible-fitting computational algorithms. In this work, we systematically analyze the quality of final models of various proteins obtained via molecular dynamics flexible fitting (MDFF) by varying the map-resolution, strength of structural restraints, and the steering forces. We find that MDFF can be extended to understand conformational changes in lower-resolution maps if larger structural restraints and lower steering forces are used to prevent overfitting. We further show that the capabilities of MDFF can be extended by combining it with an enhanced conformational sampling method, temperature-accelerated molecular dynamics (TAMD). Specifically, either TAMD can be used to generate better starting configurations for MDFF fitting or TAMD-assisted MDFF (TAMDFF) can be performed to accelerate conformational search in atomistic simulations.

Changes in lipid density induce membrane curvature

Highly curved bilayer lipid membranes make up the shell of many intra- and extracellular compartments, including organelles and vesicles. Using all-atom molecular dynamics simulations, we show that increasing the density of lipids in the bilayer membrane can induce the membrane to form a curved shape.

Cryo-electron microscopy modeling by the molecular dynamics flexible fitting method

The increasing power and popularity of cryo-electron microscopy (cryo-EM) in structural biology brought about the development of so-called hybrid methods, which permit the interpretation of cryo-EM density maps beyond their nominal resolution in terms of atomic models. The Cryo-EM Modeling Challenge 2010 is the first community effort to bring together developers of hybrid methods as well as cryo-EM experimentalists. Participating in the challenge, the molecular dynamics flexible fitting (MDFF) method was applied to a number of cryo-EM density maps. The results are described here with special emphasis on the use of symmetry-based restraints to improve the quality of atomic models derived from density maps of symmetric complexes; on a comparison of the stereochemical quality of atomic models resulting from different hybrid methods; and on application of MDFF to electron crystallography data.

Go hybrid: EM, crystallography, and beyond

A mechanistic understanding of the molecular transactions that govern cellular function requires knowledge of the dynamic organization of the macromolecular machines involved in these processes. Structural biologists employ a variety of biophysical methods to study large macromolecular complexes, but no single technique is likely to provide a complete description of the structure-function relationship of all the constituent components. Since structural studies generally only provide snapshots of these dynamic machines as they accomplish their molecular functions, combining data from many methodologies is crucial to our understanding of molecular function.

Adaptation of intracytoplasmic membranes to altered light intensity in Rhodobacter sphaeroides

The model photosynthetic bacterium Rhodobacter sphaeroides uses a network of bacteriochlorophyll (BChl)-protein complexes embedded in spherical intracytoplasmic membranes (ICM) to collect and utilise solar energy. We studied the effects of high- and low-light growth conditions, where BChl levels increased approximately four-fold from 1.6×10(6) to 6.5×10(6) molecules per cell. Most of this extra pigment is accommodated in the proliferating ICM system, which increases from approximately 274 to 1468 vesicles per cell. Thus, 16×10(6)nm(2) of specialised membrane surface area is made available for harvesting and utilising solar energy compared to 3×10(6)nm(2) under high-light conditions. Membrane mapping using atomic force microscopy revealed closely packed dimeric and monomeric reaction centre-light harvesting 1-PufX (RC-LH1-PufX) complexes in high-light ICM with room only for small clusters of LH2, whereas extensive LH2-only domains form during adaptation to low light, with the LH2/RC ratio increasing three-fold. The number of upper pigmented band (UPB) sites where membrane invagination is initiated hardly varied; 704 (5.8×10(5) BChls/cell) and 829 (4.9×10(5) BChls/cell) UPB sites per cell were estimated under high- and low-light conditions, respectively. Thus, the lower ICM content in high-light cells is a consequence of fewer ICM invaginations reaching maturity. Taking into account the relatively poor LH2-to-LH1 energy transfer in UPB membranes it is likely that high-light cells are relatively inefficient at energy trapping, but can grow well enough without the need to fully develop their photosynthetic membranes from the relatively inefficient UPB to highly efficient mature ICM.

Consensus among flexible fitting approaches improves the interpretation of cryo-EM data

Cryo-elecron microscopy (cryo-EM) can provide important structural information of large macromolecular assemblies in different conformational states. Recent years have seen an increase in structures deposited in the Protein Data Bank (PDB) by fitting a high-resolution structure into its low-resolution cryo-EM map. A commonly used protocol for accommodating the conformational changes between the X-ray structure and the cryo-EM map is rigid body fitting of individual domains. With the emergence of different flexible fitting approaches, there is a need to compare and revise these different protocols for the fitting. We have applied three diverse automated flexible fitting approaches on a protein dataset for which rigid domain fitting (RDF) models have been deposited in the PDB. In general, a consensus is observed in the conformations, which indicates a convergence from these theoretically different approaches to the most probable solution corresponding to the cryo-EM map. However, the result shows that the convergence might not be observed for proteins with complex conformational changes or with missing densities in cryo-EM map. In contrast, RDF structures deposited in the PDB can represent conformations that not only differ from the consensus obtained by flexible fitting but also from X-ray crystallography. Thus, this study emphasizes that a "consensus" achieved by the use of several automated flexible fitting approaches can provide a higher level of confidence in the modeled configurations. Following this protocol not only increases the confidence level of fitting, but also highlights protein regions with uncertain fitting. Hence, this protocol can lead to better interpretation of cryo-EM data.

Oligomerization state of photosynthetic core complexes is correlated with the dimerization affinity of a transmembrane helix

In the Rhodobacter (Rba.) species of photosynthetic purple bacteria, a single transmembrane α-helix, PufX, is found within the core complex, an essential photosynthetic macromolecular assembly that performs the absorption and the initial processing of light energy. Despite its structural simplicity, many unresolved questions surround PufX, the most important of which is its location within the photosynthetic core complex. One proposed placement of PufX is at the center of a core complex dimer, where two PufX helices associate in the membrane and form a homodimer. Inability for PufX of certain Rba. species to form a homodimer is thought to lead to monomeric core complexes. In the present study, we employ a combination of computational and experimental techniques to test the hypothesized homodimerization of PufX. We carry out a systematic investigation to measure the dimerization affinity of PufX from four Rba. species, Rba. blasticus , Rba. capsulatus , Rba. sphaeroides , and Rba. veldkampii , using a molecular dynamics-based free-energy method, as well as experimental TOXCAT assays. We found that the four PufX helices have substantially different dimerization affinities. Both computational and experimental techniques demonstrate that species with dimeric core complexes have PufX that can potentially form a homodimer, whereas the one species with monomeric core complexes has a PufX with little to no dimerization propensity. Our analysis of the helix-helix interface revealed a number of positions that may be important for PufX dimerization and the formation of a hydrogen-bond network between these GxxxG-containing helices. Our results suggest that the different oligomerization states of core complexes in various Rba. species can be attributed, among other factors, to the different propensity of its PufX helix to homodimerize.

Carotenoids are essential for normal levels of dimerisation of the RC-LH1-PufX core complex of Rhodobacter sphaeroides: characterisation of R-26 as a crtB (phytoene synthase) mutant

Carotenoids play important roles in photosynthesis where they are involved in light-harvesting, in photo-protection and in the assembly and structural stability of light-harvesting and reaction centre complexes. In order to examine the effects of carotenoids on the oligomeric state of the reaction centre-light-harvesting 1 -PufX (RC-LH1-PufX) core complex of Rhodobacter sphaeroides two carotenoid-less mutants, TC70 and R-26, were studied. Detergent fractionation showed that in the absence of carotenoids LH2 complexes do not assemble, as expected, but also that core complexes are predominantly found as monomers, although levels of the PufX polypeptide appeared to be unaffected. Analysis of R-26 membranes by electron microscopy and atomic force microscopy reveals arrays of hexagonally packed monomeric RC-LH1-PufX complexes. Transfer of the crtB gene encoding phytoene synthase to TC70 and R-26 restores the normal synthesis of carotenoids demonstrating that the R-26 mutant of Rba. sphaeroides harbours a mutation in crtB, among its other defects. The transconjugant TC70 and R-26 strains containing crtB had regained their ability to assemble wild-type levels of dimeric RC-LH1-PufX core complexes and normal energy transfer pathways were restored, demonstrating that carotenoids are essential for the normal assembly and function of both the LH2 and RC-LH1-PufX complexes in this bacterial photosystem.