Solution structure of (Cd2+)1-calbindin D9k reveals details of the stepwise structural changes along the Apo-->(Ca2+)II1-->(Ca2+)I,II2 binding pathway

Calcium Activation of the Androgen Receptor in Prostate Cells

Background

Although prostate cancer patients initially respond to androgen deprivation therapy, most patients progress to a resistant phenotype. Castration resistance is due, in part, to intratumoral and/or adrenal synthesis of androgens, overexpression or mutation of the androgen receptor (AR), stabilization of AR by chaperones, and ligand-independent activation of AR. Increasing evidence also links disruption of calcium homeostasis to progression of prostate cancer. Our previous study shows that heavy metal cadmium activates the AR through a ligand-independent mechanism. Cadmium mimics calcium in biological systems due to their similar ionic charge and radius. This study determines whether calcium activates AR and whether first- and second-generation antiandrogens block the ability of calcium to activate the receptor.

Methods

The expression of androgen-responsive genes and calcium channels was measured in prostate cells using a quantitative real-time polymerase chain reaction assay. Cell growth was measured.

Results

To ask whether calcium activates AR, prostate cells were treated with calcium in the absence and presence of the first-generation antiandrogens hydroxyflutamide and bicalutamide and the second-generation antiandrogen enzalutamide, and the expression of androgen-responsive genes and cell growth was measured. In the normal PWR-1E cells and HEK293T cells transiently expressing AR, treatment with calcium increased the expression of androgen-responsive genes by approximately 3-fold. The increase was blocked by enzalutamide but was not consistently blocked by the first-generation antiandrogens. In LNCaP cells which contain a mutant AR, treatment with calcium also increased the expression of androgen-responsive genes by approximately 3-fold, and the increase was more effectively blocked by enzalutamide than by hydroxyflutamide or bicalutamide. Treatment with calcium also increased cell growth that was blocked by enzalutamide. To ask whether dysregulation of calcium channels is associated with castration resistance, calcium channels were measured in the normal PWR-1E prostate cells, the hormone-responsive LNCaP cells, and the castration-resistant VCaP and 22RV1 cells. Compared to normal prostate cells, the hormone-responsive and hormone-resistant cells overexpressed several calcium channels.

Conclusions

The results of this study show that calcium activates AR and increases cell growth and that calcium channels are overexpressed in hormone-responsive and hormone-resistant prostate cancer cells. Taken together, the results suggest a novel role of calcium in the castration-resistant phenotype.

Mouse S100G protein exhibits properties characteristic of a calcium sensor

Bovine S100 G (calbindin D_9k, small Ca²⁺-binding protein of the EF-hand superfamily) is considered as a calcium buffer protein; i.e., the binding of Ca²⁺ practically does not change its general conformation. A set of experimental approaches has been used to study structural properties of apo- and Ca²⁺-loaded forms of mouse S100 G (81.4% identity in amino acid sequence with bovine S100 G). This analysis revealed that, in contrast to bovine S100 G, the removal of calcium ions increases α-helices content of mouse S100 G protein and enhances its accessibility to digestion by α-chymotrypsin. Furthermore, mouse apo-S100 G is characterized by a decreased surface hydrophobicity and reduced tendency for oligomerization. Such behavior is typical of calcium sensor proteins. Apo-state of mouse S100 G still has rather compact structure, which can be cooperatively unfolded by temperature and GdnHCl. Computational analysis of amino acid sequences of S100 G proteins shows that these proteins could be in a disordered state upon a removal of the bound calcium ions. The experimental data show that, although mouse apo-S100 G is flexible compared to the Ca²⁺-loaded state, the apo-form is not completely disordered and preserves some cooperatively meting structure. The origin of the unexpectedly high stability of mouse S100 G can be rationalized by an exceptionally strong association of its N- and C-terminal parts containing the EF-hands I and II, respectively.

Effect of monovalent ion binding on molecular dynamics of the S100-family calcium-binding protein calbindin D9k

Calbindin D_9k is a member of the S100 subfamily of EF-hand calcium binding proteins, and has served as an important model system for biophysical studies. The fast timescale dynamics of the calcium-free (apo) state is characterized using molecular dynamics simulations. Order parameters for the backbone NH bond vectors are determined from the simulations and compared with experimentally derived values, with a focus on the dynamics of calcium-binding site I. There is a significant discrepancy between simulated and experimental order parameters for site I residues in the case of no ion bound in site I. However, it was found in the simulations that a Na⁺ ion can bind in site I, and the resulting order parameters determined from the simulations are in excellent agreement with experiment. Comparisons are made to X-ray structures of other S100 family members in which Na⁺ ions were observed or suggested to be bound in site I. © 2019 Wiley Periodicals, Inc.

Structural Insights into the Calcium-Mediated Allosteric Transition in the C-Terminal Domain of Calmodulin from Nuclear Magnetic Resonance Measurements

Calmodulin is a two-domain signaling protein that becomes activated upon binding cooperatively two pairs of calcium ions, leading to large-scale conformational changes that expose its binding site. Despite significant advances in understanding the structural biology of calmodulin functions, the mechanistic details of the conformational transition between closed and open states have remained unclear. To investigate this transition, we used a combination of molecular dynamics simulations and nuclear magnetic resonance (NMR) experiments on the Ca(2+)-saturated E140Q C-terminal domain variant. Using chemical shift restraints in replica-averaged metadynamics simulations, we obtained a high-resolution structural ensemble consisting of two conformational states and validated such an ensemble against three independent experimental data sets, namely, interproton nuclear Overhauser enhancements, (15)N order parameters, and chemical shift differences between the exchanging states. Through a detailed analysis of this structural ensemble and of the corresponding statistical weights, we characterized a calcium-mediated conformational transition whereby the coordination of Ca(2+) by just one oxygen of the bidentate ligand E140 triggers a concerted movement of the two EF-hands that exposes the target binding site. This analysis provides atomistic insights into a possible Ca(2+)-mediated activation mechanism of calmodulin that cannot be achieved from static structures alone or from ensemble NMR measurements of the transition between conformations.

Analysis of common bean (Phaseolus vulgaris L., genotype BAT93) calmodulin cDNA using computational tools

BACKGROUND

Common bean (Phaseolus vulgaris L.) is an important part of the human diet and serves as a source of natural products. Identification and understanding of genes in P. vulgaris is important for its improvement. Characterization of expressed sequence tags (ESTs) is one of the approaches in understanding the expressed genes. For the understanding of genes expression in P. vulgaris pod-tissue, research work of ESTs generation was initiated by constructing cDNA libraries using 5-day and 20-day old bean-pod-tissues. Altogether, 5972 cDNA clones were isolated to have ESTs. While processing ESTs, we found a transcript for calmodulin (CaM) gene. It is an important gene that encodes for a calcium-binding protein and known to express in all eukaryotic cells. Hence, this study was undertaken to analyse and annotate it.

OBJECTIVE

The objective of this study was to analyze and annotate P. vulgaris CaM (PvCaM) gene cDNA and its deduced protein (amino acids) sequence.

MATERIALS AND METHODS

Both strands of PvCaM cDNA clone were sequenced using M13 forward and reverse primer to elucidate the nucleotide sequence. The cDNA sequence and deduced protein sequence were analyzed and annotated using bioinformatics tools available online. The secondary structures and three-dimensional (3D) structure of PvCaM protein were predicted using the Phyre automatic fold recognition server.

RESULTS

Results showed that PvCaM cDNA is 818 bp in length. The cDNA analysis results showed that it contains an open reading frame that encodes for 149 amino acid residues. The deduced protein sequence analysis results showed the presence of conserved domains required for CaM function. The predicted secondary structures and 3D structure are analogous to the Solanum tuberosum CaM.

CONCLUSIONS

This study analyzed and annotated PvCaM cDNA and protein. However, in order to obtain a complete understanding of PvCaM protein, further study on its expression, structure and regulation is essential.

Myoplasmic resting Ca2+ regulation by ryanodine receptors is under the control of a novel Ca2+-binding region of the receptor

Passive SR (sarcoplasmic reticulum) Ca²⁺ leak through the RyR (ryanodine receptor) plays a critical role in the mechanisms that regulate [Ca²⁺]_rest (intracellular resting myoplasmic free Ca²⁺ concentration) in muscle. This process appears to be isoform-specific as expression of either RyR1 or RyR3 confers on myotubes different [Ca²⁺]_rest. Using chimaeric RyR3–RyR1 receptors expressed in dyspedic myotubes, we show that isoform-dependent regulation of [Ca²⁺]_rest is primarily defined by a small region of the receptor encompassing amino acids 3770–4007 of RyR1 (amino acids 3620–3859 of RyR3) named as the CLR (Ca²⁺ leak regulatory) region. [Ca²⁺]_rest regulation by the CLR region was associated with alteration of RyRs’ Ca²⁺-activation profile and changes in SR Ca²⁺-leak rates. Biochemical analysis using Tb³⁺-binding assays and intrinsic tryptophan fluorescence spectroscopy of purified CLR domains revealed that this determinant of RyRs holds a novel Ca²⁺-binding domain with conformational properties that are distinctive to each isoform. Our data suggest that the CLR region provides channels with unique functional properties that modulate the rate of passive SR Ca²⁺ leak and confer on RyR1 and RyR3 distinctive [Ca²⁺]_rest regulatory properties. The identification of a new Ca²⁺-binding domain of RyRs with a key modulatory role in [Ca²⁺]_rest regulation provides new insights into Ca²⁺-mediated regulation of RyRs.

This paper reports the finding of a new class of Ca²⁺-binding domain of intracellular Ca²⁺ channels from muscle cells. This domain provides channels with distinctive properties that result in channel-specific modulation of the intracellular resting Ca²⁺ concentration.

Automated structure refinement for a protein heterodimer complex using limited EPR spectroscopic data and a rigid-body docking algorithm: a three-dimensional model for an ankyrin-CDB3 complex

We report here specialized functions incorporated recently in the rigid-body docking software toolkit TagDock to utilize electron paramagnetic resonance derived (EPR-derived) interresidue distance measurements and spin-label accessibility data. The TagDock package extensions include a custom methanethiosulfonate spin label rotamer library to enable explicit, all-atom spin-label side-chain modeling and scripts to evaluate spin-label surface accessibility. These software enhancements enable us to better utilize the biophysical data routinely available from various spin-labeling experiments. To illustrate the power and utility of these tools, we report the refinement of an ankyrin:CDB3 complex model that exhibits much improved agreement with the EPR distance measurements, compared to model structures published previously.

TagDock: an efficient rigid body docking algorithm for oligomeric protein complex model construction and experiment planning

We report here new computational tools and strategies to efficiently generate three-dimensional models for oligomeric biomolecular complexes in cases where there is limited experimental restraint data to guide the docking calculations. Our computational tools are designed to rapidly and exhaustively enumerate all geometrically possible docking poses for an oligomeric complex, rather than generate detailed, atomic-resolution models. Experimental data, such as interatomic distance measurements, are then used to select and refine docking poses that are consistent with the experimental restraints. Our computational toolkit is designed for use with sparse data sets to generate intermediate-resolution docking models, and utilizes distance difference matrix analysis to identify further restraint measurements that will provide maximum additional structural refinement. Thus, these tools can be used to help plan optimal residue positions for probe incorporation in labor-intensive biophysical experiments such as chemical cross-linking, electron paramagnetic resonance, or Förster resonance energy transfer spectroscopy studies. We present benchmark results for docking the collection of all 176 heterodimer protein complexes from the ZDOCK database, as well as a protein homodimer with recently collected experimental distance restraints, to illustrate the toolkit's capabilities and performance, and to demonstrate how distance difference matrix analysis can automatically identify and prioritize additional restraint measurements that allow us to rapidly optimize docking poses.

Metals and breast cancer

Metalloestrogens are metals that activate the estrogen receptor in the absence of estradiol. The metalloestrogens fall into two subclasses: metal/metalloid anions and bivalent cationic metals. The metal/metalloid anions include compounds such as arsenite, nitrite, selenite, and vanadate while the bivalent cations include metals such as cadmium, calcium, cobalt, copper, nickel, chromium, lead, mercury, and tin. The best studied metalloestrogen is cadmium. It is a heavy metal and a prevalent environmental contaminant with no known physiological function. This review addresses our current understanding of the mechanism by which cadmium and the bivalent cationic metals activate estrogen receptor-α. The review also summarizes the in vitro and in vivo evidence that cadmium functions as an estrogen and the potential role of cadmium in breast cancer.

Crystal structure of cardiac troponin C regulatory domain in complex with cadmium and deoxycholic acid reveals novel conformation

The amino-terminal regulatory domain of cardiac troponin C (cNTnC) plays an important role as the calcium sensor for the troponin complex. Calcium binding to cNTnC results in conformational changes that trigger a cascade of events that lead to cardiac muscle contraction. The cardiac N-terminal domain of TnC consists of two EF-hand calcium binding motifs, one of which is dysfunctional in binding calcium. Nevertheless, the defunct EF-hand still maintains a role in cNTnC function. For its structural analysis by X-ray crystallography, human cNTnC with the wild-type primary sequence was crystallized under a novel crystallization condition. The crystal structure was solved by the single-wavelength anomalous dispersion method and refined to 2.2 Å resolution. The structure displays several novel features. Firstly, both EF-hand motifs coordinate cadmium ions derived from the crystallization milieu. Secondly, the ion coordination in the defunct EF-hand motif accompanies unusual changes in the protein conformation. Thirdly, deoxycholic acid, also derived from the crystallization milieu, is bound in the central hydrophobic cavity. This is reminiscent of the interactions observed for cardiac calcium sensitizer drugs that bind to the same core region and maintain the "open" conformational state of calcium-bound cNTnC. The cadmium ion coordination in the defunct EF-hand indicates that this vestigial calcium binding site retains the structural and functional elements that allow it to coordinate a cadmium ion. However, it is a result of, or concomitant with, large and unusual structural changes in cNTnC.

Relating form and function of EF-hand calcium binding proteins

The EF hand, a helix-loop-helix structure, is one of the most common motifs found in animal genomes, and EF-hand Ca(2+)-binding proteins (EFCaBPs) are widely distributed throughout the cell. However, researchers remain confounded by a lack of understanding of how peptide sequences code for specific functions and by uncertainty about the molecular mechanisms that enable EFCaBPs to distinguish among many diverse cellular targets. Such knowledge could define the roles of EFCaBPs in health and disease and ultimately enable control or even design of Ca(2+)-dependent functions in medicine and biotechnology. In this Account, we describe our structural and biochemical research designed to understand the sequence-to-function relationship in EFCaBPs. The first structural goal was to define conformational changes induced by binding Ca(2+), and our group and others established that solution NMR spectroscopy is well suited for this task. We pinpointed residues critical to the differences in Ca(2+) response of calbindin D(9k) and calmodulin (CaM), homologous EFCaBPs from different functional classes, by using direct structure determination with site-directed mutagenesis and protein engineering. Structure combined with biochemistry provided the foundation for identifying the fundamental mechanism of cooperativity in the binding of Ca(2+) ions: this cooperativity provides EFCaBPs with the ability to detect the relatively small changes in concentration that constitute Ca(2+) signals. Using calbindin D(9k) as a model system, studies of the structure and fast time scale dynamics of each of the four ion binding states in a typical EF-hand domain provided direct evidence that site-site communication lowers the free energy cost of reorganization for binding the second ion. Our work has also extended models of how EFCaBPs interact with their cellular targets. We determined the unique dimeric architecture of S100 proteins, a specialized subfamily of EFCaBPs found exclusively in vertebrates. We described the implications for how these proteins transduce signals and went on to characterize interactions with peptide fragments of important cellular targets. Studies of the CaM homolog centrin revealed novel characteristics of its binding of Ca(2+) and its interaction with its cellular target Kar1. These results provided clear examples of how subtle differences in sequence fine-tune EFCaBPs to interact with their specific targets. The structural approach stands at a critical crossroad, shifting in emphasis from descriptive structural biochemistry to integrated biology and medicine. We present our dual-molecular-switch model for Ca(2+) regulation of gating functions of voltage-gated sodium channels in which both CaM and an intrinsic EF-hand domain serve as coupled Ca(2+) sensors. A second example involves novel EFCaBP extracellular function, that is, the role of S100A8/S100A9 heterodimer in the innate immune response to bacterial pathogens. A mechanism for the antimicrobial activity of S100A8/S100A9 was discovered. We describe interactions of S100A8/S100A9 and S100B with the cell surface receptor for advanced glycation end products. Biochemical and structural studies are now uncovering the mechanisms by which EFCaBPs work and are helping to define their biological activities, while simultaneously expanding knowledge of the roles of these proteins in normal cellular physiology and the pathology of disease.

The role of calcium in the activation of estrogen receptor-alpha

Environmental estrogen mimics, including metalloestrogens that can activate estrogen receptor-alpha (ERα), may contribute to breast cancer risk. However, the underlying mechanisms through which these molecular mimics activate the ERα are generally poorly understood. With concern to this important question, we investigated whether intracellular calcium may mediate the cross-talk between signaling pathways that activate ERα and the ligand-binding domain of ERα. MCF-7 cells treated with EGF, ATP, extracellular calcium, or caffeine to increase intracellular calcium triggered a rapid recruitment of ERα to estrogen-responsive promoters and stimulated expression of estrogen-responsive genes including pS2, complement C3, and progesterone receptor. Induction was blocked by an antiestrogen but also by the chelation of intracellular calcium. Treatment with extracellular calcium also increased the growth of MCF-7 cells through an ER-dependent mechanism. We found that EGF and extracellular calcium activated the C-terminus of ERα and the activation was blocked by the antiestrogen. Mechanistic investigations identified four potential sites on the solvent-accessible surface of the ERα ligand-binding domain as important for calcium activation of the receptor. Taken together, our results suggest that calcium mediates the cross-talk between ERα-activating signaling pathways and the ligand-binding domain of ERα providing a potential explanation for the ability of certain environmental metalloestrogens to activate the receptor.

Binding of calcium and other metal ions to the EF-hand loops of calmodulin studied by quantum chemical calculations and molecular dynamics simulations

Calcium ion binding by the four EF-hand motifs of the protein calmodulin (CaM) is a central event in Ca2+-based cellular signaling. To understand molecular details of this complex process, isolated Ca2+-binding loops can be studied, by use of both experiments and calculations. In this work, we explore the metal specificity of the four Ca2+-binding loops of CaM using density functional theory (DFT) quantum chemical calculations and molecular dynamics simulations. We study CaM complexes with the physiologically important ions of calcium (Ca2+) and magnesium (Mg2+) and also with two other ions, strontium (Sr2+) and lanthanum (La3+). The former is of interest in the area of radioactive waste bioremediation, whereas the latter is often used as a probe of Ca2+-binding sites. We obtain intrinsic metal ion-loop binding energies as well as their components: vacuum, charge-transfer, solvation, entropy, and deformation terms. A detailed analysis of the results reveals that the total binding energy depends on a delicate balance among these energy components. They, in turn, are determined by the cation's charge and size as well as the amino acid composition and flexibility of the loops and the identity of the metal-chelating residues.

A new methodology for the simulation of flexible protein–ligand interactions

A methodology has been developed for the simulation of induced fit between a ligand and its target protein. It utilizes constrained molecular dynamics where atoms determined to be immobile from difference distance matrix studies are fixed. Application of this methodology to HIV-1 reverse transcriptase (RT) as the example target protein has demonstrated its robustness. Short simulation times are sufficient to achieve good refinement of docking poses resulting from exchange of structurally dissimilar inhibitors between crystal structures.

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.

Identification and dissection of Ca(2+)-binding sites in the extracellular domain of Ca(2+)-sensing receptor

Ca(2+)-sensing receptors (CaSRs) represent a class of receptors that respond to changes in the extracellular Ca(2+) concentration ([Ca(2+)](o)) and activate multiple signaling pathways. A major barrier to advancing our understanding of the role of Ca(2+) in regulating CaSRs is the lack of adequate information about their Ca(2+)-binding locations, which is largely hindered by the lack of a solved three-dimensional structure and rapid off rates due to low Ca(2+)-binding affinities. In this paper, we have reported the identification of three potential Ca(2+)-binding sites in a modeled CaSR structure using computational algorithms based on the geometric description and surface electrostatic potentials. Mutation of the predicted ligand residues in the full-length CaSR caused abnormal responses to [Ca(2+)](o), similar to those observed with naturally occurring activating or inactivating mutations of the CaR, supporting the essential role of these predicted Ca(2+)-binding sites in the sensing capability of the CaSR. In addition, to probe the intrinsic Ca(2+)-binding properties of the predicted sequences, we engineered two predicted continuous Ca(2+)-binding sequences individually into a scaffold protein provided by a non-Ca(2+)-binding protein, CD2. We report herein the estimation of the metal-binding affinities of these predicted sites in the CaSR by monitoring aromatic-sensitized Tb(3+) fluorescence energy transfer. Removing the predicted Ca(2+)-binding ligands resulted in the loss of or significantly weakened cation binding. The potential Ca(2+)-binding residues were shown to be involved in Ca(2+)/Ln(3+) binding by high resolution NMR and site-directed mutagenesis, further validating our prediction of Ca(2+)-binding sites within the extracellular domain of the CaSR.

"Four-dimensional" protein structures: examples from metalloproteins

The fact that an object, for example, a protein, possesses a three-dimensional structure seems an obvious concept. However, when the object is flexible, the concept is less obvious. Growing experimental data over several decades show that proteins are not rigid objects, but they may sample more or less wide ranges of different conformations. To stress this concept, we propose to call the range of sampled conformations the "fourth dimension" of the protein structure. Nuclear magnetic resonance is a precious technique to define this fourth dimension. Examples of conformational heterogeneity taken from the realm of metalloproteins and their functional implications are discussed.

Assignment strategy for fast relaxing signals: complete aminoacid identification in thulium substituted calbindin D 9K

Paramagnetic proteins generally contain regions with diverse relaxation properties. Nuclei in regions far from the metal center may behave like those in diamagnetic proteins, but those closer to the metal experience rapid relaxation with accompanying line broadening. We have used a set of NMR experiments optimized to capture data from these various concentric regions in assigning the signals from a paramagnetic Calbindin D 9K derivative in which one of the two calcium ions has been replaced by thulium(III). Normal double- and triple-resonance experiments with 1H detection were used in collecting data from nuclei in the diamagnetic-like region; these approaches identified signals from fewer than 50% of the amino acid residues (those with d > 17.5 A from thulium(III)). Paramagnetism-optimized two-dimensional NMR experiments with 1H detection were used in collecting data from nuclei in the next nearer region (d > 15 A). Standard (d > 14 A) and optimized (d > 9 A) 13C direct-detection experiments were used to capture data from nuclei in the next layer. Finally nuclei closest to the metal were detected by one-dimensional 13C (d > 5 A) and one-dimensional 15N data collection (d > 4.2 A). NMR signals were assigned on the basis of through-bond correlations and, for signals closest to the metal, pseudocontact shifts. The latter were determined from chemical shift differences between assigned signals in thulium(III) and lanthanum(III) derivatives of Calbindin D 9K and they were interpreted on the basis of a structural model for the lanthanide-substituted protein. This approach yielded assignments of at least one resonance per amino acid residue, including those in the thulium(III) coordination sphere.

Principal Component Analysis of the Conformational Freedom within the EF-Hand Superfamily

A database of nonredundant structures of EF-hand domains--i.e., pairs of helix-loop-helix motifs--has been assembled, and the six angles among the four helices re-determined. A principal component analysis of these angles allows us to use two such components (PC1 and PC2) to describe the system retaining 80% of the total variance. A PC2 against PC1 plot representation allows us to represent in a compact way the full range of structural diversity of EF-hand domains, their grouping into protein families, and the variation for each family upon calcium and peptide binding.

Calcium-modulated S100 protein-phospholipid interactions. An NMR study of calbindin D9k and DPC

The cellular functions of several S100 proteins involve specific interactions with phospholipids and the cell membrane. The interactions between calbindin D(9k) (S100D) and the detergent dodecyl phosphocholine (DPC) were studied using NMR spectroscopy. In the absence of Ca(2+), the protein associates with DPC micelles. The micelle-associated state has intact helical secondary structures but no apparent tertiary fold. At neutral pH, Ca(2+)-loaded calbindin D(9k) does not associate with DPC micelles. However, a specific interaction is observed with individual DPC molecules at a site close to the linker between the two EF-hands. Binding to this site occurs only when Ca(2+) is bound to the protein. A reduction in pH in the absence of Ca(2+) increases the stability of the micelle-associated state. This along with the corresponding reduction in Ca(2+) affinity causes a transition to the micelle-associated state also in the presence of Ca(2+) when the pH is lowered. Site-specific analysis of the data indicates that calbindin D(9k) has a core of three tightly packed helices (A, B, and D), with a dynamic fourth helix (C) more loosely associated. Evidence is presented that the Ca(2+)-binding characteristics of the two EF-hands are distinctly different in a micelle environment. The role of calbindin D(9k) in the cell is discussed, along with the broader implications for the function of the S100 protein family.

Probing Site-Specific Calmodulin Calcium and Lanthanide Affinity by Grafting

Ca2+ binding is essential for the biological functions of calmodulin (CaM) as a trigger/sensor protein to regulate many biological processes in the Ca2+ -signaling cascade. A challenge in understanding the mechanism of Ca2+ signaling is to obtain site-specific information about the Ca2+ binding properties of individual Ca2+ -binding sites of EF-hand proteins, especially for CaM. In this paper, we report the first estimation of the intrinsic Ca2+ affinities of the four EF-hand loops of calmoduin (I-IV) by individually grafting into the domain 1 of CD2. Taking advantage of the Trp residues in the host protein, we first determined metal-binding affinities for Tb3+, Ca2+, and La3+ for all four grafted EF-loops using Tb3+ aromatic resonance energy transfer. EF-loop I exhibits the strongest binding affinity for Ca2+, La3+, and Tb3+, while EF-loop IV has the weakest metal-binding affinity. EF-loops I-IV of CaM have dissociation constants for Ca2+ of 34, 245, 185, and 814 microM, respectively, with the order I > III approximately equal to II > IV. These findings support a charge-ligand-balanced model in which both the number of negatively charged ligand residues and the balanced electrostatic dentate-dentate repulsion by the adjacent charged residues are two major determinants for the relative Ca2+ -binding affinities of EF-loops in CaM. Our grafting method provides a new strategy to obtain site-specific Ca2+ binding properties and a better estimation of the cooperativity and conformational change contributions of coupled EF-hand proteins.

Structure of the Mg2+-loaded C-lobe of cardiac troponin C bound to the N-domain of cardiac troponin I: comparison with the Ca2+-loaded structure

Cardiac troponin C (cTnC) is the Ca(2+)-binding component of the troponin complex and, as such, is the Ca(2+)-dependent switch in muscle contraction. This protein consists of two globular lobes, each containing a pair of EF-hand metal-binding sites, connected by a linker. In the N lobe, Ca(2+)-binding site I is inactive and Ca(2+)-binding site II is primarily responsible for initiation of muscle contraction. The C lobe contains Ca(2+)/Mg(2+)-binding sites III and IV, which bind Mg(2+) with lower affinity and play a structural as well as a secondary role in modulating the Ca(2+) signal. To understand the structural consequences of Ca(2+)/Mg(2+) exchange in the C lobe, we have determined the NMR solution structure of the Mg(2+)-loaded C lobe, cTnC(81-161), in a complex with the N domain of cardiac troponin I, cTnI(33-80), and compared it with a refined Ca(2+)-loaded structure. The overall tertiary structure of the Mg(2+)-loaded C lobe is very similar to that of the refined Ca(2+)-loaded structure as evidenced by the root-mean-square deviation of 0.94 A for all backbone atoms. While metal-dependent conformational changes are minimal, substitution of Mg(2+) for Ca(2+) is characterized by condensation of the C-terminal portion of the metal-binding loops with monodentate Mg(2+) ligation by the conserved Glu at position 12 and partial closure of the cTnI hydrophobic binding cleft around site IV. Thus, conformational plasticity in the Ca(2+)/Mg(2+)-dependent binding loops may represent a mechanism to modulate C-lobe cTnC interactions with the N domain of cTnI.

Design and characterization of helical peptides that inhibit the E6 protein of papillomavirus

The E6 protein from HPV type 16 binds proteins containing a seven-residue leucine-containing motif. Previous work demonstrated that peptides containing the consensus sequence are a mixture of alpha-helix and unstructured conformations. To design monomeric E6-binding peptides that are stable in aqueous solution, we used a protein grafting approach where the critical residues of the E6-binding motif of E6-associated protein, E6AP, LQELLGE, were incorporated into exposed helices of two stably folded peptide scaffolds. One series was built using the third zinc finger of the Sp1 protein, which contains a C-terminal helix. A second series was built using a Trp-cage scaffold, which contains an N-terminal helix. The chimeric peptides had very different activities in out-competing the E6-E6AP interaction. We characterized the peptides by circular dichroism spectroscopy and determined high-resolution structures by NMR methods. The E6-binding consensus motif was found to be helical in the high-quality structures, which had backbone root-mean-square deviations of less than 0.4 A. We have successfully grafted the E6-binding motif into two parent peptides to create ligands that have biological activity while preserving the stable, native fold of their scaffolds. The data also indicate that conformational change is common in E6-binding proteins during the formation of the complex with the viral E6 protein.

Small molecular ion adsorption on proteins and DNAs revealed by separation techniques

Ion binding is a term that assumes that the ion is included in the solvation sphere characterising the biomolecule. The binding forces are not clearly stated except for electrostatic attraction; weak forces (hydrogen bonds and Van der Waals forces) are likely involved. Many publications have dealt with ion binding to proteins and the consequences over the past 10 years, but only a few studies were performed using high-performance liquid chromatography (HPLC: ion exchange, reversed phase without the well-identified immobilised metal affinity chromatography) and capillary zone electrophoresis (CZE). This review focuses on the binding of proteins and DNAs mainly to the oxyanions (phosphate, borate, citrate) and amines used as buffers for both the HPLC eluent and the background electrolyte of CZE. Such specific ion adsorption on biomolecules is evidenced by physico-chemical characteristics such as the mobility or retention volume, closely associated with the net charge, which differ from the expected or experimental data obtained under the conditions of an indifferent electrolyte. It is shown that ion binding to proteins is a key parameter in the electrostatic repulsion between the free protein and a fouled membrane in the ultrafiltration separation of a protein mixture.

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.

Combining conformational flexibility and continuum electrostatics for calculating pK(a)s in proteins

Protein stability and function relies on residues being in their appropriate ionization states at physiological pH. In situ residue pK(a)s also provides a sensitive measure of the local protein environment. Multiconformation continuum electrostatics (MCCE) combines continuum electrostatics and molecular mechanics force fields in Monte Carlo sampling to simultaneously calculate side chain ionization and conformation. The response of protein to charges is incorporated both in the protein dielectric constant (epsilon(prot)) of four and by explicit conformational changes. The pK(a) of 166 residues in 12 proteins was determined. The root mean square error is 0.83 pH units, and >90% have errors of <1 pH units whereas only 3% have errors >2 pH units. Similar results are found with crystal and solution structures, showing that the method's explicit conformational sampling reduces sensitivity to the initial structure. The outcome also changes little with protein dielectric constant (epsilon(prot) 4-20). Multiconformation continuum electrostatics titrations show coupling of conformational flexibility and changes in ionization state. Examples are provided where ionizable side chain position (protein G), Asn orientation (lysozyme), His tautomer distribution (RNase A), and phosphate ion binding (RNase A and H) change with pH. Disallowing these motions changes the calculated pK(a).

Structure, specificity, and mode of interaction for bacterial albumin-binding modules

We have determined the solution structure of an albumin binding domain of protein G, a surface protein of group C and G streptococci. We find that it folds into a left handed three-helix bundle similar to the albumin binding domain of protein PAB from Peptostreptococcus magnus. The two domains share 59% sequence identity, are thermally very stable, and bind to the same site on human serum albumin. The albumin binding site, the first determined for this structural motif known as the GA module, comprises residues spanning the first loop to the beginning of the third helix and includes the most conserved region of GA modules. The two GA modules have different affinities for albumin from different species, and their albumin binding patterns correspond directly to the host specificity of C/G streptococci and P. magnus, respectively. These studies of the evolution, structure, and binding properties of the GA module emphasize the power of bacterial adaptation and underline ecological and medical problems connected with the use of antibiotics.

Continuum secondary structure captures protein flexibility

The DSSP program assigns protein secondary structure to one of eight states. This discrete assignment cannot describe the continuum of thermal fluctuations. Hence, a continuous assignment is proposed. Technically, the continuum results from averaging over ten discrete DSSP assignments with different hydrogen bond thresholds. The final continuous assignment for a single NMR model successfully reflected the structural variations observed between all NMR models in the ensemble. The structural variations between NMR models were verified to correlate with thermal motion; these variations were captured by the continuous assignments. Because the continuous assignment reproduces the structural variation between many NMR models from one single model, functionally important variation can be extracted from a single X-ray structure. Thus, continuous assignments of secondary structure may affect future protein structure analysis, comparison, and prediction.

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.

Ion-induced conformational and stability changes in Nereis sarcoplasmic calcium binding protein: evidence that the APO state is a molten globule

Nereis sarcoplasmic Ca(2+)-binding protein (NSCP) is a calcium buffer protein that binds Ca(2+) ions with high affinity but is also able to bind Mg(2+) ions with high positive cooperativity. We investigated the conformational and stability changes induced by the two metal ions. The thermal reversible unfolding, monitored by circular dichroism spectroscopy, shows that the thermal stability is maximum at neutral pH and increases in the order apo < Mg(2+) < Ca(2+). The stability against chemical denaturation (urea, guanidinium chloride) studied by circular dichroism or intrinsic fluorescence was found to have a similar ion dependence. To explore in more detail the structural basis of stability, we used the fluorescent probes to evaluate the hydrophobic surface exposure in the different ligation states. The apo-NSCP exhibits accessible hydrophobic surfaces, able to bind fluorescent probes, in clear contrast with denatured or Ca(2+)/Mg(2+)-bound states. Gel filtration experiments showed that, although the metal-bound NSCP has a hydrodynamic volume in agreement with the molecular mass, the volume of the apo form is considerably larger. The present results demonstrate that the apo state has many properties in common with the molten globule. The possible factors of the metal-dependent structural changes and stability are discussed.

Development of computational and graphical tools for analysis of movement and flexibility in large molecules

We developed a computer program for the calculation and display of the difference distance matrices (DDMs) of macromolecules that has the ability to compare multiple structures simultaneously. To demonstrate its use, a data set of atoms for superimposition of the HIV-1 reverse transcriptase enzyme was defined using the coordinates for the 21 available crystal structures of this enzyme and its complexes. The DDM technique for superimposition data set generation allows selection of atoms that are invariant in all structures, is free from user bias, and represents the most accurate and precise method of producing such subsets. Comparison of this technique was made against other published methods of generating superimposition data sets, and it was found that significant differences in magnitude and trends of atom movements are observed depending on which data set was used.

Peptide and metal ion-dependent association of isolated helix-loop-helix calcium binding domains: studies of thrombic fragments of calmodulin

Calmodulin (CaM), the ubiquitous, eukaryotic, bilobal calcium-binding regulatory protein, has been cleaved by thrombin to create two fragments. TM1 (1-106) and TM2 (107-148). NMR and CD results indicate that TMI and TM2 can associate in the presence of Ca2+ to form a complex similar to native CaM, even though the cleavage site is not in the linker region between two helix-loop-helix domains, but rather within an alpha-helix. Cadmium-113 NMR results show that this complex has enhanced metal-ion binding properties when compared to either TM1 or TM2 alone. This complex can bind several CaM-binding target peptides, as shown by gel bandshift assays, circular dichroism spectra, and 13C NMR spectra of biosynthetically methyl-13C-Met-labeled TM1 and TM2; moreover, gel bandshift assays show that the addition of a target peptide strengthens the interactions between TM1 and TM2 and increases the stability of the complex. Cadmium-113 NMR spectra indicate that the TM1:TM2 complex can also bind the antipsychotic drug trifluoperazine. However, in contrast to CaM:peptide complexes, the TM1:TM2:peptide complexes are disrupted by 4 M urea; moreover, TM1 and TM2 in combination are unable to activate CaM-dependent enzymes. This suggests that TM1:TM2 mixtures cannot bind target molecules as tightly as intact CaM, or perhaps that binding occurs but additional interactions with the target enzymes that are necessary for proper activation are perturbed by the proteolytic cleavage. The results presented here reflect the importance of the existence of helix-loop-helix Ca2+-binding domains in pairs in proteins such as CaM, and extend the understanding of the association of such domains in this class of proteins in general.

Solution structure and dynamics of the DNA-binding domain of the adipocyte-transcription factor FREAC-11

Transcription factors of the forkhead type share a highly conserved DNA-binding domain of about 100 amino acid residues. FREAC-11, expressed in adipocytes, belongs to this class. Here, we report on NMR studies that established the three-dimensional structure of the FREAC-11, DNA-binding domain. Although apparent similarities to the structures of other members within the forkhead family are observed, the structure also reveals some remarkable differences. Along with the complementary dynamics, the data provide insight into the fundamentals of sequence specificity within a highly conserved motif.

Structural dynamics in the C-terminal domain of calmodulin at low calcium levels

Calmodulin undergoes Ca2+-induced structural rearrangements that are intimately coupled to the regulation of numerous cellular processes. The C-terminal domain of calmodulin has previously been observed to exhibit conformational exchange in the absence of Ca2+. Here, we characterize further the conformational dynamics in the presence of low concentrations of Ca2+ using 15N spin relaxation experiments. The analysis included 1H-15N dipolar/15N chemical shift anisotropy interference cross-correlation relaxation rates to improve the description of the exchange processes, as well as the picosecond to nanosecond dynamics. Conformational transitions on microsecond to millisecond time scales were revealed by exchange contributions to the transverse auto-relaxation rates. In order to separate the effects of Ca2+ exchange from intramolecular conformational exchange processes in the apo state, transverse auto-relaxation rates were measured at different concentrations of free Ca2+. The results reveal a Ca2+-dependent contribution due mainly to exchange between the apo and (Ca2+)1 states with an apparent Ca2+ off-rate of approximately 5115 s(-1), as well as Ca2+-independent contributions due to conformational exchange within the apo state. 15N chemical shift differences estimated from the exchange data suggest that the first Ca2+ binds preferentially to loop IV. Thus, characterization of chemical exchange as a function of Ca2+ concentration has enabled the extraction of unique information on the rapidly exchanging and weakly populated (<10 %) (Ca2+)1 state that is otherwise inaccessible to direct study due to strongly cooperative Ca2+ binding. The conformational exchange within the apo state appears to involve transitions between a predominantly populated closed conformation and a smaller population of more open conformations. The picosecond to nanosecond dynamics of the apo state are typical of a well-folded protein, with reduced amplitudes of motions in the helical segments, but with significant flexibility in the Ca2+-binding loops. Comparisons with order parameters for skeletal troponin C and calbindin D9k reveal key structural and dynamical differences that correlate with the different Ca2+-binding properties of these proteins.

Determination of the metal-binding cooperativity of wild-type and mutant calbindin D9K by electrospray ionization mass spectrometry

Since the initial reports showing the ability of electrospray ionization mass spectrometry (ESI-MS) to study intact noncovalent biomolecular complexes, an increasing number of uses for this technique in studying biochemical systems is emerging. We have investigated the ability of ESI-MS to characterize the metal-binding properties of calcium (Ca2+) binding proteins by studying the incorporation of Ca2+ and cadmium (Cd2+) into wild-type and mutant calbindin D9K. ESI-MS showed that wild-type calbindin D9K binds two Ca2+ ions with similar affinities while the binding of two Cd2+ ions is sequential, as is the binding of the two Ca2+ or Cd2+ ions to the N56A mutant of calbindin. The binding of Ca2+ to the wild-type protein was clearly seen to be cooperative. These results demonstrate the potential efficacy of ESI-MS to discriminate between cooperative and independent site metal binding to metalloproteins.

Cooperative cyclic interactions involved in metal binding to pairs of sites in EF-hand proteins

This study focuses on a closed net of electron-pair donor-acceptor interactions, present in the core of all metal-bound EF-hand pairs, that link both metal ions across a short two-stranded beta-sheet. A molecular model based on the above cycle of interactions was studied using semi-empirical molecular orbital quantum mechanical methods. The calculations indicate that the interactions in the model cycle are cooperative, that is, that the interaction energy of the cyclic structure is greater than that of the sum of isolated interactions between its components. The cooperativity in this cycle can be attributed to an increase in the stability of the interactions resulting from a mutual polarisation of the associated groups. The predicted polarisation of the amide groups in the cycle is in agreement with experimental NMR 15N deshielding observed for these amide groups upon metal binding. Experimental observations of strengthening of the beta-sheet hydrogen bonds are also consistent with the model calculations. By this mechanism, the binding of the first metal ion would enhance the binding of the second metal ion, and thus, the intradomain cooperativity in cation binding of calmodulin and related EF-hand proteins can be ascribed, at least partly, to this short-range molecular mechanism.

Ca2+ binding and conformational changes in a calmodulin domain

Calcium activation of the C-terminal domain of calmodulin was studied using 1H and 15N NMR spectroscopy. The important role played by the conserved bidentate glutamate Ca2+ ligand in the binding loops is emphasized by the striking effects resulting from a mutation of this glutamic acid to a glutamine, i.e. E104Q in loop III and E140Q in loop IV. The study involves determination of Ca2+ binding constants, assignments, and structural characterizations of the apo, (Ca2+)1, and (Ca2+)2 states of the E104Q mutant and comparisons to the wild-type protein and the E140Q mutant [Evenäs et al. (1997) Biochemistry 36, 3448-3457]. NMR titration data show sequential Ca2+ binding in the E104Q mutant. The first Ca2+ binds to loop IV and the second to loop III, which is the order reverse to that observed for the E140Q mutant. In both mutants, the major structural changes occur upon Ca2+ binding to loop IV, which implies a different response to Ca2+ binding in the N- and C-terminal EF-hands. Spectral characteristics show that the (Ca2+)1 and (Ca2+)2 states of the E104Q mutant undergo global exchange on a 10-100 micros time scale between conformations seemingly similar to the closed and open structures of this domain in wild-type calmodulin, paralleling earlier observations for the (Ca2+)2 state of the E140Q mutant, indicating that both glutamic acid residues, E104 and E140, are required for stabilization of the open conformation in the (Ca2+)2 state. To verify that the NOE constraints cannot be fulfilled in a single structure, solution structures of the (Ca2+)2 state of the E104Q mutant are calculated. Within the ensemble of structures the precision is good. However, the clearly dynamic nature of the state, a large number of violated distance restraints, ill-defined secondary structural elements, and comparisons to the structures of calmodulin indicate that the ensemble does not provide a good picture of the (Ca2+)2 state of the E104Q mutant but rather represents the distance-averaged structure of at least two distinct different conformations.

Electrostatic coupling to pH-titrating sites as a source of cooperativity in protein-ligand binding

This paper describes an alternative mechanism for the cooperative binding of charged ligands to proteins. The ligand-binding sites are electrostatically coupled to protein side chains that can undergo protonation and deprotonation. The binding of one ligand alters the protein's protonation equilibrium in a manner that makes the the binding of the second ligand more favorable. This mechanism requires no conformational change to produce a cooperative effect, although it is not exclusive of conformational change. We present a theoretical description of the mechanism, and calculations on three kinds of systems: A model system containing one protonation site and two ligand-binding sites; a model system containing two protonation sites and two ligand-binding sites; and calbindin D9k, which contains two Ca2+-binding sites and 30 protonation sites. For the one-protonation-site model, it is shown that the influence of the protonation site can only be cooperative. The competition of this effect with the anticooperative effect of ligand-ligand repulsion is studied in detail. For the two-protonation site model, the effect can be either cooperative or, in special cases, anticooperative. For calbindin D9k, the calculations predict that six protonation sites in or near the ligand-binding sites make a cooperative contribution that approximately cancels the anticooperative effect of Ca2+-Ca2+ repulsion, accounting for more than half of the total cooperative effect that is needed to overcome repulsion and produce the net cooperativity observed experimentally. We argue that cooperative mechanisms of the kind described here are likely when there is more than one ligand-binding site in a protein domain.

Solution structure of the N-terminal EGF-like domain from human factor VII

Blood coagulation is initiated by Ca(2+)-dependent binding of coagulation factor VIIa (FVIIa) to its cofactor, tissue factor (TF). The TF:FVIIa complex activates factors IX and X, ultimately leading to the formation of thrombin and the coagulation of blood. FVII consists of an N-terminal gamma-carboxyglutamic-acid-containing (Gla) domain followed by two epidermal growth factor (EGF) like domains, the first of which can bind one Ca2+ ion (Kd approximately 150 microM) and a C-terminal serine protease domain. Using 1H nuclear magnetic resonance spectroscopy, we have determined the solution structure of a synthetic N-terminal EGF-like domain (EGF1) of human FVII (residues 45-85) in the absence of Ca2+. A comparison of this structure of apo EGF1 with the Ca(2+)-bound EGF1 in the complex of FVIIa and TF [Banner, D. W., et al. (1996) Nature 380, 41-46] suggests that the structural changes in the EGF1 domain upon Ca2+ binding are minor and are concentrated near the Ca(2+)-binding site, which is facing away from the TF interaction surface. Amino acid side chains that are crucial for the binding of FVII to TF show a similar conformation in both structures and are therefore unlikely to directly influence the Ca(2+)-dependent binding of FVII to TF. As Ca2+ binding to EGF1 does not lead to a conformational change in the residues constituting the interaction surface for binding to TF, our results are consistent with the idea that the altered orientation between the Gla and EGF1 domains that result from Ca2+ binding is responsible for the increased affinity of FVII/FVIIa for TF.

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

Hydrophobic core residues have a marked influence on the Ca2+-binding properties of calbindin D9k, even though there are no direct contacts between these residues and the bound Ca2+ ions. Eleven different mutants with substitutions in the hydrophobic core were produced, and their equilibrium Ca2+-binding constants measured from Ca2+ titrations in the presence of chromophoric chelators. The Ca2+-dissociation rate constants were estimated from Ca2+ titrations followed by 1H NMR1 and were measured more accurately using stopped-flow fluorescence. The parameters were measured at four KCl concentrations to assess the salt dependence of the perturbations. The high similarity between the NMR spectra of mutants and wild-type calbindin D9k suggests that the structure is largely unperturbed by the substitutions. More detailed NMR investigations of the mutant in which Val61 is substituted by Ala showed that the mutation causes only very minimal perturbations in the immediate vicinity of residue 61. Substitutions of alanines or glycines for bulky residues in the center of the core were found to have significant effects on both Ca2+ affinity and dissociation rates. These substitutions caused a reduction in affinity and an increase in off-rate. Small effects, both increases and decreases, were observed for substitutions involving residues far from the Ca2+ sites and toward the outer part of the hydrophobic core. The mutant with the substitution Phe66 --> Trp behaved differently from all other mutants, and displayed a 25-fold increase in overall affinity of binding two Ca2+ ions and a 6-fold reduction in calcium dissociation rate. A strong correlation (R = 0.94) was found between the observed Ca2+-dissociation rates and affinities, as well as between the salt dependence of the off-rate and the distance to the nearest Ca2+-coordinating atom. There was also a strong correlation (R = 0.95) between the Ca2+ affinity and stability of the Ca2+ state and a correlation (R = 0. 69) between the Ca2+ affinity and stability of the apo state, as calculated from the results in the present and preceding paper in this issue [Julenius, K., Thulin, E., Linse, S., and Finn, B. E. (1998) Biochemistry 37, 8915-8925]. The change in salt dependencies of koff and cooperativity were most pronounced for residues completely buried in the core of the protein (solvent accessible surface area approximately 0). Altogether, the results suggest that the hydrophobic core residues promote Ca2+ binding both by contributing to the preformation of the Ca2+ sites in the apo state and by preferentially stabilizing the Ca2+-bound state.

Ca2+ coordination to backbone carbonyl oxygen atoms in calmodulin and other EF-hand proteins: 15N chemical shifts as probes for monitoring individual-site Ca2+ coordination

Examination of the NMR 15N chemical shifts of a number of EF-hand proteins shows that the shift value for the amido nitrogen of the residue in position 8 of a canonical EF-hand loop (or position 10 of a pseudo EF-hand loop) provides a good indication of metal occupation of that site. The NH of the residue in position 8 is covalently bonded to the carbonyl of residue 7, the only backbone carbonyl that coordinates to the metal ion in a canonical EF-hand loop. Upon metal coordination to this carbonyl, there is an appreciable deshielding of the 15N nucleus at position 8 (+4 to +8 ppm) due to the polarization of the O(7)=C(7)-N(8) amido group and the corresponding reduction in the electron density of the nitrogen atom. This deshielding effect is effectively independent of the binding of metal to the other site of an EF-hand pair, allowing the 15N shifts to be used as probes for site-specific occupancy of metal binding sites. In addition, a Ca2+-induced change in side-chain Halpha-Calpha-Cbeta-Hbeta torsion angle for isoleucine or valine residues in position 8 can also contribute to the deshielding of the amide 15N nucleus. This conformational effect occurs only in sites I or III and takes place upon binding a Ca2+ ion to the other site of an EF-hand pair (site II or IV) regardless of whether the first site is occupied. The magnitude of this effect is in the range +5 to +7 ppm. A Ca2+ titration of 15N-labeled apo-calmodulin was performed using 2D 1H-15N HSQC NMR spectra. The changes in the 15N chemical shifts and intensities for the peaks corresponding to the NH groups of residues in position 8 of the EF-hand loops allowed the amount of metal bound at sites II, III and IV to be monitored directly at partial degrees of saturation. The peak corresponding to site I could only be monitored at the beginning and end of the titration because of line broadening effects in the intermediate region of the titration. Sites III and IV both titrate preferentially and the results demonstrate clearly that sites in either domain fill effectively in parallel, consistent with a significant positive intradomain cooperativity of calcium binding.

Solution structure and dynamics of a designed monomeric variant of the lambda Cro repressor

The solution structure of a monomeric variant of the lambda Cro repressor has been determined by multidimensional NMR. Cro K56[DGEVK] differs from wild-type Cro by the insertion of five amino acids at the center of the dimer interface. 1H and 15N resonances for 70 of the 71 residues have been assigned. Thirty-two structures were calculated by hybrid distance geometry/simulated annealing methods using 463 NOE-distance restraints, 26 hydrogen-bond, and 39 dihedral-angle restraints. The root-mean-square deviation (RMSD) from the average structure for atoms in residues 3-60 is 1.03 +/- 0.44 A for the peptide backbone and 1.6 +/- 0.73 A for all nonhydrogen atoms. The overall structure conforms very well to the original design. Although the five inserted residues form a beta hairpin as expected, this engineered turn as well as other turns in the structure are not well defined by the NMR data. Dynamics studies of backbone amides reveal T1/T2 ratios of residues in the alpha2-alpha3, beta2-beta3, and engineered turn that are reflective of chemical exchange or internal motion. The solution structure and dynamics are discussed in light of the conformational variation that has been observed in other Cro structures, and the importance of flexibility in DNA recognition.

Sequence and context dependence of EF-hand loop dynamics. An 15N relaxation study of a calcium-binding site mutant of calbindin D9k

The influence of amino acid sequence and structural context on the backbone dynamics of EF-hand calcium-binding loops was investigated using 15N spin relaxation measurements on the calcium-free state of the calbindin D9k mutant (A14D+A15Delta+P20Delta+N21G+P43M), in which the N-terminal pseudo-EF-hand loop, characteristic of S100 proteins, was engineered so as to conform with the C-terminal consensus EF-hand loop. The results were compared to a previous study of the apo state of the wild-type-like P43G calbindin D9k mutant. In the helical regions, the agreement with the P43G data is excellent, indicating that the structure and dynamics of the protein core are unaffected by the substitutions in the N-terminal loop. In the calcium-binding loops, the flexibility is drastically decreased compared to P43G, with the modified N-terminal loop showing a motional restriction comparable to that of the surrounding helixes. As in P43G, the motions in the C-terminal loop are less restricted than in the N-terminal loop. Differences in key hydrogen-bonding interactions correlate well with differences in dynamics and offer insights into the relationship between structure and dynamics of these EF-hand loops. It appears that the entire N-terminal EF-hand is built to form a rigid structure that allows calcium binding with only minor rearrangements and that the structural and dynamical properties of the entire EF-hand--rather than the loop sequence per se--is the major determinant of loop flexibility in this system.

The three-dimensional structure of Ca(2+)-bound calcyclin: implications for Ca(2+)-signal transduction by S100 proteins

BACKGROUND

Calcyclin is a member of the S100 subfamily of EF-hand Ca(2+)-binding proteins. This protein has implied roles in the regulation of cell growth and division, exhibits deregulated expression in association with cell transformation, and is found in high abundance in certain breast cancer cell lines. The novel homodimeric structural motif first identified for apo calcyclin raised the possibility that S100 proteins recognize their targets in a manner that is distinctly different from that of the prototypical EF-hand Ca2+ sensor, calmodulin. The NMR solution structure of Ca(2+)-bound calcyclin has been determined in order to identify Ca(2+)-induced structural changes and to obtain insights into the mechanism of Ca(2+)-triggered target protein recognition.

RESULTS

The three-dimensional structure of Ca(2+)-bound calcyclin was calculated with 1372 experimental constraints, and is represented by an ensemble of 20 structures that have a backbone root mean square deviation of 1.9 A for the eight helices. Ca(2+)-bound calcyclin has the same symmetric homodimeric fold as observed for the apo protein. The helical packing within the globular domains and the subunit interface also change little upon Ca2+ binding. A distinct homology was found between the Ca(2+)-bound states of the calcyclin subunit and the monomeric S100 protein calbindin D9k.

CONCLUSIONS

Only very modest Ca(2+)-induced changes are observed in the structure of calcyclin, in sharp contrast to the domain-opening that occurs in calmodulin and related Ca(2+)-sensor proteins. Thus, calcyclin, and by inference other members of the S100 family, must have a different mode for transducing Ca2+ signals and recognizing target proteins. This proposal raises significant questions concerning the purported roles of S100 proteins as Ca2+ sensors.

An interaction-based analysis of calcium-induced conformational changes in Ca2+ sensor proteins

Calcium sensor proteins translate transient increases in intracellular calcium levels into metabolic or mechanical responses, by undergoing dramatic conformational changes upon Ca2+ binding. A detailed analysis of the calcium binding-induced conformational changes in the representative calcium sensors calmodulin (CaM) and troponin C was performed to obtain insights into the underlying molecular basis for their response to the binding of calcium. Distance difference matrices, analysis of interresidue contacts, comparisons of interhelical angles, and inspection of structures using molecular graphics were used to make unbiased comparisons of the various structures. The calcium-induced conformational changes in these proteins are dominated by reorganization of the packing of the four helices within each domain. Comparison of the closed and open conformations confirms that calcium binding causes opening within each of the EF-hands. A secondary analysis of the conformation of the C-terminal domain of CaM (CaM-C) clearly shows that CaM-C occupies a closed conformation in the absence of calcium that is distinct from the semi-open conformation observed in the C-terminal EF-hand domains of myosin light chains. These studies provide insight into the structural basis for these changes and into the differential response to calcium binding of various members of the EF-hand calcium-binding protein family. Factors contributing to the stability of the Ca2+-loaded open conformation are discussed, including a new hypothesis that critical hydrophobic interactions stabilize the open conformation in Ca2+ sensors, but are absent in "non-sensor" proteins that remain closed upon Ca2+ binding. A role for methionine residues in stabilizing the open conformation is also proposed.

High-resolution polypeptide structure in a lamellar phase lipid environment from solid state NMR derived orientational constraints

BACKGROUND

Solid-state nuclear magnetic resonance (NMR) spectroscopy provides novel structural constraints from uniformly aligned samples. These orientational constraints orient specific atomic sites with respect to the magnetic field direction and the unique molecular axis of alignment. Solid-state NMR is uniquely and ideally suited for providing such structural constraints on polypeptides and proteins in a lamellar phase lipid environment. Membrane protein structure represents a great challenge for structural biologists; a new approach for characterizing high resolution three-dimensional structure in such an environment is needed.

RESULTS

The optimal use of orientational constraints for defining three-dimensional structures is demonstrated with the elucidation of the gramicidin A channel structure at high resolution. Initial structures are refined against both the experimental constraints and the CHARMM energy using a novel simulated-annealing protocol to define torsion angle solutions with an error bar of approximately +/- 5 degrees.

CONCLUSIONS

This analysis results in the determination of a high-resolution, time averaged structure of gramicidin A obtained in a lipid bilayer environment above the gel-to-liquid crystalline phase transition temperature. It is demonstrated that solid-state NMR can be used to establish polypeptide, and potentially protein, structures in such an environment. Furthermore, this high-resolution structure is demonstrated to provide new insights into polypeptide function. For the gramicidin A channel the roles of the indole groups that facilitate ion transport and details of the cation solvation environment provided by the amide oxygens are characterized.