Hydrophobic core substitutions in calbindin D9k: effects on Ca2+ binding and dissociation

Structural characterization and calcium absorption-promoting effect of sucrose-calcium chelate in Caco-2 monolayer cells and mice

Sucrose emerges as a chelating agent to form a stable sucrose-metal-ion chelate that can potentially improve metal-ion absorption. This study aimed to analyze the structure of sucrose-calcium chelate and its potential to promote calcium absorption in both Caco-2 monolayer cells and mice. The characterization results showed that calcium ions mainly chelated with hydroxyl groups in sucrose to produce sucrose-calcium chelate, altering the crystal structure of sucrose (forming polymer particles) and improving its thermal stability. Sucrose-calcium chelate dose dependently increased the amount of calcium uptake, retention, and transport in the Caco-2 monolayer cell model. Compared to CaCl2 , there was a significant improvement in the proportion of absorbed calcium utilized for transport but not retention (93.13 ± 1.75% vs. 67.67 ± 7.55%). Further treatment of calcium channel inhibitors demonstrated the active transport of sucrose-calcium chelate through Cav1.3. Cellular thermal shift assay and quantitative reverse transcription-polymerase chain reaction (qRT-PCR) assays indicated that the ability of sucrose-calcium chelate to promote calcium transport was attributed to its superior ability to bind with PMCA1b, a calcium transporter located on the basement membrane, and stimulate its gene expression compared to CaCl2 . Pharmacokinetic analysis of mice confirmed the calcium absorption-promoting effect of sucrose-calcium chelate, as evident by the higher serum calcium level (44.12 ± 1.90 mg/L vs. 37.42 ± 1.88 mmol/L) and intestinal PMCA1b gene expression than CaCl2 . These findings offer a new understanding of how sucrose-calcium chelate enhances intestinal calcium absorption and could be used as an ingredient in functional foods to treat calcium deficiency. PRACTICAL APPLICATION: The development of high-quality calcium supplements is crucial for addressing the various adverse symptoms associated with calcium deficiency. This study aimed to prepare a sucrose-calcium chelate and analyze its structure, as well as its potential to enhance calcium absorption in Caco-2 monolayer cells and mice. The results demonstrated that the sucrose-calcium chelate effectively promoted calcium absorption. Notably, its ability to enhance calcium transport was linked to its strong binding with PMCA1b, a calcium transporter located on the basement membrane, and its capacity to stimulate PMCA1b gene expression. These findings contribute to a deeper understanding of how the sucrose-calcium chelate enhances intestinal calcium absorption and suggest its potential use as an ingredient in functional foods for treating calcium deficiency.

LibGENiE - A bioinformatic pipeline for the design of information-enriched enzyme libraries

Enzymes are potent catalysts with high specificity and selectivity. To leverage nature's synthetic potential for industrial applications, various protein engineering techniques have emerged which allow to tailor the catalytic, biophysical, and molecular recognition properties of enzymes. However, the many possible ways a protein can be altered forces researchers to carefully balance between the exhaustiveness of an enzyme screening campaign and the required resources. Consequently, the optimal engineering strategy is often defined on a case-by-case basis. Strikingly, while predicting mutations that lead to an improved target function is challenging, here we show that the prediction and exclusion of deleterious mutations is a much more straightforward task as analyzed for an engineered carbonic acid anhydrase, a transaminase, a squalene-hopene cyclase and a Kemp eliminase. Combining such a pre-selection of allowed residues with advanced gene synthesis methods opens a path toward an efficient and generalizable library construction approach for protein engineering. To give researchers easy access to this methodology, we provide the website LibGENiE containing the bioinformatic tools for the library design workflow.

Insights from molecular dynamics simulations for computational protein design

A grand challenge in the field of structural biology is to design and engineer proteins that exhibit targeted functions. Although much success on this front has been achieved, design success rates remain low, an ever-present reminder of our limited understanding of the relationship between amino acid sequences and the structures they adopt. In addition to experimental techniques and rational design strategies, computational methods have been employed to aid in the design and engineering of proteins. Molecular dynamics (MD) is one such method that simulates the motions of proteins according to classical dynamics. Here, we review how insights into protein dynamics derived from MD simulations have influenced the design of proteins. One of the greatest strengths of MD is its capacity to reveal information beyond what is available in the static structures deposited in the Protein Data Bank. In this regard simulations can be used to directly guide protein design by providing atomistic details of the dynamic molecular interactions contributing to protein stability and function. MD simulations can also be used as a virtual screening tool to rank, select, identify, and assess potential designs. MD is uniquely poised to inform protein design efforts where the application requires realistic models of protein dynamics and atomic level descriptions of the relationship between dynamics and function. Here, we review cases where MD simulations was used to modulate protein stability and protein function by providing information regarding the conformation(s), conformational transitions, interactions, and dynamics that govern stability and function. In addition, we discuss cases where conformations from protein folding/unfolding simulations have been exploited for protein design, yielding novel outcomes that could not be obtained from static structures.

Electrostatics effects on Ca(2+) binding and conformational changes in EF-hand domains: Functional implications for EF-hand proteins

Mutations of Gln41 and Lys75 with nonpolar residues in the N-terminal domain of calmodulin (N-Cam) revealed the importance of solvation energetics in conformational change of Ca(2+) sensor EF-hand domains. While in general these domains have polar residues at these corresponding positions yet the extent of their conformational response to Ca(2+) binding and their Ca(2+) binding affinity can be different from N-Cam. Consequently, here we address the charge state of the polar residues at these positions. The results show that the charge state of these polar residues can affect substantially the conformational change and the Ca(2+) binding affinity of our N-Cam variants. Since all the variants kept their conformational activity in the presence of Ca(2+) suggests that the differences observed among them mainly originate from the difference in their molecular dynamics. Hence we propose that the molecular dynamics of Ca(2+) sensor EF-hand domains is a key factor in the multifunctional aspect of EF-hand proteins.

Two structural motifs within canonical EF-hand calcium-binding domains identify five different classes of calcium buffers and sensors

Proteins with EF-hand calcium-binding motifs are essential for many cellular processes, but are also associated with cancer, autism, cardiac arrhythmias, and Alzheimer's, skeletal muscle and neuronal diseases. Functionally, all EF-hand proteins are divided into two groups: (1) calcium sensors, which function to translate the signal to various responses; and (2) calcium buffers, which control the level of free Ca²⁺ ions in the cytoplasm. The borderline between the two groups is not clear, and many proteins cannot be described as definitive buffers or sensors. Here, we describe two highly-conserved structural motifs found in all known different families of the EF-hand proteins. The two motifs provide a supporting scaffold for the DxDxDG calcium binding loop and contribute to the hydrophobic core of the EF hand domain. The motifs allow more precise identification of calcium buffers and calcium sensors. Based on the characteristics of the two motifs, we could classify individual EF-hand domains into five groups: (1) Open static; (2) Closed static; (3) Local dynamic; (4) Dynamic; and (5) Local static EF-hand domains.

Functional and structural analysis of the conserved EFhd2 protein

EFhd2 is a novel protein conserved from C. elegans to H. sapiens. This novel protein was originally identified in cells of the immune and central nervous systems. However, it is most abundant in the central nervous system, where it has been found associated with pathological forms of the microtubule-associated protein tau. The physiological or pathological roles of EFhd2 are poorly understood. In this study, a functional and structural analysis was carried to characterize the molecular requirements for EFhd2's calcium binding activity. The results showed that mutations of a conserved aspartate on either EF-hand motif disrupted the calcium binding activity, indicating that these motifs work in pair as a functional calcium binding domain. Furthermore, characterization of an identified single-nucleotide polymorphisms (SNP) that introduced a missense mutation indicates the importance of a conserved phenylalanine on EFhd2 calcium binding activity. Structural analysis revealed that EFhd2 is predominantly composed of alpha helix and random coil structures and that this novel protein is thermostable. EFhd2's thermo stability depends on its N-terminus. In the absence of the N-terminus, calcium binding restored EFhd2's thermal stability. Overall, these studies contribute to our understanding on EFhd2 functional and structural properties, and introduce it into the family of canonical EF-hand domain containing proteins.

On the mechanism of protein fold-switching by a molecular sensor

Alternate frame folding (AFF) is a mechanism by which conformational change can be engineered into a protein. The protein structure switches from the wild-type fold (N) to a circularly-permuted fold (N'), or vice versa, in response to a signaling event such as ligand binding. Despite the fact that the two native states have similar structures, their interconversion involves folding and unfolding of large parts of the molecule. This rearrangement is reported by fluorescent groups whose relative proximities change as a result of the order-disorder transition. The nature of the conformational change is expected to be similar from protein to protein; thus, it may be possible to employ AFF as a general method to create optical biosensors. Toward that goal, we test basic aspects of the AFF mechanism using the AFF variant of calbindin D(9k). A simple three-state model for fold switching holds that N and N' interconvert through the unfolded state. This model predicts that the fundamental properties of the switch--calcium binding affinity, signal response (i.e., fluorescence change upon binding), and switching rate--can be controlled by altering the relative stabilities of N and N'. We find that selectively destabilizing N or N' changes the equilibrium properties of the switch (binding affinity and signal response) in accordance with the model. However, kinetic data indicate that the switching pathway does not require whole-molecule unfolding. The rate is instead limited by unfolding of a portion of the protein, possibly in concert with folding of a corresponding region.

Protein stability, flexibility and function

Proteins rely on flexibility to respond to environmental changes, ligand binding and chemical modifications. Potentially, a perturbation that changes the flexibility of a protein may interfere with its function. Millions of mutations have been performed on thousands of proteins in quests for a delineation of the molecular details of their function. Several of these mutations interfered with the binding of a specific ligand with a concomitant effect on the stability of the protein scaffold. It has been ambiguous and not straightforward to recognize if any relationships exist between the stability of a protein and the affinity for its ligand. In this review, we present examples of proteins where changes in stability results in changes in affinity and of proteins where stability and affinity are uncorrelated. We discuss the possibility for a relationship between stability and binding. From the data presented is it clear that there are specific sites (flexibility hotspots) in proteins that are important for both binding and stability. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

A Ca2+-sensing molecular switch based on alternate frame protein folding

Existing strategies for creating biosensors mainly rely on large conformational changes to transduce a binding event to an output signal. Most molecules, however, do not exhibit large-scale structural changes upon substrate binding. Here, we present a general approach (alternate frame folding, or AFF) for engineering allosteric control into ligand binding proteins. AFF can in principle be applied to any protein to establish a binding-induced conformational change, even if none exists in the natural molecule. The AFF design duplicates a portion of the amino acid sequence, creating an additional "frame" of folding. One frame corresponds to the wild-type sequence, and folding produces the normal structure. Folding in the second frame yields a circularly permuted protein. Because the two native structures compete for a shared sequence, they fold in a mutually exclusive fashion. Binding energy is used to drive the conformational change from one fold to the other. We demonstrate the approach by converting the protein calbindin D(9k) into a molecular switch that senses Ca2+. The structures of Ca2+-free and Ca2+-bound calbindin are nearly identical. Nevertheless, the AFF mechanism engineers a robust conformational change that we detect using two covalently attached fluorescent groups. Biological fluorophores can also be employed to create a genetically encoded sensor. AFF should be broadly applicable to create sensors for a variety of small molecules.

Structures and metal-ion-binding properties of the Ca2+-binding helix–loop–helix EF-hand motifs

The ‘EF-hand’ Ca2+-binding motif plays an essential role in eukaryotic cellular signalling, and the proteins containing this motif constitute a large and functionally diverse family. The EF-hand is defined by its helix–loop–helix secondary structure as well as the ligands presented by the loop to bind the Ca2+ ion. The identity of these ligands is semi-conserved in the most common (the ‘canonical’) EF-hand; however, several non-canonical EF-hands exist that bind Ca2+ by a different co-ordination mechanism. EF-hands tend to occur in pairs, which form a discrete domain so that most family members have two, four or six EF-hands. This pairing also enables communication, and many EF-hands display positive co-operativity, thereby minimizing the Ca2+ signal required to reach protein saturation. The conformational effects of Ca2+ binding are varied, function-dependent and, in some cases, minimal, but can lead to the creation of a protein target interaction site or structure formation from a molten-globule apo state. EF-hand proteins exhibit various sensitivities to Ca2+, reflecting the intrinsic binding ability of the EF-hand as well as the degree of co-operativity in Ca2+ binding to paired EF-hands. Two additional factors can influence the ability of an EF-hand to bind Ca2+: selectivity over Mg2+ (a cation with very similar chemical properties to Ca2+ and with a cytoplasmic concentration several orders of magnitude higher) and interaction with a protein target. A structural approach is used in this review to examine the diversity of family members, and a biophysical perspective provides insight into the ability of the EF-hand motif to bind Ca2+ with a wide range of affinities.

Effects of Calcium Binding on the Side-chain Methyl Dynamics of Calbindin D9k: A 2H NMR Relaxation Study

The effects of Ca(2+) binding on the side-chain methyl dynamics of calbindin D(9k) have been characterized by (2)H NMR relaxation rate measurements. Longitudinal, transverse in-phase, quadrupolar order, transverse anti-phase and double quantum relaxation rates are reported for both the apo and Ca(2+)-loaded states of the protein at two magnetic field strengths. The relatively large size of the data set allows for a detailed analysis of the underlying conformational dynamics by spectral density mapping and model-free fitting procedures. The results reveal a correlation between a methyl group's distance from the Ca(2+) binding sites and its conformational dynamics. Several methyl groups segregate into two limiting classes, one proximal and the other distal to the binding sites. Methyl groups in these two classes respond differently to Ca(2+) binding, both in terms of the timescale and amplitude of their fluctuations. Ca(2+) binding elicits a partial immobilization among methyl groups in the proximal class, which is consistent with previous studies of calbindin's backbone dynamics. The distal class, however, exhibits a trend that could not be inferred from the backbone data in that its mobility actually increases with Ca(2+) binding. We have introduced the term polar dynamics to describe this type of organization across the molecule. The trend may represent an important mechanism by which calbindin D(9k) achieves high affinity binding while minimizing the corresponding loss of conformational entropy.

Designing sequence to control protein function in an EF-hand protein

The extent of conformational change that calcium binding induces in EF-hand proteins is a key biochemical property specifying Ca(2+) sensor versus signal modulator function. To understand how differences in amino acid sequence lead to differences in the response to Ca(2+) binding, comparative analyses of sequence and structures, combined with model building, were used to develop hypotheses about which amino acid residues control Ca(2+)-induced conformational changes. These results were used to generate a first design of calbindomodulin (CBM-1), a calbindin D(9k) re-engineered with 15 mutations to respond to Ca(2+) binding with a conformational change similar to that of calmodulin. The gene for CBM-1 was synthesized, and the protein was expressed and purified. Remarkably, this protein did not exhibit any non-native-like molten globule properties despite the large number of mutations and the nonconservative nature of some of them. Ca(2+)-induced changes in CD intensity and in the binding of the hydrophobic probe, ANS, implied that CBM-1 does undergo Ca(2+) sensorlike conformational changes. The X-ray crystal structure of Ca(2+)-CBM-1 determined at 1.44 A resolution reveals the anticipated increase in hydrophobic surface area relative to the wild-type protein. A nascent calmodulin-like hydrophobic docking surface was also found, though it is occluded by the inter-EF-hand loop. The results from this first calbindomodulin design are discussed in terms of progress toward understanding the relationships between amino acid sequence, protein structure, and protein function for EF-hand CaBPs, as well as the additional mutations for the next CBM design.

Designing calcium-sensitizing mutations in the regulatory domain of cardiac troponin C

Cardiac troponin C belongs to the EF-hand superfamily of calcium-binding proteins and plays an essential role in the regulation of muscle contraction and relaxation. To follow calcium binding and exchange with the regulatory N-terminal domain (N-domain) of human cardiac troponin C, we substituted Phe at position 27 with Trp, making a fluorescent cardiac troponin C(F27W). Trp(27) accurately reported the kinetics of calcium association and dissociation of the N-domain of cardiac troponin C(F27W). To sensitize the N-domain of cardiac troponin C(F27W) to calcium, we individually substituted the hydrophobic residues Phe(20), Val(44), Met(45), Leu(48), and Met(81) with polar Gln. These mutations were designed to increase the calcium affinity of the N-domain of cardiac troponin C by facilitating the movement of helices B and C (BC unit) away from helices N, A, and D (NAD unit). As anticipated, these selected hydrophobic residue substitutions increased the calcium affinity of the regulatory domain of cardiac troponin C(F27W) approximately 2.1-15.2-fold. Surprisingly, the increased calcium affinity caused by the hydrophobic residue substitutions was largely due to faster calcium association rates (2.6-8.7-fold faster) rather than to slower calcium dissociation rates (1.2-2.9-fold slower). The regulatory N-domains of cardiac troponin C(F27W) and its mutants were also able to bind magnesium competitively and with physiologically relevant affinities (1.2-2.7 mm). The design of calcium-sensitizing cardiac troponin C mutants presented in this work enhances the understanding of how to control cation binding properties of EF-hand proteins and ultimately their structure and physiological function.

Correlation between normal modes in the 20-200 cm-1 frequency range and localized torsion motions related to certain collective motions in proteins

In certain biologically relevant collective motions, such as protein domain motions and sub-domain motions, large amplitude movements are localized in one or a few flexible regions consisting of a small number of residues. This paper explores the possible use of normal mode analysis in probing localized vibrational torsion motions in these flexible regions that may be related to certain collective motions. The normal modes of 10 structures of five proteins in different conformation (TRP repressor, calmodulin, calbindin D(9k), HIV-1 protease and troponin C), known to have shear or hinge domain or sub-domain motion, respectively, are analyzed. Our study identifies, for each structure, unique normal modes in the 20-200 cm-1 frequency range, whose corresponding motions are primarily concentrated in the region where large amplitude torsion movements of a known domain or sub-domain motion occur. This suggests possible correlation between normal modes at 20-200 cm-1 frequency range and initial fluctuational motions leading to localized collective motions in proteins, and thus the potential application of normal mode analysis in facilitating the study of biologically important localized motions in biomolecules.

The EF-hand domain: a globally cooperative structural unit

EF-hand Ca(2+)-binding proteins participate in both modulation of Ca(2+) signals and direct transduction of the ionic signal into downstream biochemical events. The range of biochemical functions of these proteins is correlated with differences in the way in which they respond to the binding of Ca(2+). The EF-hand domains of calbindin D(9k) and calmodulin are homologous, yet they respond to the binding of calcium ions in a drastically different manner. A series of comparative analyses of their structures enabled the development of hypotheses about which residues in these proteins control the calcium-induced changes in conformation. To test our understanding of the relationship between protein sequence and structure, we specifically designed the F36G mutation of the EF-hand protein calbindin D(9k) to alter the packing of helices I and II in the apoprotein. The three-dimensional structure of apo F36G was determined in solution by nuclear magnetic resonance spectroscopy and showed that the design was successful. Surprisingly, significant structural perturbations also were found to extend far from the site of mutation. The observation of such long-range effects provides clear evidence that four-helix EF-hand domains should be treated as a single globally cooperative unit. A hypothetical mechanism for how the long-range effects are transmitted is described. Our results support the concept of energetic and structural coupling of the key residues that are crucial for a protein's fold and function.

Structure of rat parvalbumin with deleted AB domain: implications for the evolution of EF hand calcium-binding proteins and possible physiological relevance

Among the EF-hand Ca(2+)-binding proteins, parvalbumin (PV) and calbindin D9k (CaB) have the function of Ca(2+) buffers. They evolved from an ancestor protein through two phylogenetic pathways, keeping one pair of EF-hands. They differ by the extra helix-loop-helix (AB domain) found in PV and by the linker between the binding sites. To investigate whether the deletion of AB in PV restores a CaB-like structure, we prepared and solved the structure of the truncated rat PV (PVratDelta37) by X-ray and NMR. PVratDelta37 keeps the PV fold, but is more compact, having a well-structured linker, which differs remarkably from CaB. PvratDelta37 has no stable apo-form, has lower affinity for Ca(2+) than full-length PV, and does not bind Mg(2+), in contrast to CaB. Structural differences of the hydrophobic core are partially responsible for lowering the calcium-binding affinity of the truncated protein. It can be concluded that the AB domain, like the linker of CaB, plays a role in structural stabilization. The AB domain of PV protects the hydrophobic core, and is required to maintain high affinity for divalent cation binding. Therefore, the AB domain possibly modulates PV buffer function.

The dimerization interface of the metastasis-associated protein S100A4 (Mts1): in vivo and in vitro studies

The S100 calcium-binding proteins are implicated in signal transduction, motility, and cytoskeletal dynamics. The three-dimensional structure of several S100 proteins revealed that the proteins form non-covalent dimers. However, the mechanism of the S100 dimerization is still obscure. In this study we characterized the dimerization of S100A4 (also named Mts1) in vitro and in vivo. Analytical ultracentrifugation revealed that apoS100A4 was present in solution as a mixture of monomers and dimers in a rapidly reversible equilibrium (K(d) = 4 +/- 2 microm). The binding of calcium promoted dimerization. Replacement of Tyr-75 by Phe resulted in the stabilization of the dimer. Helix IV is known to form the major part of the dimerization interface in homologous S100 proteins. By using the yeast two-hybrid system we showed that only a few residues of helix IV, namely Phe-72, Tyr-75, Phe-78, and Leu-79, are essential for dimerization in vivo. A homology model demonstrated that these residues form a hydrophobic cluster on helix IV. Their role is to stabilize the structure of individual subunits rather than provide specific interactions across the dimerization surface. Our mutation data showed that the specificity at the dimerization surface is not particularly stringent, which is consistent with recent data indicating that S100 proteins can form heterodimers.

Role of calcium ions in the structure and function of the di-isopropylfluorophosphatase from Loligo vulgaris

Di-isopropylfluorophosphatase (DFPase) is shown to contain two high-affinity Ca(2+)-binding sites, which are required for catalytic activity and stability. Incubation with chelating agents results in the irreversible inactivation of DFPase. From titrations with Quin 2 [2-([2-[bis(carboxymethyl)amino]-5-methylphenoxy]-methyl)-6-methoxy-8-[bis(carboxymethyl)-amino]quinoline], a lower-affinity site with dissociation constants of 21 and 840 nM in the absence and the presence of 150 mM KCl respectively was calculated. The higher-affinity site was not accessible, indicating a dissociation constant of less than 5.3 nM. Stopped-flow experiments have shown that the dissociation of bound Ca(2+) occurs in two phases, with rates of approx. 1.1 and 0.026 s(-1) corresponding to the dissociation from the low-affinity and high-affinity sites respectively. Dissociation rates depend strongly on temperature but not on ionic strength, indicating that Ca(2+) dissociation is connected with conformational changes. Limited proteolysis, CD spectroscopy, dynamic light scattering and the binding of 8-anilino-1-naphthalenesulphonic acid have been combined to give a detailed picture of the conformational changes induced on the removal of Ca(2+) from DFPase. The Ca(2+) dissociation is shown to result in a primary, at least partly reversible, step characterized by a large decrease in DFPase activity and some changes in enzyme structure and shape. This step is followed by an irreversible denaturation and aggregation of the apo-enzyme. From the temperature dependence of Ca(2+) dissociation and the denaturation results we conclude that the higher-affinity Ca(2+) site is required for stabilizing DFPase's structure, whereas the lower-affinity site is likely to fulfil a catalytic function.