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.

Clinico-radiological features, molecular spectrum, and identification of prognostic factors in developmental and epileptic encephalopathy due to inosine triphosphate pyrophosphatase (ITPase) deficiency

Developmental and epileptic encephalopathy 35 (DEE 35) is a severe neurological condition caused by biallelic variants in ITPA, encoding inosine triphosphate pyrophosphatase, an essential enzyme in purine metabolism. We delineate the genotypic and phenotypic spectrum of DEE 35, analyzing possible predictors for adverse clinical outcomes. We investigated a cohort of 28 new patients and reviewed previously described cases, providing a comprehensive characterization of 40 subjects. Exome sequencing was performed to identify underlying ITPA pathogenic variants. Brain MRI (magnetic resonance imaging) scans were systematically analyzed to delineate the neuroradiological spectrum. Survival curves according to the Kaplan-Meier method and log-rank test were used to investigate outcome predictors in different subgroups of patients. We identified 18 distinct ITPA pathogenic variants, including 14 novel variants, and two deletions. All subjects showed profound developmental delay, microcephaly, and refractory epilepsy followed by neurodevelopmental regression. Brain MRI revision revealed a recurrent pattern of delayed myelination and restricted diffusion of early myelinating structures. Congenital microcephaly and cardiac involvement were statistically significant novel clinical predictors of adverse outcomes. We refined the molecular, clinical, and neuroradiological characterization of ITPase deficiency, and identified new clinical predictors which may have a potentially important impact on diagnosis, counseling, and follow-up of affected individuals.

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.

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.

•

Spontaneous transamidation at internal sites harnessing a DogTag/DogCatcher pair

•

DogCatcher is designed and bred for high solubility and rapid reaction

•

Within protein loops DogTag can clamp on its partner faster than SpyTag003

•

Fast and faithful fluorescent labeling of an ion channel at the cell surface via DogTag

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.

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.

Fragment binding to the Nsp3 macrodomain of SARS-CoV-2 identified through crystallographic screening and computational docking

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) macrodomain within the nonstructural protein 3 counteracts host-mediated antiviral adenosine diphosphate-ribosylation signaling. This enzyme is a promising antiviral target because catalytic mutations render viruses nonpathogenic. Here, we report a massive crystallographic screening and computational docking effort, identifying new chemical matter primarily targeting the active site of the macrodomain. Crystallographic screening of 2533 diverse fragments resulted in 214 unique macrodomain-binders. An additional 60 molecules were selected from docking more than 20 million fragments, of which 20 were crystallographically confirmed. X-ray data collection to ultra-high resolution and at physiological temperature enabled assessment of the conformational heterogeneity around the active site. Several fragment hits were confirmed by solution binding using three biophysical techniques (differential scanning fluorimetry, homogeneous time-resolved fluorescence, and isothermal titration calorimetry). The 234 fragment structures explore a wide range of chemotypes and provide starting points for development of potent SARS-CoV-2 macrodomain inhibitors.

An ancestral 10-bp repeat expansion in VWA1 causes recessive hereditary motor neuropathy

The extracellular matrix comprises a network of macromolecules such as collagens, proteoglycans and glycoproteins. VWA1 (von Willebrand factor A domain containing 1) encodes a component of the extracellular matrix that interacts with perlecan/collagen VI, appears to be involved in stabilizing extracellular matrix structures, and demonstrates high expression levels in tibial nerve. Vwa1-deficient mice manifest with abnormal peripheral nerve structure/function; however, VWA1 variants have not previously been associated with human disease. By interrogating the genome sequences of 74 180 individuals from the 100K Genomes Project in combination with international gene-matching efforts and targeted sequencing, we identified 17 individuals from 15 families with an autosomal-recessive, non-length dependent, hereditary motor neuropathy and rare biallelic variants in VWA1. A single disease-associated allele p.(G25Rfs*74), a 10-bp repeat expansion, was observed in 14/15 families and was homozygous in 10/15. Given an allele frequency in European populations approaching 1/1000, the seven unrelated homozygote individuals ascertained from the 100K Genomes Project represents a substantial enrichment above expected. Haplotype analysis identified a shared 220 kb region suggesting that this founder mutation arose >7000 years ago. A wide age-range of patients (6-83 years) helped delineate the clinical phenotype over time. The commonest disease presentation in the cohort was an early-onset (mean 2.0 ± 1.4 years) non-length-dependent axonal hereditary motor neuropathy, confirmed on electrophysiology, which will have to be differentiated from other predominantly or pure motor neuropathies and neuronopathies. Because of slow disease progression, ambulation was largely preserved. Neurophysiology, muscle histopathology, and muscle MRI findings typically revealed clear neurogenic changes with single isolated cases displaying additional myopathic process. We speculate that a few findings of myopathic changes might be secondary to chronic denervation rather than indicating an additional myopathic disease process. Duplex reverse transcription polymerase chain reaction and immunoblotting using patient fibroblasts revealed that the founder allele results in partial nonsense mediated decay and an absence of detectable protein. CRISPR and morpholino vwa1 modelling in zebrafish demonstrated reductions in motor neuron axonal growth, synaptic formation in the skeletal muscles and locomotive behaviour. In summary, we estimate that biallelic variants in VWA1 may be responsible for up to 1% of unexplained hereditary motor neuropathy cases in Europeans. The detailed clinical characterization provided here will facilitate targeted testing on suitable patient cohorts. This novel disease gene may have previously evaded detection because of high GC content, consequential low coverage and computational difficulties associated with robustly detecting repeat-expansions. Reviewing previously unsolved exomes using lower QC filters may generate further diagnoses.

Fragment Binding to the Nsp3 Macrodomain of SARS-CoV-2 Identified Through Crystallographic Screening and Computational Docking

The SARS-CoV-2 macrodomain (Mac1) within the non-structural protein 3 (Nsp3) counteracts host-mediated antiviral ADP-ribosylation signalling. This enzyme is a promising antiviral target because catalytic mutations render viruses non-pathogenic. Here, we report a massive crystallographic screening and computational docking effort, identifying new chemical matter primarily targeting the active site of the macrodomain. Crystallographic screening of diverse fragment libraries resulted in 214 unique macrodomain-binding fragments, out of 2,683 screened. An additional 60 molecules were selected from docking over 20 million fragments, of which 20 were crystallographically confirmed. X-ray data collection to ultra-high resolution and at physiological temperature enabled assessment of the conformational heterogeneity around the active site. Several crystallographic and docking fragment hits were validated for solution binding using three biophysical techniques (DSF, HTRF, ITC). Overall, the 234 fragment structures presented explore a wide range of chemotypes and provide starting points for development of potent SARS-CoV-2 macrodomain inhibitors.

MichelaNglo: sculpting protein views on web pages without coding

MOTIVATION

The sharing of macromolecular structural information online by scientists is predominantly performed via 2D static images, since the embedding of interactive 3D structures in webpages is non-trivial. Whilst the technologies to do so exist, they are often only implementable with significant web coding experience.

RESULTS

Michelaɴɢʟo is an accessible and open-source web-based application that supports the generation, customization and sharing of interactive 3D macromolecular visualizations for digital media without requiring programming skills. A PyMOL file, PDB file, PDB identifier code or protein/gene name can be provided to form the basis of visualizations using the NGL JavaScript library. Hyperlinks that control the view can be added to text within the page. Protein-coding variants can be highlighted to support interpretation of their potential functional consequences. The resulting visualizations and text can be customized and shared, as well as embedded within existing websites by following instructions and using a self-contained download. Michelaɴɢʟo allows researchers to move away from static images and instead engage, describe and explain their protein to a wider audience in a more interactive fashion.

AVAILABILITY AND IMPLEMENTATION

Michelaɴɢʟo is hosted at michelanglo.sgc.ox.ac.uk. The Python code is freely available at https://github.com/thesgc/MichelaNGLo, along with documentations about its implementation.

Identification of Circulating Genomic and Metabolic Biomarkers in Intrahepatic Cholangiocarcinoma

Intrahepatic cholangiocarcinoma (ICC) is an aggressive cancer arising from the bile ducts with a need for earlier diagnosis and a greater range of treatment options. KRAS/NRAS mutations are common in ICC tumours and 6-32% of patients also have isocitrate dehydrogenase 1 and 2 (IDH1 and IDH2) gene mutations associated with metabolic changes. This feasibility study investigated sequencing circulating tumour DNA (ctDNA) combined with metabolite profiling of plasma as a method for biomarker discovery in ICC patients. Plasma was collected from four ICC patients receiving radio-embolisation and healthy controls at multiple time points. ctDNA was sequenced using Ampliseq cancer hotspot panel-v2 on Ion Torrent PGM for single nucleotide variants (SNV) detection and with Illumina whole genome sequencing for copy number variants (CNV) and further targeted examination for SNVs. Untargeted analysis of metabolites from patient and control plasma was performed using liquid chromatography coupled with high-resolution tandem mass spectrometry (LC-MS/MS). Metabolite identification was performed using multi-parameter comparisons with analysis of authentic standards, and univariate statistical analysis was performed to identify differences in metabolite abundance between patient and control samples. Recurrent somatic SNVs and CNVs were identified in ctDNA from three out of four patients that included both NRAS and IDH1 mutations linked to ICC. Plasma metabolite analysis revealed biomarker metabolites associated with ICC and in particular 2-hydroxyglutarate (2-HG) levels were elevated in both samples from the only patient showing a variant allele in IDH1. A reduction in the number of CNVs was observed with treatment. This study demonstrates that ctDNA and metabolite levels can be identified and correlated in ICC patient blood samples and differentiated from healthy controls. We conclude that combining genomic and metabolic analysis of plasma offers an effective approach to biomarker identification with potential for disease stratification and early detection studies.

Delineation of dominant and recessive forms of LZTR1-associated Noonan syndrome

Noonan syndrome (NS) is characterised by distinctive facial features, heart defects, variable degrees of intellectual disability and other phenotypic manifestations. Although the mode of inheritance is typically dominant, recent studies indicate LZTR1 may be associated with both dominant and recessive forms. Seeking to describe the phenotypic characteristics of LZTR1-associated NS, we searched for likely pathogenic variants using two approaches. First, scrutiny of exomes from 9624 patients recruited by the Deciphering Developmental Disorders (DDDs) study uncovered six dominantly-acting mutations (p.R97L; p.Y136C; p.Y136H, p.N145I, p.S244C; p.G248R) of which five arose de novo, and three patients with compound-heterozygous variants (p.R210*/p.V579M; p.R210*/p.D531N; c.1149+1G>T/p.R688C). One patient also had biallelic loss-of-function mutations in NEB, consistent with a composite phenotype. After removing this complex case, analysis of human phenotype ontology terms indicated significant phenotypic similarities (P = 0.0005), supporting a causal role for LZTR1. Second, targeted sequencing of eight unsolved NS-like cases identified biallelic LZTR1 variants in three further subjects (p.W469*/p.Y749C, p.W437*/c.-38T>A and p.A461D/p.I462T). Our study strengthens the association of LZTR1 with NS, with de novo mutations clustering around the KT1-4 domains. Although LZTR1 variants explain ~0.1% of cases across the DDD cohort, the gene is a relatively common cause of unsolved NS cases where recessive inheritance is suspected.

Evolving Accelerated Amidation by SpyTag/SpyCatcher to Analyze Membrane Dynamics

SpyTag is a peptide that forms a spontaneous amide bond with its protein partner SpyCatcher. This protein superglue is a broadly useful tool for molecular assembly, locking together biological building blocks efficiently and irreversibly in diverse architectures. We initially developed SpyTag and SpyCatcher by rational design, through splitting a domain from a Gram-positive bacterial adhesin. In this work, we established a phage-display platform to select for specific amidation, leading to an order of magnitude acceleration for interaction of the SpyTag002 variant with the SpyCatcher002 variant. We show that the 002 pair bonds rapidly under a wide range of conditions and at either protein terminus. SpyCatcher002 was fused to an intimin derived from enterohemorrhagic Escherichia coli. SpyTag002 reaction enabled specific and covalent decoration of intimin for live cell fluorescent imaging of the dynamics of the bacterial outer membrane as cells divide.

Bacterial methionine biosynthesis

Methionine is essential in all organisms, as it is both a proteinogenic amino acid and a component of the cofactor, S-adenosyl methionine. The metabolic pathway for its biosynthesis has been extensively characterized in Escherichia coli; however, it is becoming apparent that most bacterial species do not use the E. coli pathway. Instead, studies on other organisms and genome sequencing data are uncovering significant diversity in the enzymes and metabolic intermediates that are used for methionine biosynthesis. This review summarizes the different biochemical strategies that are employed in the three key steps for methionine biosynthesis from homoserine (i.e. acylation, sulfurylation and methylation). A survey is presented of the presence and absence of the various biosynthetic enzymes in 1593 representative bacterial species, shedding light on the non-canonical nature of the E. coli pathway. This review also highlights ways in which knowledge of methionine biosynthesis can be utilized for biotechnological applications. Finally, gaps in the current understanding of bacterial methionine biosynthesis are noted. For example, the paper discusses the presence of one gene (metC) in a large number of species that appear to lack the gene encoding the enzyme for the preceding step in the pathway (metB), as it is understood in E. coli. Therefore, this review aims to move the focus away from E. coli, to better reflect the true diversity of bacterial pathways for methionine biosynthesis.

New rRNA gene-based phylogenies of the Alphaproteobacteria provide perspective on major groups, mitochondrial ancestry and phylogenetic instability

Bacteria in the class Alphaproteobacteria have a wide variety of lifestyles and physiologies. They include pathogens of humans and livestock, agriculturally valuable strains, and several highly abundant marine groups. The ancestor of mitochondria also originated in this clade. Despite significant effort to investigate the phylogeny of the Alphaproteobacteria with a variety of methods, there remains considerable disparity in the placement of several groups. Recent emphasis on phylogenies derived from multiple protein-coding genes remains contentious due to disagreement over appropriate gene selection and the potential influences of systematic error. We revisited previous investigations in this area using concatenated alignments of the small and large subunit (SSU and LSU) rRNA genes, as we show here that these loci have much lower GC bias than whole genomes. This approach has allowed us to update the canonical 16S rRNA gene tree of the Alphaproteobacteria with additional important taxa that were not previously included, and with added resolution provided by concatenating the SSU and LSU genes. We investigated the topological stability of the Alphaproteobacteria by varying alignment methods, rate models, taxon selection and RY-recoding to circumvent GC content bias. We also introduce RYMK-recoding and show that it avoids some of the information loss in RY-recoding. We demonstrate that the topology of the Alphaproteobacteria is sensitive to inclusion of several groups of taxa, but it is less affected by the choice of alignment and rate methods. The majority of topologies and comparative results from Approximately Unbiased tests provide support for positioning the Rickettsiales and the mitochondrial branch within a clade. This composite clade is a sister group to the abundant marine SAR11 clade (Pelagibacterales). Furthermore, we add support for taxonomic assignment of several recently sequenced taxa. Accordingly, we propose three subclasses within the Alphaproteobacteria: the Caulobacteridae, the Rickettsidae, and the Magnetococcidae.

The origin and ecological significance of multiple branches for histidine utilization in Pseudomonas aeruginosa PAO1

Pseudomonas proliferate in a wide spectrum of harsh and variable environments. In many of these environments, amino acids, such as histidine, are a valuable source of carbon, nitrogen and energy. Here, we demonstrate that the histidine uptake and utilization (hut) pathway of Pseudomonas aeruginosa PAO1 contains two branches from the intermediate formiminoglutamate to the product glutamate. Genetic analysis revealed that the four-step route is dispensable as long as the five-step route is present (and vice versa). Mutants with deletions of either the four-step (HutE) or five-step (HutFG) branches were competed against each other and the wild-type strain to test the hypothesis of ecological redundancy; that is, that the presence of two pathways confers no benefit beyond that delivered by the individual pathways. Fitness assays performed under several environmental conditions led us to reject this hypothesis; the four-step pathway can provide an advantage when histidine is the sole carbon source. An IclR-type regulator (HutR) was identified that regulates the four-step pathway. Comparison of sequenced genomes revealed that P.aeruginosa strains and P.fluorescens Pf-5 have branched hut pathways. Phylogenetic analyses suggests that the gene encoding formiminoglutamase (hutE) was acquired by horizontal gene transfer from a Ralstonia-like ancestor. Potential barriers to inter-species transfer of the hutRE module were explored by transferring it from P.aeruginosa PAO1 to P.fluorescens SBW25. Transfer of the operon conferred the ability to utilize histidine via the four-step pathway in a single step, but the fitness cost of acquiring this new operon was found to be environment dependent.

Eumalacostracan phylogeny and total evidence: limitations of the usual suspects

BACKGROUND

The phylogeny of Eumalacostraca (Crustacea) remains elusive, despite over a century of interest. Recent morphological and molecular phylogenies appear highly incongruent, but this has not been assessed quantitatively. Moreover, 18S rRNA trees show striking branch length differences between species, accompanied by a conspicuous clustering of taxa with similar branch lengths. Surprisingly, previous research found no rate heterogeneity. Hitherto, no phylogenetic analysis of all major eumalacostracan taxa (orders) has either combined evidence from multiple loci, or combined molecular and morphological evidence.

RESULTS

We combined evidence from four nuclear ribosomal and mitochondrial loci (18S rRNA, 28S rRNA, 16S rRNA, and cytochrome c oxidase subunit I) with a newly synthesized morphological dataset. We tested the homogeneity of data partitions, both in terms of character congruence and the topological congruence of inferred trees. We also performed Bayesian and parsimony analyses on separate and combined partitions, and tested the contribution of each partition. We tested for potential long-branch attraction (LBA) using taxon deletion experiments, and with relative rate tests. Additionally we searched for molecular polytomies (spurious clades). Lastly, we investigated the phylogenetic stability of taxa, and assessed their impact on inferred relationships over the whole tree. We detected significant conflict between data partitions, especially between morphology and molecules. We found significant rate heterogeneity between species for both the 18S rRNA and combined datasets, introducing the possibility of LBA. As a test case, we showed that LBA probably affected the position of Spelaeogriphacea in the combined molecular evidence analysis. We also demonstrated that several clades, including the previously reported and surprising clade of Amphipoda plus Spelaeogriphacea, are 'supported' by zero length branches. Furthermore we showed that different sets of taxa have the greatest impact upon the relationships within molecular versus morphological trees.

CONCLUSION

Rate heterogeneity and conflict between data partitions mean that existing molecular and morphological evidence is unable to resolve a well-supported eumalacostracan phylogeny. We believe that it will be necessary to look beyond the most commonly utilized sources of data (nuclear ribosomal and mitochondrial sequences) to obtain a robust tree in the future.