The folding and evolution of multidomain proteins

Mechanism of chaperone coordination during cotranslational protein folding in bacteria

Protein folding is assisted by molecular chaperones that bind nascent polypeptides during mRNA translation. Several structurally distinct classes of chaperones promote de novo folding, suggesting that their activities are coordinated at the ribosome. We used biochemical reconstitution and structural proteomics to explore the molecular basis for cotranslational chaperone action in bacteria. We found that chaperone binding is disfavored close to the ribosome, allowing folding to precede chaperone recruitment. Trigger factor recognizes compact folding intermediates that expose an extensive unfolded surface, and dictates DnaJ access to nascent chains. DnaJ uses a large surface to bind structurally diverse intermediates and recruits DnaK to sequence-diverse solvent-accessible sites. Neither Trigger factor, DnaJ, nor DnaK destabilize cotranslational folding intermediates. Instead, the chaperones collaborate to protect incipient structure in the nascent polypeptide well beyond the ribosome exit tunnel. Our findings show how the chaperone network selects and modulates cotranslational folding intermediates.

Folding and Binding Kinetics of the Tandem of SH2 Domains from SHP2

The SH2 domains of SHP2 play a crucial role in determining the function of the SHP2 protein. While the folding and binding properties of the isolated NSH2 and CSH2 domains have been extensively studied, there is limited information about the tandem SH2 domains. This study aims to elucidate the folding and binding kinetics of the NSH2-CSH2 tandem domains of SHP2 through rapid kinetic experiments, complementing existing data on the isolated domains. The results indicate that while the domains generally fold and unfold independently, acidic pH conditions induce complex scenarios involving the formation of a misfolded intermediate. Furthermore, a comparison of the binding kinetics of isolated NSH2 and CSH2 domains with the NSH2-CSH2 tandem domains, using peptides that mimic specific portions of Gab2, suggests a dynamic interplay between NSH2 and CSH2 in binding Gab2 that modulate the microscopic association rate constant of the binding reaction. These findings, discussed in the context of previous research on the NSH2 and CSH2 domains, enhance our understanding of the function of the SH2 domain tandem of SHP2.

Navigating the complexities of multi-domain protein folding

Proteome complexity has expanded tremendously over evolutionary time, enabling biological diversification. Much of this complexity is achieved by combining a limited set of structural units into long polypeptides. This widely used evolutionary strategy poses challenges for folding of the resulting multi-domain proteins. As a consequence, their folding differs from that of small single-domain proteins, which generally fold quickly and reversibly. Co-translational processes and chaperone interactions are important aspects of multi-domain protein folding. In this review, we discuss some of the recent experimental progress toward understanding these processes.

Evaluation of Essential Dynamics and Fixed-Length Coarse Graining for Multidomain Proteins

For multiscale modeling of biomolecules, reliable coarse-grained (CG) models can offer great potential to simulate larger temporal and spatial scales than traditional all-atom (AA) models. In this study, we explore the essential dynamics coarse graining (EDCG) and fixed-length coarse graining (FLCG) approaches for constructing highly coarse-grained models for multidomain proteins (MDPs), with 1 to 10 amino acid residues per CG site. In the studies of 13 MDPs, our data indicate that both EDCG and FLCG can preserve the protein dynamics of MDPs. FLCG, which restricts an equal number of residues in each CG site, represents an excellent approximation to EDCG and a straightforward approach for coarse-graining MDPs. Furthermore, FLCG is tested with a class B G-protein-coupled receptor protein, and the agreement with prior experiments suggests its general application to various MDPs in different environments or conditions. Finally, we demonstrate another application of FLCG through progressive backmapping, showcasing the ability to recover from lower-resolution CG models (6 residues/CG site) to higher-resolution ones (1 residue/CG site). These promising outcomes underscore the broad applicability of FLCG to construct highly or ultra-coarse-grained models of complex biomolecules for multiscale simulations.

Genome-wide analysis of cellulose synthase gene superfamily in Tectona grandis L.f

This study aimed to explore Cellulose synthase gene superfamily of teak, and its evolutionary relationship with homologous genes of other woody species. The incidence of evolutionary events like gene duplication and gene loss, influence of the selection pressure, and consequent adaptive functional divergence of the duplicated TgCes gene were assessed alongside it's role in wood coloration. This study identified 39 full-length non-redundant proteins belonging to CesA and Csl gene families. TgCesA and TgCsl proteins with Cellulose synthase domain repeats indicated tandem gene duplication and probable genetic variability, enabling local adaptation. Further, multi-domain protein (MYB-like DNA-binding domain and CesA domain) with maximum introns was also identified indicating gene fusion and formation of complex protein with novel functions. Phylogenetic analysis grouped the genes into seven subfamilies (CesA, CslA, CslC, CslD, CslE, CslG, and CslM) with each undergoing gene duplication and loss along their evolutionary history. Post-species gene duplications and probable neofunctionalization were identified in TgCesA and TgCsl gene families. Each subfamily was found to be under strong purifying selection with a few or no sites under positive selection. Functional divergence analysis further revealed site-specific selective constraints in CesA and Csl genes of the teak Cellulose synthase gene family. Furthermore, protein-protein interaction network analysis identified co-expression of Cellulose synthase gene with flavonoid 3',5'-hydroxylase (F3'5'H, CYP75A), involved in the biosynthesis of xylem anthocyanin compounds, probably responsible for wood coloration. This study thus offers a foundation for future research in wood formation and wood property traits specific to teak and its provenances.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-03927-6.

Conflicting Interfacial Electrostatic Interactions as a Design Principle to Modulate Long-Range Interdomain Communication

The extent and molecular basis of interdomain communication in multidomain proteins, central to understanding allostery and function, is an open question. One simple evolutionary strategy could involve the selection of either conflicting or favorable electrostatic interactions across the interface of two closely spaced domains to tune the magnitude of interdomain connectivity. Here, we study a bilobed domain FF34 from the eukaryotic p190A RhoGAP protein to explore one such design principle that mediates interdomain communication. We find that while the individual structural units in wild-type FF34 are marginally coupled, they exhibit distinct intrinsic stabilities and low cooperativity, manifesting as slow folding. The FF3-FF4 interface harbors a frustrated network of highly conserved electrostatic interactions-a charge troika-that promotes the population of multiple, decoupled, and non-native structural modes on a rugged native landscape. Perturbing this network via a charge-reversal mutation not only enhances stability and cooperativity but also dampens the fluctuations globally and speeds up the folding rate by at least an order of magnitude. Our work highlights how a conserved but nonoptimal network of interfacial electrostatic interactions shapes the native ensemble of a bilobed protein, a feature that could be exploited in designing molecular systems with long-range connectivity and enhanced cooperativity.

N-glycosylation as a eukaryotic protective mechanism against protein aggregation

The tendency for proteins to form aggregates is an inherent part of every proteome and arises from the self-assembly of short protein segments called aggregation-prone regions (APRs). While posttranslational modifications (PTMs) have been implicated in modulating protein aggregation, their direct role in APRs remains poorly understood. In this study, we used a combination of proteome-wide computational analyses and biophysical techniques to investigate the potential involvement of PTMs in aggregation regulation. Our findings reveal that while most PTM types are disfavored near APRs, N-glycosylation is enriched and evolutionarily selected, especially in proteins prone to misfolding. Experimentally, we show that N-glycosylation inhibits the aggregation of peptides in vitro through steric hindrance. Moreover, mining existing proteomics data, we find that the loss of N-glycans at the flanks of APRs leads to specific protein aggregation in Neuro2a cells. Our findings indicate that, among its many molecular functions, N-glycosylation directly prevents protein aggregation in higher eukaryotes.

A complex network of interdomain interactions underlies the conformational ensemble of monomeric TDP-43 and modulates its phase behavior

TAR DNA-binding protein 43 (TDP-43) is a multidomain protein involved in the regulation of RNA metabolism, and its aggregates have been observed in neurodegenerative diseases, including amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). Numerous studies indicate TDP-43 can undergo liquid-liquid phase separation (LLPS) in vitro and is a component of biological condensates. Homo-oligomerization via the folded N-terminal domain (aa:1-77) and the conserved helical region (aa:319-341) of the disordered, C-terminal domain is found to be an important driver of TDP-43 phase separation. However, a comprehensive molecular view of TDP-43 phase separation, particularly regarding the nature of heterodomain interactions, is lacking due to the challenges associated with its stability and purification. Here, we utilize all-atom and coarse-grained (CG) molecular dynamics (MD) simulations to uncover the network of interdomain interactions implicated in TDP-43 phase separation. All-atom simulations uncovered the presence of transient, interdomain interactions involving flexible linkers, RNA-recognition motif (RRM) domains and a charged segment of disordered C-terminal domain (CTD). CG simulations indicate these inter-domain interactions which affect the conformational landscape of TDP-43 in the dilute phase are also prevalent in the condensed phase. Finally, sequence and surface charge distribution analysis coupled with all-atom simulations (at high salt) confirmed that the transient interdomain contacts are predominantly electrostatic in nature. Overall, our findings from multiscale simulations lead to a greater appreciation of the complex interaction network underlying the structural landscape and phase separation of TDP-43.

Tandem repeats of highly bioluminescent NanoLuc are refolded noncanonically by the Hsp70 machinery

Chaperones are a large family of proteins crucial for maintaining cellular protein homeostasis. One such chaperone is the 70 kDa heat shock protein (Hsp70), which plays a crucial role in protein (re)folding, stability, functionality, and translocation. While the key events in the Hsp70 chaperone cycle are well established, a relatively small number of distinct substrates were repetitively investigated. This is despite Hsp70 engaging with a plethora of cellular proteins of various structural properties and folding pathways. Here we analyzed novel Hsp70 substrates, based on tandem repeats of NanoLuc (Nluc), a small and highly bioluminescent protein with unique structural characteristics. In previous mechanical unfolding and refolding studies, we have identified interesting misfolding propensities of these Nluc-based tandem repeats. In this study, we further investigate these properties through in vitro bulk experiments. Similar to monomeric Nluc, engineered Nluc dyads and triads proved to be highly bioluminescent. Using the bioluminescence signal as the proxy for their structural integrity, we determined that heat-denatured Nluc dyads and triads can be efficiently refolded by the E. coli Hsp70 chaperone system, which comprises DnaK, DnaJ, and GrpE. In contrast to previous studies with other substrates, we observed that Nluc repeats can be efficiently refolded by DnaK and DnaJ, even in the absence of GrpE co-chaperone. Taken together, our study offers a new powerful substrate for chaperone research and raises intriguing questions about the Hsp70 mechanisms, particularly in the context of structurally diverse proteins.

Conformational Distribution of a Multidomain Protein Measured by Single-Pair Small-Angle X-ray Scattering

The difficulty in evaluating the conformational distribution of proteins in solution often hinders mechanistic insights. One possible strategy for visualizing conformational distribution is distance distribution measurement by single-pair small-angle X-ray scattering (SAXS), in which the scattering interference from only a specific pair of atoms in the target molecule is extracted. Despite this promising concept, with few applications in synthetic small molecules and DNA, technical difficulties have prevented its application in protein conformational studies. This study used a synthetic tag to fix the lanthanide ion at desired sites on the protein and used single-pair SAXS with contrast matching to evaluate the conformational distribution of the multidomain protein enzyme MurD. These data highlighted the broad conformational and ligand-driven distribution shifts of MurD in solution. This study proposes an important strategy in solution structural biology that targets dynamic proteins, including multidomain and intrinsically disordered proteins.

Folding pathway of a discontinuous two-domain protein

It is estimated that two-thirds of all proteins in higher organisms are composed of multiple domains, many of them containing discontinuous folds. However, to date, most in vitro protein folding studies have focused on small, single-domain proteins. As a model system for a two-domain discontinuous protein, we study the unfolding/refolding of a slow-folding double mutant of the maltose binding protein (DM-MBP) using single-molecule two- and three-color Förster Resonance Energy Transfer experiments. We observe a dynamic folding intermediate population in the N-terminal domain (NTD), C-terminal domain (CTD), and at the domain interface. The dynamic intermediate fluctuates rapidly between unfolded states and compact states, which have a similar FRET efficiency to the folded conformation. Our data reveals that the delayed folding of the NTD in DM-MBP is imposed by an entropic barrier with subsequent folding of the highly dynamic CTD. Notably, accelerated DM-MBP folding is routed through the same dynamic intermediate within the cavity of the GroEL/ES chaperone system, suggesting that the chaperonin limits the conformational space to overcome the entropic folding barrier. Our study highlights the subtle tuning and co-dependency in the folding of a discontinuous multi-domain protein.

Line Shape Analysis of 19F NMR-Monitored Chemical Denaturation of a Fold-Switching Protein RfaH Reveals Its Slow Folding Dynamics

The recent discovery of metamorphic proteins, which can switch between multiple conformations under native conditions, has challenged the well-established one sequence-one structure paradigm of protein folding. This is exemplified in the C-terminal domain of the multidomain transcription factor RfaH, which converts from an α-helical coiled-coil conformation in its autoinhibited state to a β-barrel conformation upon activation. Here, we use multisite line shape analysis of 19F NMR-monitored equilibrium chemical denaturation measurements of two 19F-labeled variants of full-length RfaH, to show that it folds/unfolds slowly on the NMR time scale, in an apparent all-or-none fashion at physiological pH and room temperature in the 3.3-4.8 M urea concentration range. The significant thermodynamic stability and slow unfolding rate (kinetic stability) are likely responsible for maintaining the closed autoinhibited state of RfaH, preventing spurious interactions with RNA polymerase (RNAP) when not functional. Our results provide a quantitative understanding of the folding-function relationship in the model fold-switching protein RfaH.

Enzyme function and evolution through the lens of bioinformatics

Enzymes have been shaped by evolution over billions of years to catalyse the chemical reactions that support life on earth. Dispersed in the literature, or organised in online databases, knowledge about enzymes can be structured in distinct dimensions, either related to their quality as biological macromolecules, such as their sequence and structure, or related to their chemical functions, such as the catalytic site, kinetics, mechanism, and overall reaction. The evolution of enzymes can only be understood when each of these dimensions is considered. In addition, many of the properties of enzymes only make sense in the light of evolution. We start this review by outlining the main paradigms of enzyme evolution, including gene duplication and divergence, convergent evolution, and evolution by recombination of domains. In the second part, we overview the current collective knowledge about enzymes, as organised by different types of data and collected in several databases. We also highlight some increasingly powerful computational tools that can be used to close gaps in understanding, in particular for types of data that require laborious experimental protocols. We believe that recent advances in protein structure prediction will be a powerful catalyst for the prediction of binding, mechanism, and ultimately, chemical reactions. A comprehensive mapping of enzyme function and evolution may be attainable in the near future.

Proteome-wide structural analysis quantifies structural conservation across distant species

Traditional evolutionary biology research mainly relies on sequence information to infer evolutionary relationships between genes or proteins. In contrast, protein structural information has long been overlooked, although structures are more conserved and closely linked to the functions than the sequences. To address this gap, we conducted a proteome-wide structural analysis using experimental and computed protein structures for organisms from the three distinct domains, including Homo sapiens (eukarya), Escherichia coli (bacteria), and Methanocaldococcus jannaschii (archaea). We reveal the distribution of structural similarity and sequence identity at the genomic level and characterize the twilight zone, where signals obtained from sequence alignment are blurred and evolutionary relationships cannot be inferred unambiguously. We find that structurally similar homologous protein pairs in the twilight zone account for ∼0.004%-0.021% of all possible protein pair combinations, which translates to ∼8%-32% of the protein-coding genes, depending on the species under comparison. In addition, by comparing the structural homologs, we show that human proteins involved in the energy supply are more similar to their E. coli homologs, whereas proteins relating to the central dogma are more similar to their M. jannaschii homologs. We also identify a bacterial GPCR homolog in the E. coli proteome that displays distinctive domain architecture. Our results shed light on the characteristics of the twilight zone and the origin of different pathways from a protein structure perspective, highlighting an exciting new frontier in evolutionary biology.

Large-scale clustering of AlphaFold2 3D models shines light on the structure and function of proteins

Two recent studies exploited ultra-fast structural aligners and deep-learning approaches to cluster the protein structure space in the AlphaFold Database. Barrio-Hernandez et al.1 and Durairaj et al.2 uncovered fascinating new protein functions and structural features previously unknown.

Retracing the evolution of a modern periplasmic binding protein

Investigating the evolution of structural features in modern multidomain proteins helps to understand their immense diversity and functional versatility. The class of periplasmic binding proteins (PBPs) offers an opportunity to interrogate one of the main processes driving diversification: the duplication and fusion of protein sequences to generate new architectures. The symmetry of their two-lobed topology, their mechanism of binding, and the organization of their operon structure led to the hypothesis that PBPs arose through a duplication and fusion event of a single common ancestor. To investigate this claim, we set out to reverse the evolutionary process and recreate the structural equivalent of a single-lobed progenitor using ribose-binding protein (RBP) as our model. We found that this modern PBP can be deconstructed into its lobes, producing two proteins that represent possible progenitor halves. The isolated halves of RBP are well folded and monomeric proteins, albeit with a lower thermostability, and do not retain the original binding function. However, the two entities readily form a heterodimer in vitro and in-cell. The x-ray structure of the heterodimer closely resembles the parental protein. Moreover, the binding function is fully regained upon formation of the heterodimer with a ligand affinity similar to that observed in the modern RBP. This highlights how a duplication event could have given rise to a stable and functional PBP-like fold and provides insights into how more complex functional structures can evolve from simpler molecular components.

Computational analysis of regulatory regions in human protein kinases

Eukaryotic proteins often feature modular domain structures comprising globular domains that are connected by linker regions and intrinsically disordered regions that may contain important functional motifs. The intramolecular interactions of globular domains and nonglobular regions can play critical roles in different aspects of protein function. However, studying these interactions and their regulatory roles can be challenging due to the flexibility of nonglobular regions, the long insertions separating interacting modules, and the transient nature of some interactions. Obtaining the experimental structures of multiple domains and functional regions is more difficult than determining the structures of individual globular domains. High-quality structural models generated by AlphaFold offer a unique opportunity to study intramolecular interactions in eukaryotic proteins. In this study, we systematically explored intramolecular interactions between human protein kinase domains (KDs) and potential regulatory regions, including globular domains, N- and C-terminal tails, long insertions, and distal nonglobular regions. Our analysis identified intramolecular interactions between human KDs and 35 different types of globular domains, exhibiting a variety of interaction modes that could contribute to orthosteric or allosteric regulation of kinase activity. We also identified prevalent interactions between human KDs and their flanking regions (N- and C-terminal tails). These interactions exhibit group-specific characteristics and can vary within each specific kinase group. Although long-range interactions between KDs and nonglobular regions are relatively rare, structural details of these interactions offer new insights into the regulation mechanisms of several kinases, such as HASPIN, MAPK7, MAPK15, and SIK1B.

Genome-wide analysis of KIX gene family for organ size regulation in soybean (Glycine max L.)

The KIX domain, conserved among various nuclear and co-activator factors, acts as a binding site that interacts with other transcriptional activators and co-activators, playing a crucial role in gene expression regulation. In plants, the KIX domain is involved in plant hormone signaling, stress response regulation, cell cycle control, and differentiation, indicating its potential relevance to crop productivity. This study aims to identify and characterize KIX domains within the soybean (Glycine max L.) genome to predict their potential role in improving crop productivity. The conservation and evolutionary history of the KIX domains were explored in 59 plant species, confirming the presence of the KIX domains in diverse plants. Specifically, 13 KIX domains were identified within the soybean genome and classified into four main groups, namely GmKIX8/9, GmMED15, GmHAC, and GmRECQL, through sequence alignment, structural analysis, and phylogenetic tree construction. Association analysis was performed between KIX domain haplotypes and soybean seed-related agronomic traits using re-sequencing data from a core collection of 422 accessions. The results revealed correlations between SNP variations observed in GmKIX8-3 and GmMED15-4 and soybean seed phenotypic traits. Additionally, transcriptome analysis confirmed significant expression of the KIX domains during the early stages of soybean seed development. This study provides the first characterization of the structural, expression, genomic haplotype, and molecular features of the KIX domain in soybean, offering a foundation for functional analysis of the KIX domain in soybean and other plants.

Ubiquitin-Dependent and Independent Proteasomal Degradation in Host-Pathogen Interactions

Ubiquitin, a small protein, is well known for tagging target proteins through a cascade of enzymatic reactions that lead to protein degradation. The ubiquitin tag, apart from its signaling role, is paramount in destabilizing the modified protein. Here, we explore the complex role of ubiquitin-mediated protein destabilization in the intricate proteolysis process by the 26S proteasome. In addition, the significance of the so-called ubiquitin-independent pathway and the role of the 20S proteasome are considered. Next, we discuss the ubiquitin-proteasome system's interplay with pathogenic microorganisms and how the microorganisms manipulate this system to establish infection by a range of elaborate pathways to evade or counteract host responses. Finally, we focus on the mechanisms that rely either on (i) hijacking the host and on delivering pathogenic E3 ligases and deubiquitinases that promote the degradation of host proteins, or (ii) counteracting host responses through the stabilization of pathogenic effector proteins.

Thermodynamic and kinetic analysis of the LAO binding protein and its isolated domains reveal non-additivity in stability, folding and function

Substrate-binding proteins (SBPs) are used by organisms from the three domains of life for transport and signalling. SBPs are composed of two domains that collectively trap ligands with high affinity and selectivity. To explore the role of the domains and the integrity of the hinge region between them in the function and conformation of SBPs, here, we describe the ligand binding, conformational stability and folding kinetics of the Lysine Arginine Ornithine (LAO) binding protein from Salmonella thiphimurium and constructs corresponding to its two independent domains. LAO is a class II SBP formed by a continuous and a discontinuous domain. Contrary to the expected behaviour based on their connectivity, the discontinuous domain shows a stable native-like structure that binds l-arginine with moderate affinity, whereas the continuous domain is barely stable and shows no detectable ligand binding. Regarding folding kinetics, studies of the entire protein revealed the presence of at least two intermediates. While the unfolding and refolding of the continuous domain exhibited only a single intermediate and simpler and faster kinetics than LAO, the folding mechanism of the discontinuous domain was complex and involved multiple intermediates. These findings suggest that in the complete protein the continuous domain nucleates folding and that its presence funnels the folding of the discontinuous domain avoiding nonproductive interactions. The strong dependence of the function, stability and folding pathway of the lobes on their covalent association is most likely the result of the coevolution of both domains as a single unit.

The New Functional Hybrid Chaperone Protein ADGroEL-SacSm

The creation of new proteins by combining natural domains is a commonly used technique in protein engineering. In this work, we have tested the possibilities and limitations of using circular homo-oligomeric Sm-like proteins as a basis for attaching other domains. Attachment to such a stable base should bring target domains together and keep them in the correct mutual orientation. We chose a circular homoheptameric Sm-like protein from Sulfolobus acidocaldarius as a stable backbone and the apical domain of the GroEL chaperone protein as the domain of study. This domain by itself, separated from the rest of the GroEL molecule, does not form an oligomeric ring. In our design, the hyperstable SacSm held the seven ADGroELs together and forced them to oligomerize. The designed hybrid protein was obtained and studied with various physical and chemical methods. Stepwise assembly and self-organization of this protein have been shown. First, the SacSm base was assembled, and then ADGroEL was folded on it. Functional testing showed that the obtained fusion protein was able to bind the same non-native proteins as the full-length GroEL chaperone. It also reduced the aggregation of a number of proteins when they were heated, which confirms its chaperone activity. Thus, the engineering path we chose made it possible to create an efficient thermostable chaperone. The result obtained shows the productivity of the way we chose for the creation and stabilization of oligomeric proteins.

An Initial Proteomic Analysis of Biogas-Related Metabolism of Euryarchaeota Consortia in Sediments from the Santiago River, México

In this paper, sediments from the Santiago River were characterized to look for an alternative source of inoculum for biogas production. A proteomic analysis of methane-processing archaea present in these sediments was carried out. The Euryarchaeota superkingdom of archaea is responsible for methane production and methane assimilation in the environment. The Santiago River is a major river in México with great pollution and exceeded recovery capacity. Its sediments could contain nutrients and the anaerobic conditions for optimal growth of Euryarchaeota consortia. Batch bioreactor experiments were performed, and a proteomic analysis was conducted with current database information. The maximum biogas production was 266 NmL·L^-1·g VS^-1, with 33.34% of methane, and for proteomics, 3206 proteins were detected from 303 species of 69 genera. Most of them are metabolically versatile members of the genera Methanosarcina and Methanosarcinales, both with 934 and 260 proteins, respectively. These results showed a diverse euryarcheotic species with high potential to methane production. Although related proteins were found and could be feeding this metabolism through the methanol and acetyl-CoA pathways, the quality obtained from the biogas suggests that this metabolism is not the main one in carbon use, possibly the sum of several conditions including growth conditions and the pollution present in these sediments.

Polyprotein synthesis: a journey from the traditional pre-translational method to modern post-translational approaches for single-molecule force spectroscopy

Polyproteins, an array of protein units of similar or differential functions in tandem, have been extensively utilized by organisms, unicellular or multicellular, as concentrators of the myriad of molecular activities. Most eukaryotic proteins, two-thirds in unicellular organisms, and more than 80% in metazoans, are polyproteins. Although the use of polyproteins continues to evolve in nature, our understanding of the structure-function-property of polyproteins is still limited. Cumbersome recombinant strategies and the lack of convenient in vitro synthetic routes of polyproteins have been rate-determining factors in the progress. However, in this review we have discussed the revolutionary journey of polyprotein synthesis with a major focus on surface-based structure-function-property studies, especially using force spectroscopy at the single-molecule level.

Domain-wise dissection of thermal stability enhancement in multidomain proteins

Stability is critical for the proper functioning of all proteins. Optimization of protein thermostability is a key step in the development of industrial enzymes and biologics. Herein, we demonstrate that multidomain proteins can be stabilized significantly using domain-based engineering followed by the recombination of the optimized domains. Domain-level analysis of designed protein variants with similar structures but different thermal profiles showed that the independent enhancement of the thermostability of a constituent domain improves the overall stability of the whole multidomain protein. The crystal structure and AlphaFold-predicted model of the designed proteins via domain-recombination provided a molecular explanation for domain-based stepwise stabilization. Our study suggests that domain-based modular engineering can minimize the sequence space for calculations in computational design and experimental errors, thereby offering useful guidance for multidomain protein engineering.

How protein topology controls allosteric regulations

Allostery is an important regulatory mechanism of protein functions. Among allosteric proteins, certain protein structure types are more observed. However, how allosteric regulation depends on protein topology remains elusive. In this study, we extracted protein topology graphs at the fold level and found that known allosteric proteins mainly contain multiple domains or subunits and allosteric sites reside more often between two or more domains of the same fold type. Only a small fraction of fold-fold combinations are observed in allosteric proteins, and homo-fold-fold combinations dominate. These analyses imply that the locations of allosteric sites including cryptic ones depend on protein topology. We further developed TopoAlloSite, a novel method that uses the kernel support vector machine to predict the location of allosteric sites on the overall protein topology based on the subgraph-matching kernel. TopoAlloSite successfully predicted known cryptic allosteric sites in several allosteric proteins like phosphopantothenoylcysteine synthetase, spermidine synthase, and sirtuin 6, demonstrating its power in identifying cryptic allosteric sites without performing long molecular dynamics simulations or large-scale experimental screening. Our study demonstrates that protein topology largely determines how its function can be allosterically regulated, which can be used to find new druggable targets and locate potential binding sites for rational allosteric drug design.

Understanding the molecular basis of folding cooperativity through a comparative analysis of a multidomain protein and its isolated domains

Although cooperativity is a well-established and general property of folding, our current understanding of this feature in multidomain folding is still relatively limited. In fact, there are contrasting results indicating that the constituent domains of a multidomain protein may either fold independently on each other or exhibit interdependent supradomain phenomena. To address this issue, here we present the comparative analysis of the folding of a tandem repeat protein, comprising two contiguous PDZ domains, in comparison to that of its isolated constituent domains. By analyzing in detail the equilibrium and kinetics of folding at different experimental conditions, we demonstrate that despite each of the PDZ domains in isolation being capable of independent folding, at variance with previously characterized PDZ tandem repeats, the full-length construct folds and unfolds as a single cooperative unit. By exploiting quantitatively, the comparison of the folding of the tandem repeat to those observed for its constituent domains, as well as by characterizing a truncated variant lacking a short autoinhibitory segment, we successfully rationalize the molecular basis of the observed cooperativity and attempt to infer some general conclusions for multidomain systems.

The role of single protein elasticity in mechanobiology

In addition to biochemical signals and genetic considerations, mechanical forces are rapidly emerging as a master regulator of human physiology. Yet the molecular mechanisms that regulate force-induced functionalities across a wide range of scales, encompassing the cell, tissue or organ levels, are comparatively not so well understood. With the advent, development and refining of single molecule nanomechanical techniques, enabling to exquisitely probe the conformational dynamics of individual proteins under the effect of a calibrated force, we have begun to acquire a comprehensive knowledge on the rich plethora of physicochemical principles that regulate the elasticity of single proteins. Here we review the major advances underpinning our current understanding of how the elasticity of single proteins regulates mechanosensing and mechanotransduction. We discuss the present limitations and future challenges of such a prolific and burgeoning field.

Direct observation of Hsp90-induced compaction in a protein chain

The chaperone heat shock protein 90 (Hsp90) is well known to undergo important conformational changes, which depend on nucleotide and substrate interactions. Conversely, how the conformations of its unstable and disordered substrates are affected by Hsp90 is difficult to address experimentally yet is central to its function. Here, using optical tweezers, we find that Hsp90 promotes local contractions in unfolded chains that drive their global compaction down to dimensions of folded states. This compaction has a gradual nature while showing small steps, is stimulated by ATP, and performs mechanical work against counteracting forces that expand the chain dimensions. The Hsp90 interactions suppress the formation of larger-scale folded, misfolded, and aggregated structures. The observations support a model in which Hsp90 alters client conformations directly by promoting local intra-chain interactions while suppressing distant ones. We conjecture that chain compaction may be central to how Hsp90 protects unstable clients and cooperates with Hsp70.

Probing Heterogeneous Folding Pathways of DNA Origami Self-Assembly at the Molecular Level with Atomic Force Microscopy

A myriad of DNA origami nanostructures have been demonstrated in various intriguing applications. In pursuit of facile yet high-yield synthesis, the mechanisms underlying DNA origami folding need to be resolved. Here, we visualize the folding processes of several multidomain DNA origami structures under ambient annealing conditions in solution using atomic force microscopy with submolecular resolution. We reveal the coexistence of diverse transitional structures that might result in the same prescribed products. Based on the experimental observations and the simulation of the energy landscapes, we propose the heterogeneity of the folding pathways of multidomain DNA origami structures. Our findings may contribute to understanding the high-yield folding mechanism of DNA origami.

Characterisation of HOIP RBR E3 ligase conformational dynamics using integrative modelling

Multidomain proteins composed of individual domains connected by flexible linkers pose a challenge for structural studies due to their intrinsic conformational dynamics. Integrated modelling approaches provide a means to characterise protein flexibility by combining experimental measurements with molecular simulations. In this study, we characterise the conformational dynamics of the catalytic RBR domain of the E3 ubiquitin ligase HOIP, which regulates immune and inflammatory signalling pathways. Specifically, we combine small angle X-ray scattering experiments and molecular dynamics simulations to generate weighted conformational ensembles of the HOIP RBR domain using two different approaches based on maximum parsimony and maximum entropy principles. Both methods provide optimised ensembles that are instrumental in rationalising observed differences between SAXS-based solution studies and available crystal structures and highlight the importance of interdomain linker flexibility.

Fuzzy supertertiary interactions within PSD-95 enable ligand binding

The scaffold protein PSD-95 links postsynaptic receptors to sites of presynaptic neurotransmitter release. Flexible linkers between folded domains in PSD-95 enable a dynamic supertertiary structure. Interdomain interactions within the PSG supramodule, formed by PDZ3, SH3, and Guanylate Kinase domains, regulate PSD-95 activity. Here we combined discrete molecular dynamics and single molecule Förster resonance energy transfer (FRET) to characterize the PSG supramodule, with time resolution spanning picoseconds to seconds. We used a FRET network to measure distances in full-length PSD-95 and model the conformational ensemble. We found that PDZ3 samples two conformational basins, which we confirmed with disulfide mapping. To understand effects on activity, we measured binding of the synaptic adhesion protein neuroligin. We found that PSD-95 bound neuroligin well at physiological pH while truncated PDZ3 bound poorly. Our hybrid structural models reveal how the supertertiary context of PDZ3 enables recognition of this critical synaptic ligand.

I-TASSER-MTD: a deep-learning-based platform for multi-domain protein structure and function prediction

Most proteins in cells are composed of multiple folding units (or domains) to perform complex functions in a cooperative manner. Relative to the rapid progress in single-domain structure prediction, there are few effective tools available for multi-domain protein structure assembly, mainly due to the complexity of modeling multi-domain proteins, which involves higher degrees of freedom in domain-orientation space and various levels of continuous and discontinuous domain assembly and linker refinement. To meet the challenge and the high demand of the community, we developed I-TASSER-MTD to model the structures and functions of multi-domain proteins through a progressive protocol that combines sequence-based domain parsing, single-domain structure folding, inter-domain structure assembly and structure-based function annotation in a fully automated pipeline. Advanced deep-learning models have been incorporated into each of the steps to enhance both the domain modeling and inter-domain assembly accuracy. The protocol allows for the incorporation of experimental cross-linking data and cryo-electron microscopy density maps to guide the multi-domain structure assembly simulations. I-TASSER-MTD is built on I-TASSER but substantially extends its ability and accuracy in modeling large multi-domain protein structures and provides meaningful functional insights for the targets at both the domain- and full-chain levels from the amino acid sequence alone.

The Wako-Saitô-Muñoz-Eaton Model for Predicting Protein Folding and Dynamics

Despite the recent advances in the prediction of protein structures by deep neutral networks, the elucidation of protein-folding mechanisms remains challenging. A promising theory for describing protein folding is a coarse-grained statistical mechanical model called the Wako-Saitô-Muñoz-Eaton (WSME) model. The model can calculate the free-energy landscapes of proteins based on a three-dimensional structure with low computational complexity, thereby providing a comprehensive understanding of the folding pathways and the structure and stability of the intermediates and transition states involved in the folding reaction. In this review, we summarize previous and recent studies on protein folding and dynamics performed using the WSME model and discuss future challenges and prospects. The WSME model successfully predicted the folding mechanisms of small single-domain proteins and the effects of amino-acid substitutions on protein stability and folding in a manner that was consistent with experimental results. Furthermore, extended versions of the WSME model were applied to predict the folding mechanisms of multi-domain proteins and the conformational changes associated with protein function. Thus, the WSME model may contribute significantly to solving the protein-folding problem and is expected to be useful for predicting protein folding, stability, and dynamics in basic research and in industrial and medical applications.

Industrial light at the end of the Iron-containing (group III) alcohol dehydrogenase tunnel

There are three prominent alcohol dehydrogenases superfamilies: Short-chain, Medium-chain, and Iron-containing alcohol dehydrogenases (FeADHs). Many members are valuable catalysts for producing industrially relevant products such as Active pharmaceutical Intermediates, Chiral synthons, Biopolymers, Biofuels and secondary metabolites. However, FeADHs are the least explored enzymes among the superfamilies for commercial tenacities. They portray a conserved structure having a 'tunnel-like' cofactor and substrate binding site with particular functions, despite representing high sequence diversity. Interestingly, phylogenetic analysis demarcates enzymes catalyzing distinct native substrates where closely related clades convert similar molecules. Further, homologs from various mesophilic and thermophilic microbes have been explored for designing a solvent and temperature resistant enzyme for industrial purposes. The review explores different Iron-containing alcohol dehydrogenases potential engineering of the enzymes and substrates helpful in manufacturing commercial products. This article is protected by copyright. All rights reserved.

Simulating domain architecture evolution

Motivation

Simulation is an essential technique for generating biomolecular data with a ‘known’ history for use in validating phylogenetic inference and other evolutionary methods. On longer time scales, simulation supports investigations of equilibrium behavior and provides a formal framework for testing competing evolutionary hypotheses. Twenty years of molecular evolution research have produced a rich repertoire of simulation methods. However, current models do not capture the stringent constraints acting on the domain insertions, duplications, and deletions by which multidomain architectures evolve. Although these processes have the potential to generate any combination of domains, only a tiny fraction of possible domain combinations are observed in nature. Modeling these stringent constraints on domain order and co-occurrence is a fundamental challenge in domain architecture simulation that does not arise with sequence and gene family simulation.

Results

Here, we introduce a stochastic model of domain architecture evolution to simulate evolutionary trajectories that reflect the constraints on domain order and co-occurrence observed in nature. This framework is implemented in a novel domain architecture simulator, DomArchov, using the Metropolis–Hastings algorithm with data-driven transition probabilities. The use of a data-driven event module enables quick and easy redeployment of the simulator for use in different taxonomic and protein function contexts. Using empirical evaluation with metazoan datasets, we demonstrate that domain architectures simulated by DomArchov recapitulate properties of genuine domain architectures that reflect the constraints on domain order and adjacency seen in nature. This work expands the realm of evolutionary processes that are amenable to simulation.

Availability and implementation

DomArchov is written in Python 3 and is available at http://www.cs.cmu.edu/~durand/DomArchov. The data underlying this article are available via the same link.

Supplementary information

Supplementary data are available at Bioinformatics online.

High-Pressure, Low-Temperature Induced Unfolding and Aggregation of Monoclonal Antibodies: Role of the Fc and Fab Fragments

The effects of high pressure and low temperature on the stability of two different monoclonal antibodies (MAbs) were examined in this work. Fluorescence and small-angle neutron scattering were used to monitor the in situ effects of pressure to infer shifts in tertiary structure and characterize aggregation prone intermediates. Partial unfolding was observed for both MAbs, to different extents, under a range of pressure/temperature conditions. Fourier transform infrared spectroscopy was also used to monitor ex situ changes in secondary structure. Preservation of native secondary structure after incubation at elevated pressures and subzero ° C temperatures was independent of the extent of tertiary unfolding and reversibility. Several combinations of pressure and temperature were also used to discern the respective contributions of the isolated Ab fragments (Fab and Fc) to unfolding and aggregation. The fragments for each antibody showed significantly different partial unfolding profiles and reversibility. There was not a simple correlation between stability of the full MAb and either the Fc or Fab fragment stabilities across all cases, demonstrating a complex relationship to full MAb unfolding and aggregation behavior. That notwithstanding, the combined use of spectroscopic and scattering techniques provides insights into MAb conformational stability and hysteresis in high-pressure, low-temperature environments.

Cooperativity and Folding Kinetics in a Multidomain Protein with Interwoven Chain Topology

Although a large percentage of eukaryotic proteomes consist of proteins with multiple domains, not much is known about their assembly mechanism, especially those with intricate native state architectures. Some have a complex topology in which the structural elements along the sequence are interwoven in such a manner that the domains cannot be separated by cutting at any location along the sequence. Such proteins are multiply connected multidomain proteins (MMPs) with the three-domain (NMP, LID, and CORE) phosphotransferase enzyme adenylate kinase (ADK) being an example. We devised a coarse-grained model to simulate ADK folding initiated by changing either the temperature or guanidinium chloride (GdmCl) concentration. The simulations reproduce the experimentally measured melting temperatures (associated with two equilibrium transitions), FRET efficiency as a function of GdmCl concentration, and the folding times quantitatively. Although the NMP domain orders independently, cooperative interactions between the LID and the CORE domains are required for complete assembly of the enzyme. Kinetic simulations show that, on the collapse time scale, multiple interconnected metastable states are populated, attesting to the folding heterogeneity. The network of kinetically connected states reveals that the CORE domain folds only after the NMP and LID domains, reflecting the interwoven nature of the chain topology.

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.

Protein folding in vitro and in the cell: From a solitary journey to a team effort

Correct protein folding is essential for the health and function of living organisms. Yet, it is not well understood how unfolded proteins reach their native state and avoid aggregation, especially within the cellular milieu. Some proteins, especially small, single-domain and apparent two-state folders, successfully attain their native state upon dilution from denaturant. Yet, many more proteins undergo misfolding and aggregation during this process, in a concentration-dependent fashion. Once formed, native and aggregated states are often kinetically trapped relative to each other. Hence, the early stages of protein life are absolutely critical for proper kinetic channeling to the folded state and for long-term solubility and function. This review summarizes current knowledge on protein folding/aggregation mechanisms in buffered solution and within the bacterial cell, highlighting early stages. Remarkably, teamwork between nascent chain, ribosome, trigger factor and Hsp70 molecular chaperones enables all proteins to overcome aggregation propensities and reach a long-lived bioactive state.

Structural and Kinetic Views of Molecular Chaperones in Multidomain Protein Folding

Despite recent developments in protein structure prediction, the process of the structure formation, folding, remains poorly understood. Notably, folding of multidomain proteins, which involves multiple steps of segmental folding, is one of the biggest questions in protein science. Multidomain protein folding often requires the assistance of molecular chaperones. Molecular chaperones promote or delay the folding of the client protein, but the detailed mechanisms are still unclear. This review summarizes the findings of biophysical and structural studies on the mechanism of multidomain protein folding mediated by molecular chaperones and explains how molecular chaperones recognize the client proteins and alter their folding properties. Furthermore, we introduce several recent studies that describe the concept of kinetics-activity relationships to explain the mechanism of functional diversity of molecular chaperones.

The Importance of Gene Duplication and Domain Repeat Expansion for the Function and Evolution of Fertilization Proteins

The process of gene duplication followed by gene loss or evolution of new functions has been studied extensively, yet the role gene duplication plays in the function and evolution of fertilization proteins is underappreciated. Gene duplication is observed in many fertilization protein families including Izumo, DCST, ZP, and the TFP superfamily. Molecules mediating fertilization are part of larger gene families expressed in a variety of tissues, but gene duplication followed by structural modifications has often facilitated their cooption into a fertilization function. Repeat expansions of functional domains within a gene also provide opportunities for the evolution of novel fertilization protein. ZP proteins with domain repeat expansions are linked to species-specificity in fertilization and TFP proteins that experienced domain duplications were coopted into a novel sperm function. This review outlines the importance of gene duplications and repeat domain expansions in the evolution of fertilization proteins.

Structure of putative tumor suppressor ALDH1L1

Putative tumor suppressor ALDH1L1, the product of natural fusion of three unrelated genes, regulates folate metabolism by catalyzing NADP+-dependent conversion of 10-formyltetrahydrofolate to tetrahydrofolate and CO2. Cryo-EM structures of tetrameric rat ALDH1L1 revealed the architecture and functional domain interactions of this complex enzyme. Highly mobile N-terminal domains, which remove formyl from 10-formyltetrahydrofolate, undergo multiple transient inter-domain interactions. The C-terminal aldehyde dehydrogenase domains, which convert formyl to CO2, form unusually large interfaces with the intermediate domains, homologs of acyl/peptidyl carrier proteins (A/PCPs), which transfer the formyl group between the catalytic domains. The 4'-phosphopantetheine arm of the intermediate domain is fully extended and reaches deep into the catalytic pocket of the C-terminal domain. Remarkably, the tetrameric state of ALDH1L1 is indispensable for catalysis because the intermediate domain transfers formyl between the catalytic domains of different protomers. These findings emphasize the versatility of A/PCPs in complex, highly dynamic enzymatic systems.

Blocking PSD95-PDZ3's amyloidogenesis through point mutations that inhibit high-temperature reversible oligomerization (RO)

The third PDZ domain of postsynaptic density protein 95 (PSD95-PDZ3; 11 kDa, 103 residues) has a propensity to form amyloid fibrils at high temperatures. At neutral pH, wild type PDZ3 is natively folded, but it exhibits a peculiar three-state thermal unfolding with a reversible oligomerization (RO) equilibrium at high temperatures, which is uncharacteristic in the unfolding of small globular protein as PDZ3 is. Here, we examined RO's role in PDZ3's amyloidogenesis at high-temperature using two variants (F340A and L342A) that suppress the high-temperature RO and five single-alanine-mutated variants, where we mutated surface-exposed hydrophobic residues to alanine. Circular Dichroism (CD), Analytical Ultracentrifuge (AUC), and other spectroscopic measurements confirmed the retention of the native structure at ambient temperature. Differential Scanning Calorimetry (DSC) was used to assess the presence or absence of the high-temperature RO, and the amyloidogenicity of the variants was measured by Thioflavin T (ThT) fluorescence and Transmission Electron Microscopy (TEM). By comparing the fraction of RO and the ThT signal, we found that mutations that suppressed the high-temperature RO strongly inhibited amyloidogenesis. On the other hand, all variants forming RO also formed amyloids under the same conditions as the wild type PDZ3.

Slowest-first protein translation scheme: Structural asymmetry and co-translational folding

Proteins are translated from the N to the C terminus, raising the basic question of how this innate directionality affects their evolution. To explore this question, we analyze 16,200 structures from the Protein Data Bank (PDB). We find remarkable enrichment of α helices at the C terminus and β strands at the N terminus. Furthermore, this α-β asymmetry correlates with sequence length and contact order, both determinants of folding rate, hinting at possible links to co-translational folding (CTF). Hence, we propose the "slowest-first" scheme, whereby protein sequences evolved structural asymmetry to accelerate CTF: the slowest of the cooperatively folding segments are positioned near the N terminus so they have more time to fold during translation. A phenomenological model predicts that CTF can be accelerated by asymmetry in folding rate, up to double the rate, when folding time is commensurate with translation time; analysis of the PDB predicts that structural asymmetry is indeed maximal in this regime. This correspondence is greater in prokaryotes, which generally require faster protein production. Altogether, this indicates that accelerating CTF is a substantial evolutionary force whose interplay with stability and functionality is encoded in secondary structure asymmetry.

A unique NLRC4 receptor from echinoderms mediates Vibrio phagocytosis via rearrangement of the cytoskeleton and polymerization of F-actin

Many members of the nucleotide-binding and oligomerization domain (NACHT)- and leucine-rich-repeat-containing protein (NLR) family play crucial roles in pathogen recognition and innate immune response regulation. In our previous work, a unique and Vibrio splendidus-inducible NLRC4 receptor comprising Ig and NACHT domains was identified from the sea cucumber Apostichopus japonicus, and this receptor lacked the CARD and LRR domains that are typical of common cytoplasmic NLRs. To better understand the functional role of AjNLRC4, we confirmed that AjNLRC4 was a bona fide membrane PRR with two transmembrane structures. AjNLRC4 was able to directly bind microbes and polysaccharides via its extracellular Ig domain and agglutinate a variety of microbes in a Ca²⁺-dependent manner. Knockdown of AjNLRC4 by RNA interference and blockade of AjNLRC4 by antibodies in coelomocytes both could significantly inhibit the phagocytic activity and elimination of V. splendidus. Conversely, overexpression of AjNLRC4 enhanced the phagocytic activity of V. splendidus, and this effect could be specifically blocked by treatment with the actin-mediated endocytosis inhibitor cytochalasin D but not other endocytosis inhibitors. Moreover, AjNLRC4-mediated phagocytic activity was dependent on the interaction between the intracellular domain of AjNLRC4 and the β-actin protein and further regulated the Arp2/3 complex to mediate the rearrangement of the cytoskeleton and the polymerization of F-actin. V. splendidus was found to be colocalized with lysosomes in coelomocytes, and the bacterial quantities were increased after injection of chloroquine, a lysosome inhibitor. Collectively, these results suggested that AjNLRC4 served as a novel membrane PRR in mediating coelomocyte phagocytosis and further clearing intracellular Vibrio through the AjNLRC4-β-actin-Arp2/3 complex-lysosome pathway.

Vibrio splendidus is ubiquitously present in marine environments and in or on many aquaculture species and is considered to be an important opportunistic pathogen that has caused serious economic losses to the aquaculture industry worldwide. Phagocytosis is the first step of pathogen clearance and is triggered by specific interactions between host pattern recognition receptors (PRRs) and pathogen-associated molecular patterns (PAMPs) from invasive bacteria. However, the mechanism that underlies receptor-mediated V. splendidus phagocytosis is poorly understood. In this study, an atypical AjNLRC4 receptor without LRR and CARD domains was found to serve as the membrane receptor for V. splendidus, not the common cytoplasmic NLRs. The Ig domain of AjNLRC4 is replaced with a conventional LRR domain to bind V. splendidus, and the intracellular domain of AjNLRC4 specifically interacts with β-actin to mediate V. splendidus endocytosis in an actin-dependent manner. Endocytic V. splendidus is ultimately degraded in phagolysosomes. Our findings will contribute to the development of novel strategies for treating V. splendidus infection by modulating the actin-dependent endocytosis pathway.

Probing ion channel functional architecture and domain recombination compatibility by massively parallel domain insertion profiling

Protein domains are the basic units of protein structure and function. Comparative analysis of genomes and proteomes showed that domain recombination is a main driver of multidomain protein functional diversification and some of the constraining genomic mechanisms are known. Much less is known about biophysical mechanisms that determine whether protein domains can be combined into viable protein folds. Here, we use massively parallel insertional mutagenesis to determine compatibility of over 300,000 domain recombination variants of the Inward Rectifier K⁺ channel Kir2.1 with channel surface expression. Our data suggest that genomic and biophysical mechanisms acted in concert to favor gain of large, structured domain at protein termini during ion channel evolution. We use machine learning to build a quantitative biophysical model of domain compatibility in Kir2.1 that allows us to derive rudimentary rules for designing domain insertion variants that fold and traffic to the cell surface. Positional Kir2.1 responses to motif insertion clusters into distinct groups that correspond to contiguous structural regions of the channel with distinct biophysical properties tuned towards providing either folding stability or gating transitions. This suggests that insertional profiling is a high-throughput method to annotate function of ion channel structural regions.

Tracing protein and proteome history with chronologies and networks: folding recapitulates evolution

INTRODUCTION

While the origin and evolution of proteins remain mysterious, advances in evolutionary genomics and systems biology are facilitating the historical exploration of the structure, function and organization of proteins and proteomes. Molecular chronologies are series of time events describing the history of biological systems and subsystems and the rise of biological innovations. Together with time-varying networks, these chronologies provide a window into the past.

AREAS COVERED

Here, we review molecular chronologies and networks built with modern methods of phylogeny reconstruction. We discuss how chronologies of structural domain families uncover the explosive emergence of metabolism, the late rise of translation, the co-evolution of ribosomal proteins and rRNA, and the late development of the ribosomal exit tunnel; events that coincided with a tendency to shorten folding time. Evolving networks described the early emergence of domains and a late 'big bang' of domain combinations.

EXPERT OPINION

Two processes, folding and recruitment appear central to the evolutionary progression. The former increases protein persistence. The later fosters diversity. Chronologically, protein evolution mirrors folding by combining supersecondary structures into domains, developing translation machinery to facilitate folding speed and stability, and enhancing structural complexity by establishing long-distance interactions in novel structural and architectural designs.