RIBFIND2: Identifying rigid bodies in protein and nucleic acid structures.

Molecular structures are often fitted into cryo-EM maps by flexible fitting. When this requires large conformational changes, identifying rigid bodies can help optimize the model-map fit. Tools for identifying rigid bodies in protein structures exist, however an equivalent for nucleic acid structures is lacking. With the increase in cryo-EM maps containing RNA and progress in RNA structure prediction, there is a need for such tools. We previously developed RIBFIND, a program for clustering protein secondary structures into rigid bodies. In RIBFIND2, this approach is extended to nucleic acid structures. RIBFIND2 can identify biologically relevant rigid bodies in important groups of complex RNA structures, capturing a wide range of dynamics, including large rigid-body movements. The usefulness of RIBFIND2-assigned rigid bodies in cryo-EM model refinement was demonstrated on three examples, with two conformations each: Group II Intron complexed IEP, Internal Ribosome Entry Site and the Processome, using cryo-EM maps at 2.7-5 Å resolution. A hierarchical refinement approach, performed on progressively smaller sets of RIBFIND2 rigid bodies, was clearly shown to have an advantage over classical all-atom refinement. RIBFIND2 is available via a web server with structure visualization and as a standalone tool.

Overview and applications of map and model validation tools in the CCP-EM software suite.

Cryogenic electron microscopy (cryo-EM) has recently been established as a powerful technique for solving macromolecular structures. Although the best resolutions achievable are improving, a significant majority of data are still resolved at resolutions worse than 3 Å, where it is non-trivial to build or fit atomic models. The map reconstructions and atomic models derived from the maps are also prone to errors accumulated through the different stages of data processing. Here, we highlight the need to evaluate both model geometry and fit to data at different resolutions. Assessment of cryo-EM structures from SARS-CoV-2 highlights a bias towards optimising the model geometry to agree with the most common conformations, compared to the agreement with data. We present the CoVal web service which provides multiple validation metrics to reflect the quality of atomic models derived from cryo-EM data of structures from SARS-CoV-2. We demonstrate that further refinement can lead to improvement of the agreement with data without the loss of geometric quality. We also discuss the recent CCP-EM developments aimed at addressing some of the current shortcomings.

ProCarbDB: a database of carbohydrate-binding proteins.

Carbohydrate-binding proteins play crucial roles across all organisms and viruses. The complexity of carbohydrate structures, together with inconsistencies in how their 3D structures are reported, has led to difficulties in characterizing the protein-carbohydrate interfaces. In order to better understand protein-carbohydrate interactions, we have developed an open-access database, ProCarbDB, which, unlike the Protein Data Bank (PDB), clearly distinguishes between the complete carbohydrate ligands and their monomeric units. ProCarbDB is a comprehensive database containing over 5200 3D X-ray crystal structures of protein-carbohydrate complexes. In ProCarbDB, the complete carbohydrate ligands are annotated and all their interactions are displayed. Users can also select any protein residue in the proximity of the ligand to inspect its interactions with the carbohydrate ligand and with other neighbouring protein residues. Where available, additional curated information on the binding affinity of the complex and the effects of mutations on the binding have also been provided in the database. We believe that ProCarbDB will be an invaluable resource for understanding protein-carbohydrate interfaces. The ProCarbDB web server is freely available at http://www.procarbdb.science/procarb.

Fragment-based discovery of a new class of inhibitors targeting mycobacterial tRNA modification.

Translational frameshift errors are often deleterious to the synthesis of functional proteins and could therefore be promoted therapeutically to kill bacteria. TrmD (tRNA-(N(1)G37) methyltransferase) is an essential tRNA modification enzyme in bacteria that prevents +1 errors in the reading frame during protein translation and represents an attractive potential target for the development of new antibiotics. Here, we describe the application of a structure-guided fragment-based drug discovery approach to the design of a new class of inhibitors against TrmD in Mycobacterium abscessus. Fragment library screening, followed by structure-guided chemical elaboration of hits, led to the rapid development of drug-like molecules with potent in vitro TrmD inhibitory activity. Several of these compounds exhibit activity against planktonic M. abscessus and M. tuberculosis as well as against intracellular M. abscessus and M. leprae, indicating their potential as the basis for a novel class of broad-spectrum mycobacterial drugs.

KMT2B-related disorders: expansion of the phenotypic spectrum and long-term efficacy of deep brain stimulation.

Heterozygous mutations in KMT2B are associated with an early-onset, progressive and often complex dystonia (DYT28). Key characteristics of typical disease include focal motor features at disease presentation, evolving through a caudocranial pattern into generalized dystonia, with prominent oromandibular, laryngeal and cervical involvement. Although KMT2B-related disease is emerging as one of the most common causes of early-onset genetic dystonia, much remains to be understood about the full spectrum of the disease. We describe a cohort of 53 patients with KMT2B mutations, with detailed delineation of their clinical phenotype and molecular genetic features. We report new disease presentations, including atypical patterns of dystonia evolution and a subgroup of patients with a non-dystonic neurodevelopmental phenotype. In addition to the previously reported systemic features, our study has identified co-morbidities, including the risk of status dystonicus, intrauterine growth retardation, and endocrinopathies. Analysis of this study cohort (n = 53) in tandem with published cases (n = 80) revealed that patients with chromosomal deletions and protein truncating variants had a significantly higher burden of systemic disease (with earlier onset of dystonia) than those with missense variants. Eighteen individuals had detailed longitudinal data available after insertion of deep brain stimulation for medically refractory dystonia. Median age at deep brain stimulation was 11.5 years (range: 4.5-37.0 years). Follow-up after deep brain stimulation ranged from 0.25 to 22 years. Significant improvement of motor function and disability (as assessed by the Burke Fahn Marsden's Dystonia Rating Scales, BFMDRS-M and BFMDRS-D) was evident at 6 months, 1 year and last follow-up (motor, P = 0.001, P = 0.004, and P = 0.012; disability, P = 0.009, P = 0.002 and P = 0.012). At 1 year post-deep brain stimulation, >50% of subjects showed BFMDRS-M and BFMDRS-D improvements of >30%. In the long-term deep brain stimulation cohort (deep brain stimulation inserted for >5 years, n = 8), improvement of >30% was maintained in 5/8 and 3/8 subjects for the BFMDRS-M and BFMDRS-D, respectively. The greatest BFMDRS-M improvements were observed for trunk (53.2%) and cervical (50.5%) dystonia, with less clinical impact on laryngeal dystonia. Improvements in gait dystonia decreased from 20.9% at 1 year to 16.2% at last assessment; no patient maintained a fully independent gait. Reduction of BFMDRS-D was maintained for swallowing (52.9%). Five patients developed mild parkinsonism following deep brain stimulation. KMT2B-related disease comprises an expanding continuum from infancy to adulthood, with early evidence of genotype-phenotype correlations. Except for laryngeal dysphonia, deep brain stimulation provides a significant improvement in quality of life and function with sustained clinical benefit depending on symptoms distribution.

Cryo-electron microscopy targets in CASP13: Overview and evaluation of results.

Structures of seven CASP13 targets were determined using cryo-electron microscopy (cryo-EM) technique with resolution between 3.0 and 4.0 Å. We provide an overview of the experimentally derived structures and describe results of the numerical evaluation of the submitted models. The evaluation is carried out by comparing coordinates of models to those of reference structures (CASP-style evaluation), as well as checking goodness-of-fit of modeled structures to the cryo-EM density maps. The performance of contributing research groups in the CASP-style evaluation is measured in terms of backbone accuracy, all-atom local geometry and similarity of inter-subunit interfaces. The results on the cryo-EM targets are compared with those on the whole set of eighty CASP13 targets. A posteriori refinement of the best models in their corresponding cryo-EM density maps resulted in structures that are very close to the reference structure, including some regions with better fit to the density.

Disruption of a Structurally Important Extracellular Element in the Glycine Receptor Leads to Decreased Synaptic Integration and Signaling Resulting in Severe Startle Disease.

Functional impairments or trafficking defects of inhibitory glycine receptors (GlyRs) have been linked to human hyperekplexia/startle disease and autism spectrum disorders. We found that a lack of synaptic integration of GlyRs, together with disrupted receptor function, is responsible for a lethal startle phenotype in a novel spontaneous mouse mutant shaky, caused by a missense mutation, Q177K, located in the extracellular ß8-ß9 loop of the GlyR α1 subunit. Recently, structural data provided evidence that the flexibility of the ß8-ß9 loop is crucial for conformational transitions during opening and closing of the ion channel and represents a novel allosteric binding site in Cys-loop receptors. We identified the underlying neuropathological mechanisms in male and female shaky mice through a combination of protein biochemistry, immunocytochemistry, and both in vivo and in vitro electrophysiology. Increased expression of the mutant GlyR α1Q177K subunit in vivo was not sufficient to compensate for a decrease in synaptic integration of α1Q177Kß GlyRs. The remaining synaptic heteromeric α1Q177Kß GlyRs had decreased current amplitudes with significantly faster decay times. This functional disruption reveals an important role for the GlyR α1 subunit ß8-ß9 loop in initiating rearrangements within the extracellular-transmembrane GlyR interface and that this structural element is vital for inhibitory GlyR function, signaling, and synaptic clustering.SIGNIFICANCE STATEMENT GlyR dysfunction underlies neuromotor deficits in startle disease and autism spectrum disorders. We describe an extracellular GlyR α1 subunit mutation (Q177K) in a novel mouse startle disease mutant shaky Structural data suggest that during signal transduction, large transitions of the ß8-ß9 loop occur in response to neurotransmitter binding. Disruption of the ß8-ß9 loop by the Q177K mutation results in a disruption of hydrogen bonds between Q177 and the ligand-binding residue R65. Functionally, the Q177K change resulted in decreased current amplitudes, altered desensitization decay time constants, and reduced GlyR clustering and synaptic strength. The GlyR ß8-ß9 loop is therefore an essential regulator of conformational rearrangements during ion channel opening and closing.

Identification of new allosteric sites and modulators of AChE through computational and experimental tools.

Allosteric sites on proteins are targeted for designing more selective inhibitors of enzyme activity and to discover new functions. Acetylcholinesterase (AChE), which is most widely known for the hydrolysis of the neurotransmitter acetylcholine, has a peripheral allosteric subsite responsible for amyloidosis in Alzheimer's disease through interaction with amyloid ß-peptide. However, AChE plays other non-hydrolytic functions. Here, we identify and characterise using computational tools two new allosteric sites in AChE, which have allowed us to identify allosteric inhibitors by virtual screening guided by structure-based and fragment hotspot strategies. The identified compounds were also screened for in vitro inhibition of AChE and three were observed to be active. Further experimental (kinetic) and computational (molecular dynamics) studies have been performed to verify the allosteric activity. These new compounds may be valuable pharmacological tools in the study of non-cholinergic functions of AChE.

Refinement of atomic models in high resolution EM reconstructions using Flex-EM and local assessment.

As the resolutions of Three Dimensional Electron Microscopic reconstructions of biological macromolecules are being improved, there is a need for better fitting and refinement methods at high resolutions and robust approaches for model assessment. Flex-EM/MODELLER has been used for flexible fitting of atomic models in intermediate-to-low resolution density maps of different biological systems. Here, we demonstrate the suitability of the method to successfully refine structures at higher resolutions (2.5-4.5Å) using both simulated and experimental data, including a newly processed map of Apo-GroEL. A hierarchical refinement protocol was adopted where the rigid body definitions are relaxed and atom displacement steps are reduced progressively at successive stages of refinement. For the assessment of local fit, we used the SMOC (segment-based Manders' overlap coefficient) score, while the model quality was checked using the Qmean score. Comparison of SMOC profiles at different stages of refinement helped in detecting regions that are poorly fitted. We also show how initial model errors can have significant impact on the goodness-of-fit. Finally, we discuss the implementation of Flex-EM in the CCP-EM software suite.

Structure-guided, target-based drug discovery - exploiting genome information from HIV to mycobacterial infections.

The use of protein crystallography in structure-guided drug discovery allows identification of potential inhibitor-binding sites and optimisation of interactions of hits and lead compounds with a target protein. An early example of this approach was the use of the structure of HIV protease in designing AIDS antivirals. More recently, use of structure-guided design with fragment-based drug discovery, which reduces the size of screening libraries by decreasing complexity, has improved ligand efficiency in drug design. Here, we discuss the use of structure-guided target identification and lead optimisation using fragment-based approaches in the development of new antimicrobials for mycobacterial infections.

DOCKSCORE: a webserver for ranking protein-protein docked poses.

BACKGROUND: Proteins interact with a variety of other molecules such as nucleic acids, small molecules and other proteins inside the cell. Structure-determination of protein-protein complexes is challenging due to several reasons such as the large molecular weights of these macromolecular complexes, their dynamic nature, difficulty in purification and sample preparation. Computational docking permits an early understanding of the feasibility and mode of protein-protein interactions. However, docking algorithms propose a number of solutions and it is a challenging task to select the native or near native pose(s) from this pool. DockScore is an objective scoring scheme that can be used to rank protein-protein docked poses. It considers several interface parameters, namely, surface area, evolutionary conservation, hydrophobicity, short contacts and spatial clustering at the interface for scoring. RESULTS: We have implemented DockScore in form of a webserver for its use by the scientific community. DockScore webserver can be employed, subsequent to docking, to perform scoring of the docked solutions, starting from multiple poses as inputs. The results, on scores and ranks for all the poses, can be downloaded as a csv file and graphical view of the interface of best ranking poses is possible. CONCLUSIONS: The webserver for DockScore is made freely available for the scientific community at: http://caps.ncbs.res.in/dockscore/ .

Genome sequencing of herb Tulsi (Ocimum tenuiflorum) unravels key genes behind its strong medicinal properties.

BACKGROUND: Krishna Tulsi, a member of Lamiaceae family, is a herb well known for its spiritual, religious and medicinal importance in India. The common name of this plant is 'Tulsi' (or 'Tulasi' or 'Thulasi') and is considered sacred by Hindus. We present the draft genome of Ocimum tenuiflurum L (subtype Krishna Tulsi) in this report. The paired-end and mate-pair sequence libraries were generated for the whole genome sequenced with the Illumina Hiseq 1000, resulting in an assembled genome of 374 Mb, with a genome coverage of 61 % (612 Mb estimated genome size). We have also studied transcriptomes (RNA-Seq) of two subtypes of O. tenuiflorum, Krishna and Rama Tulsi and report the relative expression of genes in both the varieties. RESULTS: The pathways leading to the production of medicinally-important specialized metabolites have been studied in detail, in relation to similar pathways in Arabidopsis thaliana and other plants. Expression levels of anthocyanin biosynthesis-related genes in leaf samples of Krishna Tulsi were observed to be relatively high, explaining the purple colouration of Krishna Tulsi leaves. The expression of six important genes identified from genome data were validated by performing q-RT-PCR in different tissues of five different species, which shows the high extent of urosolic acid-producing genes in young leaves of the Rama subtype. In addition, the presence of eugenol and ursolic acid, implied as potential drugs in the cure of many diseases including cancer was confirmed using mass spectrometry. CONCLUSIONS: The availability of the whole genome of O.tenuiflorum and our sequence analysis suggests that small amino acid changes at the functional sites of genes involved in metabolite synthesis pathways confer special medicinal properties to this herb.

Genome-wide survey of DNA-binding proteins in Arabidopsis thaliana: analysis of distribution and functions.

The interaction of proteins with their respective DNA targets is known to control many high-fidelity cellular processes. Performing a comprehensive survey of the sequenced genomes for DNA-binding proteins (DBPs) will help in understanding their distribution and the associated functions in a particular genome. Availability of fully sequenced genome of Arabidopsis thaliana enables the review of distribution of DBPs in this model plant genome. We used profiles of both structure and sequence-based DNA-binding families, derived from PDB and PFam databases, to perform the survey. This resulted in 4471 proteins, identified as DNA-binding in Arabidopsis genome, which are distributed across 300 different PFam families. Apart from several plant-specific DNA-binding families, certain RING fingers and leucine zippers also had high representation. Our search protocol helped to assign DNA-binding property to several proteins that were previously marked as unknown, putative or hypothetical in function. The distribution of Arabidopsis genes having a role in plant DNA repair were particularly studied and noted for their functional mapping. The functions observed to be overrepresented in the plant genome harbour DNA-3-methyladenine glycosylase activity, alkylbase DNA N-glycosylase activity and DNA-(apurinic or apyrimidinic site) lyase activity, suggesting their role in specialized functions such as gene regulation and DNA repair.

Sequence search and analysis of gene products containing RNA recognition motifs in the human genome.

BACKGROUND: Gene expression is tightly regulated at both transcriptional and post-transcriptional levels. RNA-binding proteins are involved in post-transcriptional gene regulation events. They are involved in a variety of functions such as splicing, alternative splicing, nuclear import and export of mRNA, RNA stability and translation. There are several well-characterized RNA-binding motifs present in a whole genome, such as RNA recognition motif (RRM), KH domain, zinc-fingers etc. In the present study, we have investigated human genome for the presence of RRM-containing gene products starting from RRM domains in the Pfam (Protein family database) repository. RESULTS: In Pfam, seven families are recorded to contain RRM-containing proteins. We studied these families for their taxonomic representation, sequence features (identity, length, phylogeny) and structural properties (mapping conservation on the structures). We then examined the presence of RRM-containing gene products in Homo sapiens genome and identified 928 RRM-containing gene products. These were studied for their predicted domain architectures, biological processes, involvement in pathways, disease relevance and disorder content. RRM domains were observed to occur multiple times in a single polypeptide. However, there are 56 other co-existing domains involved in different regulatory functions. Further, functional enrichment analysis revealed that RRM-containing gene products are mainly involved in biological functions such as mRNA splicing and its regulation. CONCLUSIONS: Our sequence analysis identified RRM-containing gene products in the human genome and provides insights into their domain architectures and biological functions. Since mRNA splicing and gene regulation are important in the cellular machinery, this analysis provides an early overview of genes that carry out these functions.

Clinical Phenotype in Individuals With Birk-Landau-Perez Syndrome Associated With Biallelic SLC30A9 Pathogenic Variants.

BACKGROUND AND OBJECTIVES: Birk-Landau-Perez syndrome is a genetic disorder caused by biallelic pathogenic variants in SLC30A9 presenting with a complex movement disorder, developmental regression, oculomotor abnormalities, and renal impairment. It has previously been reported in 2 families. We describe the clinical phenotype of 8 further individuals from 4 unrelated families with SLC30A9-related disease. METHOD: Following detailed clinical phenotyping, 1 family underwent research whole-genome sequencing (WGS), 1 research whole-exome sequencing, and 2 diagnostic WGS. Variants of interest were assessed for pathogenicity using in silico prediction tools, homology modeling, and, where relevant, sequencing of complementary DNA (cDNA) for splicing effect. RESULTS: In 2 unrelated families of Pakistani origin (1 consanguineous and 1 not), the same homozygous missense variant in SLC30A9 (c.1253G>T, p.Gly418Val) was identified. Family 1 included 2 affected brothers, and family 2 one affected boy. In family 3, also consanguineous, there were 4 affected siblings homozygous for the variant c.1049delCAG, pAla350del. The fourth family was nonconsanguineous: the 1 affected individual was compound heterozygous for c.1083dup, p.Val362Cysfs*5, and c.1413A>G, p.Ser471=. Despite phenotypic variability between the 4 families, all affected patients manifested with a progressive hyperkinetic movement disorder, associated with oculomotor apraxia and ptosis. None had evidence of severe renal impairment. For the novel missense variant, the conformation of the loop domain and packing of transmembrane helices are likely to be disrupted based on structure modeling. Its presence in 2 unrelated Pakistani families suggests a possible founder variant. For the synonymous variant p.Ser471=, an effect on splicing was confirmed through cDNA analysis. DISCUSSION: Pathogenic variants in SLC30A9 cause a progressive autosomal recessive neurologic syndrome associated with a complex hyperkinetic movement disorder. Our report highlights the expanding disease phenotype, which can present with a wider spectrum of severity than has previously been recognized.

PTPN11 Mosaicism Causes a Spectrum of Pigmentary and Vascular Neurocutaneous Disorders and Predisposes to Melanoma.

Phakomatosis pigmentovascularis is a diagnosis that denotes the coexistence of pigmentary and vascular birthmarks of specific types, accompanied by variable multisystem involvement, including CNS disease, asymmetrical growth, and a predisposition to malignancy. Using a tight phenotypic group and high-depth next-generation sequencing of affected tissues, we discover here clonal mosaic variants in gene PTPN11 encoding SHP2 phosphatase as a cause of phakomatosis pigmentovascularis type III or spilorosea. Within an individual, the same variant is found in distinct pigmentary and vascular birthmarks and is undetectable in blood. We go on to show that the same variants can cause either the pigmentary or vascular phenotypes alone, and drive melanoma development within pigmentary lesions. Protein structure modeling highlights that although variants lead to loss of function at the level of the phosphatase domain, resultant conformational changes promote longer ligand binding. In vitro modeling of the missense variants confirms downstream MAPK pathway overactivation and widespread disruption of human endothelial cell angiogenesis. Importantly, patients with PTPN11 mosaicism theoretically risk passing on the variant to their children as the germline RASopathy Noonan syndrome with lentigines. These findings improve our understanding of the pathogenesis and biology of nevus spilus and capillary malformation syndromes, paving the way for better clinical management.

Re-visiting protein-centric two-tier classification of existing DNA-protein complexes.

BACKGROUND: Precise DNA-protein interactions play most important and vital role in maintaining the normal physiological functioning of the cell, as it controls many high fidelity cellular processes. Detailed study of the nature of these interactions has paved the way for understanding the mechanisms behind the biological processes in which they are involved. Earlier in 2000, a systematic classification of DNA-protein complexes based on the structural analysis of the proteins was proposed at two tiers, namely groups and families. With the advancement in the number and resolution of structures of DNA-protein complexes deposited in the Protein Data Bank, it is important to revisit the existing classification. RESULTS: On the basis of the sequence analysis of DNA binding proteins, we have built upon the protein centric, two-tier classification of DNA-protein complexes by adding new members to existing families and making new families and groups. While classifying the new complexes, we also realised the emergence of new groups and families. The new group observed was where ß-propeller was seen to interact with DNA. There were 34 SCOP folds which were observed to be present in the complexes of both old and new classifications, whereas 28 folds are present exclusively in the new complexes. Some new families noticed were NarL transcription factor, Z-α DNA binding proteins, Forkhead transcription factor, AP2 protein, Methyl CpG binding protein etc. CONCLUSIONS: Our results suggest that with the increasing number of availability of DNA-protein complexes in Protein Data Bank, the number of families in the classification increased by approximately three fold. The folds present exclusively in newly classified complexes is suggestive of inclusion of proteins with new function in new classification, the most populated of which are the folds responsible for DNA damage repair. The proposed re-visited classification can be used to perform genome-wide surveys in the genomes of interest for the presence of DNA-binding proteins. Further analysis of these complexes can aid in developing algorithms for identifying DNA-binding proteins and their family members from mere sequence information.

Atomic model validation using the CCP-EM software suite.

Recently, there has been a dramatic improvement in the quality and quantity of data derived using cryogenic electron microscopy (cryo-EM). This is also associated with a large increase in the number of atomic models built. Although the best resolutions that are achievable are improving, often the local resolution is variable, and a significant majority of data are still resolved at resolutions worse than 3 Å. Model building and refinement is often challenging at these resolutions, and hence atomic model validation becomes even more crucial to identify less reliable regions of the model. Here, a graphical user interface for atomic model validation, implemented in the CCP-EM software suite, is presented. It is aimed to develop this into a platform where users can access multiple complementary validation metrics that work across a range of resolutions and obtain a summary of evaluations. Based on the validation estimates from atomic models associated with cryo-EM structures from SARS-CoV-2, it was observed that models typically favor adopting the most common conformations over fitting the observations when compared with the model agreement with data. At low resolutions, the stereochemical quality may be favored over data fit, but care should be taken to ensure that the model agrees with the data in terms of resolvable features. It is demonstrated that further re-refinement can lead to improvement of the agreement with data without the loss of geometric quality. This also highlights the need for improved resolution-dependent weight optimization in model refinement and an effective test for overfitting that would help to guide the refinement process.

The pore conformation of lymphocyte perforin.

Perforin is a pore-forming protein that facilitates rapid killing of pathogen-infected or cancerous cells by the immune system. Perforin is released from cytotoxic lymphocytes, together with proapoptotic granzymes, to bind to a target cell membrane where it oligomerizes and forms pores. The pores allow granzyme entry, which rapidly triggers the apoptotic death of the target cell. Here, we present a 4-Å resolution cryo-electron microscopy structure of the perforin pore, revealing previously unidentified inter- and intramolecular interactions stabilizing the assembly. During pore formation, the helix-turn-helix motif moves away from the bend in the central ß sheet to form an intermolecular contact. Cryo-electron tomography shows that prepores form on the membrane surface with minimal conformational changes. Our findings suggest the sequence of conformational changes underlying oligomerization and membrane insertion, and explain how several pathogenic mutations affect function.

Mycobacterium abscessus pathogenesis identified by phenogenomic analyses.

The medical and scientific response to emerging and established pathogens is often severely hampered by ignorance of the genetic determinants of virulence, drug resistance and clinical outcomes that could be used to identify therapeutic drug targets and forecast patient trajectories. Taking the newly emergent multidrug-resistant bacteria Mycobacterium abscessus as an example, we show that combining high-dimensional phenotyping with whole-genome sequencing in a phenogenomic analysis can rapidly reveal actionable systems-level insights into bacterial pathobiology. Through phenotyping of 331 clinical isolates, we discovered three distinct clusters of isolates, each with different virulence traits and associated with a different clinical outcome. We combined genome-wide association studies with proteome-wide computational structural modelling to define likely causal variants, and employed direct coupling analysis to identify co-evolving, and therefore potentially epistatic, gene networks. We then used in vivo CRISPR-based silencing to validate our findings and discover clinically relevant M. abscessus virulence factors including a secretion system, thus illustrating how phenogenomics can reveal critical pathways within emerging pathogenic bacteria.