Delineation of functionally essential protein regions for 242 neurodevelopmental genes

in Brain (2022), awac381

Neurodevelopmental disorders (NDDs), including severe pediatric epilepsy, autism, and intellectual disabilities are heterogeneous conditions in which clinical genetic testing can often identify a ... [more ▼]

Neurodevelopmental disorders (NDDs), including severe pediatric 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 enroll 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 three-dimensional (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 neurodevelopmental disorders 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. [less ▲]

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

E-print/Working paper (2022)

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 ... [more ▼]

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 Format [less ▲]

Comprehensive characterization of amino acidpositions in protein structures reveals moleculareffect of missense variants

in Proceedings of the National Academy of Sciences of the United States of America (2020)

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 ... [more ▼]

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. [less ▲]

MISCAST: MIssense variant to protein StruCture Analysis web SuiTe

in Nucleic Acids Research (2020)

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 ... [more ▼]

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. [less ▲]

Assessment of genetic variant burden in epilepsy-associated brain lesions

in European Journal of Human Genetics (2019)

It is challenging to estimate genetic variant burden across different subtypes of epilepsy. Herein, we used a comparative approach to assess the diagnostic yield and genotype-phenotype correlations in the ... [more ▼]

It is challenging to estimate genetic variant burden across different subtypes of epilepsy. Herein, we used a comparative approach to assess the diagnostic yield and genotype-phenotype correlations in the four most common brain lesions in patients with drug-resistant focal epilepsy. Targeted sequencing analysis was performed for a panel of 161 genes with a mean coverage of > 400x. Lesional tissue was histopathologically reviewed and dissected from hippocampal sclerosis (n=15), ganglioglioma (n=16), dysembryoplastic neuroepithelial tumors (n=8) and ocal cortical dysplasia type II (n=15). Peripheral blood (n=12) or surgical tissue samples histopathologically classified as lesion-free (n=42) were available for comparison. Variants were classified as pathogenic or likely pathogenic according to American College of Medical Genetics and Genomics guidelines. Overall, we identified pathogenic and likely pathogenic variants in 25.9% of patients with a mean coverage of 383x. The highest number of pathogenic/ likely pathogenic variants was observed in patients with ganglioglioma (43.75%; all somatic) and dysembryoplastic neuroepithelial tumors (37.5%; all somatic), and in 20% of cases with focal cortical dysplasia type II (13.33% somatic, 6.67% germline). Pathogenic/likely pathogenic positive genes were disorder-specific and BRAF V600E the only recurrent pathogenic variant. This study represents a reference for diagnostic yield across the four most common lesion entities in patients with drug-resistant focal epilepsy. The observed large variability in variant burden by epileptic lesion type calls for whole exome sequencing of histopathologically well characterized tissue in a diagnostic setting and in research to discover novel disease-associated genes. [less ▲]

Simple ClinVar:an interactive web server to explore and retrieve gene and disease variants 1 aggregated in ClinVar database

in Nucleic Acids Research (2019)

Clinical genetic testing has exponentially expanded in recent years, leading to an overwhelming amount of patient variants with high variability in pathogenicity and heterogeneous phenotypes. A large part ... [more ▼]

Clinical genetic testing has exponentially expanded in recent years, leading to an overwhelming amount of patient variants with high variability in pathogenicity and heterogeneous phenotypes. A large part of the variant level data are comprehensively aggregated in public databases such as ClinVar. However, the ability to explore this rich resource and answer general questions such as “How many genes inside ClinVar are associated with a specific disease? or “In which part of the protein are patient variants located?” is limited and requires advanced bioinformatics processing. Here, we present Simple ClinVar (http://simple-clinvar.broadinstitute.org/) a web- server application that is able to provide variant, gene, and disease level summary statistics based on the entire ClinVar database in a dynamic and user-friendly web-interface. Overall, our web application is able to interactively answer basic questions regarding genetic variation and its known relationships to disease. By typing a disease term of interest, t he user can identify in seconds the genes and phenotypes most frequently reported to ClinVar. Subsets of variants can then be further explored, filtered, or mapped and visualized in the corresponding protein sequences. Our website will follow ClinVar monthly releases and provide easy access to rich ClinVar resources to a broader audience including basic and clinical scientists. [less ▲]

Variant Score Ranker - a web application for intuitive missense variant prioritization

in Bioinformatics (2019)

The correct classification of missense variants as benign or pathogenic remains challenging. Pathogenic variants are expected to have higher deleterious prediction scores than benign variants in the same ... [more ▼]

The correct classification of missense variants as benign or pathogenic remains challenging. Pathogenic variants are expected to have higher deleterious prediction scores than benign variants in the same gene. However, most of the existing variant annotation tools do not reference the score range of benign population variants on gene level. Here, we present a web-application, Variant Score Ranker, which enables users to rapidly annotate variants and perform gene-specific variant score ranking on the population level. We also provide an intuitive example of how gene- and population-calibrated variant ranking scores can improve epilepsy variant prioritization. [less ▲]

Spectrum of GABAA receptor variants in epilepsy

in Current Opinion in Neurology (2019)

Purpose of review: Recent publications point to an increasingly important role of variants in genes encoding GABAA receptor subunits associated with both common and rare forms of epilepsies. The aim of ... [more ▼]

Purpose of review: Recent publications point to an increasingly important role of variants in genes encoding GABAA receptor subunits associated with both common and rare forms of epilepsies. The aim of this review is to give an overview of the current clinical phenotypes, genetic findings and pathophysiological mechanisms related to GABAA receptor variants. Recent findings: Early work showed that inherited variants in GABRG2 and GABRA1 cause relatively mild forms of monogenic epilepsies in large families. More recent studies have revealed that de novo variants in several GABAA receptor genes cause severe developmental and epileptic encephalopathies, inherited variants cause remarkably variable phenotypes within the same pedigrees ranging from asymptomatic carriers to developmental and epileptic encephalopathies, and variants in all GABAA receptor genes are enriched in common forms of epilepsy, namely rolandic epilepsy and genetic generalized epilepsy. Analyses from cellular expression systems and mouse models suggest that all variants cause a loss of GABAA receptor function resulting in GABAergic disinhibition. Summary: Genetic studies have revealed a crucial role of the GABAergic system in the underlying pathogenesis of various forms of common and rare epilepsies. Our understanding of functional consequences of GABAA receptor variants provide an opportunity to develop precision-based therapeutic strategies that are hopefully free from the side-effect burden seen with currently available GABAergic drugs. [less ▲]

Identification and Characterization of Variant Intolerant Sites across Human Protein 3-Dimensional Structures

in Biophysical Journal (2018, February 02), 114(3, Suppl. 1), 664

The functional interpretation of genetic variation in disease-associated genes is far outpaced by data generation. Existing algorithms for prediction of variant consequences do not adequately distinguish ... [more ▼]

The functional interpretation of genetic variation in disease-associated genes is far outpaced by data generation. Existing algorithms for prediction of variant consequences do not adequately distinguish pathogenic variants from benign rare variants. This lack of statistical and bioinformatics analyses, accompanied by an ever-increasing number of identified variants in biomedical research and clinical applications, has become a major challenge. Established methods to predict the functional effect of genetic variation use the degree of amino acid conservation across species in linear protein sequence alignment. More recent methods include the spatial distribution pattern of known patient and control variants. Here, we propose to combine the linear conservation and spatial constrained based scores to devise a novel score that incorporates 3-dimensional structural properties of amino acid residues, such as the solvent-accessible surface area, degree of flexibility, secondary structure propensity and binding tendency, to quantify the effect of amino acid substitutions. For this study, we develop a framework for large-scale mapping of established linear sequence-based paralog and ortholog conservation scores onto the tertiary structures of human proteins. This framework can be utilized to map the spatial distribution of mutations on solved protein structures as well as homology models. As a proof of concept, using a homology model of the human Nav1.2 voltage-gated sodium channel structure, we observe spatial clustering in distinct domains of mutations, associated with Autism Spectrum Disorder (>20 variants) and Epilepsy (>100 variants), that exert opposing effects on channel function. We are currently characterizing all variants (>300k individuals) found in ClinVar, the largest disease variant database, as well as variants identified in >140k individuals from general population. The variant mapping framework and our score, informed with structural information, will be useful in identifying structural motifs of proteins associated with disease risk. [less ▲]