The role of cardiolipin concentration and acyl chain composition on mitochondrial inner membrane molecular organization and function

Cardiolipin (CL) is a key phospholipid of the mitochondria. A loss of CL content and remodeling of CL's acyl chains is observed in several pathologies. Strong shifts in CL concentration and acyl chain composition would presumably disrupt mitochondrial inner membrane biophysical organization. However, it remains unclear in the literature as to which is the key regulator of mitochondrial membrane biophysical properties. We review the literature to discriminate the effects of CL concentration and acyl chain composition on mitochondrial membrane organization. A widely applicable theme emerges across several pathologies, including cardiovascular diseases, diabetes, Barth syndrome, and neurodegenerative ailments. The loss of CL, often accompanied by increased levels of lyso-CLs, impairs mitochondrial inner membrane organization. Modest remodeling of CL acyl chains is not a major driver of impairments and only in cases of extreme remodeling is there an influence on membrane properties.     

Dynamics of the Ethanolamine Glycerophospholipid Remodeling Network

Acyl chain remodeling in lipids is a critical biochemical process that plays a central role in disease. However, remodeling remains poorly understood, despite massive increases in lipidomic data. In this work, we determine the dynamic network of ethanolamine glycerophospholipid (PE) remodeling, using data from pulse-chase experiments and a novel bioinformatic network inference approach. The model uses a set of ordinary differential equations based on the assumptions that (1) sn1 and sn2 acyl positions are independently remodeled; (2) remodeling reaction rates are constant over time; and (3) acyl donor concentrations are constant. We use a novel fast and accurate two-step algorithm to automatically infer model parameters and their values. This is the first such method applicable to dynamic phospholipid lipidomic data. Our inference procedure closely fits experimental measurements and shows strong cross-validation across six independent experiments with distinct deuterium-labeled PE precursors, demonstrating the validity of our assumptions. In constrast, fits of randomized data or fits using random model parameters are worse. A key outcome is that we are able to robustly distinguish deacylation and reacylation kinetics of individual acyl chain types at the sn1 and sn2 positions, explaining the established prevalence of saturated and unsaturated chains in the respective positions. The present study thus demonstrates that dynamic acyl chain remodeling processes can be reliably determined from dynamic lipidomic data.     

Unremodeled and Remodeled Cardiolipin Are Functionally Indistinguishable in Yeast

After biosynthesis, an evolutionarily conserved acyl chain remodeling process generates a final highly homogeneous and yet tissue-specific molecular form of the mitochondrial lipid cardiolipin. Hence, cardiolipin molecules in different organisms, and even different tissues within the same organism, contain a distinct collection of attached acyl chains. This observation is the basis for the widely accepted paradigm that the acyl chain composition of cardiolipin is matched to the unique mitochondrial demands of a tissue. For this hypothesis to be correct, cardiolipin molecules with different acyl chain compositions should have distinct functional capacities, and cardiolipin that has been remodeled should promote cardiolipin-dependent mitochondrial processes better than its unremodeled form. However, functional disparities between different molecular forms of cardiolipin have never been established. Here, we interrogate this simple but crucial prediction utilizing the best available model to do so, Saccharomyces cerevisiae. Specifically, we compare the ability of unremodeled and remodeled cardiolipin, which differ markedly in their acyl chain composition, to support mitochondrial activities known to require cardiolipin. Surprisingly, defined changes in the acyl chain composition of cardiolipin do not alter either mitochondrial morphology or oxidative phosphorylation. Importantly, preventing cardiolipin remodeling initiation in yeast lacking TAZ1, an ortholog of the causative gene in Barth syndrome, ameliorates mitochondrial dysfunction. Thus, our data do not support the prevailing hypothesis that unremodeled cardiolipin is functionally distinct from remodeled cardiolipin, at least for the functions examined, suggesting alternative physiological roles for this conserved pathway.     

Deacylation on the matrix side of the mitochondrial inner membrane regulates cardiolipin remodeling

The mitochondrial-specific lipid cardiolipin (CL) is required for numerous processes therein. After its synthesis on the matrix-facing leaflet of the inner membrane (IM), CL undergoes acyl chain remodeling to achieve its final form. In yeast, this process is completed by the transacylase tafazzin, which associates with intermembrane space (IMS)-facing membrane leaflets. Mutations in TAZ1 result in the X-linked cardiomyopathy Barth syndrome. Amazingly, despite this clear pathophysiological association, the physiological importance of CL remodeling is unresolved. In this paper, we show that the lipase initiating CL remodeling, Cld1p, is associated with the matrix-facing leaflet of the mitochondrial IM. Thus monolysocardiolipin generated by Cld1p must be transported to IMS-facing membrane leaflets to gain access to tafazzin, identifying a previously unknown step required for CL remodeling. Additionally, we show that Cld1p is the major site of regulation in CL remodeling; and that, like CL biosynthesis, CL remodeling is augmented in growth conditions requiring mitochondrially produced energy. However, unlike CL biosynthesis, dissipation of the mitochondrial membrane potential stimulates CL remodeling, identifying a novel feedback mechanism linking CL remodeling to oxidative phosphorylation capacity.     

Statistical analysis of the processes controlling choline and ethanolamine glycerophospholipid molecular species composition

The regulation and maintenance of the cellular lipidome through biosynthetic, remodeling, and catabolic mechanisms are critical for biological homeostasis during development, health and disease. These complex mechanisms control the architectures of lipid molecular species, which have diverse yet highly regulated fatty acid chains at both the sn1 and sn2 positions. Phosphatidylcholine (PC) and phosphatidylethanolamine (PE) serve as the predominant biophysical scaffolds in membranes, acting as reservoirs for potent lipid signals and regulating numerous enzymatic processes. Here we report the first rigorous computational dissection of the mechanisms influencing PC and PE molecular architectures from high-throughput shotgun lipidomic data. Using novel statistical approaches, we have analyzed multidimensional mass spectrometry-based shotgun lipidomic data from developmental mouse heart and mature mouse heart, lung, brain, and liver tissues. We show that in PC and PE, sn1 and sn2 positions are largely independent, though for low abundance species regulatory processes may interact with both the sn1 and sn2 chain simultaneously, leading to cooperative effects. Chains with similar biochemical properties appear to be remodeled similarly. We also see that sn2 positions are more regulated than sn1, and that PC exhibits stronger cooperative effects than PE. A key aspect of our work is a novel statistically rigorous approach to determine cooperativity based on a modified Fisher's exact test using Markov Chain Monte Carlo sampling. This computational approach provides a novel tool for developing mechanistic insight into lipidomic regulation.     

In yeast,cardiolipin unsaturation level plays a key role in mitochondrial function and inner membrane integrity

Cardiolipin is the signature phospholipid of the mitochondrial inner membrane. It participates in shaping the inner membrane as well as in modulating the activity of many membrane-bound proteins. The acyl chain composition of cardiolipin is finely tuned post-biosynthesis depending on the surrounding phospholipids to produce mature or unsaturated cardiolipin. However, experimental evidence showing that immature and mature cardiolipin are functionally equivalents for mitochondria poses doubts on the relevance of cardiolipin remodeling. In this work, we studied the role of cardiolipin acyl chain composition in mitochondrial bioenergetics, including a detailed bioenergetic profile of yeast mitochondria. Cardiolipin acyl chains were modified by genetic and nutritional manipulation. We found that both the bioenergetic efficiency and osmotic stability of mitochondria are dependent on the unsaturation level of cardiolipin acyl chains. It is proposed that cardiolipin remodeling and, consequently, mature cardiolipins play an important role in mitochondrial inner membrane integrity and functionality.     

Phospholipid turnover and acyl chain remodeling in the yeast ER

The turnover of phospholipids plays an essential role in membrane lipid homeostasis by impacting both lipid head group and acyl chain composition. This review focusses on the degradation and acyl chain remodeling of the major phospholipid classes present in the ER membrane of the reference eukaryote Saccharomyces cerevisiae, i.e. phosphatidylcholine (PC), phosphatidylinositol (PI) and phosphatidylethanolamine (PE). Phospholipid turnover reactions are introduced, and the occurrence and important functions of phospholipid remodeling in higher eukaryotes are briefly summarized. After presenting an inventory of established mechanisms of phospholipid acyl chain exchange, current knowledge of phospholipid degradation and remodeling by phospholipases and acyltransferases localized to the yeast ER is summarized. PC is subject to the PC deacylation-reacylation remodeling pathway (PC-DRP) involving a phospholipase B, the recently identified glycerophosphocholine acyltransferase Gpc1p, and the broad specificity acyltransferase Ale1p. PI is post-synthetically enriched in C18:0 acyl chains by remodeling reactions involving Cst26p. PE may undergo turnover by the phospholipid: diacylglycerol acyltransferase Lro1p as first step in acyl chain remodeling. Clues as to the functions of phospholipid acyl chain remodeling are discussed.     

The influence of the acyl chain composition of cardiolipin on the stability of mitochondrial complexes; An unexpected effect of cardiolipin in α-ketoglutarate dehydrogenase and prohibitin complexes

The role of cardiolipin acyl chain composition in assembly/stabilization of mitochondrial complexes was investigated using three yeast deletion mutants (acb1Δ strain; taz1Δ strain; and acb1Δtaz1Δ strain). Deletion of the TAZ1 gene, involved in cardiolipin acyl chain remodeling, is known to increase the content of monolyso-cardiolipin (MLCL) at the expense of CL, and to decrease the unsaturation of the remaining CL. Deletion of the ACB1 gene encoding the acyl-CoA-binding protein, involved in fatty acid elongation, decreases the average length of the CL acyl chains. Furthermore, a TAZ1ACB1 double deletion mutant strain was used in this study which has both a decrease in the length of the CL acyl chains and an increase in MLCL. BN/SDS PAGE analysis revealed that cardiolipin is important for the prohibitin–m-AAA protease complex, the α-ketoglutarate dehydrogenase complex and respiratory chain supercomplexes. The results indicate that the decreased level of complexes in taz1Δ and acb1Δtaz1Δ mitochondria is due to a decreased content of CL or the presence of MLCL.     

Shotgun lipidomics reveals the temporally dependent, highly diversified cardiolipin profile in the mammalian brain: temporally coordinated postnatal diversification of cardiolipin molecular species with neuronal remodeling

Large-scale neuronal remodeling through apoptosis occurs shortly after birth in all known mammalian species. Apoptosis, in large part, depends upon critical interactions between mitochondrial membranes and cytochrome c. Herein, we examined the hypothesis that the large-scale reorganization of neuronal circuitry after birth is accompanied by profound alterations in cardiolipin (CL) content and molecular species distribution. During embryonic development, over 100 CL molecular species were identified and quantitated in murine neuronal tissues. The embryonic CL profile was notable for the presence of abundant amounts of relatively short aliphatic chains (e.g., palmitoleic and oleic acids). In sharp contrast, after birth, the CL profile contained a remarkably complex repertoire of CL molecular species, in which the signaling fatty acids (i.e., arachidonic and docosahexaenoic acids) were markedly increased. These results identify the rapid remodeling of CL in the perinatal period with resultant alterations in the physical properties of the mitochondrial membrane. The complex distribution of aliphatic chains in the neuronal CL pool is separate and distinct from that in other organs (e.g., heart, liver, etc.), where CL molecular species contain predominantly only one major type of aliphatic chain (e.g., linoleic acid). Analyses of mRNA levels by real-time quantitative polymerase chain reactions suggested that the alterations in CL content were due to the combined effects of both attenuation of de novo CL biosynthesis and decreased remodeling of CL. Collectively, these results provide a new perspective on the complexity of CL in neuronal signaling, mitochondrial bioenergetics, and apoptosis.     

Barth syndrome cells display widespread remodeling of mitochondrial complexes without affecting metabolic flux distribution

Barth syndrome (BTHS) is a rare X-linked disorder that is characterized by cardiac and skeletal myopathy, neutropenia and growth abnormalities. The disease is caused by mutations in the tafazzin (TAZ) gene encoding an enzyme involved in the acyl chain remodeling of the mitochondrial phospholipid cardiolipin (CL). Biochemically, this leads to decreased levels of mature CL and accumulation of the intermediate monolysocardiolipin (MLCL). At a cellular level, this causes mitochondrial fragmentation and reduced stability of the respiratory chain supercomplexes. However, the exact mechanism through which tafazzin deficiency leads to disease development remains unclear. We therefore aimed to elucidate the pathways affected in BTHS cells by employing proteomic and metabolic profiling assays. Complexome profiling of patient skin fibroblasts revealed significant effects for about 200 different mitochondrial proteins. Prominently, we found a specific destabilization of higher order oxidative phosphorylation (OXPHOS) supercomplexes, as well as changes in complexes involved in cristae organization and CL trafficking. Moreover, the key metabolic complexes 2-oxoglutarate dehydrogenase (OGDH) and branched-chain ketoacid dehydrogenase (BCKD) were profoundly destabilized in BTHS patient samples. Surprisingly, metabolic flux distribution assays using stable isotope tracer-based metabolomics did not show reduced flux through the TCA cycle. Overall, insights from analyzing the impact of TAZ mutations on the mitochondrial complexome provided a better understanding of the resulting functional and structural consequences and thus the pathological mechanisms leading to Barth syndrome.     

The two biosynthetic routes leading to phosphatidylcholine in yeast produce different sets of molecular species. Evidence for lipid remodeling

Phosphatidylcholine (PC), a major lipid class in the membranes of eukaryotes, is synthesized either via the triple methylation of phosphatidylethanolamine (PE) or via the CDP-choline route. To investigate whether the two biosynthetic routes contribute differently to the steady-state profile of PC species, i.e., PC molecules with specific acyl chain compositions, the pools of newly synthesized PC species were monitored by labeling Saccharomyces cerevisiae with deuterated precursors of the two routes, (methyl-D3)-methionine and (D13)-choline, respectively. Electrospray ionization tandem mass spectrometry (ESI-MS/MS) revealed that the two PC biosynthetic pathways yield different sets of PC species, with the CDP-choline route contributing most to the molecular diversity. Moreover, yeast was shown to be capable of remodeling PC by acyl chain exchange at the sn-1 position of the glycerol backbone. Remodeling was found to be required to generate the steady-state species distribution of PC. This is the first study demonstrating a functional difference between the two biosynthetic routes in yeast.     

Murine diet-induced obesity remodels cardiac and liver mitochondrial phospholipid acyl chains with differential effects on respiratory enzyme activity

Cardiac phospholipids, notably cardiolipin, undergo acyl chain remodeling and/or loss of content in aging and cardiovascular diseases, which is postulated to mechanistically impair mitochondrial function. Less is known about how diet-induced obesity influences cardiac phospholipid acyl chain composition and thus mitochondrial responses. Here we first tested if a high fat diet remodeled murine cardiac mitochondrial phospholipid acyl chain composition and consequently disrupted membrane packing, supercomplex formation and respiratory enzyme activity. Mass spectrometry analyses revealed that mice consuming a high fat diet displayed 0.8–3.3 fold changes in cardiac acyl chain remodeling of cardiolipin, phosphatidylcholine, and phosphatidylethanolamine. Biophysical analysis of monolayers constructed from mitochondrial phospholipids of obese mice showed impairment in the packing properties of the membrane compared to lean mice. However, the high fat diet, relative to the lean controls, had no influence on cardiac mitochondrial supercomplex formation, respiratory enzyme activity, and even respiration. To determine if the effects were tissue specific, we subsequently conducted select studies with liver tissue. Compared to the control diet, the high fat diet remodeled liver mitochondrial phospholipid acyl chain composition by 0.6–5.3-fold with notable increases in n-6 and n-3 polyunsaturation. The remodeling in the liver was accompanied by diminished complex I to III respiratory enzyme activity by 3.5-fold. Finally, qRT-PCR analyses demonstrated an upregulation of liver mRNA levels of tafazzin, which contributes to cardiolipin remodeling. Altogether, these results demonstrate that diet-induced obesity remodels acyl chains in the mitochondrial phospholipidome and exerts tissue specific impairments of respiratory enzyme activity.     

Cardiolipin and Its Different Properties in Mitophagy and Apoptosis

Cardiolipin (CL) is a unique dimeric phospholipid that exists almost exclusively in the inner mitochondrial membrane (IMM) in eukaryotic cells. Two chiral carbons and four fatty acyl chains in CL result in a flexible body allowing interactions with respiratory chain complexes and mitochondrial substrate carriers. Due to its high content of unsaturated fatty acids, CL is particularly prone to reactive oxygen species (ROS)-induced oxidative attacks. Under mild mitochondrial damage, CL is redistributed to the outer mitochondrial membrane (OMM) and serves as a recognition signal for dysfunctional mitochondria, which are rapidly sequestered by autophagosomes. However, peroxidation of CL is far greater in response to severe stress than under normal or mild-damage conditions. The accumulation of oxidized CL on the OMM results in recruitment of Bax and formation of the mitochondrial permeability transition pore (MPTP), which releases Cytochrome c (Cyt c) from mitochondria. Over the past decade, the significance of CL in the function of mitochondrial bioenergy has been explored. Moreover, approaches to analyzing CL have become more effective and accurate. In this review, we discuss the unique structural features of CL as well as the current understanding of CL-based molecular mechanisms of mitophagy and apoptosis.     

Evidence for a regular distribution of cholesterol in phospholipid bilayers from diphenylhexatriene fluorescence.

Cholesterol/dimyristoylphosphatidylcholine (DMPC) multilamellar vesicles were studied by steady-state fluorescence using diphenylhexatriene (DPH) as a probe. A series of dips were found in the plot of DPH fluorescence intensity versus cholesterol concentration at certain specific cholesterol concentrations. This observation indicates that there are dominant domains in which cholesterol molecules are regularly distributed on a hexagonal superlattice in the acyl chain matrix of DMPC at critical cholesterol concentrations. These concentrations can be predicted by an equation or a mathematical series, except the one at 33 mol %. These dips of DPH fluorescence intensity are temperature dependent. The excellent agreement between experimental data and calculated values as well as similar previous findings of dips and/or kinks in the excimer-over-monomer fluorescence in pyrenephosphatidylcholine/phospholipid mixtures confirm our conclusion about lateral organizations of cholesterol and acyl lipid chains in cholesterol/phospholipid multilamellar vesicles. The regular distribution model at critical concentration is consistent with the phase diagram of cholesterol/DMPC. Using the model of regular distribution, the physical origin of the liquid-disordered (Ld) phase, liquid-ordered phase (Lo), and coexistence of liquid-disordered phase and Lo phase (Lo + Ld) is discussed on the molecular level.     

Dynamic simulation of cardiolipin remodeling: greasing the wheels for an interpretative approach to lipidomics

Cardiolipin is a class of mitochondrial specific phospholipid, which is intricately involved in mitochondrial functionality. Differences in cardiolipin species exist in a variety of tissues and diseases. It has been demonstrated that the cardiolipin profile is a key modulator of the functions of many mitochondrial proteins. However, the chemical mechanism(s) leading to normal and/or pathological distribution of cardiolipin species remain elusive. Herein, we describe a novel approach for investigating the molecular mechanism of cardiolipin remodeling through a dynamic simulation. This approach applied data from shotgun lipidomic analyses of the heart, liver, brain, and lung mitochondrial lipidomes to model cardiolipin remodeling, including relative content, regiospecificity, and isomeric composition of cardiolipin species. Generated cardiolipin profiles were nearly identical to those determined by shotgun lipidomics. Importantly, the simulated isomeric compositions of cardiolipin species were further substantiated through product ion analysis. Finally, unique enzymatic activities involved in cardiolipin remodeling were assessed from the parameters used in the dynamic simulation of cardiolipin profiles. Collectively, we described, verified, and demonstrated a novel approach by integrating both lipidomic analysis and dynamic simulation to study cardiolipin biology. We believe this study provides a foundation to investigate cardiolipin metabolism and bioenergetic homeostasis in normal and disease states.     

Metabolism of phosphatidylcholine and its implications for lipid acyl chain composition in Saccharomyces cerevisiae

Phosphatidylcholine (PC) is a very abundant membrane lipid in most eukaryotes including the model organism Saccharomyces cerevisiae. Consequently, the molecular species profile of PC, i.e. the ensemble of PC molecules with acyl chains differing in number of carbon atoms and double bonds, is important in determining the physical properties of eukaryotic membranes, and should be tightly regulated. In this review current insights in the contributions of biosynthesis, turnover, and remodeling by acyl chain exchange to the maintenance of PC homeostasis at the level of the molecular species in yeast are summarized. In addition, the phospholipid class-specific changes in membrane acyl chain composition induced by PC depletion are discussed, which identify PC as key player in a novel regulatory mechanism balancing the proportions of bilayer and non-bilayer lipids in yeast.     

Distribution of Pentachlorophenol in Phospholipid Bilayers: A Molecular Dynamics Study

Molecular dynamics computer simulations of pentachlorophenol (PCP) in palmitoyl-oleoyl-phosphatidylethanolamine and palmitoyl-oleoyl-phosphatidylcholine lipid bilayers were carried out to investigate the distribution of PCP and the effects of PCP on the phospholipid bilayer structure. Starting from two extreme starting structures, including PCP molecules outside the lipid bilayer, the PCP distribution converges in simulations of up to 50 ns. PCP preferentially occupies the region between the carbonyl groups and the double bonds in the acyl chains of the lipid molecules in the bilayer. In the presence of PCP, the lipid chain order increases somewhat in both chains, and the average tilt angle of the lipid chains decreases. The increase in the lipid chain order in the presence of PCP was more pronounced in the palmitoyl-oleoyl-phosphatidylcholine bilayer compared to the palmitoyl-oleoyl-phosphatidylethanolamine bilayer. The number of trans conformations of lipid chain dihedrals does not change significantly. PCP aligns parallel to the alkyl chains of the lipid to optimize the packing in the dense ordered chain region of the bilayer. The hydroxyl group of PCP forms hydrogen bonds with both water and lipid oxygen atoms in the water/lipid interface region.     

Cellular functions of cardiolipin in yeast

Cardiolipin (CL), the signature lipid of mitochondria, plays a critical role in mitochondrial function and biogenesis. The availability of yeast mutants blocked in CL synthesis has facilitated studies of the biological role of this lipid. Perturbation of CL synthesis leads to growth defects not only during respiratory growth but also under conditions in which respiration is not essential. CL was shown to play a role in mitochondrial protein import, cell wall biogenesis, aging and apoptosis, ceramide synthesis, and translation of electron transport chain components. The genetic disorder Barth syndrome (BTHS) is caused by mutations in the tafazzin gene resulting in decreased total CL levels, accumulation of monolysocardiolipin (MLCL), and decreased unsaturated fatty acyl species of CL. The variation in clinical presentation of BTHS indicates that other physiological factors play a significant role in modifying the phenotype resulting from tafazzin deficiency. Elucidating the functions of CL is expected to shed light on the role of this important lipid in BTHS and other disorders of mitochondrial dysfunction.     

Defective remodeling of cardiolipin and phosphatidylglycerol in Barth syndrome

Cardiolipin (CL) and phosphatidylglycerol (PG) are the major polyglycerophospholipids observed in mammalian tissues. CL is exclusively found in the inner mitochondrial membrane and is required for optimal function of many of the respiratory and ATP-synthesizing enzymes. The role of CL in oxidative phosphorylation is, however, not fully understood and although reduced CL content leads to aberrant cell function, no human disorders with a primary defect in cardiolipin metabolism have been described. In this paper we present evidence that patients with the rare disorder X-linked cardioskeletal myopathy and neutropenia (Barth syndrome, MIM 302060) have a primary defect in CL and PG remodeling. We investigated phospholipid metabolism in cultured skin fibroblasts of patients and show that the biosynthesis rate of PG and CL is normal but that the CL pool size is 75% reduced, indicating accelerated degradation. Moreover, the incorporation of linoleic acid, which is the characteristic acyl side chain found in mammalian CL, into both PG and CL is significantly reduced, whereas the incorporation of other fatty acids into these phospholipids is normal. We show that this defect was only observed in Barth syndrome patients' cells and not in cells obtained from patients with primary defects in the respiratory chain, demonstrating that the observed defect is not secondary to respiratory chain dysfunction. These results imply that the G4.5 gene product, which is mutated in Barth syndrome patients, is specifically involved in the remodeling of PG and CL and for the first time identify an essential factor in this important cellular process.