Solution structure and dynamics of S100A5 in the apo and Ca2+-bound states

Binding and Functional Folding (BFF): A Physiological Framework for Studying Biomolecular Interactions and Allostery

EF-hand Ca2+-binding proteins (CBPs), such as S100 proteins (S100s) and calmodulin (CaM), are signaling proteins that undergo conformational changes upon increasing intracellular Ca2+. Upon binding Ca2+, S100 proteins and CaM interact with protein targets and induce important biological responses. The Ca2+-binding affinity of CaM and most S100s in the absence of target is weak (CaKD > 1 μM). However, upon effector protein binding, the Ca2+ affinity of these proteins increases via heterotropic allostery (CaKD < 1 μM). Because of the high number and micromolar concentrations of EF-hand CBPs in a cell, at any given time, allostery is required physiologically, allowing for (i) proper Ca2+ homeostasis and (ii) strict maintenance of Ca2+-signaling within a narrow dynamic range of free Ca2+ ion concentrations, [Ca2+]free. In this review, mechanisms of allostery are coalesced into an empirical "binding and functional folding (BFF)" physiological framework. At the molecular level, folding (F), binding and folding (BF), and BFF events include all atoms in the biomolecular complex under study. The BFF framework is introduced with two straightforward BFF types for proteins (type 1, concerted; type 2, stepwise) and considers how homologous and nonhomologous amino acid residues of CBPs and their effector protein(s) evolved to provide allosteric tightening of Ca2+ and simultaneously determine how specific and relatively promiscuous CBP-target complexes form as both are needed for proper cellular function.

Were Ancestral Proteins Less Specific?

Some have hypothesized that ancestral proteins were, on average, less specific than their descendants. If true, this would provide a universal axis along which to organize protein evolution and suggests that reconstructed ancestral proteins may be uniquely powerful tools for protein engineering. Ancestral sequence reconstruction studies are one line of evidence used to support this hypothesis. Previously, we performed such a study, investigating the evolution of peptide-binding specificity for the paralogs S100A5 and S100A6. The modern proteins appeared more specific than their last common ancestor (ancA5/A6), as each paralog bound a subset of the peptides bound by ancA5/A6. In this study, we revisit this transition, using quantitative phage display to measure the interactions of 30,533 random peptides with human S100A5, S100A6, and ancA5/A6. This unbiased screen reveals a different picture. While S100A5 and S100A6 do indeed bind to a subset of the peptides recognized by ancA5/A6, they also acquired new peptide partners outside of the set recognized by ancA5/A6. Our previous work showed that ancA5/A6 had lower specificity than its descendants when measured against biological targets; our new work shows that ancA5/A6 has similar specificity to the modern proteins when measured against a random set of peptide targets. This demonstrates that altered biological specificity does not necessarily indicate altered intrinsic specificity, and sounds a cautionary note for using ancestral reconstruction studies with biological targets as a means to infer global evolutionary trends in specificity.

SRLS Analysis of 15N-1H NMR Relaxation from the Protein S100A1: Dynamic Structure, Calcium Binding, and Related Changes in Conformational Entropy

We report on amide (N-H) NMR relaxation from the protein S100A1 analyzed with the slowly relaxing local structure (SRLS) approach. S100A1 comprises two calcium-binding "EF-hands" (helix-loop-helix motifs) connected by a linker. The dynamic structure of this protein, in both calcium-free and calcium-bound form, is described as the restricted local N-H motion coupled to isotropic protein tumbling. The restrictions are given by a rhombic potential, u (∼10 kT), the local motion by a diffusion tensor with rate constant D₂ (∼10⁹ s^-1), and principal axis tilted from the N-H bond at angle β (10-20°). This parameter combination provides a physically insightful picture of the dynamic structure of S100A1 from the N-H bond perspective. Calcium binding primarily affects the C-terminal EF-hand, among others slowing down the motion of helices III and IV approximately 10-fold. Overall, it brings about significant changes in the shape of the local potential, u, and the orientation of the local diffusion axis, β. Conformational entropy derived from u makes an unfavorable entropic contribution to the free energy of calcium binding estimated at 8.6 ± 0.5 kJ/mol. The N-terminal EF-hand undergoes moderate changes. These findings provide new insights into the calcium-binding process. The same data were analyzed previously with the extended model-free (EMF) method, which is a simple limit of SRLS. In that interpretation, the protein tumbles anisotropically. Locally, calcium binding increases ordering in the loops of S100A1 and conformational exchange (R_ex) in the helices of its N-terminal EF-hand. These are very unusual features. We show that they most likely stem from problematic data-fitting, oversimplifications inherent in EMF, and experimental imperfections. R_ex is shown to be mainly a fit parameter. By reanalyzing the experimental data with SRLS, which is largely free of these deficiencies, we obtain-as delineated above-physically-relevant structural, kinetic, geometric, and binding information.

Learning peptide recognition rules for a low-specificity protein

Many proteins interact with short linear regions of target proteins. For some proteins, however, it is difficult to identify a well-defined sequence motif that defines its target peptides. To overcome this difficulty, we used supervised machine learning to train a model that treats each peptide as a collection of easily-calculated biochemical features rather than as an amino acid sequence. As a test case, we dissected the peptide-recognition rules for human S100A5 (hA5), a low-specificity calcium binding protein. We trained a Random Forest model against a recently released, high-throughput phage display dataset collected for hA5. The model identifies hydrophobicity and shape complementarity, rather than polar contacts, as the primary determinants of peptide binding specificity in hA5. We tested this hypothesis by solving a crystal structure of hA5 and through computational docking studies of diverse peptides onto hA5. These structural studies revealed that peptides exhibit multiple binding modes at the hA5 peptide interface-all of which have few polar contacts with hA5. Finally, we used our trained model to predict new, plausible binding targets in the human proteome. This revealed a fragment of the protein α-1-syntrophin that binds to hA5. Our work helps better understand the biochemistry and biology of hA5, as well as demonstrating how high-throughput experiments coupled with machine learning of biochemical features can reveal the determinants of binding specificity in low-specificity proteins.

Zooming into the Dark Side of Human Annexin-S100 Complexes: Dynamic Alliance of Flexible Partners

Annexins and S100 proteins form two large families of Ca²⁺-binding proteins. They are quite different both structurally and functionally, with S100 proteins being small (10–12 kDa) acidic regulatory proteins from the EF-hand superfamily of Ca²⁺-binding proteins, and with annexins being at least three-fold larger (329 ± 12 versus 98 ± 7 residues) and using non-EF-hand-based mechanism for calcium binding. Members of both families have multiple biological roles, being able to bind to a large cohort of partners and possessing a multitude of functions. Furthermore, annexins and S100 proteins can interact with each other in either a Ca²⁺-dependent or Ca²⁺-independent manner, forming functional annexin-S100 complexes. Such functional polymorphism and binding indiscrimination are rather unexpected, since structural information is available for many annexins and S100 proteins, which therefore are considered as ordered proteins that should follow the classical “one protein–one structure–one function” model. On the other hand, the ability to be engaged in a wide range of interactions with multiple, often unrelated, binding partners and possess multiple functions represent characteristic features of intrinsically disordered proteins (IDPs) and intrinsically disordered protein regions (IDPRs); i.e., functional proteins or protein regions lacking unique tertiary structures. The aim of this paper is to provide an overview of the functional roles of human annexins and S100 proteins, and to use the protein intrinsic disorder perspective to explain their exceptional multifunctionality and binding promiscuity.

Friend or Foe: S100 Proteins in Cancer

S100 proteins are widely expressed small molecular EF-hand calcium-binding proteins of vertebrates, which are involved in numerous cellular processes, such as Ca²⁺ homeostasis, proliferation, apoptosis, differentiation, and inflammation. Although the complex network of S100 signalling is by far not fully deciphered, several S100 family members could be linked to a variety of diseases, such as inflammatory disorders, neurological diseases, and also cancer. The research of the past decades revealed that S100 proteins play a crucial role in the development and progression of many cancer types, such as breast cancer, lung cancer, and melanoma. Hence, S100 family members have also been shown to be promising diagnostic markers and possible novel targets for therapy. However, the current knowledge of S100 proteins is limited and more attention to this unique group of proteins is needed. Therefore, this review article summarises S100 proteins and their relation in different cancer types, while also providing an overview of novel therapeutic strategies for targeting S100 proteins for cancer treatment.

Reviewing the Crystal Structure of S100Z and Other Members of the S100 Family: Implications in Calcium-Regulated Quaternary Structure

This paper takes the cue from the previously solved crystal structure of human apo-S100Z and compares it with that of the calcium-bound S100Z from zebrafish in order to stress, for this particular S100, the significant role of the presence of calcium in promoting supramolecular assemblies with likely biological meaning. This consideration is then expanded through a wider review on analogous situations concerning all other S100s for which there is crystallographic o biochemical evidence of how the presence of calcium promotes the formation of quaternary complexes.The paper also deals with some considerations on the quality of the crystals obtained for the solved members of this family and on the need for experimental phasing for solving some of the structures where the good general sequence homology among the members of the family would have suggested molecular replacement (MR) as the easiest way to solve them.These considerations, along with the PCA analysis carried out on all the known S100s, further demonstrate that calcium plays a fundamental role in triggering quaternary structure formation for several members of this family of proteins.

Conservation of Specificity in Two Low-Specificity Proteins

Many regulatory proteins bind peptide regions of target proteins and modulate their activity. Such regulatory proteins can often interact with highly diverse target peptides. In many instances, it is not known if the peptide-binding interface discriminates targets in a biological context, or whether biological specificity is achieved exclusively through external factors such as subcellular localization. We used an evolutionary biochemical approach to distinguish these possibilities for two such low-specificity proteins: S100A5 and S100A6. We used isothermal titration calorimetry to study the binding of peptides with diverse sequence and biochemistry to human S100A5 and S100A6. These proteins bound distinct, but overlapping, sets of peptide targets. We then studied the peptide binding properties of orthologs sampled from across five amniote species. Binding specificity was conserved along all lineages, for the last 320 million years, despite the low specificity of each protein. We used ancestral sequence reconstruction to determine the binding specificity of the last common ancestor of the paralogs. The ancestor bound the entire set of peptides bound by modern S100A5 and S100A6 proteins, suggesting that paralog specificity evolved via subfunctionalization. To rule out the possibility that specificity is conserved because it is difficult to modify, we identified a single historical mutation that, when reverted in human S100A5, gave it the ability to bind an S100A6-specific peptide. These results reveal strong evolutionary constraints on peptide binding specificity. Despite being able to bind a large number of targets, the specificity of S100 peptide interfaces is likely important for the biology of these proteins.

Human S100A5 binds Ca2+ and Cu2+ independently

Background

S100A5 is a calcium binding protein found in a small subset of amniote tissues. Little is known about the biological roles of S100A5, but it may be involved in inflammation and olfactory signaling. Previous work indicated that S100A5 displays antagonism between binding of Ca²⁺ and Cu²⁺ ions-one of the most commonly cited features of the protein. We set out to characterize the interplay between Ca²⁺ and Cu²⁺ binding by S100A5 using isothermal titration calorimetry (ITC), circular dichroism spectroscopy (CD), and analytical ultracentrifugation (AUC).

Results

We found that human S100A5 is capable of binding both Cu²⁺ and Ca²⁺ ions simultaneously. The wildtype protein was extremely aggregation-prone in the presence of Cu²⁺ and Ca²⁺. A Cys-free version of S100A5, however, was not prone to precipitation or oligomerization. Mutation of the cysteines does not disrupt the binding of either Ca²⁺ or Cu²⁺ to S100A5. In the Cys-free background, we measured Ca²⁺ and Cu²⁺ binding in the presence and absence of the other metal using ITC. Saturating concentrations of Ca²⁺ or Cu²⁺ do not disrupt the binding of one another. Ca²⁺ and Cu²⁺ binding induce structural changes in S100A5, which are measurable using CD spectroscopy. We show via sedimentation velocity AUC that the wildtype protein is prone to the formation of soluble oligomers, which are not present in Cys-free samples.

Conclusions

S100A5 can bind Ca²⁺ and Cu²⁺ ions simultaneously and independently. This observation is in direct contrast to previously-reported antagonism between binding of Cu²⁺ and Ca²⁺ ions. The previous result is likely due to metal-dependent aggregation. Little is known about the biology of S100A5, so an accurate understanding of the biochemistry is necessary to make informed biological hypotheses. Our observations suggest the possibility of independent biological functions for Cu²⁺ and Ca²⁺ binding by S100A5.

Solving the crystal structure of human calcium-free S100Z: the siege and conquer of one of the last S100 family strongholds

The X-ray structure of human apo-S100Z has been solved and compared with that of the zebrafish calcium-bound S100Z, which is the closest in sequence. Human apo-S100A12, which shows only 43% sequence identity to human S100Z, has been used as template model to solve the crystallographic phase problem. Although a significant buried surface area between the two physiological dimers is present in the asymmetric unit of human apo-S100Z, the protein does not form the superhelical arrangement in the crystal as observed for the zebrafish calcium-bound S100Z and human calcium-bound S100A4. These findings further demonstrate that calcium plays a fundamental role in triggering quaternary structure formation in several S100s. Solving the X-ray structure of human apo-S100Z by standard molecular replacement procedures turned out to be a challenge and required trying different models and different software tools among which only one was successful. The model that allowed structure solution was that with one of the lowest sequence identity with the target protein among the S100 family in the apo state. Based on the previously solved zebrafish holo-S100Z, a putative human holo-S100Z structure has been then calculated through homology modeling; the differences between the experimental human apo and calculated holo structure have been compared to those existing for other members of the family.

Biophysical characterization of Ca2+-binding of S100A5 and Ca2+-induced interaction with RAGE

S100A5 is a calcium-binding protein of S100 family, which represents a major ligand to the receptor for advanced glycation end product (RAGE), a pattern recognition receptor engaged in diverse pathological processes. Here we have characterized calcium binding of S100A5 and the complex formation between S100A5 and RAGE using calorimetry and NMR spectroscopy. S100A5 binds to calcium ions in a sequential manner with the equilibrium dissociation constants (K_D) of 1.3 μM and 3.5 μM, which corresponds to the calcium-binding at the C-terminal and N-terminal EF-hands. Upon calcium binding, S100A5 interacts with the V domain of RAGE (RAGE-v) to form a heterotrimer (K_D ∼5.9 μM) that is distinct among the S100 family proteins. Chemical shift perturbation data from NMR titration experiments indicates that S100A5 employs the periphery of the dimer interface to interact with RAGE-v. Distinct binding mode and stoichiometry of RAGE against different S100 family proteins could be important to modulate diverse RAGE signaling.

Multiple Evolutionary Origins of Ubiquitous Cu2+ and Zn2+ Binding in the S100 Protein Family

The S100 proteins are a large family of signaling proteins that play critical roles in biology and disease. Many S100 proteins bind Zn2+, Cu2+, and/or Mn2+ as part of their biological functions; however, the evolutionary origins of binding remain obscure. One key question is whether divalent transition metal binding is ancestral, or instead arose independently on multiple lineages. To tackle this question, we combined phylogenetics with biophysical characterization of modern S100 proteins. We demonstrate an earlier origin for established S100 subfamilies than previously believed, and reveal that transition metal binding is widely distributed across the tree. Using isothermal titration calorimetry, we found that Cu2+ and Zn2+ binding are common features of the family: the full breadth of human S100 paralogs-as well as two early-branching S100 proteins found in the tunicate Oikopleura dioica-bind these metals with μM affinity and stoichiometries ranging from 1:1 to 3:1 (metal:protein). While binding is consistent across the tree, structural responses to binding are quite variable. Further, mutational analysis and structural modeling revealed that transition metal binding occurs at different sites in different S100 proteins. This is consistent with multiple origins of transition metal binding over the evolution of this protein family. Our work reveals an evolutionary pattern in which the overall phenotype of binding is a constant feature of S100 proteins, even while the site and mechanism of binding is evolutionarily labile.

Multilevel Changes in Protein Dynamics upon Complex Formation of the Calcium-Loaded S100A4 with a Nonmuscle Myosin IIA Tail Fragment

Dysregulation of Ca²⁺ -binding S100 proteins plays important role in various diseases. The asymmetric complex of Ca²⁺ -bound S100A4 with nonmuscle myosin IIA has high stability and highly increased Ca²⁺ affinity. Here we investigated the possible causes of this allosteric effect by NMR spectroscopy. Chemical shift-based secondary-structure analysis did not show substantial changes for the complex. Backbone dynamics revealed slow-timescale local motions in the H1 helices of homodimeric S100A4; these were less pronounced in the complex form and might be accompanied by an increase in dimer stability. Different mobilities in the Ca²⁺ -coordinating EF-hand sites indicate that they communicate by an allosteric mechanism operating through changes in protein dynamics; this must be responsible for the elevated Ca²⁺ affinity. These multilevel changes in protein dynamics as conformational adaptation allow S100A4 fine-tuning of its protein-protein interactions inside the cell during Ca²⁺ signaling.

Pentamidine blocks the interaction between mutant S100A5 and RAGE V domain and inhibits the RAGE signaling pathway

The human S100 protein family contains small, dimeric and acidic proteins that contain two EF-hand motifs and bind calcium. When S100A5 binds calcium, its conformation changes and promotes interaction with the target protein. The extracellular domain of RAGE (Receptor of Advanced Glycation End products) contain three domains: C1, C2 and V. The RAGE V domain is the target protein of S100A5 that promotes cell survival, growth and differentiation by activating several signaling pathways. Pentamidine is an apoptotic and antiparasitic drug that is used to treat or prevent pneumonia. Here, we found that pentamidine interacts with S100A5 using HSQC titration. We elucidated the interactions of S100A5 with RAGE V domain and pentamidine using fluorescence and NMR spectroscopy. We generated two binary models-the S100A5-RAGE V domain and S100A5-Pentamidine complex-and then observed that the pentamidine and RAGE V domain share a similar binding region in mS100A5. We also used the WST-1 assay to investigate the bioactivity of S100A5, RAGE V domain and pentamidine. These results indicated that pentamidine blocks the binding between S100A5 and RAGE V domain. This finding is useful for the development of new anti-proliferation drugs.

Novel interactions of the TRTK12 peptide with S100 protein family members: specificity and thermodynamic characterization

The S100 protein family consists of small, dimeric proteins that exert their biological functions in response to changing calcium concentrations. S100B is the best-studied member and has been shown to interact with more than 20 binding partners in a calcium-dependent manner. The TRTK12 peptide, derived from the consensus binding sequence for S100B, has previously been found to interact with S100A1 and has been proposed to be a general binding partner of the S100 family. To test this hypothesis and gain a better understanding of the specificity of binding for the S100 proteins, 16 members of the human S100 family were screened against this peptide and its alanine variants. Novel interactions were found with only two family members, S100P and S100A2, indicating that TRTK12 selectively interacts with a small subset of the S100 proteins. Substantial promiscuity was observed in the binding site of S100B thereby accommodating variations in the peptide sequence, while S100A1, S100A2, and S100P exhibited larger differences in the binding constants for the TRTK12 alanine variants. This suggests that single-point substitutions can be used to selectively modulate the affinity of TRTK12 peptides for individual S100 proteins. This study has important implications for the rational drug design of inhibitors for the S100 proteins, which are involved in a variety of cancers and neurodegenerative diseases.

Solution structure and dynamics of human S100A14

Human S100A14 is a member of the EF-hand calcium-binding protein family that has only recently been described in terms of its functional and pathological properties. The protein is overexpressed in a variety of tumor cells and it has been shown to trigger receptor for advanced glycation end products (RAGE)-dependent signaling in cell cultures. The solution structure of homodimeric S100A14 in the apo state has been solved at physiological temperature. It is shown that the protein does not bind calcium(II) ions and exhibits a "semi-open" conformation that thus represents the physiological structure of the S100A14. The lack of two ligands in the canonical EF-hand calcium(II)-binding site explains the negligible affinity for calcium(II) in solution, and the exposed cysteines and histidine account for the observed precipitation in the presence of zinc(II) or copper(II) ions.

Intrinsic disorder in S100 proteins

Although the members of the largest subfamily of the EF-hand proteins, S100 proteins, are evolutionarily young, their functional diversity is extremely broad, partly due to their ability to adapt to various targets. This feature is a hallmark of intrinsically disordered proteins (IDPs), but none of the S100 proteins are recognized as IDPs. S100 are predicted to be enriched in intrinsic disorder, with 62% of them being predicted to be disordered by at least one of the predictors: 31% are recognized as 'molten globules' and 15% are shown to be in extended disordered form. The disorder level of predicted disordered S100 regions is conserved compared to that of more structured regions. The central disordered stretch corresponds to the major part of pseudo EF-hand loop, helix II, hinge region, and an initial part of helix III. It contains about half of known sites of enzymatic post-translational modifications (PTMs), confirming that this region can be flexible in vivo. Most of the internal residues missing in tertiary structures belong to the hinge. Both hinge and pseudo EF-hand loop correspond to the local maxima of the PONDR® VSL2 score and are shown to be evolutionary hotspots, leading to gain of new functional properties. The action of PTMs is shown to be destabilizing, in contrast with the effect of metal-binding or S100 dimerization. Formation of the S100 heterodimers relies on the interplay between the structural rigidity of one of the S100 monomers and the flexibility of another monomer. The ordered regions dominate in the S100 homodimerization sites. Target-binding sites generally consist of distant regions, drastically differing in their disorder level. The disordered region comprising most of the hinge and the N-terminal half of helix III is virtually not involved into dimerization, being intended solely for target recognition. The structural flexibility of this region is essential for recognition of diverse target proteins. At least 86% of multiple interactions of S100 proteins with binding partners are attributed to the S100 proteins predicted to be disordered. Overall, the intrinsic disorder is inherent to many S100 proteins and is vital for activity and functional diversity of the family.

Structural characterization of human S100A16, a low-affinity calcium binder

The homodimeric structure of human S100A16 in the apo state has been obtained both in the solid state and in solution, resulting in good agreement between the structures with the exception of two loop regions. The homodimeric solution structure of human S100A16 was also calculated in the calcium(II)-bound form. Differently from most S100 proteins, the conformational rearrangement upon calcium binding is minor. This characteristic is likely to be related to the weak binding affinity of the protein for the calcium(II) ions. In turn, this is ascribed to the lack of the glutamate residue at the end of the S100-specific N-domain binding site, which in most S100 proteins provides two important side chain oxygen atoms as calcium(II) ligands. Furthermore, the presence of hydrophobic interactions stronger than for other S100 proteins, present in the closed form of S100A16 between the third and fourth helices, likely make the closed structure of the second EF-hand particularly stable, so even upon calcium(II) binding such a conformation is not disrupted.