Protein complexes are under evolutionary selection to assemble via ordered pathways

Altered assembly paths mitigate interference among paralogous complexes

Protein complexes are fundamental to all cellular processes, so understanding their evolutionary history and assembly processes is important. Gene duplication followed by divergence is considered a primary mechanism for diversifying protein complexes. Nonetheless, to what extent assembly of present-day paralogous complexes has been constrained by their long evolutionary pathways and how cross-complex interference is avoided remain unanswered questions. Subunits of protein complexes are often stabilized upon complex formation, whereas unincorporated subunits are degraded. How such cooperative stability influences protein complex assembly also remains unclear. Here, we demonstrate that subcomplexes determined by cooperative stabilization interactions serve as building blocks for protein complex assembly. We further develop a protein stability-guided method to compare the assembly processes of paralogous complexes in cellulo. Our findings support that oligomeric state and the structural organization of paralogous complexes can be maintained even if their assembly processes are rearranged. Our results indicate that divergent assembly processes by paralogous complexes not only enable the complexes to evolve new functions, but also reinforce their segregation by establishing incompatibility against deleterious hybrid assemblies.

Hallmarks and evolutionary drivers of cotranslational protein complex assembly

Recent discoveries have highlighted the prevalence of cotranslational assembly in proteomes, revealing a range of mechanisms that enables the assembly of protein complex subunits on the ribosome. Structural analyses have uncovered emergent properties that may inherently control whether a subunit undergoes cotranslational assembly. However, the evolutionary paths that have yielded such complexes over an extended timescale remain largely unclear. In this review, we reflect on historical experiments that contributed to the field, including breakthroughs that have made possible the proteome-wide detection of cotranslational assembly, and the technical challenges yet to be overcome. We introduce a simple framework that encapsulates the hallmarks of cotranslational assembly and discuss how results from new experiments are shaping our view of the mechanistic, structural and evolutionary factors driving the phenomenon.

Imperfect symmetry facilitated the evolution of specificity and high-order stoichiometry in vertebrate hemoglobin

Many proteins form paralogous multimers - molecular complexes in which evolutionarily related proteins are arranged into specific quaternary structures. Little is known about the mechanisms by which they acquired their stoichiometry (the number of total subunits in the complex) and heterospecificity (the preference of subunits for their paralogs rather than other copies of the same protein). Here we use ancestral protein reconstruction and biochemical experiments to study historical increases in stoichiometry and specificity during the evolution of vertebrate hemoglobin (Hb), a α2β2 heterotetramer that evolved from a homodimeric ancestor after a gene duplication. We show that the mechanisms for this evolutionary transition was simple. One hydrophobic substitution in subunit β after the gene duplication was sufficient to cause the ancestral dimer to homotetramerize with high affinity across a new interface. During this same interval, a single-residue deletion in subunit α at the older interface conferred specificity for the heterotetrameric form and the trans-orientation of subunits within it. These sudden transitions in stoichiometry and specificity were possible because the interfaces in Hb are isologous - involving the same surface patch on interacting subunits, rotated 180° relative to each other - but the symmetry is slightly imperfect. This architecture amplifies the impacts of individual mutations on stoichiometry and specificity, especially in higher-order complexes, and allows single substitutions to differentially affect heteromeric vs homomeric interactions. Many multimers are isologous, and symmetry in proteins is always imperfect; our findings therefore suggest that elaborate and specific molecular complexes may often evolve via simple genetic and physical mechanisms.

Are there double knots in proteins? Prediction and in vitro verification based on TrmD-Tm1570 fusion from C. nitroreducens

We have been aware of the existence of knotted proteins for over 30 years-but it is hard to predict what is the most complicated knot that can be formed in proteins. Here, we show new and the most complex knotted topologies recorded to date-double trefoil knots (31 #31). We found five domain arrangements (architectures) that result in a doubly knotted structure in almost a thousand proteins. The double knot topology is found in knotted membrane proteins from the CaCA family, that function as ion transporters, in the group of carbonic anhydrases that catalyze the hydration of carbon dioxide, and in the proteins from the SPOUT superfamily that gathers 31 knotted methyltransferases with the active site-forming knot. For each family, we predict the presence of a double knot using AlphaFold and RoseTTaFold structure prediction. In the case of the TrmD-Tm1570 protein, which is a member of SPOUT superfamily, we show that it folds in vitro and is biologically active. Our results show that this protein forms a homodimeric structure and retains the ability to modify tRNA, which is the function of the single-domain TrmD protein. However, how the protein folds and is degraded remains unknown.

Discovering optimal kinetic pathways for self-assembly using automatic differentiation

Macromolecular complexes are often composed of diverse subunits. The self-assembly of these subunits is inherently nonequilibrium and must avoid kinetic traps to achieve high yield over feasible timescales. We show how the kinetics of self-assembly benefits from diversity in subunits because it generates an expansive parameter space that naturally improves the "expressivity" of self-assembly, much like a deeper neural network. By using automatic differentiation algorithms commonly used in deep learning, we searched the parameter spaces of mass-action kinetic models to identify classes of kinetic protocols that mimic biological solutions for productive self-assembly. Our results reveal how high-yield complexes that easily become kinetically trapped in incomplete intermediates can instead be steered by internal design of rate-constants or external and active control of subunits to efficiently assemble. Internal design of a hierarchy of subunit binding rates generates self-assembly that can robustly avoid kinetic traps for all concentrations and energetics, but it places strict constraints on selection of relative rates. External control via subunit titration is more versatile, avoiding kinetic traps for any system without requiring molecular engineering of binding rates, albeit less efficiently and robustly. We derive theoretical expressions for the timescales of kinetic traps, and we demonstrate our optimization method applies not just for design but inference, extracting intersubunit binding rates from observations of yield-vs.-time for a heterotetramer. Overall, we identify optimal kinetic protocols for self-assembly as a powerful mechanism to achieve efficient and high-yield assembly in synthetic systems whether robustness or ease of "designability" is preferred.

Co-translational Assembly Pathways of Nuclear Multiprotein Complexes Involved in the Regulation of Gene Transcription

Most factors that regulate gene transcription in eukaryotic cells are multimeric, often large, protein complexes. The understanding of the biogenesis pathways of such large and heterogeneous protein assemblies, as well as the dimerization partner choice among transcription factors, is crucial to interpret and control gene expression programs and consequent cell fate decisions. Co-translational assembly (Co-TA) is thought to play key roles in the biogenesis of protein complexes by directing complex formation during protein synthesis. In this review we discuss the principles of Co-TA with a special focus for the assembly of transcription regulatory complexes. We outline the expected molecular advantages of establishing co-translational interactions, pointing at the available, or missing, evidence for each of them. We hypothesize different molecular mechanisms based on Co-TA to explain the allocation "dilemma" of paralog proteins and subunits shared by different transcription complexes. By taking as a paradigm the different assembly pathways employed by three related transcription regulatory complexes (TFIID, SAGA and ATAC), we discuss alternative Co-TA strategies for nuclear multiprotein complexes and the widespread - yet specific - use of Co-TA for the formation of nuclear complexes involved in gene transcription. Ultimately, we outlined a series of open questions which demand well-defined lines of research to investigate the principles of gene regulation that rely on the coordinated assembly of protein complexes.

An atlas of protein homo-oligomerization across domains of life

Protein structures are essential to understanding cellular processes in molecular detail. While advances in artificial intelligence revealed the tertiary structure of proteins at scale, their quaternary structure remains mostly unknown. We devise a scalable strategy based on AlphaFold2 to predict homo-oligomeric assemblies across four proteomes spanning the tree of life. Our results suggest that approximately 45% of an archaeal proteome and a bacterial proteome and 20% of two eukaryotic proteomes form homomers. Our predictions accurately capture protein homo-oligomerization, recapitulate megadalton complexes, and unveil hundreds of homo-oligomer types, including three confirmed experimentally by structure determination. Integrating these datasets with omics information suggests that a majority of known protein complexes are symmetric. Finally, these datasets provide a structural context for interpreting disease mutations and reveal coiled-coil regions as major enablers of quaternary structure evolution in human. Our strategy is applicable to any organism and provides a comprehensive view of homo-oligomerization in proteomes.

DNA-Nanostructure-Guided Assembly of Proteins into Programmable Shapes

The development of methods to synthesize artificial protein complexes with precisely controlled configurations will enable diverse biological and medical applications. Using DNA to link proteins provides programmability that can be difficult to achieve with other methods. Here, we use DNA origami as an "assembler" to guide the linking of protein-DNA conjugates using a series of oligonucleotide hybridization and displacement operations. We constructed several isomeric protein nanostructures, including a dimer, two types of trimer structures, and three types of tetramer assemblies, on a DNA origami platform by using a C3-symmetric building block composed of a protein trimer modified with DNA handles. Our approach expands the scope for the precise assembly of protein-based nanostructures and will enable the formulation of functional protein complexes with stoichiometric and geometric control.

Structure and dynamics of endogenous cardiac troponin complex in human heart tissue captured by native nanoproteomics

Protein complexes are highly dynamic entities that display substantial diversity in their assembly, post-translational modifications, and non-covalent interactions, allowing them to play critical roles in various biological processes. The heterogeneity, dynamic nature, and low abundance of protein complexes in their native states present challenges to study using conventional structural biology techniques. Here we develop a native nanoproteomics strategy for the enrichment and subsequent native top-down mass spectrometry (nTDMS) analysis of endogenous cardiac troponin (cTn) complex directly from human heart tissue. The cTn complex is enriched and purified using peptide-functionalized superparamagnetic nanoparticles under non-denaturing conditions to enable the isotopic resolution of cTn complex, revealing their complex structure and assembly. Moreover, nTDMS elucidates the stoichiometry and composition of the cTn complex, localizes Ca2+ binding domains, defines cTn-Ca2+ binding dynamics, and provides high-resolution mapping of the proteoform landscape. This native nanoproteomics strategy opens a paradigm for structural characterization of endogenous native protein complexes.

New algorithms for structure informed genome rearrangement

We define two new computational problems in the domain of perfect genome rearrangements, and propose three algorithms to solve them. The rearrangement scenarios modeled by the problems consider Reversal and Block Interchange operations, and a PQ-tree is utilized to guide the allowed operations and to compute their weights. In the first problem, [Formula: see text] ([Formula: see text]), we define the basic structure-informed rearrangement measure. Here, we assume that the gene order members of the gene cluster from which the PQ-tree is constructed are permutations. The PQ-tree representing the gene cluster is ordered such that the series of gene IDs spelled by its leaves is equivalent to that of the reference gene order. Then, a structure-informed genome rearrangement distance is computed between the ordered PQ-tree and the target gene order. The second problem, [Formula: see text] ([Formula: see text]), generalizes [Formula: see text], where the gene order members are not necessarily permutations and the structure informed rearrangement measure is extended to also consider up to [Formula: see text] and [Formula: see text] gene insertion and deletion operations, respectively, when modelling the PQ-tree informed divergence process from the reference gene order to the target gene order. The first algorithm solves [Formula: see text] in [Formula: see text] time and [Formula: see text] space, where [Formula: see text] is the maximum number of children of a node, n is the length of the string and the number of leaves in the tree, and [Formula: see text] and [Formula: see text] are the number of P-nodes and Q-nodes in the tree, respectively. If one of the penalties of [Formula: see text] is 0, then the algorithm runs in [Formula: see text] time and [Formula: see text] space. The second algorithm solves [Formula: see text] in [Formula: see text] time and [Formula: see text] space, where [Formula: see text] is the maximum number of children of a node, n is the length of the string, m is the number of leaves in the tree, [Formula: see text] and [Formula: see text] are the number of P-nodes and Q-nodes in the tree, respectively, and allowing up to [Formula: see text] deletions from the tree and up to [Formula: see text] deletions from the string. The third algorithm is intended to reduce the space complexity of the second algorithm. It solves a variant of the problem (where one of the penalties of [Formula: see text] is 0) in [Formula: see text] time and [Formula: see text] space. The algorithm is implemented as a software tool, denoted MEM-Rearrange, and applied to the comparative and evolutionary analysis of 59 chromosomal gene clusters extracted from a dataset of 1487 prokaryotic genomes.

Unidirectional gene pairs in archaea and bacteria require overlaps or very short intergenic distances for translational coupling via termination-reinitiation and often encode subunits of heteromeric complexes

Genomes of bacteria and archaea contain a much larger fraction of unidirectional (serial) gene pairs than convergent or divergent gene pairs. Many of the unidirectional gene pairs have short overlaps of -4 nt and -1 nt. As shown previously, translation of the genes in overlapping unidirectional gene pairs is tightly coupled. Two alternative models for the fate of the post-termination ribosome predict either that overlaps or very short intergenic distances are essential for translational coupling or that the undissociated post-termination ribosome can scan through long intergenic regions, up to hundreds of nucleotides. We aimed to experimentally resolve the contradiction between the two models by analyzing three native gene pairs from the model archaeon Haloferax volcanii and three native pairs from Escherichia coli. A two reporter gene system was used to quantify the reinitiation frequency, and several stop codons in the upstream gene were introduced to increase the intergenic distances. For all six gene pairs from two species, an extremely strong dependence of the reinitiation efficiency on the intergenic distance was unequivocally demonstrated, such that even short intergenic distances of about 20 nt almost completely abolished translational coupling. Bioinformatic analysis of the intergenic distances in all unidirectional gene pairs in the genomes of H. volcanii and E. coli and in 1,695 prokaryotic species representative of 49 phyla showed that intergenic distances of -4 nt or -1 nt (= short gene overlaps of 4 nt or 1 nt) were by far most common in all these groups of archaea and bacteria. A small set of genes in E. coli, but not in H. volcanii, had intergenic distances of around +10 nt. Our experimental and bioinformatic analyses clearly show that translational coupling requires short gene overlaps, whereas scanning of intergenic regions by the post-termination ribosome occurs rarely, if at all. Short overlaps are enriched among genes that encode subunits of heteromeric complexes, and co-translational complex formation requiring precise subunit stoichiometry likely confers an evolutionary advantage that drove the formation and conservation of overlapping gene pairs during evolution.

Puzzling out nuclear pore complex assembly

Nuclear pore complexes (NPCs) are sophisticated multiprotein assemblies embedded within the nuclear envelope and controlling the exchanges of molecules between the cytoplasm and the nucleus. In this review, we summarize the mechanisms by which these elaborate complexes are built from their subunits, the nucleoporins, based on our ever-growing knowledge of NPC structural organization and on the recent identification of additional features of this process. We present the constraints faced during the production of nucleoporins, their gathering into oligomeric complexes, and the formation of NPCs within nuclear envelopes, and review the cellular strategies at play, from co-translational assembly to the enrolment of a panel of cofactors. Remarkably, the study of NPCs can inform our perception of the biogenesis of multiprotein complexes in general - and vice versa.

A hierarchical assembly pathway directs the unique subunit arrangement of TRiC/CCT

How the essential eukaryotic chaperonin TRiC/CCT assembles from eight distinct subunits into a unique double-ring architecture remains undefined. We show TRiC assembly involves a hierarchical pathway that segregates subunits with distinct functional properties until holocomplex (HC) completion. A stable, likely early intermediate arises from small oligomers containing CCT2, CCT4, CCT5, and CCT7, contiguous subunits that constitute the negatively charged hemisphere of the TRiC chamber, which has weak affinity for unfolded actin. The remaining subunits CCT8, CCT1, CCT3, and CCT6, which comprise the positively charged chamber hemisphere that binds unfolded actin more strongly, join the ring individually. Unincorporated late-assembling subunits are highly labile in cells, which prevents their accumulation and premature substrate binding. Recapitulation of assembly in a recombinant system demonstrates that the subunits in each hemisphere readily form stable, noncanonical TRiC-like HCs with aberrant functional properties. Thus, regulation of TRiC assembly along a biochemical axis disfavors the formation of stable alternative chaperonin complexes.

Discovering optimal kinetic pathways for self-assembly using automatic differentiation

During self-assembly of macromolecules ranging from ribosomes to viral capsids, the formation of long-lived intermediates or kinetic traps can dramatically reduce yield of the functional products. Understanding biological mechanisms for avoiding traps and efficiently assembling is essential for designing synthetic assembly systems, but learning optimal solutions requires numerical searches in high-dimensional parameter spaces. Here, we exploit powerful automatic differentiation algorithms commonly employed by deep learning frameworks to optimize physical models of reversible self-assembly, discovering diverse solutions in the space of rate constants for 3-7 subunit complexes. We define two biologically-inspired protocols that prevent kinetic trapping through either internal design of subunit binding kinetics or external design of subunit titration in time. Our third protocol acts to recycle intermediates, mimicking energy-consuming enzymes. Preventative solutions via interface design are the most efficient and scale better with more subunits, but external control via titration or recycling are effective even for poorly evolved binding kinetics. Whilst all protocols can produce good solutions, diverse subunits always helps; these complexes access more efficient solutions when following external control protocols, and are simpler to design for internal control, as molecular interfaces do not need modification during assembly given sufficient variation in dimerization rates. Our results identify universal scaling in the cost of kinetic trapping, and provide multiple strategies for eliminating trapping and maximizing assembly yield across large parameter spaces.

Fungal COP9 signalosome assembly requires connection of two trimeric intermediates for integration of intrinsic deneddylase

The conserved eight-subunit COP9 signalosome (CSN) is required for multicellular fungal development. The CSN deneddylase cooperates with the Cand1 exchange factor to control replacements of E3 ubiquitin cullin RING ligase receptors, providing specificity to eukaryotic protein degradation. Aspergillus nidulans CSN assembles through a heptameric pre-CSN, which is activated by integration of the catalytic CsnE deneddylase. Combined genetic and biochemical approaches provided the assembly choreography within a eukaryotic cell for native fungal CSN. Interactomes of functional GFP-Csn subunit fusions in pre-CSN deficient fungal strains were compared by affinity purifications and mass spectrometry. Two distinct heterotrimeric CSN subcomplexes were identified as pre-CSN assembly intermediates. CsnA-C-H and CsnD-F-G form independently of CsnB, which connects the heterotrimers to a heptamer and enables subsequent integration of CsnE to form the enzymatically active CSN complex. Surveillance mechanisms control accurate Csn subunit amounts and correct cellular localization for sequential assembly since deprivation of Csn subunits changes the abundance and location of remaining Csn subunits.

Wiggle and Shake: Managing and Exploiting Conformational Dynamics during Proteasome Biogenesis

The 26S proteasome is the largest and most complicated protease known, and changes to proteasome assembly or function contribute to numerous human diseases. Assembly of the 26S proteasome from its ~66 individual polypeptide subunits is a highly orchestrated process requiring the concerted actions of both intrinsic elements of proteasome subunits, as well as assistance by extrinsic, dedicated proteasome assembly chaperones. With the advent of near-atomic resolution cryo-electron microscopy, it has become evident that the proteasome is a highly dynamic machine, undergoing numerous conformational changes in response to ligand binding and during the proteolytic cycle. In contrast, an appreciation of the role of conformational dynamics during the biogenesis of the proteasome has only recently begun to emerge. Herein, we review our current knowledge of proteasome assembly, with a particular focus on how conformational dynamics guide particular proteasome biogenesis events. Furthermore, we highlight key emerging questions in this rapidly expanding area.

Structure and dynamics of endogenous protein complexes in human heart tissue captured by native nanoproteomics

Protein complexes are highly dynamic entities that display substantial diversity in their assembly, post-translational modifications, and non-covalent interactions, allowing them to play critical roles in various biological processes. The heterogeneity, dynamic nature, and low abundance of protein complexes in their native states present tremendous challenges to study using conventional structural biology techniques. Here we develop a "native nanoproteomics" strategy for the native enrichment and subsequent native top-down mass spectrometry (nTDMS) of low-abundance protein complexes. Specifically, we demonstrate the first comprehensive characterization of the structure and dynamics of cardiac troponin (cTn) complexes directly from human heart tissue. The endogenous cTn complex is effectively enriched and purified using peptide-functionalized superparamagnetic nanoparticles under non-denaturing conditions to enable the isotopic resolution of cTn complexes, revealing their complex structure and assembly. Moreover, nTDMS elucidates the stoichiometry and composition of the heterotrimeric cTn complex, localizes Ca2+ binding domains (II-IV), defines cTn-Ca2+ binding dynamics, and provides high-resolution mapping of the proteoform landscape. This native nanoproteomics strategy opens a new paradigm for structural characterization of low-abundance native protein complexes.

Structure and dynamics of endogenous protein complexes in human heart tissue captured by native nanoproteomics

Protein complexes are highly dynamic entities that display substantial diversity in their assembly, post-translational modifications, and non-covalent interactions, allowing them to play critical roles in various biological processes. The heterogeneity, dynamic nature, and low abundance of protein complexes in their native states present tremendous challenges to study using conventional structural biology techniques. Here we develop a "native nanoproteomics" strategy for the native enrichment and subsequent native top-down mass spectrometry (nTDMS) of low-abundance protein complexes. Specifically, we demonstrate the first comprehensive characterization of the structure and dynamics of cardiac troponin (cTn) complexes directly from human heart tissue. The endogenous cTn complex is effectively enriched and purified using peptide-functionalized superparamagnetic nanoparticles under non-denaturing conditions to enable the isotopic resolution of cTn complexes, revealing their complex structure and assembly. Moreover, nTDMS elucidates the stoichiometry and composition of the heterotrimeric cTn complex, localizes Ca2+ binding domains (II-IV), defines cTn-Ca2+ binding dynamics, and provides high-resolution mapping of the proteoform landscape. This native nanoproteomics strategy opens a new paradigm for structural characterization of low-abundance native protein complexes.

A highly contiguous genome assembly reveals sources of genomic novelty in the symbiotic fungus Rhizophagus irregularis

The root systems of most plant species are aided by the soil-foraging capacities of symbiotic arbuscular mycorrhizal (AM) fungi of the Glomeromycotina subphylum. Despite recent advances in our knowledge of the ecology and molecular biology of this mutualistic symbiosis, our understanding of the AM fungi genome biology is just emerging. Presented here is a close to T2T genome assembly of the model AM fungus Rhizophagus irregularis DAOM197198, achieved through Nanopore long-read DNA sequencing and Hi-C data. This haploid genome assembly of R. irregularis, alongside short- and long-read RNA-Sequencing data, was used to produce a comprehensive annotation catalog of gene models, repetitive elements, small RNA loci, and DNA cytosine methylome. A phylostratigraphic gene age inference framework revealed that the birth of genes associated with nutrient transporter activity and transmembrane ion transport systems predates the emergence of Glomeromycotina. While nutrient cycling in AM fungi relies on genes that existed in ancestor lineages, a burst of Glomeromycotina-restricted genetic innovation is also detected. Analysis of the chromosomal distribution of genetic and epigenetic features highlights evolutionarily young genomic regions that produce abundant small RNAs, suggesting active RNA-based monitoring of genetic sequences surrounding recently evolved genes. This chromosome-scale view of the genome of an AM fungus genome reveals previously unexplored sources of genomic novelty in an organism evolving under an obligate symbiotic life cycle.

Architecture and genomic arrangement of the MurE-MurF bacterial cell wall biosynthesis complex

Peptidoglycan (PG) is a central component of the bacterial cell wall, and the disruption of its biosynthetic pathway has been a successful antibacterial strategy for decades. PG biosynthesis is initiated in the cytoplasm through sequential reactions catalyzed by Mur enzymes that have been suggested to associate into a multimembered complex. This idea is supported by the observation that in many eubacteria, mur genes are present in a single operon within the well conserved dcw cluster, and in some cases, pairs of mur genes are fused to encode a single, chimeric polypeptide. We performed a vast genomic analysis using >140 bacterial genomes and mapped Mur chimeras in numerous phyla, with Proteobacteria carrying the highest number. MurE-MurF, the most prevalent chimera, exists in forms that are either directly associated or separated by a linker. The crystal structure of the MurE-MurF chimera from Bordetella pertussis reveals a head-to-tail, elongated architecture supported by an interconnecting hydrophobic patch that stabilizes the positions of the two proteins. Fluorescence polarization assays reveal that MurE-MurF interacts with other Mur ligases via its central domains with KDs in the high nanomolar range, backing the existence of a Mur complex in the cytoplasm. These data support the idea of stronger evolutionary constraints on gene order when encoded proteins are intended for association, establish a link between Mur ligase interaction, complex assembly and genome evolution, and shed light on regulatory mechanisms of protein expression and stability in pathways of critical importance for bacterial survival.

Folding of heterologous proteins in bacterial cell factories: Cellular mechanisms and engineering strategies

The expression of correctly folded and functional heterologous proteins is important in many biotechnological production processes, whether it is enzymes, biopharmaceuticals or biosynthetic pathways for production of sustainable chemicals. For industrial applications, bacterial platform organisms, such as E. coli, are still broadly used due to the availability of tools and proven suitability at industrial scale. However, expression of heterologous proteins in these organisms can result in protein aggregation and low amounts of functional protein. This review provides an overview of the cellular mechanisms that can influence protein folding and expression, such as co-translational folding and assembly, chaperone binding, as well as protein quality control, across different model organisms. The knowledge of these mechanisms is then linked to different experimental methods that have been applied in order to improve functional heterologous protein folding, such as codon optimization, fusion tagging, chaperone co-production, as well as strain and protein engineering strategies.

Building Quantitative Bridges between Dynamics and Sequences of SARS-CoV-2 Main Protease and a Diverse Set of Thirty-Two Proteins

Proteases are major drug targets for many viral diseases. However, mutations can render several antiprotease drugs inefficient rapidly even though these mutations may not alter protein structures significantly. Understanding relations between quickly mutating residues, protease structures, and the dynamics of the proteases is crucial for designing potent drugs. Due to this reason, we studied relations between the evolutionary information on residues in the amino acid sequences and protein dynamics for SARS-CoV-2 main protease. More precisely, we analyzed three dynamical quantities (Schlitter entropy, root-mean-square fluctuations, and dynamical flexibility index) and their relation to the amino acid conservation extracted from multiple sequence alignments of the main protease. We showed that a quantifiable similarity can be built between a sequence-based quantity called Jensen-Shannon conservation and those three dynamical quantities. We validated this similarity for a diverse set of 32 different proteins, other than the SARS-CoV-2 main protease. We believe that establishing these kinds of quantitative bridges will have larger implications for all viral proteases as well as all proteins.

Simulation toolkits at the molecular scale for trans-scale thermal signaling

Thermogenesis is a physiological activity of releasing heat that originates from intracellular biochemical reactions. Recent experimental studies discovered that externally applied heat changes intracellular signaling locally, resulting in global changes in cell morphology and signaling. Therefore, we hypothesize an inevitable contribution of thermogenesis in modulating biological system functions throughout the spatial scales from molecules to individual organisms. One key issue examining the hypothesis, namely, the "trans-scale thermal signaling," resides at the molecular scale on the amount of heat released via individual reactions and by which mechanism the heat is employed for cellular function operations. This review introduces atomistic simulation tool kits for studying the mechanisms of thermal signaling processes at the molecular scale that even state-of-the-art experimental methodologies of today are hardly accessible. We consider biological processes and biomolecules as potential heat sources in cells, such as ATP/GTP hydrolysis and multiple biopolymer complex formation and disassembly. Microscopic heat release could be related to mesoscopic processes via thermal conductivity and thermal conductance. Additionally, theoretical simulations to estimate these thermal properties in biological membranes and proteins are introduced. Finally, we envisage the future direction of this research field.

Simple mechanisms for the evolution of protein complexity

Proteins are tiny models of biological complexity: specific interactions among their many amino acids cause proteins to fold into elaborate structures, assemble with other proteins into higher-order complexes, and change their functions and structures upon binding other molecules. These complex features are classically thought to evolve via long and gradual trajectories driven by persistent natural selection. But a growing body of evidence from biochemistry, protein engineering, and molecular evolution shows that naturally occurring proteins often exist at or near the genetic edge of multimerization, allostery, and even new folds, so just one or a few mutations can trigger acquisition of these properties. These sudden transitions can occur because many of the physical properties that underlie these features are present in simpler proteins as fortuitous by-products of their architecture. Moreover, complex features of proteins can be encoded by huge arrays of sequences, so they are accessible from many different starting points via many possible paths. Because the bridges to these features are both short and numerous, random chance can join selection as a key factor in explaining the evolution of molecular complexity.

Predicting the structure of large protein complexes using AlphaFold and Monte Carlo tree search

AlphaFold can predict the structure of single- and multiple-chain proteins with very high accuracy. However, the accuracy decreases with the number of chains, and the available GPU memory limits the size of protein complexes which can be predicted. Here we show that one can predict the structure of large complexes starting from predictions of subcomponents. We assemble 91 out of 175 complexes with 10-30 chains from predicted subcomponents using Monte Carlo tree search, with a median TM-score of 0.51. There are 30 highly accurate complexes (TM-score ≥0.8, 33% of complete assemblies). We create a scoring function, mpDockQ, that can distinguish if assemblies are complete and predict their accuracy. We find that complexes containing symmetry are accurately assembled, while asymmetrical complexes remain challenging. The method is freely available and accesible as a Colab notebook https://colab.research.google.com/github/patrickbryant1/MoLPC/blob/master/MoLPC.ipynb .

RL-MLZerD: Multimeric protein docking using reinforcement learning

Numerous biological processes in a cell are carried out by protein complexes. To understand the molecular mechanisms of such processes, it is crucial to know the quaternary structures of the complexes. Although the structures of protein complexes have been determined by biophysical experiments at a rapid pace, there are still many important complex structures that are yet to be determined. To supplement experimental structure determination of complexes, many computational protein docking methods have been developed; however, most of these docking methods are designed only for docking with two chains. Here, we introduce a novel method, RL-MLZerD, which builds multiple protein complexes using reinforcement learning (RL). In RL-MLZerD a multi-chain assembly process is considered as a series of episodes of selecting and integrating pre-computed pairwise docking models in a RL framework. RL is effective in correctly selecting plausible pairwise models that fit well with other subunits in a complex. When tested on a benchmark dataset of protein complexes with three to five chains, RL-MLZerD showed better modeling performance than other existing multiple docking methods under different evaluation criteria, except against AlphaFold-Multimer in unbound docking. Also, it emerged that the docking order of multi-chain complexes can be naturally predicted by examining preferred paths of episodes in the RL computation.

Cotranslational Mechanisms of Protein Biogenesis and Complex Assembly in Eukaryotes

The formation of protein complexes is crucial to most biological functions. The cellular mechanisms governing protein complex biogenesis are not yet well understood, but some principles of cotranslational and posttranslational assembly are beginning to emerge. In bacteria, this process is favored by operons encoding subunits of protein complexes. Eukaryotic cells do not have polycistronic mRNAs, raising the question of how they orchestrate the encounter of unassembled subunits. Here we review the constraints and mechanisms governing eukaryotic co- and posttranslational protein folding and assembly, including the influence of elongation rate on nascent chain targeting, folding, and chaperone interactions. Recent evidence shows that mRNAs encoding subunits of oligomeric assemblies can undergo localized translation and form cytoplasmic condensates that might facilitate the assembly of protein complexes. Understanding the interplay between localized mRNA translation and cotranslational proteostasis will be critical to defining protein complex assembly in vivo.

Large protein complex interfaces have evolved to promote cotranslational assembly

Assembly pathways of protein complexes should be precise and efficient to minimise misfolding and unwanted interactions with other proteins in the cell. One way to achieve this efficiency is by seeding assembly pathways during translation via the cotranslational assembly of subunits. While recent evidence suggests that such cotranslational assembly is widespread, little is known about the properties of protein complexes associated with the phenomenon. Here, using a combination of proteome-specific protein complex structures and publicly available ribosome profiling data, we show that cotranslational assembly is particularly common between subunits that form large intermolecular interfaces. To test whether large interfaces have evolved to promote cotranslational assembly, as opposed to cotranslational assembly being a non-adaptive consequence of large interfaces, we compared the sizes of first and last translated interfaces of heteromeric subunits in bacterial, yeast, and human complexes. When considering all together, we observe the N-terminal interface to be larger than the C-terminal interface 54% of the time, increasing to 64% when we exclude subunits with only small interfaces, which are unlikely to cotranslationally assemble. This strongly suggests that large interfaces have evolved as a means to maximise the chance of successful cotranslational subunit binding.

Evolution of homo-oligomerization of methionine S-adenosyltransferases is replete with structure-function constrains

Homomers are prevalent in bacterial proteomes, particularly among core metabolic enzymes. Homomerization is often key to function and regulation, and interfaces that facilitate the formation of homomeric enzymes are subject to intense evolutionary change. However, our understanding of the molecular mechanisms that drive evolutionary variation in homomeric complexes is still lacking. How is the diversification of protein interfaces linked to variation in functional regulation and structural integrity of homomeric complexes? To address this question, we studied quaternary structure evolution of bacterial methionine S-adenosyltransferases (MATs)-dihedral homotetramers formed along a large and conserved dimeric interface harboring two active sites, and a small, recently evolved, interdimeric interface. Here, we show that diversity in the physicochemical properties of small interfaces is directly linked to variability in the kinetic stability of MAT quaternary complexes and in modes of their functional regulation. Specifically, hydrophobic interactions within the small interface of Escherichia coli MAT render the functional homotetramer kinetically stable yet impose severe aggregation constraints on complex assembly. These constraints are alleviated by electrostatic interactions that accelerate dimer-dimer assembly. In contrast, Neisseria gonorrhoeae MAT adopts a nonfunctional dimeric state due to the low hydrophobicity of its small interface and the high flexibility of its active site loops, which perturbs small interface integrity. Remarkably, in the presence of methionine and ATP, N. gonorrhoeae MAT undergoes substrate-induced assembly into a functional tetrameric state. We suggest that evolution acts on the interdimeric interfaces of MATs to tailor the regulation of their activity and stability to unique organismal needs.

STAT3 in medulloblastoma: a key transcriptional regulator and potential therapeutic target

Medulloblastoma is the most common malignant brain tumor of childhood accounting for about 60% of all pediatric embryonal tumors. Despite improvements in the overall survival rate, this tumor still lacks an efficient, reliable, and less toxic therapeutic approach. Characterization of the molecular mechanisms involved in medulloblastoma initiation and progression is a crucial step for the development of effective therapies. Signal transducer and activator of transcription 3 is a convergence point for several signaling cascades that are implicated in medulloblastoma tumorigenesis. Accumulated evidence has revealed the pivotal role of signal transducer and activator of transcription 3 in medulloblastoma pathogenesis such as proliferation, survival, angiogenesis, and immunosuppression as well as maintenance, drug resistance, and recurrence. In this review, we focus on the role of signal transducer and activator of transcription 3 in medulloblastoma tumorigenesis and discuss the recent advances of signal transducer and activator of transcription 3 inhibition as a promising developed strategy for medulloblastoma therapy.

Biogenesis of a bacterial metabolosome for propanediol utilization

Bacterial metabolosomes are a family of protein organelles in bacteria. Elucidating how thousands of proteins self-assemble to form functional metabolosomes is essential for understanding their significance in cellular metabolism and pathogenesis. Here we investigate the de novo biogenesis of propanediol-utilization (Pdu) metabolosomes and characterize the roles of the key constituents in generation and intracellular positioning of functional metabolosomes. Our results demonstrate that the Pdu metabolosome undertakes both "Shell first" and "Cargo first" assembly pathways, unlike the β-carboxysome structural analog which only involves the "Cargo first" strategy. Shell and cargo assemblies occur independently at the cell poles. The internal cargo core is formed through the ordered assembly of multiple enzyme complexes, and exhibits liquid-like properties within the metabolosome architecture. Our findings provide mechanistic insight into the molecular principles driving bacterial metabolosome assembly and expand our understanding of liquid-like organelle biogenesis.

Rethinking the MtInhA tertiary and quaternary structure flexibility: a molecular dynamics view

Flexibility and function are related properties in the study of protein dynamics. Flexibility reflects in the conformational potential of proteins and thus in their functionalities. The presence of interactions between protein-ligands and protein-protein complexes, substrates, and environmental changes can alter protein plasticity, acting from the rearrangement of the side chains of amino acids to the folding/unfolding of large structural motifs. To evaluate the effects of the flexibility in protein systems, we defined the enzyme 2-trans-enoyl-ACP (CoA) reductase from Mycobacterium tuberculosis, or MtInhA, as our target system. MtInhA is biologically active as a tetramer in solution; however, computational studies commonly use the monomer justifying the independence of its active sites due to their distances. However, differences in flexibility between tertiary and quaternary structures could present impact on the size of the active site, influencing the drug discovery process. In this study, we investigated the influence of flexibility restrictions in A- and B-loops of the MtInhA in order to suggest a monomeric structure that describes the conformational behavior of the tetrameric system. Overall, we observed that simulations where restrictions were applied to the A- and B-loops present a more similar behavior to the native structure when compared to unrestricted simulations. Therefore, our work presents a monomeric model of MtInhA, which has conformational characteristics of the biologically active structure. Thus, the data obtained in this work can be applied to the MtInhA system for the generation of more reliable flexible models for molecular docking experiments, and also for the performance of longer simulations by molecular dynamics and with a lower computational cost.

Computational prediction of heteromeric protein complex disassembly order using hybrid Monte Carlo/molecular dynamics simulation

The physicochemical entities comprising the biological phenomena in the cell form a network of biochemical reactions and the activity of such a network is regulated by multimeric protein complexes. Mass spectroscopy (MS) experiments and multimeric protein docking simulations based on structural bioinformatics techniques have revealed the molecular-level stoichiometry and static configuration of subcomplexes in their bound forms, thus revealing the subcomplex population and formation orders. Meanwhile, these methodologies are not designed to straightforwardly examine the temporal dynamics of multimeric protein assembly and disassembly, essential physicochemical properties to understand the functional expression mechanisms of proteins in the biological environment. To address this problem, we have developed an atomistic simulation in the framework of the hybrid Monte Carlo/molecular dynamics (hMC/MD) method and succeeded in observing the disassembly of a homomeric pentamer of the serum amyloid P component protein in an experimentally consistent order. In this study, we improved the hMC/MD method to examine the disassembly processes of the tryptophan synthase tetramer, a paradigmatic heteromeric protein complex in MS studies. We employed the likelihood-based selection scheme to determine a dissociation-prone subunit pair at every hMC/MD simulation cycle and achieved highly reliable predictions of the disassembly orders without a priori knowledge of the MS experiments and structural bioinformatics simulations. The success rate for the experimentally-observed disassembly order is over 0.9. We similarly succeeded in reliable predictions for three other tetrameric protein complexes. These achievements indicate the potential applicability of our hMC/MD approach as a general-purpose methodology to obtain microscopic and physicochemical insights into multimeric protein complex formation.

Surface-induced Dissociation Mass Spectrometry as a Structural Biology Tool

Native mass spectrometry (nMS) is evolving into a workhorse for structural biology. The plethora of online and offline preparation, separation, and purification methods as well as numerous ionization techniques combined with powerful new hybrid ion mobility and mass spectrometry systems has illustrated the great potential of nMS for structural biology. Fundamental to the progression of nMS has been the development of novel activation methods for dissociating proteins and protein complexes to deduce primary, secondary, tertiary, and quaternary structure through the combined use of multiple MS/MS technologies. This review highlights the key features and advantages of surface collisions (surface-induced dissociation, SID) for probing the connectivity of subunits within protein and nucleoprotein complexes and, in particular, for solving protein structure in conjunction with complementary techniques such as cryo-EM and computational modeling. Several case studies highlight the significant role SID, and more generally nMS, will play in structural elucidation of biological assemblies in the future as the technology becomes more widely adopted. Cases are presented where SID agrees with solved crystal or cryoEM structures or provides connectivity maps that are otherwise inaccessible by "gold standard" structural biology techniques.

Co-Translational Folding of Multi-Domain Proteins

The majority of proteins in nature are composed of multiple domains connected in a single polypeptide. How these long sequences fold into functional structures without forming toxic misfolds or aggregates is poorly understood. Their folding is inextricably linked to protein synthesis and interactions with cellular machinery, making mechanistic studies challenging. Recent progress has revealed critical features of multi-domain protein folding in isolation and in the context of translation by the ribosome. In this review, we discuss challenges and progress in understanding multi-domain protein folding, and highlight how molecular interactions shape folding and misfolding pathways. With the development of new approaches and model systems, the stage is now set for mechanistically exploring the folding of large multi-domain proteins.

Quality control of protein complex assembly by the ubiquitin-proteasome system

The majority of human proteins operate as multimeric complexes with defined compositions and distinct architectures. How the assembly of these complexes is surveyed and how defective complexes are recognized is just beginning to emerge. In eukaryotes, over 600 E3 ubiquitin ligases form part of the ubiquitin-proteasome system (UPS) which detects structural characteristics in its target proteins and selectively induces their degradation. The UPS has recently been shown to oversee key quality control steps during the assembly of protein complexes. We review recent findings on how E3 ubiquitin ligases regulate protein complex assembly and highlight unanswered questions relating to their mechanism of action.

Co-translational assembly orchestrates competing biogenesis pathways

During the co-translational assembly of protein complexes, a fully synthesized subunit engages with the nascent chain of a newly synthesized interaction partner. Such events are thought to contribute to productive assembly, but their exact physiological relevance remains underexplored. Here, we examine structural motifs contained in nucleoporins for their potential to facilitate co-translational assembly. We experimentally test candidate structural motifs and identify several previously unknown co-translational interactions. We demonstrate by selective ribosome profiling that domain invasion motifs of beta-propellers, coiled-coils, and short linear motifs may act as co-translational assembly domains. Such motifs are often contained in proteins that are members of multiple complexes (moonlighters) and engage with closely related paralogs. Surprisingly, moonlighters and paralogs assemble co-translationally in only some but not all of the relevant biogenesis pathways. Our results highlight the regulatory complexity of assembly pathways.

The properties of human disease mutations at protein interfaces

The assembly of proteins into complexes and their interactions with other biomolecules are often vital for their biological function. While it is known that mutations at protein interfaces have a high potential to be damaging and cause human genetic disease, there has been relatively little consideration for how this varies between different types of interfaces. Here we investigate the properties of human pathogenic and putatively benign missense variants at homomeric (isologous and heterologous), heteromeric, DNA, RNA and other ligand interfaces, and at different regions in proteins with respect to those interfaces. We find that different types of interfaces vary greatly in their propensity to be associated with pathogenic mutations, with homomeric heterologous and DNA interfaces being particularly enriched in disease. We also find that residues that do not directly participate in an interface, but are close in three-dimensional space, show a significant disease enrichment. Finally, we observe that mutations at different types of interfaces tend to have distinct property changes when undergoing amino acid substitutions associated with disease, and that this is linked to substantial variability in their identification by computational variant effect predictors.

Nearly all proteins interact with other molecules as part of their biological function. For example, proteins can interact with other copies of the same type of protein, with different proteins, with DNA, or with small ligand molecules. Many mutations at protein interfaces, the regions of proteins that interact with other molecules, are known to cause human genetic disease. In this study, we first investigate how different types of protein interfaces have different tendencies to be associated with disease. We also show that the closer a mutation is to an interface, the more likely it is to cause disease. Finally, we study how mutations at different types of interfaces tend to be associated with different changes in amino acid properties, which appears to influence our ability to computationally predict the effects of mutations. Ultimately, we hope that consideration of protein interface properties will eventually improve our ability to identify new disease-causing mutations.

Sequence Divergence and Functional Specializations of the Ancient Spliceosomal SF3b: Implications in Flexibility and Adaptations of the Multi-Protein Complex

Multi-protein assemblies are complex molecular systems that perform highly sophisticated biochemical functions in an orchestrated manner. They are subject to changes that are governed by the evolution of individual components. We performed a comparative analysis of the ancient and functionally conserved spliceosomal SF3b complex, to recognize molecular signatures that contribute to sequence divergence and functional specializations. For this, we recognized homologous sequences of individual SF3b proteins distributed across 10 supergroups of eukaryotes and identified all seven protein components of the complex in 578 eukaryotic species. Using sequence and structural analysis, we establish that proteins occurring on the surface of the SF3b complex harbor more sequence variation than the proteins that lie in the core. Further, we show through protein interface conservation patterns that the extent of conservation varies considerably between interacting partners. When we analyze phylogenetic distributions of individual components of the complex, we find that protein partners that are known to form independent subcomplexes are observed to share similar profiles, reaffirming the link between differential conservation of interface regions and their inter-dependence. When we extend our analysis to individual protein components of the complex, we find taxa-specific variability in molecular signatures of the proteins. These trends are discussed in the context of proline-rich motifs of SF3b4, functional and drug binding sites of SF3b1. Further, we report key protein-protein interactions between SF3b1 and SF3b6 whose presence is observed to be lineage-specific across eukaryotes. Together, our studies show the association of protein location within the complex and subcomplex formation patterns with the sequence conservation of SF3b proteins. In addition, our study underscores evolutionarily flexible elements that appear to confer adaptive features in individual components of the multi-protein SF3b complexes and may contribute to its functional adaptability.

Gene utility recapitulates chromosomal aberrancies in advanced stage neuroblastoma

Neuroblastoma (NB) is the most common extracranial solid tumor in children. Although only a few recurrent somatic mutations have been identified, chromosomal abnormalities, including the loss of heterozygosity (LOH) at the chromosome 1p and gains of chromosome 17q, are often seen in the high-risk cases. The biological basis and evolutionary forces that drive such genetic abnormalities remain enigmatic. Here, we conceptualize the Gene Utility Model (GUM) that seeks to identify genes driving biological signaling via their collective gene utilities and apply it to understand the impact of those differentially utilized genes on constraining the evolution of NB karyotypes. By employing a computational process-guided flow algorithm to model gene utility in protein–protein networks that built based on transcriptomic data, we conducted several pairwise comparative analyses to uncover genes with differential utilities in stage 4 NBs with distinct classification. We then constructed a utility karyotype by mapping these differentially utilized genes to their respective chromosomal loci. Intriguingly, hotspots of the utility karyotype, to certain extent, can consistently recapitulate the major chromosomal abnormalities of NBs and also provides clues to yet identified predisposition sites. Hence, our study not only provides a new look, from a gene utility perspective, into the known chromosomal abnormalities detected by integrative genomic sequencing efforts, but also offers new insights into the etiology of NB and provides a framework to facilitate the identification of novel therapeutic targets for this devastating childhood cancer.

Assembly status transition offers an avenue for activity modulation of a supramolecular enzyme

Nature has evolved many supramolecular proteins assembled in certain, sometimes even seemingly oversophisticated, morphological manners. The rationale behind such evolutionary efforts is often poorly understood. Here, we provide atomic-resolution insights into how the dynamic building of a structurally complex enzyme with higher order symmetry offers amenability to intricate regulation. We have established the functional coupling between enzymatic activity and protein morphological states of glutamine synthetase (GS), an old multi-subunit enzyme essential for cellular nitrogen metabolism. Cryo-EM structure determination of GS in both the catalytically active and inactive assembly states allows us to reveal an unanticipated self-assembly-induced disorder-order transition paradigm, in which the remote interactions between two subcomplex entities significantly rigidify the otherwise structurally fluctuating active sites, thereby regulating activity. We further show in vivo evidences that how the enzyme morphology transitions could be modulated by cellular factors on demand. Collectively, our data present an example of how assembly status transition offers an avenue for activity modulation, and sharpens our mechanistic understanding of the complex functional and regulatory properties of supramolecular enzymes.

Structural Analysis of Saccharomyces cerevisiae Dihydroorotase Reveals Molecular Insights into the Tetramerization Mechanism

Dihydroorotase (DHOase), a dimetalloenzyme containing a carbamylated lysine within the active site, is a member of the cyclic amidohydrolase family, which also includes allantoinase (ALLase), dihydropyrimidinase (DHPase), hydantoinase, and imidase. Unlike most known cyclic amidohydrolases, which are tetrameric, DHOase exists as a monomer or dimer. Here, we report and analyze two crystal structures of the eukaryotic Saccharomyces cerevisiae DHOase (ScDHOase) complexed with malate. The structures of different DHOases were also compared. An asymmetric unit of these crystals contained four crystallographically independent ScDHOase monomers. ScDHOase shares structural similarity with Escherichia coli DHOase (EcDHOase). Unlike EcDHOase, ScDHOase can form tetramers, both in the crystalline state and in solution. In addition, the subunit-interacting residues of ScDHOase for dimerization and tetramerization are significantly different from those of other DHOases. The tetramerization pattern of ScDHOase is also different from those of DHPase and ALLase. Based on sequence analysis and structural evidence, we identify two unique helices (α6 and α10) and a loop (loop 7) for tetramerization, and discuss why the residues for tetramerization in ScDHOase are not necessarily conserved among DHOases.

The Hsp90 cochaperone TTT promotes cotranslational maturation of PIKKs prior to complex assembly

Phosphatidylinositol 3-kinase-related kinases (PIKKs) are a family of kinases that control fundamental processes, including cell growth, DNA damage repair, and gene expression. Although their regulation and activities are well characterized, little is known about how PIKKs fold and assemble into active complexes. Previous work has identified a heat shock protein 90 (Hsp90) cochaperone, the TTT complex, that specifically stabilizes PIKKs. Here, we describe a mechanism by which TTT promotes their de novo maturation in fission yeast. We show that TTT recognizes newly synthesized PIKKs during translation. Although PIKKs form multimeric complexes, we find that they do not engage in cotranslational assembly with their partners. Rather, our findings suggest a model by which TTT protects nascent PIKK polypeptides from misfolding and degradation because PIKKs acquire their native state after translation is terminated. Thus, PIKK maturation and assembly are temporally segregated, suggesting that the biogenesis of large complexes requires both dedicated chaperones and cotranslational interactions between subunits.

Native mass spectrometry analyses of chaperonin complex TRiC/CCT reveal subunit N-terminal processing and re-association patterns

The eukaryotic chaperonin TRiC/CCT is a large ATP-dependent complex essential for cellular protein folding. Its subunit arrangement into two stacked eight-membered hetero-oligomeric rings is conserved from yeast to man. A recent breakthrough enables production of functional human TRiC (hTRiC) from insect cells. Here, we apply a suite of mass spectrometry techniques to characterize recombinant hTRiC. We find all subunits CCT1-8 are N-terminally processed by combinations of methionine excision and acetylation observed in native human TRiC. Dissociation by organic solvents yields primarily monomeric subunits with a small population of CCT dimers. Notably, some dimers feature non-canonical inter-subunit contacts absent in the initial hTRiC. This indicates individual CCT monomers can promiscuously re-assemble into dimers, and lack the information to assume the specific interface pairings in the holocomplex. CCT5 is consistently the most stable subunit and engages in the greatest number of non-canonical dimer pairings. These findings confirm physiologically relevant post-translational processing and function of recombinant hTRiC and offer quantitative insight into the relative stabilities of TRiC subunits and interfaces, a key step toward reconstructing its assembly mechanism. Our results also highlight the importance of assigning contacts identified by native mass spectrometry after solution dissociation as canonical or non-canonical when investigating multimeric assemblies.

Unraveling Genome Evolution Throughout Visual Analysis: The XCout Portal

Due to major breakthroughs in sequencing technologies throughout the last decades, the time and cost per sequencing experiment have reduced drastically, overcoming the data generation barrier during the early genomic era. Such a shift has encouraged the scientific community to develop new computational methods that are able to compare large genomic sequences, thus enabling large-scale studies of genome evolution. The field of comparative genomics has proven itself invaluable for studying the evolutionary mechanisms and the forces driving genome evolution. In this line, a full genome comparison study between 2 species requires a quadratic number of comparisons in terms of the number of sequences (around 400 chromosome comparisons in the case of mammalian genomes); however, when studying conserved syntenies or evolutionary rearrangements, many sequence comparisons can be skipped for not all will contain significant signals. Subsequently, the scientific community has developed fast heuristics to perform multiple pairwise comparisons between large sequences to determine whether significant sets of conserved similarities exist. The data generation problem is no longer an issue, yet the limitations have shifted toward the analysis of such massive data. Therefore, we present XCout, a Web-based visual analytics application for multiple genome comparisons designed to improve the analysis of large-scale evolutionary studies using novel techniques in Web visualization. XCout enables to work on hundreds of comparisons at once, thus reducing the time of the analysis by identifying significant signals between chromosomes across multiple species. Among others, XCout introduces several techniques to aid in the analysis of large-scale genome rearrangements, particularly (1) an interactive heatmap interface to display comparisons using automatic color scales based on similarity thresholds to ease detection at first sight, (2) an overlay system to detect individual signal contributions between chromosomes, (3) a tracking tool to trace conserved blocks across different species to perform evolutionary studies, and (4) a search engine to search annotations throughout different species.

Mechanisms of Cotranslational Protein Maturation in Bacteria

Growing cells invest a significant part of their biosynthetic capacity into the production of proteins. To become functional, newly-synthesized proteins must be N-terminally processed, folded and often translocated to other cellular compartments. A general strategy is to integrate these protein maturation processes with translation, by cotranslationally engaging processing enzymes, chaperones and targeting factors with the nascent polypeptide. Precise coordination of all factors involved is critical for the efficiency and accuracy of protein synthesis and cellular homeostasis. This review provides an overview of the current knowledge on cotranslational protein maturation, with a focus on the production of cytosolic proteins in bacteria. We describe the role of the ribosome and the chaperone network in protein folding and how the dynamic interplay of all cotranslationally acting factors guides the sequence of cotranslational events. Finally, we discuss recent data demonstrating the coupling of protein synthesis with the assembly of protein complexes and end with a brief discussion of outstanding questions and emerging concepts in the field of cotranslational protein maturation.

Protein Complexes Form a Basis for Complex Hybrid Incompatibility

Proteins are the workhorses of the cell and execute many of their functions by interacting with other proteins forming protein complexes. Multi-protein complexes are an admixture of subunits, change their interaction partners, and modulate their functions and cellular physiology in response to environmental changes. When two species mate, the hybrid offspring are usually inviable or sterile because of large-scale differences in the genetic makeup between the two parents causing incompatible genetic interactions. Such reciprocal-sign epistasis between inter-specific alleles is not limited to incompatible interactions between just one gene pair; and, usually involves multiple genes. Many of these multi-locus incompatibilities show visible defects, only in the presence of all the interactions, making it hard to characterize. Understanding the dynamics of protein-protein interactions (PPIs) leading to multi-protein complexes is better suited to characterize multi-locus incompatibilities, compared to studying them with traditional approaches of genetics and molecular biology. The advances in omics technologies, which includes genomics, transcriptomics, and proteomics can help achieve this end. This is especially relevant when studying non-model organisms. Here, we discuss the recent progress in the understanding of hybrid genetic incompatibility; omics technologies, and how together they have helped in characterizing protein complexes and in turn multi-locus incompatibilities. We also review advances in bioinformatic techniques suitable for this purpose and propose directions for leveraging the knowledge gained from model-organisms to identify genetic incompatibilities in non-model organisms.

A Supramolecular Assembly of Hemoproteins Formed in a Star-Shaped Structure via Heme-Heme Pocket Interactions

Proteins have been used as building blocks to provide various supramolecular structures in efforts to develop nano-biomaterials possessing broad biological functionalities. A series of unique structures have been obtained from the engineering of hemoproteins which contain the iron porphyrin known as heme, as a prosthetic group. This work in developing assembling systems is extended using cytochrome b₅₆₂, a small electron transfer hemoprotein engineered to include an externally-attached heme moiety. The engineered units, which form a one-dimensional assembly via interprotein heme–heme pocket interactions, are conjugated to an apo-form of hexameric tyrosine-coordinated hemoprotein (apoHTHP) to provide a branching unit promoting the assembly of a star-shaped structure. The incorporation of the heme moiety attached to the protein surface of cytochrome b₅₆₂ into apoHTHP can be accelerated by elevating the reaction temperature to generate a new assembly. The formation of a new larger assembly structure was confirmed by size exclusion chromatography. The ratio of the heme-containing units in the assemblies was analyzed by UV-Vis spectroscopy and the population of protein units estimated from SDS PAGE suggests the presence of plausible star-shaped structures, which are supported by hydrodynamic diameter data obtained by dynamic light scattering.