Reduced creatine-stimulated respiration in doxorubicin challenged mitochondria: particular sensitivity of the heart

Doxorubicin (DXR) belongs to the most efficient anticancer drugs. However, its use is limited by a risk of cardiotoxicity, which is not completely understood. Recently, we have shown that DXR impairs essential properties of purified mitochondrial creatine kinase (MtCK), with cardiac isoenzyme (sMtCK) being particularly sensitive. In this study we assessed the effects of DXR on respiration of isolated structurally and functionally intact heart mitochondria, containing sMtCK, in the presence and absence of externally added creatine (Cr), and compared these effects with the response of brain mitochondria expressing uMtCK, the ubiquitous, non-muscle MtCK isoenzyme. DXR impaired respiration of isolated heart mitochondria already after short-term exposure (minutes), affecting both ADP- and Cr-stimulated respiration. During a first short time span (minutes to 1 h), detachment of MtCK from membranes occurred, while a decrease of MtCK activity related to oxidative damage was only observed after longer exposure (several hours). The early inhibition of Cr-stimulated respiration, in addition to impairment of components of the respiratory chain involves a partial disturbance of functional coupling between MtCK and ANT, likely due to interaction of DXR with cardiolipin leading to competitive inhibition of MtCK/membrane binding. The relevance of these findings for the regulation of mitochondrial energy production in the heart, as well as the obvious differences of DXR action in the heart as compared to brain tissue, is discussed.     

Reduced creatine-stimulated respiration in doxorubicin challenged mitochondria: Particular sensitivity of the heart

Doxorubicin (DXR) belongs to the most efficient anticancer drugs. However, its use is limited by a risk of cardiotoxicity, which is not completely understood. Recently, we have shown that DXR impairs essential properties of purified mitochondrial creatine kinase (MtCK), with cardiac isoenzyme (sMtCK) being particularly sensitive. In this study we assessed the effects of DXR on respiration of isolated structurally and functionally intact heart mitochondria, containing sMtCK, in the presence and absence of externally added creatine (Cr), and compared these effects with the response of brain mitochondria expressing uMtCK, the ubiquitous, non-muscle MtCK isoenzyme. DXR impaired respiration of isolated heart mitochondria already after short-term exposure (minutes), affecting both ADP- and Cr-stimulated respiration. During a first short time span (minutes to 1 h), detachment of MtCK from membranes occurred, while a decrease of MtCK activity related to oxidative damage was only observed after longer exposure (several hours). The early inhibition of Cr-stimulated respiration, in addition to impairment of components of the respiratory chain involves a partial disturbance of functional coupling between MtCK and ANT, likely due to interaction of DXR with cardiolipin leading to competitive inhibition of MtCK/membrane binding. The relevance of these findings for the regulation of mitochondrial energy production in the heart, as well as the obvious differences of DXR action in the heart as compared to brain tissue, is discussed.     

C-terminal lysines determine phospholipid interaction of sarcomeric mitochondrial creatine kinase

High affinity interaction between octameric mitochondrial creatine kinase (MtCK) and the phospholipid cardiolipin in the inner mitochondrial membrane plays an important role in metabolite channeling between MtCK and inner membrane adenylate translocator, which itself is tightly bound to cardiolipin. Three C-terminal basic residues revealed as putative cardiolipin anchors in the x-ray structures of MtCK and corresponding to lysines in human sarcomeric MtCK (sMtCK) were exchanged by in vitro mutagenesis (K369A/E, K379Q/A/E, K380Q/A/E) to yield double and triple mutants. sMtCK proteins were bacterially expressed, purified to homogeneity, and verified for structural integrity by enzymatic activity, gel filtration chromatography, and CD spectroscopy. Interaction with cardiolipin and other acidic phospholipids was quantitatively analyzed by light scattering, surface plasmon resonance, and fluorescence spectroscopy. All mutant sMtCKs showed a strong decrease in vesicle cross-linking, membrane affinity, binding capacity, membrane ordering capability, and binding-induced changes in protein structure as compared with wild type. These effects did not depend on the nature of the replacing amino acid but on the number of exchanged lysines. They were moderate for Lys-379/Lys-380 double mutants but pronounced for triple mutants, with a 30-fold lower membrane affinity and an entire lack of alterations in protein structure compared with wild-type sMtCK. However, even triple mutants partially maintained an increased order of cardiolipin-containing membranes. Thus, the three C-terminal lysines determine high affinity sMtCK/cardiolipin interaction and its effects on MtCK structure, whereas low level binding and some effect on membrane fluidity depend on other structural components. These results are discussed in regard to MtCK microcompartments and evolution.     

Mitochondrial kinases and their molecular interaction with cardiolipin

Mitochondrial isoforms of creatine kinase (MtCK) and nucleoside diphosphate kinase (NDPK-D) are not phylogenetically related but share functionally important properties. They both use mitochondrially generated ATP with the ultimate goal of maintaining proper nucleotide pools, are located in the intermembrane/cristae space, have symmetrical oligomeric structures, and show high affinity binding to anionic phospholipids, in particular cardiolipin. The structural basis and functional consequences of the cardiolipin interaction have been studied and are discussed in detail in this review. They mainly result in a functional interaction of MtCK and NDPK-D with inner membrane adenylate translocator, probably by forming proteolipid complexes. These interactions allow for privileged exchange of metabolites (channeling) that ultimately regulate mitochondrial respiration. Further functions of the MtCK/membrane interaction include formation of cardiolipin membrane patches, stabilization of mitochondria and a role in apoptotic signaling, as well as in case of both kinases, a role in facilitating lipid transfer between two membranes. Finally, disturbed cardiolipin interactions of MtCK, NDPK-D and other proteins like cytochrome c and truncated Bid are discussed more generally in the context of apoptosis and necrosis.     

Cardiolipin dynamics and binding to conserved residues in the mitochondrial ADP/ATP carrier

Cardiolipin in eukaryotes is found in the mitochondrial inner membrane, where it interacts with membrane proteins and, although not essential, is necessary for the optimal activity of a number of proteins. One of them is the mitochondrial ADP/ATP carrier, which imports ADP into the mitochondrion and exports ATP. In the crystal structures, cardiolipin is bound to three equivalent sites of the ADP/ATP carrier, but its role is unresolved. Conservation of residues at these cardiolipin binding sites across other members of the mitochondrial carrier superfamily indicates cardiolipin binding is likely to be important for the function of all mitochondrial carriers. Multiscale simulations were performed in a cardiolipin-containing membrane to investigate the dynamics of cardiolipin around the yeast and bovine ADP/ATP carriers in a lipid bilayer and the properties of the cardiolipin-binding sites. In coarse-grain simulations, cardiolipin molecules bound to the carriers for longer periods of time than phosphatidylcholine and phosphatidylethanolamine lipids—with timescales in the tens of microseconds. Three long-lived cardiolipin binding sites overlapped with those in the crystal structures of the carriers. Other shorter-lived cardiolipin interaction sites were identified in both membrane leaflets. However, the timescales of the interactions were of the same order as phosphatidylcholine and phosphatidylethanolamine, suggesting that these sites are not specific for cardiolipin binding. The calculation of lipid binding times and the overlap of the cardiolipin binding sites between the structures and simulations demonstrate the potential of multiscale simulations to investigate the dynamics and behavior of lipids interacting with membrane proteins.     

Mitochondrial creatine kinase in contact sites: interaction with porin and adenine nucleotide translocase, role in permeability transition and sensitivity to oxidative damage

The creatine/phosphocreatine circuit provides an efficient energy buffering and transport system in a variety of cells with high and fluctuating energy requirements. It connects sites of energy production (mitochondria, glycolysis) with sites of energy consumption (various cellular ATPases). The cellular creatine/phosphocreatine pool is linked to the ATP/ADP pool by the action of different isoforms of creatine kinase located at distinct subcellular compartments. Octameric mitochondrial creatine kinase (MtCK), together with porin and adenine nucleotide translocase, forms a microcompartment at contact sites between inner and outer mitochondrial membranes and facilitates the production and export into the cytosol of phosphocreatine. MtCK is probably in direct protein-protein contact with outer membrane porin, whereas interaction with inner membrane adenine nucleotide translocase is rather mediated by acidic phopholipids (like cardiolipin) present in significant amounts in the inner membrane. Octamer-dimer transitions of MtCK as well as different creatine kinase substrates have a profound influence on controlling mitochondrial permeability transition (MPT). Inactivation by reactive oxygen species of MtCK and destabilization of its octameric structure are factors that contribute to impairment of energy homeostasis and facilitated opening of the MPT pore, which eventually lead to tissue damage during periods of ischemia/reperfusion.     

Novel lipid transfer property of two mitochondrial proteins that bridge the inner and outer membranes

This study provides evidence of a novel function for mitochondrial creatine kinase (MtCK) and nucleoside diphosphate kinase (NDPK-D). Both are basic peripheral membrane proteins with symmetrical homo-oligomeric structure, which in the case of MtCK was already shown to allow crossbridging of lipid bilayers. Here, different lipid dilution assays clearly demonstrate that both kinases also facilitate lipid transfer from one bilayer to another. Lipid transfer occurs between liposomes mimicking the lipid composition of mitochondrial contact sites, containing 30 mol % cardiolipin, but transfer does not occur when cardiolipin is replaced by phosphatidylglycerol. Ubiquitous MtCK, but not NDPK-D, shows some specificity in the nature of the lipids transferred and it is not active with phosphatidylcholine alone. MtCK can undergo reversible oligomerization between dimeric and octameric forms, but only the octamer can bridge membranes and promote lipid transfer. Cytochrome c, another basic mitochondrial protein known to bind to anionic membranes but not crosslinking them, is also incapable of promoting lipid transfer. The lipid transfer process does not involve vesicle fusion or loss of the internal contents of the liposomes.     

Lipid Fingerprints and Cofactor Dynamics of Light-Harvesting Complex II in Different Membranes

Plant light-harvesting complex II (LHCII) is the key antenna complex for plant photosynthesis. We present coarse-grained molecular dynamics simulations of monomeric and trimeric LHCII in a realistic thylakoid membrane environment based on the Martini force field. The coarse-grained protein model has been optimized with respect to atomistic reference simulations. Our simulations provide detailed insights in the thylakoid lipid fingerprint of LHCII which compares well with experimental data from membrane protein purification. Comparing the monomer and trimeric LHCII reveals a stabilizing effect of trimerization on the chromophores as well as the protein. Moreover, the average chromophore distance shortens in the trimer leading to stronger excitonic couplings. When changing the native thylakoid environment to a model membrane the protein flexibility remains constant, whereas the chromophore flexibility is reduced. Overall, the presented LHCII model lays the foundation to investigate the μs dynamics of this key antenna protein of plants.     

Metabolic control analysis of integrated energy metabolism in permeabilized cardiomyocytes - experimental study

The main focus of this research was to apply Metabolic Control Analysis to quantitative investigation of the regulation of respiration by components of the Mitochondrial Interactosome (MI, a supercomplex consisting of ATP Synthasome, mitochondrial creatine kinase (MtCK), voltage dependent anion channel (VDAC), and tubulin) in permeabilized cardiomyocytes. Flux control coefficients (FCC) were measured using two protocols: 1) with direct ADP activation, and 2) with MtCK activation by creatine (Cr) in the presence of ATP and pyruvate kinase-phosphoenolpyruvate system. The results show that the metabolic control is much stronger in the latter case: the sum of the measured FCC is 2.7 versus 0.74 (ADP activation). This is consistent with previous data showing recycling of ADP and ATP inside the MI due to the functional coupling between MtCK and ANT and limited permeability of VDAC for these compounds, PCr being the major energy carrier between the mitochondria and ATPases. In physiological conditions, when the MI is activated, the key sites of regulation of respiration in mitochondria are MtCK (FCC = 0.93), adenine nucleotide translocase ANT (FCC = 0.95) and CoQ cytochrome c oxidoreductase (FCC = 0.4). These results show clearly that under the physiological conditions the energy transfer from mitochondria to the cytoplasm is regulated by the MI supercomplex and is very sensitive to metabolic signals.     

High efficiency of energy flux controls within mitochondrial interactosome in cardiac intracellular energetic units

The aim of our study was to analyze a distribution of metabolic flux controls of all mitochondrial complexes of ATP-Synthasome and mitochondrial creatine kinase (MtCK) in situ in permeabilized cardiac cells. For this we used their specific inhibitors to measure flux control coefficients (C(vi)(JATP)) in two different systems: A) direct stimulation of respiration by ADP and B) activation of respiration by coupled MtCK reaction in the presence of MgATP and creatine. In isolated mitochondria the C(vi)(JATP) were for system A: Complex I - 0.19, Complex III - 0.06, Complex IV 0.18, adenine nucleotide translocase (ANT) - 0.11, ATP synthase - 0.01, Pi carrier - 0.20, and the sum of C(vi)(JATP) was 0.75. In the presence of 10mM creatine (system B) the C(vi)(JATP) were 0.38 for ANT and 0.80 for MtCK. In the permeabilized cardiomyocytes inhibitors had to be added in much higher final concentration, and the following values of C(vi)(JATP) were determined for condition A and B, respectively: Complex I - 0.20 and 0.64, Complex III - 0.41 and 0.40, Complex IV - 0.40 and 0.49, ANT - 0.20 and 0.92, ATP synthase - 0.065 and 0.38, Pi carrier - 0.06 and 0.06, MtCK 0.95. The sum of C(vi)(JATP) was 1.33 and 3.84, respectively. Thus, C(vi)(JATP) were specifically increased under conditions B only for steps involved in ADP turnover and for Complex I in permeabilized cardiomyocytes within Mitochondrial Interactosome, a supercomplex consisting of MtCK, ATP-Synthasome, voltage dependent anion channel associated with tubulin βII which restricts permeability of the mitochondrial outer membrane.     

Octamers of mitochondrial creatine kinase isoenzymes differ in stability and membrane binding

Octamer stability and membrane binding of mitochondrial creatine kinase (MtCK) are important for proper functioning of the enzyme and were suggested as targets for regulatory mechanisms. A quantitative analysis of these properties, using fluorescence spectroscopy, gel filtration, and surface plasmon resonance, revealed substantial differences between the two types of MtCK isoenzymes, sarcomeric (sMtCK) and ubiquitous (uMtCK). As compared with human and chicken sMtCK, human uMtCK showed a 23-34 times slower octamer dissociation rate, a reduced reoctamerization rate and a superior octamer stability as deduced from the octamer/dimer ratios at thermodynamic equilibrium. Octamer stability of sMtCK increased with temperature up to 30 degrees C, indicating a substantial contribution of hydrophobic interactions, while it decreased in the case of uMtCK, indicating the presence of additional polar dimer/dimer interactions. These conclusions are consistent with the recently solved x-ray structure of the human uMtCK (Eder, M., Fritz-Wolf, K., Kabsch, W., Wallimann, T., and Schlattner, U. (2000) Proteins 39, 216-225). When binding to 16% cardiolipin membranes, sMtCK showed slightly faster on-rates and higher affinities than uMtCK. However, human uMtCK was able to recruit the highest number of binding sites on the vesicle surface. The observed divergence of ubiquitous and sarcomeric MtCK is discussed with respect to their molecular structures and the possible physiological implications.     

Multiscale methods for macromolecular simulations

In this article we review the key modeling tools available for simulating biomolecular systems. We consider recent developments and representative applications of mixed quantum mechanics/molecular mechanics (QM/MM), elastic network models (ENMs), coarse-grained molecular dynamics, and grid-based tools for calculating interactions between essentially rigid protein assemblies. We consider how the different length scales can be coupled, both in a sequential fashion (e.g. a coarse-grained or grid model using parameterization from MD simulations), and via concurrent approaches, where the calculations are performed together and together control the progression of the simulation. We suggest how the concurrent coupling approach familiar in the context of QM/MM calculations can be generalized, and describe how this has been done in the CHARMM macromolecular simulation package.     

Coarse-grained simulations of protein-protein association: an energy landscape perspective

Understanding protein-protein association is crucial in revealing the molecular basis of many biological processes. Here, we describe a theoretical simulation pipeline to study protein-protein association from an energy landscape perspective. First, a coarse-grained model is implemented and its applications are demonstrated via molecular dynamics simulations for several protein complexes. Second, an enhanced search method is used to efficiently sample a broad range of protein conformations. Third, multiple conformations are identified and clustered from simulation data and further projected on a three-dimensional globe specifying protein orientations and interacting energies. Results from several complexes indicate that the crystal-like conformation is favorable on the energy landscape even if the landscape is relatively rugged with metastable conformations. A closer examination on molecular forces shows that the formation of associated protein complexes can be primarily electrostatics-driven, hydrophobics-driven, or a combination of both in stabilizing specific binding interfaces. Taken together, these results suggest that the coarse-grained simulations and analyses provide an alternative toolset to study protein-protein association occurring in functional biomolecular complexes.     

Molecular interactions between mitochondrial membrane proteins and the C-terminal domain of PB1-F2: an in silico approach

PB1-F2 is a recently described influenza A viral protein that induces apoptosis by binding with two mitochondrial membrane proteins, i.e. VDAC1 (outer membrane) and ANT3 (inner membrane). Knowledge of this binding mechanism could provide insights that would aid in the design of novel inhibitors against this protein. Therefore, to better understand these interactions, we have undertaken this study to model the PB1-F2 protein of the highly pathogenic influenza A virus subtype H5N1. Moreover, a model of human ANT3 was also established. The dynamics of the molecular interactions between the C-terminal region of PB1-F2 protein and VDAC1 and ANT3 were expounded by employing an in silico approach. Our results suggest the involvement of 12 amino acids of PB1-F2 protein, which form hydrophobic contacts with 22 amino acids of VDAC1. Of these, Leu64, Arg75 and Val76 were found to be crucial for mitochondrial targetting. In the case of the PB1-F2-ANT3 complex, 14 amino acids of ANT3 were found to make hydrophobic contacts with 9 amino acids of PB1-F2. Furthermore, two hydrogen bonds were predicted in both complexes PB1-F2/VDAC1 and PB1-F2/ANT3. This study reveals the molecular interactions required for PB1-F2-induced apoptosis and suggests a hypothetical model for future study.     

Mitochondrial creatine kinase and mitochondrial outer membrane porin show a direct interaction that is modulated by calcium

Mitochondrial creatine kinase (MtCK) co-localizes with mitochondrial porin (voltage-dependent anion channel) and adenine nucleotide translocator in mitochondrial contact sites. A specific, direct protein-protein interaction between MtCK and mitochondrial porin was demonstrated using surface plasmon resonance spectroscopy. This interaction was independent of the immobilized binding partner (porin reconstituted in liposomes or MtCK) or the analyzed isoform (chicken sarcomeric MtCK or human ubiquitous MtCK, human recombinant porin, or purified bovine porin). Increased ionic strength reduced the binding of MtCK to porin, suggesting predominantly ionic interactions. By contrast, micromolar concentrations of Ca(2+) increased the amount of bound MtCK, indicating a physiological regulation of complex formation. No interaction of MtCK with reconstituted adenine nucleotide translocator was detectable in our experimental setup. The relevance of these findings for structure and function of mitochondrial contact sites is discussed.     

Multiscale modeling of lipids and lipid bilayers

A multiscale modeling approach is applied for simulations of lipids and lipid assemblies on mesoscale. First, molecular dynamics simulation of initially disordered system of lipid molecules in water within all-atomic model was carried out. On the next stage, structural data obtained from the molecular dynamics (MD) simulation were used to build a coarse-grained (ten sites) lipid model, with effective interaction potentials computed by the inverse Monte Carlo method. Finally, several simulations of the coarse-grained model on longer length- and time-scale were performed, both within Monte Carlo and molecular dynamics simulations: a periodical sample of lipid molecules ordered in bilayer, a free sheet of such bilayer without periodic boundary conditions, formation of vesicle from a plain membrane, process of self-assembly of lipids randomly dispersed in volume. It was shown that the coarse-grained model, developed exclusively from all-atomic simulation data, reproduces well all the basic features of lipids in water solution.     

PIP2-induced membrane binding of the vinculin tail competes with its other binding partners

Vinculin plays a key role during the first phase of focal adhesion formation and interacts with the plasma membrane through specific binding of its tail domain to the lipid phosphatidylinositol 4,5-bisphosphate (PIP₂). Our understanding of the PIP₂-vinculin interaction has been hampered by contradictory biochemical and structural data. Here, we used a multiscale molecular dynamics simulation approach, in which unbiased coarse-grained molecular dynamics were used to generate starting structures for subsequent microsecond-long all-atom simulations. This allowed us to map the interaction of the vinculin tail with PIP₂-enriched membranes in atomistic detail. In agreement with experimental data, we have shown that membrane binding is sterically incompatible with the intramolecular interaction between vinculin’s head and tail domain. Our simulations further confirmed biochemical and structural results, which identified two positively charged surfaces, the basic collar and the basic ladder, as the main PIP₂ interaction sites. By introducing a valency-disaggregated binding network analysis, we were able to map the protein-lipid interactions in unprecedented detail. In contrast to the basic collar, in which PIP₂ is specifically recognized by an up to hexavalent binding pocket, the basic ladder forms a series of low-valency binding sites. Importantly, many of these PIP₂ binding residues are also involved in maintaining vinculin in a closed, autoinhibited conformation. These findings led us to propose a molecular mechanism for the coupling between vinculin activation and membrane binding. Finally, our refined binding site suggests an allosteric relationship between PIP₂ and F-actin binding that disfavors simultaneous interaction with both ligands, despite nonoverlapping binding sites.