Defining the core proteome of the chloroplast envelope membranes

High-throughput protein localization studies require multiple strategies. Mass spectrometric analysis of defined cellular fractions is one of the complementary approaches to a diverse array of cell biological methods. In recent years, the protein content of different cellular (sub-)compartments was approached. Despite of all the efforts made, the analysis of membrane fractions remains difficult, in that the dissection of the proteomes of the envelope membranes of chloroplasts or mitochondria is often not reliable because sample purity is not always warranted. Moreover, proteomic studies are often restricted to single (model) species, and therefore limited in respect to differential individual evolution. In this study we analyzed the chloroplast envelope proteomes of different plant species, namely, the individual proteomes of inner and outer envelope (OE) membrane of Pisum sativum and the mixed envelope proteomes of Arabidopsis thaliana and Medicago sativa. The analysis of all three species yielded 341 identified proteins in total, 247 of them being unique. 39 proteins were genuine envelope proteins found in at least two species. Based on this and previous envelope studies we defined the core envelope proteome of chloroplasts. Comparing the general overlap of the available six independent studies (including ours) revealed only a number of 27 envelope proteins. Depending on the stringency of applied selection criteria we found 231 envelope proteins, while less stringent criteria increases this number to 649 putative envelope proteins. Based on the latter we provide a map of the outer and inner envelope core proteome, which includes many yet uncharacterized proteins predicted to be involved in transport, signaling, and response. Furthermore, a foundation for the functional characterization of yet unidentified functions of the inner and OE for further analyses is provided

A systems approach delivers a functional microRNA catalog and expanded targets for seizure suppression in temporal lobe epilepsy

Temporal lobe epilepsy is the most common drug-resistant form of epilepsy in adults. The reorganization of neural networks and the gene expression landscape underlying pathophysiologic network behavior in brain structures such as the hippocampus has been suggested to be controlled, in part, by microRNAs. To systematically assess their significance, we sequenced Argonaute-loaded microRNAs to define functionally engaged microRNAs in the hippocampus of three different animal models in two species and at six time points between the initial precipitating insult through to the establishment of chronic epilepsy. We then selected commonly up-regulated microRNAs for a functional in vivo therapeutic screen using oligonucleotide inhibitors. Argonaute sequencing generated 1.44 billion small RNA reads of which up to 82% were microRNAs, with over 400 unique microRNAs detected per model. Approximately half of the detected microRNAs were dysregulated in each epilepsy model. We prioritized commonly up-regulated microRNAs that were fully conserved in humans and designed custom antisense oligonucleotides for these candidate targets. Antiseizure phenotypes were observed upon knockdown of miR-10a-5p, miR-21a-5p, and miR-142a-5p and electrophysiological analyses indicated broad safety of this approach. Combined inhibition of these three microRNAs reduced spontaneous seizures in epileptic mice. Proteomic data, RNA sequencing, and pathway analysis on predicted and validated targets of these microRNAs implicated derepressed TGF-\u3b2 signaling as a shared seizure-modifying mechanism. Correspondingly, inhibition of TGF-\u3b2 signaling occluded the antiseizure effects of the antagomirs. Together, these results identify shared, dysregulated, and functionally active microRNAs during the pathogenesis of epilepsy which represent therapeutic antiseizure targets

Gene variant effects across sodium channelopathies predict function and guide precision therapy

Pathogenic variants in the voltage-gated sodium channel gene family lead to early onset epilepsies, neurodevelopmental disorders, skeletal muscle channelopathies, peripheral neuropathies and cardiac arrhythmias. Disease-associated variants have diverse functional effects ranging from complete loss-of-function to marked gain-of-function. Therapeutic strategy is likely to depend on functional effect. Experimental studies offer important insights into channel function but are resource intensive and only performed in a minority of cases. Given the evolutionarily conserved nature of the sodium channel genes, we investigated whether similarities in biophysical properties between different voltage-gated sodium channels can predict function and inform precision treatment across sodium channelopathies. We performed a systematic literature search identifying functionally assessed variants in any of the nine voltage-gated sodium channel genes until 28 April 2021. We included missense variants that had been electrophysiologically characterized in mammalian cells in whole-cell patch-clamp recordings. We performed an alignment of linear protein sequences of all sodium channel genes and correlated variants by their overall functional effect on biophysical properties. Of 951 identified records, 437 sodium channel-variants met our inclusion criteria and were reviewed for functional properties. Of these, 141 variants were epilepsy-associated (SCN1/2/3/8A), 79 had a neuromuscular phenotype (SCN4/9/10/11A), 149 were associated with a cardiac phenotype (SCN5/10A) and 68 (16%) were considered benign. We detected 38 missense variant pairs with an identical disease-associated variant in a different sodium channel gene. Thirty-five out of 38 of those pairs resulted in similar functional consequences, indicating up to 92% biophysical agreement between corresponding sodium channel variants (odds ratio = 11.3; 95% confidence interval = 2.8 to 66.9; P < 0.001). Pathogenic missense variants were clustered in specific functional domains, whereas population variants were significantly more frequent across non-conserved domains (odds ratio = 18.6; 95% confidence interval = 10.9-34.4; P < 0.001). Pore-loop regions were frequently associated with loss-of-function variants, whereas inactivation sites were associated with gain-of-function (odds ratio = 42.1, 95% confidence interval = 14.5-122.4; P < 0.001), whilst variants occurring in voltage-sensing regions comprised a range of gain- and loss-of-function effects. Our findings suggest that biophysical characterisation of variants in one SCN-gene can predict channel function across different SCN-genes where experimental data are not available. The collected data represent the first gain- versus loss-of-function topological map of SCN proteins indicating shared patterns of biophysical effects aiding variant analysis and guiding precision therapy. We integrated our findings into a free online webtool to facilitate functional sodium channel gene variant interpretation (http://SCN-viewer.broadinstitute.org).Peer reviewe

Molecular dynamics simulations reveal molecular mechanisms for the gain and loss of function effects of four SCN2A variants

Genetic variants in the SCN2A gene cause a wide spectrum of neurodevelopmental disorders. The current binary classification of pathogenic missense variants as either gain or loss of function does not correlate with the diverse clinical presentations observed in affected individuals. We hypothesize that missense variants lead to a spectrum of disease mechanisms and that computational extensive molecular dynamics (MD) simulations can uncover these. We have selected four variants within the SCN2A gene based on their high prevalence among SCN2A patients and distinct molecular read-outs, classified as either gain or loss of function based upon in vitro demonstrated increase or decrease in ion flux. We designed five systems, employed 500ns MD simulations with 3 replicas each, to investigate a comprehensive range of biophysical properties exhibited by disease-causing variants at the atomic level. All patient variant structure models were less tightly packed and exhibited a greater range of conformations compared to the wild type, suggesting that all mutants undergo structural changes. In addition, we also studied the possible variant specific disease mechanisms responsible for all four variants. These findings go beyond the binary classification of gain and loss of function obtained from electrophysiological studies. We identified two different mechanisms resulting in a loss of function and two different mechanisms resulting in a gain of function, therefore providing greater detail and understanding of the variant effect, which also suggests that categorization of a variant solely as a gain or loss of function oversimplifies the complex nature of these mechanisms. These molecular dynamics (MD) studies serve as the initial step in examining numerous variants in sodium ion channel genes, and with adequate computing resources, this approach could be applied to analyze all missense variants

Conserved patterns across ion channels correlate with variant pathogenicity and clinical phenotypes

Clinically identified genetic variants in ion channels can be benign or cause disease by increasing or decreasing the protein function. Consequently, therapeutic decision-making is challenging without molecular testing of each variant. Our biophysical knowledge of ion channel structures and function is just emerging, and it is currently not well understood which amino acid residues cause disease when mutated.We sought to systematically identify biological properties associated with variant pathogenicity across all major voltage and ligand-gated ion channel families. We collected and curated 3,049 pathogenic variants from hundreds of neurodevelopmental and other disorders and 12,546 population variants for 30 ion channel or channel subunits for which a high-quality protein structure was available. Using a wide range of bioinformatics approaches, we computed 163 structural features and tested them for pathogenic variant enrichment. We developed a novel 3D spatial distance scoring approach that enables comparisons of pathogenic and population variant distribution across protein structures.We discovered and independently replicated that several pore residue properties and proximity to the pore axis were most significantly enriched for pathogenic variants compared to population variants. Using our 3D scoring approach, we showed that the strongest pathogenic variant enrichment was observed for pore-lining residues and alpha-helix residues within 5Å distance from the pore axis center and not involved in gating. Within the subset of residues located at the pore, the hydrophobicity of the pore was the feature most strongly associated with variant pathogenicity. We also found an association between the identified properties and both clinical phenotypes and functional in vitro assays for voltage-gated sodium channels (SCN1A, SCN2A, SCN8A) and N-methyl-D-aspartate (NMDA) receptor (GRIN1, GRIN2A, GRIN2B) encoding genes. In an independent expert-curated dataset of 1,422 neurodevelopmental disorder pathogenic patient variants and 679 electrophysiological experiments, we show that pore axis distance is associated with seizure age of onset and cognitive performance as well as differential gain vs. loss-of-channel function.In summary, we identified biological properties associated with ion-channel malfunction and show that these are correlated with in vitro functional read-outs and clinical phenotypes in patients with neurodevelopmental disorders. Our results suggest that clinical decision support algorithms that predict variant pathogenicity and function are feasible in the future

Conserved patterns across ion channels correlate with variant pathogenicity and clinical phenotypes 2022.03.23.485339

Clinically identified genetic variants in ion channels can be benign or cause disease by increasing or decreasing the protein function. Consequently, therapeutic decision-making is challenging without molecular testing of each variant. Our biophysical knowledge of ion channel structures and function is just emerging, and it is currently not well understood which amino acid residues cause disease when mutated.We sought to systematically identify biological properties associated with variant pathogenicity across all major voltage and ligand-gated ion channel families. We collected and curated 3,049 pathogenic variants from hundreds of neurodevelopmental and other disorders and 12,546 population variants for 30 ion channel or channel subunits for which a high-quality protein structure was available. Using a wide range of bioinformatics approaches, we computed 163 structural features and tested them for pathogenic variant enrichment. We developed a novel 3D spatial distance scoring approach that enables comparisons of pathogenic and population variant distribution across protein structures.We discovered and independently replicated that several pore residue properties and proximity to the pore axis were most significantly enriched for pathogenic variants compared to population variants. Using our novel 3D scoring approach, we showed that the strongest pathogenic variant enrichment was observed for pore-lining residues and alpha-helix residues within 5 A distance from the pore axis center and not involved in gating. Within the subset of residues located at the pore, the hydrophobicity of the pore was the feature most strongly associated with variant pathogenicity. We also found an association between the identified properties and both clinical phenotypes and fucntional in vitro assays for voltage-gated sodium channels (SCN1A, SCN2A, SCN8A) and N-methyl-D-aspartate (NMDA) receptor (GRIN1, GRIN2A, GRIN2B) encoding genes. In an independent expert-curated dataset of 1,422 neurodevelopmental disorder pathogenic patient variants, and 679 electrophysiological experiments that pore axis distance is associated with seizure age of onset and cognitive performance as well as differential gain vs. loss-of-channel function.In summary, we identified biological properties associated with ion-channel malfunction and show that these are correlated with in vitro functional read-outs and clinical phenotypes in patients with neurodevelopmental disorders. Our results suggest that clinical decision support algorithms that predict variant pathogenicity and function are feasible in the future.Competing Interest StatementThe authors have declared no competing interest.DSSPDictionary of Protein Secondary StructuregnomADGenome aggregation DatabaseGoFGain of functionGRIN genesGRIN1, GRIN2A. GRIN2BHGMDHuman Gene Mutation DatabaseNMDA receptorN-methyl-D-aspartate receptorGABA receptorGamma-aminobutyric acid receptorLoFLoss of functionSCN genesSCN1A, SCN2A, SCN8AVCFVariant Call Forma

The association of late-acting snoRNPs with human pre-ribosomal complexes requires the RNA helicase DDX21

Translation fidelity and efficiency require multiple ribosomal (r)RNA modifications that are mostly mediated by small nucleolar (sno)RNPs during ribosome production. Overlapping basepairing of snoRNAs with pre-rRNAs often necessitates sequential and efficient association and dissociation of the snoRNPs, however, how such hierarchy is established has remained unknown so far. Here, we identify several late-acting snoRNAs that bind pre-40S particles in human cells and show that their association and function in pre-40S complexes is regulated by the RNA helicase DDX21. We map DDX21 crosslinking sites on pre-rRNAs and show their overlap with the basepairing sites of the affected snoRNAs. While DDX21 activity is required for recruitment of the late-acting snoRNAs SNORD56 and SNORD68, earlier snoRNAs are not affected by DDX21 depletion. Together, these observations provide an understanding of the timing and ordered hierarchy of snoRNP action in pre-40S maturation and reveal a novel mode of regulation of snoRNP function by an RNA helicase in human cells

Antagonizing Increased miR-135a Levels at the Chronic Stage of Experimental TLE Reduces Spontaneous Recurrent Seizures

Mesial temporal lobe epilepsy (mTLE) is a chronic neurological disease characterized by recurrent seizures. The antiepileptic drugs currently available to treat mTLE are ineffective in one-third of patients and lack disease-modifying effects. miRNAs, a class of small noncoding RNAs which control gene expression at the post-transcriptional level, play a key role in the pathogenesis of mTLE and other epilepsies. Although manipulation of miRNAs at acute stages has been reported to reduce subsequent spontaneous seizures, it is uncertain whether targeting miRNAs at chronic stages of mTLE can also reduce seizures. Furthermore, the functional role and downstream targets of most epilepsy-associated miRNAs remain poorly understood. Here, we show that miR-135a is selectively upregulated within neurons in epileptic brain and report that targeting miR-135a in vivo using antagomirs after onset of spontaneous recurrent seizures can reduce seizure activity at the chronic stage of experimental mTLE in male mice. Further, by using an unbiased approach combining immunoprecipitation and RNA sequencing, we identify several novel neuronal targets of miR-135a, including Mef2a Mef2 proteins are key regulators of excitatory synapse density. Mef2a and miR-135a show reciprocal expression regulation in human (of both sexes) and experimental TLE, and miR-135a regulates dendritic spine number and type through Mef2. Together, our data show that miR-135a is target for reducing seizure activity in chronic epilepsy, and that deregulation of miR-135a in epilepsy may alter Mef2a expression and thereby affect synaptic function and plasticity.SIGNIFICANCE STATEMENT miRNAs are post-transcriptional regulators of gene expression with roles in the pathogenesis of epilepsy. However, the precise mechanism of action and therapeutic potential of most epilepsy-associated miRNAs remain poorly understood. Our study reveals dramatic upregulation of the key neuronal miRNA miR-135a in both experimental and human mesial temporal lobe epilepsy. Silencing miR-135a in experimental temporal lobe epilepsy reduces seizure activity at the spontaneous recurrent seizure stage. These data support the exciting possibility that miRNAs can be targeted to combat seizures after spontaneous seizure activity has been established. Further, by using unbiased approaches novel neuronal targets of miR-135a, including members of the Mef2 protein family, are identified that begin to explain how deregulation of miR-135a may contribute to epilepsy

Antagonizing Increased miR-135a Levels at the Chronic Stage of Experimental TLE Reduces Spontaneous Recurrent Seizures

Mesial temporal lobe epilepsy (mTLE) is a chronic neurological disease characterized by recurrent seizures. The antiepileptic drugs currently available to treat mTLE are ineffective in one-third of patients and lack disease-modifying effects. miRNAs, a class of small noncoding RNAs which control gene expression at the post-transcriptional level, play a key role in the pathogenesis of mTLE and other epilepsies. Although manipulation of miRNAs at acute stages has been reported to reduce subsequent spontaneous seizures, it is uncertain whether targeting miRNAs at chronic stages of mTLE can also reduce seizures. Furthermore, the functional role and downstream targets of most epilepsy-associated miRNAs remain poorly understood. Here, we show that miR-135a is selectively upregulated within neurons in epileptic brain and report that targeting miR-135a in vivo using antagomirs after onset of spontaneous recurrent seizures can reduce seizure activity at the chronic stage of experimental mTLE in male mice. Further, by using an unbiased approach combining immunoprecipitation and RNA sequencing, we identify several novel neuronal targets of miR-135a, including Mef2a Mef2 proteins are key regulators of excitatory synapse density. Mef2a and miR-135a show reciprocal expression regulation in human (of both sexes) and experimental TLE, and miR-135a regulates dendritic spine number and type through Mef2. Together, our data show that miR-135a is target for reducing seizure activity in chronic epilepsy, and that deregulation of miR-135a in epilepsy may alter Mef2a expression and thereby affect synaptic function and plasticity.SIGNIFICANCE STATEMENT miRNAs are post-transcriptional regulators of gene expression with roles in the pathogenesis of epilepsy. However, the precise mechanism of action and therapeutic potential of most epilepsy-associated miRNAs remain poorly understood. Our study reveals dramatic upregulation of the key neuronal miRNA miR-135a in both experimental and human mesial temporal lobe epilepsy. Silencing miR-135a in experimental temporal lobe epilepsy reduces seizure activity at the spontaneous recurrent seizure stage. These data support the exciting possibility that miRNAs can be targeted to combat seizures after spontaneous seizure activity has been established. Further, by using unbiased approaches novel neuronal targets of miR-135a, including members of the Mef2 protein family, are identified that begin to explain how deregulation of miR-135a may contribute to epilepsy

Antagonizing Increased miR-135a Levels at the Chronic Stage of Experimental TLE Reduces Spontaneous Recurrent Seizures

Mesial temporal lobe epilepsy (mTLE) is a chronic neurological disease characterized by recurrent seizures. The antiepileptic drugs currently available to treat mTLE are ineffective in one-third of patients and lack disease-modifying effects. miRNAs, a class of small noncoding RNAs which control gene expression at the post-transcriptional level, play a key role in the pathogenesis of mTLE and other epilepsies. Although manipulation of miRNAs at acute stages has been reported to reduce subsequent spontaneous seizures, it is uncertain whether targeting miRNAs at chronic stages of mTLE can also reduce seizures. Furthermore, the functional role and downstream targets of most epilepsy-associated miRNAs remain poorly understood. Here, we show that miR-135a is selectively upregulated within neurons in epileptic brain and report that targeting miR-135a in vivo using antagomirs after onset of spontaneous recurrent seizures can reduce seizure activity at the chronic stage of experimental mTLE in male mice. Further, by using an unbiased approach combining immunoprecipitation and RNA sequencing, we identify several novel neuronal targets of miR-135a, including Mef2a Mef2 proteins are key regulators of excitatory synapse density. Mef2a and miR-135a show reciprocal expression regulation in human (of both sexes) and experimental TLE, and miR-135a regulates dendritic spine number and type through Mef2. Together, our data show that miR-135a is target for reducing seizure activity in chronic epilepsy, and that deregulation of miR-135a in epilepsy may alter Mef2a expression and thereby affect synaptic function and plasticity.SIGNIFICANCE STATEMENT miRNAs are post-transcriptional regulators of gene expression with roles in the pathogenesis of epilepsy. However, the precise mechanism of action and therapeutic potential of most epilepsy-associated miRNAs remain poorly understood. Our study reveals dramatic upregulation of the key neuronal miRNA miR-135a in both experimental and human mesial temporal lobe epilepsy. Silencing miR-135a in experimental temporal lobe epilepsy reduces seizure activity at the spontaneous recurrent seizure stage. These data support the exciting possibility that miRNAs can be targeted to combat seizures after spontaneous seizure activity has been established. Further, by using unbiased approaches novel neuronal targets of miR-135a, including members of the Mef2 protein family, are identified that begin to explain how deregulation of miR-135a may contribute to epilepsy