Three-dimensional missense tolerance ratio analysis

Distances from ligands as main predictive features for pathogenicity and functional effect of variants in NMDA receptors

Genetic variants in genes GRIN1 , GRIN2A , GRIN2B , and GRIN2D , which encode subunits of the N-methyl-D-aspartate receptor (NMDAR), have been associated with severe and heterogeneous neurologic diseases. Missense variants in these genes can result in gain or loss of the NMDAR function, requiring opposite therapeutic treatments. Computational methods that predict pathogenicity and molecular functional effects are therefore crucial for accurate diagnosis and therapeutic applications. We assembled missense variants: 201 from patients, 631 from general population, and 159 characterized by electrophysiological readouts showing whether they can enhance or reduce the receptor function. This includes new functional data from 47 variants reported here, for the first time. We found that pathogenic/benign variants and variants that increase/decrease the channel function were distributed unevenly on the protein structure, with spatial proximity to ligands bound to the agonist and antagonist binding sites being key predictive features. Leveraging distances from ligands, we developed two independent machine learning-based predictors for NMDAR missense variants: a pathogenicity predictor which outperforms currently available predictors (AUC=0.945, MCC=0.726), and the first binary predictor of molecular function (increase or decrease) (AUC=0.809, MCC=0.523). Using these, we reclassified variants of uncertain significance in the ClinVar database and refined a previous genome-informed epidemiological model to estimate the birth incidence of molecular mechanism-defined GRIN disorders. Our findings demonstrate that distance from ligands is an important feature in NMDARs that can enhance variant pathogenicity prediction and enable functional prediction. Further studies with larger numbers of phenotypically and functionally characterized variants will enhance the potential clinical utility of this method.

The landscape of regional missense mutational intolerance quantified from 125,748 exomes

Missense variants can have a range of functional impacts depending on factors such as the specific amino acid substitution and location within the gene. To interpret their deleteriousness, studies have sought to identify regions within genes that are specifically intolerant of missense variation 1-12 . Here, we leverage the patterns of rare missense variation in 125,748 individuals in the Genome Aggregation Database (gnomAD) 13 against a null mutational model to identify transcripts that display regional differences in missense constraint. Missense-depleted regions are enriched for ClinVar 14 pathogenic variants, de novo missense variants from individuals with neurodevelopmental disorders (NDDs) 15,16 , and complex trait heritability. Following ClinGen calibration recommendations for the ACMG/AMP guidelines, we establish that regions with less than 20% of their expected missense variation achieve moderate support for pathogenicity. We create a missense deleteriousness metric (MPC) that incorporates regional constraint and outperforms other deleteriousness scores at stratifying case and control de novo missense variation, with a strong enrichment in NDDs. These results provide additional tools to aid in missense variant interpretation.

A unified analysis of evolutionary and population constraint in protein domains highlights structural features and pathogenic sites

Protein evolution is constrained by structure and function, creating patterns in residue conservation that are routinely exploited to predict structure and other features. Similar constraints should affect variation across individuals, but it is only with the growth of human population sequencing that this has been tested at scale. Now, human population constraint has established applications in pathogenicity prediction, but it has not yet been explored for structural inference. Here, we map 2.4 million population variants to 5885 protein families and quantify residue-level constraint with a new Missense Enrichment Score (MES). Analysis of 61,214 structures from the PDB spanning 3661 families shows that missense depleted sites are enriched in buried residues or those involved in small-molecule or protein binding. MES is complementary to evolutionary conservation and a combined analysis allows a new classification of residues according to a conservation plane. This approach finds functional residues that are evolutionarily diverse, which can be related to specificity, as well as family-wide conserved sites that are critical for folding or function. We also find a possible contrast between lethal and non-lethal pathogenic sites, and a surprising clinical variant hot spot at a subset of missense enriched positions.

De novo GRIN variants in M3 helix associated with neurological disorders control channel gating of NMDA receptor

N-methyl-D-aspartate receptors (NMDARs) are members of the glutamate receptor family and participate in excitatory postsynaptic transmission throughout the central nervous system. Genetic variants in GRIN genes encoding NMDAR subunits are associated with a spectrum of neurological disorders. The M3 transmembrane helices of the NMDAR couple directly to the agonist-binding domains and form a helical bundle crossing in the closed receptors that occludes the pore. The M3 functions as a transduction element whose conformational change couples ligand binding to opening of an ion conducting pore. In this study, we report the functional consequences of 48 de novo missense variants in GRIN1, GRIN2A, and GRIN2B that alter residues in the M3 transmembrane helix. These de novo variants were identified in children with neurological and neuropsychiatric disorders including epilepsy, developmental delay, intellectual disability, hypotonia and attention deficit hyperactivity disorder. All 48 variants in M3 for which comprehensive testing was completed produce a gain-of-function (28/48) compared to loss-of-function (9/48); 11 variants had an indeterminant phenotype. This supports the idea that a key structural feature of the M3 gate exists to stabilize the closed state so that agonist binding can drive channel opening. Given that most M3 variants enhance channel gating, we assessed the potency of FDA-approved NMDAR channel blockers on these variant receptors. These data provide new insight into the structure-function relationship of the NMDAR gate, and suggest that variants within the M3 transmembrane helix produce a gain-of-function.

Clinical features, functional consequences, and rescue pharmacology of missense GRID1 and GRID2 human variants

GRID1 and GRID2 encode the enigmatic GluD1 and GluD2 proteins, which form tetrameric receptors that play important roles in synapse organization and development of the central nervous system. Variation in these genes has been implicated in neurodevelopmental phenotypes. We evaluated GRID1 and GRID2 human variants from the literature, ClinVar, and clinical laboratories and found that many of these variants reside in intolerant domains, including the amino terminal domain of both GRID1 and GRID2. Other conserved regions, such as the M3 transmembrane domain, show different intolerance between GRID1 and GRID2. We introduced these variants into GluD1 and GluD2 cDNA and performed electrophysiological and biochemical assays to investigate the mechanisms of dysfunction of GRID1/2 variants. One variant in the GRID1 distal amino terminal domain resides at a position predicted to interact with Cbln2/Cbln4, and the variant disrupts complex formation between GluD1 and Cbln2, which could perturb its role in synapse organization. We also discovered that, like the lurcher mutation (GluD2-A654T), other rare variants in the GRID2 M3 domain create constitutively active receptors that share similar pathogenic phenotypes. We also found that the SCHEMA schizophrenia M3 variant GluD1-A650T produced constitutively active receptors. We tested a variety of compounds for their ability to inhibit constitutive currents of GluD receptor variants and found that pentamidine potently inhibited GluD2-T649A constitutive channels (IC50 50 nM). These results identify regions of intolerance to variation in the GRID genes, illustrate the functional consequences of GRID1 and GRID2 variants, and suggest how these receptors function normally and in disease.

Therapeutic potential of N-methyl-D-aspartate receptor modulators in psychiatry

N-methyl-D-aspartate (NMDA) receptors mediate a slow component of excitatory synaptic transmission, are widely distributed throughout the central nervous system, and regulate synaptic plasticity. NMDA receptor modulators have long been considered as potential treatments for psychiatric disorders including depression and schizophrenia, neurodevelopmental disorders such as Rett Syndrome, and neurodegenerative conditions such as Alzheimer's disease. New interest in NMDA receptors as therapeutic targets has been spurred by the findings that certain inhibitors of NMDA receptors produce surprisingly rapid and robust antidepressant activity by a novel mechanism, the induction of changes in the brain that well outlast the presence of drug in the body. These findings are driving research into an entirely new paradigm for using NMDA receptor antagonists in a host of related conditions. At the same time positive allosteric modulators of NMDA receptors are being pursued for enhancing synaptic function in diseases that feature NMDA receptor hypofunction. While there is great promise, developing the therapeutic potential of NMDA receptor modulators must also navigate the potential significant risks posed by the use of such agents. We review here the emerging pharmacology of agents that target different NMDA receptor subtypes, offering new avenues for capturing the therapeutic potential of targeting this important receptor class.

Clinical and functional consequences of GRIA variants in patients with neurological diseases

AMPA receptors are members of the glutamate receptor family and mediate a fast component of excitatory synaptic transmission at virtually all central synapses. Thus, their functional characteristics are a critical determinant of brain function. We evaluate intolerance of each GRIA gene to genetic variation using 3DMTR and report here the functional consequences of 52 missense variants in GRIA1-4 identified in patients with various neurological disorders. These variants produce changes in agonist EC50, response time course, desensitization, and/or receptor surface expression. We predict that these functional and localization changes will have important consequences for circuit function, and therefore likely contribute to the patients' clinical phenotype. We evaluated the sensitivity of variant receptors to AMPAR-selective modulators including FDA-approved drugs to explore potential targeted therapeutic options.

De novo CLPTM1 variants with reduced GABAA R current response in patients with epilepsy

OBJECTIVE

To investigate the clinical features and potential pathogenesis mechanism of de novo CLPTM1 variants associated with epilepsy.

METHODS

Identify de novo genetic variants associated with epilepsy by reanalyzing trio-based whole-exome sequencing data. We analyzed the clinical characteristics of patients with these variants and performed functional in vitro studies in cells expressing mutant complementary DNA for these variants using whole-cell voltage-clamp current recordings and outside-out patch-clamp recordings from transiently transfected human embryonic kidney (HEK) cells.

RESULTS

Two de novo missense variants related to epilepsy were identified in the CLPTM1 gene. Functional studies indicated that CLPTM1-p.R454H and CLPTM1-p.R568Q variants reduced the γ-aminobutyric acid A receptor (GABAA R) current response amplitude recorded under voltage clamp compared to the wild-type receptors. These variants also reduced the charge transfer and altered the time course of desensitization and deactivation following rapid removal of GABA. The surface expression of the GABAA R γ2 subunit from the CLPTM1-p.R568Q group was significantly reduced compared to CLPTM1-WT.

SIGNIFICANCE

This is the first report of functionally relevant variants within the CLPTM1 gene. Patch-clamp recordings showed that these de novo CLPTM1 variants reduce GABAA R currents and charge transfer, which should promote excitation and hypersynchronous activity. This study may provide insights into the molecular mechanisms of the CLPTM1 variants underlying the patients' phenotypes, as well as for exploring potential therapeutic targets for epilepsy.

Classification of missense variants in the N-methyl-d-aspartate receptor GRIN gene family as gain- or loss-of-function

Advances in sequencing technology have generated a large amount of genetic data from patients with neurological conditions. These data have provided diagnosis of many rare diseases, including a number of pathogenic de novo missense variants in GRIN genes encoding N-methyl-d-aspartate receptors (NMDARs). To understand the ramifications for neurons and brain circuits affected by rare patient variants, functional analysis of the variant receptor is necessary in model systems. For NMDARs, this functional analysis needs to assess multiple properties in order to understand how variants could impact receptor function in neurons. One can then use these data to determine whether the overall actions will increase or decrease NMDAR-mediated charge transfer. Here, we describe an analytical and comprehensive framework by which to categorize GRIN variants as either gain-of-function (GoF) or loss-of-function (LoF) and apply this approach to GRIN2B variants identified in patients and the general population. This framework draws on results from six different assays that assess the impact of the variant on NMDAR sensitivity to agonists and endogenous modulators, trafficking to the plasma membrane, response time course and channel open probability. We propose to integrate data from multiple in vitro assays to arrive at a variant classification, and suggest threshold levels that guide confidence. The data supporting GoF and LoF determination are essential to assessing pathogenicity and patient stratification for clinical trials as personalized pharmacological and genetic agents that can enhance or reduce receptor function are advanced. This approach to functional variant classification can generalize to other disorders associated with missense variants.

Differential responses of disease-related GRIN variants located in pore-forming M2 domain of N-methyl-D-aspartate receptor to FDA-approved inhibitors

N-methyl-D-aspartate receptors (NMDAR), ionotropic glutamate receptors, mediate a slow component of excitatory synaptic transmission in the central nervous system and play a key role in normal brain function and development. Genetic variations in GRIN genes encoding NMDAR subunits that alter the receptor's functional characteristics are associated with a wide range of neurological and neuropsychiatric conditions. Pathological GRIN variants located in the M2 re-entrant loop lining the channel pore cause significant functional changes, the most consequential alteration being a reduction in voltage-dependent Mg2+ inhibition. Voltage-dependent Mg2+ block is a unique feature of NMDAR biology whereby channel activation requires both ligand binding and postsynaptic membrane depolarization. Thus, loss of NMDAR Mg2+ block will have a profound impact on synaptic function and plasticity. Here, we choose 11 missense variants within the GRIN1, GRIN2A, and GRIN2B genes that alter residues located in the M2 loop and significantly reduce Mg2+ inhibition. Each variant was evaluated for tolerance to genetic variation using the 3-dimensional structure and assessed for functional rescue pharmacology via electrophysiological recordings. Three FDA-approved NMDAR drugs-memantine, dextromethorphan, and ketamine-were chosen based on their ability to bind near the M2 re-entrant loop, potentially rectifying dysregulated NMDAR function by supplementing the reduced voltage-dependent Mg2+ block. These results provide insight of structural determinants of FDA-approved NMDAR drugs at their binding sites in the channel pore and may further define conditions necessary for the use of such agents as potential rescue pharmacology.

Functional effects of disease-associated variants reveal that the S1-M1 linker of the NMDA receptor critically controls channel opening

The short pre-M1 helix within the S1-M1 linker (also referred to as the pre-M1 linker) between the agonist-binding domain (ABD, S1) and the M1 transmembrane helix of the NMDA receptor (NMDAR) is devoid of missense variants within the healthy population but is a locus for de novo pathogenic variants associated with neurological disorders. Several de novo variants within this helix have been identified in patients presenting early in life with intellectual disability, developmental delay, and/or epilepsy. In this study, we evaluated functional properties for twenty variants within the pre-M1 linker in GRIN1, GRIN2A, and GRIN2B genes, including six novel missense variants. The effects of pre-M1 variants on agonist potency, sensitivity to endogenous allosteric modulators, response time course, channel open probability, and surface expression were assessed. Our data indicated that virtually all of the variants evaluated altered channel function, and multiple variants had profound functional consequences, which may contribute to the neurological conditions in the patients harboring the variants in this region. These data strongly suggest that the residues within the pre-M1 helix play a key role in channel gating and are highly intolerant to genetic variation.

Functional Assessment of Cancer-Linked Mutations in Sensitive Regions of Regulators of G Protein Signaling Predicted by Three-Dimensional Missense Tolerance Ratio Analysis

Regulators of G protein signaling (RGS) proteins modulate G protein-coupled receptor (GPCR) signaling by acting as negative regulators of G proteins. Genetic variants in RGS proteins are associated with many diseases, including cancers, although the impact of these mutations on protein function is uncertain. Here we analyze the RGS domains of 15 RGS protein family members using a novel bioinformatic tool that measures the missense tolerance ratio (MTR) using a three-dimensional (3D) structure (3DMTR). Subsequent permutation analysis can define the protein regions that are most significantly intolerant (P < 0.05) in each dataset. We further focused on RGS14, RGS10, and RGS4. RGS14 exhibited seven significantly tolerant and seven significantly intolerant residues, RGS10 had six intolerant residues, and RGS4 had eight tolerant and six intolerant residues. Intolerant and tolerant-control residues that overlap with pathogenic cancer mutations reported in the COSMIC cancer database were selected to define the functional phenotype. Using complimentary cellular and biochemical approaches, proteins were tested for effects on GPCR-Gα activation, Gα binding properties, and downstream cAMP levels. Identified intolerant residues with reported cancer-linked mutations RGS14-R173C/H and RGS4-K125Q/E126K, and tolerant RGS14-S127P and RGS10-S64T resulted in a loss-of-function phenotype in GPCR-G protein signaling activity. In downstream cAMP measurement, tolerant RGS14-D137Y and RGS10-S64T and intolerant RGS10-K89M resulted in change of function phenotypes. These findings show that 3DMTR identified intolerant residues that overlap with cancer-linked mutations cause phenotypic changes that negatively impact GPCR-G protein signaling and suggests that 3DMTR is a potentially useful bioinformatics tool for predicting functionally important protein residues. SIGNIFICANCE STATEMENT: Human genetic variant/mutation information has expanded rapidly in recent years, including cancer-linked mutations in regulator of G protein signaling (RGS) proteins. However, experimental testing of the impact of this vast catalogue of mutations on protein function is not feasible. We used the novel bioinformatics tool three-dimensional missense tolerance ratio (3DMTR) to define regions of genetic intolerance in RGS proteins and prioritize which cancer-linked mutants to test. We found that 3DMTR more accurately classifies loss-of-function mutations in RGS proteins than other databases thereby offering a valuable new research tool.

Compendium of proteins containing segments that exhibit zero-tolerance to amino acid variation in humans

Genetic missense tolerance ratio (MTR) analysis systematically evaluates all possible segments in a given protein-encoding transcript found in the human population. This method scores each segment for the number of observed missense variants versus the number of silent mutations in that same segment. An MTR score of 0 indicates that no missense mutations are observed within a given segment. This is indicative of evolutionary purifying selection, which excludes mutations in that segment from the general human population. Here, we conducted MTR analysis on each of the roughly 20,000 protein-encoding human genes. It was seen that there are 257 genes with at least one 31-residue encoding segment with MTR = 0 (1.3% of all human genes). The proteins encoded by these 257 genes were tabulated along with information regarding the sequence location of each intolerant segment, the likely function of the protein, and so forth. The most functionally-enriched family among these proteins is a collection of several dozen proteins that are directly involved in RNA splicing. Some of the other proteins with zero-tolerance segments have thus far escaped significant characterization. Indeed, while a number of these proteins have previously been genetically linked to human disorders, many have not. We hypothesize that this compendium of human proteins with zero-tolerance segments can be used to complement disease mutation data as a pointer to genes and proteins that are associated with interesting and underexplored human biology.

The 3D mutational constraint on amino acid sites in the human proteome

Quantification of the tolerance of protein sites to genetic variation has become a cornerstone of variant interpretation. We hypothesize that the constraint on missense variation at individual amino acid sites is largely shaped by direct interactions with 3D neighboring sites. To quantify this constraint, we introduce a framework called COntact Set MISsense tolerance (or COSMIS) and comprehensively map the landscape of 3D mutational constraint on 6.1 million amino acid sites covering 16,533 human proteins. We show that 3D mutational constraint is pervasive and that the level of constraint is strongly associated with disease relevance both at the site and the protein level. We demonstrate that COSMIS performs significantly better at variant interpretation tasks than other population-based constraint metrics while also providing structural insight into the functional roles of constrained sites. We anticipate that COSMIS will facilitate the interpretation of protein-coding variation in evolution and prioritization of sites for mechanistic investigation.

Integration of Protein Structure and Population-Scale DNA Sequence Data for Disease Gene Discovery and Variant Interpretation

The experimental and computational techniques for capturing information about protein structures and genetic variation within the human genome have advanced dramatically in the past 20 years, generating extensive new data resources. In this review, we discuss these advances, along with new approaches for determining the impact a genetic variant has on protein function. We focus on the potential of new methods that integrate human genetic variation into protein structures to discover relationships to disease, including the discovery of mutational hotspots in cancer-related proteins, the localization of protein-altering variants within protein regions for common complex diseases, and the assessment of variants of unknown significance for Mendelian traits. We expect that approaches that integrate these data sources will play increasingly important roles in disease gene discovery and variant interpretation. Expected final online publication date for the Annual Review of Biomedical Data Science, Volume 5 is August 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Opportunities for Precision Treatment of GRIN2A and GRIN2B Gain-of-Function Variants in Triheteromeric N-Methyl-D-Aspartate Receptors

N-methyl-D-aspartate receptors (NMDARs) are tetrameric assemblies of two glutamate N-methyl-D-aspartate receptor subunits, GluN1 and two GluN2, that mediate excitatory synaptic transmission in the central nervous system. Four genes (GRIN2A-D) encode four distinct GluN2 subunits (GluN2A-D). Thus, NMDARs can be diheteromeric assemblies of two GluN1 plus two identical GluN2 subunits, or triheteromeric assemblies of two GluN1 subunits plus two different GluN2 subunits. An increasing number of de novo GRIN variants have been identified in patients with neurologic conditions and with GRIN2A and GRIN2B harboring the vast majority (> 80%) of variants in these cases. These variants produce a wide range of effects on NMDAR function depending upon its subunit subdomain location and additionally on the subunit composition of diheteromeric versus triheteromeric NMDARs. Increasing evidence implicates triheteromeric GluN1/GluN2A/GluN2B receptors as a major component of the NMDAR pool in the adult cortex and hippocampus. Here, we explore the ability of GluN2A- and GluN2B-selective inhibitors to reduce excess current flow through triheteromeric GluN1/GluN2A/GluN2B receptors that contain one copy of GRIN2A or GRIN2B gain-of-function variants. Our data reveal a broad range of sensitivities for variant-containing triheteromeric receptors to subunit-selective inhibitors, with some variants still showing strong sensitivity to inhibitors, whereas others are relatively insensitive. Most variants, however, retain sensitivity to non-selective channel blockers and the competitive antagonist D-(-)-2-amino-5-phosphonopentanoic acid. These results suggest that with comprehensive analysis, certain disease-related GRIN2A and GRIN2B variants can be identified as potential targets for subunit-selective modulation and potential therapeutic gain. SIGNIFICANCE STATEMENT: Triheteromeric NMDA receptors that contain one copy each of the GluN2A and GluN2B subunits show intermediate sensitivity to GluN2A- and GluN2B-selective inhibitors, making these compounds candidates for attenuating overactive, GRIN variant-containing NMDA receptors associated with neurological conditions. We show that functional evaluation of variant properties with inhibitor pharmacology can support selection of a subset of variants for which GluN2 subunit-selective agents remain effective inhibitors of variant-containing triheteromeric NMDA receptors.

Structure, Function, and Pharmacology of Glutamate Receptor Ion Channels

Many physiologic effects of l-glutamate, the major excitatory neurotransmitter in the mammalian central nervous system, are mediated via signaling by ionotropic glutamate receptors (iGluRs). These ligand-gated ion channels are critical to brain function and are centrally implicated in numerous psychiatric and neurologic disorders. There are different classes of iGluRs with a variety of receptor subtypes in each class that play distinct roles in neuronal functions. The diversity in iGluR subtypes, with their unique functional properties and physiologic roles, has motivated a large number of studies. Our understanding of receptor subtypes has advanced considerably since the first iGluR subunit gene was cloned in 1989, and the research focus has expanded to encompass facets of biology that have been recently discovered and to exploit experimental paradigms made possible by technological advances. Here, we review insights from more than 3 decades of iGluR studies with an emphasis on the progress that has occurred in the past decade. We cover structure, function, pharmacology, roles in neurophysiology, and therapeutic implications for all classes of receptors assembled from the subunits encoded by the 18 ionotropic glutamate receptor genes. SIGNIFICANCE STATEMENT: Glutamate receptors play important roles in virtually all aspects of brain function and are either involved in mediating some clinical features of neurological disease or represent a therapeutic target for treatment. Therefore, understanding the structure, function, and pharmacology of this class of receptors will advance our understanding of many aspects of brain function at molecular, cellular, and system levels and provide new opportunities to treat patients.