Crystallographic fragment screening delivers diverse chemical scaffolds for Zika virus NS2B-NS3 protease inhibitor development.

Zika virus (ZIKV) infections cause microcephaly in new-borns and Guillain-Barre syndrome in adults raising a significant global public health concern, yet no vaccines or antiviral drugs have been developed to prevent or treat ZIKV infections. The viral protease NS3 and its co-factor NS2B are essential for the cleavage of the Zika polyprotein precursor into individual structural and non-structural proteins and is therefore an attractive drug target. Generation of a robust crystal system of co-expressed NS2B-NS3 protease has enabled us to perform a crystallographic fragment screening campaign with 1076 fragments. 48 binders with diverse chemical scaffolds were identified in the active site of the protease, with another 6 fragment hits observed in a potential allosteric binding site. Our work provides potential starting points for the development of potent NS2B-NS3 protease inhibitors. Furthermore, we have structurally characterized a potential allosteric binding pocket, identifying opportunities for allosteric inhibitor development.

De novo variants in DENND5B cause a neurodevelopmental disorder.

The Rab family of guanosine triphosphatases (GTPases) includes key regulators of intracellular transport and membrane trafficking targeting specific steps in exocytic, endocytic, and recycling pathways. DENND5B (Rab6-interacting Protein 1B-like protein, R6IP1B) is the longest isoform of DENND5, an evolutionarily conserved DENN domain-containing guanine nucleotide exchange factor (GEF) that is highly expressed in the brain. Through exome sequencing and international matchmaking platforms, we identified five de novo variants in DENND5B in a cohort of five unrelated individuals with neurodevelopmental phenotypes featuring cognitive impairment, dysmorphism, abnormal behavior, variable epilepsy, white matter abnormalities, and cortical gyration defects. We used biochemical assays and confocal microscopy to assess the impact of DENND5B variants on protein accumulation and distribution. Then, exploiting fluorescent lipid cargoes coupled to high-content imaging and analysis in living cells, we investigated whether DENND5B variants affected the dynamics of vesicle-mediated intracellular transport of specific cargoes. We further generated an in silico model to investigate the consequences of DENND5B variants on the DENND5B-RAB39A interaction. Biochemical analysis showed decreased protein levels of DENND5B mutants in various cell types. Functional investigation of DENND5B variants revealed defective intracellular vesicle trafficking, with significant impairment of lipid uptake and distribution. Although none of the variants affected the DENND5B-RAB39A interface, all were predicted to disrupt protein folding. Overall, our findings indicate that DENND5B variants perturb intracellular membrane trafficking pathways and cause a complex neurodevelopmental syndrome with variable epilepsy and white matter involvement.

Open science discovery of potent noncovalent SARS-CoV-2 main protease inhibitors.

We report the results of the COVID Moonshot, a fully open-science, crowdsourced, and structure-enabled drug discovery campaign targeting the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) main protease. We discovered a noncovalent, nonpeptidic inhibitor scaffold with lead-like properties that is differentiated from current main protease inhibitors. Our approach leveraged crowdsourcing, machine learning, exascale molecular simulations, and high-throughput structural biology and chemistry. We generated a detailed map of the structural plasticity of the SARS-CoV-2 main protease, extensive structure-activity relationships for multiple chemotypes, and a wealth of biochemical activity data. All compound designs (>18,000 designs), crystallographic data (>490 ligand-bound x-ray structures), assay data (>10,000 measurements), and synthesized molecules (>2400 compounds) for this campaign were shared rapidly and openly, creating a rich, open, and intellectual property-free knowledge base for future anticoronavirus drug discovery.

Structural and non-coding variants increase the diagnostic yield of clinical whole genome sequencing for rare diseases.

BACKGROUND: Whole genome sequencing is increasingly being used for the diagnosis of patients with rare diseases. However, the diagnostic yields of many studies, particularly those conducted in a healthcare setting, are often disappointingly low, at 25-30%. This is in part because although entire genomes are sequenced, analysis is often confined to in silico gene panels or coding regions of the genome. METHODS: We undertook WGS on a cohort of 122 unrelated rare disease patients and their relatives (300 genomes) who had been pre-screened by gene panels or arrays. Patients were recruited from a broad spectrum of clinical specialties. We applied a bioinformatics pipeline that would allow comprehensive analysis of all variant types. We combined established bioinformatics tools for phenotypic and genomic analysis with our novel algorithms (SVRare, ALTSPLICE and GREEN-DB) to detect and annotate structural, splice site and non-coding variants. RESULTS: Our diagnostic yield was 43/122 cases (35%), although 47/122 cases (39%) were considered solved when considering novel candidate genes with supporting functional data into account. Structural, splice site and deep intronic variants contributed to 20/47 (43%) of our solved cases. Five genes that are novel, or were novel at the time of discovery, were identified, whilst a further three genes are putative novel disease genes with evidence of causality. We identified variants of uncertain significance in a further fourteen candidate genes. The phenotypic spectrum associated with RMND1 was expanded to include polymicrogyria. Two patients with secondary findings in FBN1 and KCNQ1 were confirmed to have previously unidentified Marfan and long QT syndromes, respectively, and were referred for further clinical interventions. Clinical diagnoses were changed in six patients and treatment adjustments made for eight individuals, which for five patients was considered life-saving. CONCLUSIONS: Genome sequencing is increasingly being considered as a first-line genetic test in routine clinical settings and can make a substantial contribution to rapidly identifying a causal aetiology for many patients, shortening their diagnostic odyssey. We have demonstrated that structural, splice site and intronic variants make a significant contribution to diagnostic yield and that comprehensive analysis of the entire genome is essential to maximise the value of clinical genome sequencing.

Biallelic MED27 variants lead to variable ponto-cerebello-lental degeneration with movement disorders.

MED27 is a subunit of the Mediator multiprotein complex, which is involved in transcriptional regulation. Biallelic MED27 variants have recently been suggested to be responsible for an autosomal recessive neurodevelopmental disorder with spasticity, cataracts and cerebellar hypoplasia. We further delineate the clinical phenotype of MED27-related disease by characterizing the clinical and radiological features of 57 affected individuals from 30 unrelated families with biallelic MED27 variants. Using exome sequencing and extensive international genetic data sharing, 39 unpublished affected individuals from 18 independent families with biallelic missense variants in MED27 have been identified (29 females, mean age at last follow-up 17 ± 12.4 years, range 0.1-45). Follow-up and hitherto unreported clinical features were obtained from the published 12 families. Brain MRI scans from 34 cases were reviewed. MED27-related disease manifests as a broad phenotypic continuum ranging from developmental and epileptic-dyskinetic encephalopathy to variable neurodevelopmental disorder with movement abnormalities. It is characterized by mild to profound global developmental delay/intellectual disability (100%), bilateral cataracts (89%), infantile hypotonia (74%), microcephaly (62%), gait ataxia (63%), dystonia (61%), variably combined with epilepsy (50%), limb spasticity (51%), facial dysmorphism (38%) and death before reaching adulthood (16%). Brain MRI revealed cerebellar atrophy (100%), white matter volume loss (76.4%), pontine hypoplasia (47.2%) and basal ganglia atrophy with signal alterations (44.4%). Previously unreported 39 affected individuals had seven homozygous pathogenic missense MED27 variants, five of which were recurrent. An emerging genotype-phenotype correlation was observed. This study provides a comprehensive clinical-radiological description of MED27-related disease, establishes genotype-phenotype and clinical-radiological correlations and suggests a differential diagnosis with syndromes of cerebello-lental neurodegeneration and other subtypes of 'neuro-MEDopathies'.

The prevalence and phenotypic range associated with biallelic PKDCC variants.

PKDCC encodes a component of Hedgehog signalling required for normal chondrogenesis and skeletal development. Although biallelic PKDCC variants have been implicated in rhizomelic shortening of limbs with variable dysmorphic features, this association was based on just two patients. In this study, data from the 100 000 Genomes Project was used in conjunction with exome sequencing and panel-testing results accessed via international collaboration to assemble a cohort of eight individuals from seven independent families with biallelic PKDCC variants. The allelic series included six frameshifts, a previously described splice-donor site variant and a likely pathogenic missense variant observed in two families that was supported by in silico structural modelling. Database queries suggested that the prevalence of this condition is between 1 of 127 and 1 of 721 in clinical cohorts with skeletal dysplasia of unknown aetiology. Clinical assessments, combined with data from previously published cases, indicate a predominantly upper limb involvement. Micrognathia, hypertelorism and hearing loss appear to be commonly co-occurring features. In conclusion, this study strengthens the link between biallelic inactivation of PKDCC and rhizomelic limb-shortening and will enable clinical testing laboratories to better interpret variants in this gene.

Iterative computational design and crystallographic screening identifies potent inhibitors targeting the Nsp3 macrodomain of SARS-CoV-2.

The nonstructural protein 3 (NSP3) of the severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2) contains a conserved macrodomain enzyme (Mac1) that is critical for pathogenesis and lethality. While small-molecule inhibitors of Mac1 have great therapeutic potential, at the outset of the COVID-19 pandemic, there were no well-validated inhibitors for this protein nor, indeed, the macrodomain enzyme family, making this target a pharmacological orphan. Here, we report the structure-based discovery and development of several different chemical scaffolds exhibiting low- to sub-micromolar affinity for Mac1 through iterations of computer-aided design, structural characterization by ultra-high-resolution protein crystallography, and binding evaluation. Potent scaffolds were designed with in silico fragment linkage and by ultra-large library docking of over 450 million molecules. Both techniques leverage the computational exploration of tangible chemical space and are applicable to other pharmacological orphans. Overall, 160 ligands in 119 different scaffolds were discovered, and 153 Mac1-ligand complex crystal structures were determined, typically to 1 Å resolution or better. Our analyses discovered selective and cell-permeable molecules, unexpected ligand-mediated conformational changes within the active site, and key inhibitor motifs that will template future drug development against Mac1.

Biallelic PRMT7 pathogenic variants are associated with a recognizable syndromic neurodevelopmental disorder with short stature, obesity, and craniofacial and digital abnormalities.

PURPOSE: Protein arginine methyltransferase 7 (PRMT7) is a member of a family of enzymes that catalyzes the methylation of arginine residues on several protein substrates. Biallelic pathogenic PRMT7 variants have previously been associated with a syndromic neurodevelopmental disorder characterized by short stature, brachydactyly, intellectual developmental disability, and seizures. To our knowledge, no comprehensive study describes the detailed clinical characteristics of this syndrome. Thus, we aim to delineate the phenotypic spectrum of PRMT7-related disorder. METHODS: We assembled a cohort of 51 affected individuals from 39 different families, gathering clinical information from 36 newly described affected individuals and reviewing data of 15 individuals from the literature. RESULTS: The main clinical characteristics of the PRMT7-related syndrome are short stature, mild to severe developmental delay/intellectual disability, hypotonia, brachydactyly, and distinct facial morphology, including bifrontal narrowing, prominent supraorbital ridges, sparse eyebrows, short nose with full/broad nasal tip, thin upper lip, full and everted lower lip, and a prominent or squared-off jaw. Additional variable findings include seizures, obesity, nonspecific magnetic resonance imaging abnormalities, eye abnormalities (i.e., strabismus or nystagmus), and hearing loss. CONCLUSION: This study further delineates and expands the molecular, phenotypic spectrum and natural history of PRMT7-related syndrome characterized by a neurodevelopmental disorder with skeletal, growth, and endocrine abnormalities.

Iterative computational design and crystallographic screening identifies potent inhibitors targeting the Nsp3 Macrodomain of SARS-CoV-2.

The nonstructural protein 3 (NSP3) of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) contains a conserved macrodomain enzyme (Mac1) that is critical for pathogenesis and lethality. While small molecule inhibitors of Mac1 have great therapeutic potential, at the outset of the COVID-19 pandemic there were no well-validated inhibitors for this protein nor, indeed, the macrodomain enzyme family, making this target a pharmacological orphan. Here, we report the structure-based discovery and development of several different chemical scaffolds exhibiting low- to sub-micromolar affinity for Mac1 through iterations of computer-aided design, structural characterization by ultra-high resolution protein crystallography, and binding evaluation. Potent scaffolds were designed with in silico fragment linkage and by ultra-large library docking of over 450 million molecules. Both techniques leverage the computational exploration of tangible chemical space and are applicable to other pharmacological orphans. Overall, 160 ligands in 119 different scaffolds were discovered, and 152 Mac1-ligand complex crystal structures were determined, typically to 1 Å resolution or better. Our analyses discovered selective and cell-permeable molecules, unexpected ligand-mediated protein dynamics within the active site, and key inhibitor motifs that will template future drug development against Mac1.

Venus: Elucidating the Impact of Amino Acid Variants on Protein Function Beyond Structure Destabilisation.

Exploring the functional effect of a non-synonymous coding variant at the protein level requires multiple pieces of information to be interpreted appropriately. This is particularly important when embarking on the study of a potentially pathogenic variant linked to a rare or monogenic disease. Whereas accurate protein stability predictions alone are generally informative, other effects, such as disruption of post-translational modifications or weakened ligand binding, may also contribute to the disease phenotype. Furthermore, consideration of nearby variants that are found in the healthy population may strengthen or refute a given mechanistic hypothesis. Whilst there are several bioinformatics tools available that score a genetic variant in terms of deleteriousness, there is no single tool that assembles multiple effects of a variant on the encoded protein, beyond structural stability, and presents them on the structure for inspection. Venus is a web application which, given a protein substitution, rapidly estimates the predicted effect on protein stability of the variant, flags if the variant affects a post-translational modification site, a predicted linear motif or known annotation, and determines the effect on protein stability of variants which affect nearby residues and have been identified in healthy populations. Venus is built upon Michelanglo and the results can be exported to it, allowing them to be annotated and shared with other researchers. Venus is freely accessible at https://venus.cmd.ox.ac.uk and its source code is openly available at https://github.com/CMD-Oxford/Michelanglo-and-Venus.

Biallelic variants in SLC38A3 encoding a glutamine transporter cause epileptic encephalopathy.

The solute carrier (SLC) superfamily encompasses >400 transmembrane transporters involved in the exchange of amino acids, nutrients, ions, metals, neurotransmitters and metabolites across biological membranes. SLCs are highly expressed in the mammalian brain; defects in nearly 100 unique SLC-encoding genes (OMIM: https://www.omim.org) are associated with rare Mendelian disorders including developmental and epileptic encephalopathy and severe neurodevelopmental disorders. Exome sequencing and family-based rare variant analyses on a cohort with neurodevelopmental disorders identified two siblings with developmental and epileptic encephalopathy and a shared deleterious homozygous splicing variant in SLC38A3. The gene encodes SNAT3, a sodium-coupled neutral amino acid transporter and a principal transporter of the amino acids asparagine, histidine, and glutamine, the latter being the precursor for the neurotransmitters GABA and glutamate. Additional subjects with a similar developmental and epileptic encephalopathy phenotype and biallelic predicted-damaging SLC38A3 variants were ascertained through GeneMatcher and collaborations with research and clinical molecular diagnostic laboratories. Untargeted metabolomic analysis was performed to identify novel metabolic biomarkers. Ten individuals from seven unrelated families from six different countries with deleterious biallelic variants in SLC38A3 were identified. Global developmental delay, intellectual disability, hypotonia, and absent speech were common features while microcephaly, epilepsy, and visual impairment were present in the majority. Epilepsy was drug-resistant in half. Metabolomic analysis revealed perturbations of glutamate, histidine, and nitrogen metabolism in plasma, urine, and CSF of selected subjects, potentially representing biomarkers of disease. Our data support the contention that SLC38A3 is a novel disease gene for developmental and epileptic encephalopathy and illuminate the likely pathophysiology of the disease as perturbations in glutamine homeostasis.

DogCatcher allows loop-friendly protein-protein ligation.

There are many efficient ways to connect proteins at termini. However, connecting at a loop is difficult because of lower flexibility and variable environment. Here, we have developed DogCatcher, a protein that forms a spontaneous isopeptide bond with DogTag peptide. DogTag/DogCatcher was generated initially by splitting a Streptococcus pneumoniae adhesin. We optimized DogTag/DogCatcher through rational design and evolution, increasing reaction rate by 250-fold and establishing millimolar solubility of DogCatcher. When fused to a protein terminus, DogTag/DogCatcher reacts slower than SpyTag003/SpyCatcher003. However, inserted in loops of a fluorescent protein or enzyme, DogTag reacts much faster than SpyTag003. Like many membrane proteins, the ion channel TRPC5 has no surface-exposed termini. DogTag in a TRPC5 extracellular loop allowed normal calcium flux and specific covalent labeling on cells in 1 min. DogTag/DogCatcher reacts under diverse conditions, at nanomolar concentrations, and to 98% conversion. Loop-friendly ligation should expand the toolbox for creating protein architectures.

Recurrent de novo missense variants in GNB2 can cause syndromic intellectual disability.

PURPOSE: Binding proteins (G-proteins) mediate signalling pathways involved in diverse cellular functions and comprise Gα and Gßγ units. Human diseases have been reported for all five Gß proteins. A de novo missense variant in GNB2 was recently reported in one individual with developmental delay/intellectual disability (DD/ID) and dysmorphism. We aim to confirm GNB2 as a neurodevelopmental disease gene, and elucidate the GNB2-associated neurodevelopmental phenotype in a patient cohort. METHODS: We discovered a GNB2 variant in the index case via exome sequencing and sought individuals with GNB2 variants via international data-sharing initiatives. In silico modelling of the variants was assessed, along with multiple lines of evidence in keeping with American College of Medical Genetics and Genomics guidelines for interpretation of sequence variants. RESULTS: We identified 12 unrelated individuals with five de novo missense variants in GNB2, four of which are recurrent: p.(Ala73Thr), p.(Gly77Arg), p.(Lys89Glu) and p.(Lys89Thr). All individuals have DD/ID with variable dysmorphism and extraneurologic features. The variants are located at the universally conserved shared interface with the Gα subunit, which modelling suggests weaken this interaction. CONCLUSION: Missense variants in GNB2 cause a congenital neurodevelopmental disorder with variable syndromic features, broadening the spectrum of multisystem phenotypes associated with variants in genes encoding G-proteins.

Expanding the phenotypic spectrum of BCS1L-related mitochondrial disease.

OBJECTIVE: To delineate the full phenotypic spectrum of BCS1L-related disease, provide better understanding of the genotype-phenotype correlations and identify reliable prognostic disease markers. METHODS: We performed a retrospective multinational cohort study of previously unpublished patients followed in 15 centres from 10 countries. Patients with confirmed biallelic pathogenic BCS1L variants were considered eligible. Clinical, laboratory, neuroimaging and genetic data were analysed. Patients were stratified into different groups based on the age of disease onset, whether homozygous or compound heterozygous for the c.232A>G (p.Ser78Gly) variant, and those with other pathogenic BCS1L variants. RESULTS: Thirty-three patients were included. We found that growth failure, lactic acidosis, tubulopathy, hepatopathy and early death were more frequent in those with disease onset within the first month of life. In those with onset after 1 month, neurological features including movement disorders and seizures were more frequent. Novel phenotypes, particularly involving movement disorder, were identified in this group. The presence of the c.232A>G (p.Ser78Gly) variant was associated with significantly worse survival and exclusively found in those with disease onset within the first month of life, whilst other pathogenic BCS1L variants were more frequent in those with later symptom onset. INTERPRETATION: The phenotypic spectrum of BCS1L-related disease comprises a continuum of clinical features rather than a set of separate syndromic clinical identities. Age of onset defines BCS1L-related disease clinically and early presentation is associated with poor prognosis. Genotype correlates with phenotype in the presence of the c.232A>G (p.Ser78Gly) variant.

Biallelic PI4KA variants cause neurological, intestinal and immunological disease.

Phosphatidylinositol 4-kinase IIIα (PI4KIIIα/PI4KA/OMIM:600286) is a lipid kinase generating phosphatidylinositol 4-phosphate (PI4P), a membrane phospholipid with critical roles in the physiology of multiple cell types. PI4KIIIα's role in PI4P generation requires its assembly into a heterotetrameric complex with EFR3, TTC7 and FAM126. Sequence alterations in two of these molecular partners, TTC7 (encoded by TTC7A or TCC7B) and FAM126, have been associated with a heterogeneous group of either neurological (FAM126A) or intestinal and immunological (TTC7A) conditions. Here we show that biallelic PI4KA sequence alterations in humans are associated with neurological disease, in particular hypomyelinating leukodystrophy. In addition, affected individuals may present with inflammatory bowel disease, multiple intestinal atresia and combined immunodeficiency. Our cellular, biochemical and structural modelling studies indicate that PI4KA-associated phenotypical outcomes probably stem from impairment of PI4KIIIα-TTC7-FAM126's organ-specific functions, due to defective catalytic activity or altered intra-complex functional interactions. Together, these data define PI4KA gene alteration as a cause of a variable phenotypical spectrum and provide fundamental new insight into the combinatorial biology of the PI4KIIIα-FAM126-TTC7-EFR3 molecular complex.

Non-coding region variants upstream of MEF2C cause severe developmental disorder through three distinct loss-of-function mechanisms.

Clinical genetic testing of protein-coding regions identifies a likely causative variant in only around half of developmental disorder (DD) cases. The contribution of regulatory variation in non-coding regions to rare disease, including DD, remains very poorly understood. We screened 9,858 probands from the Deciphering Developmental Disorders (DDD) study for de novo mutations in the 5' untranslated regions (5' UTRs) of genes within which variants have previously been shown to cause DD through a dominant haploinsufficient mechanism. We identified four single-nucleotide variants and two copy-number variants upstream of MEF2C in a total of ten individual probands. We developed multiple bespoke and orthogonal experimental approaches to demonstrate that these variants cause DD through three distinct loss-of-function mechanisms, disrupting transcription, translation, and/or protein function. These non-coding region variants represent 23% of likely diagnoses identified in MEF2C in the DDD cohort, but these would all be missed in standard clinical genetics approaches. Nonetheless, these variants are readily detectable in exome sequence data, with 30.7% of 5' UTR bases across all genes well covered in the DDD dataset. Our analyses show that non-coding variants upstream of genes within which coding variants are known to cause DD are an important cause of severe disease and demonstrate that analyzing 5' UTRs can increase diagnostic yield. We also show how non-coding variants can help inform both the disease-causing mechanism underlying protein-coding variants and dosage tolerance of the gene.

A human importin-ß-related disorder: Syndromic thoracic aortic aneurysm caused by bi-allelic loss-of-function variants in IPO8.

Importin 8, encoded by IPO8, is a ubiquitously expressed member of the importin-ß protein family that translocates cargo molecules such as proteins, RNAs, and ribonucleoprotein complexes into the nucleus in a RanGTP-dependent manner. Current knowledge of the cargoes of importin 8 is limited, but TGF-ß signaling components such as SMAD1-4 have been suggested to be among them. Here, we report that bi-allelic loss-of-function variants in IPO8 cause a syndromic form of thoracic aortic aneurysm (TAA) with clinical overlap with Loeys-Dietz and Shprintzen-Goldberg syndromes. Seven individuals from six unrelated families showed a consistent phenotype with early-onset TAA, motor developmental delay, connective tissue findings, and craniofacial dysmorphic features. A C57BL/6N Ipo8 knockout mouse model recapitulates TAA development from 8-12 weeks onward in both sexes but most prominently shows ascending aorta dilatation with a propensity for dissection in males. Compliance assays suggest augmented passive stiffness of the ascending aorta in male Ipo8-/- mice throughout life. Immunohistological investigation of mutant aortic walls reveals elastic fiber disorganization and fragmentation along with a signature of increased TGF-ß signaling, as evidenced by nuclear pSmad2 accumulation. RT-qPCR assays of the aortic wall in male Ipo8-/- mice demonstrate decreased Smad6/7 and increased Mmp2 and Ccn2 (Ctgf) expression, reinforcing a role for dysregulation of the TGF-ß signaling pathway in TAA development. Because importin 8 is the most downstream TGF-ß-related effector implicated in TAA pathogenesis so far, it offers opportunities for future mechanistic studies and represents a candidate drug target for TAA.