Evolution of oligomeric state through geometric coupling of protein interfaces

Shifting the substrate scope of dimeric pyranose oxidase from monosaccharide to glycoside preference through oligomeric state modification

Pyranose oxidase (POx) and C-glycoside oxidase (CGOx) are FAD-dependent oxidoreductases belonging to the glucose-methanol-choline oxidoreductase superfamily and share the same sequence space. Despite a shared structural fold, these two members possess homologous domains that enable (arm and head domain) or disable (insertion-1 domain and barrel-shaped bottom) oligomerization. POxs with a higher oligomerization state (dimeric or tetrameric) exclusively catalyze the oxidation of monosaccharides (d-glucose, d-xylose). In contrast, the monomeric state of POxs/CGOxs is observed to prefer glycosides (homoorientin, phlorizin) and has low activity with free monosaccharides. We aimed to engineer dimeric POx from Kitasatospora aureofaciens (KaPOx) to form a functional monomer, and monomeric POx/CGOx from Streptomyces canus (ScPOx) to a dimeric structure. Deletion of the head and arm domains of the KaPOx subunit resulted in enzyme variants with a less hydrophobic surface, thus affecting its oligomerization. These monomeric KaPOx variants KaPOx_xal and KaPOx_xalh resembled monomeric wild-type POxs/CGOxs and preferred glycosides as substrates over monosaccharides with catalytic efficiencies for phlorizin being 24 × 106 higher compared to those for d-xylose. The wild-type dimeric KaPOx showed no activity towards glycosides. We hypothesize that KaPOx_xalh is unable to react with monosaccharides because the introduced mutations alter the positions of monosaccharide-binding residues. The inability of KaPOx to react with glycosides is likely caused by steric hindrance and the inaccessibility of the active site to bulky glycosides due to dimerization. The attempt to engineer ScPOx into a dimeric structure failed at the stage of soluble expression, likely due to exposed hydrophobic patches and aggregation.

The Dynamical Asymmetry in SARS-CoV2 Protease Reveals the Exchange Between Catalytic Activity and Stability in Homodimers

The molecular approach to understanding the mechanisms of emerging diseases, like COVID-19, has largely accelerated the search for successful therapeutical strategies. In this work, we present an extensive molecular dynamics (MD) analysis of two forms of the SARS-CoV-2 main protease MPro. We analyzed the free form (apo) and compared the results with those coming from the (holo) form bound to the inhibitor Boceprevir, an FDA-approved drug repurposed for COVID-19 therapy. We applied Dynamic Cross Correlation (DCC) analysis to the MD simulations to trace the concerted motion patterns within the protein structure. Although symmetric, the homodimer in the bound form showed clearly asymmetric dynamical behavior. In particular, the presence of concerted motions was detected in the protomer where the expulsion of the substrate from the active site happened. Such behavior was not observed in the same time lapses in the apo form. These results highlight a sort of 'symmetry breaking', making a symmetric structure to display functional induced asymmetric behavior in response to a perturbation. This highly coordinated dynamics in response to an external cue confirms the character of 'complex molecular machines' of biopolymers.

The mineralocorticoid receptor forms higher order oligomers upon DNA binding

The prevailing model of steroid hormone nuclear receptor function assumes ligand-induced homodimer formation followed by binding to DNA hormone response elements (HREs). This model has been challenged by evidence showing that the glucocorticoid receptor (GR) forms tetramers upon ligand and DNA binding, which then drive receptor-mediated gene transactivation and transrepression. GR and the closely-related mineralocorticoid receptors (MR) interact to transduce corticosteroid hormone signaling, but whether they share the same quaternary arrangement is unknown. Here, we used a fluorescence imaging technique, Number & Brightness, to study oligomerization in a cell system allowing real-time analysis of receptor-DNA interactions. Agonist-bound MR forms tetramers in the nucleoplasm and higher order oligomers upon binding to HREs. Antagonists form intermediate-size quaternary arrangements, suggesting that large oligomers are essential for function. Divergence between MR and GR quaternary structure is driven by different functionality of known and new multimerization interfaces, which does not preclude formation of heteromers. Thus, influencing oligomerization may be important to selectively modulate corticosteroid signaling.

Experimental modification in thermal stability of oligomers by alanine substitution and site saturation mutagenesis of interfacial residues

For certain industrial applications, the stability of protein oligomers is important. In this study, we demonstrated an efficient method to improve the thermal stability of oligomers using the trimeric protein chloramphenicol acetyltransferase (CAT) as the model. We substituted all interfacial residues of CAT with alanine to detect residues critical for oligomer stability. Mutation of six of the forty-nine interfacial residues enhanced oligomer thermal stability. Site saturation mutagenesis was performed on these six residues to optimize the side chains. About 15% of mutations enhanced thermal stability by more than 0.5 °C and most did not disrupt activity of CAT. Certain combinations of mutations further improved thermal stability and resistance against heat treatment. The quadruple mutant, H17V/N34S/F134A/D157C, retained the same activity as the wild-type after heat treatment at 9 °C higher temperature than the wild-type CAT. Furthermore, combinations with only alanine substitutions also improved thermal stability, suggesting the method we developed can be used for rapid modification of industrially important proteins.

Diversity of function and higher-order structure within HWE sensor histidine kinases

Integral to the protein structure/function paradigm, oligomeric state is typically conserved along with function across evolution. However, notable exceptions such as the hemoglobins show how evolution can alter oligomerization to enable new regulatory mechanisms. Here, we examine this linkage in histidine kinases (HKs), a large class of widely distributed prokaryotic environmental sensors. While the majority of HKs are transmembrane homodimers, members of the HWE/HisKA2 family can deviate from this architecture as exemplified by our finding of a monomeric soluble HWE/HisKA2 HK (EL346, a photosensing light-oxygen-voltage [LOV]-HK). To further explore the diversity of oligomerization states and regulation within this family, we biophysically and biochemically characterized multiple EL346 homologs and found a range of HK oligomeric states and functions. Three LOV-HK homologs are primarily dimeric with differing structural and functional responses to light, while two Per-ARNT-Sim-HKs interconvert between differentially active monomers and dimers, suggesting dimerization might control enzymatic activity for these proteins. Finally, we examined putative interfaces in a dimeric LOV-HK, finding that multiple regions contribute to dimerization. Our findings suggest the potential for novel regulatory modes and oligomeric states beyond those traditionally defined for this important family of environmental sensors.

The mineralocorticoid receptor forms higher order oligomers upon DNA binding

The prevailing model of steroid hormone nuclear receptor function assumes ligand-induced homodimer formation followed by binding to DNA hormone response elements (HREs). This model has been challenged by evidence showing that the glucocorticoid receptor (GR) forms tetramers upon ligand and DNA binding, which then drive receptor-mediated gene transactivation and transrepression. GR and the closely-related mineralocorticoid receptors (MR) interact to transduce corticosteroid hormone signaling, but whether they share the same quaternary arrangement is unknown. Here, we used a fluorescence imaging technique, Number & Brightness, to study oligomerization in a cell system allowing real-time analysis of receptor-DNA interactions. Agonist-bound MR forms tetramers in the nucleoplasm and higher order oligomers upon binding to HREs. Antagonists form intermediate quaternary arrangements, suggesting that large oligomers are essential for function. Divergence between MR and GR quaternary structure is driven by different functionality of known and new multimerization interfaces, which does not preclude formation of heteromers. Thus, influencing oligomerization may be important to selectively modulate corticosteroid signaling.

Emerging maps of allosteric regulation in cellular networks

Allosteric regulation is classically defined as action at a distance, where a perturbation outside of a protein active site affects function. While this definition has motivated many studies of allosteric mechanisms at the level of protein structure, translating these insights to the allosteric regulation of entire cellular processes - and their crosstalk - has received less attention, despite the broad importance of allostery for cellular regulation foreseen by Jacob and Monod. Here, we revisit an evolutionary model for the widespread emergence of allosteric regulation in colocalized proteins, describe supporting evidence, and discuss emerging advances in mapping allostery in cellular networks that link precise and often allosteric perturbations at the molecular level to functional changes at the pathway and systems levels.

Structures of distantly related interacting protein homologs are less divergent than non-interacting homologs

Homologous proteins can display high structural variation due to evolutionary divergence at low sequence identity. This classical inverse relationship between sequence identity and structural similarity, established many years ago, has remained true between homologous proteins of known structure over time. However, a large number of heteromeric proteins also exist in the structural data bank, where the interacting subunits belong to the same fold and maintain low sequence identity between themselves. It is not clear if there is any selection pressure to deviate from the inverse sequence-structure relationship for such interacting distant homologs, in comparison to pairs of homologs which are not known to interact. We examined 12,824 fold pairs of interacting homologs of known structure, which includes both heteromers and multi-domain proteins. These were compared with monomeric proteins, resulting in 26,082 fold pairs as a dataset of non-interacting homologous systems. Interacting homologs were found to retain higher structural similarity than non-interacting homologs at diminishing sequence identity in a statistically significant manner. Interacting homologs are more similar in their 3D structures than non-interacting homologs and have a preference towards symmetric association. There appears to be a structural constraint between remote homologs due to this commitment.

Structural plasticity enables evolution and innovation of RuBisCO assemblies

Oligomerization is a core structural feature that defines the form and function of many proteins. Most proteins form molecular complexes; however, there remains a dearth of diversity-driven structural studies investigating the evolutionary trajectory of these assemblies. Ribulose-1,5-bisphosphate carboxylase-oxygenase (RuBisCO) is one such enzyme that adopts multiple assemblies, although the origins and distribution of its different oligomeric states remain cryptic. Here, we retrace the evolution of ancestral and extant form II RuBisCOs, revealing a complex and diverse history of oligomerization. We structurally characterize a newly discovered tetrameric RuBisCO, elucidating how solvent-exposed surfaces can readily adopt new interactions to interconvert or give rise to new oligomeric states. We further use these principles to engineer and demonstrate how changes in oligomerization can be mediated by relatively few mutations. Our findings yield insight into how structural plasticity may give rise to new oligomeric states.

Oligomeric interactions maintain active-site structure in a noncooperative enzyme family

The evolutionary benefit accounting for widespread conservation of oligomeric structures in proteins lacking evidence of intersubunit cooperativity remains unclear. Here, crystal and cryo-EM structures, and enzymological data, demonstrate that a conserved tetramer interface maintains the active-site structure in one such class of proteins, the short-chain dehydrogenase/reductase (SDR) superfamily. Phylogenetic comparisons support a significantly longer polypeptide being required to maintain an equivalent active-site structure in the context of a single subunit. Oligomerization therefore enhances evolutionary fitness by reducing the metabolic cost of enzyme biosynthesis. The large surface area of the structure-stabilizing oligomeric interface yields a synergistic gain in fitness by increasing tolerance to activity-enhancing yet destabilizing mutations. We demonstrate that two paralogous SDR superfamily enzymes with different specificities can form mixed heterotetramers that combine their individual enzymological properties. This suggests that oligomerization can also diversify the functions generated by a given metabolic investment, enhancing the fitness advantage provided by this architectural strategy.

Exploring the therapeutic potential of modern and ancestral phenylalanine/tyrosine ammonia-lyases as supplementary treatment of hereditary tyrosinemia

Phenylalanine/tyrosine ammonia-lyases (PAL/TALs) have been approved by the FDA for treatment of phenylketonuria and may harbour potential for complementary treatment of hereditary tyrosinemia Type I. Herein, we explore ancestral sequence reconstruction as an enzyme engineering tool to enhance the therapeutic potential of PAL/TALs. We reconstructed putative ancestors from fungi and compared their catalytic activity and stability to two modern fungal PAL/TALs. Surprisingly, most putative ancestors could be expressed as functional tetramers in Escherichia coli and thus retained their ability to oligomerize. All ancestral enzymes displayed increased thermostability compared to both modern enzymes, however, the increase in thermostability was accompanied by a loss in catalytic turnover. One reconstructed ancestral enzyme in particular could be interesting for further drug development, as its ratio of specific activities is more favourable towards tyrosine and it is more thermostable than both modern enzymes. Moreover, long-term stability assessment showed that this variant retained substantially more activity after prolonged incubation at 25 °C and 37 °C, as well as an increased resistance to incubation at 60 °C. Both of these factors are indicative of an extended shelf-life of biopharmaceuticals. We believe that ancestral sequence reconstruction has potential for enhancing the properties of enzyme therapeutics, especially with respect to stability. This work further illustrates that resurrection of putative ancestral oligomeric proteins is feasible and provides insight into the extent of conservation of a functional oligomerization surface area from ancestor to modern enzyme.

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.

Infinite Assembly of Folded Proteins in Evolution, Disease, and Engineering

Mutations and changes in a protein's environment are well known for their potential to induce misfolding and aggregation, including amyloid formation. Alternatively, such perturbations can trigger new interactions that lead to the polymerization of folded proteins. In contrast to aggregation, this process does not require misfolding and, to highlight this difference, we refer to it as agglomeration. This term encompasses the amorphous assembly of folded proteins as well as the polymerization in one, two, or three dimensions. We stress the remarkable potential of symmetric homo-oligomers to agglomerate even by single surface point mutations, and we review the double-edged nature of this potential: how aberrant assemblies resulting from agglomeration can lead to disease, but also how agglomeration can serve in cellular adaptation and be exploited for the rational design of novel biomaterials.

Structural and dynamical description of the enzymatic reaction of a phosphohexomutase

Enzymes are known to adopt various conformations at different points along their catalytic cycles. Here, we present a comprehensive analysis of 15 isomorphous, high resolution crystal structures of the enzyme phosphoglucomutase from the bacterium Xanthomonas citri. The protein was captured in distinct states critical to function, including enzyme-substrate, enzyme-product, and enzyme-intermediate complexes. Key residues in ligand recognition and regions undergoing conformational change are identified and correlated with the various steps of the catalytic reaction. In addition, we use principal component analysis to examine various subsets of these structures with two goals: (1) identifying sites of conformational heterogeneity through a comparison of room temperature and cryogenic structures of the apo-enzyme and (2) a priori clustering of the enzyme-ligand complexes into functionally related groups, showing sensitivity of this method to structural features difficult to detect by traditional methods. This study captures, in a single system, the structural basis of diverse substrate recognition, the subtle impact of covalent modification, and the role of ligand-induced conformational change in this representative enzyme of the α-D-phosphohexomutase superfamily.

Inferring and Using Protein Quaternary Structure Information from Crystallographic Data

A precise knowledge of the quaternary structure of proteins is essential to illuminate both their function and their evolution. The major part of our knowledge on quaternary structure is inferred from X-ray crystallography data, but this inference process is hard and error-prone. The difficulty lies in discriminating fortuitous protein contacts, which make up the lattice of protein crystals, from biological protein contacts that exist in the native cellular environment. Here, we review methods devised to discriminate between both types of contacts and describe resources for downloading protein quaternary structure information and identifying high-confidence quaternary structures. The use of high-confidence datasets of quaternary structures will be critical for the analysis of structural, functional, and evolutionary properties of proteins.

Modeling protein quaternary structure of homo- and hetero-oligomers beyond binary interactions by homology

Cellular processes often depend on interactions between proteins and the formation of macromolecular complexes. The impairment of such interactions can lead to deregulation of pathways resulting in disease states, and it is hence crucial to gain insights into the nature of macromolecular assemblies. Detailed structural knowledge about complexes and protein-protein interactions is growing, but experimentally determined three-dimensional multimeric assemblies are outnumbered by complexes supported by non-structural experimental evidence. Here, we aim to fill this gap by modeling multimeric structures by homology, only using amino acid sequences to infer the stoichiometry and the overall structure of the assembly. We ask which properties of proteins within a family can assist in the prediction of correct quaternary structure. Specifically, we introduce a description of protein-protein interface conservation as a function of evolutionary distance to reduce the noise in deep multiple sequence alignments. We also define a distance measure to structurally compare homologous multimeric protein complexes. This allows us to hierarchically cluster protein structures and quantify the diversity of alternative biological assemblies known today. We find that a combination of conservation scores, structural clustering, and classical interface descriptors, can improve the selection of homologous protein templates leading to reliable models of protein complexes.

Reconstructing Ancient Proteins to Understand the Causes of Structure and Function

A central goal in biochemistry is to explain the causes of protein sequence, structure, and function. Mainstream approaches seek to rationalize sequence and structure in terms of their effects on function and to identify function's underlying determinants by comparing related proteins to each other. Although productive, both strategies suffer from intrinsic limitations that have left important aspects of many proteins unexplained. These limits can be overcome by reconstructing ancient proteins, experimentally characterizing their properties, and retracing their evolution through time. This approach has proven to be a powerful means for discovering how historical changes in sequence produced the functions, structures, and other physical/chemical characteristics of modern proteins. It has also illuminated whether protein features evolved because of functional optimization, historical constraint, or blind chance. Here we review recent studies employing ancestral protein reconstruction and show how they have produced new knowledge not only of molecular evolutionary processes but also of the underlying determinants of modern proteins' physical, chemical, and biological properties.

Physical and molecular bases of protein thermal stability and cold adaptation

The molecular bases of thermal and cold stability and adaptation, which allow proteins to remain folded and functional in the temperature ranges in which their host organisms live and grow, are still only partially elucidated. Indeed, both experimental and computational studies fail to yield a fully precise and global physical picture, essentially because all effects are context-dependent and thus quite intricate to unravel. We present a snapshot of the current state of knowledge of this highly complex and challenging issue, whose resolution would enable large-scale rational protein design.

Coevolution analysis of Hepatitis C virus genome to identify the structural and functional dependency network of viral proteins

A novel computational approach of coevolution analysis allowed us to reconstruct the protein-protein interaction network of the Hepatitis C Virus (HCV) at the residue resolution. For the first time, coevolution analysis of an entire viral genome was realized, based on a limited set of protein sequences with high sequence identity within genotypes. The identified coevolving residues constitute highly relevant predictions of protein-protein interactions for further experimental identification of HCV protein complexes. The method can be used to analyse other viral genomes and to predict the associated protein interaction networks.

Comparative analyses of quaternary arrangements in homo-oligomeric proteins in superfamilies: Functional implications

A comprehensive analysis of the quaternary features of distantly related homo-oligomeric proteins is the focus of the current study. This study has been performed at the levels of quaternary state, symmetry, and quaternary structure. Quaternary state and quaternary structure refers to the number of subunits and spatial arrangements of subunits, respectively. Using a large dataset of available 3D structures of biologically relevant assemblies, we show that only 53% of the distantly related homo-oligomeric proteins have the same quaternary state. Considering these homologous homo-oligomers with the same quaternary state, conservation of quaternary structures is observed only in 38% of the pairs. In 36% of the pairs of distantly related homo-oligomers with different quaternary states the larger assembly in a pair shows high structural similarity with the entire quaternary structure of the related protein with lower quaternary state and it is referred as "Russian doll effect." The differences in quaternary state and structure have been suggested to contribute to the functional diversity. Detailed investigations show that even though the gross functions of many distantly related homo-oligomers are the same, finer level differences in molecular functions are manifested by differences in quaternary states and structures. Comparison of structures of biological assemblies in distantly and closely related homo-oligomeric proteins throughout the study differentiates the effects of sequence divergence on the quaternary structures and function. Knowledge inferred from this study can provide insights for improved protein structure classification and function prediction of homo-oligomers. Proteins 2016; 84:1190-1202. © 2016 Wiley Periodicals, Inc.

Structural basis for the antipolymer activity of Hb ζ2βs2 trapped in a tense conformation

The phenotypical severity of sickle-cell disease (SCD) can be mitigated by modifying mutant hemoglobin S (Hb S, Hb α₂β^s₂) to contain embryonic ζ-globin in place of adult α-globin subunits (Hb ζ₂β^s₂). Crystallographical analyses of liganded Hb ζζ₂β^s₂, though, demonstrate a tense (T-state) quaternary structure that paradoxically predicts its participation in--rather than its exclusion from--pathological deoxyHb S polymers. We resolved this structure-function conundrum by examining the effects of α→ζ exchange on the characteristics of specific amino acids that mediate sickle polymer assembly. Superposition analyses of the β^s subunits of T-state deoxyHb α₂β^s₂ and T-state CO-liganded Hb ζ₂β^s₂ reveal significant displacements of both mutant β^sVal6 and conserved β-chain contact residues, predicting weakening of corresponding polymer-stabilizing interactions. Similar comparisons of the α- and ζ-globin subunits implicate four amino acids that are either repositioned or undergo non-conservative substitution, abrogating critical polymer contacts. CO-Hb ζ₂β^s₂ additionally exhibits a unique trimer-of-heterotetramers crystal packing that is sustained by novel intermolecular interactions involving the pathological β^sVal6, contrasting sharply with the classical double-stranded packing of deoxyHb S. Finally, the unusually large buried solvent-accessible surface area for CO-Hb ζ₂β^s₂ suggests that it does not co-assemble with deoxyHb S in vivo. In sum, the antipolymer activities of Hb ζ₂β^s₂ appear to arise from both repositioning and replacement of specific α- and β^s-chain residues, favoring an alternate T-state solution structure that is excluded from pathological deoxyHb S polymers. These data account for the antipolymer activity of Hb ζ₂β^s₂, and recommend the utility of SCD therapeutics that capitalize on α-globin exchange strategies.

The value of protein structure classification information-Surveying the scientific literature

The Structural Classification of Proteins (SCOP) and Class, Architecture, Topology, Homology (CATH) databases have been valuable resources for protein structure classification for over 20 years. Development of SCOP (version 1) concluded in June 2009 with SCOP 1.75. The SCOPe (SCOP-extended) database offers continued development of the classic SCOP hierarchy, adding over 33,000 structures. We have attempted to assess the impact of these two decade old resources and guide future development. To this end, we surveyed recent articles to learn how structure classification data are used. Of 571 articles published in 2012-2013 that cite SCOP, 439 actually use data from the resource. We found that the type of use was fairly evenly distributed among four top categories: A) study protein structure or evolution (27% of articles), B) train and/or benchmark algorithms (28% of articles), C) augment non-SCOP datasets with SCOP classification (21% of articles), and D) examine the classification of one protein/a small set of proteins (22% of articles). Most articles described computational research, although 11% described purely experimental research, and a further 9% included both. We examined how CATH and SCOP were used in 158 articles that cited both databases: while some studies used only one dataset, the majority used data from both resources. Protein structure classification remains highly relevant for a diverse range of problems and settings.

Protein-Protein Interfaces in Viral Capsids Are Structurally Unique

Viral capsids exhibit elaborate and symmetrical architectures of defined sizes and remarkable mechanical properties not seen with cellular macromolecular complexes. Given the uniqueness of the higher-order organization of viral capsid proteins in the virosphere, we explored the question of whether the patterns of protein-protein interactions within viral capsids are distinct from those in generic protein complexes. Our comparative analysis involving a non-redundant set of 551 inter-subunit interfaces in viral capsids from VIPERdb and 20,014 protein-protein interfaces in non-capsid protein complexes from the Protein Data Bank found 418 generic protein-protein interfaces that share similar physicochemical patterns with some protein-protein interfaces in the capsid set, using the program PCalign we developed for comparing protein-protein interfaces. This overlap in the structural space of protein-protein interfaces is significantly small, with a p-value <0.0001, based on a permutation test on the total set of protein-protein interfaces. Furthermore, the generic protein-protein interfaces that bear similarity in their spatial and chemical arrangement with capsid ones are mostly small in size with fewer than 20 interfacial residues, which results from the relatively limited choices of natural design for small interfaces rather than having significant biological implications in terms of functional relationships. We conclude based on this study that protein-protein interfaces in viral capsids are non-representative of patterns in the smaller, more compact cellular protein complexes. Our finding highlights the design principle of building large biological containers from repeated, self-assembling units and provides insights into specific targets for antiviral drug design for improved efficacy.

The Energetics of Streptococcal Enolase Octamer Formation: The Quantitative Contributions of the Last Eight Amino Acids at the Carboxy-Terminus

The enolase produced by Streptococcus pyogenes is a homo-octamer whose overall shape resembles that of a donut. The octamer is best described as a tetramer of dimers. As such, it contains two types of interfaces. The first is common to almost all enolases as most enolases that have been studied are dimers. The second is unique to the octamers and includes residues near the carboxy-terminus. The primary sequence of the enolase contains 435 residues with an added 19 as an N-terminal hexahistine tag. We have systematically truncated the carboxy-terminus, individually removing the first 8 residues. This gave rise to a series of eight structures containing respectively, 435, 434, 433, 432, 431, 430, 429 and 427 residues. The truncations cause the protein to gradually dissociate from octamers to enzymatically inactive monomers with very small amounts of intermediate tetramers and dimers. We have evaluated the contributions of the missing residues to the monomer/octamer equilibrium using a combination of analytical ultracentrifugation and activity assays. For the dissociation reaction, octamer ⇐⇒ 8 monomer truncation of all eight C-terminal residues resulted in a diminution in the standard Gibbs energy of dissociation of about 59 kJ/mole of octamer relative to the full length protein. Considering that this change is spread over eight subunits, this translates to a change in standard Gibbs interaction energy of less than 8 kJ/mole of monomer distributed over the eight monomers. The resulting proteins, containing 434, 433, 432, 431, 430, 429 and 427 residues per monomer, showed intermediate free energies of dissociation. Finally, three other mutations were introduced into our reference protein to establish how they influenced the equilibrium. The main importance of this work is it shows that for homo-multimeric proteins a small change in the standard Gibbs interaction energy between subunits can have major physiological effects.

Archaeal actin from a hyperthermophile forms a single-stranded filament

The prokaryotic origins of the actin cytoskeleton have been firmly established, but it has become clear that the bacterial actins form a wide variety of different filaments, different both from each other and from eukaryotic F-actin. We have used electron cryomicroscopy (cryo-EM) to examine the filaments formed by the protein crenactin (a crenarchaeal actin) from Pyrobaculum calidifontis, an organism that grows optimally at 90 °C. Although this protein only has ∼ 20% sequence identity with eukaryotic actin, phylogenetic analyses have placed it much closer to eukaryotic actin than any of the bacterial homologs. It has been assumed that the crenactin filament is double-stranded, like F-actin, in part because it would be hard to imagine how a single-stranded filament would be stable at such high temperatures. We show that not only is the crenactin filament single-stranded, but that it is remarkably similar to each of the two strands in F-actin. A large insertion in the crenactin sequence would prevent the formation of an F-actin-like double-stranded filament. Further, analysis of two existing crystal structures reveals six different subunit-subunit interfaces that are filament-like, but each is different from the others in terms of significant rotations. This variability in the subunit-subunit interface, seen at atomic resolution in crystals, can explain the large variability in the crenactin filaments observed by cryo-EM and helps to explain the variability in twist that has been observed for eukaryotic actin filaments.

Tracing the evolution of the p53 tetramerization domain

The tetrameric transcription factors p53, p63, and p73 evolved from a common ancestor and play key roles in tumor suppression and development. Surprisingly, p63 and p73 require a second helix in their tetramerization domain for the formation of stable tetramers that is absent in human p53, raising questions about the evolutionary processes leading to diversification. Here we determined the crystal structure of the zebrafish p53 tetramerization domain, which contains a second helix, reminiscent of p63 and p73, combined with p53-like features. Through comprehensive phylogenetic analyses, we systematically traced the evolution of vertebrate p53 family oligomerization domains back to the beginning of multicellular life. We provide evidence that their last common ancestor also had an extended p63/p73-like domain and pinpoint evolutionary events that shaped this domain during vertebrate radiation. Domain compaction and transformation of a structured into a flexible, intrinsically disordered region may have contributed to the expansion of the human p53 interactome.

Biophysics of protein evolution and evolutionary protein biophysics

The study of molecular evolution at the level of protein-coding genes often entails comparing large datasets of sequences to infer their evolutionary relationships. Despite the importance of a protein's structure and conformational dynamics to its function and thus its fitness, common phylogenetic methods embody minimal biophysical knowledge of proteins. To underscore the biophysical constraints on natural selection, we survey effects of protein mutations, highlighting the physical basis for marginal stability of natural globular proteins and how requirement for kinetic stability and avoidance of misfolding and misinteractions might have affected protein evolution. The biophysical underpinnings of these effects have been addressed by models with an explicit coarse-grained spatial representation of the polypeptide chain. Sequence-structure mappings based on such models are powerful conceptual tools that rationalize mutational robustness, evolvability, epistasis, promiscuous function performed by 'hidden' conformational states, resolution of adaptive conflicts and conformational switches in the evolution from one protein fold to another. Recently, protein biophysics has been applied to derive more accurate evolutionary accounts of sequence data. Methods have also been developed to exploit sequence-based evolutionary information to predict biophysical behaviours of proteins. The success of these approaches demonstrates a deep synergy between the fields of protein biophysics and protein evolution.

Quantitative and Systems-Based Approaches for Deciphering Bacterial Membrane Interactome and Gene Function

High-throughput genomic and proteomic methods provide a concise description of the molecular constituents of a cell, whereas systems biology strives to understand the way these components function as a whole. Recent developments, such as genome editing technologies and protein epitope-tagging coupled with high-sensitivity mass-spectrometry, allow systemic studies to be performed at an unprecedented scale. Available methods can be successfully applied to various goals, both expanding fundamental knowledge and solving applied problems. In this review, we discuss the present state and future of bacterial cell envelope interactomics, with a specific focus on host-pathogen interactions and drug target discovery. Both experimental and computational methods will be outlined together with examples of their practical implementation.

Evolution of oligomeric state through allosteric pathways that mimic ligand binding

Evolution and design of protein complexes are almost always viewed through the lens of amino acid mutations at protein interfaces. We showed previously that residues not involved in the physical interaction between proteins make important contributions to oligomerization by acting indirectly or allosterically. In this work, we sought to investigate the mechanism by which allosteric mutations act, using the example of the PyrR family of pyrimidine operon attenuators. In this family, a perfectly sequence-conserved helix that forms a tetrameric interface is exposed as solvent-accessible surface in dimeric orthologs. This means that mutations must be acting from a distance to destabilize the interface. We identified 11 key mutations controlling oligomeric state, all distant from the interfaces and outside ligand-binding pockets. Finally, we show that the key mutations introduce conformational changes equivalent to the conformational shift between the free versus nucleotide-bound conformations of the proteins.

Structure, dynamics, assembly, and evolution of protein complexes

The assembly of individual proteins into functional complexes is fundamental to nearly all biological processes. In recent decades, many thousands of homomeric and heteromeric protein complex structures have been determined, greatly improving our understanding of the fundamental principles that control symmetric and asymmetric quaternary structure organization. Furthermore, our conception of protein complexes has moved beyond static representations to include dynamic aspects of quaternary structure, including conformational changes upon binding, multistep ordered assembly pathways, and structural fluctuations occurring within fully assembled complexes. Finally, major advances have been made in our understanding of protein complex evolution, both in reconstructing evolutionary histories of specific complexes and in elucidating general mechanisms that explain how quaternary structure tends to evolve. The evolution of quaternary structure occurs via changes in self-assembly state or through the gain or loss of protein subunits, and these processes can be driven by both adaptive and nonadaptive influences.

A new model for pore formation by cholesterol-dependent cytolysins

Cholesterol Dependent Cytolysins (CDCs) are important bacterial virulence factors that form large (200–300 Å) membrane embedded pores in target cells. Currently, insights from X-ray crystallography, biophysical and single particle cryo-Electron Microscopy (cryo-EM) experiments suggest that soluble monomers first interact with the membrane surface via a C-terminal Immunoglobulin-like domain (Ig; Domain 4). Membrane bound oligomers then assemble into a prepore oligomeric form, following which the prepore assembly collapses towards the membrane surface, with concomitant release and insertion of the membrane spanning subunits. During this rearrangement it is proposed that Domain 2, a region comprising three β-strands that links the pore forming region (Domains 1 and 3) and the Ig domain, must undergo a significant yet currently undetermined, conformational change. Here we address this problem through a systematic molecular modeling and structural bioinformatics approach. Our work shows that simple rigid body rotations may account for the observed collapse of the prepore towards the membrane surface. Support for this idea comes from analysis of published cryo-EM maps of the pneumolysin pore, available crystal structures and molecular dynamics simulations. The latter data in particular reveal that Domains 1, 2 and 4 are able to undergo significant rotational movements with respect to each other. Together, our data provide new and testable insights into the mechanism of pore formation by CDCs.

Pore formation is central to the ability of cholesterol dependent cytolysins (CDCs) to act as important bacterial virulence factors. Secreted by numerous pathogens the toxins assemble into a circular ring and then perforate the target membrane to form the largest self-assembling proteinaceous pores known. In this paper we investigated computationally the conformational properties of the CDC molecule and deduced a new structural model of pore formation and membrane insertion that reconciles all experimental data. The mechanism of membrane perforation by CDCs put forward here reveals concerted and unsuspected domains motion of large amplitude, which conflicts with the currently proposed model. The work presented here procures a plausible structural mechanism of CDC oligomeric transition and furthers our understanding of pore formation by these important toxins.

An overview of recent advances in structural bioinformatics of protein-protein interactions and a guide to their principles

Rich data bearing on the structural and evolutionary principles of protein-protein interactions are paving the way to a better understanding of the regulation of function in the cell. This is particularly the case when these interactions are considered in the framework of key pathways. Knowledge of the interactions may provide insights into the mechanisms of crucial 'driver' mutations in oncogenesis. They also provide the foundation toward the design of protein-protein interfaces and inhibitors that can abrogate their formation or enhance them. The main features to learn from known 3-D structures of protein-protein complexes and the extensive literature which analyzes them computationally and experimentally include the interaction details which permit undertaking structure-based drug discovery, the evolution of complexes and their interactions, the consequences of alterations such as post-translational modifications, ligand binding, disease causing mutations, host pathogen interactions, oligomerization, aggregation and the roles of disorder, dynamics, allostery and more to the protein and the cell. This review highlights some of the recent advances in these areas, including design, inhibition and prediction of protein-protein complexes. The field is broad, and much work has been carried out in these areas, making it challenging to cover it in its entirety. Much of this is due to the fast increase in the number of molecules whose structures have been determined experimentally and the vast increase in computational power. Here we provide a concise overview.

Protein flexibility facilitates quaternary structure assembly and evolution

The flexibility of individual proteins aids their evolutionary recruitment into complexes with increasing numbers of distinct subunits.

The intrinsic flexibility of proteins allows them to undergo large conformational fluctuations in solution or upon interaction with other molecules. Proteins also commonly assemble into complexes with diverse quaternary structure arrangements. Here we investigate how the flexibility of individual protein chains influences the assembly and evolution of protein complexes. We find that flexibility appears to be particularly conducive to the formation of heterologous (i.e., asymmetric) intersubunit interfaces. This leads to a strong association between subunit flexibility and homomeric complexes with cyclic and asymmetric quaternary structure topologies. Similarly, we also observe that the more nonhomologous subunits that assemble together within a complex, the more flexible those subunits tend to be. Importantly, these findings suggest that subunit flexibility should be closely related to the evolutionary history of a complex. We confirm this by showing that evolutionarily more recent subunits are generally more flexible than evolutionarily older subunits. Finally, we investigate the very different explorations of quaternary structure space that have occurred in different evolutionary lineages. In particular, the increased flexibility of eukaryotic proteins appears to enable the assembly of heteromeric complexes with more unique components.

Proteins often interact with other proteins and assemble into complexes. Here we show that the flexibility of individual proteins is important for their recruitment to complexes, as it facilitates the formation of asymmetric interfaces between different subunits. The role of flexibility becomes increasingly important as a greater number of distinct proteins are packed together within a single complex: the more distinct subunits, the more flexible those subunits need to be. A consequence of this is that, when a protein complex gains a new subunit during evolution, the newer subunit will tend to be more flexible than the older subunits. This suggests that we may be able to partially reconstruct the evolutionary history of a protein complex by considering the flexibility of its subunits. We also find that the types of protein complexes an organism forms are closely related to the flexibility of its proteins, with eukaryotic species, and particularly animals, using their increased flexibility to assemble complexes involving more distinct components.

The active form of E6-associated protein (E6AP)/UBE3A ubiquitin ligase is an oligomer

Employing 125I-polyubiquitin chain formation as a functional readout of ligase activity, biochemical and biophysical evidence demonstrates that catalytically active E6-associated protein (E6AP)/UBE3A is an oligomer. Based on an extant structure previously discounted as an artifact of crystal packing forces, we propose that the fully active form of E6AP is a trimer, analysis of which reveals a buried surface of 7508Å2 and radially symmetric interacting residues that are conserved within the Hect (homologous to E6AP C terminus) ligase superfamily. An absolutely conserved interaction between Phe(727) and a hydrophobic pocket present on the adjacent subunit is critical for trimer stabilization because mutation disrupts the oligomer and decreases kcat 62-fold but fails to affect E2 ubiquitin binding or subsequent formation of the Hect domain Cys(820) ubiquitin thioester catalytic intermediate. Exogenous N-acetylphenylalanylamide reversibly antagonizes Phe(727)-dependent trimer formation and catalytic activity (Ki12 mM), as does a conserved-helical peptide corresponding to residues 474–490 of E6A Pisoform 1 (Ki22M) reported to bind the hydrophobic pocket of other Hect ligases, presumably blocking Phe(727) intercalation and trimer formation. Conversely, oncogenic human papillomavirus-16/18 E6 protein significantly enhances E6AP catalytic activity by promoting trimer formation (Kactivation 1.5 nM) through the ability of E6 to form homodimers. Recombinant E6 protein additionally rescues the kcat defect of the Phe(727) mutation and that of a specific loss-of-function Angelman syndrome mutation that promotes trimer destabilization. The present findings codify otherwise disparate observations regarding the mechanism of E6AP and related Hect ligases in addition to suggesting therapeutic approaches for modulating ligase activity.

Thermal and chemical stability of two homologous POZ/BTB domains of KCTD proteins characterized by a different oligomeric organization

POZ/BTB domains are widespread modules detected in a variety of different biological contexts. Here, we report a biophysical characterization of the POZ/BTB of KCTD6, a protein that is involved in the turnover of the muscle small ankyrin-1 isoform 5 and, in combination with KCTD11, in the ubiquitination and degradation of HDAC1. The analyses show that the domain is a tetramer made up by subunits with the expected α/β structure. A detailed investigation of its stability, carried out in comparison with the homologous pentameric POZ/BTB domain isolated from KCTD5, highlights a number of interesting features, which are shared by the two domains despite their different organization. Their thermal/chemical denaturation curves are characterized by a single and sharp inflection point, suggesting that the denaturation of the two domains is a cooperative two-state process. Furthermore, both domains present a significant content of secondary structure in their denatured state and a reversible denaturation process. We suggest that the ability of these domains to fold and unfold reversibly, a property that is somewhat unexpected for these oligomeric assemblies, may have important implications for their biological function. Indeed, these properties likely favor the formation of heteromeric associations that may be essential for the intricate regulation of the processes in which these proteins are involved.

Structural templates for modeling homodimers

Oligomeric proteins are more abundant in nature than monomeric proteins, and involved in all biological processes. In the absence of an experimental structure, their subunits can be modeled from their sequence like monomeric proteins, but reliable procedures to build the oligomeric assembly are scarce. Template-based methods, which start from known protein structures, are commonly applied to model subunits. We present a method to model homodimers that relies on a structural alignment of the subunits, and test it on a set of 511 target structures recently released by the Protein Data Bank, taking as templates the earlier released structures of 3108 homodimeric proteins (H-set), and 2691 monomeric proteins that form dimer-like assemblies in crystals (M-set). The structural alignment identifies a H-set template for 97% of the targets, and in half of the cases, it yields a correct model of the dimer geometry and residue-residue contacts in the target. It also identifies a M-set template for most of the targets, and some of the crystal dimers are very similar to the target homodimers. The procedure efficiently detects homology at low levels of sequence identities, and points to erroneous quaternary structures in the Protein Data Bank. The high coverage of the target set suggests that the content of the Protein Data Bank already approaches the structural diversity of protein assemblies in nature, and that template-based methods should become the choice method for modeling oligomeric as well as monomeric proteins.

Structural, evolutionary, and assembly principles of protein oligomerization

In the protein universe, 30-50% of proteins self-assemble to form symmetrical complexes consisting of multiple copies of themselves, called homomers. The prevalence of homomers motivates us to review many of their properties. In Section 1, we describe the methods and challenges associated with quaternary structure inference-these methods are indeed at the basis of any analysis on homomers. In Section 2, we describe the morphological properties of homomers, as well as the database 3DComplex, which provides a taxonomy for both homomeric and heteromeric protein complexes. In Section 3, we review interface properties of homomeric complexes. In Section 4, we then present recent findings on the evolution of homomer interfaces, which we link in Section 5 to the evolution of homomers as entire entities. In Section 6, we discuss mechanisms involved in their assembly and how these mechanisms can be linked to evolution.

Evolutionary, physicochemical, and functional mechanisms of protein homooligomerization

Protein homooligomers afford several important benefits for the cell; they mediate and regulate gene expression, activity of many enzymes, ion channels, receptors, and cell-cell adhesion processes. The evolutionary and physical mechanisms of oligomer formation are very diverse and are not well understood. Certain homooligomeric states may be conserved within protein subfamilies and between different subfamilies, therefore providing the specificity to particular substrates while minimizing interactions with unwanted partners. In addition, transitions between different oligomeric states may regulate protein activity and support the switch between different pathways. In this chapter, we summarize the biological importance of homooligomeric assemblies, physicochemical properties of their interfaces, experimental methods for their identification, their evolution, and role in human diseases.

On the binding affinity of macromolecular interactions: daring to ask why proteins interact

Interactions between proteins are orchestrated in a precise and time-dependent manner, underlying cellular function. The binding affinity, defined as the strength of these interactions, is translated into physico-chemical terms in the dissociation constant (K_d), the latter being an experimental measure that determines whether an interaction will be formed in solution or not. Predicting binding affinity from structural models has been a matter of active research for more than 40 years because of its fundamental role in drug development. However, all available approaches are incapable of predicting the binding affinity of protein–protein complexes from coordinates alone. Here, we examine both theoretical and experimental limitations that complicate the derivation of structure–affinity relationships. Most work so far has concentrated on binary interactions. Systems of increased complexity are far from being understood. The main physico-chemical measure that relates to binding affinity is the buried surface area, but it does not hold for flexible complexes. For the latter, there must be a significant entropic contribution that will have to be approximated in the future. We foresee that any theoretical modelling of these interactions will have to follow an integrative approach considering the biology, chemistry and physics that underlie protein–protein recognition.