RESUMEN
The EMDataResource Ligand Model Challenge aimed to assess the reliability and reproducibility of modeling ligands bound to protein and protein-nucleic acid complexes in cryogenic electron microscopy (cryo-EM) maps determined at near-atomic (1.9-2.5 Å) resolution. Three published maps were selected as targets: Escherichia coli beta-galactosidase with inhibitor, SARS-CoV-2 virus RNA-dependent RNA polymerase with covalently bound nucleotide analog and SARS-CoV-2 virus ion channel ORF3a with bound lipid. Sixty-one models were submitted from 17 independent research groups, each with supporting workflow details. The quality of submitted ligand models and surrounding atoms were analyzed by visual inspection and quantification of local map quality, model-to-map fit, geometry, energetics and contact scores. A composite rather than a single score was needed to assess macromolecule+ligand model quality. These observations lead us to recommend best practices for assessing cryo-EM structures of liganded macromolecules reported at near-atomic resolution.
Asunto(s)
Microscopía por Crioelectrón , Modelos Moleculares , Microscopía por Crioelectrón/métodos , Ligandos , SARS-CoV-2 , COVID-19/virología , Escherichia coli , beta-Galactosidasa/química , beta-Galactosidasa/metabolismo , Conformación Proteica , Reproducibilidad de los ResultadosRESUMEN
Cyclooxygenase-2 (COX-2) is an enzyme that plays a pivotal role in peripheral inflammation and pain via the prostaglandin pathway. In the central nervous system (CNS), COX-2 is implicated in neurodegenerative and psychiatric disorders as a potential therapeutic target and biomarker. However, clinical studies with COX-2 have yielded inconsistent results, partly due to limited mechanistic understanding of how COX-2 activity relates to CNS pathology. Therefore, developing COX-2 positron emission tomography (PET) radiotracers for human neuroimaging is of interest. This study introduces [11C]BRD1158, which is a potent and uniquely fast-binding, selective COX-2 PET radiotracer. [11C]BRD1158 was developed by prioritizing potency at COX-2, isoform selectivity over COX-1, fast binding kinetics, and free fraction in the brain. Evaluated through in vivo PET neuroimaging in rodent models with human COX-2 overexpression, [11C]BRD1158 demonstrated high brain uptake, fast target-engagement, functional reversibility, and excellent specific binding, which is advantageous for human imaging applications. Lastly, post-mortem samples from Huntington's disease (HD) patients and preclinical HD mouse models showed that COX-2 levels were elevated specifically in disease-affected brain regions, primarily from increased expression in microglia. These findings indicate that COX-2 holds promise as a novel clinical marker of HD onset and progression, one of many potential applications of [11C]BRD1158 human PET.
RESUMEN
The EMDataResource Ligand Model Challenge aimed to assess the reliability and reproducibility of modeling ligands bound to protein and protein/nucleic-acid complexes in cryogenic electron microscopy (cryo-EM) maps determined at near-atomic (1.9-2.5 Å) resolution. Three published maps were selected as targets: E. coli beta-galactosidase with inhibitor, SARS-CoV-2 RNA-dependent RNA polymerase with covalently bound nucleotide analog, and SARS-CoV-2 ion channel ORF3a with bound lipid. Sixty-one models were submitted from 17 independent research groups, each with supporting workflow details. We found that (1) the quality of submitted ligand models and surrounding atoms varied, as judged by visual inspection and quantification of local map quality, model-to-map fit, geometry, energetics, and contact scores, and (2) a composite rather than a single score was needed to assess macromolecule+ligand model quality. These observations lead us to recommend best practices for assessing cryo-EM structures of liganded macromolecules reported at near-atomic resolution.
RESUMEN
Human genetic studies have revealed rare missense and protein-truncating variants in GRIN2A, encoding for the GluN2A subunit of the NMDA receptors, that confer significant risk for schizophrenia (SCZ). Mutations in GRIN2A are also associated with epilepsy and developmental delay/intellectual disability (DD/ID). However, it remains enigmatic how alterations to the same protein can result in diverse clinical phenotypes. Here, we performed functional characterization of human GluN1/GluN2A heteromeric NMDA receptors that contain SCZ-linked GluN2A variants, and compared them to NMDA receptors with GluN2A variants associated with epilepsy or DD/ID. Our findings demonstrate that SCZ-associated GRIN2A variants were predominantly loss-of-function (LoF), whereas epilepsy and DD/ID-associated variants resulted in both gain- and loss-of-function phenotypes. We additionally show that M653I and S809R, LoF GRIN2A variants associated with DD/ID, exert a dominant-negative effect when co-expressed with a wild-type GluN2A, whereas E58Ter and Y698C, SCZ-linked LoF variants, and A727T, an epilepsy-linked LoF variant, do not. These data offer a potential mechanism by which SCZ/epilepsy and DD/ID-linked variants can cause different effects on receptor function and therefore result in divergent pathological outcomes.
Asunto(s)
Epilepsia , Trastornos del Neurodesarrollo , Esquizofrenia , Humanos , Epilepsia/genética , Mutación , Trastornos del Neurodesarrollo/genética , Receptores de N-Metil-D-Aspartato/genética , Receptores de N-Metil-D-Aspartato/metabolismo , Esquizofrenia/genéticaRESUMEN
Recent advances in AI-based methods have revolutionized the field of structural biology. Concomitantly, high-throughput sequencing and functional genomics technologies have enabled the detection and generation of variants at an unprecedented scale. However, efficient tools and resources are needed to link these two disparate data types - to "map" variants onto protein structures, to better understand how the variation causes disease and thereby design therapeutics. Here we present the Genomics 2 Proteins Portal (G2P; g2p.broadinstitute.org/): a human proteome-wide resource that maps 19,996,443 genetic variants onto 42,413 protein sequences and 77,923 structures, with a comprehensive set of structural and functional features. Additionally, the G2P portal generalizes the capability of linking genomics to proteins beyond databases by allowing users to interactively upload protein residue-wise annotations (variants, scores, etc.) as well as the protein structure to establish the connection. The portal serves as an easy-to-use discovery tool for researchers and scientists to hypothesize the structure-function relationship between natural or synthetic variations and their molecular phenotype.
RESUMEN
Neurodevelopmental disorders (NDDs), including severe paediatric epilepsy, autism and intellectual disabilities are heterogeneous conditions in which clinical genetic testing can often identify a pathogenic variant. For many of them, genetic therapies will be tested in this or the coming years in clinical trials. In contrast to first-generation symptomatic treatments, the new disease-modifying precision medicines require a genetic test-informed diagnosis before a patient can be enrolled in a clinical trial. However, even in 2022, most identified genetic variants in NDD genes are 'variants of uncertain significance'. To safely enrol patients in precision medicine clinical trials, it is important to increase our knowledge about which regions in NDD-associated proteins can 'tolerate' missense variants and which ones are 'essential' and will cause a NDD when mutated. In addition, knowledge about functionally indispensable regions in the 3D structure context of proteins can also provide insights into the molecular mechanisms of disease variants. We developed a novel consensus approach that overlays evolutionary, and population based genomic scores to identify 3D essential sites (Essential3D) on protein structures. After extensive benchmarking of AlphaFold predicted and experimentally solved protein structures, we generated the currently largest expert curated protein structure set for 242 NDDs and identified 14 377 Essential3D sites across 189 gene disorders associated proteins. We demonstrate that the consensus annotation of Essential3D sites improves prioritization of disease mutations over single annotations. The identified Essential3D sites were enriched for functional features such as intermembrane regions or active sites and discovered key inter-molecule interactions in protein complexes that were otherwise not annotated. Using the currently largest autism, developmental disorders, and epilepsies exome sequencing studies including >360 000 NDD patients and population controls, we found that missense variants at Essential3D sites are 8-fold enriched in patients. In summary, we developed a comprehensive protein structure set for 242 NDDs and identified 14 377 Essential3D sites in these. All data are available at https://es-ndd.broadinstitute.org for interactive visual inspection to enhance variant interpretation and development of mechanistic hypotheses for 242 NDDs genes. The provided resources will enhance clinical variant interpretation and in silico drug target development for NDD-associated genes and encoded proteins.
Asunto(s)
Discapacidad Intelectual , Trastornos del Neurodesarrollo , Humanos , Niño , Trastornos del Neurodesarrollo/genética , Pruebas Genéticas , Mutación/genética , Discapacidad Intelectual/genética , Mutación MissenseRESUMEN
PPM1D encodes a serine/threonine phosphatase that regulates numerous pathways including the DNA damage response and p53. Activating mutations and amplification of PPM1D are found across numerous cancer types. GSK2830371 is a potent and selective allosteric inhibitor of PPM1D, but its mechanism of binding and inhibition of catalytic activity are unknown. Here we use computational, biochemical and functional genetic studies to elucidate the molecular basis of GSK2830371 activity. These data confirm that GSK2830371 binds an allosteric site of PPM1D with high affinity. By further incorporating data from hydrogen deuterium exchange mass spectrometry and sedimentation velocity analytical ultracentrifugation, we demonstrate that PPM1D exists in an equilibrium between two conformations that are defined by the movement of the flap domain, which is required for substrate recognition. A hinge region was identified that is critical for switching between the two conformations and was directly implicated in the high-affinity binding of GSK2830371 to PPM1D. We propose that the two conformations represent active and inactive forms of the protein reflected by the position of the flap, and that binding of GSK2830371 shifts the equilibrium to the inactive form. Finally, we found that C-terminal truncating mutations proximal to residue 400 result in destabilization of the protein via loss of a stabilizing N- and C-terminal interaction, consistent with the observation from human genetic data that nearly all PPM1D mutations in cancer are truncating and occur distal to residue 400. Taken together, our findings elucidate the mechanism by which binding of a small molecule to an allosteric site of PPM1D inhibits its activity and provides insights into the biology of PPM1D.
Asunto(s)
Neoplasias , Proteína Fosfatasa 2C , Sitio Alostérico , Aminopiridinas/farmacología , Dipéptidos/farmacología , Humanos , Mutación , Neoplasias/tratamiento farmacológico , Neoplasias/enzimología , Neoplasias/genética , Conformación Proteica , Proteína Fosfatasa 2C/antagonistas & inhibidores , Proteína Fosfatasa 2C/química , Proteína Fosfatasa 2C/genética , Proteína Fosfatasa 2C/metabolismo , Serina/genética , Serina/metabolismo , Relación Estructura-ActividadRESUMEN
All proteomes contain both proteins and polypeptide segments that don't form a defined three-dimensional structure yet are biologically active-called intrinsically disordered proteins and regions (IDPs and IDRs). Most of these IDPs/IDRs lack useful functional annotation limiting our understanding of their importance for organism fitness. Here we characterized IDRs using protein sequence annotations of functional sites and regions available in the UniProt knowledgebase ("UniProt features": active site, ligand-binding pocket, regions mediating protein-protein interactions, etc.). By measuring the statistical enrichment of twenty-five UniProt features in 981 IDRs of 561 human proteins, we identified eight features that are commonly located in IDRs. We then collected the genetic variant data from the general population and patient-based databases and evaluated the prevalence of population and pathogenic variations in IDPs/IDRs. We observed that some IDRs tolerate 2 to 12-times more single amino acid-substituting missense mutations than synonymous changes in the general population. However, we also found that 37% of all germline pathogenic mutations are located in disordered regions of 96 proteins. Based on the observed-to-expected frequency of mutations, we categorized 34 IDRs in 20 proteins (DDX3X, KIT, RB1, etc.) as intolerant to mutation. Finally, using statistical analysis and a machine learning approach, we demonstrate that mutation-intolerant IDRs carry a distinct signature of functional features. Our study presents a novel approach to assign functional importance to IDRs by leveraging the wealth of available genetic data, which will aid in a deeper understating of the role of IDRs in biological processes and disease mechanisms.
Asunto(s)
Proteínas Intrínsecamente Desordenadas , Secuencia de Aminoácidos , Variación Genética/genética , Humanos , Proteínas Intrínsecamente Desordenadas/química , Conformación Proteica , Proteoma/genéticaRESUMEN
PURPOSE: Pathogenic variants in GABRB3 have been associated with a spectrum of phenotypes from severe developmental disorders and epileptic encephalopathies to milder epilepsy syndromes and mild intellectual disability (ID). In this study, we analyzed a large cohort of individuals with GABRB3 variants to deepen the phenotypic understanding and investigate genotype-phenotype correlations. METHODS: Through an international collaboration, we analyzed electro-clinical data of unpublished individuals with variants in GABRB3, and we reviewed previously published cases. All missense variants were mapped onto the 3-dimensional structure of the GABRB3 subunit, and clinical phenotypes associated with the different key structural domains were investigated. RESULTS: We characterized 71 individuals with GABRB3 variants, including 22 novel subjects, expressing a wide spectrum of phenotypes. Interestingly, phenotypes correlated with structural locations of the variants. Generalized epilepsy, with a median age at onset of 12 months, and mild-to-moderate ID were associated with variants in the extracellular domain. Focal epilepsy with earlier onset (median: age 4 months) and severe ID were associated with variants in both the pore-lining helical transmembrane domain and the extracellular domain. CONCLUSION: These genotype-phenotype correlations will aid the genetic counseling and treatment of individuals affected by GABRB3-related disorders. Future studies may reveal whether functional differences underlie the phenotypic differences.
Asunto(s)
Epilepsia , Discapacidad Intelectual , Epilepsia/genética , Estudios de Asociación Genética , Humanos , Discapacidad Intelectual/genética , Mutación , Fenotipo , Receptores de GABA-A/genéticaRESUMEN
PRMT5 and its substrate adaptor proteins (SAPs), pICln and Riok1, are synthetic lethal dependencies in MTAP-deleted cancer cells. SAPs share a conserved PRMT5 binding motif (PBM) which mediates binding to a surface of PRMT5 distal to the catalytic site. This interaction is required for methylation of several PRMT5 substrates, including histone and spliceosome complexes. We screened for small molecule inhibitors of the PRMT5-PBM interaction and validated a compound series which binds to the PRMT5-PBM interface and directly inhibits binding of SAPs. Mode of action studies revealed the formation of a covalent bond between a halogenated pyridazinone group and cysteine 278 of PRMT5. Optimization of the starting hit produced a lead compound, BRD0639, which engages the target in cells, disrupts PRMT5-RIOK1 complexes, and reduces substrate methylation. BRD0639 is a first-in-class PBM-competitive inhibitor that can support studies of PBM-dependent PRMT5 activities and the development of novel PRMT5 inhibitors that selectively target these functions.
Asunto(s)
Proteínas Adaptadoras Transductoras de Señales/antagonistas & inhibidores , Descubrimiento de Drogas , Proteína-Arginina N-Metiltransferasas/antagonistas & inhibidores , Piridazinas/farmacología , Proteínas Adaptadoras Transductoras de Señales/metabolismo , Relación Dosis-Respuesta a Droga , Humanos , Modelos Moleculares , Estructura Molecular , Proteína-Arginina N-Metiltransferasas/metabolismo , Piridazinas/síntesis química , Piridazinas/química , Relación Estructura-ActividadRESUMEN
Background: Oculocutaneous albinism (OCA) is a Mendelian disorder characterized by hypopigmentation of the skin, hair, and eyes, hypoplastic fovea, and low vision, known to be caused by mutations in the Tyrosinase (TYR) gene. Among the known TYR variants, some reduce but do not completely eliminate tyrosinase activity, allowing residual production of melanin and resulting in a contradictory assignment as either pathogenic or benign, preventing a precise clinical diagnostic.Materials and Methods: In the present work, we performed Whole Exome Sequencing and subsequent Sanger sequencing in a young male clinically diagnosed with OCA.Results: Whole-exome sequencing analysis revealed the identification of two variants in trans in TYR. The first, corresponds to a known pathogenic variant G47D, while the second S192Y, was considered a polymorphism due to its relatively high frequency in the European population.Conclusion: The lack of other pathogenic variants in TYR, the reported reduced enzymatic activity (ca. 40% respect to wt) for S192Y, together with the structural in-silico analysis strongly suggest that both reported variants are jointly disease-causing and that S192Y should be considered as likely pathogenic, especially when it is found in trans with a null variant.
Asunto(s)
Albinismo Oculocutáneo/genética , Monofenol Monooxigenasa/genética , Mutación Missense/genética , Polimorfismo de Nucleótido Simple/genética , Adolescente , Albinismo Oculocutáneo/diagnóstico , Secuencia de Aminoácidos , Estudios de Asociación Genética , Predisposición Genética a la Enfermedad , Genotipo , Humanos , Masculino , Datos de Secuencia Molecular , Linaje , Secuenciación del ExomaRESUMEN
Genetic variation of the 16p11.2 deletion locus containing the KCTD13 gene and of CUL3 is linked with autism. This genetic connection suggested that substrates of a CUL3-KCTD13 ubiquitin ligase may be involved in disease pathogenesis. Comparison of Kctd13 mutant (Kctd13 -/- ) and wild-type neuronal ubiquitylomes identified adenylosuccinate synthetase (ADSS), an enzyme that catalyzes the first step in adenosine monophosphate (AMP) synthesis, as a KCTD13 ligase substrate. In Kctd13 -/- neurons, there were increased levels of succinyl-adenosine (S-Ado), a metabolite downstream of ADSS. Notably, S-Ado levels are elevated in adenylosuccinate lyase deficiency, a metabolic disorder with autism and epilepsy phenotypes. The increased S-Ado levels in Kctd13 -/- neurons were decreased by treatment with an ADSS inhibitor. Lastly, functional analysis of human KCTD13 variants suggests that KCTD13 variation may alter ubiquitination of ADSS. These data suggest that succinyl-AMP metabolites accumulate in Kctd13 -/- neurons, and this observation may have implications for our understanding of 16p11.2 deletion syndrome.
RESUMEN
Advances in gene discovery have identified genetic variants in the solute carrier family 6 member 1 gene as a monogenic cause of neurodevelopmental disorders, including epilepsy with myoclonic atonic seizures, autism spectrum disorder and intellectual disability. The solute carrier family 6 member 1 gene encodes for the GABA transporter protein type 1, which is responsible for the reuptake of the neurotransmitter GABA, the primary inhibitory neurotransmitter in the central nervous system, from the extracellular space. GABAergic inhibition is essential to counterbalance neuronal excitation, and when significantly disrupted, it negatively impacts brain development leading to developmental differences and seizures. Aggregation of patient variants and observed clinical manifestations expand understanding of the genotypic and phenotypic spectrum of this disorder. Here, we assess genetic and phenotypic features in 116 individuals with solute carrier family 6 member 1 variants, the vast majority of which are likely to lead to GABA transporter protein type 1 loss-of-function. The knowledge acquired will guide therapeutic decisions and the development of targeted therapies that selectively enhance transporter function and may improve symptoms. We analysed the longitudinal and cell type-specific expression of solute carrier family 6 member 1 in humans and localization of patient and control missense variants in a novel GABA transporter protein type 1 protein structure model. In this update, we discuss the progress made in understanding and treating solute carrier family 6 member 1-related disorders thus far, through the concerted efforts of clinicians, scientists and family support groups.
RESUMEN
Interpretation of the colossal number of genetic variants identified from sequencing applications is one of the major bottlenecks in clinical genetics, with the inference of the effect of amino acid-substituting missense variations on protein structure and function being especially challenging. Here we characterize the three-dimensional (3D) amino acid positions affected in pathogenic and population variants from 1,330 disease-associated genes using over 14,000 experimentally solved human protein structures. By measuring the statistical burden of variations (i.e., point mutations) from all genes on 40 3D protein features, accounting for the structural, chemical, and functional context of the variations' positions, we identify features that are generally associated with pathogenic and population missense variants. We then perform the same amino acid-level analysis individually for 24 protein functional classes, which reveals unique characteristics of the positions of the altered amino acids: We observe up to 46% divergence of the class-specific features from the general characteristics obtained by the analysis on all genes, which is consistent with the structural diversity of essential regions across different protein classes. We demonstrate that the function-specific 3D features of the variants match the readouts of mutagenesis experiments for BRCA1 and PTEN, and positively correlate with an independent set of clinically interpreted pathogenic and benign missense variants. Finally, we make our results available through a web server to foster accessibility and downstream research. Our findings represent a crucial step toward translational genetics, from highlighting the impact of mutations on protein structure to rationalizing the variants' pathogenicity in terms of the perturbed molecular mechanisms.
Asunto(s)
Mutación Missense/genética , Proteínas/química , Proteínas/genética , Secuencia de Aminoácidos , Proteína BRCA1/química , Proteína BRCA1/genética , Biología Computacional/métodos , Humanos , Aprendizaje Automático , Modelos Moleculares , Mutación Missense/fisiología , Fosfohidrolasa PTEN/química , Fosfohidrolasa PTEN/genética , Conformación Proteica , Proteínas/fisiologíaRESUMEN
Malfunctions of voltage-gated sodium and calcium channels (encoded by SCNxA and CACNA1x family genes, respectively) have been associated with severe neurologic, psychiatric, cardiac, and other diseases. Altered channel activity is frequently grouped into gain or loss of ion channel function (GOF or LOF, respectively) that often corresponds not only to clinical disease manifestations but also to differences in drug response. Experimental studies of channel function are therefore important, but laborious and usually focus only on a few variants at a time. On the basis of known gene-disease mechanisms of 19 different diseases, we inferred LOF (n = 518) and GOF (n = 309) likely pathogenic variants from the disease phenotypes of variant carriers. By training a machine learning model on sequence- and structure-based features, we predicted LOF or GOF effects [area under the receiver operating characteristics curve (ROC) = 0.85] of likely pathogenic missense variants. Our LOF versus GOF prediction corresponded to molecular LOF versus GOF effects for 87 functionally tested variants in SCN1/2/8A and CACNA1I (ROC = 0.73) and was validated in exome-wide data from 21,703 cases and 128,957 controls. We showed respective regional clustering of inferred LOF and GOF nucleotide variants across the alignment of the entire gene family, suggesting shared pathomechanisms in the SCNxA/CACNA1x family genes.
Asunto(s)
Canales de Calcio , Preparaciones Farmacéuticas , Mutación Missense/genética , Fenotipo , SodioRESUMEN
Prion disease is a rapidly progressive neurodegenerative disorder caused by misfolding and aggregation of the prion protein (PrP), and there are currently no therapeutic options. PrP ligands could theoretically antagonize prion formation by protecting the native protein from misfolding or by targeting it for degradation, but no validated small-molecule binders have been discovered to date. We deployed a variety of screening methods in an effort to discover binders of PrP, including 19F-observed and saturation transfer difference (STD) NMR spectroscopy, differential scanning fluorimetry (DSF), DNA-encoded library selection, and in silico screening. A single benzimidazole compound was confirmed in concentration-response, but affinity was very weak (Kd > 1 mm), and it could not be advanced further. The exceptionally low hit rate observed here suggests that PrP is a difficult target for small-molecule binders. Whereas orthogonal binder discovery methods could yield high-affinity compounds, non-small-molecule modalities may offer independent paths forward against prion disease.
Asunto(s)
Bencimidazoles/farmacología , Enfermedades por Prión/tratamiento farmacológico , Proteínas Priónicas/antagonistas & inhibidores , Bibliotecas de Moléculas Pequeñas/farmacología , Bencimidazoles/química , Descubrimiento de Drogas , Evaluación Preclínica de Medicamentos , Humanos , Espectroscopía de Resonancia Magnética , Enfermedades por Prión/metabolismo , Proteínas Priónicas/metabolismo , Bibliotecas de Moléculas Pequeñas/químicaRESUMEN
Human genome sequencing efforts have greatly expanded, and a plethora of missense variants identified both in patients and in the general population is now publicly accessible. Interpretation of the molecular-level effect of missense variants, however, remains challenging and requires a particular investigation of amino acid substitutions in the context of protein structure and function. Answers to questions like 'Is a variant perturbing a site involved in key macromolecular interactions and/or cellular signaling?', or 'Is a variant changing an amino acid located at the protein core or part of a cluster of known pathogenic mutations in 3D?' are crucial. Motivated by these needs, we developed MISCAST (missense variant to protein structure analysis web suite; http://miscast.broadinstitute.org/). MISCAST is an interactive and user-friendly web server to visualize and analyze missense variants in protein sequence and structure space. Additionally, a comprehensive set of protein structural and functional features have been aggregated in MISCAST from multiple databases, and displayed on structures alongside the variants to provide users with the biological context of the variant location in an integrated platform. We further made the annotated data and protein structures readily downloadable from MISCAST to foster advanced offline analysis of missense variants by a wide biological community.
Asunto(s)
Mutación Missense , Conformación Proteica , Programas Informáticos , Humanos , Internet , Proteínas/química , Proteínas/genéticaRESUMEN
OBJECTIVE: Voltage-gated sodium channels (SCNs) share similar amino acid sequence, structure, and function. Genetic variants in the four human brain-expressed SCN genes SCN1A/2A/3A/8A have been associated with heterogeneous epilepsy phenotypes and neurodevelopmental disorders. To better understand the biology of seizure susceptibility in SCN-related epilepsies, our aim was to determine similarities and differences between sodium channel disorders, allowing us to develop a broader perspective on precision treatment than on an individual gene level alone. METHODS: We analyzed genotype-phenotype correlations in large SCN-patient cohorts and applied variant constraint analysis to identify severe sodium channel disease. We examined temporal patterns of human SCN expression and correlated functional data from in vitro studies with clinical phenotypes across different sodium channel disorders. RESULTS: Comparing 865 epilepsy patients (504 SCN1A, 140 SCN2A, 171 SCN8A, four SCN3A, 46 copy number variation [CNV] cases) and analysis of 114 functional studies allowed us to identify common patterns of presentation. All four epilepsy-associated SCN genes demonstrated significant constraint in both protein truncating and missense variation when compared to other SCN genes. We observed that age at seizure onset is related to SCN gene expression over time. Individuals with gain-of-function SCN2A/3A/8A missense variants or CNV duplications share similar characteristics, most frequently present with early onset epilepsy (<3 months), and demonstrate good response to sodium channel blockers (SCBs). Direct comparison of corresponding SCN variants across different SCN subtypes illustrates that the functional effects of variants in corresponding channel locations are similar; however, their clinical manifestation differs, depending on their role in different types of neurons in which they are expressed. SIGNIFICANCE: Variant function and location within one channel can serve as a surrogate for variant effects across related sodium channels. Taking a broader view on precision treatment suggests that in those patients with a suspected underlying genetic epilepsy presenting with neonatal or early onset seizures (<3 months), SCBs should be considered.
Asunto(s)
Síndromes Epilépticos/genética , Canal de Sodio Activado por Voltaje NAV1.1/genética , Canal de Sodio Activado por Voltaje NAV1.2/genética , Canal de Sodio Activado por Voltaje NAV1.3/genética , Canal de Sodio Activado por Voltaje NAV1.6/genética , Canales de Sodio/genética , Edad de Inicio , Trastorno del Espectro Autista/genética , Trastorno del Espectro Autista/fisiopatología , Niño , Preescolar , Codón sin Sentido , Variaciones en el Número de Copia de ADN , Electroencefalografía , Síndromes Epilépticos/tratamiento farmacológico , Síndromes Epilépticos/fisiopatología , Femenino , Mutación con Ganancia de Función , Eliminación de Gen , Duplicación de Gen , Expresión Génica , Regulación del Desarrollo de la Expresión Génica , Genotipo , Humanos , Lactante , Recién Nacido , Mutación con Pérdida de Función , Masculino , Mutación Missense , Canal de Sodio Activado por Voltaje NAV1.1/metabolismo , Canal de Sodio Activado por Voltaje NAV1.2/metabolismo , Canal de Sodio Activado por Voltaje NAV1.3/metabolismo , Canal de Sodio Activado por Voltaje NAV1.6/metabolismo , Trastornos del Neurodesarrollo/genética , Trastornos del Neurodesarrollo/fisiopatología , Fenotipo , Bloqueadores de los Canales de Sodio/uso terapéutico , Canales de Sodio/metabolismoRESUMEN
The four voltage-gated sodium channels SCN1/2/3/8A have been associated with heterogeneous types of developmental disorders, each presenting with disease specific temporal and cell type specific gene expression. Using single-cell RNA sequencing transcriptomic data from humans and mice, we observe that SCN1A is predominantly expressed in inhibitory neurons. In contrast, SCN2/3/8A are profoundly expressed in excitatory neurons with SCN2/3A starting prenatally, followed by SCN1/8A neonatally. In contrast to previous observations from low resolution RNA screens, we observe that all four genes are expressed in both excitatory and inhibitory neurons, however, exhibit differential expression strength. These findings provide molecular evidence, at single-cell resolution, to support the hypothesis that the excitatory/inhibitory (E/I) neuronal expression ratios of sodium channels are important regulatory mechanisms underlying brain homeostasis and neurological diseases. Modulating the E/I expression balance within cell types of sodium channels could serve as a potential strategy to develop targeted treatment for NaV-associated neuronal developmental disorders.
Asunto(s)
Encéfalo/metabolismo , Discapacidades del Desarrollo/metabolismo , Neuronas/metabolismo , Canales de Sodio Activados por Voltaje/metabolismo , Animales , Discapacidades del Desarrollo/genética , Humanos , Ratones , Canales de Sodio Activados por Voltaje/genéticaRESUMEN
The human RNA helicase DDX6 is an essential component of membrane-less organelles called processing bodies (PBs). PBs are involved in mRNA metabolic processes including translational repression via coordinated storage of mRNAs. Previous studies in human cell lines have implicated altered DDX6 in molecular and cellular dysfunction, but clinical consequences and pathogenesis in humans have yet to be described. Here, we report the identification of five rare de novo missense variants in DDX6 in probands presenting with intellectual disability, developmental delay, and similar dysmorphic features including telecanthus, epicanthus, arched eyebrows, and low-set ears. All five missense variants (p.His372Arg, p.Arg373Gln, p.Cys390Arg, p.Thr391Ile, and p.Thr391Pro) are located in two conserved motifs of the RecA-2 domain of DDX6 involved in RNA binding, helicase activity, and protein-partner binding. We use functional studies to demonstrate that the first variants identified (p.Arg373Gln and p.Cys390Arg) cause significant defects in PB assembly in primary fibroblast and model human cell lines. These variants' interactions with several protein partners were also disrupted in immunoprecipitation assays. Further investigation via complementation assays included the additional variants p.Thr391Ile and p.Thr391Pro, both of which, similarly to p.Arg373Gln and p.Cys390Arg, demonstrated significant defects in P-body assembly. Complementing these molecular findings, modeling of the variants on solved protein structures showed distinct spatial clustering near known protein binding regions. Collectively, our clinical and molecular data describe a neurodevelopmental syndrome associated with pathogenic missense variants in DDX6. Additionally, we suggest DDX6 join the DExD/H-box genes DDX3X and DHX30 in an emerging class of neurodevelopmental disorders involving RNA helicases.