Mechanisms of protein oligomerization, the critical role of insertions and deletions in maintaining different oligomeric states

Deciphering Pathways and Thermodynamics of Protein Assembly Using Native Mass Spectrometry

Protein oligomerization regulates many critical physiological processes, and its dysregulation can contribute to dysfunction and diseases. Elucidating the assembly pathways and quantifying their underlying thermodynamic and kinetic parameters are crucial for a comprehensive understanding of biological processes and for advancing therapeutics targeting abnormal protein oligomerization. Established binding assays, with limited mass precision, often rely on simplified models for data interpretation. In contrast, high-resolution native mass spectrometry (nMS) can directly determine the stoichiometry of biomolecular complexes in vitro. However, quantification is hindered by the fact that the relative abundances of gas-phase ions generally do not reflect solution concentrations due to nonuniform response factors. Recently, slow mixing mode (SLOMO)-nMS, which can quantify the relative response factors of interacting species, has been demonstrated to reliably measure the affinity (Kd) of binary biomolecular complexes. Here, we introduce an extended form of SLOMO-nMS that enables simultaneous quantification of the thermodynamics in multistep association reactions. Application of this method to homo-oligomerization of concanavalin A and insulin confirmed the reliability of the assay and uncovered details about the assembly processes that had previously resisted elucidation. Results acquired using SLOMO-nMS implemented with charge detection shed new light on the binding of recombinant human angiotensin-converting enzyme 2 and the SARS-CoV-2 spike protein. Importantly, new assembly pathways were uncovered, and the affinities of these interactions, which regulate host cell infection, were quantified. Together, these findings highlight the tremendous potential of SLOMO-nMS to accelerate the characterization of protein assembly pathways and thermodynamics and, in so doing, enhance fundamental biological understanding and facilitate therapeutic development. https://orcid.org/0000-0002-3389-7112.

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.

Characterization of a unique pH-dependent amylosucrase from Deinococcus cellulosilyticus

The amylosucrase (ASase, EC 2.4.1.4) utilizes sucrose as the sole substrate to catalyze multifunctional reactions. It can naturally synthesize α-1,4-linked glucans such as amylose as well as sucrose isomers with more favorable properties than sucrose with a lower intestinal digestibility and non-cariogenic properties. The amino acid sequence of the asase gene from Deinococcus cellulosilyticus (DceAS) exhibits low homology with those of other ASases from other Deinococcus species. In this study, we cloned and expressed DceAS and demonstrated its high activity at pH 6 and pH 8 and maintained stability. It showed higher polymerization activity at pH 6 than at pH 8, but similar isomerization activity and produced more turanose and trehalulose at pH 6 than at pH 8 and produced more isomaltulose at pH 8. Furthermore, the molecular weight of DceAS was 226.6 kDa at pH 6 and 145.5 kDa at pH 8, indicating that it existed as a trimer and dimer, respectively under those conditions. Additionally, circular dichroism spectra showed that the DceAS secondary structure was different at pH 6 and pH 8. These differences in reaction products at different pHs can be harnessed to naturally produce sucrose alternatives that are more beneficial to human health.

Multi-oligomeric and catalytically compromised serine acetyltransferase and cysteine regulatory complex of Mycobacterium tuberculosis

l-cysteine, a primary building block of mycothiol, plays an essential role in the defense mechanism of Mycobacterium tuberculosis (Mtb). However, it is unclear how Mtb regulates cysteine biosynthesis as no study has reported the cysteine regulatory complex (CRC) in Mtb. Serine acetyltransferase (SAT) and cysteine synthase (CS) interact to form CRC. Although MtCS has been characterized well, minimal information is available on MtSAT, which synthesizes, O-acetylserine (OAS), the precursor of cysteine. This study fills the gap and provides experimental evidence for the presence of MtCRC and a non-canonical multi-oligomeric MtSAT. We employed multiple analytical methods to characterize the oligomeric and kinetic properties of MtSAT and MtCRC. Results show that MtSAT, lacking >75 N-terminal amino acids exists in three different assembly states; trimer, hexamer, and dodecamer, compared to the single hexameric state of SAT of other bacteria. While hexamers display the highest catalytic turnover, the trimer is the least active. The predominance of trimers at low physiologically relevant concentrations suggests that MtSAT displays the lowest catalytic potential known. Further, the catalytic potential of MtSAT is also significantly reduced in CRC state, in contrast to enhanced activity of SAT in CRC of other organisms. Our study provides insights into multi-oligomeric MtSAT with reduced catalytic potential and demonstrates that both MtSAT and MtCS of Mycobacterium interact to form CRC, although with altered catalytic properties. We discuss our results in light of the altered biochemistry of the last step of canonical sulfate-dependent cysteine biosynthesis of Mycobacterium.

Rationally Engineering pH Adaptation of Acid-Induced Arginine Decarboxylase from Escherichia coli to Alkaline Environments to Efficiently Biosynthesize Putrescine

Acid-induced arginine decarboxylase AdiA is a typical homo-oligomeric protein biosynthesizing alkaline nylon monomer putrescine. However, upon loss of the AdiA decamer oligomeric state at neutral and alkaline conditions the activity also diminishes, obstructing the whole-cell biosynthesis of alkaline putrescine. Here, a structure cohesion strategy is proposed to change the pH adaptation of AdiA to alkaline environments based on the rational engineering of meridional and latitudinal oligomerization interfaces. After integrating substitutions of E467K at the latitudinal interface and H736E at the meridional channel interface, the structural stability of AdiA decamer and its substrate transport efficiency at neutral and alkaline conditions are improved. Finally, E467K_H736E is well adapted to neutral and alkaline environments (pH 7.0-9.0), and its enzymatic activity is 35-fold higher than that of wild AdiA at pH 8.0. Using E467K_H736E in the putrescine synthesis pathway, the titer of putrescine is up to 128.9 g·L-1 with a conversion of 0.94 mol·mol-1 in whole-cell catalysis. Additionally, the neutral pH adaptation of lysine decarboxylase, with a decamer structure similar to AdiA, is also improved using this cohesion strategy, providing an option for pH-adaptation engineering of other oligomeric decarboxylases.

Structural insights into the ubiquitylation strategy of the oligomeric CRL2FEM1B E3 ubiquitin ligase

Cullin-RING E3 ubiquitin ligase (CRL) family members play critical roles in numerous biological processes and diseases including cancer and Alzheimer's disease. Oligomerization of CRLs has been reported to be crucial for the regulation of their activities. However, the structural basis for its regulation and mechanism of its oligomerization are not fully known. Here, we present cryo-EM structures of oligomeric CRL2FEM1B in its unneddylated state, neddylated state in complex with BEX2 as well as neddylated state in complex with FNIP1/FLCN. These structures reveal that asymmetric dimerization of N8-CRL2FEM1B is critical for the ubiquitylation of BEX2 while FNIP1/FLCN is ubiquitylated by monomeric CRL2FEM1B. Our data present an example of the asymmetric homo-dimerization of CRL. Taken together, this study sheds light on the ubiquitylation strategy of oligomeric CRL2FEM1B according to substrates with different scales.

Mechanistic insights into glycoside 3-oxidases involved in C-glycoside metabolism in soil microorganisms

C-glycosides are natural products with important biological activities but are recalcitrant to degradation. Glycoside 3-oxidases (G3Oxs) are recently identified bacterial flavo-oxidases from the glucose-methanol-coline (GMC) superfamily that catalyze the oxidation of C-glycosides with the concomitant reduction of O2 to H2O2. This oxidation is followed by C-C acid/base-assisted bond cleavage in two-step C-deglycosylation pathways. Soil and gut microorganisms have different oxidative enzymes, but the details of their catalytic mechanisms are largely unknown. Here, we report that PsG3Ox oxidizes at 50,000-fold higher specificity (kcat/Km) the glucose moiety of mangiferin to 3-keto-mangiferin than free D-glucose to 2-keto-glucose. Analysis of PsG3Ox X-ray crystal structures and PsG3Ox in complex with glucose and mangiferin, combined with mutagenesis and molecular dynamics simulations, reveal distinctive features in the topology surrounding the active site that favor catalytically competent conformational states suitable for recognition, stabilization, and oxidation of the glucose moiety of mangiferin. Furthermore, their distinction to pyranose 2-oxidases (P2Oxs) involved in wood decay and recycling is discussed from an evolutionary, structural, and functional viewpoint.

Diacylglycerol-dependent hexamers of the SNARE-assembling chaperone Munc13-1 cooperatively bind vesicles

Munc13-1 is essential for vesicle docking and fusion at the active zone of synapses. Here, we report that Munc13-1 self-assembles into molecular clusters within diacylglycerol-rich microdomains present in phospholipid bilayers. Although the copy number of Munc13-1 molecules in these clusters has a broad distribution, a systematic Poisson analysis shows that this is most likely the result of two molecular species: monomers and mainly hexameric oligomers. Each oligomer is able to capture one vesicle independently. Hexamers have also been observed in crystals of Munc13-1 that form between opposed phospholipid bilayers [K. Grushin, R. V. Kalyana Sundaram, C. V. Sindelar, J. E. Rothman, Proc. Natl. Acad. Sci. U.S.A. 119, e2121259119 (2022)]. Mutations targeting the contacts stabilizing the crystallographic hexagons also disrupt the isolated hexamers, suggesting they are identical. Additionally, these mutations also convert vesicle binding from a cooperative to progressive mode. Our study provides an independent approach showing that Munc13-1 can form mainly hexamers on lipid bilayers each capable of vesicle capture.

A morpheein equilibrium regulates catalysis in phosphoserine phosphatase SerB2 from Mycobacterium tuberculosis

Mycobacterium tuberculosis phosphoserine phosphatase MtSerB2 is of interest as a new antituberculosis target due to its essential metabolic role in L-serine biosynthesis and effector functions in infected cells. Previous works indicated that MtSerB2 is regulated through an oligomeric transition induced by L-Ser that could serve as a basis for the design of selective allosteric inhibitors. However, the mechanism underlying this transition remains highly elusive due to the lack of experimental structural data. Here we describe a structural, biophysical, and enzymological characterisation of MtSerB2 oligomerisation in the presence and absence of L-Ser. We show that MtSerB2 coexists in dimeric, trimeric, and tetrameric forms of different activity levels interconverting through a conformationally flexible monomeric state, which is not observed in two near-identical mycobacterial orthologs. This morpheein behaviour exhibited by MtSerB2 lays the foundation for future allosteric drug discovery and provides a starting point to the understanding of its peculiar multifunctional moonlighting properties.

Salmonella YqiC exerts its function through an oligomeric state

Protein oligomerization occurs frequently both in vitro and in vivo, with specific functionalities associated with different oligomeric states. The YqiC protein from Salmonella Typhimurium forms a homotrimer through its C-terminal coiled-coil domain, and the protein is closely linked to the colonization and invasion of the bacteria to the host cells. To elucidate the importance of the oligomeric state of YqiC in vivo and its relation with bacterial infection, we mutated crucial residues in YqiC's coiled-coil region and confirmed the loss of trimer formation using chemical crosslinking and size exclusion chromatography coupled with multiple angle light scattering (SEC-MALS) techniques. The yqiC-knockout strain complemented with mutant YqiC showed significantly reduced colonization and invasion of Salmonella to host cells, demonstrating the critical role of YqiC oligomerization in bacterial pathogenesis. Furthermore, we conducted a protein-protein interaction study of YqiC using a pulled-down assay coupled with mass spectrometry analysis to investigate the protein's role in bacterial virulence. The results reveal that YqiC interacts with subunits of Complex II of the electron transport chain (SdhA and SdhB) and the β-subunit of F0 F1 -ATP synthase. These interactions suggest that YqiC may modulate the energy production of Salmonella and subsequently affect the assembly of crucial virulence factors, such as flagella. Overall, our findings provide new insights into the molecular mechanisms of YqiC's role in S. Typhimurium pathogenesis and suggest potential therapeutic targets for bacterial infections.

What's the defect? Using mass defects to study oligomerization of membrane proteins and peptides in nanodiscs with native mass spectrometry

Many membrane proteins form functional complexes that are either homo- or hetero-oligomeric. However, it is challenging to characterize membrane protein oligomerization in intact lipid bilayers, especially for polydisperse mixtures. Native mass spectrometry of membrane proteins and peptides inserted in lipid nanodiscs provides a unique method to study the oligomeric state distribution and lipid preferences of oligomeric assemblies. To interpret these complex spectra, we developed novel data analysis methods using macromolecular mass defect analysis. Here, we provide an overview of how mass defect analysis can be used to study oligomerization in nanodiscs, discuss potential limitations in interpretation, and explore strategies to resolve these ambiguities. Finally, we review recent work applying this technique to studying formation of antimicrobial peptide, amyloid protein, and viroporin complexes with lipid membranes.

Designing Cyclic-Constrained Peptides to Inhibit Human Phosphoglycerate Dehydrogenase

Although loop epitopes at protein-protein binding interfaces often play key roles in mediating oligomer formation and interaction specificity, their binding sites are underexplored as drug targets owing to their high flexibility, relatively few hot spots, and solvent accessibility. Prior attempts to develop molecules that mimic loop epitopes to disrupt protein oligomers have had limited success. In this study, we used structure-based approaches to design and optimize cyclic-constrained peptides based on loop epitopes at the human phosphoglycerate dehydrogenase (PHGDH) dimer interface, which is an obligate homo-dimer with activity strongly dependent on the oligomeric state. The experimental validations showed that these cyclic peptides inhibit PHGDH activity by directly binding to the dimer interface and disrupting the obligate homo-oligomer formation. Our results demonstrate that loop epitope derived cyclic peptides with rationally designed affinity-enhancing substitutions can modulate obligate protein homo-oligomers, which can be used to design peptide inhibitors for other seemingly intractable oligomeric proteins.

A mutation-induced drug resistance database (MdrDB)

Mutation-induced drug resistance is a significant challenge to the clinical treatment of many diseases, as structural changes in proteins can diminish drug efficacy. Understanding how mutations affect protein-ligand binding affinities is crucial for developing new drugs and therapies. However, the lack of a large-scale and high-quality database has hindered the research progresses in this area. To address this issue, we have developed MdrDB, a database that integrates data from seven publicly available datasets, which is the largest database of its kind. By integrating information on drug sensitivity and cell line mutations from Genomics of Drug Sensitivity in Cancer and DepMap, MdrDB has substantially expanded the existing drug resistance data. MdrDB is comprised of 100,537 samples of 240 proteins (which encompass 5119 total PDB structures), 2503 mutations, and 440 drugs. Each sample brings together 3D structures of wild type and mutant protein-ligand complexes, binding affinity changes upon mutation (ΔΔG), and biochemical features. Experimental results with MdrDB demonstrate its effectiveness in significantly enhancing the performance of commonly used machine learning models when predicting ΔΔG in three standard benchmarking scenarios. In conclusion, MdrDB is a comprehensive database that can advance the understanding of mutation-induced drug resistance, and accelerate the discovery of novel chemicals.

The role of native cysteine residues in the oligomerization of KCNQ1 channels

KCNQ1, the major component of the slow-delayed rectifier potassium channel, is responsible for repolarization of cardiac action potential. Mutations in this channel can lead to a variety of diseases, most notably long QT syndrome. It is currently unknown how many of these mutations change channel function and structure on a molecular level. Since tetramerization is key to proper function and structure of the channel, it is likely that mutations modify the stability of KCNQ1 oligomers. Presently, the C-terminal domain of KCNQ1 has been noted as the driving force for oligomer formation. However, truncated versions of this protein lacking the C-terminal domain still tetramerize. Therefore, we explored the role of native cysteine residues in a truncated construct of human KCNQ1, amino acids 100-370, by blocking potential interactions of cysteines with a nitroxide based spin label. Mobility of the spin labels was investigated with continuous wave electron paramagnetic resonance (CW-EPR) spectroscopy. The oligomerization state was examined by gel electrophoresis. The data provide information on tetramerization of human KCNQ1 without the C-terminal domain. Specifically, how blocking the side chains of native cysteines residues reduces oligomerization. A better understanding of tetramer formation could provide improved understanding of the molecular etiology of long QT syndrome and other diseases related to KCNQ1.

Structural Dynamics Predominantly Determine the Adaptability of Proteins to Amino Acid Deletions

The insertion or deletion (indel) of amino acids has a variety of effects on protein function, ranging from disease-forming changes to gaining new functions. Despite their importance, indels have not been systematically characterized towards protein engineering or modification goals. In the present work, we focus on deletions composed of multiple contiguous amino acids (mAA-dels) and their effects on the protein (mutant) folding ability. Our analysis reveals that the mutant retains the native fold when the mAA-del obeys well-defined structural dynamics properties: localization in intrinsically flexible regions, showing low resistance to mechanical stress, and separation from allosteric signaling paths. Motivated by the possibility of distinguishing the features that underlie the adaptability of proteins to mAA-dels, and by the rapid evaluation of these features using elastic network models, we developed a positive-unlabeled learning-based classifier that can be adopted for protein design purposes. Trained on a consolidated set of features, including those reflecting the intrinsic dynamics of the regions where the mAA-dels occur, the new classifier yields a high recall of 84.3% for identifying mAA-dels that are stably tolerated by the protein. The comparative examination of the relative contribution of different features to the prediction reveals the dominant role of structural dynamics in enabling the adaptation of the mutant to mAA-del without disrupting the native fold.

Subunit association, and thermal and chemical unfolding of Cucurbitaceae phloem exudate lectins. A review

Detailed characterization of protein (un)folding intermediates is crucial for understanding the (un)folding pathway, aggregation, stability and their functional properties. In recent years, stress-inducible lectins are being investigated with much interest. In plants phloem proteins PP1 and PP2 are major components of the phloem fluid. While PP1 is a structural protein, PP2 exhibits lectin activity, and was proposed to play key roles in wound sealing, anti-pathogenic activity, and transportation of various molecules including RNA within the plant. Cucurbitaceae fruits contain high concentrations of PP2 lectins, which recognize chitooligosaccharides with high specificity. Although the presence of PP2 lectins in the phloem exudate of Cucurbitaceae species was documented over 40 years ago, so far only a few proteins from this family have been purified and characterized in detail. This review summarizes the results of biophysical studies aimed at investigating the oligomeric status of these lectins, their thermal stability, structural perturbations caused by changes in pH and addition of chaotropic agents and characterization of intermediates observed in the unfolding process. The implications of these results in the functional roles played by PP2 type lectins in their native environment are discussed. Finally, perspectives for future biophysical research on these proteins are given.

Structural evolution of an amphibian-specific globin: A computational evolutionary biochemistry approach

Studies on the globin family are continuously revealing insights into the mechanisms of gene and protein evolution. The rise of a new globin gene type in Pelobatoidea and Neobatrachia (Amphibia:Anura) from an α-globin precursor provides the opportunity to investigate the genetic and physical mechanisms underlying the origin of new protein structural and functional properties. This amphibian-specific globin (globin A/GbA) discovered in the heart of Rana catesbeiana is a monomer. As the ancestral oligomeric state of α-globins is a homodimer, we inferred that the ancestral state was lost somewhere in the GbA lineage. Here, we combined computational molecular evolution with structural bioinformatics to determine the extent to which the loss of the homodimeric state is pervasive in the GbA clade. We also characterized the loci of GbA genes in Bufo bufo. We found two GbA clades in Neobatrachia. One was deleted in Ranidae, but retained and expanded to yield a new globin cluster in Bufonidae species. Loss of the ancestral oligomeric state seems to be pervasive in the GbA clade. However, a taxonomic sampling that includes more Pelobatoidea, as well as early Neobatrachia, lineages would be necessary to determine the oligomeric state of the last common ancestor of all GbA. The evidence presented here points out a possible loss of oligomerization in Pelobatoidea GbA as a result of amino acid substitutions that weaken the homodimeric state. In contrast, the loss of oligomerization in both Neobatrachia GbA clades was linked to independent deletions that disrupted many packing contacts at the homodimer interface.

Assembly and disassembly dynamics of nonmuscle myosin II control endosomal fission

Endocytic vesicular trafficking requires merging of two lipid bilayers, but how the two lipid bilayers can come close together during fusion and fission in endocytic trafficking is not well explored. Here, we establish that knocking down nonmuscle myosin IIs (NM IIs) by small interfering RNA (siRNA) or inhibition of their activities by (-) blebbistatin causes the formation of a ring-like assembly of early endosomes (raEE). Inhibition of NM II assembly by an inhibitor of regulatory light-chain (RLC) kinase results in the formation of raEE, whereas inhibition of NM II disassembly by inhibitors of heavy chain kinases, protein kinase C (PKC) and casein kinase 2 (CK2), causes the dispersion of early endosomes. The raEEs retain EEA1, Rab7, and LAMP2 markers. Overexpression of an assembly incompetent form, RLC-AA, and disassembly incompetent form, NMHCIIB-S6A or NMHCIIA-1916A, induces such defects, respectively. Altogether, these data support that NM II assembly and disassembly dynamics participate in endocytic trafficking by regulating fission to maintain the size of early endosomes.

Insertions and Deletions (Indels): A Missing Piece of the Protein Engineering Jigsaw

Over the years, protein engineers have studied nature and borrowed its tricks to accelerate protein evolution in the test tube. While there have been considerable advances, our ability to generate new proteins in the laboratory is seemingly limited. One explanation for these shortcomings may be that insertions and deletions (indels), which frequently arise in nature, are largely overlooked during protein engineering campaigns. The profound effect of indels on protein structures, by way of drastic backbone alterations, could be perceived as "saltation" events that bring about significant phenotypic changes in a single mutational step. Should we leverage these effects to accelerate protein engineering and gain access to unexplored regions of adaptive landscapes? In this Perspective, we describe the role played by indels in the functional diversification of proteins in nature and discuss their untapped potential for protein engineering, despite their often-destabilizing nature. We hope to spark a renewed interest in indels, emphasizing that their wider study and use may prove insightful and shape the future of protein engineering by unlocking unique functional changes that substitutions alone could never achieve.

Unravelling the regulation pathway of photosynthetic AB-GAPDH

Oxygenic phototrophs perform carbon fixation through the Calvin–Benson cycle. Different mechanisms adjust the cycle and the light-harvesting reactions to rapid environmental changes. Photosynthetic glyceraldehyde 3-phosphate dehydrogenase (GAPDH) is a key enzyme in the cycle. In land plants, different photosynthetic GAPDHs exist: the most abundant isoform is formed by A2B2 heterotetramers and the least abundant by A4 homotetramers. Regardless of the subunit composition, GAPDH is the major consumer of photosynthetic NADPH and its activity is strictly regulated. While A4-GAPDH is regulated by CP12, AB-GAPDH is autonomously regulated through the C-terminal extension (CTE) of its B subunits. Reversible inhibition of AB-GAPDH occurs via the oxidation of a cysteine pair located in the CTE and the substitution of NADP(H) with NAD(H) in the cofactor-binding site. These combined conditions lead to a change in the oligomerization state and enzyme inhibition. SEC–SAXS and single-particle cryo-EM analysis were applied to reveal the structural basis of this regulatory mechanism. Both approaches revealed that spinach (A2B2) n -GAPDH oligomers with n = 1, 2, 4 and 5 co-exist in a dynamic system. B subunits mediate the contacts between adjacent tetramers in A4B4 and A8B8 oligomers. The CTE of each B subunit penetrates into the active site of a B subunit of the adjacent tetramer, which in turn moves its CTE in the opposite direction, effectively preventing the binding of the substrate 1,3-bisphosphoglycerate in the B subunits. The whole mechanism is made possible, and eventually controlled, by pyridine nucleotides. In fact, NAD(H), by removing NADP(H) from A subunits, allows the entrance of the CTE into the active site of the B subunit, hence stabilizing inhibited oligomers.

MVsim is a toolset for quantifying and designing multivalent interactions

Arising through multiple binding elements, multivalency can specify the avidity, duration, cooperativity, and selectivity of biomolecular interactions, but quantitative prediction and design of these properties has remained challenging. Here we present MVsim, an application suite built around a configurational network model of multivalency to facilitate the quantification, design, and mechanistic evaluation of multivalent binding phenomena through a simple graphical user interface. To demonstrate the utility and versatility of MVsim, we first show that both monospecific and multispecific multivalent ligand-receptor interactions, with their noncanonical binding kinetics, can be accurately simulated. Further, to illustrate the conceptual insights into multivalent systems that MVsim can provide, we apply it to quantitatively predict the ultrasensitivity and performance of multivalent-encoded protein logic gates, evaluate the inherent programmability of multispecificity for selective receptor targeting, and extract rate constants of conformational switching for the SARS-CoV-2 spike protein and model its binding to ACE2 as well as multivalent inhibitors of this interaction. MVsim and instructional tutorials are freely available at https://sarkarlab.github.io/MVsim/ .

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.

How gene duplication diversifies the landscape of protein oligomeric state and function

Oligomeric proteins are central to cellular life and the duplication and divergence of their genes is a key driver of evolutionary innovations. The duplication of a gene coding for an oligomeric protein has numerous possible outcomes, which motivates questions on the relationship between structural and functional divergence. How do protein oligomeric states diversify after gene duplication? In the simple case of duplication of a homo-oligomeric protein gene, what properties can influence the fate of descendant paralogs toward forming independent homomers or maintaining their interaction as a complex? Furthermore, how are functional innovations associated with the diversification of oligomeric states? Here, we review recent literature and present specific examples in an attempt to illustrate and answer these questions.

Trigger factor both holds and folds its client proteins

ATP-independent chaperones like trigger factor are generally assumed to play passive roles in protein folding by acting as holding chaperones. Here we show that trigger factor plays a more active role. Consistent with a role as an aggregation inhibiting chaperone, we find that trigger factor rapidly binds to partially folded glyceraldehyde 3-phosphate dehydrogenase (GAPDH) and prevents it from non-productive self-association by shielding oligomeric interfaces. In the traditional view of holding chaperone action, trigger factor would then be expected to transfer its client to a chaperone foldase system for complete folding. Unexpectedly, we noticed that GAPDH folds into a monomeric but otherwise rather native-like intermediate state while trigger factor-bound. Upon release from trigger factor, the mostly folded monomeric GAPDH rapidly self-associates into its native tetramer and acquires enzymatic activity without needing additional folding factors. The mechanism we propose here for trigger factor bridges the holding and folding activities of chaperone function.

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.

A peptide interfering with the dimerization of oncogenic KITENIN protein and its stability suppresses colorectal tumour progression

The stability of a protein, as well as its function and versatility, can be enhanced through oligomerization. KITENIN (KAI1 C-terminal interacting tetraspanin) is known to promote the malignant progression of colorectal cancer (CRC). How KITENIN maintains its structural integrity and stability are largely unknown, however. Here we investigated the mechanisms regulating the stability of KITENIN with the aim of developing therapeutics blocking its oncogenic functions. We found that KITENIN formed a homo-oligomeric complex and that the intracellular C-terminal domain (KITENIN-CTD) was needed for this oligomerization. Expression of the KITENIN-CTD alone interfered with the formation of the KITENIN homodimer, and the amino acid sequence from 463 to 471 within the KITENIN-CTD was the most effective. This sequence coupled with a cell-penetrating peptide was named a KITENIN dimerization-interfering peptide (KDIP). We next studied the mechanisms by which KDIP affected the stability of KITENIN. The KITENIN-interacting protein myosin-X (Myo10), which has oncogenic activity in several cancers, functioned as an effector to stabilize the KITENIN homodimer in the cis formation. Treatment with KDIP resulted in the disintegration of the homodimer via downregulation of Myo10, which led to increased binding of RACK1 to the exposed RACK1-interacting motif (463-471 aa), and subsequent autophagy-dependent degradation of KITENIN and reduced CRC cell invasion. Intravenous injection of KDIP significantly reduced the tumour burden in a syngeneic mouse tumour model and colorectal liver metastasis in an intrasplenic hepatic metastasis model. Collectively, our present results provide a new cancer therapeutic peptide for blocking colorectal liver metastasis, which acts by inducing the downregulation of Myo10 and specifically targeting the stability of the oncogenic KITENIN protein.

Structural and biophysical characterization of transcription factor HNF-1A as a tool to study MODY3 diabetes variants

Hepatocyte nuclear factor 1A (HNF-1A) is a transcription factor expressed in several embryonic and adult tissues, modulating the expression of numerous target genes. Pathogenic variants in the HNF1A gene are known to cause maturity-onset diabetes of the young 3 (MODY3 or HNF1A MODY), a disease characterized by dominant inheritance, age of onset before 25 to 35 years of age, and pancreatic β-cell dysfunction. A precise diagnosis can alter management of this disease, as insulin can be exchanged with sulfonylurea tablets and genetic counseling differs from polygenic forms of diabetes. Therefore, more knowledge on the mechanisms of HNF-1A function and the level of pathogenicity of the numerous HNF1A variants is required for precise diagnostics. Here, we structurally and biophysically characterized an HNF-1A protein containing both the DNA-binding domain and the dimerization domain, and determined the folding and DNA-binding capacity of two established MODY3 HNF-1A variant proteins (P112L, R263C) and one variant of unknown significance (N266S). All three variants showed reduced functionality compared to the WT protein. Furthermore, while the R263C and N266S variants displayed reduced binding to an HNF-1A target promoter, we found the P112L variant was unstable in vitro and in cells. Our results support and mechanistically explain disease causality for these investigated variants and present a novel approach for the dissection of structurally unstable and DNA-binding defective variants. This study indicates that structural and biochemical investigation of HNF-1A is a valuable tool in reliable variant classification needed for precision diabetes diagnostics and management.

Conserved Molecular Signatures in the Spike, Nucleocapsid, and Polymerase Proteins Specific for the Genus Betacoronavirus and Its Different Subgenera

The genus Betacoronavirus, consisting of four main subgenera (Embecovirus, Merbecovirus, Nobecovirus, and Sarbecovirus), encompasses all clinically significant coronaviruses (CoVs), including SARS, MERS, and the SARS-CoV-2 virus responsible for current COVID-19 pandemic. Very few molecular characteristics are known that are specific for the genus Betacoronavirus or its different subgenera. In this study, our analyses of the sequences of four essential proteins of CoVs, viz., spike, nucleocapsid, envelope, and RNA-dependent RNA polymerase (RdRp), identified ten novel molecular signatures consisting of conserved signature indels (CSIs) in these proteins which are specific for the genus Betacoronavirus or its subgenera. Of these CSIs, two 14-aa-conserved deletions found within the heptad repeat motifs 1 and 2 of the spike protein are specific for all betacoronaviruses, except for their shared presence in the highly infectious avian coronavirus. Six additional CSIs present in the nucleocapsid protein and one CSI in the RdRp protein are distinctive characteristics of either the Merbecovirus, Nobecovirus, or Sarbecovirus subgenera. In addition, a 4-aa insert is present in the spike protein, which is uniquely shared by all viruses from the subgenera Merbecovirus, Nobecovirus, and Sarbecovirus, but absent in Embecovirus and all other genera of CoVs. This molecular signature provides evidence that viruses from the three subgenera sharing this CSI are more closely related to each other, and they evolved after the divergence of embecoviruses and other CoVs. As all CSIs specific for different groups of CoVs are flanked by conserved regions, their sequences provide novel means for identifying the above groups of CoVs and for developing novel diagnostic tests. Furthermore, our analyses of the structures of the spike and nucleocapsid proteins show that all identified CSIs are localized in the surface-exposed loops of these protein. It is postulated that these surface loops, through their interactions with other cellular proteins/ligands, play important roles in the biology/pathology of these viruses.

The tetrameric structure of Plasmodium falciparum phosphoglycerate mutase is critical for optimal enzymatic activity

The glycolytic enzyme phosphoglycerate mutase (PGM) is of utmost importance for overall cellular metabolism and has emerged as a novel therapeutic target in cancer cells. This enzyme is also conserved in the rapidly proliferating malarial parasite Plasmodium falciparum, which have a similar metabolic framework as cancer cells and rely on glycolysis as the sole energy-yielding process during intraerythrocytic development. There is no redundancy among the annotated PGM enzymes in Plasmodium, and PfPGM1 is absolutely required for the parasite survival as evidenced by conditional knockdown in our study. A detailed comparison of PfPGM1 with its counterparts followed by in-depth structure-function analysis revealed unique attributes of this parasitic protein. Here, we report for the first time the importance of oligomerization for the optimal functioning of the enzyme in vivo, as earlier studies in eukaryotes only focused on the effects in vitro. We show that single point mutation of the amino acid residue W68 led to complete loss of tetramerization and diminished catalytic activity in vitro. Additionally, ectopic expression of the WT PfPGM1 protein enhanced parasite growth, whereas the monomeric form of PfPGM1 failed to provide growth advantage. Furthermore, mutation of the evolutionarily conserved residue K100 led to a drastic reduction in enzymatic activity. The indispensable nature of this parasite enzyme highlights the potential of PfPGM1 as a therapeutic target against malaria, and targeting the interfacial residues critical for oligomerization can serve as a focal point for promising drug development strategies that may not be restricted to malaria only.

Conserved Signatures in Protein Sequences Reliably Demarcate Different Clades of Rodents/Glires Species and Consolidate Their Evolutionary Relationships

The grandorder Glires, consisting of the orders Rodentia and Lagomorpha, encompasses a significant portion of the extant mammalian species including Rat, Mouse, Squirrel, Guinea pig and Beaver. Glires species play an important role in the ecosystem and provide valuable animal models for genetic studies and animal testing. Thus, it is important to reliably determine their evolutionary relationships and identify molecular characteristics that are specific for different species groups within the Glires. In this work, we have constructed a phylogenetic tree for >30 genome sequenced Glires species based on concatenated sequences of 25 conserved proteins. In this tree, members of different orders, suborders, and families within Glires formed strongly supported clades, and their interrelationships were also generally reliably resolved. In parallel, we conducted comparative analyses on more than 1500 protein sequences from Glires species to identify highly conserved molecular markers. These markers were comprised of conserved signature indels (CSIs) in proteins, which are specific for different Rodentia/Glires clades. Of the 41 novel CSIs identified in this work, some are specific for the entire Glires, Rodentia, or Lagomorpha clades, whereas many others reliably demarcate different family/suborder level clades of Rodentia (viz. Myomorpha, Castorimorpha, Sciuromorpha, Hystricomorpha, and Muroidea). Additionally, some of the CSIs also provide information regarding the interrelationships among Rodentia subgroups. Our analysis has also identified one CSI that is commonly shared by the Glires and Scandentia species (tree shrew), however, its evolutionary significance is unclear. Several of the identifed rodents-specific CSIs are present in conserved disease-related proteins. Thus, they provide novel molecular markers for genetic and biochemical studies on the functions of these proteins.

Cryo-EM structure of acylpeptide hydrolase reveals substrate selection by multimerization and a multi-state serine-protease triad

The first structure of tetrameric mammalian acylaminoacyl peptidase, an enzyme that functions as an upstream regulator of the proteasome through the removal of terminal N-acetylated residues from its protein substrates, was determined by cryo-EM and further elucidated by MD simulations. Self-association results in a toroid-shaped quaternary structure, guided by an amyloidogenic β-edge and unique inserts. With a Pro introduced into its central β-sheet, sufficient conformational freedom is awarded to the segment containing the catalytic Ser587 that the serine protease catalytic triad alternates between active and latent states. Active site flexibility suggests that the dual function of catalysis and substrate selection are fulfilled by a novel mechanism: substrate entrance is regulated by flexible loops creating a double-gated channel system, while binding of the substrate to the active site is required for stabilization of the catalytic apparatus – as a second filter before hydrolysis. The structure not only underlines that within the family of S9 proteases homo-multimerization acts as a crucial tool for substrate selection, but it will also allow drug design targeting of the ubiquitin-proteasome system.

The structure of tetrameric mammalian acylaminoacyl peptidase – a key upstream regulator of the proteasome – was determined by cryo-EM (and elucidated by MD), showing a “shutters-and-channels” substrate selection apparatus created by oligomerization.

Conserved molecular signatures in the spike protein provide evidence indicating the origin of SARS-CoV-2 and a Pangolin-CoV (MP789) by recombination(s) between specific lineages of Sarbecoviruses

Both SARS-CoV-2 and SARS coronaviruses (CoVs) are members of the subgenus Sarbecovirus. To understand the origin of SARS-CoV-2, sequences for the spike and nucleocapsid proteins from sarbecoviruses were analyzed to identify molecular markers consisting of conserved inserts or deletions (termed CSIs) that are specific for either a particular clade of Sarbecovirus or are commonly shared by two or more clades of these viruses. Three novel CSIs in the N-terminal domain (NTD) of the spike protein S1-subunit (S1-NTD) are uniquely shared by SARS-CoV-2, Bat-CoV-RaTG13 and most pangolin CoVs (SARS-CoV-2r clade). Three other sarbecoviruses viz. bat-CoVZXC21, -CoVZC45 and -PrC31 (forming CoVZC/PrC31 clade), and a pangolin-CoV_MP789 also contain related CSIs in the same positions. In contrast to the S1-NTD, both SARS and SARS-CoV-2r viruses contain two large CSIs in the S1-C-terminal domain (S1-CTD) that are absent in the CoVZC/PrC31 clade. One of these CSIs, consisting of a 12 aa insert, is also present in the RShSTT clade (Cambodia-CoV strains). Sequence similarity studies show that the S1-NTD of SARS-CoV-2r viruses is most similar to the CoVZC/PrC31 clade, whereas their S1-CTD exhibits highest similarity to the RShSTT- (and the SARS-related) CoVs. Results from the shared presence of CSIs and sequence similarity studies on different CoV lineages support the inference that the SARS-CoV-2r cluster of viruses has originated by a genetic recombination between the S1-NTD of the CoVZC/PrC31 clade of CoVs and the S1-CTD of RShSTT/SARS viruses, respectively. We also present compelling evidence, based on the shared presence of CSIs and sequence similarity studies, that the pangolin-CoV_MP789, whose receptor-binding domain is most similar to the SARS-CoV-2 virus, has resulted from another independent recombination event involving the S1-NTD of the CoVZC/PrC31 CoVs and the S1-CTD of an unidentified SARS-CoV-2r related virus. The SARS-CoV-2 virus involved in this latter recombination event is postulated to be most similar to the SARS-CoV-2. Several other CSIs reported here are specific for other clusters of sarbecoviruses including a clade consisting of bat-SARS-CoVs (BM48-31/BGR/2008 and SARS_BtKY72). Structural mapping studies show that the identified CSIs form distinct loops/patches on the surface of the spike protein. It is hypothesized that these novel loops/patches on the spike protein, through their interactions with other host components, should play important roles in the biology/pathology of SARS-CoV-2 virus. Lastly, the CSIs specific for different clades of sarbecoviruses including SARS-CoV-2r clade provide novel means for the identification of these viruses and other potential applications.

Evolution of the folding landscape of effector caspases

Caspases are a family of cysteinyl proteases that control programmed cell death and maintain homeostasis in multicellular organisms. The caspase family is an excellent model to study protein evolution because all caspases are produced as zymogens (procaspases [PCPs]) that must be activated to gain full activity; the protein structures are conserved through hundreds of millions of years of evolution; and some allosteric features arose with the early ancestor, whereas others are more recent evolutionary events. The apoptotic caspases evolved from a common ancestor (CA) into two distinct subfamilies: monomers (initiator caspases) or dimers (effector caspases). Differences in activation mechanisms of the two subfamilies, and their oligomeric forms, play a central role in the regulation of apoptosis. Here, we examine changes in the folding landscape by characterizing human effector caspases and their CA. The results show that the effector caspases unfold by a minimum three-state equilibrium model at pH 7.5, where the native dimer is in equilibrium with a partially folded monomeric (PCP-7, CA) or dimeric (PCP-6) intermediate. In comparison, the unfolding pathway of PCP-3 contains both oligomeric forms of the intermediate. Overall, the data show that the folding landscape was first established with the CA and was retained for >650 million years. Partially folded monomeric or dimeric intermediates in the ancestral ensemble provide mechanisms for evolutionary changes that affect stability of extant caspases. The conserved folding landscape allows for the fine-tuning of enzyme stability in a species-dependent manner while retaining the overall caspase–hemoglobinase fold.

Structure-function relationship of terpenoid glycosyltransferases from plants

Covering: up to 2021Terpenoids are physiologically active substances that are of great importance to humans. Their physicochemical properties are modified by glycosylation, in terms of polarity, volatility, solubility and reactivity, and their bioactivities are altered accordingly. Significant scientific progress has been made in the functional study of glycosylated terpenes and numerous plant enzymes involved in regio- and enantioselective glycosylation have been characterized, a reaction that remains chemically challenging. Crucial clues to the mechanism of terpenoid glycosylation were recently provided by the first crystal structures of a diterpene glycosyltransferase UGT76G1. Here, we review biochemically characterized terpenoid glycosyltransferases, compare their functions and primary structures, discuss their acceptor and donor substrate tolerance and product specificity, and elaborate features of the 3D structures of the first terpenoid glycosyltransferases from plants.

MVsim : a toolset for quantifying and designing multivalent interactions

Arising through multiple binding elements, multivalency can specify the avidity, duration, cooperativity, and selectivity of biomolecular interactions, but quantitative prediction and design of these properties has remained challenging. Here we present MVsim , an application suite built around a configurational network model of multivalency to facilitate the quantification, design, and mechanistic evaluation of multivalent binding phenomena through a simple graphical user interface. To demonstrate the utility and versatility of MVsim , we first show that both monospecific and multispecific multivalent ligand-receptor interactions, with their noncanonical binding kinetics, can be accurately simulated. We then quantitatively predict the ultrasensitivity and performance of multivalent-encoded protein logic gates, evaluate the inherent programmability of multispecificity for selective receptor targeting, and extract rate constants of conformational switching for the SARS-CoV-2 spike protein and model its binding to ACE2 as well as multivalent inhibitors of this interaction. MVsim is freely available at https://sarkarlab.github.io/MVsim/ .

Differential scanning fluorimetry (DSF) screen to identify inhibitors of Hsp60 protein-protein interactions

There are relatively few methods available for discovering inhibitors of the protein-protein interactions (PPIs) that hold together homo-oligomers. We envisioned that Differential Scanning Fluorimetry (DSF) might be a versatile way to discover this type of inhibitor because oligomers are often more thermally stable than monomers. Using the homo-heptameric chaperonin, Hsp60, as a model, we screened ∼5000 diverse compounds in 384-well plates by DSF, revealing molecules that partially inhibited oligomerization. Because DSF does not require protein labeling or structural information, we propose that it could be a versatile way to uncover PPI inhibitors.

Design of artificial metalloenzymes with multiple inorganic elements: The more the merrier

A large fraction of metalloenzymes harbors multiple metal-centers that are electronically and/or functionally arranged within their proteinaceous environments. To explore the orchestration of inorganic and biochemical components and to develop bioinorganic catalysts and materials, we have described selected examples of artificial metalloproteins having multiple metallocofactors that were grouped depending on their initial protein scaffolds, the nature of introduced inorganic moieties, and the method used to propagate the number of metal ions within a protein. They demonstrated that diverse inorganic moieties can be selectively grafted and modulated in protein environments, providing a retrosynthetic bottom-up approach in the design of versatile and proficient biocatalysts and biomimetic model systems to explore fundamental questions in bioinorganic chemistry.

Emerging evidence for kindlin oligomerization and its role in regulating kindlin function

Integrin-mediated cell-extracellular matrix (ECM) interactions play crucial roles in a broad range of physiological and pathological processes. Kindlins are important positive regulators of integrin activation. The FERM-domain-containing kindlin family comprises three members, kindlin-1, kindlin-2 and kindlin-3 (also known as FERMT1, FERMT2 and FERMT3), which share high sequence similarity (identity >50%), as well as domain organization, but exhibit diverse tissue-specific expression patterns and cellular functions. Given the significance of kindlins, analysis of their atomic structures has been an attractive field for decades. Recently, the structures of kindlin and its β-integrin-bound form have been obtained, which greatly advance our understanding of the molecular functions that involve kindlins. In particular, emerging evidence indicates that oligomerization of kindlins might affect their integrin binding and focal adhesion localization, positively or negatively. In this Review, we presented an update on the recent progress of obtaining kindlin structures, and discuss the implication for integrin activation based on kindlin oligomerization, as well as the possible regulation of this process.

Dimerization of PHGDH via the catalytic unit is essential for its enzymatic function

Human D-3-phosphoglycerate dehydrogenase (PHGDH), a key enzyme in de novo serine biosynthesis, is amplified in various cancers and serves as a potential target for anticancer drug development. To facilitate this process, more information is needed on the basic biochemistry of this enzyme. For example, PHGDH was found to form tetramers in solution and the structure of its catalytic unit (sPHGDH) was solved as a dimer. However, how the oligomeric states affect PHGDH enzyme activity remains elusive. We studied the dependence of PHGDH enzymatic activity on its oligomeric states. We found that sPHGDH forms a mixture of monomers and dimers in solution with a dimer dissociation constant of ∼0.58 μM, with the enzyme activity depending on the dimer content. We computationally identified hotspot residues at the sPHGDH dimer interface. Single-point mutants at these sites disrupt dimer formation and abolish enzyme activity. Molecular dynamics simulations showed that dimer formation facilitates substrate binding and maintains the correct conformation required for enzyme catalysis. We further showed that the full-length PHGDH exists as a dynamic mixture of monomers, dimers, and tetramers in solution with enzyme concentration-dependent activity. Mutations that can completely disrupt the sPHGDH dimer show different abilities to interrupt the full-length PHGDH tetramer. Among them, E108A and I121A can also disrupt the oligomeric structures of the full-length PHGDH and abolish its enzyme activity. Our study indicates that disrupting the oligomeric structure of PHGDH serves as a novel strategy for PHGDH drug design and the hotspot residues identified can guide the design process.

Structural Characterization of Receptor-Receptor Interactions in the Allosteric Modulation of G Protein-Coupled Receptor (GPCR) Dimers

G protein-coupled receptor (GPCR) oligomerization, while contentious, continues to attract the attention of researchers. Numerous experimental investigations have validated the presence of GPCR dimers, and the relevance of dimerization in the effectuation of physiological functions intensifies the attractiveness of this concept as a potential therapeutic target. GPCRs, as a single entity, have been the main source of scrutiny for drug design objectives for multiple diseases such as cancer, inflammation, cardiac, and respiratory diseases. The existence of dimers broadens the research scope of GPCR functions, revealing new signaling pathways that can be targeted for disease pathogenesis that have not previously been reported when GPCRs were only viewed in their monomeric form. This review will highlight several aspects of GPCR dimerization, which include a summary of the structural elucidation of the allosteric modulation of class C GPCR activation offered through recent solutions to the three-dimensional, full-length structures of metabotropic glutamate receptor and γ-aminobutyric acid B receptor as well as the role of dimerization in the modification of GPCR function and allostery. With the growing influence of computational methods in the study of GPCRs, we will also be reviewing recent computational tools that have been utilized to map protein-protein interactions (PPI).

Multimerization variants as potential drivers of neofunctionalization

Whole-genome duplications are common during evolution, creating genetic redundancy that can enable cellular innovations. Novel protein-protein interactions provide a route to diversified gene functions, but, at present, there is limited proteome-scale knowledge on the extent to which variability in protein complex formation drives neofunctionalization. Here, we used protein correlation profiling to test for variability in apparent mass among thousands of orthologous proteins isolated from diverse species and cell types. Variants in protein complex size were unexpectedly common, in some cases appearing after relatively recent whole-genome duplications or an allopolyploidy event. In other instances, variants such as those in the carbonic anhydrase orthologous group reflected the neofunctionalization of ancient paralogs that have been preserved in extant species. Our results demonstrate that homo- and heteromer formation have the potential to drive neofunctionalization in diverse classes of enzymes, signaling, and structural proteins.

Protein phosphorylation regulates maize endosperm starch synthase IIa activity and protein-protein interactions

Starch synthesis is an elaborate process employing several isoforms of starch synthases (SSs), starch branching enzymes (SBEs) and debranching enzymes (DBEs). In cereals, some starch biosynthetic enzymes can form heteromeric complexes whose assembly is controlled by protein phosphorylation. Previous studies suggested that SSIIa forms a trimeric complex with SBEIIb, SSI, in which SBEIIb is phosphorylated. This study investigates the post-translational modification of SSIIa, and its interactions with SSI and SBEIIb in maize amyloplast stroma. SSIIa, immunopurified and shown to be free from other soluble starch synthases, was shown to be readily phosphorylated, affecting V_max but with minor effects on substrate K_d and K_m values, resulting in a 12-fold increase in activity compared with the dephosphorylated enzyme. This ATP-dependent stimulation of activity was associated with interaction with SBEIIb, suggesting that the availability of glucan branching limits SSIIa and is enhanced by physical interaction of the two enzymes. Immunoblotting of maize amyloplast extracts following non-denaturing polyacrylamide gel electrophoresis identified multiple bands of SSIIa, the electrophoretic mobilities of which were markedly altered by conditions that affected protein phosphorylation, including protein kinase inhibitors. Separation of heteromeric enzyme complexes by GPC, following alteration of protein phosphorylation states, indicated that such complexes are stable and may partition into larger and smaller complexes. The results suggest a dual role for protein phosphorylation in promoting association and dissociation of SSIIa-containing heteromeric enzyme complexes in the maize amyloplast stroma, providing new insights into the regulation of starch biosynthesis in plants.

A conserved cysteine-based redox mechanism sustains TFEB/HLH-30 activity under persistent stress

Mammalian TFEB and TFE3, as well as their ortholog in Caenorhabditis elegans HLH-30, play an important role in mediating cellular response to a variety of stress conditions, including nutrient deprivation, oxidative stress, and pathogen infection. In this study, we identify a novel mechanism of TFEB/HLH-30 regulation through a cysteine-mediated redox switch. Under stress conditions, TFEB-C212 undergoes oxidation, allowing the formation of intermolecular disulfide bonds that result in TFEB oligomerization. TFEB oligomers display increased resistance to mTORC1-mediated inactivation and are more stable under prolonged stress conditions. Mutation of the only cysteine residue present in HLH-30 (C284) significantly reduced its activity, resulting in developmental defects and increased pathogen susceptibility in worms. Therefore, cysteine oxidation represents a new type of TFEB post-translational modification that functions as a molecular switch to link changes in redox balance with expression of TFEB/HLH-30 target genes.

Supra-Molecular Assemblies of ORAI1 at Rest Precede Local Accumulation into Puncta after Activation

The Ca²⁺ selective channel ORAI1 and endoplasmic reticulum (ER)-resident STIM proteins form the core of the channel complex mediating store operated Ca²⁺ entry (SOCE). Using liquid phase electron microscopy (LPEM), the distribution of ORAI1 proteins was examined at rest and after SOCE-activation at nanoscale resolution. The analysis of over seven hundred thousand ORAI1 positions revealed a number of ORAI1 channels had formed STIM-independent distinct supra-molecular clusters. Upon SOCE activation and in the presence of STIM proteins, a fraction of ORAI1 assembled in micron-sized two-dimensional structures, such as the known puncta at the ER plasma membrane contact zones, but also in divergent structures such as strands, and ring-like shapes. Our results thus question the hypothesis that stochastically migrating single ORAI1 channels are trapped at regions containing activated STIM, and we propose instead that supra-molecular ORAI1 clusters fulfill an amplifying function for creating dense ORAI1 accumulations upon SOCE-activation.

Estimating Change in Foldability Due to Multipoint Deletions in Protein Structures

Insertions/deletions of amino acids in the protein backbone potentially result in altered structural/functional specifications. They can either contribute positively to the evolutionary process or can result in disease conditions. Despite being the second most prevalent form of protein modification, there are no databases or computational frameworks that delineate harmful multipoint deletions (MPD) from beneficial ones. We introduce a positive unlabeled learning-based prediction framework (PROFOUND) that utilizes fold-level attributes, environment-specific properties, and deletion site-specific properties to predict the change in foldability arising from such MPDs, both in the non-loop and loop regions of protein structures. In the absence of any protein structure dataset to study MPDs, we introduce a dataset with 153 MPD instances that lead to native-like folded structures and 7650 unlabeled MPD instances whose effect on the foldability of the corresponding proteins is unknown. PROFOUND on 10-fold cross-validation on our newly introduced dataset reports a recall of 82.2% (86.6%) and a fall out rate (FR) of 14.2% (20.6%), corresponding to MPDs in the protein loop (non-loop) region. The low FR suggests that the foldability in proteins subject to MPDs is not random and necessitates unique specifications of the deleted region. In addition, we find that additional evolutionary attributes contribute to higher recall and lower FR. The first of a kind foldability prediction system owing to MPD instances and the newly introduced dataset will potentially aid in novel protein engineering endeavors.

Concentration- and pH-Dependent Oligomerization of the Thrombin-Derived C-Terminal Peptide TCP-25

Peptide oligomerization dynamics affects peptide structure, activity, and pharmacodynamic properties. The thrombin C-terminal peptide, TCP-25 (GKYGFYTHVFRLKKWIQKVIDQFGE), is currently in preclinical development for improved wound healing and infection prevention. It exhibits turbidity when formulated at pH 7.4, particularly at concentrations of 0.3 mM or more. We used biochemical and biophysical approaches to explore whether the peptide self-associates and forms oligomers. The peptide showed a dose-dependent increase in turbidity as well as α-helical structure at pH 7.4, a phenomenon not observed at pH 5.0. By analyzing the intrinsic tryptophan fluorescence, we demonstrate that TCP-25 is more stable at high concentrations (0.3 mM) when exposed to high temperatures or a high concentration of denaturant agents, which is compatible with oligomer formation. The denaturation process was reversible above 100 µM of peptide. Dynamic light scattering demonstrated that TCP-25 oligomerization is sensitive to changes in pH, time, and temperature. Computational modeling with an active 18-mer region of TCP-25 showed that the peptide can form pH-dependent higher-order end-to-end oligomers and micelle-like structures, which is in agreement with the experimental data. Thus, TCP-25 exhibits pH- and temperature-dependent dynamic changes involving helical induction and reversible oligomerization, which explains the observed turbidity of the pharmacologically developed formulation.

Human Cellular Retinol Binding Protein II Forms a Domain-Swapped Trimer Representing a Novel Fold and a New Template for Protein Engineering

Domain-swapping is a mechanism for evolving new protein structure from extant scaffolds, and has been an efficient protein-engineering strategy for tailoring functional diversity. However, domain swapping can only be exploited if it can be controlled, especially in cases where various folds can coexist. Herein, we describe the structure of a domain-swapped trimer of the iLBP family member hCRBPII, and suggest a mechanism for domain-swapped trimerization. It is further shown that domain-swapped trimerization can be favored by strategic installation of a disulfide bond, thus demonstrating a strategy for fold control. We further show the domain-swapped trimer to be a useful protein design template by installing a high-affinity metal binding site through the introduction of a single mutation, taking advantage of its threefold symmetry. Together, these studies show how nature can promote oligomerization, stabilize a specific oligomer, and generate new function with minimal changes to the protein sequence.

Polymeric IgA paraprotein on agarose gel electrophoresis immunofixation identifies a unique subset of IgA myeloma patients

OBJECTIVES

On agarose gel electrophoresis, IgA paraprotein has a tendency to polymerise and form multiple bands on immunofixation. We decided to investigate if there are any differences in clinical parameters with monomeric vs. polymeric IgA paraprotein multiple myeloma (MM) patients.

METHODS

During an 18-month retrospective and prospective cross-sectional audit review period at one Australian laboratory, we identified 92 IgA MM patients that were divided up according to monomeric or polymeric IgA paraproteins based on their appearance on IFE. Medical and pathology records were reviewed for demographic details, and laboratory data to examine for end-organ manifestations of MM.

RESULTS

After correcting for age, polymeric IgA MM patients had a greater degree of proteinuria and hence, higher incidence of hypogammaglobulinaemia. The patients tended to have a higher quantity of paraprotein as well. No difference in mortality was seen.

CONCLUSIONS

Our study is the first to stratify IgA MM patients according to the laboratory appearance of their paraprotein and may hold important prognostic and predictive clues for these patients.