DeepSub: Utilizing Deep Learning for Predicting the Number of Subunits in Homo-Oligomeric Protein Complexes.

The molecular weight (MW) of an enzyme is a critical parameter in enzyme-constrained models (ecModels). It is determined by two factors: the presence of subunits and the abundance of each subunit. Although the number of subunits (NS) can potentially be obtained from UniProt, this information is not readily available for most proteins. In this study, we addressed this gap by extracting and curating subunit information from the UniProt database to establish a robust benchmark dataset. Subsequently, we propose a novel model named DeepSub, which leverages the protein language model and Bi-directional Gated Recurrent Unit (GRU), to predict NS in homo-oligomers solely based on protein sequences. DeepSub demonstrates remarkable accuracy, achieving an accuracy rate as high as 0.967, surpassing the performance of QUEEN. To validate the effectiveness of DeepSub, we performed predictions for protein homo-oligomers that have been reported in the literature but are not documented in the UniProt database. Examples include homoserine dehydrogenase from Corynebacterium glutamicum, Matrilin-4 from Mus musculus and Homo sapiens, and the Multimerins protein family from M. musculus and H. sapiens. The predicted results align closely with the reported findings in the literature, underscoring the reliability and utility of DeepSub.

Analysis of Surface Expression of NMDAR Subunits in Primary Hippocampal Neurons.

The expression and activity of ionotropic glutamate receptors control signal transduction at the excitatory synapses in the CNS. The NMDAR comprises two obligatory GluN1 subunits and two GluN2 or GluN3 subunits in different combinations. Each GluN subunit consists of four domains: the extracellular amino-terminal and agonist-binding domains, the transmembrane domain, and the intracellular C-terminal domain (CTD). The CTD interaction with various classes of intracellular proteins is critical for trafficking and synaptic localization of NMDARs. Amino acid mutations or the inclusion of premature stop codons in the CTD could contribute to the emergence of neurodevelopmental and neuropsychiatric disorders. Here, we describe the method of preparing primary hippocampal neurons and lentiviral particles expressing GluN subunits that can be used as a model to study cell surface expression and synaptic localization of NMDARs. We also show a simple method of fluorescence immunostaining of eGFP-tagged GluN2 subunits and subsequent microscopy technique and image analysis to study the effects of disease-associated mutations in the CTDs of GluN2A and GluN2B subunits.

Functional Analysis of NMDAR Subunit Components in Postsynaptic Currents of Identified Cells and Synapses in Brain Slices.

Epilepsy is one of the most represented neurological diseases worldwide. However, in many cases, the precise molecular mechanisms of epileptogenesis and ictiogenesis are unknown. Because of their important role in synaptic function and neuronal excitability, NMDA receptors are implicated in various epileptogenic mechanisms. Most of these are subunit specific and require a precise analysis of the subunit composition of the NMDARs implicated. Here, we describe an express electrophysiological method to analyze the contribution of NMDAR subunits to spontaneous postsynaptic activity in identified cells in brain slices using patch clamp whole cell recordings.

Selective Cell-Surface Expression of Triheteromeric NMDA Receptors.

NMDA-type ionotropic glutamate receptors are critically involved in many brain functions and are implicated in a variety of brain disorders. Seven NMDA receptor subunits exist (GluN1, GluN2A-D, and GluN3A-B) that assemble into tetrameric receptor subtypes with distinct functional properties and physiological roles. The majority NMDA receptors are composed of two GluN1 and two GluN2 subunits, which can assemble into four diheteromeric receptors subtypes composed of GluN1 and one type of GluN2 subunit (e.g., GluN1/2A), and presumably also six triheteromeric receptor subtypes composed of GluN1 and two different GluN2 subunits (e.g., GluN1/2A/2B). Furthermore, the GluN1 subunit exists as eight splice variants (e.g., GluN1-1a and GluN1-1b isoforms), and two different GluN1 isoforms can co-assemble to also form triheteromeric NMDA receptors (e.g., GluN1-1a/1b/2A). Here, we describe a method to faithfully express triheteromeric NMDA receptors in heterologous expression systems by controlling the identity of two of the four subunits. This method overcomes the problem that co-expression of three different NMDA receptor subunits generates two distinct diheteromeric receptor subtypes as well as one triheteromeric receptor subtype, thereby confounding studies that require a homogenous population of triheteromeric NMDA receptors. The method has been applied to selectively express recombinant triheteromeric GluN1/2A/2B, GluN1/2A/2C, GluN1/2B/2D, GluN1-1a/GluN1-1b/2A, GluN1-1a/GluN1-1b/2B receptors with negligible co-expression of the respective diheteromeric receptor subtypes. This method therefore enables quantitative evaluation of functional and pharmacological properties of triheteromeric NMDA receptors, some of which are abundant NMDA receptor subtypes in the adult brain.

Characterizing Functional Contributions of Specific GluN2 Subunits to Individual Postsynaptic NMDAR Responses Using Biophysical Parameters.

The NMDAR is a heterotetramer composed of two GluN1 subunits and two GluN2 and/or GluN3 subunits, with the GluN2 subunits exhibiting significant diversity in their structure and function. Recent studies have highlighted the importance of characterizing the specific roles of each GluN2 subunit across central nervous system regions and developmental stages, as well as their unique contributions to NMDAR-mediated signaling and plasticity. Understanding the distinct functions of GluN2 subunits is critical for the development of targeted therapeutic strategies for NMDAR-related disorders. However, measuring the functional contribution of individual GluN2 subtypes in ex vivo slices is challenging. Conventionally, pharmacological or genetic approaches are used, but, in many cases, this is not possible or is restricted to population-level NMDAR responses. Here, we describe a technique for using biophysical properties of miniature synaptic NMDAR responses as a proxy to measure the functional contribution of specific GluN2-NMDAR subunits to individual synapses within a neuron.

CORVET-specific subunit levels determine the balance between HOPS/CORVET endosomal tethering complexes.

The closely related endolysosomal tethering complexes HOPS and CORVET play pivotal roles in the homo- and heterotypic fusion of early and late endosomes, respectively, and HOPS also mediates the fusion of lysosomes with incoming vesicles including late endosomes and autophagosomes. These heterohexameric complexes share their four core subunits that assemble with additional two, complex-specific subunits. These features and the similar structure of the complexes could allow the formation of hybrid complexes, and the complex specific subunits may compete for binding to the core. Indeed, our biochemical analyses revealed the overlap of binding sites for HOPS-specific VPS41 and CORVET-specific VPS8 on the shared core subunit VPS18. We found that the overexpression of CORVET-specific VPS8 or Tgfbrap1 decreased the amount of core proteins VPS11 and VPS18 that are assembled with HOPS-specific subunits VPS41 or VPS39, indicating reduced amount of assembled HOPS. In line with this, we observed the elevation of both lipidated, autophagosome-associated LC3 protein and the autophagic cargo p62 in these cells, suggesting impaired autophagosome-lysosome fusion. In contrast, overexpression of HOPS-specific VPS39 or VPS41 did not affect the level of assembled CORVET or autophagy. VPS8 or Tgfbrap1 overexpression also induced Cathepsin D accumulation, suggesting that HOPS-dependent biosynthetic delivery of lysosomal hydrolases is perturbed, too. These indicate that CORVET-specific subunit levels fine-tune HOPS assembly and activity in vivo.

Cryo-EM Structure of Porphyromonas gingivalis RNA Polymerase.

Porphyromonas gingivalis, an anaerobic CFB (Cytophaga, Fusobacterium, and Bacteroides) group bacterium, is the keystone pathogen of periodontitis and has been implicated in various systemic diseases. Increased antibiotic resistance and lack of effective antibiotics necessitate a search for new intervention strategies. Here we report a 3.5 Å resolution cryo-EM structure of P. gingivalis RNA polymerase (RNAP). The structure displays new structural features in its ω subunit and multiple domains in ß and ß' subunits, which differ from their counterparts in other bacterial RNAPs. Superimpositions with E. coli RNAP holoenzyme and initiation complex further suggest that its ω subunit may contact the σ4 domain, thereby possibly contributing to the assembly and stabilization of initiation complexes. In addition to revealing the unique features of P. gingivalis RNAP, our work offers a framework for future studies of transcription regulation in this important pathogen, as well as for structure-based drug development.

How an intrinsically disordered regulatory subunit assembles a PP1:eIF2 complex.

Fatalska et al.1 use an interdisciplinary strategy to elucidate how an intrinsically disordered regulatory subunit of protein phosphatase 1 binds trimeric eIF2 and positions the phosphatase-substrate complex for dephosphorylation. As validation, they show that a disease mutation abolishes the interaction.

Structural basis of Integrator-dependent RNA polymerase II termination.

The Integrator complex can terminate RNA polymerase II (Pol II) in the promoter-proximal region of genes. Previous work has shed light on how Integrator binds to the paused elongation complex consisting of Pol II, the DRB sensitivity-inducing factor (DSIF) and the negative elongation factor (NELF) and how it cleaves the nascent RNA transcript1, but has not explained how Integrator removes Pol II from the DNA template. Here we present three cryo-electron microscopy structures of the complete Integrator-PP2A complex in different functional states. The structure of the pre-termination complex reveals a previously unresolved, scorpion-tail-shaped INTS10-INTS13-INTS14-INTS15 module that may use its 'sting' to open the DSIF DNA clamp and facilitate termination. The structure of the post-termination complex shows that the previously unresolved subunit INTS3 and associated sensor of single-stranded DNA complex (SOSS) factors prevent Pol II rebinding to Integrator after termination. The structure of the free Integrator-PP2A complex in an inactive closed conformation2 reveals that INTS6 blocks the PP2A phosphatase active site. These results lead to a model for how Integrator terminates Pol II transcription in three steps that involve major rearrangements.

Architecture and activation of human muscle phosphorylase kinase.

The study of phosphorylase kinase (PhK)-regulated glycogen metabolism has contributed to the fundamental understanding of protein phosphorylation; however, the molecular mechanism of PhK remains poorly understood. Here we present the high-resolution cryo-electron microscopy structures of human muscle PhK. The 1.3-megadalton PhK α4ß4γ4δ4 hexadecamer consists of a tetramer of tetramer, wherein four αßγδ modules are connected by the central ß4 scaffold. The α- and ß-subunits possess glucoamylase-like domains, but exhibit no detectable enzyme activities. The α-subunit serves as a bridge between the ß-subunit and the γδ subcomplex, and facilitates the γ-subunit to adopt an autoinhibited state. Ca2+-free calmodulin (δ-subunit) binds to the γ-subunit in a compact conformation. Upon binding of Ca2+, a conformational change occurs, allowing for the de-inhibition of the γ-subunit through a spring-loaded mechanism. We also reveal an ADP-binding pocket in the ß-subunit, which plays a role in allosterically enhancing PhK activity. These results provide molecular insights of this important kinase complex.

Structure and assembly of a bacterial gasdermin pore.

In response to pathogen infection, gasdermin (GSDM) proteins form membrane pores that induce a host cell death process called pyroptosis1-3. Studies of human and mouse GSDM pores have revealed the functions and architectures of assemblies comprising 24 to 33 protomers4-9, but the mechanism and evolutionary origin of membrane targeting and GSDM pore formation remain unknown. Here we determine a structure of a bacterial GSDM (bGSDM) pore and define a conserved mechanism of pore assembly. Engineering a panel of bGSDMs for site-specific proteolytic activation, we demonstrate that diverse bGSDMs form distinct pore sizes that range from smaller mammalian-like assemblies to exceptionally large pores containing more than 50 protomers. We determine a cryo-electron microscopy structure of a Vitiosangium bGSDM in an active 'slinky'-like oligomeric conformation and analyse bGSDM pores in a native lipid environment to create an atomic-level model of a full 52-mer bGSDM pore. Combining our structural analysis with molecular dynamics simulations and cellular assays, our results support a stepwise model of GSDM pore assembly and suggest that a covalently bound palmitoyl can leave a hydrophobic sheath and insert into the membrane before formation of the membrane-spanning ß-strand regions. These results reveal the diversity of GSDM pores found in nature and explain the function of an ancient post-translational modification in enabling programmed host cell death.

Exploring the structural assembly of rice ADP-glucose pyrophosphorylase subunits using MD simulation.

ADP-glucose pyrophosphorylase plays a pivotal role as an allosteric enzyme, essential for starch biosynthesis in plants. The higher plant AGPase comparises of a pair of large and a pair of small subunits to form a heterotetrameric complex. Growing evidence indicates that each subunit plays a distinct role in regulating the underlying mechanism of starch biosynthesis. In the rice genome, there are four large subunit genes (OsL1-L4) and three small subunit genes (OsS1, OsS2a, and OsS2b). While the structural assembly of cytosolic rice AGPase subunits (OsL2:OsS2b) has been elucidated, there is currently no such documented research available for plastidial rice AGPases (OsL1:OsS1). In this study, we employed protein modeling and MD simulation approaches to gain insights into the structural association of plastidial rice AGPase subunits. Our results demonstrate that the heterotetrameric association of OsL1:OsS1 is very similar to that of cytosolic OsL2:OsS2b and potato AGPase heterotetramer (StLS:StSS). Moreover, the yeast-two-hybrid results on OsL1:OsS1, which resemble StLS:StSS, suggest a differential protein assembly for OsL2:OsS2b. Thus, the regulatory and catalytic mechanisms for plastidial AGPases (OsL1:OsS1) could be different in rice culm and developing endosperm compared to those of OsL2:OsS2b, which are predominantly found in rice endosperm.

Cryo-EM structures of RAD51 assembled on nucleosomes containing a DSB site.

RAD51 is the central eukaryotic recombinase required for meiotic recombination and mitotic repair of double-strand DNA breaks (DSBs)1,2. However, the mechanism by which RAD51 functions at DSB sites in chromatin has remained elusive. Here we report the cryo-electron microscopy structures of human RAD51-nucleosome complexes, in which RAD51 forms ring and filament conformations. In the ring forms, the N-terminal lobe domains (NLDs) of RAD51 protomers are aligned on the outside of the RAD51 ring, and directly bind to the nucleosomal DNA. The nucleosomal linker DNA that contains the DSB site is recognized by the L1 and L2 loops-active centres that face the central hole of the RAD51 ring. In the filament form, the nucleosomal DNA is peeled by the RAD51 filament extension, and the NLDs of RAD51 protomers proximal to the nucleosome bind to the remaining nucleosomal DNA and histones. Mutations that affect nucleosome-binding residues of the RAD51 NLD decrease nucleosome binding, but barely affect DNA binding in vitro. Consistently, yeast Rad51 mutants with the corresponding mutations are substantially defective in DNA repair in vivo. These results reveal an unexpected function of the RAD51 NLD, and explain the mechanism by which RAD51 associates with nucleosomes, recognizes DSBs and forms the active filament in chromatin.

Cryo-EM structure of the Slo1 potassium channel with the auxiliary γ1 subunit suggests a mechanism for depolarization-independent activation.

Mammalian Ca2+-dependent Slo K+ channels can stably associate with auxiliary γ subunits which fundamentally alter their behavior. By a so far unknown mechanism, the four γ subunits reduce the need for voltage-dependent activation and, thereby, allow Slo to open independently of an action potential. Here, using cryo-EM, we reveal how the transmembrane helix of γ1/LRRC26 binds and presumably stabilizes the activated voltage-sensor domain of Slo1. The activation is further enhanced by an intracellular polybasic stretch which locally changes the charge gradient across the membrane. Our data provide a possible explanation for Slo1 regulation by the four γ subunits and also their different activation efficiencies. This suggests a novel activation mechanism of voltage-gated ion channels by auxiliary subunits.

The NADPH oxidase 2 subunit p47phox binds to the WAVE regulatory complex and p22phox in a mutually exclusive manner.

The actin cytoskeleton and reactive oxygen species (ROS) both play crucial roles in various cellular processes. Previous research indicated a direct interaction between two key components of these systems: the WAVE1 subunit of the WAVE regulatory complex (WRC), which promotes actin polymerization and the p47phox subunit of the NADPH oxidase 2 complex (NOX2), which produces ROS. Here, using carefully characterized recombinant proteins, we find that activated p47phox uses its dual Src homology 3 domains to bind to multiple regions within the WAVE1 and Abi2 subunits of the WRC, without altering WRC's activity in promoting Arp2/3-mediated actin polymerization. Notably, contrary to previous findings, p47phox uses the same binding pocket to interact with both the WRC and the p22phox subunit of NOX2, albeit in a mutually exclusive manner. This observation suggests that when activated, p47phox may separately participate in two distinct processes: assembling into NOX2 to promote ROS production and engaging with WRC to regulate the actin cytoskeleton.

Structure of human phagocyte NADPH oxidase in the activated state.

Phagocyte NADPH oxidase, a protein complex with a core made up of NOX2 and p22 subunits, is responsible for transferring electrons from intracellular NADPH to extracellular oxygen1. This process generates superoxide anions that are vital for killing pathogens1. The activation of phagocyte NADPH oxidase requires membrane translocation and the binding of several cytosolic factors2. However, the exact mechanism by which cytosolic factors bind to and activate NOX2 is not well understood. Here we present the structure of the human NOX2-p22 complex activated by fragments of three cytosolic factors: p47, p67 and Rac1. The structure reveals that the p67-Rac1 complex clamps onto the dehydrogenase domain of NOX2 and induces its contraction, which stabilizes the binding of NADPH and results in a reduction of the distance between the NADPH-binding domain and the flavin adenine dinucleotide (FAD)-binding domain. Furthermore, the dehydrogenase domain docks onto the bottom of the transmembrane domain of NOX2, which reduces the distance between FAD and the inner haem. These structural rearrangements might facilitate the efficient transfer of electrons between the redox centres in NOX2 and lead to the activation of phagocyte NADPH oxidase.

Structural basis for intra- and intermolecular interactions on RAD9 subunit of 9-1-1 checkpoint clamp implies functional 9-1-1 regulation by RHINO.

Eukaryotic DNA clamp is a trimeric protein featuring a toroidal ring structure that binds DNA on the inside of the ring and multiple proteins involved in DNA transactions on the outside. Eukaryotes have two types of DNA clamps: the replication clamp PCNA and the checkpoint clamp RAD9-RAD1-HUS1 (9-1-1). 9-1-1 activates the ATR-CHK1 pathway in DNA damage checkpoint, regulating cell cycle progression. Structure of 9-1-1 consists of two moieties: a hetero-trimeric ring formed by PCNA-like domains of three subunits and an intrinsically disordered C-terminal region of the RAD9 subunit, called RAD9 C-tail. The RAD9 C-tail interacts with the 9-1-1 ring and disrupts the interaction between 9-1-1 and DNA, suggesting a negative regulatory role for this intramolecular interaction. In contrast, RHINO, a 9-1-1 binding protein, interacts with both RAD1 and RAD9 subunits, positively regulating checkpoint activation by 9-1-1. This study presents a biochemical and structural analysis of intra- and inter-molecular interactions on the 9-1-1 ring. Biochemical analysis indicates that RAD9 C-tail binds to the hydrophobic pocket on the PCNA-like domain of RAD9, implying that the pocket is involved in multiple protein-protein interactions. The crystal structure of the 9-1-1 ring in complex with a RHINO peptide reveals that RHINO binds to the hydrophobic pocket of RAD9, shedding light on the RAD9-binding motif. Additionally, the study proposes a structural model of the 9-1-1-RHINO quaternary complex. Together, these findings provide functional insights into the intra- and inter-molecular interactions on the front side of RAD9, elucidating the roles of RAD9 C-tail and RHINO in checkpoint activation.

Chemical shift assignments of the ACID domain of MED25, a subunit of the mediator complex in Arabidopsis thaliana.

Mediator complex is a key component that bridges various transcription activators and RNA polymerase during eukaryotic transcription initiation. The Arabidopsis thaliana Med25 (aMed25), a subunit of the Mediator complex, plays important roles in regulating hormone signaling, biotic and abiotic stress responses and plant development by interacting with a variety of transcription factors through its activator-interacting domain (ACID). However, the recognition mechanism of aMed25-ACID for various transcription factors remains unknown. Here, we report the nearly complete 1H, 13C, and 15N backbone and side chain resonance assignments of aMED25-ACID (residues 551-681). TALOS-N analysis revealed that aMED25-ACID structure is comprised of three α-helices and seven ß-strands, which lacks the C-terminal α-helix existing in the human MED25-ACID. This study lays a foundation for further research on the structure-function relationship of aMED25-ACID.

Coregulation of NDC80 Complex Subunits Determines the Fidelity of the Spindle-Assembly Checkpoint and Mitosis.

NDC80 complex (NDC80C) is composed of four subunits (SPC24, SPC25, NDC80, and NUF2) and is vital for kinetochore-microtubule (KT-MT) attachment during mitosis. Paradoxically, NDC80C also functions in the activation of the spindle-assembly checkpoint (SAC). This raises an interesting question regarding how mitosis is regulated when NDC80C levels are compromised. Using a degron-mediated depletion system, we found that acute silencing of SPC24 triggered a transient mitotic arrest followed by mitotic slippage. SPC24-deficient cells were unable to sustain SAC activation despite the loss of KT-MT interaction. Intriguingly, our results revealed that other subunits of the NDC80C were co-downregulated with SPC24 at a posttranslational level. Silencing any individual subunit of NDC80C likewise reduced the expression of the entire complex. We found that the SPC24-SPC25 and NDC80-NUF2 subcomplexes could be individually stabilized using ectopically expressed subunits. The synergism of SPC24 downregulation with drugs that promote either mitotic arrest or mitotic slippage further underscored the dual roles of NDC80C in KT-MT interaction and SAC maintenance. The tight coordinated regulation of NDC80C subunits suggests that targeting individual subunits could disrupt mitotic progression and provide new avenues for therapeutic intervention. IMPLICATIONS: These results highlight the tight coordinated regulation of NDC80C subunits and their potential as targets for antimitotic therapies.

Coevolution of RNA and protein subunits in RNase P and RNase MRP, two RNA processing enzymes.

RNase P and RNase mitochondrial RNA processing (MRP) are ribonucleoproteins (RNPs) that consist of a catalytic RNA and a varying number of protein cofactors. RNase P is responsible for precursor tRNA maturation in all three domains of life, while RNase MRP, exclusive to eukaryotes, primarily functions in rRNA biogenesis. While eukaryotic RNase P is associated with more protein cofactors and has an RNA subunit with fewer auxiliary structural elements compared to its bacterial cousin, the double-anchor precursor tRNA recognition mechanism has remarkably been preserved during evolution. RNase MRP shares evolutionary and structural similarities with RNase P, preserving the catalytic core within the RNA moiety inherited from their common ancestor. By incorporating new protein cofactors and RNA elements, RNase MRP has established itself as a distinct RNP capable of processing ssRNA substrates. The structural information on RNase P and MRP helps build an evolutionary trajectory, depicting how emerging protein cofactors harmonize with the evolution of RNA to shape different functions for RNase P and MRP. Here, we outline the structural and functional relationship between RNase P and MRP to illustrate the coevolution of RNA and protein cofactors, a key driver for the extant, diverse RNP world.