Structure and Elasticity of Mitochondrial Membranes: A Molecular Dynamics Simulation Study

Mitochondria are known as the powerhouse of the cell because they produce energy in the form of adenosine triphosphate. They also have other crucial functions such as regulating apoptosis, calcium homeostasis, and reactive oxygen species production. To perform these diverse functions, mitochondria adopt specific structures and frequently undergo dynamic shape changes, indicating that their mechanical properties play an essential role in their functions. To gain a detailed understanding at the molecular level of the structure and mechanical properties of mitochondria, we carry out atomistic molecular dynamics simulations for three inner mitochondrial membranes and three outer mitochondrial membrane models. These models take into account variations in cardiolipin and cholesterol concentrations as well as the symmetry/asymmetry between the two leaflets. Our simulations allow us to calculate various structural quantities and the bending, twisting, and tilting elastic moduli of the membrane models. Our results indicate that the structures of the inner and outer mitochondrial membranes are quite similar and do not depend much on the variation in lipid compositions. However, the bending modulus of the membranes increases with increasing concentrations of cardiolipin or cholesterol but decreases with a membrane asymmetry. Notably, we found that the dipole potential of the membrane increases with an increasing cardiolipin concentration. Finally, possible roles of cardiolipin in regulating ion and proton currents and maintaining the cristate are discussed in some details.

FA Sliding as the Mechanism for the ANT1-Mediated Fatty Acid Anion Transport in Lipid Bilayers

Mitochondrial adenine nucleotide translocase (ANT) exchanges ADP for ATP to maintain energy production in the cell. Its protonophoric function in the presence of long-chain fatty acids (FA) is also recognized. Our previous results imply that proton/FA transport can be best described with the FA cycling model, in which protonated FA transports the proton to the mitochondrial matrix. The mechanism by which ANT1 transports FA anions back to the intermembrane space remains unclear. Using a combined approach involving measurements of the current through the planar lipid bilayers reconstituted with ANT1, site-directed mutagenesis and molecular dynamics simulations, we show that the FA anion is first attracted by positively charged arginines or lysines on the matrix side of ANT1 before moving along the positively charged protein-lipid interface and binding to R79, where it is protonated. We show that R79 is also critical for the competitive binding of ANT1 substrates (ADP and ATP) and inhibitors (carboxyatractyloside and bongkrekic acid). The binding sites are well conserved in mitochondrial SLC25 members, suggesting a general mechanism for transporting FA anions across the inner mitochondrial membrane.

Molecular docking and molecular dynamics simulation study on the toxicity mechanism of bongkrekic acid

BKA belongs to gram-negative brevibacterium. It can cause poisoning in humans or animals and can be fatal in severe cases. There are few investigations on toxic mechanisms of BKA because of foodborne factors. MD simulations were used to study the stability and intermolecular interactions of BKA and ANT complexes to reveal the mechanism of BKA in this paper. BKA blocked ANT protein translocation mainly through Van der Waals force, hydrophobic and hydrogen bonding interactions by the MD simulations. The conformational flexibility of the complex system during different simulation times indicated that BKA affected the conformational changes of ANT through strong interactions of hydrogen bonds with active domain residues Gln-93, Tyr-196, Arg-287 and Arg-245. The results of binding free energy, principal component analysis, hydrophobic interactions and root-mean-square fluctuation showed that the prominent binding force of Tyr-196 with C₂₆ of BKA was significant to the toxicity. The active site interactions analysis indicated that the essential positively charged polar amino acids which play a crucial role within the active site of the ANT protein undergo conformational changes with BKA as the branch point.

Structure, substrate binding, and symmetry of the mitochondrial ADP/ATP carrier in its matrix-open state

The mitochondrial ADP/ATP carrier (AAC) performs the first and last step in oxidative phosphorylation by exchanging ADP and ATP across the mitochondrial inner membrane. Its optimal function has been shown to be dependent on cardiolipins (CLs), unique phospholipids located almost exclusively in the mitochondrial membrane. In addition, AAC exhibits an enthralling threefold pseudosymmetry, a unique feature of members of the SLC25 family. Recently, its conformation poised for binding of ATP was solved by x-ray crystallography referred to as the matrix state. Binding of the substrate leads to conformational changes that export of ATP to the mitochondrial intermembrane space. In this contribution, we investigate the influence of CLs on the structure, substrate-binding properties, and structural symmetry of the matrix state, employing microsecond-scale molecular dynamics simulations. Our findings demonstrate that CLs play a minor stabilizing role on the AAC structure. The interdomain salt bridges and hydrogen bonds forming the cytoplasmic network and tyrosine braces, which ensure the integrity of the global AAC scaffold, highly benefit from the presence of CLs. Under these conditions, the carrier is found to be organized in a more compact structure in its interior, as revealed by analyses of the electrostatic potential, measure of the AAC cavity aperture, and the substrate-binding assays. Introducing a convenient structure-based symmetry metric, we quantified the structural threefold pseudosymmetry of AAC, not only for the crystallographic structure, but also for conformational states of the carrier explored in the molecular dynamics simulations. Our results suggest that CLs moderately contribute to preserve the pseudosymmetric structure of AAC.

Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design

We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.

Mitochondrial Uncoupling Proteins (UCP1-UCP3) and Adenine Nucleotide Translocase (ANT1) Enhance the Protonophoric Action of 2,4-Dinitrophenol in Mitochondria and Planar Bilayer Membranes

2,4-Dinitrophenol (DNP) is a classic uncoupler of oxidative phosphorylation in mitochondria which is still used in “diet pills”, despite its high toxicity and lack of antidotes. DNP increases the proton current through pure lipid membranes, similar to other chemical uncouplers. However, the molecular mechanism of its action in the mitochondria is far from being understood. The sensitivity of DNP’s uncoupling action in mitochondria to carboxyatractyloside, a specific inhibitor of adenine nucleotide translocase (ANT), suggests the involvement of ANT and probably other mitochondrial proton-transporting proteins in the DNP’s protonophoric activity. To test this hypothesis, we investigated the contribution of recombinant ANT1 and the uncoupling proteins UCP1-UCP3 to DNP-mediated proton leakage using the well-defined model of planar bilayer lipid membranes. All four proteins significantly enhanced the protonophoric effect of DNP. Notably, only long-chain free fatty acids were previously shown to be co-factors of UCPs and ANT1. Using site-directed mutagenesis and molecular dynamics simulations, we showed that arginine 79 of ANT1 is crucial for the DNP-mediated increase of membrane conductance, implying that this amino acid participates in DNP binding to ANT1.

Computational studies of the mitochondrial carrier family SLC25. Present status and future perspectives

The members of the mitochondrial carrier family, also known as solute carrier family 25 (SLC25), are transmembrane proteins involved in the translocation of a plethora of small molecules between the mitochondrial intermembrane space and the matrix. These transporters are characterized by three homologous domains structure and a transport mechanism that involves the transition between different conformations. Mutations in regions critical for these transporters’ function often cause several diseases, given the crucial role of these proteins in the mitochondrial homeostasis. Experimental studies can be problematic in the case of membrane proteins, in particular concerning the characterization of the structure–function relationships. For this reason, computational methods are often applied in order to develop new hypotheses or to support/explain experimental evidence. Here the computational analyses carried out on the SLC25 members are reviewed, describing the main techniques used and the outcome in terms of improved knowledge of the transport mechanism. Potential future applications on this protein family of more recent and advanced in silico methods are also suggested.

ANT1 Activation and Inhibition Patterns Support the Fatty Acid Cycling Mechanism for Proton Transport

Adenine nucleotide translocase (ANT) is a well-known mitochondrial exchanger of ATP against ADP. In contrast, few studies have shown that ANT also mediates proton transport across the inner mitochondrial membrane. The results of these studies are controversial and lead to different hypotheses about molecular transport mechanisms. We hypothesized that the H⁺-transport mediated by ANT and uncoupling proteins (UCP) has a similar regulation pattern and can be explained by the fatty acid cycling concept. The reconstitution of purified recombinant ANT1 in the planar lipid bilayers allowed us to measure the membrane current after the direct application of transmembrane potential ΔΨ, which would correspond to the mitochondrial states III and IV. Experimental results reveal that ANT1 does not contribute to a basal proton leak. Instead, it mediates H⁺ transport only in the presence of long-chain fatty acids (FA), as already known for UCPs. It depends on FA chain length and saturation, implying that FA’s transport is confined to the lipid-protein interface. Purine nucleotides with the preference for ATP and ADP inhibited H⁺ transport. Specific inhibitors of ATP/ADP transport, carboxyatractyloside or bongkrekic acid, also decreased proton transport. The H⁺ turnover number was calculated based on ANT1 concentration determined by fluorescence correlation spectroscopy and is equal to 14.6 ± 2.5 s⁻¹. Molecular dynamic simulations revealed a large positively charged area at the protein/lipid interface that might facilitate FA anion’s transport across the membrane. ANT’s dual function—ADP/ATP and H⁺ transport in the presence of FA—may be important for the regulation of mitochondrial membrane potential and thus for potential-dependent processes in mitochondria. Moreover, the expansion of proton-transport modulating drug targets to ANT1 may improve the therapy of obesity, cancer, steatosis, cardiovascular and neurodegenerative diseases.

Molecular Dynamics Simulations of Mitochondrial Uncoupling Protein 2

Molecular dynamics (MD) simulations of uncoupling proteins (UCP), a class of transmembrane proteins relevant for proton transport across inner mitochondrial membranes, represent a complicated task due to the lack of available structural data. In this work, we use a combination of homology modelling and subsequent microsecond molecular dynamics simulations of UCP2 in the DOPC phospholipid bilayer, starting from the structure of the mitochondrial ATP/ADP carrier (ANT) as a template. We show that this protocol leads to a structure that is impermeable to water, in contrast to MD simulations of UCP2 structures based on the experimental NMR structure. We also show that ATP binding in the UCP2 cavity is tight in the homology modelled structure of UCP2 in agreement with experimental observations. Finally, we corroborate our results with conductance measurements in model membranes, which further suggest that the UCP2 structure modeled from ANT protein possesses additional key functional elements, such as a fatty acid-binding site at the R60 region of the protein, directly related to the proton transport mechanism across inner mitochondrial membranes.

Function-related asymmetry of the specific cardiolipin binding sites on the mitochondrial ADP/ATP carrier

The ADP/ATP carrier (AAC) transports matrix ATP and cytosolic ADP across the inner mitochondrial membrane (IMM). It is well known that cardiolipin (CL) plays an important role in regulating the function of AAC, yet the underlying mechanism still remains elusive. AAC is composed of three homologous domains, and three specific CL binding sites are located at the domain-domain interfaces near the matrix side. Here we report an in-depth investigation on the dynamic properties of the bound CL within the three specific sites through all-atom molecular dynamics simulations of up to 13 μs in total. Our results highlight the importance of the basic and polar residues in CL binding. The basic residues from the linker helix and/or the [Y/W/F][K/R]G motif enable the bound CL to form an intra-domain binding mode, and the canonical inter-domain binding mode only forms when these basic residues are occupied by an additional phospholipid. Of special significance, differences in the basic and polar residues lead to remarkable asymmetry among the three specific CL binding sites. We found that the bound CL at the interface of domains 2 and 3 predominantly adopts inter-domain binding mode, while CLs at the other two sites have much more intra-domain populations. This is consistent with the asymmetric crystal structure of the matrix state (m-state) AAC which implies an asymmetric transport mechanism. The dynamic equilibrium between the inter-domain and intra-domain binding modes observed in our simulations could be highly important for the bound CLs to adapt to the movements during state transitions.

Functional Oligomeric Forms of Uncoupling Protein 2: Strong Evidence for Asymmetry in Protein and Lipid Bilayer Systems

Stoichiometry of uncoupling proteins (UCPs) and their coexistence as functional monomeric and associated forms in lipid membranes remain intriguing open questions. In this study, tertiary and quaternary structures of UCP2 were analyzed experimentally and through molecular dynamics (MD) simulations. UCP2 was overexpressed in the inner membrane of Escherichia coli, then purified and reconstituted in lipid vesicles. Structure and proton transport function of UCP2 were characterized by circular dichroism (CD) spectroscopy and fluorescence methods. Findings suggest a tetrameric functional form for UCP2. MD simulations conclude that tetrameric UCP2 is a dimer of dimers, is more stable than its monomeric and dimeric forms, is asymmetrical and induces asymmetry in the membrane's lipid structure, and a biphasic on-off switch between the dimeric units is its possible mode of transport. MD simulations also show that the water density inside the UCP2 monomer is asymmetric, with the cytoplasmic side having a higher water density and a wider radius. In contrast, the structurally comparable adenosine 5'-diphosphate (ADP)/adenosine 5'-triphosphate (ATP) carrier (AAC1) did not form tetramers, implying that tetramerization cannot be generalized to all mitochondrial carriers.