Molecular Basis of S100A1 Activation and Target Regulation Within Physiological Cytosolic Ca2+ Levels

Calcium signaling in mitochondrial intermembrane space

The mitochondrial intermembrane space (IMS) is a highly protected compartment, second only to the matrix. It is a crucial bridge, coordinating mitochondrial activities with cellular processes such as metabolites, protein, lipid, and ion exchange. This regulation influences signaling pathways for metabolic activities and cellular homeostasis. The IMS harbors various proteins critical for initiating apoptotic cascades and regulating reactive oxygen species production by controlling the respiratory chain. Calcium (Ca2+), a key intracellular secondary messenger, enter the mitochondrial matrix via the IMS, regulating mitochondrial bioenergetics, ATP production, modulating cell death pathways. IMS acts as a regulatory site for Ca2+ entry due to the presence of different Ca2+ sensors such as MICUs, solute carriers (SLCs); ion exchangers (LETM1/SCaMCs); S100A1, mitochondrial glycerol-3-phosphate dehydrogenase, and EFHD1, each with unique Ca2+ binding motifs and spatial localizations. This review primarily emphasizes the role of these IMS-localized Ca2+ sensors concerning their spatial localization, mechanism, and molecular functions. Additionally, we discuss how these sensors contribute to the progression and pathogenesis of various human health conditions and diseases.

Transcriptomic analysis of immune-related genes in Pacific white snook (Centropomus viridis) gills infected with the monogenean parasite Rhabdosynochus viridisi

The parasite Rhabdosynochus viridisi (Platyhelminthes: Monogenea) infects the Pacific white snook Centropomus viridis gills and can cause adverse effects in the aquaculture industry. The immune responses of Pacific white snook to monogenean infections are poorly understood. Thus, this study aimed to identify differentially expressed genes (DEGs) in the gills of Pacific white snook juveniles experimentally infected with R. viridisi, emphasizing immune-related genes and pathways activated or suppressed during the infection. RNA sequencing was performed on the gills of uninfected (control) and infected fish. The algorithm Seq2Fun was selected without a reference transcriptome to map the reads to transcripts of fishes available from a database for gene orthologs (EcoOmics) and obtain the counting table. The ExpressAnalyst software was used for differential expression and functional analyses. A total of 20,106 transcripts were found, and 1430 (7 %) were differentially expressed genes (DEGs) between infected and control groups. We identified 860 (60 %) downregulated and 570 (40 %) upregulated genes. Thirteen canonical pathways after the Kyoto Encyclopedia of Genes and Genomes (KEGG) database were overrepresented, and most of the DEGs were downregulated, suggesting the inactivation of these pathways. The functions of most of the DEGs with higher fold change found in this study are poorly understood in fish. Even though the well-known pro-inflammatory cytokines remained unchanged in infected gills of C. viridis, and transforming growth factor β (tgfβ) was downregulated, interleukin-17 ligands il17d and il17a/f1, as well as C-X-C motif chemokine receptor 2 (cxcr2) genes were upregulated, indicating that the infection with R. viridisi promotes Th17-like immunity. Overexpression of plasma B cell activity markers such as immunoglobulin light chain-like genes and the v-set pre-B cell surrogate light chain 3 (vpreb3) was also detected in this study. The possible implications of DEGs related to calcium imbalance, hypoxia adaptation, hemostasis, and immunity are discussed. These results will support future studies to improve the prevention and treatment of monogenean infections in finfish aquaculture.

S100A1's single cysteine is an indispensable redox switch for the protection against diastolic calcium waves in cardiomyocytes

The EF-hand calcium (Ca2+) sensor protein S100A1 combines inotropic with antiarrhythmic potency in cardiomyocytes (CMs). Oxidative posttranslational modification (ox-PTM) of S100A1's conserved, single-cysteine residue (C85) via reactive nitrogen species (i.e., S-nitrosylation or S-glutathionylation) has been proposed to modulate conformational flexibility of intrinsically disordered sequence fragments and to increase the molecule's affinity toward Ca2+. Considering the unknown biological functional consequence, we aimed to determine the impact of the C85 moiety of S100A1 as a potential redox switch. We first uncovered that S100A1 is endogenously glutathionylated in the adult heart in vivo. To prevent glutathionylation of S100A1, we generated S100A1 variants that were unresponsive to ox-PTMs. Overexpression of wild-type (WT) and C85-deficient S100A1 protein variants in isolated CM demonstrated equal inotropic potency, as shown by equally augmented Ca2+ transient amplitudes under basal conditions and β-adrenergic receptor (βAR) stimulation. However, in contrast, ox-PTM defective S100A1 variants failed to protect against arrhythmogenic diastolic sarcoplasmic reticulum (SR) Ca2+ waves and ryanodine receptor 2 (RyR2) hypernitrosylation during βAR stimulation. Despite diastolic performance failure, C85-deficient S100A1 protein variants exerted similar Ca2+-dependent interaction with the RyR2 than WT-S100A1. Dissecting S100A1's molecular structure-function relationship, our data indicate for the first time that the conserved C85 residue potentially acts as a redox switch that is indispensable for S100A1's antiarrhythmic but not its inotropic potency in CMs. We, therefore, propose a model where C85's ox-PTM determines S100A1's ability to beneficially control diastolic but not systolic RyR2 activity.NEW & NOTEWORTHY S100A1 is an emerging candidate for future gene-therapy treatment of human chronic heart failure. We aimed to study the significance of the conserved single-cysteine 85 (C85) residue in cardiomyocytes. We show that S100A1 is endogenously glutathionylated in the heart and demonstrate that this is dispensable to increase systolic Ca2+ transients, but indispensable for mediating S100A1's protection against sarcoplasmic reticulum (SR) Ca2+ waves, which was dependent on the ryanodine receptor 2 (RyR2) nitrosylation status.

Deciphering S100B Allosteric Signaling: The Role of a Peptide Target, TRTK-12, as an Ensemble Modulator

Allostery is an essential biological phenomenon in which perturbation at one site in a biomolecule elicits a functional response at a distal location(s). It is integral to biological processes, such as cellular signaling, metabolism, and transcription regulation. Understanding allostery is also crucial for rational drug discovery. In this work, we focus on an allosteric S100B protein that belongs to the S100 class of EF-hand Ca2+-binding proteins. The Ca2+-binding affinity of S100B is modulated allosterically by TRTK-12 peptide binding 25 Å away from the Ca2+-binding site. We investigated S100B allostery by carrying out nuclear magnetic resonance (NMR) measurements along with microsecond-long molecular dynamics (MD) simulations on S100B/Ca2+ with/without TRTK-12 at different NaCl salt concentrations. NMR HSQC results show that TRTK-12 reorganizes how S100B/Ca2+ responds to different salt concentrations at both orthosteric and allosteric sites. The MD data suggest that TRTK-12 breaks the dynamic aromatic and hydrogen-bond interactions (not observed in X-ray crystallographic structures) between the hinge/helix and Ca2+-binding EF-hand loop of the two subunits in the homodimeric protein. This triggers rearrangement in the protein network architectures and leads to allosteric communication. Finally, computational studies of S100B at distinct ionic strengths suggest that ligand-bound species are more robust to the changing environment relative to the S100B/Ca2+ complex.

Differential Expression of MED12-Associated Coding RNA Transcripts in Uterine Leiomyomas

Recent studies have demonstrated that somatic MED12 mutations in exon 2 occur at a frequency of up to 80% and have a functional role in leiomyoma pathogenesis. The objective of this study was to elucidate the expression profile of coding RNA transcripts in leiomyomas, with and without these mutations, and their paired myometrium. Next-generation RNA sequencing (NGS) was used to systematically profile the differentially expressed RNA transcripts from paired leiomyomas (n = 19). The differential analysis indicated there are 394 genes differentially and aberrantly expressed only in the mutated tumors. These genes were predominantly involved in the regulation of extracellular constituents. Of the differentially expressed genes that overlapped in the two comparison groups, the magnitude of change in gene expression was greater for many genes in tumors bearing MED12 mutations. Although the myometrium did not express MED12 mutations, there were marked differences in the transcriptome landscape of the myometrium from mutated and non-mutated specimens, with genes regulating the response to oxygen-containing compounds being most altered. In conclusion, MED12 mutations have profound effects on the expression of genes pivotal to leiomyoma pathogenesis in the tumor and the myometrium which could alter tumor characteristics and growth potential.

Myofilament-associated proteins with intrinsic disorder (MAPIDs) and their resolution by computational modeling

The cardiac sarcomere is a cellular structure in the heart that enables muscle cells to contract. Dozens of proteins belong to the cardiac sarcomere, which work in tandem to generate force and adapt to demands on cardiac output. Intriguingly, the majority of these proteins have significant intrinsic disorder that contributes to their functions, yet the biophysics of these intrinsically disordered regions (IDRs) have been characterized in limited detail. In this review, we first enumerate these myofilament-associated proteins with intrinsic disorder (MAPIDs) and recent biophysical studies to characterize their IDRs. We secondly summarize the biophysics governing IDR properties and the state-of-the-art in computational tools toward MAPID identification and characterization of their conformation ensembles. We conclude with an overview of future computational approaches toward broadening the understanding of intrinsic disorder in the cardiac sarcomere.

Molecular Insights into the Calcium Binding in Troponin C through a Molecular Dynamics Study

Calcium-binding proteins play critical roles in various biological processes such as signal transduction, cell growth, and transcription factor regulation. Ion binding and target binding of Ca²⁺-binding proteins are highly related. Therefore, understanding the ion binding mechanism will benefit the relevant inhibitor design toward the Ca²⁺-binding proteins. The EF-hand is the typical ion binding motif in Ca²⁺-binding proteins. Previous studies indicate that the ion binding affinity of the EF-hand increases with the peptide length, but this mechanism has not been fully understood. Herein, using molecular dynamics simulations, thermodynamic integration calculations, and molecular mechanics Poisson-Boltzmann surface area analysis, we systematically investigated four Ca²⁺-binding peptides containing the EF-hand loop in site III of rabbit skeletal troponin C. These four peptides have 13, 21, 26, and 34 residues. Our simulations reproduced the observed trend that the ion binding affinity increases with the peptide length. Our results implied that the E-helix motif preceding the EF-hand loop, likely the Phe99 residue in particular, plays a significant role in this regulation. The E-helix has a significant impact on the backbone and side-chain conformations of the Asp103 residue, rigidifying important hydrogen bonds in the EF-hand and decreasing the solvent exposure of the Ca²⁺ ion, hence leading to more favorable Ca²⁺ binding in longer peptides. The present study provides molecular insights into the ion binding in the EF-hand and establishes an important step toward elucidating the responses of Ca²⁺-binding proteins toward the ion and target availability.

S100A1 expression is increased in spinal cord injury and promotes inflammation, oxidative stress and apoptosis of PC12 cells induced by LPS via ERK signaling

Spinal cord injury (SCI) is a severe neurological disorder and the molecular mechanisms leading to its poor prognosis remain to be elucidated. S100A1, a mediator of Ca2+ handling of sarcoplasmic reticulum and mitochondrial function, operates as an endogenous danger signal (alarmin) associated with inflammatory response and tissue injury. The aim of the present study was to investigate the expression and biological effects of S100A1 in SCI. A rat model of SCI and a PC12 cell model of lipopolysaccharide (LPS)‑induced inflammation were established to examine S100A1 expression at the mRNA and protein levels. The inflammation level, which was mediated by S100A1, was determined based on inflammatory factor (IL‑1β, IL‑6 and TNF‑α) and anti‑inflammatory factor (IL‑10) expression. The effects of S100A1 on cellular oxidation and anti‑oxidation levels were observed by detecting the levels of reactive oxygen species, superoxide dismutase, catalase activities and nuclear factor erythroid 2‑related factor 2 expression. The protein levels of Bax, Bcl2 and cleaved caspase‑3 were used for the evaluation of the effects of S100A1 on apoptosis. Phosphorylated (p‑)ERK1/2 expression was used to evaluate the effects of S100A1 on ERK signaling. The results revealed that S100A1 expression was significantly upregulated in vivo and in vitro in the PC12 cell model of LPS‑inflammation. The silencing and overexpression of S100A1 helped alleviate and aggravate LPS‑induced inflammation, oxidative stress and apoptosis levels, respectively. S100A1 was found to regulate the ERK signaling pathway positively. An inhibitor of ERK signaling (MK‑8353) partially abolished the promoting effects of the overexpression of S100A1 on inflammation, oxidative stress damage and apoptosis. In conclusion, S100A1 expression was elevated in model of SCI and in the PC12 cell model of LPS‑induced inflammation. Furthermore, the overexpression/silencing S100A1 aggravated/mitigated the inflammation, oxidative stress damage and the apoptosis of LPS‑stimulated PC12 cells via the ERK signaling pathway. The present study revealed the mechanism of S100A1 in SCI, which provided a new theoretic reference for future research on SCI.

Overexpression of S100A1 in Osteosarcoma Inhibits Tumor Proliferation and Progression

Background

Osteosarcoma is the most common primary malignant tumor of bone. Abnormal expression of S100A1 protein is closely related to the occurrence and development of malignant tumors. However, S100A1 in osteosarcoma has not been studied.

Methods

All osteosarcoma tissues were collected from patients who received surgical therapy at the Affiliated Hospital of Inner Mongolia Medical University, China in 2020. QRT-PCR and western blot assays were used to detect the expression of S100A1 in osteosarcoma tissues and cells. The negative effect of S100A1 on osteosarcoma cell growth was confirmed by vitro and vivo experiments.

Results

S100A1 inhibited the growth of osteosarcoma cells in vitro. Overexpression of S100A1 may inhibit the proliferation of osteosarcoma cells by preventing the activation of AKT signaling pathway by western blot assay. Finally, animal experiments confirmed that overexpression of S100A1 could inhibit the proliferation of osteosarcoma cells. Overexpression of S100A1 obtained better survival benefit in mice.

Conclusion

Our findings provided a new insight to the treatment of osteosarcoma. It also raised the possibility that S100A1 could be used in targeted therapies for osteosarcoma.

Non-Canonical Interaction between Calmodulin and Calcineurin Contributes to the Differential Regulation of Plant-Derived Calmodulins on Calcineurin

Calmodulin (CaM) serves as an important Ca2+ signaling hub that regulates many protein signaling pathways. Recently, it was demonstrated that plant CaM homologues can regulate mammalian targets, often in a manner that opposes the impact of the mammalian CaM (mCaM). However, the molecular basis of how CaM homologue mutations differentially impact target activation is unclear. To understand these mechanisms, we examined two CaM isoforms found in soybean plants that differentially regulate a mammalian target, calcineurin (CaN). These CaM isoforms, sCaM-1 and sCaM-4, share >90 and ∼78% identity with the mCaM, respectively, and activate CaN with comparable or reduced activity relative to mCaM. We used molecular dynamics (MD) simulations and fluorometric assays of CaN-dependent dephosphorylation of MUF-P to probe whether calcium and protein-protein binding interactions are altered by plant CaMs relative to mCaM as a basis for differential CaN regulation. In the presence of CaN, we found that the two sCaMs' Ca2+ binding properties, such as their predicted coordination of Ca2+ and experimentally measured EC50 [Ca2+] values are comparable to mCaM. Furthermore, the binding of CaM to the CaM binding region (CaMBR) in CaN is comparable among the three CaMs, as evidenced by MD-predicted binding energies and experimentally measured EC50 [CaM] values. However, mCaM and sCaM-1 exhibited binding with a secondary region of CaN's regulatory domain that is weakened for sCaM-4. We speculate that this secondary interaction affects the turnover rate (kcat) of CaN based on our modeling of enzyme activity, which is consistent with our experimental data. Together, our data describe how plant-derived CaM variants alter CaN activity through enlisting interactions other than those directly influencing Ca2+ binding and canonical CaMBR binding, which may additionally play a role in the differential regulation of other mammalian targets.

Computational Analysis of Binding Interactions between the Ryanodine Receptor Type 2 and Calmodulin

Mutations in the cardiac ryanodine receptor type 2 (RyR2) have been linked to a variety of cardiac arrhythmias, such as catecholaminergic polymorphic ventricular tachycardia (CPVT). RyR2 is regulated by calmodulin (CaM), and mutations that disrupt their interaction can cause aberrant calcium release, leading to an arrhythmia. It was recently shown that increasing the RyR2-CaM binding affinity could rescue a defective CPVT-related RyR2 channel to near wild-type behavior. However, the interactions that determine the binding affinity at the RyR2-CaM binding interface are not well understood. In this study, we identify the key domains and interactions, including several new interactions, involved in the binding of CaM to RyR2. Also, our comparison between the wild-type and V3599K mutant suggests how the RyR2-CaM binding affinity can be increased via a change in the central and N-terminal lobe binding contacts for CaM. This computational approach provides new insights into the effect of a mutation at the RyR2-CaM binding interface, and it may find utility in drug design for the future treatment of cardiac arrhythmias.

Structural Changes beyond the EF-Hand Contribute to Apparent Calcium Binding Affinities: Insights from Parvalbumins

Members of the parvalbumin (PV) family of calcium (Ca2+) binding proteins (CBPs) share a relatively high level of sequence similarity. However, their Ca2+ affinities and selectivities against competing ions like Mg2+ can widely vary. We conducted molecular dynamics simulations of several α-parvalbumin (αPV) constructs with micromolar to nanomolar Ca2+ affinities to identify structural and dynamic features that contribute to their binding of ions. Specifically, we examined a D94S/G98E construct with a lower Ca2+ affinity (≈-18 kcal/mol) relative to the wild type (WT) (≈-22 kcal/mol) and an S55D/E59D variant with enhanced affinity (≈-24 kcal/mol). Additionally, we also examined the binding of Mg2+ to these isoforms, which is much weaker than Ca2+. We used mean spherical approximation (MSA) theory to evaluate ion binding thermodynamics within the proteins' EF-hand domains to account for the impact of ions' finite sizes and the surrounding electrolyte composition. While the MSA scores differentiated Mg2+ from Ca2+, they did not indicate that Ca2+ binding affinities at the binding loop differed between the PV isoforms. Instead, molecular mechanics generalized Born surface area (MM/GBSA) approximation energies, which we used to quantify the thermodynamic cost of structural rearrangement of the proteins upon binding ions, indicated that S55D/E59D αPV favored Ca2+ binding by -20 kcal/mol relative to WT versus 30 kcal/mol for D94S/G98E αPV. Meanwhile, Mg2+ binding was favored for the S55D/E59D αPV and D94S/G98E αPV variants by -18.32 and -1.65 kcal/mol, respectively. These energies implicate significant contributions to ion binding beyond oxygen coordination at the binding loop, which stemmed from changes in α-helicity, β-sheet character, and hydrogen bonding. Hence, Ca2+ affinity and selectivity against Mg2+ are emergent properties stemming from both local effects within the proteins' ion binding sites as well as non-local contributions elsewhere. Our findings broaden our understanding of the molecular bases governing αPV ion binding that are likely shared by members of the broad family of CBPs.