Hexamerization: explaining the original sin of IgG-mediated complement activation in acute lung injury.

Although antibody-mediated lung damage is a major factor in transfusion-related acute lung injury (ALI), autoimmune lung disease (for example, coatomer subunit α [COPA] syndrome), and primary graft dysfunction following lung transplantation, the mechanism by which antigen-antibody complexes activate complement to induce lung damage remains unclear. In this issue of the JCI, Cleary and colleagues utilized several approaches to demonstrate that IgG forms hexamers with MHC class I alloantibodies. This hexamerization served as a key pathophysiological mechanism in alloimmune lung injury models and was mediated through the classical pathway of complement activation. Additionally, the authors provided avenues for exploring therapeutics for this currently hard-to-treat clinical entity that has several etiologies but a potentially focused mechanism.

Asymmetric allostery in estrogen receptor-α homodimers drives responses to the ensemble of estrogens in the hormonal milieu.

The estrogen receptor-α (ER) is thought to function only as a homodimer but responds to a variety of environmental, metazoan, and therapeutic estrogens at subsaturating doses, supporting binding mixtures of ligands as well as dimers that are only partially occupied. Here, we present a series of flexible ER ligands that bind to receptor dimers with individual ligand poses favoring distinct receptor conformations-receptor conformational heterodimers-mimicking the binding of two different ligands. Molecular dynamics simulations showed that the pairs of different ligand poses changed the correlated motion across the dimer interface to generate asymmetric communication between the dimer interface, the ligands, and the surface binding sites for epigenetic regulatory proteins. By examining the binding of the same ligand in crystal structures of ER in the agonist vs. antagonist conformers, we also showed that these allosteric signals are bidirectional. The receptor conformer can drive different ligand binding modes to support agonist vs. antagonist activity profiles, a revision of ligand binding theory that has focused on unidirectional signaling from the ligand to the coregulator binding site. We also observed differences in the allosteric signals between ligand and coregulator binding sites in the monomeric vs. dimeric receptor, and when bound by two different ligands, states that are physiologically relevant. Thus, ER conformational heterodimers integrate two different ligand-regulated activity profiles, representing different modes for ligand-dependent regulation of ER activity.

The step-by-step assembly mechanism of secreted flavivirus NS1 tetramer and hexamer captured at atomic resolution.

Flaviviruses encode a conserved, membrane-associated nonstructural protein 1 (NS1) with replication and immune evasion functions. The current knowledge of secreted NS1 (sNS1) oligomers is based on several low-resolution structures, thus hindering the development of drugs and vaccines against flaviviruses. Here, we revealed that recombinant sNS1 from flaviviruses exists in a dynamic equilibrium of dimer-tetramer-hexamer states. Two DENV4 hexameric NS1 structures and several tetrameric NS1 structures from multiple flaviviruses were solved at atomic resolution by cryo-EM. The stacking of the tetrameric NS1 and hexameric NS1 is facilitated by the hydrophobic ß-roll and connector domains. Additionally, a triacylglycerol molecule located within the central cavity may play a role in stabilizing the hexamer. Based on differentiated interactions between the dimeric NS1, two distinct hexamer models (head-to-head and side-to-side hexamer) and the step-by-step assembly mechanisms of NS1 dimer into hexamer were proposed. We believe that our study sheds light on the understanding of the NS1 oligomerization and contributes to NS1-based therapies.

Study of individual domains contributing to MALT1 dimerization in BCL10-independent and dependent assembly.

The CARMA-BCL10-MALT1 (CBM) signalosome functions as a pivotal supramolecular module, integrating diverse receptor-induced signaling pathways to regulate BCL10-dependent NF-kB activation in innate and adaptive immunity. Conversely, the API2-MALT1 fusion protein in t(11; 18)(q21; q21) MALT lymphoma constitutively induces BCL10-independent NF-kB activation. MALT1 dimer formation is indispensable for the requisite proteolytic activity and is critical for NF-kB activation regulation in both scenarios. However, the molecular assembly of MALT1 individual domains in CBM activation remains elusive. Here we report the crystal structure of the MALT1 death domain (DD) at a resolution of 2.1 Å, incorporating reconstructed residues in previously disordered loops 1 and 2. Additionally, we observe a conformational regulation element (CRE) regulating stem-helix formation in NLRPs pyrin (PYD) within the MALT1 DD structure. The structure reveals a stem-helix-mediated dimer further corroborated in solution. To elucidate how the BCL10 filament facilitates MALT1 dimerization, we reconstitute a BCL10-CARD-MALT1-DD-IG1-IG2 complex model. We propose a N+7 rule for BCL10-dependent MALT1 dimerization via the IG1-IG2 domain and for MALT1-dependent cleavage in trans. Biochemical data further indicates concentration-dependent dimerization of the MALT1 IG1-IG2 domain, facilitating MALT1 dimerization in BCL10-independent manner. Our findings provide a structural and biochemical foundation for understanding MALT1 dimeric mechanisms, shedding light on potential BCL10-independent MALT1 dimer formation and high-order BCL10-MALT1 assembly.

Rescuing tri-heteromeric NMDA receptor function: the potential of pregnenolone-sulfate in loss-of-function GRIN2B variants.

N-methyl-D-aspartate receptors (NMDARs emerging from GRIN genes) are tetrameric receptors that form diverse channel compositions in neurons, typically consisting of two GluN1 subunits combined with two GluN2(A-D) subunits. During prenatal stages, the predominant channels are di-heteromers with two GluN1 and two GluN2B subunits due to the high abundance of GluN2B subunits. Postnatally, the expression of GluN2A subunits increases, giving rise to additional subtypes, including GluN2A-containing di-heteromers and tri-heteromers with GluN1, GluN2A, and GluN2B subunits. The latter  emerge as the major receptor subtype at mature synapses in the hippocampus. Despite extensive research on purely di-heteromeric receptors containing two identical GRIN variants, the impact of a single variant on the function of other channel forms, notably tri-heteromers, is lagging. In this study, we systematically investigated the effects of two de novo GRIN2B variants (G689C and G689S) in pure, mixed di- and tri-heteromers. Our findings reveal that incorporating a single variant in mixed di-heteromers or tri-heteromers exerts a dominant negative effect on glutamate potency, although 'mixed' channels show improved potency compared to pure variant-containing di-heteromers. We show that a single variant within a receptor complex does not impair the response of all receptor subtypes to the positive allosteric modulator pregnenolone-sulfate (PS), whereas spermine completely fails to potentiate tri-heteromers containing GluN2A and -2B-subunits. We examined PS on primary cultured hippocampal neurons transfected with the variants, and observed a positive impact over current amplitudes and synaptic activity. Together, our study supports previous observations showing that mixed di-heteromers exhibit improved glutamate potency and extend these findings towards the exploration of the effect of Loss-of-Function variants over tri-heteromers. Notably, we provide an initial and crucial demonstration of the beneficial effects of GRIN2B-relevant potentiators on tri-heteromers. Our results underscore the significance of studying how different variants affect distinct receptor subtypes, as these effects cannot be inferred solely from observations made on pure di-heteromers. Overall, this study contributes to ongoing efforts to understand the pathophysiology of GRINopathies and provides insights into potential treatment strategies.

AlphaFold2 modeling and molecular dynamics simulations of an intrinsically disordered protein.

We use AlphaFold2 (AF2) to model the monomer and dimer structures of an intrinsically disordered protein (IDP), Nvjp-1, assisted by molecular dynamics (MD) simulations. We observe relatively rigid dimeric structures of Nvjp-1 when compared with the monomer structures. We suggest that protein conformations from multiple AF2 models and those from MD trajectories exhibit a coherent trend: the conformations of an IDP are deviated from each other and the conformations of a well-folded protein are consistent with each other. We use a residue-residue interaction network (RIN) derived from the contact map which show that the residue-residue interactions in Nvjp-1 are mainly transient; however, those in a well-folded protein are mainly persistent. Despite the variation in 3D shapes, we show that the AF2 models of both disordered and ordered proteins exhibit highly consistent profiles of the pLDDT (predicted local distance difference test) scores. These results indicate a potential protocol to justify the IDPs based on multiple AF2 models and MD simulations.

Cyanobacterial α-carboxysome carbonic anhydrase is allosterically regulated by the Rubisco substrate RuBP.

Cyanobacterial CO2 concentrating mechanisms (CCMs) sequester a globally consequential proportion of carbon into the biosphere. Proteinaceous microcompartments, called carboxysomes, play a critical role in CCM function, housing two enzymes to enhance CO2 fixation: carbonic anhydrase (CA) and Rubisco. Despite its importance, our current understanding of the carboxysomal CAs found in α-cyanobacteria, CsoSCA, remains limited, particularly regarding the regulation of its activity. Here, we present a structural and biochemical study of CsoSCA from the cyanobacterium Cyanobium sp. PCC7001. Our results show that the Cyanobium CsoSCA is allosterically activated by the Rubisco substrate ribulose-1,5-bisphosphate and forms a hexameric trimer of dimers. Comprehensive phylogenetic and mutational analyses are consistent with this regulation appearing exclusively in cyanobacterial α-carboxysome CAs. These findings clarify the biologically relevant oligomeric state of α-carboxysomal CAs and advance our understanding of the regulation of photosynthesis in this globally dominant lineage.

Assessing nanobody interaction with SARS-CoV-2 Nsp9.

The interaction between SARS-CoV-2 non-structural protein Nsp9 and the nanobody 2NSP90 was investigated by NMR spectroscopy using the paramagnetic perturbation methodology PENELOP (Paramagnetic Equilibrium vs Nonequilibrium magnetization Enhancement or LOss Perturbation). The Nsp9 monomer is an essential component of the replication and transcription complex (RTC) that reproduces the viral gRNA for subsequent propagation. Therefore preventing Nsp9 recruitment in RTC would represent an efficient antiviral strategy that could be applied to different coronaviruses, given the Nsp9 relative invariance. The NMR results were consistent with a previous characterization suggesting a 4:4 Nsp9-to-nanobody stoichiometry with the occurrence of two epitope pairs on each of the Nsp9 units that establish the inter-dimer contacts of Nsp9 tetramer. The oligomerization state of Nsp9 was also analyzed by molecular dynamics simulations and both dimers and tetramers resulted plausible. A different distribution of the mapped epitopes on the tetramer surface with respect to the former 4:4 complex could also be possible, as well as different stoichiometries of the Nsp9-nanobody assemblies such as the 2:2 stoichiometry suggested by the recent crystal structure of the Nsp9 complex with 2NSP23 (PDB ID: 8dqu), a nanobody exhibiting essentially the same affinity as 2NSP90. The experimental NMR evidence, however, ruled out the occurrence in liquid state of the relevant Nsp9 conformational change observed in the same crystal structure.

Reconstitution and characterization of BRAF in complex with 14-3-3 and KRAS4B on nanodiscs.

RAF kinases are key components of the RAS-MAPK signaling pathway, which drives cell growth and is frequently overactivated in cancer. Upstream signaling activates the small GTPase RAS, which recruits RAF to the cell membrane, driving a transition of the latter from an auto-inhibited monomeric conformation to an active dimer. Despite recent progress, mechanistic details underlying RAF activation remain unclear, particularly the role of RAS and the membrane in mediating this conformational rearrangement of RAF together with 14-3-3 to permit RAF kinase domain dimerization. Here, we reconstituted an active complex of dimeric BRAF, a 14-3-3 dimer and two KRAS4B on a nanodisc bilayer and verified that its assembly is GTP-dependent. Biolayer interferometry (BLI) was used to compare the binding affinities of monomeric versus dimeric full-length BRAF:14-3-3 complexes for KRAS4B-conjugated nanodiscs (RAS-ND) and to investigate the effects of membrane lipid composition and spatial density of KRAS4B on binding. 1,2-Dioleoyl-sn-glycero-3-phospho-L-serine (DOPS) and higher KRAS4B density enhanced the interaction of BRAF:14-3-3 with RAS-ND to different degrees depending on BRAF oligomeric state. We utilized our reconstituted system to dissect the effects of KRAS4B and the membrane on the kinase activity of monomeric and dimeric BRAF:14-3-3 complexes, finding that KRAS4B or nanodiscs alone were insufficient to stimulate activity, whereas RAS-ND increased activity of both states of BRAF. The reconstituted assembly of full-length BRAF with 14-3-3 and KRAS on a cell-free, defined lipid bilayer offers a more holistic biophysical perspective to probe regulation of this multimeric signaling complex at the membrane surface.

Post-translational modification-dependent oligomerization switch in regulation of global transcription and DNA damage repair during genotoxic stress.

Mechanisms of functional cross-talk between global transcriptional repression and efficient DNA damage repair during genotoxic stress are poorly known. In this study, using human AF9 as representative of Super Elongation Complex (SEC) components, we delineate detailed mechanisms of these processes. Mechanistically, we describe that Poly-Serine domain-mediated oligomerization is pre-requisite for AF9 YEATS domain-mediated TFIID interaction-dependent SEC recruitment at the promoter-proximal region for release of paused RNA polymerase II. Interestingly, during genotoxic stress, CaMKII-mediated phosphorylation-dependent nuclear export of AF9-specific deacetylase HDAC5 enhances concomitant PCAF-mediated acetylation of K339 residue. This causes monomerization of AF9 and reduces TFIID interaction for transcriptional downregulation. Furthermore, the K339 acetylation-dependent enhanced AF9-DNA-PKc interaction leads to phosphorylation at S395 residue which reduces AF9-SEC interaction resulting in transcriptional downregulation and efficient repair of DNA damage. After repair, nuclear re-entry of HDAC5 reduces AF9 acetylation and restores its TFIID and SEC interaction to restart transcription.

Structures of H5N1 influenza polymerase with ANP32B reveal mechanisms of genome replication and host adaptation.

Avian influenza A viruses (IAVs) pose a public health threat, as they are capable of triggering pandemics by crossing species barriers. Replication of avian IAVs in mammalian cells is hindered by species-specific variation in acidic nuclear phosphoprotein 32 (ANP32) proteins, which are essential for viral RNA genome replication. Adaptive mutations enable the IAV RNA polymerase (FluPolA) to surmount this barrier. Here, we present cryo-electron microscopy structures of monomeric and dimeric avian H5N1 FluPolA with human ANP32B. ANP32B interacts with the PA subunit of FluPolA in the monomeric form, at the site used for its docking onto the C-terminal domain of host RNA polymerase II during viral transcription. ANP32B acts as a chaperone, guiding FluPolA towards a ribonucleoprotein-associated FluPolA to form an asymmetric dimer-the replication platform for the viral genome. These findings offer insights into the molecular mechanisms governing IAV genome replication, while enhancing our understanding of the molecular processes underpinning mammalian adaptations in avian-origin FluPolA.

Cyanovirin-N Binding to N-Acetyl-d-glucosamine Requires Carbohydrate-Binding Sites on Two Different Protomers.

Cyanovirin-N (CV-N) binds high-mannose oligosaccharides on enveloped viruses with two carbohydrate-binding sites, one bearing high affinity and one low affinity to Manα(1-2)Man moieties. A tandem repeat of two CV-N molecules (CVN2) was tested for antiviral activity against human immunodeficiency virus type I (HIV-1) by using a domain-swapped dimer. CV-N was shown to bind N-acetylmannosamine (ManNAc) and N-acetyl-d-glucosamine (GlcNAc) when the carbohydrate-binding sites in CV-N were free to interact with these monosaccharides independently. CVN2 recognized ManNAc at a Kd of 1.4 µM and bound this sugar in solution, regardless of the lectin making amino acid side chain contacts on the targeted viral glycoproteins. An interdomain cross-contacting residue Glu41, which has been shown to be hydrogen bonding with dimannose, was substituted in the monomeric CV-N. The amide derivative of glucose, GlcNAc, achieved similar high affinity to the new variant CVN-E41T as high-mannose N-glycans, but binding to CVN2 in the nanomolar range with four binding sites involved or binding to the monomeric CVN-E41A. A stable dimer was engineered and expressed from the alanine-to-threonine-substituted monomer to confirm binding to GlcNAc. In summary, low-affinity binding was achieved by CVN2 to dimannosylated peptide or GlcNAc with two carbohydrate-binding sites of differing affinities, mimicking biological interactions with the respective N-linked glycans of interest and cross-linking of carbohydrates on human T cells for lymphocyte activation.

Trimeric Bet v 1-specific nanobodies cause strong suppression of IgE binding.

Background: Around 20% of the population in Northern and Central Europe is affected by birch pollen allergy, with the major birch pollen allergen Bet v 1 as the main elicitor of allergic reactions. Together with its cross-reactive allergens from related trees and foods, Bet v 1 causes an impaired quality of life. Hence, new treatment strategies were elaborated, demonstrating the effectiveness of blocking IgG antibodies on Bet v 1-induced IgE-mediated reactions. A recent study provided evidence for the first time that Bet v 1-specific nanobodies reduce patients´ IgE binding to Bet v 1. In order to increase the potential to outcompete IgE recognition of Bet v 1 and to foster cross-reactivity and cross-protection, we developed Bet v 1-specific nanobody trimers and evaluated their capacity to suppress polyclonal IgE binding to corresponding allergens and allergen-induced basophil degranulation. Methods: Nanobody trimers were engineered by adding isoleucine zippers, thus enabling trimeric formation. Trimers were analyzed for their cross-reactivity, binding kinetics to Bet v 1, and related allergens, and patients' IgE inhibition potential. Finally, their efficacy to prevent basophil degranulation was investigated. Results: Trimers showed enhanced recognition of cross-reactive allergens and increased efficiency to reduce IgE-allergen binding compared to nanobody monomers. Furthermore, trimers displayed slow dissociation rates from allergens and suppressed allergen-induced mediator release. Conclusion: We generated high-affine nanobody trimers that target Bet v 1 and related allergens. Trimers blocked IgE-allergen interaction by competing with IgE for allergen binding. They inhibited IgE-mediated release of biological mediators, demonstrating a promising potential to prevent allergic reactions caused by Bet v 1 and relatives.

Chemokine receptor hetero-oligomers regulate monocyte chemotaxis.

It is known that stress influences immune cell function. The underlying molecular mechanisms are unclear. We recently reported that many chemokine receptors (CRs) heteromerize with α1-adrenoceptors (α1-ARs) through which CRs are regulated. Here, we show that arginine vasopressin receptor 1A (AVPR1A) heteromerizes with all human CRs, except chemokine (C-X-C motif) receptor (CXCR)1, in recombinant systems and that such heteromers are detectable in THP-1 cells and human monocytes. We demonstrate that ligand-free AVPR1A differentially regulates the efficacy of CR partners to mediate chemotaxis and that AVPR1A ligands disrupt AVPR1A:CR heteromers, which enhances chemokine (C-C motif) receptor (CCR)1-mediated chemotaxis and inhibits CCR2-, CCR8-, and CXCR4-mediated chemotaxis. Using bioluminescence resonance energy transfer to monitor G protein activation and CRISPR/Cas9 gene-edited THP-1 cells lacking AVPR1A or α1B-AR, we show that CRs that share the propensity to heteromerize with α1B/D-ARs and AVPR1A exist and function within interdependent hetero-oligomeric complexes through which the efficacy of CRs to mediate chemotaxis is controlled. Our findings suggest that hetero-oligomers composed of CRs, α1B/D-ARs, and AVPR1A may enable stress hormones to regulate immune cell trafficking.

Secreted dengue virus NS1 from infection is predominantly dimeric and in complex with high-density lipoprotein.

Severe dengue infections are characterized by endothelial dysfunction shown to be associated with the secreted nonstructural protein 1 (sNS1), making it an attractive vaccine antigen and biotherapeutic target. To uncover the biologically relevant structure of sNS1, we obtained infection-derived sNS1 (isNS1) from dengue virus (DENV)-infected Vero cells through immunoaffinity purification instead of recombinant sNS1 (rsNS1) overexpressed in insect or mammalian cell lines. We found that isNS1 appeared as an approximately 250 kDa complex of NS1 and ApoA1 and further determined the cryoEM structures of isNS1 and its complex with a monoclonal antibody/Fab. Indeed, we found that the major species of isNS1 is a complex of the NS1 dimer partially embedded in a high-density lipoprotein (HDL) particle. Crosslinking mass spectrometry studies confirmed that the isNS1 interacts with the major HDL component ApoA1 through interactions that map to the NS1 wing and hydrophobic domains. Furthermore, our studies demonstrated that the sNS1 in sera from DENV-infected mice and a human patient form a similar complex as isNS1. Our results report the molecular architecture of a biological form of sNS1, which may have implications for the molecular pathogenesis of dengue.

Allosteric coupling asymmetry mediates paradoxical activation of BRAF by type II inhibitors.

The type II class of RAF inhibitors currently in clinical trials paradoxically activate BRAF at subsaturating concentrations. Activation is mediated by induction of BRAF dimers, but why activation rather than inhibition occurs remains unclear. Using biophysical methods tracking BRAF dimerization and conformation, we built an allosteric model of inhibitor-induced dimerization that resolves the allosteric contributions of inhibitor binding to the two active sites of the dimer, revealing key differences between type I and type II RAF inhibitors. For type II inhibitors the allosteric coupling between inhibitor binding and BRAF dimerization is distributed asymmetrically across the two dimer binding sites, with binding to the first site dominating the allostery. This asymmetry results in efficient and selective induction of dimers with one inhibited and one catalytically active subunit. Our allosteric models quantitatively account for paradoxical activation data measured for 11 RAF inhibitors. Unlike type II inhibitors, type I inhibitors lack allosteric asymmetry and do not activate BRAF homodimers. Finally, NMR data reveal that BRAF homodimers are dynamically asymmetric with only one of the subunits locked in the active αC-in state. This provides a structural mechanism for how binding of only a single αC-in inhibitor molecule can induce potent BRAF dimerization and activation.

The complete assembly of human LAT1-4F2hc complex provides insights into its regulation, function and localisation.

The LAT1-4F2hc complex (SLC7A5-SLC3A2) facilitates uptake of essential amino acids, hormones and drugs. Its dysfunction is associated with many cancers and immune/neurological disorders. Here, we apply native mass spectrometry (MS)-based approaches to provide evidence of super-dimer formation (LAT1-4F2hc)2. When combined with lipidomics, and site-directed mutagenesis, we discover four endogenous phosphatidylethanolamine (PE) molecules at the interface and C-terminus of both LAT1 subunits. We find that interfacial PE binding is regulated by 4F2hc-R183 and is critical for regulation of palmitoylation on neighbouring LAT1-C187. Combining native MS with mass photometry (MP), we reveal that super-dimerization is sensitive to pH, and modulated by complex N-glycans on the 4F2hc subunit. We further validate the dynamic assemblies of LAT1-4F2hc on plasma membrane and in the lysosome. Together our results link PTM and lipid binding with regulation and localisation of the LAT1-4F2hc super-dimer.

Exploring Structural Insights of Aß42 and α-Synuclein Monomers and Heterodimer: A Comparative Study Using Implicit and Explicit Solvent Simulations.

Protein misfolding, aggregation, and fibril formation play a central role in the development of severe neurological disorders, including Alzheimer's and Parkinson's diseases. The structural stability of mature fibrils in these diseases is of great importance, as organisms struggle to effectively eliminate amyloid plaques. To address this issue, it is crucial to investigate the early stages of fibril formation when monomers aggregate into small, toxic, and soluble oligomers. However, these structures are inherently disordered, making them challenging to study through experimental approaches. Recently, it has been shown experimentally that amyloid-ß 42 (Aß42) and α-synuclein (α-Syn) can coassemble. This has motivated us to investigate the interaction between their monomers as a first step toward exploring the possibility of forming heterodimeric complexes. In particular, our study involves the utilization of various Amber and CHARMM force-fields, employing both implicit and explicit solvent models in replica exchange and conventional simulation modes. This comprehensive approach allowed us to assess the strengths and weaknesses of these solvent models and force fields in comparison to experimental and theoretical findings, ensuring the highest level of robustness. Our investigations revealed that Aß42 and α-Syn monomers can indeed form stable heterodimers, and the resulting heterodimeric model exhibits stronger interactions compared to the Aß42 dimer. The binding of α-Syn to Aß42 reduces the propensity of Aß42 to adopt fibril-prone conformations and induces significant changes in its conformational properties. Notably, in AMBER-FB15 and CHARMM36m force fields with the use of explicit solvent, the presence of Aß42 significantly increases the ß-content of α-Syn, consistent with the experiments showing that Aß42 triggers α-Syn aggregation. Our analysis clearly shows that although the use of implicit solvent resulted in too large compactness of monomeric α-Syn, structural properties of monomeric Aß42 and the heterodimer were preserved in explicit-solvent simulations. We anticipate that our study sheds light on the interaction between α-Syn and Aß42 proteins, thus providing the atom-level model required to assess the initial stage of aggregation mechanisms related to Alzheimer's and Parkinson's diseases.

Deciphering the Cofilin Oligomers via Intermolecular Disulfide Bond Formation: A Coarse-Grained Molecular Dynamics Approach to Understanding Cofilin's Regulation on Actin Filaments.

Cofilin, a key actin-binding protein, orchestrates the dynamics of the actomyosin network through its actin-severing activity and by promoting the recycling of actin monomers. Recent experiments suggest that cofilin forms functionally distinct oligomers via thiol post-translational modifications (PTMs) that promote actin nucleation and assembly. Despite these advances, the structural conformations of cofilin oligomers that modulate actin activity remain elusive because there are combinatorial ways to oxidize thiols in cysteines to form disulfide bonds rapidly. This study employs molecular dynamics simulations to investigate human cofilin 1 as a case study for exploring cofilin dimers via disulfide bond formation. Utilizing a biasing scheme in simulations, we focus on analyzing dimer conformations conducive to disulfide bond formation. Additionally, we explore potential PTMs arising from the examined conformational ensemble. Using the free energy profiling, our simulations unveil a range of probable cofilin dimer structures not represented in current Protein Data Bank entries. These candidate dimers are characterized by their distinct population distributions and relative free energies. Of particular note is a dimer featuring an interface between cysteines 139 and 147 residues, which demonstrates stable free energy characteristics and intriguingly symmetrical geometry. In contrast, the experimentally proposed dimer structure exhibits a less stable free energy profile. We also evaluate frustration quantification based on the energy landscape theory in the protein-protein interactions at the dimer interfaces. Notably, the 39-39 dimer configuration emerges as a promising candidate for forming cofilin tetramers, as substantiated by frustration analysis. Additionally, docking simulations with actin filaments further evaluate the stability of these cofilin dimer-actin complexes. Our findings thus offer a computational framework for understanding the role of thiol PTM of cofilin proteins in regulating oligomerization, and the subsequent cofilin-mediated actin dynamics in the actomyosin network.

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.