Predicted mechanistic impacts of human protein missense variants.

Genome sequencing efforts have led to the discovery of tens of millions of protein missense variants found in the human population with the majority of these having no annotated role and some likely contributing to trait variation and disease. Sequence-based artificial intelligence approaches have become highly accurate at predicting variants that are detrimental to the function of proteins but they do not inform on mechanisms of disruption. Here we combined sequence and structure-based methods to perform proteome-wide prediction of deleterious variants with information on their impact on protein stability, protein-protein interactions and small-molecule binding pockets. AlphaFold2 structures were used to predict approximately 100,000 small-molecule binding pockets and stability changes for over 200 million variants. To inform on protein-protein interfaces we used AlphaFold2 to predict structures for nearly 500,000 protein complexes. We illustrate the value of mechanism-aware variant effect predictions to study the relation between protein stability and abundance and the structural properties of interfaces underlying trans protein quantitative trait loci (pQTLs). We characterised the distribution of mechanistic impacts of protein variants found in patients and experimentally studied example disease linked variants in FGFR1.

Domain shuffling of a highly mutable ligand-binding fold drives adhesin generation across the bacterial kingdom.

Bacterial fibrillar adhesins are specialized extracellular polypeptides that promote the attachment of bacteria to the surfaces of other cells or materials. Adhesin-mediated interactions are critical for the establishment and persistence of stable bacterial populations within diverse environmental niches and are important determinants of virulence. The fibronectin (Fn)-binding fibrillar adhesin CshA, and its paralogue CshB, play important roles in host colonization by the oral commensal and opportunistic pathogen Streptococcus gordonii. As paralogues are often catalysts for functional diversification, we have probed the early stages of structural and functional divergence in Csh proteins by determining the X-ray crystal structure of the CshB adhesive domain NR2 and characterizing its Fn-binding properties in vitro. Despite sharing a common fold, CshB_NR2 displays an ~1.7-fold reduction in Fn-binding affinity relative to CshA_NR2. This correlates with reduced electrostatic charge in the Fn-binding cleft. Complementary bioinformatic studies reveal that homologues of CshA/B_NR2 domains are widely distributed in both Gram-positive and Gram-negative bacteria, where they are found housed within functionally cryptic multi-domain polypeptides. Our findings are consistent with the classification of Csh adhesins and their relatives as members of the recently defined polymer adhesin domain (PAD) family of bacterial proteins.

Cryo-EM structures of human RNA polymerase I.

RNA polymerase I (Pol I) specifically synthesizes ribosomal RNA. Pol I upregulation is linked to cancer, while mutations in the Pol I machinery lead to developmental disorders. Here we report the cryo-EM structure of elongating human Pol I at 2.7 Å resolution. In the exit tunnel, we observe a double-stranded RNA helix that may support Pol I processivity. Our structure confirms that human Pol I consists of 13 subunits with only one subunit forming the Pol I stalk. Additionally, the structure of human Pol I in complex with the initiation factor RRN3 at 3.1 Å resolution reveals stalk flipping upon RRN3 binding. We also observe an inactivated state of human Pol I bound to an open DNA scaffold at 3.3 Å resolution. Lastly, the high-resolution structure of human Pol I allows mapping of disease-related mutations that can aid understanding of disease etiology.

Discovery of fibrillar adhesins across bacterial species.

BACKGROUND: Fibrillar adhesins are long multidomain proteins that form filamentous structures at the cell surface of bacteria. They are an important yet understudied class of proteins composed of adhesive and stalk domains that mediate interactions of bacteria with their environment. This study aims to characterize fibrillar adhesins in a wide range of bacterial phyla and to identify new fibrillar adhesin-like proteins to improve our understanding of host-bacteria interactions. RESULTS: Through careful literature and computational searches, we identified 82 stalk and 27 adhesive domain families in fibrillar adhesins. Based on the presence of these domains in the UniProt Reference Proteomes database, we identified and analysed 3,542 fibrillar adhesin-like proteins across species of the most common bacterial phyla. We further enumerate the adhesive and stalk domain combinations found in nature and demonstrate that fibrillar adhesins have complex and variable domain architectures, which differ across species. By analysing the domain architecture of fibrillar adhesins, we show that in Gram positive bacteria, adhesive domains are mostly positioned at the N-terminus and cell surface anchors at the C-terminus of the protein, while their positions are more variable in Gram negative bacteria. We provide an open repository of fibrillar adhesin-like proteins and domains to enable further studies of this class of bacterial surface proteins. CONCLUSION: This study provides a domain-based characterization of fibrillar adhesins and demonstrates that they are widely found in species across the main bacterial phyla. We have discovered numerous novel fibrillar adhesins and improved our understanding of pathogenic adhesion and invasion mechanisms.

Periscope Proteins are variable-length regulators of bacterial cell surface interactions.

Changes at the cell surface enable bacteria to survive in dynamic environments, such as diverse niches of the human host. Here, we reveal "Periscope Proteins" as a widespread mechanism of bacterial surface alteration mediated through protein length variation. Tandem arrays of highly similar folded domains can form an elongated rod-like structure; thus, variation in the number of domains determines how far an N-terminal host ligand binding domain projects from the cell surface. Supported by newly available long-read genome sequencing data, we propose that this class could contain over 50 distinct proteins, including those implicated in host colonization and biofilm formation by human pathogens. In large multidomain proteins, sequence divergence between adjacent domains appears to reduce interdomain misfolding. Periscope Proteins break this "rule," suggesting that their length variability plays an important role in regulating bacterial interactions with host surfaces, other bacteria, and the immune system.

Cryo-EM structures of human RNA polymerase III in its unbound and transcribing states.

RNA polymerase III (Pol III) synthesizes transfer RNAs and other short, essential RNAs. Human Pol III misregulation is linked to tumor transformation, neurodegenerative and developmental disorders, and increased sensitivity to viral infections. Here, we present cryo-electron microscopy structures at 2.8 to 3.3 Å resolution of transcribing and unbound human Pol III. We observe insertion of the TFIIS-like subunit RPC10 into the polymerase funnel, providing insights into how RPC10 triggers transcription termination. Our structures resolve elements absent from Saccharomyces cerevisiae Pol III such as the winged-helix domains of RPC5 and an iron-sulfur cluster, which tethers the heterotrimer subcomplex to the core. The cancer-associated RPC7α isoform binds the polymerase clamp, potentially interfering with Pol III inhibition by tumor suppressor MAF1, which may explain why overexpressed RPC7α enhances tumor transformation. Finally, the human Pol III structure allows mapping of disease-related mutations and may contribute to the development of inhibitors that selectively target Pol III for therapeutic interventions.

The thrombospondin module 1 domain of the matricellular protein CCN3 shows an atypical disulfide pattern and incomplete CWR layers.

The members of the CCN (Cyr61/CTGF/Nov) family are a group of matricellular regulatory proteins that are essential to a wide range of functional pathways in cell signalling. Through interacting with extracellular matrix components and growth factors via one of their four domains, the CCN proteins are involved in critical biological processes such as angiogenesis, cell proliferation, bone development, fibrogenesis and tumorigenesis. Here, the crystal structure of the thrombospondin module 1 (TSP1) domain of CCN3 (previously known as Nov) is presented, which shares a similar three-stranded fold with the thrombospondin type 1 repeats of thrombospondin-1 and spondin-1, but with variations in the disulfide connectivity. Moreover, the CCN3 TSP1 domain lacks the typical π-stacked ladder of charged and aromatic residues on one side of the domain that is seen in other TSP1 domains. Using conservation analysis among orthologous domains, it is shown that a charged cluster in the centre of the domain is the most conserved site and this cluster is predicted to be a potential functional epitope for heparan sulfate binding. This variant TSP1 domain has also been used to revise the sequence determinants of TSP1 domains and to derive improved Pfam sequence profiles for the identification of novel TSP1 domains in more than 10 000 proteins across diverse phyla.

Defining the remarkable structural malleability of a bacterial surface protein Rib domain implicated in infection.

Streptococcus groups A and B cause serious infections, including early onset sepsis and meningitis in newborns. Rib domain-containing surface proteins are found associated with invasive strains and elicit protective immunity in animal models. Yet, despite their apparent importance in infection, the structure of the Rib domain was previously unknown. Structures of single Rib domains of differing length reveal a rare case of domain atrophy through deletion of 2 core antiparallel strands, resulting in the loss of an entire sheet of the ß-sandwich from an immunoglobulin-like fold. Previously, observed variation in the number of Rib domains within these bacterial cell wall-attached proteins has been suggested as a mechanism of immune evasion. Here, the structure of tandem domains, combined with molecular dynamics simulations and small angle X-ray scattering, suggests that variability in Rib domain number would result in differential projection of an N-terminal host-colonization domain from the bacterial surface. The identification of 2 further structures where the typical B-D-E immunoglobulin ß-sheet is replaced with an α-helix further confirms the extensive structural malleability of the Rib domain.

Assessment of protein assembly prediction in CASP13.

We present the assembly category assessment in the 13th edition of the CASP community-wide experiment. For the second time, protein assemblies constitute an independent assessment category. Compared to the last edition we see a clear uptake in participation, more oligomeric targets released, and consistent, albeit modest, improvement of the predictions quality. Looking at the tertiary structure predictions, we observe that ignoring the oligomeric state of the targets hinders modeling success. We also note that some contact prediction groups successfully predicted homomeric interfacial contacts, though it appears that these predictions were not used for assembly modeling. Homology modeling with sizeable human intervention appears to form the basis of the assembly prediction techniques in this round of CASP. Future developments should see more integrated approaches where subunits are modeled in the context of the assemblies they form.

Tandem domain swapping: determinants of multidomain protein misfolding.

Tandem homologous domains in proteins are susceptible to misfolding through the formation of domain swaps, non-native conformations involving the exchange of equivalent structural elements between adjacent domains. Cutting-edge biophysical experiments have recently allowed the observation of tandem domain swapping events at the single molecule level. In addition, computer simulations have shed light into the molecular mechanisms of domain swap formation and serve as the basis for methods to systematically predict them. At present, the number of studies on tandem domain swaps is still small and limited to a few domain folds, but they offer important insights into the folding and evolution of multidomain proteins with applications in the field of protein design.

Analyzing the symmetrical arrangement of structural repeats in proteins with CE-Symm.

Many proteins fold into highly regular and repetitive three dimensional structures. The analysis of structural patterns and repeated elements is fundamental to understand protein function and evolution. We present recent improvements to the CE-Symm tool for systematically detecting and analyzing the internal symmetry and structural repeats in proteins. In addition to the accurate detection of internal symmetry, the tool is now capable of i) reporting the type of symmetry, ii) identifying the smallest repeating unit, iii) describing the arrangement of repeats with transformation operations and symmetry axes, and iv) comparing the similarity of all the internal repeats at the residue level. CE-Symm 2.0 helps the user investigate proteins with a robust and intuitive sequence-to-structure analysis, with many applications in protein classification, functional annotation and evolutionary studies. We describe the algorithmic extensions of the method and demonstrate its applications to the study of interesting cases of protein evolution.

BioJava 5: A community driven open-source bioinformatics library.

BioJava is an open-source project that provides a Java library for processing biological data. The project aims to simplify bioinformatic analyses by implementing parsers, data structures, and algorithms for common tasks in genomics, structural biology, ontologies, phylogenetics, and more. Since 2012, we have released two major versions of the library (4 and 5) that include many new features to tackle challenges with increasingly complex macromolecular structure data. BioJava requires Java 8 or higher and is freely available under the LGPL 2.1 license. The project is hosted on GitHub at https://github.com/biojava/biojava. More information and documentation can be found online on the BioJava website (http://www.biojava.org) and tutorial (https://github.com/biojava/biojava-tutorial). All inquiries should be directed to the GitHub page or the BioJava mailing list (http://lists.open-bio.org/mailman/listinfo/biojava-l).

TADOSS: computational estimation of tandem domain swap stability.

SUMMARY: Proteins with highly similar tandem domains have shown an increased propensity for misfolding and aggregation. Several molecular explanations have been put forward, such as swapping of adjacent domains, but there is a lack of computational tools to systematically analyze them. We present the TAndem DOmain Swap Stability predictor (TADOSS), a method to computationally estimate the stability of tandem domain-swapped conformations from the structures of single domains, based on previous coarse-grained simulation studies. The tool is able to discriminate domains susceptible to domain swapping and to identify structural regions with high propensity to form hinge loops. TADOSS is a scalable method and suitable for large scale analyses. AVAILABILITY AND IMPLEMENTATION: Source code and documentation are freely available under an MIT license on GitHub at https://github.com/lafita/tadoss. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Automated evaluation of quaternary structures from protein crystals.

A correct assessment of the quaternary structure of proteins is a fundamental prerequisite to understanding their function, physico-chemical properties and mode of interaction with other proteins. Currently about 90% of structures in the Protein Data Bank are crystal structures, in which the correct quaternary structure is embedded in the crystal lattice among a number of crystal contacts. Computational methods are required to 1) classify all protein-protein contacts in crystal lattices as biologically relevant or crystal contacts and 2) provide an assessment of how the biologically relevant interfaces combine into a biological assembly. In our previous work we addressed the first problem with our EPPIC (Evolutionary Protein Protein Interface Classifier) method. Here, we present our solution to the second problem with a new method that combines the interface classification results with symmetry and topology considerations. The new algorithm enumerates all possible valid assemblies within the crystal using a graph representation of the lattice and predicts the most probable biological unit based on the pairwise interface scoring. Our method achieves 85% precision (ranging from 76% to 90% for different oligomeric types) on a new dataset of 1,481 biological assemblies with consensus of PDB annotations. Although almost the same precision is achieved by PISA, currently the most popular quaternary structure assignment method, we show that, due to the fundamentally different approach to the problem, the two methods are complementary and could be combined to improve biological assembly assignments. The software for the automatic assessment of protein assemblies (EPPIC version 3) has been made available through a web server at http://www.eppic-web.org.

Assessment of protein assembly prediction in CASP12.

We present the results of the first independent assessment of protein assemblies in CASP. A total of 1624 oligomeric models were submitted by 108 predictor groups for the 30 oligomeric targets in the CASP12 edition. We evaluated the accuracy of oligomeric predictions by comparison to their reference structures at the interface patch and residue contact levels. We find that interface patches are more reliably predicted than the specific residue contacts. Whereas none of the 15 hard oligomeric targets have successful predictions for the residue contacts at the interface, six have models with resemblance in the interface patch. Successful predictions of interface patch and contacts exist for all targets suitable for homology modeling, with at least one group improving over the best available template for each target. However, the participation in protein assembly prediction is low and uneven. Three human groups are closely ranked at the top by overall performance, but a server outperforms all other predictors for targets suitable for homology modeling. The state of the art of protein assembly prediction methods is in development and has apparent room for improvement, especially for assemblies without templates.

Biological and functional relevance of CASP predictions.

Our goal is to answer the question: compared with experimental structures, how useful are predicted models for functional annotation? We assessed the functional utility of predicted models by comparing the performances of a suite of methods for functional characterization on the predictions and the experimental structures. We identified 28 sites in 25 protein targets to perform functional assessment. These 28 sites included nine sites with known ligand binding (holo-sites), nine sites that are expected or suggested by experimental authors for small molecule binding (apo-sites), and Ten sites containing important motifs, loops, or key residues with important disease-associated mutations. We evaluated the utility of the predictions by comparing their microenvironments to the experimental structures. Overall structural quality correlates with functional utility. However, the best-ranked predictions (global) may not have the best functional quality (local). Our assessment provides an ability to discriminate between predictions with high structural quality. When assessing ligand-binding sites, most prediction methods have higher performance on apo-sites than holo-sites. Some servers show consistently high performance for certain types of functional sites. Finally, many functional sites are associated with protein-protein interaction. We also analyzed biologically relevant features from the protein assemblies of two targets where the active site spanned the protein-protein interface. For the assembly targets, we find that the features in the models are mainly determined by the choice of template.