Delineation of functionally essential protein regions for 242 neurodevelopmental genes.

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.

Characterization of intrinsically disordered regions in proteins informed by human genetic diversity.

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.

Comprehensive characterization of amino acid positions in protein structures reveals molecular effect of missense variants.

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.

Structural mapping of GABRB3 variants reveals genotype-phenotype correlations.

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.

MISCAST: MIssense variant to protein StruCture Analysis web SuiTe.

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.

Multimodal small-molecule screening for human prion protein binders.

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.

Biological concepts in human sodium channel epilepsies and their relevance in clinical practice.

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.

A dynamic checkpoint in oxidative lesion discrimination by formamidopyrimidine-DNA glycosylase.

In contrast to proteins recognizing small-molecule ligands, DNA-dependent enzymes cannot rely solely on interactions in the substrate-binding centre to achieve their exquisite specificity. It is widely believed that substrate recognition by such enzymes involves a series of conformational changes in the enzyme-DNA complex with sequential gates favoring cognate DNA and rejecting nonsubstrates. However, direct evidence for such mechanism is limited to a few systems. We report that discrimination between the oxidative DNA lesion, 8-oxoguanine (oxoG) and its normal counterpart, guanine, by the repair enzyme, formamidopyrimidine-DNA glycosylase (Fpg), likely involves multiple gates. Fpg uses an aromatic wedge to open the Watson-Crick base pair and everts the lesion into its active site. We used molecular dynamics simulations to explore the eversion free energy landscapes of oxoG and G by Fpg, focusing on structural and energetic details of oxoG recognition. The resulting energy profiles, supported by biochemical analysis of site-directed mutants disturbing the interactions along the proposed path, show that Fpg selectively facilitates eversion of oxoG by stabilizing several intermediate states, helping the rapidly sliding enzyme avoid full extrusion of every encountered base for interrogation. Lesion recognition through multiple gating intermediates may be a common theme in DNA repair enzymes.

Active destabilization of base pairs by a DNA glycosylase wedge initiates damage recognition.

Formamidopyrimidine-DNA glycosylase (Fpg) excises 8-oxoguanine (oxoG) from DNA but ignores normal guanine. We combined molecular dynamics simulation and stopped-flow kinetics with fluorescence detection to track the events in the recognition of oxoG by Fpg and its mutants with a key phenylalanine residue, which intercalates next to the damaged base, changed to either alanine (F110A) or fluorescent reporter tryptophan (F110W). Guanine was sampled by Fpg, as evident from the F110W stopped-flow traces, but less extensively than oxoG. The wedgeless F110A enzyme could bend DNA but failed to proceed further in oxoG recognition. Modeling of the base eversion with energy decomposition suggested that the wedge destabilizes the intrahelical base primarily through buckling both surrounding base pairs. Replacement of oxoG with abasic (AP) site rescued the activity, and calculations suggested that wedge insertion is not required for AP site destabilization and eversion. Our results suggest that Fpg, and possibly other DNA glycosylases, convert part of the binding energy into active destabilization of their substrates, using the energy differences between normal and damaged bases for fast substrate discrimination.

Kinetic and structural insights into the binding of histone deacetylase 1 and 2 (HDAC1, 2) inhibitors.

The structure-activity and structure-kinetic relationships of a series of novel and selective ortho-aminoanilide inhibitors of histone deacetylases (HDACs) 1 and 2 are described. Different kinetic and thermodynamic selectivity profiles were obtained by varying the moiety occupying an 11Å channel leading to the Zn(2+) catalytic pocket of HDACs 1 and 2, two paralogs with a high degree of structural similarity. The design of these novel inhibitors was informed by two ligand-bound crystal structures of truncated hHDAC2. BRD4884 and BRD7232 possess kinetic selectivity for HDAC1 versus HDAC2. We demonstrate that the binding kinetics of HDAC inhibitors can be tuned for individual isoforms in order to modulate target residence time while retaining functional activity and increased histone H4K12 and H3K9 acetylation in primary mouse neuronal cell culture assays. These chromatin modifiers, with tuned binding kinetic profiles, can be used to define the relation between target engagement requirements and the pharmacodynamic response of HDACs in different disease applications.

Ensemble-based docking using biased molecular dynamics.

Proteins are dynamic molecules, and understanding their movements, especially as they relate to molecular recognition and protein-ligand interactions, poses a significant challenge to structure-based drug discovery. In most instances, protein flexibility is underrepresented in computer-aided drug design due to uncertainties on how it should be accurately modeled as well as the computational cost associated with attempting to incorporate flexibility in the calculations. One approach that aims to address these issues is ensemble-based docking. With this technique, ligands are docked to an ensemble of rigid protein conformations. Molecular dynamics (MD) simulations can be used to generate the ensemble of protein conformations for the subsequent docking. Here we present a novel approach that uses biased-MD simulations to generate the docking ensemble. The MD simulations are biased toward an initial protein-ligand X-ray complex structure. The biasing maintains some of the original crystallographic pocket-ligand information and thereby enhances sampling of the more relevant conformational space of the protein. Resulting trajectories are clustered to select a representative set of protein conformations, and ligands are docked to that reduced set of conformations. Cross-docking to this ensemble and then selecting the lowest scoring pose enables reliable identification of the correct binding mode. Various levels of biasing are investigated, and the method is validated for cyclin-dependent kinase 2 and factor Xa.

Fragment-based lead discovery and design.

Fragment-based lead discovery and design has and continues to show increasing promise in drug discovery. In this article, the current state of the art in terms of hot-spot characterization, fragment screening techniques, and fragment-based design is discussed. Three overall fragment-based lead generation strategies are explored and involve the chemical biology characterization of biological targets via fragment screening, fragment screening as a complementary approach to high-throughput screening of drug-like compounds, and direct fragment-based drug discovery, respectively. The evolution and development of fragment libraries is described. With an emphasis on computational approaches and the strategies applied at AstraZeneca, the review illustrates how integration of data from one regime can inform the design of experiments in the other, ultimately leading to the discovery of high quality chemical matter.

Genomics 2 Proteins portal: A resource and discovery tool for linking genetic screening outputs to protein sequences and structures.

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.

Differential functional consequences of GRIN2A mutations associated with schizophrenia and neurodevelopmental disorders.

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.

A Fast-Binding, Functionally Reversible, COX-2 Radiotracer for CNS PET Imaging.

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.

Allosteric inhibition of PPM1D serine/threonine phosphatase via an altered conformational state.

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.

Energetic preference of 8-oxoG eversion pathways in a DNA glycosylase.

Base eversion is a fundamental process in the biochemistry of nucleic acids, allowing proteins engaged in DNA repair and epigenetic modifications to access target bases in DNA. Crystal structures reveal end points of these processes, but not the pathways involved in the dynamic process of base recognition. To elucidate the pathway taken by 8-oxoguanine during base excision repair by Fpg, we calculated free energy surfaces during eversion of the damaged base through the major and minor grooves. The minor groove pathway and free energy barrier (6-7 kcal/mol) are consistent with previously reported results (Qi, Y.; Spong, M. C.; Nam, K.; Banerjee, A.; Jiralerspong, S.; Karplus, M.; Verdine, G. L. Nature 2009, 462, 762.) However, eversion of 8-oxoG through the major groove encounters a significantly lower barrier (3-4 kcal/mol) more consistent with experimentally determined rates of enzymatic sliding during lesion search (Blainey, P. C.; van Oijent, A. M.; Banerjee, A.; Verdine, G. L.; Xie, X. S. Proc. Natl. Acad. Sci. U.S.A. 2006, 103, 5752.). Major groove eversion has been suggested for other glycosylases, suggesting that in addition to function, dynamics of base eversion may also be conserved.

Oculocutaneous albinism type 1B associated with a functionally significant tyrosinase gene polymorphism detected with Whole Exome Sequencing.

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.

Regulation of purine metabolism connects KCTD13 to a metabolic disorder with autistic features.

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.

Discovery of a First-in-Class Inhibitor of the PRMT5-Substrate Adaptor Interaction.

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.