Mutational burden of XPNPEP3 leads to defects in mitochondrial complex I and cilia in NPHPL1.

Nephronophthisis-like nephropathy-1 (NPHPL1) is a rare ciliopathy, caused by mutations of XPNPEP3. Despite a well-described monogenic etiology, the pathogenesis of XPNPEP3 associated with mitochondrial and ciliary function remains elusive. Here, we identified novel compound heterozygous mutations in NPHPL1 patients with renal lesion only or with extra bone cysts together. Patient-derived lymphoblasts carrying c.634G>A and c.761G>T together exhibit elevated mitochondrial XPNPEP3 levels via the reduction of mRNA degradation, leading to mitochondrial dysfunction in both urine tubular epithelial cells and lymphoblasts from patient. Mitochondrial XPNPEP3 was co-immunoprecipitated with respiratory chain complex I and was required for the stability and activity of complex I. Deletion of Xpnpep3 in mice resulted in lower activity of complex I, elongated primary cilium, and predisposition to tubular dilation and fibrosis under stress. Our findings provide valuable insights into the mitochondrial functions involved in the pathogenesis of NPHP.

A novel NF2 splicing mutant causes neurofibromatosis type 2 via liquid-liquid phase separation with large tumor suppressor and Hippo pathway.

Neurofibromatosis type 2 is an autosomal dominant multiple neoplasia syndrome and is usually caused by mutations in the neurofibromin 2 (NF2) gene, which encodes a tumor suppressor and initiates the Hippo pathway. However, the mechanism by which NF2 functions in the Hippo pathway isn't fully understood. Here we identified a NF2 c.770-784del mutation from a neurofibromatosis type 2 family. MD simulations showed that this mutation significantly changed the structure of the F3 module of the NF2-FERM domain. Functional assays indicated that the NF2 c.770-784del variant formed LLPS in the cytoplasm with LATS to restrain LATS plasma membrane localization and inactivated the Hippo pathway. Besides, this deletion partly caused a skipping of exon 8 and reduced the protein level of NF2, collectively promoting proliferation and tumorigenesis of meningeal cells. We identified an unrecognized mechanism of LLPS and splicing skipping for the NF2-induced Hippo pathway, which provided new insight into the pathogenesis of neurofibromatosis type 2.

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.

Heteroplasmic and homoplasmic m.616T>C in mitochondria tRNAPhe promote isolated chronic kidney disease and hyperuricemia.

Inherited kidney diseases are the fifth most common cause of end-stage renal disease (ESRD). Mitochondrial dysfunction plays a vital role in the progression of inherited kidney diseases, while mitochondrial-transfer RNA (mt-tRNA) variants and their pathogenic contributions to kidney disease remain largely unclear. In this study, we identified the pathogenic mt-tRNAPhe 616T>C mutation in 3 families and documented that m.616T>C showed a high pathogenic threshold, with both heteroplasmy and homoplasmy leading to isolated chronic kidney disease and hyperuricemia without hematuria, proteinuria, or renal cyst formation. Moreover, 1 proband with homoplamic m.616T>C presented ESRD as a child. No symptoms of nervous system evolvement were observed in these families. Lymphoblast cells bearing m.616T>C exhibited swollen mitochondria, underwent active mitophagy, and showed respiratory deficiency, leading to reduced mitochondrial ATP production, diminished membrane potential, and overproduction of mitochondrial ROS. Pathogenic m.616T>C abolished a highly conserved base pair (A31-U39) in the anticodon stem-loop which altered the structure of mt-tRNAPhe, as confirmed by a decreased melting temperature and slower electrophoretic mobility of the mutant tRNA. Furthermore, the unstable structure of mt-tRNAPhe contributed to a shortage of steady-state mt-tRNAPhe and enhanced aminoacylation efficiency, which resulted in impaired mitochondrial RNA translation and a significant decrease in mtDNA-encoded polypeptides. Collectively, these findings provide potentially new insights into the pathogenesis underlying inherited kidney disease caused by mitochondrial variants.

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.

Long-term prognosis of 35 patients with methionine adenosyltransferase deficiency based on newborn screening in China.

Methionine adenosyltransferase deficiency (MATD) is a rare metabolic disorder caused by mono- or biallelic MAT1A mutations that are not yet well understood. Of the 4,065,644 neonates screened between November 2010 and December 2021, 35 individuals have been diagnosed with an estimated incidence of 1: 116,161 by a cutoff value of methionine 82.7 µmol/L and follow-up over 11 years. MATD patients with autosomal recessive (AR) type had higher clinical and genetic heterogeneity than those with autosomal dominant (AD) type. Fifteen unrelated AD patients harbored one well-known dominant variant, c.791 G>A or c.776 C>T, and were clinically unaffected with a mean plasma methionine (Met) value <300 µmol/L. Twenty AR cases have unique genotypes and presented a wide range of clinical abnormalities from asymptomatic to white matter lesions. Of them, 10 AR patients displayed severe manifestations, such as verbal difficulty, motor delay, development delay, and white matter lesions, with mean Met >500 µmol/L and thereby were treated with a methionine-restricted diet alone or in combination with betaine, folate, or vitamin B6, and were healthy finally. Neurological abnormalities were evidenced in two patients (P16 and P27) with Met values >800 µmol/L by MRI scan. Neurological abnormalities were reversed here by liver transplantation or by the determination of S-adenosylmethionine supplementation. Additionally, 38 variants of MAT1A were distributed within patients and carriers, of which 24 were novel and mostly predicted to be damaged. Our findings with an extensive clinical and genetic dataset provided new insights into its diagnosis and treatment and will be helpful for its optimal management in the future.

The low contagiousness and new A958D mutation of SARS-CoV-2 in children: An observational cohort study.

AIMS: To explore the contagiousness and new SARS-CoV-2 mutations in pediatric COVID-19. METHODS: This cohort study enrolled all pediatric patients admitted to 8 hospitals in Zhejiang Province of China between 21 January and 29 February 2020, their family members and close-contact classmates. Epidemiological, demographic, clinical and laboratory data were collected. Bioinformatics was used to analyze the features of SARS-CoV-2. Individuals were divided into 3 groups by the first-generation case: Groups 1 (unclear), 2 (adult), and 3 (child). The secondary attack rate (SAR) and R0 were compared among the groups. RESULTS: The infection rate among 211 individuals was 64% (135/211). The SAR in Groups 2 and 3 was 71% (73/103) and 3% (1/30), respectively; the median R0 in Groups 2 and 3 was 2 (range: 1-8) and 0 (range: 0-1), respectively. Compared with adult cases, the SAR and R0 of pediatric cases were significantly lower (p<0.05). We obtained SARS-CoV-2 sequences from the same infant's throat and fecal samples at a two-month interval and found that the new spike protein A958D mutation detected in the stool improved thermostability theoretically. CONCLUSIONS: Children have lower ability to spread SARS-CoV-2. The new A958D mutation is a potential reason for its long residence in the intestine.

A deafness-associated tRNA mutation caused pleiotropic effects on the m1G37 modification, processing, stability and aminoacylation of tRNAIle and mitochondrial translation.

Defects in the posttranscriptional modifications of mitochondrial tRNAs have been linked to human diseases, but their pathophysiology remains elusive. In this report, we investigated the molecular mechanism underlying a deafness-associated tRNAIle 4295A>G mutation affecting a highly conserved adenosine at position 37, 3' adjacent to the tRNA's anticodon. Primer extension and methylation activity assays revealed that the m.4295A>G mutation introduced a tRNA methyltransferase 5 (TRMT5)-catalyzed m1G37 modification of tRNAIle. Molecular dynamics simulations suggested that the m.4295A>G mutation affected tRNAIle structure and function, supported by increased melting temperature, conformational changes and instability of mutated tRNA. An in vitro processing experiment revealed that the m.4295A>G mutation reduced the 5' end processing efficiency of tRNAIle precursors, catalyzed by RNase P. We demonstrated that cybrid cell lines carrying the m.4295A>G mutation exhibited significant alterations in aminoacylation and steady-state levels of tRNAIle. The aberrant tRNA metabolism resulted in the impairment of mitochondrial translation, respiratory deficiency, decreasing membrane potentials and ATP production, increasing production of reactive oxygen species and promoting autophagy. These demonstrated the pleiotropic effects of m.4295A>G mutation on tRNAIle and mitochondrial functions. Our findings highlighted the essential role of deficient posttranscriptional modifications in the structure and function of tRNA and their pathogenic consequence of deafness.

Phosphorothioate Substitutions in RNA Structure Studied by Molecular Dynamics Simulations, QM/MM Calculations, and NMR Experiments.

Phosphorothioates (PTs) are important chemical modifications of the RNA backbone where a single nonbridging oxygen of the phosphate is replaced with a sulfur atom. PT can stabilize RNAs by protecting them from hydrolysis and is commonly used as a tool to explore their function. It is, however, unclear what basic physical effects PT has on RNA stability and electronic structure. Here, we present molecular dynamics (MD) simulations, quantum mechanical (QM) calculations, and NMR spectroscopy measurements, exploring the effects of PT modifications in the structural context of the neomycin-sensing riboswitch (NSR). The NSR is the smallest biologically functional riboswitch with a well-defined structure stabilized by a U-turn motif. Three of the signature interactions of the U-turn: an H-bond, an anion-π interaction, and a potassium binding site; are formed by RNA phosphates, making the NSR an ideal model for studying how PT affects RNA structure and dynamics. By comparing with high-level QM calculations, we reveal the distinct physical properties of the individual interactions facilitated by the PT. The sulfur substitution, besides weakening the direct H-bond interaction, reduces the directionality of H-bonding while increasing its dispersion and induction components. It also reduces the induction and increases the dispersion component of the anion-π stacking. The sulfur force-field parameters commonly employed in the literature do not reflect these distinctions, leading to the unsatisfactory description of PT in simulations of the NSR. We show that it is not possible to accurately describe the PT interactions using one universal set of van der Waals sulfur parameters and provide suggestions for improving the force-field performance.

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.

O-GlcNAcylation of PGK1 coordinates glycolysis and TCA cycle to promote tumor growth.

Many cancer cells display enhanced glycolysis and suppressed mitochondrial metabolism. This phenomenon, known as the Warburg effect, is critical for tumor development. However, how cancer cells coordinate glucose metabolism through glycolysis and the mitochondrial tricarboxylic acid (TCA) cycle is largely unknown. We demonstrate here that phosphoglycerate kinase 1 (PGK1), the first ATP-producing enzyme in glycolysis, is reversibly and dynamically modified with O-linked N-acetylglucosamine (O-GlcNAc) at threonine 255 (T255). O-GlcNAcylation activates PGK1 activity to enhance lactate production, and simultaneously induces PGK1 translocation into mitochondria. Inside mitochondria, PGK1 acts as a kinase to inhibit pyruvate dehydrogenase (PDH) complex to reduce oxidative phosphorylation. Blocking T255 O-GlcNAcylation of PGK1 decreases colon cancer cell proliferation, suppresses glycolysis, enhances the TCA cycle, and inhibits tumor growth in xenograft models. Furthermore, PGK1 O-GlcNAcylation levels are elevated in human colon cancers. This study highlights O-GlcNAcylation as an important signal for coordinating glycolysis and the TCA cycle to promote tumorigenesis.

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 mechanism of Bisphenol A on androgen receptor antagonism.

Bisphenol A (BPA), one of the highest production volume chemicals, is a typical endocrine-disrupting chemical (EDC) that exhibits antiandrogenic activity. However, how BPA antagonizes androgen effects remains ambiguous. In this study, the in silico and in vitro assays were carried out to explore the molecular mechanism(s) of BPA on androgen receptor (AR) antagonism. In reporter gene assay, BPA caused a significant antagonistic effect on 5α-dihydrotestosterone (DHT)-induced AR transcriptional activity at concentrations of 10-9 M-10-5 M. The results of molecular docking and molecular dynamics simulations indicated the availability of BPA binding to the ligand binding domain of AR. BPA treatment prevented the inhibition of receptor degradation caused by DHT binding to AR. BPA exposure also abolished DHT-dependent dissociation of AR from its co-chaperone, 90-kDa heat shock protein (Hsp90), and resulted in the blockage of DHT-induced AR nuclear translocation. This is the first report to show that BPA inhibited the DHT-induced stabilization of AR and the DHT-induced dissociation of AR-Hsp90 complex. This study provided new evidence for further understanding the precise mechanisms of antagonism of BPA on AR.

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.

Contribution of mitochondrial ND1 3394T>C mutation to the phenotypic manifestation of Leber's hereditary optic neuropathy.

Mitochondrial DNA (mtDNA) mutations have been associated with Leber's hereditary optic neuropathy (LHON) and their pathophysiology remains poorly understood. In this study, we investigated the pathophysiology of a LHON susceptibility allele (m.3394T>C, p.30Y>H) in the Mitochondrial (MT)-ND1 gene. The incidence of m.3394T>C mutation was 2.7% in the cohort of 1741 probands with LHON. Extremely low penetrances of LHON were observed in 26 pedigrees carrying only m.3394T>C mutation, while 21 families bearing m.3394T>C, together with m.11778G>A or m.14484T>C mutation, exhibited higher penetrance of LHON than those in families carrying single mtDNA mutation(s). The m.3394T>C mutation disrupted the specific electrostatic interactions between Y30 of p.MT-ND1 with the sidechain of E4 and backbone carbonyl group of M1 of NDUFA1 (NADH dehydrogenase [ubiquinone] 1 alpha subcomplex subunit 1) of complex I, thereby altering the structure and function of complex I. We demonstrated that these cybrids bearing only m.3394T>C mutation caused mild mitochondrial dysfunctions and those harboring both m.3394T>C and m.11778G>A mutations exhibited greater mitochondrial dysfunctions than cybrids carrying only m.11778G>A mutation. In particular, the m.3394T>C mutation altered the stability of p.MT-ND1 and complex I assembly. Furthermore, the m.3394T>C mutation decreased the activities of mitochondrial complexes I, diminished mitochondrial ATP levels and membrane potential and increased the production of reactive oxygen species in the cybrids. These m.3394T>C mutation-induced alterations aggravated mitochondrial dysfunctions associated with the m.11778G>A mutation. These resultant biochemical defects contributed to higher penetrance of LHON in these families carrying both mtDNA mutations. Our findings provide new insights into the pathophysiology of LHON arising from the synergy between mitochondrial ND1 and ND4 mutations.