Function-Related Asymmetry of the Interactions between Matrix Loops and Conserved Sequence Motifs in the Mitochondrial ADP/ATP Carrier.

The ADP/ATP carrier (AAC) plays a central role in oxidative metabolism by exchanging ATP and ADP across the inner mitochondrial membrane. Previous experiments have shown the involvement of the matrix loops of AAC in its function, yet potential mechanisms remain largely elusive. One obstacle is the limited information on the structural dynamics of the matrix loops. In the current work, unbiased all-atom molecular dynamics (MD) simulations were carried out on c-state wild-type AAC and mutants. Our results reveal that: (1) two ends of a matrix loop are tethered through interactions between the residue of triplet 38 (Q38, D143 and Q240) located at the C-end of the odd-numbered helix and residues of the [YF]xG motif located before the N-end of the short matrix helix in the same domain; (2) the initial progression direction of a matrix loop is determined by interactions between the negatively charged residue of the [DE]G motif located at the C-end of the short matrix helix and the capping arginine (R30, R139 and R236) in the previous domain; (3) the two chemically similar residues D and E in the highly conserved [DE]G motif are actually quite different; (4) the N-end of the M3 loop is clamped by the [DE]G motif and the capping arginine of domain 2 from the two sides, which strengthens interactions between domain 2 and domain 3; and (5) a highly asymmetric stable core exists within domains 2 and 3 at the m-gate level. Moreover, our results help explain almost all extremely conserved residues within the matrix loops of the ADP/ATP carriers from a structural point of view. Taken together, the current work highlights asymmetry in the three matrix loops and implies a close relationship between asymmetry and ADP/ATP transport.

Mechanistic insights into multiple-step transport of mitochondrial ADP/ATP carrier.

The ADP/ATP carrier (AAC) is crucial for mitochondrial functions by importing ADP and exporting ATP across the inner mitochondrial membrane. However, the mechanism of highly specific ADP recognition and transport by AAC remains largely elusive. In this work, spontaneous ADP binding process to the ground c-state AAC was investigated through rigorous molecular dynamics simulations of over 31 microseconds in total. With improved simulation strategy, we have successfully identified a highly specific ADP binding site in the upper region of the cavity, and this site exhibits selectivity for ADP over ATP based on free-energy calculations. Sequence analyses on adenine nucleotide transporters also suggest that this subgroup uses the upper region of the cavity, rather than the previously proposed central binding site located at the bottom of the cavity to discriminate their substrates. Identification of the new site unveils the unusually high substrate specificity of AAC and explains the dependence of transport on the flexibility between anti and syn glycosidic conformers of ADP. Moreover, this new site together with the central site supports early biochemical findings. In light of these early findings, our simulations described a multi-step model in which the carrier uses different sites for substrate attraction, recognition and conformational transition. These results provide new insights into the transport mechanism of AAC and other adenine nucleotide transporters.

The effects of cardiolipin on the structural dynamics of the mitochondrial ADP/ATP carrier in its cytosol-open state.

Cardiolipin (CL) has been shown to play a crucial role in regulating the function of proteins in the inner mitochondrial membrane. As the most abundant protein of the inner mitochondrial membrane, the ADP/ATP carrier (AAC) has long been the model of choice to study CL-protein interactions, and specifically bound CLs have been identified in a variety of crystal structures of AAC. However, how CL binding affects the structural dynamics of AAC in atomic detail remains largely elusive. Here we compared all-atom molecular dynamics simulations on bovine AAC1 in lipid bilayers with and without CLs. Our results show that on the current microsecond simulation time scale: 1) CL binding does not significantly affect overall stability of the carrier or structural symmetry at the matrix-gate level; 2) pocket volumes of the carrier and interactions involved in the matrix-gate network become more heterogeneous in parallel simulations with membranes containing CLs; 3) CL binding consistently strengthens backbone hydrogen bonds within helix H2 near the matrix side; and 4) CLs play a consistent stabilizing role on the domain 1-2 interface through binding with the R30:R71:R151 stacking structure and fixing the M2 loop in a defined conformation. CL is necessary for the formation of this stacking structure, and this structure in turn forms a very stable CL binding site. Such a delicate equilibrium suggests the strictly conserved R30:R71:R151stacking structure of AACs could function as a switch under regulation of CLs. Taken together, these results shed new light on the CL-mediated modulation of AAC function.

Investigating the Broad Matrix-Gate Network in the Mitochondrial ADP/ATP Carrier through Molecular Dynamics Simulations.

The mitochondrial ADP/ATP carrier (AAC) exports ATP and imports ADP through alternating between cytosol-open (c-) and matrix-open (m-) states. The salt bridge networks near the matrix side (m-gate) and cytosol side (c-gate) are thought to be crucial for state transitions, yet our knowledge on these networks is still limited. In the current work, we focus on more conserved m-gate network in the c-state AAC. All-atom molecular dynamics (MD) simulations on a variety of mutants and the CATR-AAC complex have revealed that: (1) without involvement of other positive residues, the charged residues from the three Px[DE]xx[KR] motifs only are prone to form symmetrical inter-helical network; (2) R235 plays a determinant role for the asymmetry in m-gate network of AAC; (3) R235 significantly strengthens the interactions between H3 and H5; (4) R79 exhibits more significant impact on m-gate than R279; (5) CATR promotes symmetry in m-gate mainly through separating R234 from D231 and fixing R79; (6) vulnerability of the H2-H3 interface near matrix side could be functionally important. Our results provide new insights into the highly conserved yet variable m-gate network in the big mitochondrial carrier family.

Leber's Hereditary Optic Neuropathy Arising From the Synergy Between ND1 3635G>A Mutation and Mitochondrial YARS2 Mutations.

Purpose: To investigate the mechanism underlying the synergic interaction between Leber's hereditary optic neuropathy (LHON)-associated ND1 and mitochondrial tyrosyl-tRNA synthetase (YARS2) mutations. Methods: Molecular dynamics simulation and differential scanning fluorimetry were used to evaluate the structure and stability of proteins. The impact of ND1 3635G>A and YARS2 p.G191V mutations on the oxidative phosphorylation machinery was evaluated using blue native gel electrophoresis and enzymatic activities assays. Assessment of reactive oxygen species (ROS) production in cell lines was performed by flow cytometry with MitoSOX Red reagent. Analysis of effect of mutations on autophagy was undertaken via flow cytometry for autophagic flux. Results: Members of one Chinese family bearing both the YARS2 p.191Gly>Val and m.3635G>A mutations exhibited much higher penetrance of optic neuropathy than those pedigrees carrying only the m.3635G>A mutation. The m.3635G>A (p.Ser110Asn) mutation altered the ND1 structure and function, whereas the p.191Gly>Val mutation affected the stability of YARS2. Lymphoblastoid cell lines harboring both m.3635G>A and p.191Gly>Val mutations revealed more reductions in the levels of mitochondrion-encoding ND1 and CO2 than cells bearing only the m.3635G>A mutation. Strikingly, both m.3635G>A and p.191Gly>Val mutations exhibited decreases in the nucleus-encoding subunits of complex I and IV. These deficiencies manifested greater defects in the stability and activities of complex I and complex IV and overproduction of ROS and promoted greater autophagy in cell lines harboring both m.3635G>A and p.191Gly>Val mutations compared with cells bearing only the m.3635G>A mutation. Conclusions: Our findings provide new insights into the pathophysiology of LHON arising from the synergy between ND1 3635G>A mutation and mitochondrial YARS2 mutations.

An animal model for mitochondrial tyrosyl-tRNA synthetase deficiency reveals links between oxidative phosphorylation and retinal function.

Mitochondria maintain a distinct pool of ribosomal machinery, including tRNAs and tRNAs activating enzymes, such as mitochondrial tyrosyl-tRNA synthetase (YARS2). Mutations in YARS2, which typically lead to the impairment of mitochondrial protein synthesis, have been linked to an array of human diseases including optic neuropathy. However, the lack of YARS2 mutation animal model makes us difficult to elucidate the pathophysiology underlying YARS2 deficiency. To explore this system, we generated YARS2 knockout (KO) HeLa cells and zebrafish using CRISPR/Cas9 technology. We observed the aberrant tRNATyr aminoacylation overall and reductions in the levels in mitochondrion- and nucleus-encoding subunits of oxidative phosphorylation system (OXPHOS), which were especially pronounced effects in the subunits of complex I and complex IV. These deficiencies manifested the decreased levels of intact supercomplexes overall. Immunoprecipitation assays showed that YARS2 bound to specific subunits of complex I and complex IV, suggesting the posttranslational stabilization of OXPHOS. Furthermore, YARS2 ablation caused defects in the stability and activities of OXPHOS complexes. These biochemical defects could be rescued by the overexpression of YARS2 cDNA in the YARS2KO cells. In zebrafish, the yars2KO larva conferred deficient COX activities in the retina, abnormal mitochondrial morphology, and numbers in the photoreceptor and retinal ganglion cells. The zebrafish further exhibited the retinal defects affecting both rods and cones. Vision defects in yars2KO zebrafish recapitulated the clinical phenotypes in the optic neuropathy patients carrying the YARS2 mutations. Our findings highlighted the critical role of YARS2 in the stability and activity of OXPHOS and its pathological consequence in vision impairments.

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.

Novel ACADVL variants resulting in mitochondrial defects in long-chain acyl-CoA dehydrogenase deficiency.

The pathogenesis of very-long-chain acyl-CoA dehydrogenase (VLCAD) deficiency is highly heterogeneous and still unclear. Additional novel variants have been recently detected in the population. The molecular and cellular effects of these previously unreported variants are still poorly understood and require further characterization. To address this problem, we have evaluated the various functions and biochemical consequences of six novel missense variants that lead to mild VLCAD deficiency. Marked deficiencies in fatty acid oxidation (FAO) and other mitochondrial defects were observed in cells carrying one of these six variants (c.541C>T, c.863T>G, c.895A>G, c.1238T>C, c.1276G>A, and c.1505T>A), including reductions in mitochondrial respiratory-chain function and adenosine triphosphate (ATP) production, and increased levels of mitochondrial reactive oxygen species (ROS). Intriguingly, higher apoptosis levels were found in cells carrying the mutant VLCAD under glucose-limited stress. Moreover, the stability of the mutant homodimer was disturbed, and major conformational changes in each mutant VLCAD structure were predicted by molecular dynamics (MD) simulation. The data presented here may provide valuable information for improving management of diagnosis and treatment of VLCAD deficiency and for a better understanding of the general molecular bases of disease variability.

Complex I mutations synergize to worsen the phenotypic expression of Leber's hereditary optic neuropathy.

Leber's hereditary optic neuropathy (LHON) is a maternal inheritance of eye disease because of the mitochondrial DNA (mtDNA) mutations. We previously discovered a 3866T>C mutation within the gene for the ND1 subunit of complex I as possibly amplifying disease progression for patients bearing the disease-causing 11778G>A mutation within the gene for the ND4 subunit of complex I. However, whether and how the ND1 mutation exacerbates the ND4 mutation were unknown. In this report, we showed that four Chinese families bearing both m.3866T>C and m.11778G>A mutations exhibited higher penetrances of LHON than 6 Chinese pedigrees carrying only the m.3866T>C mutation or families harboring only the m.11778G>A mutation. The protein structure analysis revealed that the m.3866T>C (I187T) and m.11778G>A (R340H) mutations destabilized the specific interactions with other residues of ND1 and ND4, thereby altering the structure and function of complex I. Cellular data obtained using cybrids, constructed by transferring mitochondria from the Chinese families into mtDNA-less (ρ°) cells, demonstrated that the mutations perturbed the stability, assembly, and activity of complex I, leading to changes in mitochondrial ATP levels and membrane potential and increasing the production of reactive oxygen species. These mitochondrial dysfunctions promoted the apoptotic sensitivity of cells and decreased mitophagy. Cybrids bearing only the m.3866T>C mutation displayed mild mitochondrial dysfunctions, whereas those harboring both m.3866T>C and m.11778G>A mutations exhibited greater mitochondrial dysfunctions. These suggested that the m.3866T>C mutation acted in synergy with the m.11778G>A mutation, aggravating mitochondrial dysfunctions and contributing to higher penetrance of LHON in these families carrying both mtDNA mutations.

Overexpression of mitochondrial histidyl-tRNA synthetase restores mitochondrial dysfunction caused by a deafness-associated tRNAHis mutation.

The deafness-associated m.12201T>C mutation affects the A5-U68 base-pairing within the acceptor stem of mitochondrial tRNAHis The primary defect in this mutation is an alteration in tRNAHis aminoacylation. Here, we further investigate the molecular mechanism of the deafness-associated tRNAHis 12201T>C mutation and test whether the overexpression of the human mitochondrial histidyl-tRNA synthetase gene (HARS2) in cytoplasmic hybrid (cybrid) cells carrying the m.12201T>C mutation reverses mitochondrial dysfunctions. Using molecular dynamics simulations, we demonstrate that the m.12201T>C mutation perturbs the tRNAHis structure and function, supported by decreased melting temperature, conformational changes, and instability of mutated tRNA. We show that the m.12201T>C mutation-induced alteration of aminoacylation tRNAHis causes mitochondrial translational defects and respiratory deficiency. We found that the transfer of HARS2 into the cybrids carrying the m.12201T>C mutation raises the levels of aminoacylated tRNAHis from 56.3 to 75.0% but does not change the aminoacylation of other tRNAs. Strikingly, HARS2 overexpression increased the steady-state levels of tRNAHis and of noncognate tRNAs, including tRNAAla, tRNAGln, tRNAGlu, tRNALeu(UUR), tRNALys, and tRNAMet, in cells bearing the m.12201T>C mutation. This improved tRNA metabolism elevated the efficiency of mitochondrial translation, activities of oxidative phosphorylation complexes, and respiration capacity. Furthermore, HARS2 overexpression markedly increased mitochondrial ATP levels and membrane potential and reduced production of reactive oxygen species in cells carrying the m.12201T>C mutation. These results indicate that HARS2 overexpression corrects the mitochondrial dysfunction caused by the tRNAHis mutation. These findings provide critical insights into the pathophysiology of mitochondrial disease and represent a step toward improved therapeutic interventions for mitochondrial disorders.

Contribution of a mitochondrial tyrosyl-tRNA synthetase mutation to the phenotypic expression of the deafness-associated tRNASer(UCN) 7511A>G mutation.

Nuclear modifier genes have been proposed to modify the phenotypic expression of mitochondrial DNA mutations. Using a targeted exome-sequencing approach, here we found that the p.191Gly>Val mutation in mitochondrial tyrosyl-tRNA synthetase 2 (YARS2) interacts with the tRNASer(UCN) 7511A>G mutation in causing deafness. Strikingly, members of a Chinese family bearing both the YARS2 p.191Gly>Val and m.7511A>G mutations displayed much higher penetrance of deafness than those pedigrees carrying only the m.7511A>G mutation. The m.7511A>G mutation changed the A4:U69 base-pairing to G4:U69 pairing at the aminoacyl acceptor stem of tRNASer(UCN) and perturbed tRNASer(UCN) structure and function, including an increased melting temperature, altered conformation, instability, and aberrant aminoacylation of mutant tRNA. Using lymphoblastoid cell lines derived from symptomatic and asymptomatic members of these Chinese families and control subjects, we show that cell lines harboring only the m.7511A>G or p.191Gly>Val mutation revealed relatively mild defects in tRNASer(UCN) or tRNATyr metabolism, respectively. However, cell lines harboring both m.7511A>G and p.191Gly>Val mutations displayed more severe defective aminoacylations and lower tRNASer(UCN) and tRNATyr levels, aberrant aminoacylation, and lower levels of other tRNAs, including tRNAThr, tRNALys, tRNALeu(UUR), and tRNASer(AGY), than those in the cell lines carrying only the m.7511A>G or p.191Gly>Val mutation. Furthermore, mutant cell lines harboring both m.7511A>G and p.191Gly>Val mutations exhibited greater decreases in the levels of mitochondrial translation, respiration, and mitochondrial ATP and membrane potentials, along with increased production of reactive oxygen species. Our findings provide molecular-level insights into the pathophysiology of maternally transmitted deafness arising from the synergy between tRNASer(UCN) and mitochondrial YARS mutations.

Molecular dynamics simulations on apo ADP/ATP carrier shed new lights on the featured motif of the mitochondrial carriers.

The ADP/ATP carrier (AAC) is a transporter responsible for the equal molar exchange of cytosolic ADP and ATP synthesized within mitochondrial matrix across the mitochondrial membrane. Its primary structure consists of three homologous repeats, and each repeat contains a conserved motif that is shared by all members of the mitochondrial carrier family (MCF). Although these MCF motif residues cluster together in the crystal structure of AAC, detailed analyses on the interactions among the motif residues are still limited. In the present study, all-atom molecular dynamics (MD) simulations of up to 10 µs have been carried out on AAC, and interactions and structural dynamics of the MCF motif residues have been specifically investigated. Our simulations have revealed: i) a very asymmetrical electrostatic network at the bottom of the pocket of apo AAC, ii) the asymmetrical interactions between the Pro kink region and the [YWF][KR] G motif in three repeats, iii) the role of the conserved Arg residues in stabilizing the C-ends of the odd-numbered helices, iv) the structural change of the [YWF][KR] G motif and its potential involvement in substrate translocation process. Our results highlight the asymmetry of the MCF residues in the three repeats, which might contribute to the ability of the carriers to transport the asymmetrical substrates. Our observations provide microscopic basis for further research on the translocation mechanism of mitochondrial carriers.