Structural symmetry and protein function

Functional asymmetry in processivity clamp proteins

Symmetric homo-oligomeric proteins comprising multiple copies of identical subunits are abundant in all domains of life. To fulfill their biological function, these complexes undergo conformational changes, binding events, or posttranslational modifications, leading to loss of symmetry. Processivity clamp proteins that encircle DNA and play multiple roles in DNA replication and repair are archetypical homo-oligomeric symmetric protein complexes. The symmetrical nature of processivity clamps enables simultaneous interactions with multiple protein binding partners; such interactions result in asymmetric changes that facilitate the transition between clamp loading and DNA replication and between DNA replication and repair. The ring-shaped processivity clamps are opened and loaded onto DNA by clamp-loader complexes via asymmetric intermediates with one of the intermonomer interfaces disrupted, undergo spontaneous opening events, and bind heterogeneous partners. Eukaryotic clamp proteins are subject to ubiquitylation, SUMOylation, and acetylation, affecting their biological functions. There is increasing evidence of the functional asymmetry of the processivity clamp proteins from structural, biophysical, and computational studies. Here, we review the symmetry and asymmetry of processivity clamps and their roles in regulating the various functions of the clamps.

Cryo-EM reveals mechanisms of natural RNA multivalency

Homo-oligomerization of biological macromolecules leads to functional assemblies that are critical to understanding various cellular processes. However, RNA quaternary structures have rarely been reported. Comparative genomics analysis has identified RNA families containing hundreds of sequences that adopt conserved secondary structures and likely fold into complex three-dimensional structures. In this study, we used cryo-electron microscopy (cryo-EM) to determine structures from four RNA families, including ARRPOF and OLE forming dimers and ROOL and GOLLD forming hexameric, octameric, and dodecameric nanostructures, at 2.6- to 4.6-angstrom resolutions. These homo-oligomeric assemblies reveal a plethora of structural motifs that contribute to RNA multivalency, including kissing-loop, palindromic base-pairing, A-stacking, metal ion coordination, pseudoknot, and minor-groove interactions. These results provide the molecular basis of intermolecular interactions driving RNA multivalency with potential functional relevance.

Mass Photometry

Mass photometry (MP) is a technology for the mass measurement of biological macromolecules in solution. Its mass accuracy and resolution have transformed label-free optical detection into a quantitative measurement, enabling the identification of distinct species in a mixture and the characterization of their relative abundances. Its applicability to a variety of biomolecules, including polypeptides, nucleic acids, lipids, and sugars, coupled with the ability to quantify heterogeneity, interaction energies, and kinetics, has driven the rapid and widespread adoption of MP across the life sciences community. These applications have been largely orthogonal to those traditionally associated with microscopy, such as detection, imaging, and tracking, instead focusing on the constituents of biomolecular complexes and their change with time. Here, we present an overview of the origins of MP, its current applications, and future improvements that will further expand its scope.

Plasticity and Co-Factor-Dependent Structural Changes in the RecA Nucleoprotein Filament Studied by Small-Angle X-Ray Scattering (SAXS) Measurements and Molecular Modeling

Structural analyses of protein filaments formed by self-assembly, such as actin, tubulin, or recombinase filaments, have suffered for decades from technical issues due to difficulties in crystallization, their large size, or the dynamic behavior inherent to their cellular function. The advent of cryo-electron microscopy has finally enabled us to obtain structures at different stages of the existence of these filaments. However, these structures correspond to frozen states, and the possibility of observations in solution is still lacking, especially for filaments characterized by a high plasticity, such as the RecA protein for homologous recombination. Here, we use a combination of SAXS measurements and integrative modeling to generate the solution structure of two known forms of the RecA nucleoprotein filament, previously characterized by electron microscopy and resolved by X-ray crystallography. The two forms differ in the cofactor bound to RecA-RecA interfaces, either ATP or ADP. Cooperative transition from one form to the other has been observed during single-molecule experiments by pulling on the filament but also in solution by modifying solvent conditions. We first compare the SAXS data against known structural information. While the crystal structure of the ATP form matches well with the SAXS data, we deduce from the SAXS profiles of the ADP-form values of the pitch (72.0 Å) and the number of monomers per turn (6.4) that differ with respect to the crystal structure (respectively, 82.7 Å and 6.0). We then monitor the transition between the two states driven by the addition of magnesium, and we show this transition occurs with 0.3 mM Mg 2+ ions with a high cooperativity.

Structured Random Binding: a minimal model of protein-protein interactions

We describe Structured Random Binding (SRB), a minimal model of protein-protein interactions rooted in the statistical physics of disordered systems. In this model, nonspecific binding is a generic consequence of the interaction between random proteins, exhibiting a phase transition from a high temperature state where nonspecific complexes are transient and lack well-defined interaction interfaces, to a low temperature state where the complex structure is frozen and a definite interaction interface is present. Numerically, weakly-bound nonspecific complexes can evolve into tightly-bound, highly specific complexes, but only if the structural correlation length along the peptide backbone is short; moreover, evolved tightly-bound homodimers favor the same interface structure that is predominant in real protein homodimers.

Native Capillary Nanogel Electrophoresis Assay of Inhibitors of Neuraminidases Derived from H1N1 and H5N1 Influenza A Pandemics

Tetrameric neuraminidases cleave the end-capping sialylated monomer from oligosaccharide ligands at the surface of a host cell infected by the influenza A virus. This cleavage releases the replicated virions from the host cell, making drugs that inhibit neuraminidase function effective to treat influenza A infections. A capillary electrophoresis separation-based assay is reported that maintains the native structure of tetrameric viral neuraminidases derived from H1N1 or H5N1 influenza A pandemics which convert, in-real time, a substrate that mimics 6'-sialyllated threonine-linked glycans on human cells. The assay integrates the enzyme reaction with the separation and is operated using a background electrolyte containing 100 mM NaCl with a thermally reversible nanogel in a 10 μm inner diameter fused silica capillary. In addition to defining the 0.4 nL reaction zone maintained at 37 °C, the nanogel medium resolves the substrate from contaminants as well as the substrate from the product before and after the enzymatic conversion. The enzyme activity is quantifiable based on the percent conversion observed in the presence of a range of inhibitor concentrations. For 1918 H1N1 (A/Brevig Mission/1/18) neuraminidase, the inhibition constant of the transition state analog 2,3-dehydro-2-deoxy-N-acetylneuraminic acid (DANA) is 3.5 ± 0.8 μM (n = 5). The inhibition constants for oseltamivir acid (inhibiting compound of Tamiflu) and peramivir (Rapivab) are 18.2 ± 0.5 nM (n = 3) and 67 ± 8 nM (n = 3), respectively. For 2004 H5N1 (A/Vietnam/1203/2004) neuraminidase, which contained a foreign tetramerization domain to maintain the structure, the inhibition constant for peramivir is 5.4 nM.

DNA Tetrahedron Mass-Tagged Probe Set for the Programmed Detection of Protein Trimers by Point-to-Point Recognition and Induced Self-Assembly in Living Cells

Multimeric proteins normally perform different biological functions from their monomer components. Thus, precise recognition and quantitative detection of multimeric proteins can benefit a better understanding of complex biological processes and their roles in disease diagnosis and treatment. The challenge herein is to distinguish the multimeric proteins containing identical monomer components and recognize all the monomers in a multimeric protein on spatial scales. This situation is likely to become more significant for homomultimeric proteins. In this study, a DNA polyhedron mass-tagged probe set strategy was developed for the programmed detection of multimeric proteins in living cells. The probe set comprised recognition and displacement probes, a DNA polyhedron probe, and a mass-tagged probe. After point-to-point recognition of each monomer in the target protein complex by recognition and displacement probes, the DNA polyhedron probe could integrate the information on all the protein monomers by carrying out induced self-assembly via a cascaded toehold-mediated strand-displacement (TMSD) reaction. Afterward, the mass-tagged probe collected the integrated information, and the mass tag in the probe was released by ultraviolet (UV) irradiation and detected by mass spectrometry (MS). Using the tmTNF-α homotrimer as an example, its expression levels in different breast cancer cell lines were ultimately determined using this probe set containing a DNA tetrahedron probe. This study is among the first to quantitatively detect multimeric proteins in living cells. Using a similar strategy, more DNA polyhedron mass-tagged probe sets can be developed for the detection of higher-order multimeric proteins.

Rapid and accurate prediction of protein homo-oligomer symmetry using Seq2Symm

The majority of proteins must form higher-order assemblies to perform their biological functions, yet few machine learning models can accurately and rapidly predict the symmetry of assemblies involving multiple copies of the same protein chain. Here, we address this gap by finetuning several classes of protein foundation models, to predict homo-oligomer symmetry. Our best model named Seq2Symm, which utilizes ESM2, outperforms existing template-based and deep learning methods achieving an average AUC-PR of 0.47, 0.44 and 0.49 across homo-oligomer symmetries on three held-out test sets compared to 0.24, 0.24 and 0.25 with template-based search. Seq2Symm uses a single sequence as input and can predict at the rate of ~80,000 proteins/hour. We apply this method to 5 proteomes and ~3.5 million unlabeled protein sequences, showing its promise to be used in conjunction with downstream computationally intensive all-atom structure generation methods such as RoseTTAFold2 and AlphaFold2-multimer. Code, datasets, model are available at: https://github.com/microsoft/seq2symm .

Heligeom: A web resource to generate, analyze, and visualize filament architectures based on pairwise association geometries of biological macromolecules

At the subcellular level, macromolecules self-assemble to form molecular machinery in which the assembly modes play critical roles: the structural integrity of cell walls that allows mechanical growth, the maintenance and repair of the genetic material, membrane flow control, protein chaperoning, and ATP production, to cite just a few examples. As molecular modeling expands its scope to such systems, structural biologists are faced with the difficulty of understanding the structure and dynamics of these supramolecular assemblies. We present Heligeom, a webserver that offers a simple and efficient means for analyzing and constructing oligomeric assemblies based on user-provided structures of two interacting units. The input 3D coordinates may result from structure determination, simulations, docking trials, or deep-learning tools such as AlphaFold. For a given interface, Heligeom outputs the mathematical helical parameters of the corresponding oligomeric form, including axis, pitch, handedness, number of monomers per turn, etc. The server also allows leveraging these parameters to construct oligomers of specified size, visualizing them interactively, and downloading them as PDB files. For subunits (protomers) having multiple interaction geometries, the different interfaces and their global assembly geometry can be visualized and compared. Heligeom thus allows explicitly linking protomer-protomer interfaces to the oligomeric architecture, illuminating possible sources of plasticity in protein filaments such as mutations or thermal, mechanical, or chemical perturbations. Heligeom thus constitutes an intuitive tool to accompany integrative modeling of oligomeric filamentous assemblies. Examples of its application at different structural levels are presented. Heligeom webserver can be accessed at https://heligeom.galaxy.ibpc.fr.

Physicochemical features of subunit interfaces and their role in self-assembly across the ferritin superfamily

Ferritins are ubiquitous and play a critical role in iron homeostasis. They are classified into four main subfamilies: classical, bacterial, bacterioferritin, and Dps. These are characterized by subunits with a four-helical bundle domain and interact through three distinct regions-one antiparallel interface (IntA) and two perpendicular interfaces (IntB and IntC), collectively forming a cage-like structure. Here, we attempt to characterize the variability of these interfaces across subfamilies. We found that IntA is essential for the dimeric unit assembly and is likely to assemble first, followed by the smaller interfaces of IntB and IntC (in any order), which are crucial for cage formation. These interfaces are unique in that they are less packed, although chemically stable, and their size lies between that of protein-protein complex and obligate homodimers. This study provides a detailed exploration of the ferritin interfaces, offering insights into their assembly and their importance as carrier proteins.

Hierarchical design of pseudosymmetric protein nanocages

Discrete protein assemblies ranging from hundreds of kilodaltons to hundreds of megadaltons in size are a ubiquitous feature of biological systems and perform highly specialized functions1,2. Despite remarkable recent progress in accurately designing new self-assembling proteins, the size and complexity of these assemblies has been limited by a reliance on strict symmetry3. Here, inspired by the pseudosymmetry observed in bacterial microcompartments and viral capsids, we developed a hierarchical computational method for designing large pseudosymmetric self-assembling protein nanomaterials. We computationally designed pseudosymmetric heterooligomeric components and used them to create discrete, cage-like protein assemblies with icosahedral symmetry containing 240, 540 and 960 subunits. At 49, 71 and 96 nm diameter, these nanocages are the largest bounded computationally designed protein assemblies generated to date. More broadly, by moving beyond strict symmetry, our work substantially broadens the variety of self-assembling protein architectures that are accessible through design.

Symmetry facilitated the evolution of heterospecificity and high-order stoichiometry in vertebrate hemoglobin

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

Structural basis of nearest-neighbor cooperativity in the ring-shaped gene regulatory protein TRAP from protein engineering and cryo-EM

The homo-dodecameric ring-shaped trp RNA binding attenuation protein (TRAP) from Alkalihalobacillus halodurans (Aha) binds up to twelve tryptophan ligands (Trp) and becomes activated to bind a specific sequence in the 5' leader region of the trp operon mRNA, thereby downregulating biosynthesis of Trp. Thermodynamic measurements of Trp binding have revealed a range of cooperative behavior for different TRAP variants, even if the averaged apparent affinities for Trp have been found to be similar. Proximity between the ligand binding sites, and the ligand-coupled disorder-to-order transition has implicated nearest-neighbor interactions in cooperativity. To establish a solid basis for describing nearest-neighbor cooperativity in TRAP, we engineered variants constructed with two subunits connected by a flexible linker (dTRAP). We mutated the binding sites of alternating protomers such that only every other site was competent for Trp binding (WT-Mut dTRAP). Ligand binding monitored by NMR, calorimetry, and native mass spectrometry revealed strong cooperativity in dTRAP containing adjacent binding-competent sites, but a severe binding defect when the wild-type sites were separated by mutated sites. Cryo-EM experiments of dTRAP in its ligand-free apo state, and both dTRAP and WT-Mut dTRAP in the presence of Trp, revealed progressive stabilization of loops that gate the Trp binding site and participate in RNA binding. These studies provide important insights into the thermodynamic and structural basis for the observed ligand binding cooperativity in TRAP. Such insights can be useful for understanding allosteric control networks and for the development of those with defined ligand sensitivity and regulatory control.

Comparative Analysis of Symmetry Parameters in the E2 Inner Core of the Pyruvate Dehydrogenase Complex

Recent advances in cryo-electron microscopy (cryo-EM) have facilitated the high-resolution structural determination of macromolecular complexes in their native states, providing valuable insights into their dynamic behaviors. However, insufficient understanding or experience with the cryo-EM image processing parameters can result in the loss of biological meaning. In this paper, we investigate the dihydrolipoyl acetyltransferase (E2) inner core complex of the pyruvate dehydrogenase complex (PDC) and reconstruct the 3D maps using five different symmetry parameters. The results demonstrate that the reconstructions yield structurally identical 3D models even at a near-atomic structure. This finding underscores a crucial message for researchers engaging in single-particle analysis (SPA) with relatively user-friendly and convenient image processing software. This approach helps reduce the risk of missing critical biological details, such as the dynamic properties of macromolecules.

Symmetry facilitated the evolution of heterospecificity and high-order stoichiometry in vertebrate hemoglobin

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

Frequent transitions in self-assembly across the evolution of a central metabolic enzyme

Many enzymes assemble into homomeric protein complexes comprising multiple copies of one protein. Because structural form is usually assumed to follow function in biochemistry, these assemblies are thought to evolve because they provide some functional advantage. In many cases, however, no specific advantage is known and, in some cases, quaternary structure varies among orthologs. This has led to the proposition that self-assembly may instead vary neutrally within protein families. The extent of such variation has been difficult to ascertain because quaternary structure has until recently been difficult to measure on large scales. Here, we employ mass photometry, phylogenetics, and structural biology to interrogate the evolution of homo-oligomeric assembly across the entire phylogeny of prokaryotic citrate synthases - an enzyme with a highly conserved function. We discover a menagerie of different assembly types that come and go over the course of evolution, including cases of parallel evolution and reversions from complex to simple assemblies. Functional experiments in vitro and in vivo indicate that evolutionary transitions between different assemblies do not strongly influence enzyme catalysis. Our work suggests that enzymes can wander relatively freely through a large space of possible assembly states and demonstrates the power of characterizing structure-function relationships across entire phylogenies.

Cryo-EM Structure of AAA + ATPase Thorase Reveals Novel Helical Filament Formation

The AAA+ ( A TPases a ssociated with a variety of cellular a ctivities) ATPase, Thorase, also known as ATAD1, plays multiple roles in synaptic plasticity, mitochondrial quality control and mTOR signaling through disassembling protein complexes like AMPAR and mTORC1 in an ATP-dependent manner. The Oligomerization of Thorase is crucial for its disassembly and remodeling functions. We show that wild-type Thorase forms long helical filaments in vitro , dependent on ATP binding but not hydrolysis. We report the Cryogenic Electron Microscopy (cryo-EM) structure of the Thorase filament at a resolution of 4 Å, revealing the dimeric arrangement of the basic repeating unit that is formed through a distinct interface compared to the hexameric MSP1/ATAD1E193Q assembly. Structure-guided mutagenesis confirms the role of critical amino acid residues required for filament formation, oligomerization and disassembly of mTORC1 protein complex. Together, our data reveals a novel filament structure of Thorase and provides critical information that elucidates the mechanism underlying Thorase filament formation and Thorase-mediated disassembly of the mTORC1 complex.

Dock2D: Synthetic Data for the Molecular Recognition Problem

Predicting the physical interaction of proteins is a cornerstone problem in computational biology. New classes of learning-based algorithms are actively being developed, and are typically trained end-to-end on protein complex structures extracted from the Protein Data Bank. These training datasets tend to be large and difficult to use for prototyping and, unlike image or natural language datasets, they are not easily interpretable by non-experts. We present Dock2D-IP and Dock2D-IF, two "toy" datasets that can be used to select algorithms predicting protein-protein interactions-or any other type of molecular interactions. Using two-dimensional shapes as input, each example from Dock2D-IP ("interaction pose") describes the interaction pose of two shapes known to interact and each example from Dock2D-IF ("interaction fact") describes whether two shapes form a stable complex or not, regardless of how they bind. We propose a number of baseline solutions to the problem and show that the same underlying energy function can be learned either by solving the interaction pose task (formulated as an energy-minimization "docking" problem) or the fact-of-interaction task (formulated as a binding free energy estimation problem).

Symmetry-adapted Markov state models of closing, opening, and desensitizing in α 7 nicotinic acetylcholine receptors

α7 nicotinic acetylcholine receptors (nAChRs) are homopentameric ligand-gated ion channels with critical roles in the nervous system. Recent studies have resolved and functionally annotated closed, open, and desensitized states of these receptors, providing insight into ion permeation and lipid binding. However, the process by which α7 nAChRs transition between states remains unclear. To understand gating and lipid modulation, we generated two ensembles of molecular dynamics simulations of apo α7 nAChRs, with or without cholesterol. Using symmetry-adapted Markov state modeling, we developed a five-state gating model. Free energies recapitulated functional behavior, with the closed state dominating in absence of agonist. Open-to-nonconducting transition rates corresponded to experimental open durations. Cholesterol relatively stabilized the desensitized state, and reduced open-desensitized barriers. These results establish plausible asymmetric transition pathways between states, define lipid modulation effects on the α7 nAChR conformational cycle, and provide an ensemble of structural models applicable to rational design of lipidic pharmaceuticals.

Catalytic mechanism study of ATP-citrate lyase during citryl-CoA synthesis process

ATP-citrate lyase (ACLY) is a critical metabolic enzyme and promising target for drug development. The structure determinations of ACLY have revealed its homotetramer states with various subunit symmetries, but catalytic mechanism of ACLY tetramer and the importance of subunit symmetry have not been clarified. Here, we constructed the free energy landscape of ACLY tetramer with arbitrary subunit symmetries and investigated energetic and conformational coupling of subunits during citryl-CoA synthesis process. The optimal conformational pathway indicates that ACLY tetramer encounters three critical conformational barriers and undergoes a loss of rigid-D2 symmetry to gain an energetic advantage. Energetic coupling of conformational changes and biochemical reactions suggests that these biological events are not independent but rather coupled with each other, showing a comparable energy barrier to the experimental data for the rate-limiting step. These findings could contribute to further research on catalytic mechanism, functional modulation, and inhibitor design of ACLY.

Stable monomers in the ancestral sequence reconstruction of the last opisthokont common ancestor of dimeric triosephosphate isomerase

Function and structure are strongly coupled in obligated oligomers such as Triosephosphate isomerase (TIM). In animals and fungi, TIM monomers are inactive and unstable. Previously, we used ancestral sequence reconstruction to study TIM evolution and found that before these lineages diverged, the last opisthokonta common ancestor of TIM (LOCATIM) was an obligated oligomer that resembles those of extant TIMs. Notably, calorimetric evidence indicated that ancestral TIM monomers are more structured than extant ones. To further increase confidence about the function, structure, and stability of the LOCATIM, in this work, we applied two different inference methodologies and the worst plausible case scenario for both of them, to infer four sequences of this ancestor and test the robustness of their physicochemical properties. The extensive biophysical characterization of the four reconstructed sequences of LOCATIM showed very similar hydrodynamic and spectroscopic properties, as well as ligand-binding energetics and catalytic parameters. Their 3D structures were also conserved. Although differences were observed in melting temperature, all LOCATIMs showed reversible urea-induced unfolding transitions, and for those that reached equilibrium, high conformational stability was estimated (ΔGTot = 40.6-46.2 kcal/mol). The stability of the inactive monomeric intermediates was also high (ΔGunf = 12.6-18.4 kcal/mol), resembling some protozoan TIMs rather than the unstable monomer observed in extant opisthokonts. A comparative analysis of the 3D structure of ancestral and extant TIMs shows a correlation between the higher stability of the ancestral monomers with the presence of several hydrogen bonds located in the "bottom" part of the barrel.

Structural insights into rapamycin-induced oligomerization of a FRB-FKBP fusion protein

Inducible dimerization systems, such as rapamycin-induced dimerization of FK506 binding protein (FKBP) and FKBP-rapamycin binding (FRB) domain, are widely employed chemical biology tools to manipulate cellular functions. We previously advanced an inducible dimerization system into an inducible oligomerization system by developing a bivalent fusion protein, FRB-FKBP, which forms large oligomers upon rapamycin addition and can be used to manipulate cells. However, the oligomeric structure of FRB-FKBP remains unclear. Here, we report that FRB-FKBP forms a rotationally symmetric trimer in crystals, but a larger oligomer in solution, primarily tetramers and pentamers, which maintain similar inter-subunit contacts as in the crystal trimer. These findings expand the applications of the FRB-FKBP oligomerization system in diverse biological events.

De novo design of allosterically switchable protein assemblies

Allosteric modulation of protein function, wherein the binding of an effector to a protein triggers conformational changes at distant functional sites, plays a central part in the control of metabolism and cell signalling1-3. There has been considerable interest in designing allosteric systems, both to gain insight into the mechanisms underlying such 'action at a distance' modulation and to create synthetic proteins whose functions can be regulated by effectors4-7. However, emulating the subtle conformational changes distributed across many residues, characteristic of natural allosteric proteins, is a significant challenge8,9. Here, inspired by the classic Monod-Wyman-Changeux model of cooperativity10, we investigate the de novo design of allostery through rigid-body coupling of peptide-switchable hinge modules11 to protein interfaces12 that direct the formation of alternative oligomeric states. We find that this approach can be used to generate a wide variety of allosterically switchable systems, including cyclic rings that incorporate or eject subunits in response to peptide binding and dihedral cages that undergo effector-induced disassembly. Size-exclusion chromatography, mass photometry13 and electron microscopy reveal that these designed allosteric protein assemblies closely resemble the design models in both the presence and absence of peptide effectors and can have ligand-binding cooperativity comparable to classic natural systems such as haemoglobin14. Our results indicate that allostery can arise from global coupling of the energetics of protein substructures without optimized side-chain-side-chain allosteric communication pathways and provide a roadmap for generating allosterically triggerable delivery systems, protein nanomachines and cellular feedback control circuitry.

Structure, dimeric conformation, and coenzyme versatility of p-hydroxybenzoate hydroxylase from Arthrobacter sp. PAMC25564

p-Hydroxybenzoate hydroxylase (PHBH) catalyzes the ortho-hydroxylation of 4-hydroxybenzoate (4-HB) to protocatechuate (PCA). PHBHs are commonly known as homodimers, and the prediction of pyridine nucleotide binding and specificity remains an ongoing focus in this field. Therefore, our study aimed to determine the dimerization interface in AspPHBH from Arthrobacter sp. PAMC25564 and identify the canonical pyridine nucleotide-binding residues, along with coenzyme specificity, through site-directed mutagenesis. The results confirm a functional dimeric assembly from a tetramer that appeared in the crystallographic asymmetric unit identical to that established in previous studies. Furthermore, AspPHBH exhibits coenzyme versatility, utilizing both NADH and NADPH, with a preference for NADH. Rational engineering experiments demonstrated that targeted mutations in coenzyme surrounding residues profoundly impact NADPH binding, leading to nearly abrogated enzymatic activity compared to that of NADH. R50, R273, and S166 emerged as significant residues for NAD(P)H binding, having a near-fatal impact on NADPH binding compared to NADH. Likewise, the E44 residue plays a critical role in determining coenzyme specificity. Overall, our findings contribute to the fundamental understanding of the determinants of PHBH's active dimeric conformation, coenzyme binding and specificity holding promise for biotechnological advancements.

Hallmarks and evolutionary drivers of cotranslational protein complex assembly

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

Molecular Modeling of Single- and Double-Hydrocarbon-Stapled Coiled-Coil Inhibitors against Bcr-Abl: Toward a Treatment Strategy for CML

The chimeric oncoprotein Bcr-Abl is the causative agent of virtually all chronic myeloid leukemias and a subset of acute lymphoblastic leukemias. As a result of the so-called Philadelphia chromosome translocation t(9;22), Bcr-Abl manifests as a constitutively active tyrosine kinase, which promotes leukemogenesis by activation of cell cycle signaling pathways. Constitutive and oncogenic activation is mediated by an N-terminal coiled-coil oligomerization domain in Bcr (Bcr-CC), presenting a therapeutic target for inhibition of Bcr-Abl activity toward the treatment of Bcr-Abl+ leukemias. Previously, we demonstrated that a rationally designed Bcr-CC mutant, CCmut3, exerts a dominant negative effect upon Bcr-Abl activity by preferential oligomerization with Bcr-CC. Moreover, we have shown that conjugation to a leukemia-specific cell-penetrating peptide (CPP-CCmut3) improves intracellular delivery and activity. However, our full-length CPP-CCmut3 construct (81 aa) is encumbered by an intrinsically high degree of conformational variability and susceptibility to proteolytic degradation relative to traditional small-molecule therapeutics. Here, we iterate a new generation of CCmut3 inhibitors against Bcr-CC-mediated Bcr-Abl assembly designed to address these constraints through incorporation of all-hydrocarbon staples spanning i and i + 7 positions in α-helix 2 (CPP-CCmut3-st). We utilize computational modeling and biomolecular simulation to evaluate single- and double-stapled CCmut3 candidates in silico for dynamics and binding energetics. We further model a truncated system characterized by the deletion of α-helix 1 and the flexible loop linker, which are known to impart high conformational variability. To study the impact of the N-terminal cyclic CPP toward model stability and inhibitor activity, we also model the full-length and truncated systems devoid of the CPP, with a cyclized CPP, and with an open-configuration CPP, for a total of six systems that comprise our library. From this library, we present lead-stapled peptide candidates to be synthesized and evaluated experimentally as our next iteration of inhibitors against Bcr-Abl.

Frequent transitions in self-assembly across the evolution of a central metabolic enzyme

Many enzymes assemble into homomeric protein complexes comprising multiple copies of one protein. Because structural form is usually assumed to follow function in biochemistry, these assemblies are thought to evolve because they provide some functional advantage. In many cases, however, no specific advantage is known and, in some cases, quaternary structure varies among orthologs. This has led to the proposition that self-assembly may instead vary neutrally within protein families. The extent of such variation has been difficult to ascertain because quaternary structure has until recently been difficult to measure on large scales. Here, we employ mass photometry, phylogenetics, and structural biology to interrogate the evolution of homo-oligomeric assembly across the entire phylogeny of prokaryotic citrate synthases - an enzyme with a highly conserved function. We discover a menagerie of different assembly types that come and go over the course of evolution, including cases of parallel evolution and reversions from complex to simple assemblies. Functional experiments in vitro and in vivo indicate that evolutionary transitions between different assemblies do not strongly influence enzyme catalysis. Our work suggests that enzymes can wander relatively freely through a large space of possible assemblies and demonstrates the power of characterizing structure-function relationships across entire phylogenies.

Cracking the Code: Enhancing Molecular Tools for Progress in Nanobiotechnology

Nature continually refines its processes for optimal efficiency, especially within biological systems. This article explores the collaborative efforts of researchers worldwide, aiming to mimic nature's efficiency by developing smarter and more effective nanoscale technologies and biomaterials. Recent advancements highlight progress and prospects in leveraging engineered nucleic acids and proteins for specific tasks, drawing inspiration from natural functions. The focus is developing improved methods for characterizing, understanding, and reprogramming these materials to perform user-defined functions, including personalized therapeutics, targeted drug delivery approaches, engineered scaffolds, and reconfigurable nanodevices. Contributions from academia, government agencies, biotech, and medical settings offer diverse perspectives, promising a comprehensive approach to broad nanobiotechnology objectives. Encompassing topics from mRNA vaccine design to programmable protein-based nanocomputing agents, this work provides insightful perspectives on the trajectory of nanobiotechnology toward a future of enhanced biomimicry and technological innovation.

Translation-independent association of mRNAs encoding protomers of the 5-HT 2A -mGlu2 receptor complex in living cells

The serotonin 2A receptor (5-HT 2A R) and the metabotropic glutamate 2 receptor (mGluR2) form heteromeric G protein-coupled receptor (GPCR) complexes through a direct physical interaction. Co-translational association of mRNAs encoding subunits of heteromeric ion channels has been reported, but whether complex assembly of GPCRs occurs during translation remains unknown. Our in vitro data reveal evidence of co-translational modulation in 5-HT 2A R and mGluR2 mRNAs following siRNA-mediated knockdown. Interestingly, immunoprecipitation of either 5-HT 2A R or mGluR2, using an antibody targeting epitope tags at their N-terminus, results in detection of both transcripts associated with ribonucleoprotein complexes containing RPS24. Additionally, we demonstrate that the mRNA transcripts of 5-HT 2A R and mGluR2 associate autonomously of their respective encoded proteins. Validation of this translation-independent association is extended ex vivo using mouse frontal cortex samples. Together, these findings provide mechanistic insights into the co-translational assembly of GPCR heteromeric complexes, unraveling regulatory processes governing protein-protein interactions and complex formation.

MoLPC2: improved prediction of large protein complex structures and stoichiometry using Monte Carlo Tree Search and AlphaFold2

MOTIVATION

Today, the prediction of structures of large protein complexes solely from their sequence information requires prior knowledge of the stoichiometry of the complex. To address this challenge, we have enhanced the Monte Carlo Tree Search algorithms in MoLPC to enable the assembly of protein complexes while simultaneously predicting their stoichiometry.

RESULTS

In MoLPC2, we have improved the predictions by allowing sampling alternative AlphaFold predictions. Using MoLPC2, we accurately predicted the structures of 50 out of 175 nonredundant protein complexes (TM-score ≥ 0.8) without knowing the stoichiometry. MoLPC2 provides new opportunities for predicting protein complex structures without stoichiometry information.

AVAILABILITY AND IMPLEMENTATION

MoLPC2 is freely available at https://github.com/hychim/molpc2. A notebook is also available from the repository for easy use.

Review and Comparative Analysis of Methods and Advancements in Predicting Protein Complex Structure

Protein complexes perform diverse biological functions, and obtaining their three-dimensional structure is critical to understanding and grasping their functions. In many cases, it's not just two proteins interacting to form a dimer; instead, multiple proteins interact to form a multimer. Experimentally resolving protein complex structures can be quite challenging. Recently, there have been efforts and methods that build upon prior predictions of dimer structures to attempt to predict multimer structures. However, in comparison to monomeric protein structure prediction, the accuracy of protein complex structure prediction remains relatively low. This paper provides an overview of recent advancements in efficient computational models for predicting protein complex structures. We introduce protein-protein docking methods in detail and summarize their main ideas, applicable modes, and related information. To enhance prediction accuracy, other critical protein-related information is also integrated, such as predicting interchain residue contact, utilizing experimental data like cryo-EM experiments, and considering protein interactions and non-interactions. In addition, we comprehensively review computational approaches for end-to-end prediction of protein complex structures based on artificial intelligence (AI) technology and describe commonly used datasets and representative evaluation metrics in protein complexes. Finally, we analyze the formidable challenges faced in current protein complex structure prediction tasks, including the structure prediction of heteromeric complex, disordered regions in complex, antibody-antigen complex, and RNA-related complex, as well as the evaluation metrics for complex assessment. We hope that this work will provide comprehensive knowledge of complex structure predictions to contribute to future advanced predictions.

Photobleaching FRET-FLIM-ICS for quaternary structure quantification on cells. Theory and simulations

The oligomerization of proteins is an important biological control mechanism and has several functions in activity and stability of enzymes, structural proteins, ion channels and transcription factors. The determination of the relevant oligomeric states in terms of geometry (spatial extent), oligomer size (monomer or dimer or oligomer) and affinity (amounts of monomer, dimer and oligomer) is a challenging biophysical problem. Förster resonance energy transfer and fluorescence fluctuation spectroscopy are powerful tools that are sensitive to proximity and oligomerization respectively. Here it is proposed to combine image-based lifetime-detected Forster resonance energy transfer with image correlation spectroscopy and photobleaching to determine distances, oligomer sizes and oligomer distributions. Simulations for simple oligomeric forms illustrate the potential to improve the discrimination between different quaternary states in the cellular milieu.

Biochemical evidence for conformational variants in the anti-viral and pro-metastatic protein IFITM1

Interferon induced transmembrane proteins (IFITMs) play a dual role in the restriction of RNA viruses and in cancer progression, yet the mechanism of their action remains unknown. Currently, there is no data about the basic biochemical features or biophysical properties of the IFITM1 protein. In this work, we report on description and biochemical characterization of three conformational variants/oligomeric species of recombinant IFITM1 protein derived from an Escherichia coli expression system. The protein was extracted from the membrane fraction, affinity purified, and separated by size exclusion chromatography where two distinct oligomeric species were observed in addition to the expected monomer. These species remained stable upon re-chromatography and were designated as "dimer" and "oligomer" according to their estimated molecular weight. The dimer was found to be less stable compared to the oligomer using circular dichroism thermal denaturation and incubation with a reducing agent. A two-site ELISA and HDX mass spectrometry suggested the existence of structural motif within the N-terminal part of IFITM1 which might be significant in oligomer formation. Together, these data show the unusual propensity of recombinant IFITM1 to naturally assemble into very stable oligomeric species whose study might shed light on IFITM1 anti-viral and pro-oncogenic functions in cells.

Structural and functional analysis of the key enzyme responsible for the degradation of ochratoxin A in the Alcaligenes genus

The potential to degrade ochratoxin A (OTA), a highly poisonous mycotoxin, was investigated in cultures from Alcaligenes-type strains. Genome sequence analyses from different Alcaligenes species have permitted us to demonstrate a direct, causal link between the gene coding a known N-acyl-L-amino acid amidohydrolase from A. faecalis (AfOTH) and the OTA-degrading activity of this bacterium. In agreement with this finding, we found the gene coding AfOTH in two additional species included in the Alcaligenes genus, namely, A. pakistanensis, and A. aquatilis, which also degraded OTA. Notably, A. faecalis subsp. faecalis DSM 30030T was able to transform OTα, the product of OTA hydrolysis. AfOTH from A. faecalis subsp. phenolicus DSM 16503T was recombinantly over-produced and enzymatically characterized. AfOTH is a Zn2+-containing metalloenzyme that possesses structural features and conserved residues identified in the M20D family of enzymes. AfOTH is a tetramer in solution that shows both aminoacylase and carboxypeptidase activities. Using diverse potential substrates, namely, N-acetyl-L-amino acids and carbobenzyloxy-L-amino acids, a marked preference towards C-terminal Phe and Tyr residues could be deduced. The structural basis for this specificity has been determined by in silico molecular docking analyses. The amidase activity of AfOTH on C-terminal Phe residues structurally supports its OTA and OTB degradation activity.

Collision induced unfolding and molecular dynamics simulations of norovirus capsid dimers reveal strain-specific stability profiles

Collision induced unfolding (CIU) is a method used with ion mobility mass spectrometry to examine protein structures and their stability. Such experiments yield information about higher order protein structures, yet are unable to provide details about the underlying processes. That information can however be provided using molecular dynamics simulations. Here, we investigate the gas-phase unfolding of norovirus capsid dimers from the Norwalk and Kawasaki strains by employing molecular dynamics simulations over a range of temperatures, representing different levels of activation, together with CIU experiments. The dimers have highly similar structures, but their CIU reveals different stability that can be explained by the different dynamics that arises in response to the activation seen in the simulations, including a part of the sequence with previously observed strain-specific dynamics in solution. Our findings show how similar protein variants can be examined using mass spectrometric techniques in conjunction with atomistic molecular dynamics simulations to reveal differences in stability as well as differences in how and where unfolding takes place upon activation.

Rapid and accurate prediction of protein homo-oligomer symmetry with Seq2Symm

The majority of proteins must form higher-order assemblies to perform their biological functions. Despite the importance of protein quaternary structure, there are few machine learning models that can accurately and rapidly predict the symmetry of assemblies involving multiple copies of the same protein chain. Here, we address this gap by training several classes of protein foundation models, including ESM-MSA, ESM2, and RoseTTAFold2, to predict homo-oligomer symmetry. Our best model named Seq2Symm, which utilizes ESM2, outperforms existing template-based and deep learning methods. It achieves an average PR-AUC of 0.48 and 0.44 across homo-oligomer symmetries on two different held-out test sets compared to 0.32 and 0.23 for the template-based method. Because Seq2Symm can rapidly predict homo-oligomer symmetries using a single sequence as input (~ 80,000 proteins/hour), we have applied it to 5 entire proteomes and ~ 3.5 million unlabeled protein sequences to identify patterns in protein assembly complexity across biological kingdoms and species.

EncoMPASS: An encyclopedia of membrane proteins analyzed by structure and symmetry

Protein structure determination and prediction, active site detection, and protein sequence alignment techniques all exploit information about protein structure and structural relationships. For membrane proteins, however, there is limited agreement among available online tools for highlighting and mapping such structural similarities. Moreover, no available resource provides a systematic overview of quaternary and internal symmetries, and their orientation relative to the membrane, despite the fact that these properties can provide key insights into membrane protein function and evolution. Here, we describe the Encyclopedia of Membrane Proteins Analyzed by Structure and Symmetry (EncoMPASS), a database for relating integral membrane proteins of known structure from the points of view of sequence, structure, and symmetry. EncoMPASS is accessible through a web interface, and its contents can be easily downloaded. This allows the user not only to focus on specific proteins, but also to study general properties of the structure and evolution of membrane proteins.

Protein language model-embedded geometric graphs power inter-protein contact prediction

Accurate prediction of contacting residue pairs between interacting proteins is very useful for structural characterization of protein-protein interactions. Although significant improvement has been made in inter-protein contact prediction recently, there is still a large room for improving the prediction accuracy. Here we present a new deep learning method referred to as PLMGraph-Inter for inter-protein contact prediction. Specifically, we employ rotationally and translationally invariant geometric graphs obtained from structures of interacting proteins to integrate multiple protein language models, which are successively transformed by graph encoders formed by geometric vector perceptrons and residual networks formed by dimensional hybrid residual blocks to predict inter-protein contacts. Extensive evaluation on multiple test sets illustrates that PLMGraph-Inter outperforms five top inter-protein contact prediction methods, including DeepHomo, GLINTER, CDPred, DeepHomo2, and DRN-1D2D_Inter, by large margins. In addition, we also show that the prediction of PLMGraph-Inter can complement the result of AlphaFold-Multimer. Finally, we show leveraging the contacts predicted by PLMGraph-Inter as constraints for protein-protein docking can dramatically improve its performance for protein complex structure prediction.

Emergence of fractal geometries in the evolution of a metabolic enzyme

Fractals are patterns that are self-similar across multiple length-scales1. Macroscopic fractals are common in nature2-4; however, so far, molecular assembly into fractals is restricted to synthetic systems5-12. Here we report the discovery of a natural protein, citrate synthase from the cyanobacterium Synechococcus elongatus, which self-assembles into Sierpiński triangles. Using cryo-electron microscopy, we reveal how the fractal assembles from a hexameric building block. Although different stimuli modulate the formation of fractal complexes and these complexes can regulate the enzymatic activity of citrate synthase in vitro, the fractal may not serve a physiological function in vivo. We use ancestral sequence reconstruction to retrace how the citrate synthase fractal evolved from non-fractal precursors, and the results suggest it may have emerged as a harmless evolutionary accident. Our findings expand the space of possible protein complexes and demonstrate that intricate and regulatable assemblies can evolve in a single substitution.

Computational analysis of PD-L1 dimerization mechanism induced by small molecules and potential dynamical properties

The design of small molecule inhibitors that target the programmed death ligand-1 (PD-L1) is a forefront issue in immune checkpoint blocking therapy. Small-molecule inhibitors have been shown to exert therapeutic effects by inducing dimerization of the PD-L1 protein, however, the specific mechanisms underlying this dimerization process remain largely unexplored. Furthermore, there is a notable lack of comparative studies examining the binding modes of structurally diverse inhibitors. In view of the research gaps, this work employed molecular dynamics simulations to meticulously examine the interactions between two distinct types of inhibitors and PD-L1 in both monomeric and dimeric forms, and predicted the dimerization mechanism. The results revealed that inhibitors initially bind to a PD-L1 monomer, subsequently attracting another monomer to form a dimer. Notably, symmetric inhibitors observed superior binding efficiency compared to other inhibitors. Key residues, including Ile54, Tyr56, Met115 and Tyr123 played a leading role in binding. Structurally, symmetric inhibitors were capable of thoroughly engaging the binding pocket, promoting a more symmetrical formation of PD-L1 dimers. Furthermore, symmetric inhibitors formed more extensive hydrophobic interactions with protein residues. The insights garnered from this research are expected to significantly contribute to the rational design and optimization of small molecule inhibitors targeting PD-L1.

Understanding the complex formation of falstatin; an endogenous macromolecular inhibitor of falcipains

Proteolytic activity constitutes a fundamental process essential for the survival of the malaria parasite and is thus highly regulated. Falstatin, a protease inhibitor of Plasmodium falciparum, tightly regulates the activity of cysteine hemoglobinases, falcipain-2 and 3 (FP2, FP3), by inhibiting FP2 through a single surface exposed loop. However, the multimeric nature of falstatin and its interaction with FP2 remained unexplored. Here we report that the N-terminal falstatin region is highly disordered, and needs chaperone activity (heat-shock protein 70, HSP70) for its folding. Protein-protein interaction assays showed a significant interaction between falstatin and HSP70. Further, characterization of the falstatin multimer through a series of biophysical techniques identified the formation of a falstatin decamer, which was extremely thermostable. Computational analysis of the falstatin decamer showed the presence of five falstatin dimers, with each dimer aligned in a head-to-tail orientation. Further, the falstatin C-terminal region was revealed to be primarily involved in the oligomerization process. Stoichiometric analysis of the FP2-falstatin multimer showed the formation of a heterooligomeric complex in a 1:1 ratio, with the participation of ten subunits of each protein. Taken together, our results report a novel protease-inhibitor complex and strengthens our understanding of the regulatory mechanisms of major plasmodium hemoglobinases.

Each big journey starts with a first step: Importance of oligomerization

Protein oligomers, widely found in nature, have significant physiological and pathological functions. They are classified into three groups based on their function and toxicity. Significant advancements are being achieved in the development of functional oligomers, with a focus on various applications and their engineering. The antimicrobial peptides oligomers play roles in death of bacterial and cancer cells. The predominant pathogenic species in neurodegenerative disorders, as shown by recent results, are amyloid oligomers, which are the main subject of this chapter. They are generated throughout the aggregation process, serving as both intermediates in the subsequent aggregation pathways and ultimate products. Some of them may possess potent cytotoxic properties and through diverse mechanisms cause cellular impairment, and ultimately, the death of cells and disease progression. Information regarding their structure, formation mechanism, and toxicity is limited due to their inherent instability and structural variability. This chapter aims to provide a concise overview of the current knowledge regarding amyloid oligomers.

Synonymous Mutations Can Alter Protein Dimerization Through Localized Interface Misfolding Involving Self-entanglements

Synonymous mutations in messenger RNAs (mRNAs) can reduce protein-protein binding substantially without changing the protein's amino acid sequence. Here, we use coarse-grain simulations of protein synthesis, post-translational dynamics, and dimerization to understand how synonymous mutations can influence the dimerization of two E. coli homodimers, oligoribonuclease and ribonuclease T. We synthesize each protein from its wildtype, fastest- and slowest-translating synonymous mRNAs in silico and calculate the ensemble-averaged interaction energy between the resulting dimers. We find synonymous mutations alter oligoribonuclease's dimer properties. Relative to wildtype, the dimer interaction energy becomes 4% and 10% stronger, respectively, when translated from its fastest- and slowest-translating mRNAs. Ribonuclease T dimerization, however, is insensitive to synonymous mutations. The structural and kinetic origin of these changes are misfolded states containing non-covalent lasso-entanglements, many of which structurally perturb the dimer interface, and whose probability of occurrence depends on translation speed. These entangled states are kinetic traps that persist for long time scales. Entanglements cause altered dimerization energies for oligoribonuclease, as there is a large association (odds ratio: 52) between the co-occurrence of non-native self-entanglements and weak-binding dimer conformations. Simulated at all-atom resolution, these entangled structures persist for long timescales, indicating the conclusions are independent of model resolution. Finally, we show that regions of the protein we predict to have changes in entanglement are also structurally perturbed during refolding, as detected by limited-proteolysis mass spectrometry. Thus, non-native changes in entanglement at dimer interfaces is a mechanism through which oligomer structure and stability can be altered.

Molecular basis for curvature formation in SepF polymerization

The self-assembly of proteins into curved structures plays an important role in many cellular processes. One good example of this phenomenon is observed in the septum-forming protein (SepF), which forms polymerized structures with uniform curvatures. SepF is essential for regulating the thickness of the septum during bacteria cell division. In Bacillus subtilis, SepF polymerization involves two distinct interfaces, the β-β and α-α interfaces, which define the assembly unit and contact interfaces, respectively. However, the mechanism of curvature formation in this step is not yet fully understood. In this study, we employed solid-state NMR (SSNMR) to compare the structures of cyclic wild-type SepF assemblies with linear assemblies resulting from a mutation of G137 on the β-β interface. Our results demonstrate that while the sequence differences arise from the internal assembly unit, the dramatic changes in the shape of the assemblies depend on the α-α interface between the units. We further provide atomic-level insights into how the angular variation of the α2 helix on the α-α interface affects the curvature of the assemblies, using a combination of SSNMR, cryo-electron microscopy, and simulation methods. Our findings shed light on the shape control of protein assemblies and emphasize the importance of interhelical contacts in retaining curvature.

Simple visualization of submicroscopic protein clusters with a phase-separation-based fluorescent reporter

Protein clustering plays numerous roles in cell physiology and disease. However, protein oligomers can be difficult to detect because they are often too small to appear as puncta in conventional fluorescence microscopy. Here, we describe a fluorescent reporter strategy that detects protein clusters with high sensitivity called CluMPS (clusters magnified by phase separation). A CluMPS reporter detects and visually amplifies even small clusters of a binding partner, generating large, quantifiable fluorescence condensates. We use computational modeling and optogenetic clustering to demonstrate that CluMPS can detect small oligomers and behaves rationally according to key system parameters. CluMPS detected small aggregates of pathological proteins where the corresponding GFP fusions appeared diffuse. CluMPS also detected and tracked clusters of unmodified and tagged endogenous proteins, and orthogonal CluMPS probes could be multiplexed in cells. CluMPS provides a powerful yet straightforward approach to observe higher-order protein assembly in its native cellular context. A record of this paper's transparent peer review process is included in the supplemental information.

An atlas of protein homo-oligomerization across domains of life

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

Assembly Requirements for the Construction of Large-Scale Binary Protein Structures

The precise assembly of multiple biomacromolecules into well-defined structures and materials is of great importance for various biomedical and nanobiotechnological applications. In this study, we investigate the assembly requirements for two-component materials using charged protein nanocages as building blocks. To achieve this, we designed several variants of ferritin nanocages to determine the surface characteristics necessary for the formation of large-scale binary three-dimensional (3D) assemblies. These nanocage variants were employed in protein crystallization experiments and macromolecular crystallography analyses, complemented by computational methods. Through the screening of nanocage variant combinations at various ionic strengths, we identified three essential features for successful assembly: (1) the presence of a favored crystal contact region, (2) the presence of a charged patch not involved in crystal contacts, and (3) sufficient distinctiveness between the nanocages. Surprisingly, the absence of noncrystal contact mediating patches had a detrimental effect on the assemblies, highlighting their unexpected importance. Intriguingly, we observed the formation of not only binary structures but also both negatively and positively charged unitary structures under previously exclusively binary conditions. Overall, our findings will inform future design strategies by providing some design rules, showcasing the utility of supercharging symmetric building blocks in facilitating the assembly of biomacromolecules into large-scale binary 3D assemblies.

Structural and biochemical investigations of a HEAT-repeat protein involved in the cytosolic iron-sulfur cluster assembly pathway

Iron-sulfur clusters are essential for life and defects in their biosynthesis lead to human diseases. The mechanism of cluster assembly and delivery to cytosolic and nuclear client proteins via the cytosolic iron-sulfur cluster assembly (CIA) pathway is not well understood. Here we report cryo-EM structures of the HEAT-repeat protein Met18 from Saccharomyces cerevisiae, a key component of the CIA targeting complex (CTC) that identifies cytosolic and nuclear client proteins and delivers a mature iron-sulfur cluster. We find that in the absence of other CTC proteins, Met18 adopts tetrameric and hexameric states. Using mass photometry and negative stain EM, we show that upon the addition of Cia2, these higher order oligomeric states of Met18 disassemble. We also use pulldown assays to identify residues of critical importance for Cia2 binding and recognition of the Leu1 client, many of which are buried when Met18 oligomerizes. Our structures show conformations of Met18 that have not been previously observed in any Met18 homolog, lending support to the idea that a highly flexible Met18 may be key to how the CTC is able to deliver iron-sulfur clusters to client proteins of various sizes and shapes, i.e. Met18 conforms to the dimensions needed.

Accurate prediction of protein assembly structure by combining AlphaFold and symmetrical docking

AlphaFold can predict the structures of monomeric and multimeric proteins with high accuracy but has a limit on the number of chains and residues it can fold. Here we show that a combination of AlphaFold and all-atom symmetric docking simulations enables highly accurate prediction of the structure of complex symmetrical assemblies. We present a method to predict the structure of complexes with cubic - tetrahedral, octahedral and icosahedral - symmetry from sequence. Focusing on proteins where AlphaFold can make confident predictions on the subunit structure, 27 cubic systems were assembled with a median TM-score of 0.99 and a DockQ score of 0.72. 21 had TM-scores of above 0.9 and were categorized as acceptable- to high-quality according to DockQ. The resulting models are energetically optimized and can be used for detailed studies of intermolecular interactions in higher-order symmetrical assemblies. The results demonstrate how explicit treatment of structural symmetry can significantly expand the size and complexity of AlphaFold predictions.