Insight into the Properties of Cardiolipin Containing Bilayers from Molecular Dynamics Simulations, Using a Hybrid All-Atom/United-Atom Force Field

Curvature sensing lipid dynamics in a mitochondrial inner membrane model

Membrane curvature is essential for many cellular structures and processes, and factors such as leaflet asymmetry, lipid composition, and proteins all play important roles. Cardiolipin is the signature lipid of mitochondrial membranes and is essential for maintaining the highly curved shapes of the inner mitochondrial membrane (IMM) and the spatial arrangement of membrane proteins. In this study, we investigate the partitioning behavior of various lipids present in the IMM using coarse-grained molecular dynamics simulations. This study explores curved bilayer systems containing phosphatidylcholine (PC), phosphatidylethanolamine (PE), and cardiolipin (CDL) in binary and ternary component mixtures. Curvature properties such as mean and Gaussian curvatures, as well as the distribution of lipids into the various curved regions of the cristae models, are quantified. Overall, this work represents an advance beyond previous studies on lipid curvature sensing by simulating these systems in a geometry that has the morphological features and scales of curvature consistent with regions of the IMM. We find that CDL has a stronger preference for accumulating in regions of negative curvature than PE lipids, in agreement with previous results. Furthermore, we find lipid partitioning propensity is dominated by sensitivity to mean curvature, while there is a weaker correlation with Gaussian curvature.

Domain motions, dimerization, and membrane interactions of the murine guanylate binding protein 2

Guanylate-binding proteins (GBPs) are a group of GTPases that are induced by interferon-[Formula: see text] and are crucial components of cell-autonomous immunity against intracellular pathogens. Here, we examine murine GBP2 (mGBP2), which we have previously shown to be an essential effector protein for the control of Toxoplasma gondii replication, with its recruitment through the membrane of the parasitophorous vacuole and its involvement in the destruction of this membrane likely playing a role. The overall aim of our work is to provide a molecular-level understanding of the mutual influences of mGBP2 and the parasitophorous vacuole membrane. To this end, we performed lipid-binding assays which revealed that mGBP2 has a particular affinity for cardiolipin. This observation was confirmed by fluorescence microscopy using giant unilamellar vesicles of different lipid compositions. To obtain an understanding of the protein dynamics and how this is affected by GTP binding, mGBP2 dimerization, and membrane binding, assuming that each of these steps are relevant for the function of the protein, we carried out standard as well as replica exchange molecular dynamics simulations with an accumulated simulation time of more than 30 μs. The main findings from these simulations are that mGBP2 features a large-scale hinge motion in its M/E domain, which is present in each of the studied protein states. When bound to a cardiolipin-containing membrane, this hinge motion is particularly pronounced, leading to an up and down motion of the M/E domain on the membrane, which did not occur on a membrane without cardiolipin. Our prognosis is that this up and down motion has the potential to destroy the membrane following the formation of supramolecular mGBP2 complexes on the membrane surface.

Cardiolipin prevents pore formation in phosphatidylglycerol bacterial membrane models

Several antimicrobial peptides, including magainin and the human cathelicidin LL-37, act by forming pores in bacterial membranes. Bacteria such as Staphylococcus aureus modify their membrane's cardiolipin composition to resist such types of perturbations that compromise their membrane stability. Here, we used molecular dynamic simulations to quantify the role of cardiolipin on the formation of pores in simple bacterial-like membrane models composed of phosphatidylglycerol and cardiolipin mixtures. Cardiolipin modified the structure and ordering of the lipid bilayer, making it less susceptible to mechanical changes. Accordingly, the free-energy barrier for the formation of a transmembrane pore and its kinetic instability augmented by increasing the cardiolipin concentration. This is attributed to the unfavorable positioning of cardiolipin near the formed pore, due to its small polar head and bulky hydrophobic body. Overall, our study demonstrates how cardiolipin prevents membrane-pore formation and this constitutes a plausible mechanism used by bacteria to act against stress perturbations and, thereby, gain resistance to antimicrobial agents.

Predicting the Structure and Dynamics of Membrane Protein GerAB from Bacillus subtilis

Bacillus subtilis forms dormant spores upon nutrient depletion. Germinant receptors (GRs) in spore’s inner membrane respond to ligands such as L-alanine, and trigger spore germination. In B. subtilis spores, GerA is the major GR, and has three subunits, GerAA, GerAB, and GerAC. L-Alanine activation of GerA requires all three subunits, but which binds L-alanine is unknown. To date, how GRs trigger germination is unknown, in particular due to lack of detailed structural information about B subunits. Using homology modelling with molecular dynamics (MD) simulations, we present structural predictions for the integral membrane protein GerAB. These predictions indicate that GerAB is an α-helical transmembrane protein containing a water channel. The MD simulations with free L-alanine show that alanine binds transiently to specific sites on GerAB. These results provide a starting point for unraveling the mechanism of L-alanine mediated signaling by GerAB, which may facilitate early events in spore germination.

Unraveling the Role of Drug-Lipid Interactions in NSAIDs-Induced Cardiotoxicity

Cardiovascular (CV) toxicity is nowadays recognized as a class effect of non-aspirin nonsteroidal anti-inflammatory drugs (NSAIDs). However, their mechanisms of cardiotoxicity are not yet well understood, since different compounds with similar action mechanisms exhibit distinct cardiotoxicity. For instance, diclofenac (DIC) is among the most cardiotoxic compounds, while naproxen (NAP) is associated with low CV risk. In this sense, this study aimed to unravel the role of drug-lipid interactions in NSAIDs-induced cardiotoxicity. For that, DIC and NAP interactions with lipid bilayers as model systems of cell and mitochondrial membranes were characterized by derivative spectrophotometry, fluorometric leakage assays, and synchrotron X-ray scattering. Both DIC and NAP were found to have the ability to permeabilize the membrane models, as well as to alter the bilayers’ structure. The NSAIDs-induced modifications were dependent on the lipid composition of the membrane model, the three-dimensional structure of the drug, as well as the drug:lipid molar ratio tested. Altogether, this work supports the hypothesis that NSAIDs-lipid interactions, in particular at the mitochondrial level, may be another key step among the mechanisms underlying NSAIDs-induced cardiotoxicity.

Cardiolipin-Dependent Properties of Model Mitochondrial Membranes from Molecular Simulations

Cardiolipin is an anionic lipid found in the mitochondrial membranes of eukaryotes ranging from unicellular microorganisms to metazoans. This unique lipid contributes to various mitochondrial functions, including metabolism, mitochondrial membrane fusion and/or fission dynamics, and apoptosis. However, differences in cardiolipin content between the two mitochondrial membranes, as well as dynamic fluctuations in cardiolipin content in response to stimuli and cellular signaling events, raise questions about how cardiolipin concentration affects mitochondrial membrane structure and dynamics. Although cardiolipin’s structural and dynamic roles have been extensively studied in binary mixtures with other phospholipids, the biophysical properties of cardiolipin in higher number lipid mixtures are still not well resolved. Here, we used molecular dynamics simulations to investigate the cardiolipin-dependent properties of ternary lipid bilayer systems that mimic the major components of mitochondrial membranes. We found that changes to cardiolipin concentration only resulted in minor changes to bilayer structural features but that the lipid diffusion was significantly affected by those alterations. We also found that cardiolipin position along the bilayer surfaces correlated to negative curvature deflections, consistent with the induction of negative curvature stress in the membrane monolayers. This work contributes to a foundational understanding of the role of cardiolipin in altering the properties in ternary lipid mixtures composed of the major mitochondrial phospholipids, providing much-needed insights to help understand how cardiolipin concentration modulates the biophysical properties of mitochondrial membranes.

Molecular Dynamics Analysis of Cardiolipin and Monolysocardiolipin on Bilayer Properties

The mitochondrial lipid cardiolipin (CL) contributes to the spatial protein organization and morphological character of the inner mitochondrial membrane. Monolysocardiolipin (MLCL), an intermediate species in the CL remodeling pathway, is enriched in the multisystem disease Barth syndrome. Despite the medical relevance of MLCL, a detailed molecular description that elucidates the structural and dynamic differences between CL and MLCL has not been conducted. To this end, we performed comparative atomistic molecular dynamics studies on bilayers consisting of pure CL or MLCL to elucidate similarities and differences in their molecular and bulk bilayer properties. We describe differential headgroup dynamics and hydrogen bonding patterns between the CL variants and show an increased cohesiveness of MLCL's solvent interfacial region, which may have implications for protein interactions. Finally, using the coarse-grained Martini model, we show that substitution of MLCL for CL in bilayers mimicking mitochondrial composition induces drastic differences in bilayer mechanical properties and curvature-dependent partitioning behavior. Together, the results of this work reveal differences between CL and MLCL at the molecular and mesoscopic levels that may underpin the pathomechanisms of defects in cardiolipin remodeling.

Computational Modeling of Realistic Cell Membranes

Cell membranes contain a large variety of lipid types and are crowded with proteins, endowing them with the plasticity needed to fulfill their key roles in cell functioning. The compositional complexity of cellular membranes gives rise to a heterogeneous lateral organization, which is still poorly understood. Computational models, in particular molecular dynamics simulations and related techniques, have provided important insight into the organizational principles of cell membranes over the past decades. Now, we are witnessing a transition from simulations of simpler membrane models to multicomponent systems, culminating in realistic models of an increasing variety of cell types and organelles. Here, we review the state of the art in the field of realistic membrane simulations and discuss the current limitations and challenges ahead.

Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.

Mitochondrial ADP/ATP Carrier in Dodecylphosphocholine Binds Cardiolipins with Non-native Affinity

Biophysical investigation of membrane proteins generally requires their extraction from native sources using detergents, a step that can lead, possibly irreversibly, to protein denaturation. The propensity of dodecylphosphocholine (DPC), a detergent widely utilized in NMR studies of membrane proteins, to distort their structure has been the subject of much controversy. It has been recently proposed that the binding specificity of the yeast mitochondrial ADP/ATP carrier (yAAC3) toward cardiolipins is preserved in DPC, thereby suggesting that DPC is a suitable environment in which to study membrane proteins. In this communication, we used all-atom molecular dynamics simulations to investigate the specific binding of cardiolipins to yAAC3. Our data demonstrate that the interaction interface observed in a native-like environment differs markedly from that inferred from an NMR investigation in DPC, implying that in this detergent, the protein structure is distorted. We further investigated yAAC3 solubilized in DPC and in the milder dodecylmaltoside with thermal-shift assays. The loss of thermal transition observed in DPC confirms that the protein is no longer properly folded in this environment.

Buckling Under Pressure: Curvature-Based Lipid Segregation and Stability Modulation in Cardiolipin-Containing Bilayers

Mitochondrial metabolic function is affected by the morphology and protein organization of the mitochondrial inner membrane. Cardiolipin (CL) is a unique tetra-acyl lipid that is involved in the maintenance of the highly curved shape of the mitochondrial inner membrane as well as spatial organization of the proteins necessary for respiration and oxidative phosphorylation. Cardiolipin has been suggested to self-organize into lipid domains due to its inverted conical molecular geometry, though the driving forces for this organization are not fully understood. In this work, we use coarse-grained molecular dynamics simulations to study the mechanical properties and lipid dynamics in heterogeneous bilayers both with and without CL, as a function of membrane curvature. We find that incorporation of CL increases bilayer deformability and that CL becomes highly enriched in regions of high negative curvature. We further show that another mitochondrial inverted conical lipid, phosphatidylethanolamine (PE), does not partition or increase the deformability of the membrane in a significant manner. Therefore, CL appears to possess some unique characteristics that cannot be inferred simply from molecular geometry considerations.

Phosphoethanolamine addition to the Heptose I of the Lipopolysaccharide modifies the inner core structure and has an impact on the binding of Polymyxin B to the Escherichia coli outer membrane

Phosphoethanolamine (pEtN) decoration of E. coli Lipopolysaccharide (LPS) provides resistance to the antimicrobial Polymyxin B (PolB). While EptA and EptB enzymes catalyze the addition of pEtN to the Lipid A and Kdo (pEtN-Kdo-Lipid A), EptC catalyzes the pEtN addition to the Heptose I (pEtN-Hept_I). In this study, we investigated the contribution of pEtN-Hept_I to PolB resistance using eptA/eptB and eptC deficient E. coli K12 and its wild-type parent strains. These mutations were shown to decrease the antimicrobial activity of PolB on cells grown under pEtN-addition inducing conditions. Furthermore, the 1-N-phenylnapthylamine uptake assay revealed that in vivo PolB has a reduced OM-permeabilizing activity on the ΔeptA/eptB strain compared with the ΔeptC strain. In vitro, the changes in size and zeta potential of LPS-vesicles indicate that pEtN-Hept_I reduce the PolB binding, but in a minor extent than pEtN-Kdo-Lipid A. Molecular dynamics analysis revealed the structural basis of the PolB resistance promoted by pEtN-Hept_I, which generate a new hydrogen-bonding networks and a denser inner core region. Altogether, the experimental and theoretical assays shown herein indicate that pEtN-Hept_I addition promote an LPS conformational rearrangement, that could act as a shield by hindering the accession of PolB to inner LPS-targets moieties.

Mutations at the Qo Site of the Cytochrome bc1 Complex Strongly Affect Oxygen Binding

The homodimeric bc₁ protein complex is embedded in membranes of mitochondria and photosynthetic bacteria, where it transports protons across the membrane to maintain an electrostatic potential used to drive ATP synthesis as part of the respiratory or photosynthetic pathways. The reaction cycle of the bc₁ complex is driven by series of redox processes involving substrate molecules from the membrane, but occasional side reactions between an intermediate semiquinone substrate and molecular oxygen are suspected to be a source of toxic superoxide, which is believed to be a factor in aging. The present investigation employs molecular dynamics simulations to study the effect of mutations in the Q_o binding sites of the bc₁ complex on the ability of oxygen molecules to migrate to and bind at various locations within the complex. It is found that the mutations strongly affect the ability of oxygen to bind at the Q_o sites, and moreover, different behavior of the two monomers of the bc₁ complex is observed. The conformational differences at the Q_o sites of the two monomers are studied in detail and discussed. The anionic form of semiquinone was identified as leading to the greatest opportunity for side reactions with oxygen.

Mechanism of the Primary Charge Transfer Reaction in the Cytochrome bc1 Complex

The bc₁ complex is a critical enzyme for the ATP production in photosynthesis and cellular respiration. Its biochemical function relies on the so-called Q-cycle, which is well established and operates via quinol substrates that bind inside the protein complex. Despite decades of research, the quinol-protein interaction, which initiates the Q-cycle, has not yet been completely described. Furthermore, the initial charge transfer reactions of the Q-cycle lack a physical description. The present investigation utilizes classical molecular dynamics simulations in tandem with quantum density functional theory calculations, to provide a complete and consistent quantitative description of the primary events that occur within the bc₁ complex upon quinol binding. In particular, the electron and proton transfer reactions that trigger the Q-cycle in the bc₁ complex from Rhodobacter capsulatus are studied. The coupled nature of these charge transfer reactions was revealed by obtaining the transition energy path connecting configurations of the Q_o-site prior and after the transfers. The analysis of orbitals and partial charge distribution of the different states of the Q_o-site has further supported the conclusion. Finally, key structural elements of the bc₁ complex that trigger the charge transfer reactions were established, manifesting the importance of the environment in the process, which is furthermore evidenced by free energy calculations.

Spontaneous Binding of Molecular Oxygen at the Qo-Site of the bc1 Complex Could Stimulate Superoxide Formation

A key part of the respiratory and photosynthetic pathways is the bc1 protein complex embedded in the inner membrane of mitochondria and the plasma membrane of photosynthetic bacteria. The protein complex pumps protons across the membrane to maintain an electrostatic potential, which is in turn used to drive ATP synthesis. This molecular machinery, however, is suspected to be a source of superoxide, which is toxic to the cell, even in minuscular quantities, and believed to be a factor in aging. Through molecular dynamics simulations, we investigate here the migration of molecular oxygen in the bc1 complex in order to identify possible reaction sites that could lead to superoxide formation. It is found, in particular, that oxygen penetrates spontaneously the Qo binding site of the bc1 complex in the presence of an intermediate semiquinone radical, thus making the Qo-site a strong candidate for being a center of superoxide production.

Role of charged lipids in membrane structures - Insight given by simulations

Lipids and proteins are the main components of cell membranes. It is becoming increasingly clear that lipids, in addition to providing an environment for proteins to work in, are in many cases also able to modulate the structure and function of those proteins. Particularly charged lipids such as phosphatidylinositols and phosphatidylserines are involved in several examples of such effects. Molecular dynamics simulations have proved an invaluable tool in exploring these aspects. This so-called computational microscope can provide both complementing explanations for the experimental results and guide experiments to fruitful directions. In this paper, we review studies that have utilized molecular dynamics simulations to unravel the roles of charged lipids in membrane structures. We focus on lipids as active constituents of the membranes, affecting both general membrane properties as well as non-lipid membrane components, mainly proteins. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

The ionization properties of cardiolipin and its variants in model bilayers

The anionic phospholipid cardiolipin has an unusual dimeric structure with a two-phosphate headgroup and four acyl chains. Cardiolipin is present in energy-transducing membranes that maintain electrochemical gradients, including most bacterial plasma membranes and the mitochondrial inner membrane, where it mediates respiratory complex assembly and activation, among many other roles. Dysfunctional biogenesis of cardiolipin is implicated in the pathogenesis of several diseases including Barth syndrome. Because cardiolipin is a dominant anionic lipid in energy-conserving membranes, its headgroup is a major contributor to surface charge density and the bilayer electrostatic profile. However, the proton dissociation behavior of its headgroup remains controversial. In one model, the pKa values of the phosphates differ by several units and the headgroup exists as a monoanion at physiological pH. In another model, both phosphates ionize as strong acids with low pKa values and the headgroup exists in dianionic form at physiological pH. Using independent electrokinetic and spectroscopic approaches, coupled with analysis using Gouy-Chapman-Stern formalism, we have analyzed the ionization properties of cardiolipin within biologically relevant lipid bilayer model systems. We show that both phosphates of the cardiolipin headgroup show strong ionization behavior with low pKa values. Moreover, cardiolipin variants lacking structural features proposed to be required to maintain disparate pKa values--namely the secondary hydroxyl on the central glycerol or a full complement of four acyl chains--were shown to have ionization behavior identical to intact cardiolipin. Hence, these results indicate that within the physiological pH range, the cardiolipin headgroup is fully ionized as a dianion. We discuss the implications of these results with respect to the role of cardiolipin in defining membrane surface potential, activating respiratory complexes, and modulating membrane curvature.

Binding of methylene blue onto Langmuir monolayers representing cell membranes may explain its efficiency as photosensitizer in photodynamic therapy

We provide evidence for the electrostatic interactions between the cationic photosensitizer methylene blue (MB) and cell membrane models represented by neat and mixed Langmuir monolayers of dioleylphosphatidylcholine (DOPC) and 1,1',2,2'-tetraoleoylcardiolipin (CL). From surface pressure measurements, MB was found to adsorb strongly and expand CL-containing monolayers, while it caused an apparent decreasing in molecular area on neat DOPC monolayer. The binding site of MB could be inferred from data with the surface-specific polarization-modulated infrared reflection-absorption spectroscopy (PM-IRRAS) technique, where changes induced by MB were observed in the vibrational modes of the phosphate groups of both CL and DOPC. The incorporation of MB also affected the carbonyl groups and the packing of the alkyl chains, thus indicating that MB binding site favors singlet oxygen generation close to the double bonds in the alkyl chains, an important requirement for photodynamic efficiency. Significantly, the data presented here demonstrate that MB may act in membranes composed by PCs, such as mammalian plasma membranes, and in those containing CL, as in bacterial and inner mitochondrial membranes.

Structural and mechanical properties of cardiolipin lipid bilayers determined using neutron spin echo, small angle neutron and X-ray scattering, and molecular dynamics simulations

The detailed structural and mechanical properties of a tetraoleoyl cardiolipin (TOCL) bilayer were determined using neutron spin echo (NSE) spectroscopy, small angle neutron and X-ray scattering (SANS and SAXS, respectively), and molecular dynamics (MD) simulations. We used MD simulations to develop a scattering density profile (SDP) model, which was then utilized to jointly refine SANS and SAXS data. In addition to commonly reported lipid bilayer structural parameters, component distributions were obtained, including the volume probability, electron density and neutron scattering length density. Of note, the distance between electron density maxima DHH (39.4 Å) and the hydrocarbon chain thickness 2DC (29.1 Å) of TOCL bilayers were both found to be larger than the corresponding values for dioleoyl phosphatidylcholine (DOPC) bilayers. Conversely, TOCL bilayers have a smaller overall bilayer thickness DB (36.7 Å), primarily due to their smaller headgroup volume per phosphate. SDP analysis yielded a lipid area of 129.8 Å(2), indicating that the cross-sectional area per oleoyl chain in TOCL bilayers (i.e., 32.5 Å(2)) is smaller than that for DOPC bilayers. Multiple sets of MD simulations were performed with the lipid area constrained at different values. The calculated surface tension versus lipid area resulted in a lateral area compressibility modulus KA of 342 mN m(-1), which is slightly larger compared to DOPC bilayers. Model free comparison to experimental scattering data revealed the best simulated TOCL bilayer from which detailed molecular interactions were determined. Specifically, Na(+) cations were found to interact most strongly with the glycerol hydroxyl linkage, followed by the phosphate and backbone carbonyl oxygens. Inter- and intra-lipid interactions were facilitated by hydrogen bonding between the glycerol hydroxyl and phosphate oxygen, but not with the backbone carbonyl. Finally, analysis of the intermediate scattering functions from NSE spectroscopy measurements of TOCL bilayers yielded a bending modulus KC of 1.06 × 10(-19) J, which was larger than that observed in DOPC bilayers. Our results show the physicochemical properties of cardiolin bilayers that may be important in explaining their functionality in the inner mitochondrial membrane.

Identification of ubiquinol binding motifs at the Qo-site of the cytochrome bc1 complex

Enzymes of the bc₁ complex family power the biosphere through their central role in respiration and photosynthesis. These enzymes couple the oxidation of quinol molecules by cytochrome c to the transfer of protons across the membrane, to generate a proton-motive force that drives ATP synthesis. Key for the function of the bc₁ complex is the initial redox process that involves a bifurcated electron transfer in which the two electrons from a quinol substrate are passed to different electron acceptors in the bc₁ complex. The electron transfer is coupled to proton transfer. The overall mechanism of quinol oxidation by the bc₁ complex is well enough characterized to allow exploration at the atomistic level, but details are still highly controversial. The controversy stems from the uncertain binding motifs of quinol at the so-called Q_o active site of the bc₁ complex. Here we employ a combination of classical all atom molecular dynamics and quantum chemical calculations to reveal the binding modes of quinol at the Q_o-site of the bc₁ complex from Rhodobacter capsulatus. The calculations suggest a novel configuration of amino acid residues responsible for quinol binding and support a mechanism for proton-coupled electron transfer from quinol to iron–sulfur cluster through a bridging hydrogen bond from histidine that stabilizes the reaction complex.

CHARMM-GUI Membrane Builder toward realistic biological membrane simulations

CHARMM-GUI Membrane Builder, http://www.charmm-gui.org/input/membrane, is a web-based user interface designed to interactively build all-atom protein/membrane or membrane-only systems for molecular dynamics simulations through an automated optimized process. In this work, we describe the new features and major improvements in Membrane Builder that allow users to robustly build realistic biological membrane systems, including (1) addition of new lipid types, such as phosphoinositides, cardiolipin (CL), sphingolipids, bacterial lipids, and ergosterol, yielding more than 180 lipid types, (2) enhanced building procedure for lipid packing around protein, (3) reliable algorithm to detect lipid tail penetration to ring structures and protein surface, (4) distance-based algorithm for faster initial ion displacement, (5) CHARMM inputs for P21 image transformation, and (6) NAMD equilibration and production inputs. The robustness of these new features is illustrated by building and simulating a membrane model of the polar and septal regions of E. coli membrane, which contains five lipid types: CL lipids with two types of acyl chains and phosphatidylethanolamine lipids with three types of acyl chains. It is our hope that CHARMM-GUI Membrane Builder becomes a useful tool for simulation studies to better understand the structure and dynamics of proteins and lipids in realistic biological membrane environments.

Biophysical significance of the inner mitochondrial membrane structure on the electrochemical potential of mitochondria

The available literature supports the hypothesis that the morphology of the inner mitochondrial membrane is regulated by different energy states, that the three-dimensional morphology of cristae is dynamic, and that both are related to biochemical function. Examination of the correlation between the inner mitochondrial membrane (IMM) structure and mitochondrial energetic function is critical to an understanding of the links between mesoscale morphology and function in progressive mitochondrial dysfunction such as aging, neurodegeneration, and disease. To investigate this relationship, we develop a model to examine the effects of three-dimensional IMM morphology on the electrochemical potential of mitochondria. The two-dimensional axisymmetric finite element method is used to simulate mitochondrial electric potential and proton concentration distribution. This simulation model demonstrates that the proton motive force (Δp) produced on the membranes of cristae can be higher than that on the inner boundary membrane. The model also shows that high proton concentration in cristae can be induced by the morphology-dependent electric potential gradient along the outer side of the IMM. Furthermore, simulation results show that a high Δp is induced by the large surface-to-volume ratio of an individual crista, whereas a high capacity for ATP synthesis can primarily be achieved by increasing the surface area of an individual crista. The mathematical model presented here provides compelling support for the idea that morphology at the mesoscale is a significant driver of mitochondrial function.