Association of cytochrome c with membrane-bound cytochrome c oxidase proceeds parallel to the membrane rather than in bulk solution

Cryo-EM structure and kinetics reveal electron transfer by 2D diffusion of cytochrome c in the yeast III-IV respiratory supercomplex

Energy conversion in aerobic organisms involves an electron current from low-potential donors, such as NADH and succinate, to dioxygen through the membrane-bound respiratory chain. Electron transfer is coupled to transmembrane proton transport, which maintains the electrochemical proton gradient used to produce ATP and drive other cellular processes. Electrons are transferred from respiratory complexes III to IV (CIII and CIV) by water-soluble cytochrome (cyt.) c In Saccharomyces cerevisiae and some other organisms, these complexes assemble into larger CIII₂CIV_1/2 supercomplexes, the functional significance of which has remained enigmatic. In this work, we measured the kinetics of the S. cerevisiae supercomplex cyt. c-mediated QH₂:O₂ oxidoreductase activity under various conditions. The data indicate that the electronic link between CIII and CIV is confined to the surface of the supercomplex. Single-particle electron cryomicroscopy (cryo-EM) structures of the supercomplex with cyt. c show the positively charged cyt. c bound to either CIII or CIV or along a continuum of intermediate positions. Collectively, the structural and kinetic data indicate that cyt. c travels along a negatively charged patch on the supercomplex surface. Thus, rather than enhancing electron transfer rates by decreasing the distance that cyt. c must diffuse in three dimensions, formation of the CIII₂CIV_1/2 supercomplex facilitates electron transfer by two-dimensional (2D) diffusion of cyt. c This mechanism enables the CIII₂CIV_1/2 supercomplex to increase QH₂:O₂ oxidoreductase activity and suggests a possible regulatory role for supercomplex formation in the respiratory chain.

Structure and Mechanism of Respiratory III-IV Supercomplexes in Bioenergetic Membranes

In the final steps of energy conservation in aerobic organisms, free energy from electron transfer through the respiratory chain is transduced into a proton electrochemical gradient across a membrane. In mitochondria and many bacteria, reduction of the dioxygen electron acceptor is catalyzed by cytochrome c oxidase (complex IV), which receives electrons from cytochrome bc1 (complex III), via membrane-bound or water-soluble cytochrome c. These complexes function independently, but in many organisms they associate to form supercomplexes. Here, we review the structural features and the functional significance of the nonobligate III2IV1/2 Saccharomyces cerevisiae mitochondrial supercomplex as well as the obligate III2IV2 supercomplex from actinobacteria. The analysis is centered around the Q-cycle of complex III, proton uptake by CytcO, as well as mechanistic and structural solutions to the electronic link between complexes III and IV.

Peripheral Membrane Proteins: Promising Therapeutic Targets across Domains of Life

Membrane proteins can be classified into two main categories—integral and peripheral membrane proteins—depending on the nature of their membrane interaction. Peripheral membrane proteins are highly unique amphipathic proteins that interact with the membrane indirectly, using electrostatic or hydrophobic interactions, or directly, using hydrophobic tails or GPI-anchors. The nature of this interaction not only influences the location of the protein in the cell, but also the function. In addition to their unique relationship with the cell membrane, peripheral membrane proteins often play a key role in the development of human diseases such as African sleeping sickness, cancer, and atherosclerosis. This review will discuss the membrane interaction and role of periplasmic nitrate reductase, CymA, cytochrome c, alkaline phosphatase, ecto-5’-nucleotidase, acetylcholinesterase, alternative oxidase, type-II NADH dehydrogenase, and dihydroorotate dehydrogenase in certain diseases. The study of these proteins will give new insights into their function and structure, and may ultimately lead to ground-breaking advances in the treatment of severe diseases.

The conformational feasibility for the formation of reaching dimer in ASV and HIV integrase: a molecular dynamics study

Retroviral integrases are reported to form alternate dimer assemblies like the core-core dimer and reaching dimer. The core-core dimer is stabilized predominantly by an extensive interface between two catalytic core domains. The reaching dimer is stabilized by N-terminal domains that reach to form intermolecular interfaces with the other subunit's core and C-terminal domains (CTD), as well as CTD-CTD interactions. In this study, molecular dynamics (MD), Brownian dynamics (BD) simulations, and free energy analyses, were performed to elucidate determinants for the stability of the reaching dimer forms of full-length Avian Sarcoma Virus (ASV) and Human Immunodeficiency Virus (HIV) IN, and to examine the role of the C-tails (the last ~16-18 residues at the C-termini) in their structural dynamics. The dynamics of an HIV reaching dimer derived from small angle X-ray scattering and protein crosslinking data, was compared with the dynamics of a core-core dimer model derived from combining the crystal structures of two-domain fragments. The results showed that the core domains in the ASV reaching dimer express free dynamics, whereas those in the HIV reaching dimer are highly stable. BD simulations suggest a higher rate of association for the HIV core-core dimer than the reaching dimer. The predicted stability of these dimers was therefore ranked in the following order: ASV reaching dimer < HIV reaching dimer < composite core-core dimer. Analyses of MD trajectories have suggested residues that are critical for intermolecular contacts in each reaching dimer. Tests of these predictions and insights gained from these analyses could reveal a potential pathway for the association and dissociation of full-length IN multimers.

SDA 7: A modular and parallel implementation of the simulation of diffusional association software

The simulation of diffusional association (SDA) Brownian dynamics software package has been widely used in the study of biomacromolecular association. Initially developed to calculate bimolecular protein-protein association rate constants, it has since been extended to study electron transfer rates, to predict the structures of biomacromolecular complexes, to investigate the adsorption of proteins to inorganic surfaces, and to simulate the dynamics of large systems containing many biomacromolecular solutes, allowing the study of concentration-dependent effects. These extensions have led to a number of divergent versions of the software. In this article, we report the development of the latest version of the software (SDA 7). This release was developed to consolidate the existing codes into a single framework, while improving the parallelization of the code to better exploit modern multicore shared memory computer architectures. It is built using a modular object-oriented programming scheme, to allow for easy maintenance and extension of the software, and includes new features, such as adding flexible solute representations. We discuss a number of application examples, which describe some of the methods available in the release, and provide benchmarking data to demonstrate the parallel performance.

On calculation of the electrostatic potential of a phosphatidylinositol phosphate-containing phosphatidylcholine lipid membrane accounting for membrane dynamics

Many signaling events require the binding of cytoplasmic proteins to cell membranes by recognition of specific charged lipids, such as phosphoinositol-phosphates. As a model for a protein-membrane binding site, we consider one charged phosphoinositol phosphate (PtdIns(3)P) embedded in a phosphatidylcholine bilayer. As the protein-membrane binding is driven by electrostatic interactions, continuum solvent models require an accurate representation of the electrostatic potential of the phosphoinositol phosphate-containing membrane. We computed and analyzed the electrostatic potentials of snapshots taken at regular intervals from molecular dynamics simulations of the bilayer. We observe considerable variation in the electrostatic potential of the bilayer both along a single simulation and between simulations performed with the GAFF or CHARMM c36 force fields. However, we find that the choice of GAFF or CHARMM c36 parameters has little effect on the electrostatic potential of a given configuration of the bilayer with a PtdIns(3)P embedded in it. From our results, we propose a remedian averaging method for calculating the electrostatic potential of a membrane system that is suitable for simulations of protein-membrane binding with a continuum solvent model.

Single molecule activity measurements of cytochrome P450 oxidoreductase reveal the existence of two discrete functional states

Electron transfer between membrane spanning oxidoreductase enzymes controls vital metabolic processes. Here we studied for the first time with single molecule resolution the function of P450 oxidoreductase (POR), the canonical membrane spanning activator of all microsomal cytochrome P450 enzymes. Measurements and statistical analysis of individual catalytic turnover cycles shows POR to sample at least two major functional states. This phenotype may underlie regulatory interactions with different cytochromes P450 but to date has remained masked in bulk kinetics. To ensure that we measured the inherent behavior of POR, we reconstituted the full length POR in "native like" membrane patches, nanodiscs. Nanodisc reconstitution increased stability by ∼2-fold as compared to detergent solubilized POR and showed significantly increased activity at biologically relevant ionic strength conditions, highlighting the importance of studying POR function in a membrane environment. This assay paves the way for studying the function of additional membrane spanning oxidoreductases with single molecule resolution.

Studying the role of lipid rafts on protein receptor bindings with cellular automata

It is widely accepted that lipid rafts promote receptor clustering and thereby facilitate signaling transduction. The role of lipid rafts in inducing and promoting receptor accumulation within the cell membrane has been explored by several computational and experimental studies. However, it remains unclear whether lipid rafts influence the recruitment and binding of proteins from the cytosol as well. To provide an answer to this question a spatial membrane model has been developed based on cellular automata. Our results indicate that lipid rafts indeed influence protein receptor bindings. In particular processes with slow dissociation and binding kinetics are promoted by lipid rafts, whereas fast binding processes are slightly hampered. However, the impact depends on a variety of parameters, such as the size and mobility of the lipid rafts, the induced slow down of receptors within rafts, and also the dissociation and binding kinetics of the cytosolic proteins. Thus, for any individual signaling pathway the influence of lipid rafts on protein binding might be different. To facilitate analyzing this influence given a specific pathway, our approach has been generalized into LiRaM, a modeling and simulation tool for lipid rafts models.

Finding a needle in a haystack: the role of electrostatics in target lipid recognition by PH domains

Interactions between protein domains and lipid molecules play key roles in controlling cell membrane signalling and trafficking. The pleckstrin homology (PH) domain is one of the most widespread, binding specifically to phosphatidylinositol phosphates (PIPs) in cell membranes. PH domains must locate specific PIPs in the presence of a background of approximately 20% anionic lipids within the cytoplasmic leaflet of the plasma membrane. We investigate the mechanism of such recognition via a multiscale procedure combining Brownian dynamics (BD) and molecular dynamics (MD) simulations of the GRP1 PH domain interacting with phosphatidylinositol (3,4,5)-trisphosphate (PI(3,4,5)P₃). The interaction of GRP1-PH with PI(3,4,5)P₃ in a zwitterionic bilayer is compared with the interaction in bilayers containing different levels of anionic ‘decoy’ lipids. BD simulations reveal both translational and orientational electrostatic steering of the PH domain towards the PI(3,4,5)P₃-containing anionic bilayer surface. There is a payoff between non-PIP anionic lipids attracting the PH domain to the bilayer surface in a favourable orientation and their role as ‘decoys’, disrupting the interaction of GRP1-PH with the PI(3,4,5)P₃ molecule. Significantly, approximately 20% anionic lipid in the cytoplasmic leaflet of the bilayer is nearly optimal to both enhance orientational steering and to localise GRP1-PH proximal to the surface of the membrane without sacrificing its ability to locate PI(3,4,5)P₃ within the bilayer plane. Subsequent MD simulations reveal binding to PI(3,4,5)P₃, forming protein-phosphate contacts comparable to those in X-ray structures. These studies demonstrate a computational framework which addresses lipid recognition within a cell membrane environment, offering a link between structural and cell biological characterisation.

Cell signalling pathways are crucial for many biological processes including cell proliferation and survival. Signalling is governed by a complex network of interactions within the cell, and disruption of signalling can lead to a variety of human diseases. Often, a key event in the signalling cascade is the reversible recruitment of peripheral membrane proteins to the surface of the cell membrane, where they then bind to a specific lipid in order to perform their function. However, it is not clear how these proteins locate their target lipid in the complex multi-lipid environment of the plasma membrane. Here, we have used a combination of computational techniques to simulate the association of a signalling protein with the surface of the cell membrane. We demonstrate that the mechanism of membrane binding is dependent upon the lipid composition of the lipid bilayer, and the results show that orientational and positional steering of the protein is optimised when the anionic lipid content of our model membrane matches the physiological composition observed in cells.

Many-particle Brownian and Langevin Dynamics Simulations with the Brownmove package

Background

Brownian Dynamics (BD) is a coarse-grained implicit-solvent simulation method that is routinely used to investigate binary protein association dynamics, but due to its efficiency in handling large simulation volumes and particle numbers it is well suited to also describe many-protein scenarios as they often occur in biological cells.

Results

Here we introduce our "brownmove" simulation package which was designed to handle many-particle problems with varying particle numbers and allows for a very flexible definition of rigid and flexible protein and polymer models. Both a Brownian and a Langevin dynamics (LD) propagation scheme can be used and hydrodynamic interactions are treated efficiently with our recently introduced TEA-HI ansatz [Geyer, Winter, JCP 130 (2009) 114905]. With simulations of constrained polymers and flexible models of spherical proteins we demonstrate that it is crucial to include hydrodynamics when multi-bead models are used in BD or LD simulations. Only then both the translational and the rotational diffusion coefficients and the timescales of the internal dynamics can be reproduced correctly. In the third example project we show how constant density boundary conditions [Geyer et al, JCP 120 (2004) 4573] can be used to set up a non-equilibrium simulation of diffusional transport across an array of fixed obstacles. Finally, we demonstrate how the agglomeration dynamics of multiple particles with attractive patches can be analysed conveniently with the help of a dynamic interaction network.

Conclusions

Combining BD and LD propagation, fast hydrodynamics, a flexible protein model, and interfaces for "open" simulation settings, our freely available "brownmove" simulation package constitutes a new platform for coarse-grained many-particle simulations of biologically relevant diffusion and transport processes.

Analysis of Leigh syndrome mutations in the yeast SURF1 homolog reveals a new member of the cytochrome oxidase assembly factor family

Three missense SURF1 mutations identified in patients with Leigh syndrome (LS) were evaluated in the yeast homolog Shy1 protein. Introduction of two of the Leigh mutations, F(249)T and Y(344)D, in Shy1 failed to significantly attenuate the function of Shy1 in cytochrome c oxidase (CcO) biogenesis as seen with the human mutations. In contrast, a G(137)E substitution in Shy1 results in a nonfunctional protein conferring a CcO deficiency. The G(137)E Shy1 mutant phenocopied shy1Delta cells in impaired Cox1 hemylation and low mitochondrial copper. A genetic screen for allele-specific suppressors of the G(137)E Shy1 mutant revealed Coa2, Cox10, and a novel factor designated Coa4. Coa2 and Cox10 are previously characterized CcO assembly factors. Coa4 is a twin CX(9)C motif mitochondrial protein localized in the intermembrane space and associated with the inner membrane. Cells lacking Coa4 are depressed in CcO activity but show no impairment in Cox1 maturation or formation of the Shy1-stabilized Cox1 assembly intermediate. To glean insights into the functional role of Coa4 in CcO biogenesis, an unbiased suppressor screen of coa4Delta cells was conducted. Respiratory function of coa4Delta cells was restored by the overexpression of CYC1 encoding cytochrome c. Cyc1 is known to be important at an ill-defined step in the assembly and/or stability of CcO. This new link to Coa4 may begin to further elucidate the role of Cyc1 in CcO biogenesis.