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.

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.

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.

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.

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 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.

A coronary artery disease-associated tRNAThr mutation altered mitochondrial function, apoptosis and angiogenesis.

The tissue specificity of mitochondrial tRNA mutations remains largely elusive. In this study, we demonstrated the deleterious effects of tRNAThr 15927G>A mutation that contributed to pathogenesis of coronary artery disease. The m.15927G>A mutation abolished the highly conserved base-pairing (28C-42G) of anticodon stem of tRNAThr. Using molecular dynamics simulations, we showed that the m.15927G>A mutation caused unstable tRNAThr structure, supported by decreased melting temperature and slower electrophoretic mobility of mutated tRNA. Using cybrids constructed by transferring mitochondria from a Chinese family carrying the m.15927G>A mutation and a control into mitochondrial DNA (mtDNA)-less human umbilical vein endothelial cells, we demonstrated that the m.15927G>A mutation caused significantly decreased efficiency in aminoacylation and steady-state levels of tRNAThr. The aberrant tRNAThr metabolism yielded variable decreases in mtDNA-encoded polypeptides, respiratory deficiency, diminished membrane potential and increased the production of reactive oxygen species. The m.15927G>A mutation promoted the apoptosis, evidenced by elevated release of cytochrome c into cytosol and increased levels of apoptosis-activated proteins: caspases 3, 7, 9 and PARP. Moreover, the lower wound healing cells and perturbed tube formation were observed in mutant cybrids, indicating altered angiogenesis. Our findings provide new insights into the pathophysiology of coronary artery disease, which is manifested by tRNAThr mutation-induced alterations.

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.

Biochemical Evidence for a Nuclear Modifier Allele (A10S) in TRMU (Methylaminomethyl-2-thiouridylate-methyltransferase) Related to Mitochondrial tRNA Modification in the Phenotypic Manifestation of Deafness-associated 12S rRNA Mutation.

Nuclear modifier gene(s) was proposed to modulate the phenotypic expression of mitochondrial DNA mutation(s). Our previous investigations revealed that a nuclear modifier allele (A10S) in TRMU (methylaminomethyl-2-thiouridylate-methyltransferase) related to tRNA modification interacts with 12S rRNA 1555A→G mutation to cause deafness. The A10S mutation resided at a highly conserved residue of the N-terminal sequence. It was hypothesized that the A10S mutation altered the structure and function of TRMU, thereby causing mitochondrial dysfunction. Using molecular dynamics simulations, we showed that the A10S mutation introduced the Ser10 dynamic electrostatic interaction with the Lys106 residue of helix 4 within the catalytic domain of TRMU. The Western blotting analysis displayed the reduced levels of TRMU in mutant cells carrying the A10S mutation. The thermal shift assay revealed the Tm value of mutant TRMU protein, lower than that of the wild-type counterpart. The A10S mutation caused marked decreases in 2-thiouridine modification of U34 of tRNALys, tRNAGlu and tRNAGln However, the A10S mutation mildly increased the aminoacylated efficiency of tRNAs. The altered 2-thiouridine modification worsened the impairment of mitochondrial translation associated with the m.1555A→G mutation. The defective translation resulted in the reduced activities of mitochondrial respiration chains. The respiratory deficiency caused the reduction of mitochondrial ATP production and elevated the production of reactive oxidative species. As a result, mutated TRMU worsened mitochondrial dysfunctions associated with m.1555A→G mutation, exceeding the threshold for expressing a deafness phenotype. Our findings provided new insights into the pathophysiology of maternally inherited deafness that was manifested by interaction between mtDNA mutation and nuclear modifier gene.

Site-specific Disruption of the Oct4/Sox2 Protein Interaction Reveals Coordinated Mesendodermal Differentiation and the Epithelial-Mesenchymal Transition.

Although the Oct4/Sox2 complex is crucial for maintaining the pluripotency of stem cells, the molecular basis underlying its regulation during lineage-specific differentiation remains unknown. Here, we revealed that the highly conserved Oct4/Lys-156 is important for maintaining the stability of the Oct4 protein and the intermolecular salt bridge between Oct4/Lys-151 and Sox2/Asp-107 that contributes to the Oct4/Sox2 interaction. Post-translational modifications at Lys-156 and K156N, a somatic mutation detected in bladder cancer patients, both impaired the Lys-151-Asp-107 salt bridge and the Oct4/Sox2 interaction. When produced as a recombinant protein or overexpressed in pluripotent stem cells, Oct4/K156N, with reduced binding to Sox2, significantly down-regulated the stemness genes that are cooperatively controlled by the Oct4/Sox2 complex and specifically up-regulated the mesendodermal genes and the SNAIL family genes that promote the epithelial-mesenchymal transition. Thus, we conclude that Oct4/Lys-156-modulated Oct4/Sox2 interaction coordinately controls the epithelial-mesenchymal transition and mesendoderm specification induced by specific differentiation signals.

Cholesterol-ß1 AR interaction versus cholesterol-ß2 AR interaction.

Two 8-µs all-atom molecular dynamics simulations have been performed on the two highly homologous G protein-coupled receptor (GPCR) subtypes, ß1 - and ß2 -adrenergic receptors, which were embedded in a lipid bilayer with randomly dispersed cholesterol molecules. During the simulations, cholesterol molecules accumulate to different surface regions of the two receptors, suggesting the subtype specificity of cholesterol-ß-adrenergic receptor interaction and providing some clues to the physiological difference of the two subtypes. Meanwhile, comparison between the two receptors in interacting with cholesterols shed some new light on general determinants of cholesterol binding to GPCRs. Our results indicate that although the concave surface, charged residues and aromatic residues are important, neither of these stabilizing factors is indispensable for a cholesterol interaction site. Different combinations of these factors lead to the diversified binding modes of cholesterol binding to the receptors. Our long-time simulations, for the first time, revealed the pathway of a cholesterol molecule entering the consensus cholesterol motif (CCM) site, and the binding process of cholesterol to CCM is accompanied by a side chain flipping of the conserved Trp4.50. Moreover, the simulation results suggest that the I-/V-/L-rich region on the extracellular parts of helix 6 might be an alternatively conserved cholesterol-binding site for the class-A GPCRs.

Relative stability of different DNA guanine quadruplex stem topologies derived using large-scale quantum-chemical computations.

We provide theoretical predictions of the intrinsic stability of different arrangements of guanine quadruplex (G-DNA) stems. Most computational studies of nucleic acids have applied Molecular Mechanics (MM) approaches using simple pairwise-additive force fields. The principle limitation of such calculations is the highly approximate nature of the force fields. In this study, we for the first time apply accurate QM computations (DFT-D3 with large atomic orbital basis sets) to essentially complete DNA building blocks, seven different folds of the cation-stabilized two-quartet G-DNA stem, each having more than 250 atoms. The solvent effects are approximated by COSMO continuum solvent. We reveal sizable differences between MM and QM descriptions of relative energies of different G-DNA stems, which apparently reflect approximations of the DNA force field. Using the QM energy data, we propose correction to earlier free energy estimates of relative stabilities of different parallel, hybrid, and antiparallel G-stem folds based on classical simulations. The new energy ranking visibly improves the agreement between theory and experiment. We predict the 5'-anti-anti-3' GpG dinucleotide step to be the most stable one, closely followed by the 5'-syn-anti-3' step. The results are in good agreement with known experimental structures of 2-, 3-, and 4-quartet G-DNA stems. Besides providing specific results for G-DNA, our study highlights basic limitations of force field modeling of nucleic acids. Although QM computations have their own limitations, mainly the lack of conformational sampling and the approximate description of the solvent, they can substantially improve the quality of calculations currently relying exclusively on force fields.

Molecular dynamics simulations of G-DNA and perspectives on the simulation of nucleic acid structures.

The article reviews the application of biomolecular simulation methods to understand the structure, dynamics and interactions of nucleic acids with a focus on explicit solvent molecular dynamics simulations of guanine quadruplex (G-DNA and G-RNA) molecules. While primarily dealing with these exciting and highly relevant four-stranded systems, where recent and past simulations have provided several interesting results and novel insight into G-DNA structure, the review provides some general perspectives on the applicability of the simulation techniques to nucleic acids.

Explaining the varied glycosidic conformational, G-tract length and sequence preferences for anti-parallel G-quadruplexes.

Guanine-rich DNA sequences tend to form four-stranded G-quadruplex structures. Characteristic glycosidic conformational patterns along the G-strands, such as the 5'-syn-anti-syn-anti pattern observed with the Oxytricha nova telomeric G-quadruplexes, have been well documented. However, an explanation for these featured glycosidic patterns has not emerged. This work presents MD simulation and free energetic analyses for simplified two-quartet [d(GG)](4) models and suggests that the four base pair step patterns show quite different relative stabilities: syn-anti > anti-anti > anti-syn > syn-syn. This suggests the following rule: when folding, anti-parallel G-quadruplexes tend to maximize the number of syn-anti steps and avoid the unfavorable anti-syn and syn-syn steps. This rule is consistent with most of the anti-parallel G-quadruplex structures in the Protein Databank (PDB). Structural polymorphisms of G-quadruplexes relate to these glycosidic conformational patterns and the lengths of the G-tracts. The folding topologies of G2- and G4-tracts are not very polymorphic because each strand tends to populate the stable syn-anti repeat. G3-tracts, on the other hand, cannot present this repeating pattern on each G-tract. This leads to smaller energy differences between different geometries and helps explain the extreme structural polymorphism of the human telomeric G-quadruplexes.

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.

Conformational transition and energy landscape of ErbB4 activated by neuregulin1ß: one microsecond molecular dynamics simulations.

ErbB4, a receptor tyrosine kinase of the ErbB family, plays crucial roles in cell growth and differentiation, especially in the development of the heart and nervous system. Ligand binding to its extracellular region could modulate the activation process. To understand the mechanism of ErbB4 activation induced by ligand binding, we performed one microsecond molecular dynamics (MD) simulations on the ErbB4 extracellular region (ECR) with and without its endogenous ligand neuregulin1ß (NRG1ß). The conformational transition of the ECR-ErbB4/NRG1ß complex from a tethered inactive conformation to an extended active-like form has been observed, while such large and function-related conformational change has not been seen in the simulation on the ECR-ErbB4, suggesting that ligand binding is indeed the active inducing force for the conformational transition and further dimerization. On the basis of MD simulations and principal component analysis, we constructed a rough energy landscape for the conformational transition of ECR-ErbB4/NRG1ß complex, suggesting that the conformational change from the inactive state to active-like state involves a stable conformation. The energy barrier for the tether opening was estimated as ~2.7 kcal/mol, which is very close to the experimental value (1-2 kcal/mol) reported for ErbB1. On the basis of the simulation results, an atomic mechanism for the ligand-induced activation of ErbB4 was postulated. The present MD simulations provide a new insight into the conformational changes underlying the activation of ErbB4.

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.

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.

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.

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.