Structural and evolutionary versatility in protein complexes with uneven stoichiometry

Lignin developmental patterns and Casparian strip as apoplastic barriers: A review

Lignin and Casparian strips are two essential components of plant cells that play critical roles in plant development regulate nutrients and water across the plants cell. Recent studies have extensively investigated lignin diversity and Casparian strip formation, providing valuable insights into plant physiology. This review presents the established lignin biosynthesis pathway, as well as the developmental patterns of lignin and Casparian strip and transcriptional network associated with Casparian strip formation. It describes the biochemical and genetic mechanisms that regulate lignin biosynthesis and deposition in different plants cell types and tissues. Additionally, the review highlights recent studies that have uncovered novel lignin biosynthesis genes and enzymatic pathways, expanding our understanding of lignin diversity. This review also discusses the developmental patterns of Casparian strip in roots and their role in regulating nutrient and water transport, focusing on recent genetic and molecular studies that have identified regulators of Casparian strip formation. Previous research has shown that lignin biosynthesis genes also play a role in Casparian strip formation, suggesting that these processes are interconnected. In conclusion, this comprehensive overview provides insights into the developmental patterns of lignin diversity and Casparian strip as apoplastic barriers. It also identifies future research directions, including the functional characterization of novel lignin biosynthesis genes and the identification of additional regulators of Casparian strip formation. Overall, this review enhances our understanding of the complex and interconnected processes that drive plant growth, pathogen defense, regulation and development.

Dominant negative variants and cotranslational assembly of macromolecular complexes

Pathogenic variants occurring in protein-coding regions underlie human genetic disease through various mechanisms. They can lead to a loss of function (LOF) such as in recessive conditions or in dominant conditions due to haploinsufficiency. Dominant-negative (DN) effects, counteracting the activity of the normal gene-product, and gain of function (GOF) are also mechanisms driving dominance. Here, I discuss a few papers on these specific mechanisms. In short, there is accumulating evidence pointing to differences between LOF versus non-LOF variants (DN and GOF). The latter are thought to have milder effects on protein structure and, as expected, DN variants are enriched at protein interfaces. This tendency to cluster in 3D space can help improve the ability of computational tools to predict the pathogenicity of DN variants, which is currently a challenging issue. More recent results support the hypothesis whereby cotranslational assembly of macromolecular complexes can buffer deleterious consequences of variants that would otherwise lead to DN effects (DNEs). Indeed, subunits the variants of which are responsible for DNEs tend to elude cotranslational assembly, thus poisoning complexes involving wild-type subunits. The constraints explaining why the buffering of DNEs is not universal require further investigation.

Co-evolution at protein-protein interfaces guides inference of stoichiometry of oligomeric protein complexes by de novo structure prediction

The quaternary structure with specific stoichiometry is pivotal to the specific function of protein complexes. However, determining the structure of many protein complexes experimentally remains a major bottleneck. Structural bioinformatics approaches, such as the deep learning algorithm Alphafold2-multimer (AF2-multimer), leverage the co-evolution of amino acids and sequence-structure relationships for accurate de novo structure and contact prediction. Pseudo-likelihood maximization direct coupling analysis (plmDCA) has been used to detect co-evolving residue pairs by statistical modeling. Here, we provide evidence that combining both methods can be used for de novo prediction of the quaternary structure and stoichiometry of a protein complex. We achieve this by augmenting the existing AF2-multimer confidence metrics with an interpretable score to identify the complex with an optimal fraction of native contacts of co-evolving residue pairs at intermolecular interfaces. We use this strategy to predict the quaternary structure and non-trivial stoichiometries of Bacillus subtilis spore germination protein complexes with unknown structures. Co-evolution at intermolecular interfaces may therefore synergize with AI-based de novo quaternary structure prediction of structurally uncharacterized bacterial protein complexes.

Molecular determinants of protein evolvability

The plethora of biological functions that sustain life is rooted in the remarkable evolvability of proteins. An emerging view highlights the importance of a protein's initial state in dictating evolutionary success. A deeper comprehension of the mechanisms that govern the evolvability of these initial states can provide invaluable insights into protein evolution. In this review, we describe several molecular determinants of protein evolvability, unveiled by experimental evolution and ancestral sequence reconstruction studies. We further discuss how genetic variation and epistasis can promote or constrain functional innovation and suggest putative underlying mechanisms. By establishing a clear framework for these determinants, we provide potential indicators enabling the forecast of suitable evolutionary starting points and delineate molecular mechanisms in need of deeper exploration.

Buffering of genetic dominance by allele-specific protein complex assembly

Protein complex assembly often occurs while subunits are being translated, resulting in complexes whose subunits were translated from the same mRNA in an allele-specific manner. It has thus been hypothesized that such cotranslational assembly may counter the assembly-mediated dominant-negative effect, whereby co-assembly of mutant and wild-type subunits "poisons" complex activity. Here, we show that cotranslationally assembling subunits are much less likely to be associated with autosomal dominant relative to recessive disorders, and that subunits with dominant-negative disease mutations are significantly depleted in cotranslational assembly compared to those associated with loss-of-function mutations. We also find that complexes with known dominant-negative effects tend to expose their interfaces late during translation, lessening the likelihood of cotranslational assembly. Finally, by combining complex properties with other features, we trained a computational model for predicting proteins likely to be associated with non-loss-of-function disease mechanisms, which we believe will be of considerable utility for protein variant interpretation.

TTT (Tel2-Tti1-Tti2) Complex, the Co-Chaperone of PIKKs and a Potential Target for Cancer Chemotherapy

The heterotrimeric Tel2-Tti1-Tti2 or TTT complex is essential for cell viability and highly observed in eukaryotes. As the co-chaperone of ATR, ATM, DNA-PKcs, mTOR, SMG1, and TRRAP, the phosphatidylinositol 3-kinase-related kinases (PIKKs) and a group of large proteins of 300-500 kDa, the TTT plays crucial roles in genome stability, cell proliferation, telomere maintenance, and aging. Most of the protein kinases in the kinome are targeted by co-chaperone Cdc37 for proper folding and stability. Like Cdc37, accumulating evidence has established the mechanism by which the TTT interacts with chaperone Hsp90 via R2TP (Rvb1-Rvb2-Tah1-Pih1) complex or other proteins for co-translational maturation of the PIKKs. Recent structural studies have revealed the α-solenoid structure of the TTT and its interactions with the R2TP complex, which shed new light on the co-chaperone mechanism and provide new research opportunities. A series of mutations of the TTT have been identified that cause disease syndrome with neurodevelopmental defects, and misregulation of the TTT has been shown to contribute to myeloma, colorectal, and non-small-cell lung cancers. Surprisingly, Tel2 in the TTT complex has recently been found to be a target of ivermectin, an antiparasitic drug that has been used by millions of patients. This discovery provides mechanistic insight into the anti-cancer effect of ivermectin and thus promotes the repurposing of this Nobel-prize-winning medicine for cancer chemotherapy. Here, we briefly review the discovery of the TTT complex, discuss the recent studies, and describe the perspectives for future investigation.

A Modular Scaffold for Controlling Transcriptional Activation via Homomeric Protein-Protein Interactions

Protein-protein interactions (PPIs) have been extensively utilized in synthetic biology to construct artificial gene networks. However, synthetic regulation of gene expression by PPIs in E. coli has largely relied upon repressors, and existing PPI-controlled transcriptional activators have generally been employed with heterodimeric interactions. Here we report a highly modular, PPI-dependent transcriptional activator, cCadC, that was designed to be compatible with homomeric interactions. We describe the process of engineering cCadC from a transmembrane protein into a soluble cytosolic regulator. We then show that gene transcription by cCadC can be controlled by homomeric PPIs and furthermore discriminates between dimeric and higher-order interactions. Finally, we demonstrate that cCadC activity can be placed under small molecule regulation using chemically induced dimerization or ligand dependent stabilization. This work should greatly expand the scope of PPIs that can be employed in artificial gene circuits in E. coli and complements the existing repertoire of tools for transcriptional regulation in synthetic biology.

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.

Are transient protein-protein interactions more dispensable?

Protein-protein interactions (PPIs) are key drivers of cell function and evolution. While it is widely assumed that most permanent PPIs are important for cellular function, it remains unclear whether transient PPIs are equally important. Here, we estimate and compare dispensable content among transient PPIs and permanent PPIs in human. Starting with a human reference interactome mapped by experiments, we construct a human structural interactome by building three-dimensional structural models for PPIs, and then distinguish transient PPIs from permanent PPIs using several structural and biophysical properties. We map common mutations from healthy individuals and disease-causing mutations onto the structural interactome, and perform structure-based calculations of the probabilities for common mutations (assumed to be neutral) and disease mutations (assumed to be mildly deleterious) to disrupt transient PPIs and permanent PPIs. Using Bayes' theorem we estimate that a similarly small fraction (<~20%) of both transient and permanent PPIs are completely dispensable, i.e., effectively neutral upon disruption. Hence, transient and permanent interactions are subject to similarly strong selective constraints in the human interactome.

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.

Studying the Structures of Relaxed and Fuzzy Interactions: The Diverse World of S100 Complexes

S100 proteins are small, dimeric, Ca²⁺-binding proteins of considerable interest due to their associations with cancer and rheumatic and neurodegenerative diseases. They control the functions of numerous proteins by forming protein–protein complexes with them. Several of these complexes were found to display “fuzzy” properties. Examining these highly flexible interactions, however, is a difficult task, especially from a structural biology point of view. Here, we summarize the available in vitro techniques that can be deployed to obtain structural information about these dynamic complexes. We also review the current state of knowledge about the structures of S100 complexes, focusing on their often-asymmetric nature.

Modeling and Structure Determination of Homo-Oligomeric Proteins: An Overview of Challenges and Current Approaches

Protein homo-oligomerization is a very common phenomenon, and approximately half of proteins form homo-oligomeric assemblies composed of identical subunits. The vast majority of such assemblies possess internal symmetry which can be either exploited to help or poses challenges during structure determination. Moreover, aspects of symmetry are critical in the modeling of protein homo-oligomers either by docking or by homology-based approaches. Here, we first provide a brief overview of the nature of protein homo-oligomerization. Next, we describe how the symmetry of homo-oligomers is addressed by crystallographic and non-crystallographic symmetry operations, and how biologically relevant intermolecular interactions can be deciphered from the ordered array of molecules within protein crystals. Additionally, we describe the most important aspects of protein homo-oligomerization in structure determination by NMR. Finally, we give an overview of approaches aimed at modeling homo-oligomers using computational methods that specifically address their internal symmetry and allow the incorporation of other experimental data as spatial restraints to achieve higher model reliability.

Heterodimer Formation of the Homodimeric ABC Transporter OpuA

Many proteins have a multimeric structure and are composed of two or more identical subunits. While this can be advantageous for the host organism, it can be a challenge when targeting specific residues in biochemical analyses. In vitro splitting and re-dimerization to circumvent this problem is a tedious process that requires stable proteins. We present an in vivo approach to transform homodimeric proteins into apparent heterodimers, which then can be purified using two-step affinity-tag purification. This opens the door to both practical applications such as smFRET to probe the conformational dynamics of homooligomeric proteins and fundamental research into the mechanism of protein multimerization, which is largely unexplored for membrane proteins. We show that expression conditions are key for the formation of heterodimers and that the order of the differential purification and reconstitution of the protein into nanodiscs is important for a functional ABC-transporter complex.

Analysing the yeast complexome-the Complex Portal rising to the challenge

The EMBL-EBI Complex Portal is a knowledgebase of macromolecular complexes providing persistent stable identifiers. Entries are linked to literature evidence and provide details of complex membership, function, structure and complex-specific Gene Ontology annotations. Data are freely available and downloadable in HUPO-PSI community standards and missing entries can be requested for curation. In collaboration with Saccharomyces Genome Database and UniProt, the yeast complexome, a compendium of all known heteromeric assemblies from the model organism Saccharomyces cerevisiae, was curated. This expansion of knowledge and scope has led to a 50% increase in curated complexes compared to the previously published dataset, CYC2008. The yeast complexome is used as a reference resource for the analysis of complexes from large-scale experiments. Our analysis showed that genes coding for proteins in complexes tend to have more genetic interactions, are co-expressed with more genes, are more multifunctional, localize more often in the nucleus, and are more often involved in nucleic acid-related metabolic processes and processes where large machineries are the predominant functional drivers. A comparison to genetic interactions showed that about 40% of expanded co-complex pairs also have genetic interactions, suggesting strong functional links between complex members.

The near-symmetry of protein oligomers: NMR-derived structures

The majority of oligomeric proteins form clusters which have rotational or dihedral symmetry. Despite the many advantages of symmetric packing, protein oligomers are only nearly symmetric, and the origin of this phenomenon is still in need to be fully explored. Here we apply near-symmetry analyses by the Continuous Symmetry Measures methodology of protein homomers to their natural state, namely their structures in solution. NMR-derived structural data serves us for that purpose. We find that symmetry deviations of proteins are by far higher in solution, compared to the crystalline state; that much of the symmetry distortion is due to amino acids along the interface between the subunits; that the distortions are mainly due to hydrophilic amino acids; and that distortive oligomerization processes such as the swap-domain mechanism can be identified by the symmetry analysis. Most of the analyses were carried out on distorted C2-symmetry dimers, but C3 and D2 cases were analyzed as well. Our NMR analysis supports the idea that the crystallographic B-factor represents non-classical crystals, in which different conformers pack in the crystal, perhaps from the conformers which the NMR analysis provides.

Keeping the Proportions of Protein Complex Components in Check

How do cells maintain relative proportions of protein complex components? Advances in quantitative, genome-wide measurements have begun to shed light onto the roles of protein synthesis and degradation in establishing the precise proportions in living cells: on the one hand, ribosome profiling studies indicate that proteins are already produced in the correct relative proportions. On the other hand, proteomic studies found that many complexes contain subunits that are made in excess and subsequently degraded. Here, we discuss these seemingly contradictory findings, emerging principles, and remaining open questions. We conclude that establishing precise protein levels involves both coordinated synthesis and post-translational fine-tuning via protein degradation.

Thermokinetic Analysis of Protein Subunit Exchange by Variable-Temperature Native Mass Spectrometry

Many protein complexes are assembled from a varying number of subunits, which are continuously exchanging with diverse time scales. This structural dynamics is considered to be important for many regulatory and sensory adaptation processes that occur in vivo. We have developed an accurate method for monitoring protein subunit exchange by using native electrospray ionization mass spectrometry (ESI-MS), exemplified here for an extremely stable Rad50 zinc hook (Hk) dimer assembly, Zn(Hk)₂. The method has two steps: appropriate protein/peptide mutation and native ESI-MS analysis using a variable-temperature sample inlet. In this work, two Hk mutants were produced, mixed with wild-type Hk, and measured at three different temperatures. A thermokinetic analysis of heterodimer formation allowed us to determine the enthalpy, entropy, and Gibbs free energy of activation for subunit exchange, showing that the reaction is slow and associated with a high enthalpic barrier, consistent with the exceptionally high stability of the Zn(Hk)₂ assembly.

The interdimeric interface controls function and stability of Ureaplasma urealiticum methionine S-adenosyltransferase

Methionine S-adenosyltransferases (MATs) are predominantly homotetramers, comprised of dimers of dimers. The larger, highly conserved intradimeric interface harbors two active sites, making the dimer the obligatory functional unit. However, functionality of the smaller, more diverged, and recently evolved interdimeric interface is largely unknown. Here, we show that the interdimeric interface of Ureaplasmaurealiticum MAT has evolved to control the catalytic activity and structural integrity of the homotetramer in response to product accumulation. When all four active sites are occupied with the product, S-adenosylmethionine (SAM), binding of four additional SAM molecules to the interdimeric interface prompts a ∼45° shift in the dimer orientation and a concomitant ∼60% increase in the interface area. This rearrangement inhibits the enzymatic activity by locking the flexible active site loops in a closed state and renders the tetramer resistant to proteolytic degradation. Our findings suggest that the interdimeric interface of MATs is subject to rapid evolutionary changes that tailor the molecular properties of the entire homotetramer to the specific needs of the organism.

Ligand Binding Site Structure Influences the Evolution of Protein Complex Function and Topology

It has been suggested that the evolution of protein complexes is significantly influenced by stochastic, non-adaptive processes. Using ligand binding as a proxy of function, we show that the structure of ligand-binding sites significantly influences the evolution of protein complexes. We show that homomers with multi-chain binding sites (MBSs) evolve new functions slower than monomers or other homomers, and those binding cofactors and metals have more conserved quaternary structure than other homomers. Moreover, the ligands and ligand-binding pockets of homologous MBS homomers are more similar than monomers and other homomers. Our results suggest strong evolutionary selection for quaternary structure in cofactor-binding MBS homomers, whereas neutral processes are more important in complexes with single-chain binding sites. They also have pharmacological implications, suggesting that complexes with single-chain binding sites are better targets for selective drugs, whereas MBS homomers are good candidates for broad-spectrum antibiotic and multitarget drug design.

•

Ligand binding site structure significantly influences protein function evolution

•

MBS homomers have more similar ligand binding pockets than monomers and other homomers

•

Cofactor and metal-binding MBS homomers have more conserved QS than other homomers

•

MBS homomers are promising targets for developing antibiotics and multitarget drugs

Evolution of interface binding strengths in simplified model of protein quaternary structure

The self-assembly of proteins into protein quaternary structures is of fundamental importance to many biological processes, and protein misassembly is responsible for a wide range of proteopathic diseases. In recent years, abstract lattice models of protein self-assembly have been used to simulate the evolution and assembly of protein quaternary structure, and to provide a tractable way to study the genotype-phenotype map of such systems. Here we generalize these models by representing the interfaces as mutable binary strings. This simple change enables us to model the evolution of interface strengths, interface symmetry, and deterministic assembly pathways. Using the generalized model we are able to reproduce two important results established for real protein complexes: The first is that protein assembly pathways are under evolutionary selection to minimize misassembly. The second is that the assembly pathway of a complex mirrors its evolutionary history, and that both can be derived from the relative strengths of interfaces. These results demonstrate that the generalized lattice model offers a powerful new idealized framework to facilitate the study of protein self-assembly processes and their evolution.

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

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

Many protein structures show a great deal of regularity. Even within single polypeptide chains, about 25% of proteins contain self-similar repeating structures, which can be organized in ring-like symmetric arrangements or linear open repeats. The repeats are often related, and thus comparing the sequence and structure of repeats can give an idea as to the early evolutionary history of a protein family. Additionally, the conservation and divergence of repeats can lead to insights about the function of the proteins. This work describes CE-Symm 2.0, a tool for the analysis of protein symmetry. The method automatically detects internal symmetry in protein structures and produces a multiple alignment of structural repeats. The algorithm is able to detect the geometric relationships between the repeats, including cyclic, dihedral, and polyhedral symmetries, translational repeats, and cases where multiple symmetry operators are applicable in a hierarchical manner. These complex relationships can then be visualized in a graphical interface as a complete structure, as a superposition of repeats, or as a multiple alignment of the protein sequence. CE-Symm 2.0 can be systematically used for the automatic detection of internal symmetry in protein structures, or as an interactive tool for the analysis of structural repeats.

Principles and characteristics of biological assemblies in experimentally determined protein structures

More than half of all structures in the PDB are assemblies of two or more proteins, including both homooligomers and heterooligomers. Structural information on these assemblies comes from X-ray crystallography, NMR, and cryo-EM spectroscopy. The correct assembly in an X-ray structure is often ambiguous, and computational methods have been developed to identify the most likely biologically relevant assembly based on physical properties of assemblies and sequence conservation in interfaces. Taking advantage of the large number of structures now available, some of the most recent methods have relied on similarity of interfaces and assemblies across structures of homologous proteins.

Investigation of protein quaternary structure via stoichiometry and symmetry information

The Protein Data Bank (PDB) is the single worldwide archive of experimentally-determined three-dimensional (3D) structures of proteins and nucleic acids. As of January 2017, the PDB housed more than 125,000 structures and was growing by more than 11,000 structures annually. Since the 3D structure of a protein is vital to understand the mechanisms of biological processes, diseases, and drug design, correct oligomeric assembly information is of critical importance. Unfortunately, the biologically relevant oligomeric form of a 3D structure is not directly obtainable by X-ray crystallography, whilst in solution methods (NMR or single particle EM) it is known from the experiment. Instead, this information may be provided by the PDB Depositor as metadata coming from additional experiments, be inferred by sequence-sequence comparisons with similar proteins of known oligomeric state, or predicted using software, such as PISA (Proteins, Interfaces, Structures and Assemblies) or EPPIC (Evolutionary Protein Protein Interface Classifier). Despite significant efforts by professional PDB Biocurators during data deposition, there remain a number of structures in the archive with incorrect quaternary structure descriptions (or annotations). Further investigation is, therefore, needed to evaluate the correctness of quaternary structure annotations. In this study, we aim to identify the most probable oligomeric states for proteins represented in the PDB. Our approach evaluated the performance of four independent prediction methods, including text mining of primary publications, inference from homologous protein structures, and two computational methods (PISA and EPPIC). Aggregating predictions to give consensus results outperformed all four of the independent prediction methods, yielding 83% correct, 9% wrong, and 8% inconclusive predictions, when tested with a well-curated benchmark dataset. We have developed a freely-available web-based tool to make this approach accessible to researchers and PDB Biocurators (http://quatstruct.rcsb.org/).

Automated evaluation of quaternary structures from protein crystals

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

X-ray diffraction experiments are the main experimental technique to reveal the detailed atomic 3-dimensional structure of proteins. In these experiments, proteins are packed into crystals, an environment that is far away from their native solution environment. Determining which parts of the structure reflect the protein’s state in the cell rather than being artifacts of the crystal environment can be a difficult task. How the different protein subunits assemble together in solution is known as the quaternary structure. Finding the correct quaternary structure is important both to understand protein oligomerization and for the understanding of protein-protein interactions at large. Here we present a new method to automatically determine the quaternary structure of proteins given their crystal structure. We provide a theoretical basis for properties that correct protein assemblies should possess, and provide a systematic evaluation of all possible assemblies according to these properties. The method provides a guidance to the experimental structural biologist as well as to structural bioinformaticians analyzing protein structures in bulk. Assemblies are provided for all proteins in the Protein Data Bank through a public website and database that is updated weekly as new structures are released.

Assessing the impact of the 4CL enzyme complex on the robustness of monolignol biosynthesis using metabolic pathway analysis

Lignin is a polymer present in the secondary cell walls of all vascular plants. It is a known barrier to pulping and the extraction of high-energy sugars from cellulosic biomass. The challenge faced with predicting outcomes of transgenic plants with reduced lignin is due in part to the presence of unique protein-protein interactions that influence the regulation and metabolic flux in the pathway. Yet, it is unclear why certain plants have evolved to create these protein complexes. In this study, we use mathematical models to investigate the role that the protein complex, formed specifically between Ptr4CL3 and Ptr4CL5 enzymes, have on the monolignol biosynthesis pathway. The role of this Ptr4CL3-Ptr4CL5 enzyme complex on the steady state flux distribution was quantified by performing Monte Carlo simulations. The effect of this complex on the robustness and the homeostatic properties of the pathway were identified by performing sensitivity and stability analyses, respectively. Results from these robustness and stability analyses suggest that the monolignol biosynthetic pathway is resilient to mild perturbations in the presence of the Ptr4CL3-Ptr4CL5 complex. Specifically, the presence of Ptr4CL3-Ptr4CL5 complex increased the stability of the pathway by 22%. The robustness in the pathway is maintained due to the presence of multiple enzyme isoforms as well as the presence of alternative pathways resulting from the presence of the Ptr4CL3-Ptr4CL5 complex.

Signalling assemblies: the odds of symmetry

The assembly of proteins into complexes is fundamental to nearly all biological signalling processes. Symmetry is a dominant feature of the structures of experimentally determined protein complexes, observed in the vast majority of homomers and many heteromers. However, some asymmetric structures exist, and asymmetry also often forms transiently, intractable to traditional structure determination methods. Here, we explore the role of protein complex symmetry and asymmetry in cellular signalling, focusing on receptors, transcription factors and transmembrane channels, among other signalling assemblies. We highlight a recurrent tendency for asymmetry to be crucial for signalling function, often being associated with activated states. We conclude with a discussion of how consideration of protein complex symmetry and asymmetry has significant potential implications and applications for pharmacology and human disease.

Functional determinants of protein assembly into homomeric complexes

Approximately half of proteins with experimentally determined structures can interact with other copies of themselves and assemble into homomeric complexes, the overwhelming majority of which (>96%) are symmetric. Although homomerisation is often assumed to a functionally beneficial result of evolutionary selection, there has been little systematic analysis of the relationship between homomer structure and function. Here, utilizing the large numbers of structures and functional annotations now available, we have investigated how proteins that assemble into different types of homomers are associated with different biological functions. We observe that homomers from different symmetry groups are significantly enriched in distinct functions, and can often provide simple physical and geometrical explanations for these associations in regards to substrate recognition or physical environment. One of the strongest associations is the tendency for metabolic enzymes to form dihedral complexes, which we suggest is closely related to allosteric regulation. We provide a physical explanation for why allostery is related to dihedral complexes: it allows for efficient propagation of conformational changes across isologous (i.e. symmetric) interfaces. Overall we demonstrate a clear relationship between protein function and homomer symmetry that has important implications for understanding protein evolution, as well as for predicting protein function and quaternary structure.

Regulation, evolution and consequences of cotranslational protein complex assembly

Most proteins assemble into complexes, which are involved in almost all cellular processes. Thus it is crucial for cell viability that mechanisms for correct assembly exist. The timing of assembly plays a key role in determining the fate of the protein: if the protein is allowed to diffuse into the crowded cellular milieu, it runs the risk of forming non-specific interactions, potentially leading to aggregation or other deleterious outcomes. It is therefore expected that strong regulatory mechanisms should exist to ensure efficient assembly. In this review we discuss the cotranslational assembly of protein complexes and discuss how it occurs, ways in which it is regulated, potential disadvantages of cotranslational interactions between proteins and the implications for the inheritance of dominant-negative genetic disorders.

Co-translational assembly of protein complexes

The interaction of biological macromolecules is a fundamental attribute of cellular life. Proteins, in particular, often form stable complexes with one another. Although the importance of protein complexes is widely recognized, we still have only a very limited understanding of the mechanisms underlying their assembly within cells. In this article, we review the available evidence for one such mechanism, namely the coupling of protein complex assembly to translation at the polysome. We discuss research showing that co-translational assembly can occur in both prokaryotic and eukaryotic organisms and can have important implications for the correct functioning of the complexes that result. Co-translational assembly can occur for both homomeric and heteromeric protein complexes and for both proteins that are translated directly into the cytoplasm and those that are translated into or across membranes. Finally, we discuss the properties of proteins that are most likely to be associated with co-translational assembly.

Operon Gene Order Is Optimized for Ordered Protein Complex Assembly

The assembly of heteromeric protein complexes is an inherently stochastic process in which multiple genes are expressed separately into proteins, which must then somehow find each other within the cell. Here, we considered one of the ways by which prokaryotic organisms have attempted to maximize the efficiency of protein complex assembly: the organization of subunit-encoding genes into operons. Using structure-based assembly predictions, we show that operon gene order has been optimized to match the order in which protein subunits assemble. Exceptions to this are almost entirely highly expressed proteins for which assembly is less stochastic and for which precisely ordered translation offers less benefit. Overall, these results show that ordered protein complex assembly pathways are of significant biological importance and represent a major evolutionary constraint on operon gene organization.

•

Operon-encoded subunits tend to be encoded by neighboring genes and form large interfaces

•

Operon gene order is often optimized for the order of protein complex assembly

•

Exceptions are mostly highly expressed proteins for which assembly is less stochastic

Principles of assembly reveal a periodic table of protein complexes

Structural insights into protein complexes have had a broad impact on our understanding of biological function and evolution. In this work, we sought a comprehensive understanding of the general principles underlying quaternary structure organization in protein complexes. We first examined the fundamental steps by which protein complexes can assemble, using experimental and structure-based characterization of assembly pathways. Most assembly transitions can be classified into three basic types, which can then be used to exhaustively enumerate a large set of possible quaternary structure topologies. These topologies, which include the vast majority of observed protein complex structures, enable a natural organization of protein complexes into a periodic table. On the basis of this table, we can accurately predict the expected frequencies of quaternary structure topologies, including those not yet observed. These results have important implications for quaternary structure prediction, modeling, and engineering.

Understanding the fabric of protein crystals: computational classification of biological interfaces and crystal contacts

Modern structural biology still draws the vast majority of information from crystallography, a technique where the objects being investigated are embedded in a crystal lattice. Given the complexity and variety of those objects, it becomes fundamental to computationally assess which of the interfaces in the lattice are biologically relevant and which are simply crystal contacts. Since the mid-1990s, several approaches have been applied to obtain high-accuracy classification of crystal contacts and biological protein–protein interfaces. This review provides an overview of the concepts and main approaches to protein interface classification: thermodynamic estimation of interface stability, evolutionary approaches based on conservation of interface residues, and co-occurrence of the interface across different crystal forms. Among the three categories, evolutionary approaches offer the strongest promise for improvement, thanks to the incessant growth in sequence knowledge. Importantly, protein interface classification algorithms can also be used on multimeric structures obtained using other high-resolution techniques or for protein assembly design or validation purposes. A key issue linked to protein interface classification is the identification of the biological assembly of a crystal structure and the analysis of its symmetry. Here, we highlight the most important concepts and problems to be overcome in assembly prediction. Over the next few years, tools and concepts of interface classification will probably become more frequently used and integrated in several areas of structural biology and structural bioinformatics. Among the main challenges for the future are better addressing of weak interfaces and the application of interface classification concepts to prediction problems like protein–protein docking.

Supplementary information: Supplementary data are available at Bioinformatics online.

Contact:guido.capitani@psi.ch