Calcium-induced structural changes and domain autonomy in calmodulin

Noncovalent Peptide Assembly Enables Crystalline, Permutable, and Reactive Thiol Frameworks

Though thiols are exceptionally versatile, their high reactivity has also hindered the synthesis and characterization of well-defined thiol-containing porous materials. Leveraging the mild conditions of the noncovalent peptide assembly, we readily synthesized and characterized a number of frameworks with thiols displayed at many unique positions and in several permutations. Importantly, nearly all assemblies were structurally determined using single-crystal X-ray diffraction to reveal their rich sequence-structure landscape and the cooperative noncovalent interactions underlying their assembly. These observations and supporting molecular dynamics calculations enabled rational engineering by the positive and negative design of noncovalent interactions. Furthermore, the thiol-containing frameworks undergo diverse single-crystal-to-single-crystal reactions, including toxic metal ion coordination (e.g., Cd2+, Pb2+, and Hg2+), selective uptake of Hg2+ ions, and redox transformations. Notably, we find a framework that supports thiol-nitrosothiol interconversion, which is applicable for biocompatible nitric oxide delivery. The modularity, ease of synthesis, functionality, and well-defined nature of these peptide-based thiol frameworks are expected to accelerate the design of complex materials with reactive active sites.

Half-calcified calmodulin promotes basal activity and inactivation of the L-type calcium channel CaV1.2

The L-type Ca2+ channel CaV1.2 controls gene expression, cardiac contraction, and neuronal activity. Calmodulin (CaM) governs CaV1.2 open probability (Po) and Ca2+-dependent inactivation (CDI) but the mechanisms remain unclear. Here, we present electrophysiological data that identify a half Ca2+-saturated CaM species (Ca2/CaM) with Ca2+ bound solely at the third and fourth EF-hands (EF3 and EF4) under resting Ca2+ concentrations (50-100 nM) that constitutively preassociates with CaV1.2 to promote Po and CDI. We also present an NMR structure of a complex between the CaV1.2 IQ motif (residues 1644-1665) and Ca2/CaM12', a calmodulin mutant in which Ca2+ binding to EF1 and EF2 is completely disabled. We found that the CaM12' N-lobe does not interact with the IQ motif. The CaM12' C-lobe bound two Ca2+ ions and formed close contacts with IQ residues I1654 and Y1657. I1654A and Y1657D mutations impaired CaM binding, CDI, and Po, as did disabling Ca2+ binding to EF3 and EF4 in the CaM34 mutant when compared to WT CaM. Accordingly, a previously unappreciated Ca2/CaM species promotes CaV1.2 Po and CDI, identifying Ca2/CaM as an important mediator of Ca signaling.

Interdomain Dynamics via Paramagnetic NMR on the Highly Flexible Complex Calmodulin/Munc13-1

Paramagnetic NMR constraints are very useful to study protein interdomain motion, but their interpretation is not always straightforward. On the example of the particularly flexible complex Calmodulin/Munc13-1, we present a new approach to characterize this motion with pseudocontact shifts and residual dipolar couplings. Using molecular mechanics, we sampled the conformational space of the complex and used a genetic algorithm to find ensembles that are in agreement with the data. We used the Bayesian information criterion to determine the ideal ensemble size. This way, we were able to make an accurate, unambiguous, reproducible model of the interdomain motion of Calmodulin/Munc13-1 without prior knowledge about the domain orientation from crystallography.

High-affinity tamoxifen analogues retain extensive positional disorder when bound to calmodulin

Using a combination of NMR and fluorescence measurements, we have investigated the structure and dynamics of the complexes formed between calcium-loaded calmodulin (CaM) and the potent breast cancer inhibitor idoxifene, a derivative of tamoxifen. High-affinity binding (K d ∼ 300  nM) saturates with a 2 : 1 idoxifene : CaM complex. The complex is an ensemble where each idoxifene molecule is predominantly in the vicinity of one of the two hydrophobic patches of CaM but, in contrast with the lower-affinity antagonists TFP, J-8, and W-7, does not substantially occupy the hydrophobic pocket. At least four idoxifene orientations per domain of CaM are necessary to satisfy the intermolecular nuclear Overhauser effect (NOE) restraints, and this requires that the idoxifene molecules switch rapidly between positions. The CaM molecule is predominantly in the form where the N and C-terminal domains are in close proximity, allowing for the idoxifene molecules to contact both domains simultaneously. Hence, the 2 : 1 idoxifene : CaM complex illustrates how high-affinity binding occurs without the loss of extensive positional dynamics.

Tribbles Pseudokinase 2 (TRIB2) Regulates Expression of Binding Partners in Bovine Granulosa Cells

Members of the Tribbles (TRIB) family of pseudokinases are critical components of intracellular signal transduction pathways in physiological and pathological processes. TRIBs, including TRIB2, have been previously shown as signaling mediators and scaffolding proteins regulating numerous cellular events such as proliferation, differentiation and cell death through protein stability and activity. However, the signaling network associated with TRIB2 and its binding partners in granulosa cells during ovarian follicular development is not fully defined. We previously reported that TRIB2 is differentially expressed in growing dominant follicles while downregulated in ovulatory follicles following the luteinizing hormone (LH) surge or human chorionic gonadotropin (hCG) injection. In the present study, we used the yeast two-hybrid screening system and in vitro coimmunoprecipitation assays to identify and confirm TRIB2 interactions in granulosa cells (GCs) of dominant ovarian follicles (DFs), which yielded individual candidate binding partners including calmodulin 1 (CALM1), inhibin subunit beta A (INHBA), inositol polyphosphate phosphatase-like 1 (INPPL1), 5'-nucleotidase ecto (NT5E), stearoyl-CoA desaturase (SCD), succinate dehydrogenase complex iron sulfur subunit B (SDHB) and Ras-associated protein 14 (RAB14). Further analyses showed that all TRIB2 binding partners are expressed in GCs of dominant follicles but are differentially regulated throughout the different stages of follicular development. CRISPR/Cas9-driven inhibition along with pQE-driven overexpression of TRIB2 showed that TRIB2 differently regulates expression of binding partners, which reveals the importance of TRIB2 in the control of gene expression linked to various biological processes such as proliferation, differentiation, cell migration, apoptosis, calcium signaling and metabolism. These data provide a larger view of potential TRIB2-regulated signal transduction pathways in GCs and provide strong evidence that TRIB2 may act as a regulator of target genes during ovarian follicular development.

Interactions of a Signal Transduction Protein Investigated by Fluorescence Stopped-Flow Kinetics

To understand cellular processes such as biochemical pathways and signaling networks, we need to understand binding and reaction rates of often competing reactions, their dependence on cellular concentrations of participating molecules, and the regulation of these rates through allostery, posttranslational modifications, or other mechanisms. To do so, we break these systems down into their elementary steps, which are almost invariably either unimolecular or bimolecular reactions that frequently occur on sub-second, often sub-millisecond, time scales. Rapid mixing techniques, which generally achieve mixing in less than 2 ms, are generally suitable for the study of such reactions. The application of these techniques to the study of enzyme mechanisms is described in several excellent texts (Cornish-Bowden, Fundamentals of enzyme kinetics, 1995; Gutfreund, Kinetics for the life sciences. Receptors, transmitters and catalysis, 1995); flow techniques are used to study individual steps by monitoring the approach to equilibrium (the pre-steady state) under single turnover conditions.The individual steps in complex biochemical reaction schemes determine how fast systems can respond to incoming signals and adapt to changed conditions [1, 2]. This chapter is concerned with in vitro techniques that have been developed to study fast reactions in solution, and we present the study of various interactions of calmodulin as an example. The kinetic information obtained with these techniques is indispensable for understanding the dynamics of biochemical processes and complements the static structural and thermodynamic information available from X-ray crystallography, NMR, and equilibrium binding studies.

Analyzing In Silico the Relationship Between the Activation of the Edema Factor and Its Interaction With Calmodulin

Molecular dynamics (MD) simulations have been recorded on the complex between the edema factor (EF) of Bacilllus anthracis and calmodulin (CaM), starting from a structure with the orthosteric inhibitor adefovir bound in the EF catalytic site. The starting structure has been destabilized by alternately suppressing different co-factors, such as adefovir ligand or ions, revealing several long-distance correlations between the conformation of CaM, the geometry of the CaM/EF interface, the enzymatic site and the overall organization of the complex. An allosteric communication between CaM/EF interface and the EF catalytic site, highlighted by these correlations, was confirmed by several bioinformatics approaches from the literature. A network of hydrogen bonds and stacking interactions extending from the helix V of of CaM, and the residues of the switches A, B and C, and connecting to catalytic site residues, is a plausible candidate for the mediation of allosteric communication. The greatest variability in volume between the different MD conditions was also found for cavities present at the EF/CaM interface and in the EF catalytic site. The similarity between the predictions from literature and the volume variability might introduce the volume variability as new descriptor of allostery.

Conformational ensembles of non-peptide ω-conotoxin mimetics and Ca+2 ion binding to human voltage-gated N-type calcium channel Cav2.2

Chronic neuropathic pain is the most complex and challenging clinical problem of a population that sets a major physical and economic burden at the global level. Ca²⁺-permeable channels functionally orchestrate the processing of pain signals. Among them, N-type voltage-gated calcium channels (VGCC) hold prominent contribution in the pain signal transduction and serve as prime targets for synaptic transmission block and attenuation of neuropathic pain. Here, we present detailed in silico analysis to comprehend the underlying conformational changes upon Ca²⁺ ion passage through Ca_v2.2 to differentially correlate subtle transitions induced via binding of a conopeptide-mimetic alkylphenyl ether-based analogue MVIIA. Interestingly, pronounced conformational changes were witnessed at the proximal carboxyl-terminus of Ca_v2.2 that attained an upright orientation upon Ca⁺² ion permeability. Moreover, remarkable changes were observed in the architecture of channel tunnel. These findings illustrate that inhibitor binding to Ca_v2.2 may induce more narrowing in the pore size as compared to Ca²⁺ binding through modulating the hydrophilicity pattern at the selectivity region. A significant reduction in the tunnel volume at the selectivity filter and its enhancement at the activation gate of Ca⁺²-bound Ca_v2.2 suggests that ion binding modulates the outward splaying of pore-lining S6 helices to open the voltage gate. Overall, current study delineates dynamic conformational ensembles in terms of Ca⁺² ion and MVIIA-associated structural implications in the Ca_v2.2 that may help in better therapeutic intervention to chronic and neuropathic pain management.

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.

Neurogranin stimulates Ca2+/calmodulin-dependent kinase II by suppressing calcineurin activity at specific calcium spike frequencies

Calmodulin sits at the center of molecular mechanisms underlying learning and memory. Its complex and sometimes opposite influences, mediated via the binding to various proteins, are yet to be fully understood. Calcium/calmodulin-dependent protein kinase II (CaMKII) and calcineurin (CaN) both bind open calmodulin, favoring Long-Term Potentiation (LTP) or Depression (LTD) respectively. Neurogranin binds to the closed conformation of calmodulin and its impact on synaptic plasticity is less clear. We set up a mechanistic computational model based on allosteric principles to simulate calmodulin state transitions and its interactions with calcium ions and the three binding partners mentioned above. We simulated calcium spikes at various frequencies and show that neurogranin regulates synaptic plasticity along three modalities. At low spike frequencies, neurogranin inhibits the onset of LTD by limiting CaN activation. At intermediate frequencies, neurogranin facilitates LTD, but limits LTP by precluding binding of CaMKII with calmodulin. Finally, at high spike frequencies, neurogranin promotes LTP by enhancing CaMKII autophosphorylation. While neurogranin might act as a calmodulin buffer, it does not significantly preclude the calmodulin opening by calcium. On the contrary, neurogranin synchronizes the opening of calmodulin's two lobes and promotes their activation at specific frequencies. Neurogranin suppresses basal CaN activity, thus increasing the chance of CaMKII trans-autophosphorylation at high-frequency calcium spikes. Taken together, our study reveals dynamic regulatory roles played by neurogranin on synaptic plasticity, which provide mechanistic explanations for opposing experimental findings.

Understanding Calcium-Dependent Conformational Changes in S100A1 Protein: A Combination of Molecular Dynamics and Gene Expression Study in Skeletal Muscle

The S100A1 protein, involved in various physiological activities through the binding of calcium ions (Ca²⁺), participates in several protein-protein interaction (PPI) events after Ca²⁺-dependent activation. The present work investigates Ca²⁺-dependent conformational changes in the helix-EF hand-helix using the molecular dynamics (MD) simulation approach that facilitates the understanding of Ca²⁺-dependent structural and dynamic distinctions between the apo and holo forms of the protein. Furthermore, the process of ion binding by inserting Ca²⁺ into the bulk of the apo structure was simulated by molecular dynamics. Expectations of the simulation were demonstrated using cluster analysis and a variety of structural metrics, such as interhelical angle estimation, solvent accessible surface area, hydrogen bond analysis, and contact analysis. Ca²⁺ triggered a rise in the interhelical angles of S100A1 on the binding site and solvent accessible surface area. Significant configurational regulations were observed in the holo protein. The findings would contribute to understanding the molecular basis of the association of Ca²⁺ with the S100A1 protein, which may be an appropriate study to understand the Ca²⁺-mediated conformational changes in the protein target. In addition, we investigated the expression profile of S100A1 in myoblast differentiation and muscle regeneration. These data showed that S100A1 is expressed in skeletal muscles. However, the expression decreases with time during the process of myoblast differentiation.

Magnesium and calcium reveal different chelating effects in a steroid compound: A model study of prednisolone using NMR spectroscopy

In this work, we used high resolution NMR spectroscopy to investigate metal cation chelation by the steroidal drug 1,4-pregnadiene-11β,17α,21-triol-3,20-dione (Prednisolone; abbreviated as Prd). Prd/MgCl₂ and Prd/CaCl₂ mixtures were prepared at eight different molar ratios. Using two-dimensional ¹H/¹³C heteronuclear correlation spectroscopy, we were able to resolve most of the ¹H signals, except those at 1.4-1.55 ppm, where signals for H15β, H16α and Me-19 are superimposed. The chelation sites were determined by the cation concentration-dependent ¹³C signals. Both ring A and ring D of Prd were found to be involved in Mg²⁺ chelation, whereas only ring A was involved in Ca²⁺ chelation. The dihedral angles deduced from the ³J_H-H coupling constants indicated that ring D of Prd might undergo rather small, but different, distortions in the presence of Mg²⁺ and Ca²⁺. Additionally, using the continuous variation method, we deduced that the stoichiometric ratios of the Prd/Mg²⁺ and Prd/Ca²⁺ complexes were 1:1 and 2:1, respectively. All of the evidence led us to conclude that the Prd/Mg²⁺ and Prd/Ca²⁺ complexes are mediated by different chelating mechanisms.

Structural Biology and Molecular Modeling to Analyze the Entry of Bacterial Toxins and Virulence Factors into Host Cells

Understanding the functions and mechanisms of biological systems is an outstanding challenge. One way to overcome it is to combine together several approaches such as molecular modeling and experimental structural biology techniques. Indeed, the interplay between structural and dynamical properties of the system is crucial to unravel the function of molecular machinery’s. In this review, we focus on how molecular simulations along with structural information can aid in interpreting biological data. Here, we examine two different cases: (i) the endosomal translocation toxins (diphtheria, tetanus, botulinum toxins) and (ii) the activation of adenylyl cyclase inside the cytoplasm (edema factor, CyA, ExoY).

How calcium ion binding induces the conformational transition of the calmodulin N-terminal domain—an atomic level characterization

The Ca²⁺binding and triggering conformation transition of nCaM were detected in unbiased molecular dynamics simulations.

Structural dynamics of calmodulin-ryanodine receptor interactions: electron paramagnetic resonance using stereospecific spin labels

We have used electron paramagnetic resonance, with rigid and stereospecific spin labels, to resolve structural states in calmodulin (CaM), as affected by binding of Ca and a CaM-binding peptide (RyRp) derived from the ryanodine receptor (RyR), the Ca channel that triggers muscle contraction. CaM mutants containing a pair of cysteines in the N-lobe and/or C-lobe were engineered and labeled with a stereospecifically bound bifunctional spin label (BSL). RyRp was synthesized with and without TOAC (a stereospecifically attached spin-labeled amino acid) substituted for a single amino acid near the N-terminus. Intramolecular DEER distance measurements of doubly-labeled BSL-CaM revealed that CaM exists in dynamic equilibrium among multiple states, consistent with open, closed, and compact structural models. Addition of RyRp shifted the equilibrium partially toward the compact state in the absence of Ca, and completely toward the compact state in the presence of Ca, supporting a conformational selection model. Inter-protein distance measurements show that Ca stabilizes the compact state primarily by inducing ordered binding of the CaM N-lobe to RyRp, while only slightly affecting the C-lobe. The results provide insight into the structural mechanism of CaM-mediated RyR regulation, while demonstrating the power of using two types of rigidly and stereospecifically bound spin labels.

Cyanylated Cysteine Reports Site-Specific Changes at Protein-Protein-Binding Interfaces Without Perturbation

To investigate the cyanylated cysteine vibrational probe group’s ability to report on binding-induced changes along a protein–protein interface, the probe group was incorporated at several sites in a peptide of the calmodulin (CaM)-binding domain of skeletal muscle myosin light chain kinase. Isothermal titration calorimetry was used to determine the binding thermodynamics between calmodulin and each peptide. For all probe positions, the binding affinity was nearly identical to that of the unlabeled peptide. The CN stretching infrared band was collected for each peptide free in solution and bound to calmodulin. Binding-induced shifts in the IR spectral frequencies were correlated with estimated solvent accessibility based on molecular dynamics simulations. This work generally suggests (1) that site-specific incorporation of this vibrational probe group does not cause major perturbations to its local structural environment and (2) that this small probe group might be used quite broadly to map dynamic protein-binding interfaces. However, site-specific perturbations due to artificial labeling groups can be somewhat unpredictable and should be evaluated on a site-by-site basis through complementary measurements. A fully quantitative, simulation-based interpretation of the rich probe IR spectra is still needed but appears to be possible given recent advances in simulation techniques.

Probing Activated and Non-Activated Single Calmodulin Molecules under a Piconewton Compressive Force

Interrogating the protein structure-function inter-relationship under a piconewton force manipulation has been highly promising and informative. Although protein conformational changes under pulling force manipulations have been extensively studied, protein conformational changes under a compressive force have not been explored in detail. Using our home-modified sensitive and high signal-to-noise atomic force microscopy (AFM) approach, we have applied a piconewton compressive force, manipulating a Calmodulin (CaM) molecule to characterize two different forms of CaM, the Ca²⁺-ligated activated form and the Ca²⁺ free non-activated form (apo-CaM). We observed sudden and spontaneous structural rupture of apo-CaM under compressive force applied by an AFM tip, though no such events were recorded in the case of Ca²⁺-ligated activated CaM form. The sudden spontaneous structural rupture under a piconewton force compression has never been reported before, which presents an unexplored function that is likely important for protein-protein interactions and cell signaling functions.

Coordination to lanthanide ions distorts binding site conformation in calmodulin

The Ca²⁺-sensing protein calmodulin (CaM) is a popular model of biological ion binding since it is both experimentally tractable and essential to survival in all eukaryotic cells. CaM modulates hundreds of target proteins and is sensitive to complex patterns of Ca²⁺ exposure, indicating that it functions as a sophisticated dynamic transducer rather than a simple on/off switch. Many details of this transduction function are not well understood. Fourier transform infrared (FTIR) spectroscopy, ultrafast 2D infrared (2D IR) spectroscopy, and electronic structure calculations were used to probe interactions between bound metal ions (Ca²⁺ and several trivalent lanthanide ions) and the carboxylate groups in CaM's EF-hand ion-coordinating sites. Since Tb³⁺ is commonly used as a luminescent Ca²⁺ analog in studies of protein-ion binding, it is important to characterize distinctions between the coordination of Ca²⁺ and the lanthanides in CaM. Although functional assays indicate that Tb³⁺ fully activates many Ca²⁺-dependent proteins, our FTIR spectra indicate that Tb³⁺, La³⁺, and Lu³⁺ disrupt the bidentate coordination geometry characteristic of the CaM binding sites' strongly conserved position 12 glutamate residue. The 2D IR spectra indicate that, relative to the Ca²⁺-bound form, lanthanide-bound CaM exhibits greater conformational flexibility and larger structural fluctuations within its binding sites. Time-dependent 2D IR lineshapes indicate that binding sites in Ca²⁺-CaM occupy well-defined configurations, whereas binding sites in lanthanide-bound-CaM are more disordered. Overall, the results show that binding to lanthanide ions significantly alters the conformation and dynamics of CaM's binding sites.

Bilobal architecture is a requirement for calmodulin signaling to CaV1.3 channels

Calmodulin (CaM) regulation of voltage-gated calcium (Ca_V) channels is a powerful Ca²⁺ feedback mechanism that adjusts Ca²⁺ influx, affording rich mechanistic insights into Ca²⁺ decoding. CaM possesses a dual-lobed architecture, a salient feature of the myriad Ca²⁺-sensing proteins, where two homologous lobes that recognize similar targets hint at redundant signaling mechanisms. Here, by tethering CaM lobes, we demonstrate that bilobal architecture is obligatory for signaling to Ca_V channels. With one lobe bound, Ca_V carboxy tail rearranges itself, resulting in a preinhibited configuration precluded from Ca²⁺ feedback. Reconstitution of two lobes, even as separate molecules, relieves preinhibition and restores Ca²⁺ feedback. Ca_V channels thus detect the coincident binding of two Ca²⁺-free lobes to promote channel opening, a molecular implementation of a logical NOR operation that processes spatiotemporal Ca²⁺ signals bifurcated by CaM lobes. Overall, a unified scheme of Ca_V channel regulation by CaM now emerges, and our findings highlight the versatility of CaM to perform exquisite Ca²⁺ computations.

Conformational Ensembles of Calmodulin Revealed by Nonperturbing Site-Specific Vibrational Probe Groups

Seven native residues on the regulatory protein calmodulin, including three key methionine residues, were replaced (one by one) by the vibrational probe amino acid cyanylated cysteine, which has a unique CN stretching vibration that reports on its local environment. Almost no perturbation was caused by this probe at any of the seven sites, as reported by CD spectra of calcium-bound and apo calmodulin and binding thermodynamics for the formation of a complex between calmodulin and a canonical target peptide from skeletal muscle myosin light chain kinase measured by isothermal titration. The surprising lack of perturbation suggests that this probe group could be applied directly in many protein–protein binding interfaces. The infrared absorption bands for the probe groups reported many dramatic changes in the probes’ local environments as CaM went from apo- to calcium-saturated to target peptide-bound conditions, including large frequency shifts and a variety of line shapes from narrow (interpreted as a rigid and invariant local environment) to symmetric to broad and asymmetric (likely from multiple coexisting and dynamically exchanging structures). The fast intrinsic time scale of infrared spectroscopy means that the line shapes report directly on site-specific details of calmodulin’s variable structural distribution. Though quantitative interpretation of the probe line shapes depends on a direct connection between simulated ensembles and experimental data that does not yet exist, formation of such a connection to data such as that reported here would provide a new way to evaluate conformational ensembles from data that directly contains the structural distribution. The calmodulin probe sites developed here will also be useful in evaluating the binding mode of calmodulin with many uncharacterized regulatory targets.

Using single molecule force spectroscopy to facilitate a rational design of Ca2+-responsive β-roll peptide-based hydrogels

Mechanical stability of Ca²⁺-responsive β-roll peptides (RTX) is largely responsible for the Ca²⁺-dependent mechanical properties of the RTX-based hydrogels.

Arabidopsis Calmodulin-Like Proteins, CML15 and CML16 Possess Biochemical Properties Distinct from Calmodulin and Show Non-overlapping Tissue Expression Patterns

Calcium ions are used as ubiquitous, key second messengers in cells across eukaryotic taxa. In plants, calcium signal transduction is involved in a wide range of cellular processes from abiotic and biotic stress responses to development and growth. Calcium signals are detected by calcium sensor proteins, of which calmodulin (CaM), is the most evolutionarily conserved and well-studied. These sensors regulate downstream targets to propagate the information in signaling pathways. Plants possess a large family of calcium sensors related to CaM, termed CaM-like (CMLs), that are not found in animals and remain largely unstudied at the structural and functional level. Here, we investigated the biochemical properties and gene promoter activity of two closely related members of the Arabidopsis CML family, CML15 and CML16. Biochemical characterization of recombinant CML15 and CML16 indicated that they possess properties consistent with their predicted roles as calcium sensors. In the absence of calcium, CML15 and CML16 display greater intrinsic hydrophobicity than CaM. Both CMLs displayed calcium-dependent and magnesium-independent conformational changes that expose hydrophobic residues, but the degree of hydrophobic exposure was markedly less than that observed for CaM. Isothermal titration calorimetry indicated two and three calcium-binding sites for CML15 and CML16, respectively, with affinities expected to be within a physiological range. Both CML15 and CML16 bound calcium with high affinity in the presence of excess magnesium. Promoter-reporter analysis demonstrated that the CML16 promoter is active across a range of Arabidopsis tissues and developmental stages, whereas the CML15 promoter activity is very restricted and was observed only in floral tissues, specifically anthers and pollen. Collectively, our data indicate that these CMLs behave biochemically like calcium sensors but with properties distinct from CaM and likely have non-overlapping roles in floral development. We discuss our findings in the broader context of calcium sensors and signaling in Arabidopsis.

Biophysical characterization of the calmodulin-like domain of Plasmodium falciparum calcium dependent protein kinase 3

Calcium dependent protein kinases are unique to plants and certain parasites and comprise an N-terminal segment and a kinase domain that is regulated by a C-terminal calcium binding domain. Since the proteins are not found in man they are potential drug targets. We have characterized the calcium binding lobes of the regulatory domain of calcium dependent protein kinase 3 from the malaria parasite Plasmodium falciparum. Despite being structurally similar, the two lobes differ in several other regards. While the monomeric N-terminal lobe changes its structure in response to calcium binding and shows global dynamics on the sub-millisecond time-scale both in its apo and calcium bound states, the C-terminal lobe could not be prepared calcium-free and forms dimers in solution. If our results can be generalized to the full-length protein, they suggest that the C-terminal lobe is calcium bound even at basal levels and that activation is caused by the structural reorganization associated with binding of a single calcium ion to the N-terminal lobe.

Nanopore-Based Measurements of Protein Size, Fluctuations, and Conformational Changes

Proteins are structurally dynamic macromolecules, and it is challenging to quantify the conformational properties of their native state in solution. Nanopores can be efficient tools to study proteins in a solution environment. In this method, an electric field induces electrophoretic and/or electro-osmotic transport of protein molecules through a nanopore slightly larger than the protein molecule. High-bandwidth ion current measurement is used to detect the transit of each protein molecule. First, our measurements reveal a correlation between the mean current blockade amplitude and the radius of gyration for each protein. Next, we find a correlation between the shape of the current signal amplitude distributions and the protein fluctuation as obtained from molecular dynamics simulations. Further, the magnitude of the structural fluctuations, as probed by experiments and simulations, correlates with the ratio of α-helix to β-sheet content. We highlight the resolution of our measurements by resolving two states of calmodulin, a canonical protein that undergoes a conformational change in response to calcium binding.

Molecular Modeling of the Catalytic Domain of CyaA Deepened the Knowledge of Its Functional Dynamics

Although CyaA has been studied for over three decades and revealed itself to be a very good prototype for developing various biotechnological applications, only a little is known about its functional dynamics and about the conformational landscape of this protein. Molecular dynamics simulations helped to clarify the view on these points in the following way. First, the model of interaction between AC and calmodulin (CaM) has evolved from an interaction centered on the surface between C-CaM hydrophobic patch and the α helix H of AC, to a more balanced view, in which the C-terminal tail of AC along with the C-CaM Calcium loops play an important role. This role has been confirmed by the reduction of the affinity of AC for calmodulin in the presence of R338, D360 and N347 mutations. In addition, enhanced sampling studies have permitted to propose a representation of the conformational space for the isolated AC. It remains to refine this representation using structural low resolution information measured on the inactive state of AC. Finally, due to a virtual screening study on another adenyl cyclase from Bacillus anthracis, weak inhibitors of AC have been discovered.

Incorporating an allosteric regulatory site in an antibody through backbone design

Allosteric regulation underlies living cells' ability to sense changes in nutrient and signaling-molecule concentrations, but the ability to computationally design allosteric regulation into non-allosteric proteins has been elusive. Allosteric-site design is complicated by the requirement to encode the relative stabilities of active and inactive conformations of the same protein in the presence and absence of both ligand and effector. To address this challenge, we used Rosetta to design the backbone of the flexible heavy-chain complementarity-determining region 3 (HCDR3), and used geometric matching and sequence optimization to place a Zn²⁺ -coordination site in a fluorescein-binding antibody. We predicted that due to HCDR3's flexibility, the fluorescein-binding pocket would configure properly only upon Zn²⁺ application. We found that regulation by Zn²⁺ was reversible and sensitive to the divalent ion's identity, and came at the cost of reduced antibody stability and fluorescein-binding affinity. Fluorescein bound at an order of magnitude higher affinity in the presence of Zn²⁺ than in its absence, and the increase in fluorescein affinity was due almost entirely to faster fluorescein on-rate, suggesting that Zn²⁺ preorganized the antibody for fluorescein binding. Mutation analysis demonstrated the extreme sensitivity of Zn²⁺ regulation on the atomic details in and around the metal-coordination site. The designed antibody could serve to study how allosteric regulation evolved from non-allosteric binding proteins, and suggests a way to designing molecular sensors for environmental and biomedical targets.

Computational comparison of a calcium-dependent jellyfish protein (apoaequorin) and calmodulin-cholesterol in short-term memory maintenance

Memory reconsolidation and maintenance depend on calcium channels and on calcium/calmodulin-dependent kinases regulating protein turnover in the hippocampus. Ingestion of a jellyfish protein, apoaequorin, reportedly protects and/or improves verbal learning in adults and is currently widely advertised for use by the elderly. Apoaequorin is a member of the EF-hand calcium binding family of proteins that includes calmodulin. Calmodulin-1 (148 residues) differs from Apoaequorin (195 residues) in that it contains four rather than three Ca²⁺-binding sites and three rather than four cholesterol-binding (CRAC, CARC) domains. All three cholesterol-binding CARC domains in calmodulin have a high interaction affinity for cholesterol compared to only two high affinity CARC domains in apoaequorin. Both calmodulin and apoaequorin can form dimers with a potential of eight bound Ca²⁺ ions and six high affinity-bound cholesterol molecules in calmodulin with six bound Ca²⁺ ions and a mixed population of eight cholesterols bound to both CARC and CRAC domains in apoaqueorin. MEMSAT-SVM analysis indicates that both calmodulin and apoaqueorin have a pore-lining region. The Peptide-Cutter algorithm predicts that calmodulin-1 contains 11 trypsin-specific cleavage sites (compared to 21 in apoaqueorin), four of which are potentially blocked by cholesterol and three are within the Ca-binding domains and/or the pore-lining region. Three are clustered between the third and fourth Ca²⁺-binding sites. Only calmodulin pore-lining regions contain Ca²⁺ binding sites and as dimers may insert into the plasma membrane of neural cells and act as Ca²⁺ channels. In a dietary supplement, bound cholesterol may protect both apoaequorin and calmodulin from proteolysis in the gut as well as facilitate uptake across the blood-brain barrier. Our results suggest that a physiological calmodulin-cholesterol complex, not cholesterol-free jellyfish protein, may better serve as a dietary supplement to facilitate memory maintenance.

Molecular Basis of S100A1 Activation at Saturating and Subsaturating Calcium Concentrations

The S100A1 protein mediates a wide variety of physiological processes through its binding of calcium (Ca(2+)) and endogenous target proteins. S100A1 presents two Ca(2+)-binding domains: a high-affinity "canonical" EF (cEF) hand and a low-affinity "pseudo" EF (pEF) hand. Accumulating evidence suggests that both Ca(2+)-binding sites must be saturated to stabilize an open state conducive to peptide recognition, yet the pEF hand's low affinity limits Ca(2+) binding at normal physiological concentrations. To understand the molecular basis of Ca(2+) binding and open-state stabilization, we performed 100 ns molecular dynamics simulations of S100A1 in the apo/holo (Ca(2+)-free/bound) states and a half-saturated state, for which only the cEF sites are Ca(2+)-bound. Our simulations indicate that the pattern of oxygen coordination about Ca(2+) in the cEF relative to the pEF site contributes to the former's higher affinity, whereas Ca(2+) binding strongly reshapes the protein's conformational dynamics by disrupting β-sheet coupling between EF hands. Moreover, modeling of the half-saturated configuration suggests that the open state is unstable and reverts toward a closed state in the absence of the pEF Ca(2+) ion. These findings indicate that Ca(2+) binding at the cEF site alone is insufficient to stabilize opening; thus, posttranslational modification of the protein may be required for target peptide binding at subsaturating intracellular Ca(2+) levels.

Single-Molecule FRET States, Conformational Interchange, and Conformational Selection by Dye Labels in Calmodulin

We investigate the roles of measurement time scale and the nature of the fluorophores in the FRET states measured for calmodulin, a calcium signaling protein known to undergo pronounced conformational changes. The measured FRET distributions depend markedly on the measurement time scale (nanosecond or microsecond). Comparison of FRET distributions measured by donor fluorescence decay with FRET distributions recovered from single-molecule burst measurements binned over time scales of 90 μs to 1 ms reveals conformational averaging over the intervening time regimes. We find further that, particularly in the presence of saturating Ca(2+), the nature of the measured single-molecule FRET distribution depends markedly on the identity of the FRET pair. The results suggest interchange between conformational states on time scales of hundreds of microseconds or less. Interaction with a fluorophore such as the dye Texas Red alters both the nature of the measured FRET distributions and the dynamics of conformational interchange. The results further suggest that the fluorophore may not be merely a benign reporter of protein conformations in FRET studies, but may in fact alter the conformational landscape.

Conformational heterogeneity of the calmodulin binding interface

Calmodulin (CaM) is a ubiquitous Ca(2+) sensor and a crucial signalling hub in many pathways aberrantly activated in disease. However, the mechanistic basis of its ability to bind diverse signalling molecules including G-protein-coupled receptors, ion channels and kinases remains poorly understood. Here we harness the high resolution of molecular dynamics simulations and the analytical power of Markov state models to dissect the molecular underpinnings of CaM binding diversity. Our computational model indicates that in the absence of Ca(2+), sub-states in the folded ensemble of CaM's C-terminal domain present chemically and sterically distinct topologies that may facilitate conformational selection. Furthermore, we find that local unfolding is off-pathway for the exchange process relevant for peptide binding, in contrast to prior hypotheses that unfolding might account for binding diversity. Finally, our model predicts a novel binding interface that is well-populated in the Ca(2+)-bound regime and, thus, a candidate for pharmacological intervention.

Ca-Dependent Folding of Human Calumenin

Human calumenin (hCALU) is a six EF-hand protein belonging to the CREC family. As other members of the family, it is localized in the secretory pathway and regulates the activity of SERCA2a and of the ryanodine receptor in the endoplasmic reticulum (ER). We have studied the effects of Ca²⁺ binding to the protein and found it to attain a more compact structure upon ion binding. Circular Dichroism (CD) measurements suggest a major rearrangement of the protein secondary structure, which reversibly switches from disordered at low Ca²⁺ concentrations to predominantly alpha-helical when Ca²⁺ is added. SAXS experiments confirm the transition from an unfolded to a compact structure, which matches the structural prediction of a trilobal fold. Overall our experiments suggest that calumenin is a Ca²⁺ sensor, which folds into a compact structure, capable of interacting with its molecular partners, when Ca²⁺ concentration within the ER reaches the millimolar range.

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.

Calmodulin and STIM proteins: Two major calcium sensors in the cytoplasm and endoplasmic reticulum

The calcium (Ca(2+)) ion is a universal signalling messenger which plays vital physiological roles in all eukaryotes. To decode highly regulated intracellular Ca(2+) signals, cells have evolved a number of sensor proteins that are ideally adapted to respond to a specific range of Ca(2+) levels. Among many such proteins, calmodulin (CaM) is a multi-functional cytoplasmic Ca(2+) sensor with a remarkable ability to interact with and regulate a plethora of structurally diverse target proteins. CaM achieves this 'multi-talented' functionality through two EF-hand domains, each with an independent capacity to bind targets, and an adaptable flexible linker. By contrast, stromal interaction molecule-1 and -2 (STIMs) have evolved for a specific role in endoplasmic reticulum (ER) Ca(2+) sensing using EF-hand machinery analogous to CaM; however, whereas CaM structurally adjusts to dissimilar binding partners, STIMs use the EF-hand machinery to self-regulate the stability of the Ca(2+) sensing domain. The molecular mechanisms underlying the Ca(2+)-dependent signal transduction by CaM and STIMs have revealed a remarkable repertoire of actions and underscore the flexibility of nature in molecular evolution and adaption to discrete Ca(2+) levels. Recent genomic sequencing efforts have uncovered a number of disease-associated mutations in both CaM and STIM1. This article aims to highlight the most recent key structural and functional findings in the CaM and STIM fields, and discusses how these two Ca(2+) sensor proteins execute their biological functions.

Calcium binding properties of calcium dependent protein kinase 1 (CaCDPK1) from Cicer arietinum

Calcium plays a crucial role as a secondary messenger in all aspects of plant growth, development and survival. Calcium dependent protein kinases (CDPKs) are the major calcium decoders, which couple the changes in calcium level to an appropriate physiological response. The mechanism by which calcium regulates CDPK protein is not well understood. In this study, we investigated the interactions of Ca(2+) ions with the CDPK1 isoform of Cicer arietinum (CaCDPK1) using a combination of biophysical tools. CaCDPK1 has four different EF hands as predicted by protein sequence analysis. The fluorescence emission spectrum of CaCDPK1 showed quenching with a 5 nm red shift upon addition of calcium, indicating conformational changes in the tertiary structure. The plot of changes in intensity against calcium concentrations showed a biphasic curve with binding constants of 1.29 μM and 120 μM indicating two kinds of binding sites. Isothermal calorimetric (ITC) titration with CaCl2 also showed a biphasic curve with two binding constants of 0.027 μM and 1.7 μM. Circular dichroism (CD) spectra showed two prominent peaks at 208 and 222 nm indicating that CaCDPK1 is a α-helical rich protein. Calcium binding further increased the α-helical content of CaCDPK1 from 75 to 81%. Addition of calcium to CaCDPK1 also increased fluorescence of 8-anilinonaphthalene-1-sulfonic acid (ANS) indicating exposure of hydrophobic surfaces. Thus, on the whole this study provides evidence for calcium induced conformational changes, exposure of hydrophobic surfaces and heterogeneity of EF hands in CaCDPK1.

Allosteric activation of Bordetella pertussis adenylyl cyclase by calmodulin: molecular dynamics and mutagenesis studies

Adenylyl cyclase (AC) toxin is an essential toxin that allows Bordetella pertussis to invade eukaryotic cells, where it is activated after binding to calmodulin (CaM). Based on the crystal structure of the AC catalytic domain in complex with the C-terminal half of CaM (C-CaM), our previous molecular dynamics simulations (Selwa, E., Laine, E., and Malliavin, T. (2012) Differential role of calmodulin and calcium ions in the stabilization of the catalytic domain of adenyl cyclase CyaA from Bordetella pertussis. Proteins 80, 1028–1040) suggested that three residues (i.e. Arg(338), Asn(347), and Asp(360)) might be important for stabilizing the AC/CaM interaction. These residues belong to a loop-helix-loop motif at the C-terminal end of AC, which is located at the interface between CaM and the AC catalytic loop. In the present study, we conducted the in silico and in vitro characterization of three AC variants, where one (Asn(347); ACm1A), two (Arg(338) and Asp(360); ACm2A), or three residues (Arg(338), Asn(347), and Asp(360); ACm3A) were substituted with Ala. Biochemical studies showed that the affinities of ACm1A and ACm2A for CaM were not affected significantly, whereas that of ACm3A was reduced dramatically. To understand the effects of these modifications, molecular dynamics simulations were performed based on the modified proteins. The molecular dynamics trajectories recorded for the ACm3AC-CaM complex showed that the calcium-binding loops of C-CaM exhibited large fluctuations, which could be related to the weakened interaction between ACm3A and its activator. Overall, our results suggest that the loop-helix-loop motif at the C-terminal end of AC is crucial during CaM binding for stabilizing the AC catalytic loop in an active configuration.

Modulation of calmodulin lobes by different targets: an allosteric model with hemiconcerted conformational transitions

Calmodulin is a calcium-binding protein ubiquitous in eukaryotic cells, involved in numerous calcium-regulated biological phenomena, such as synaptic plasticity, muscle contraction, cell cycle, and circadian rhythms. It exibits a characteristic dumbell shape, with two globular domains (N- and C-terminal lobe) joined by a linker region. Each lobe can take alternative conformations, affected by the binding of calcium and target proteins. Calmodulin displays considerable functional flexibility due to its capability to bind different targets, often in a tissue-specific fashion. In various specific physiological environments (e.g. skeletal muscle, neuron dendritic spines) several targets compete for the same calmodulin pool, regulating its availability and affinity for calcium. In this work, we sought to understand the general principles underlying calmodulin modulation by different target proteins, and to account for simultaneous effects of multiple competing targets, thus enabling a more realistic simulation of calmodulin-dependent pathways. We built a mechanistic allosteric model of calmodulin, based on an hemiconcerted framework: each calmodulin lobe can exist in two conformations in thermodynamic equilibrium, with different affinities for calcium and different affinities for each target. Each lobe was allowed to switch conformation on its own. The model was parameterised and validated against experimental data from the literature. In spite of its simplicity, a two-state allosteric model was able to satisfactorily represent several sets of experiments, in particular the binding of calcium on intact and truncated calmodulin and the effect of different skMLCK peptides on calmodulin's saturation curve. The model can also be readily extended to include multiple targets. We show that some targets stabilise the low calcium affinity T state while others stabilise the high affinity R state. Most of the effects produced by calmodulin targets can be explained as modulation of a pre-existing dynamic equilibrium between different conformations of calmodulin's lobes, in agreement with linkage theory and MWC-type models.

Conformational frustration in calmodulin-target recognition

Calmodulin (CaM) is a primary calcium (Ca(2+) )-signaling protein that specifically recognizes and activates highly diverse target proteins. We explored the molecular basis of target recognition of CaM with peptides representing the CaM-binding domains from two Ca(2+) -CaM-dependent kinases, CaMKI and CaMKII, by employing experimentally constrained molecular simulations. Detailed binding route analysis revealed that the two CaM target peptides, although similar in length and net charge, follow distinct routes that lead to a higher binding frustration in the CaM-CaMKII complex than in the CaM-CaMKI complex. We discovered that the molecular origin of the binding frustration is caused by intermolecular contacts formed with the C-domain of CaM that need to be broken before the formation of intermolecular contacts with the N-domain of CaM. We argue that the binding frustration is important for determining the kinetics of the recognition process of proteins involving large structural fluctuations.

Electron paramagnetic resonance spectroscopy of nitroxide-labeled calmodulin

Calmodulin (CaM) is a highly conserved calcium-binding protein consisting of two homologous domains, each of which contains two EF-hands, that is known to bind well over 300 proteins and peptides. In most cases the (Ca(2+))(4-)form of CaM leads to the activation of a key regulatory enzyme or protein in a myriad of biological processes. Using the nitroxide spin-labeling reagent, 3-(2-iodoacetamido)-2,2,5,5-tetramethyl-1-pyrrolidinyl oxyl, bovine brain CaM was modified at 2-3 methionines with retention of activity as judged by the activation of cyclic nucleotide phosphodiesterase. X-band electron paramagnetic resonance (EPR) spectroscopy was used to measure the spectral changes upon addition of Ca(2+) to the apo-form of spin-labeled protein. A significant loss of spectral intensity, arising primarily from reductions in the heights of the low, intermediate, and high field peaks, accompanied Ca(2+) binding. The midpoint of the Ca(2+)-mediated transition determined by EPR occurred at a higher Ca(2+) concentration than that measured with circular dichroic spectroscopy and enzyme activation. Recent data have indicated that the transition from the apo-state of CaM to the fully saturated form, [(Ca(2+))(4-)CaM], contains a compact intermediate corresponding to [(Ca(2+))(2-)CaM], and the present results suggest that the spin probes are reporting on Ca(2+) binding to the last two sites in the N-terminal domain, i.e. for the [(Ca(2+))(2)-CaM] → [(Ca(2+))(4-)CaM] transition in which the compact structure becomes more extended. EPR of CaM, spin-labeled at methionines, offers a different approach for studying Ca(2+)-mediated conformational changes and may emerge as a useful technique for monitoring interactions with target proteins.

Structure of myosin-1c tail bound to calmodulin provides insights into calcium-mediated conformational coupling

Class I myosins can sense cellular mechanical forces and function as tension-sensitive anchors or transporters. How mechanical load is transduced from the membrane-binding tail to the force-generating head in myosin-1 is unknown. Here we determined the crystal structure of the entire tail of mouse myosin-1c in complex with apocalmodulin, showing that myosin-1c adopts a stable monomer conformation suited for force transduction. The lever-arm helix and the C-terminal extended PH domain of the motor are coupled by a stable post-IQ domain bound to calmodulin in a highly unusual mode. Ca(2+) binding to calmodulin induces major conformational changes in both IQ motifs and the post-IQ domain and increases flexibility of the myosin-1c tail. Our study provides a structural blueprint for the neck and tail domains of myosin-1 and expands the target binding modes of the master Ca(2+)-signal regulator calmodulin.

Characterization of calmodulin-Fas death domain interaction: an integrated experimental and computational study

The Fas death receptor-activated death-inducing signaling complex (DISC) regulates apoptosis in many normal and cancer cells. Qualitative biochemical experiments demonstrate that calmodulin (CaM) binds to the death domain of Fas. The interaction between CaM and Fas regulates Fas-mediated DISC formation. A quantitative understanding of the interaction between CaM and Fas is important for the optimal design of antagonists for CaM or Fas to regulate the CaM-Fas interaction, thus modulating Fas-mediated DISC formation and apoptosis. The V254N mutation of the Fas death domain (Fas DD) is analogous to an identified mutant allele of Fas in lpr-cg mice that have a deficiency in Fas-mediated apoptosis. In this study, the interactions of CaM with the Fas DD wild type (Fas DD WT) and with the Fas DD V254N mutant were characterized using isothermal titration calorimetry (ITC), circular dichroism spectroscopy (CD), and molecular dynamics (MD) simulations. ITC results reveal an endothermic binding characteristic and an entropy-driven interaction of CaM with Fas DD WT or with Fas DD V254N. The Fas DD V254N mutation decreased the association constant (Ka) for CaM-Fas DD binding from (1.79 ± 0.20) × 10(6) to (0.88 ± 0.14) × 10(6) M(-1) and slightly increased a standard state Gibbs free energy (ΔG°) for CaM-Fas DD binding from -8.87 ± 0.07 to -8.43 ± 0.10 kcal/mol. CD secondary structure analysis and MD simulation results did not show significant secondary structural changes of the Fas DD caused by the V254N mutation. The conformational and dynamical motion analyses, the analyses of hydrogen bond formation within the CaM binding region, the contact numbers of each residue, and the electrostatic potential for the CaM binding region based on MD simulations demonstrated changes caused by the Fas DD V254N mutation. These changes caused by the Fas DD V254N mutation could affect the van der Waals interactions and electrostatic interactions between CaM and Fas DD, thereby affecting CaM-Fas DD interactions. Results from this study characterize CaM-Fas DD interactions in a quantitative way, providing structural and thermodynamic evidence of the role of the Fas DD V254N mutation in the CaM-Fas DD interaction. Furthermore, the results could help to identify novel strategies for regulating CaM-Fas DD interactions and Fas DD conformation and thus to modulate Fas-mediated DISC formation and thus Fas-mediated apoptosis.

Effect of the Brugada syndrome mutation A39V on calmodulin regulation of Cav1.2 channels

Background

The L-type calcium channel Cav1.2 is important for brain and heart function. The ubiquitous calcium sensing protein calmodulin (CaM) regulates calcium dependent gating of Cav1.2 channels by reducing calcium influx, a process known as calcium-dependent inactivation (CDI). Dissecting the calcium-dependence of CaM in this process has benefited greatly from the use of mutant CaM molecules which are unable to bind calcium to their low affinity (N-lobe) and high affinity (C-lobe) binding sites. Unlike CDI, it is unknown whether CaM can modulate the activation gating of Cav1.2 channels.

Results

We examined a Cav1.2 point mutant in the N-terminus region of the channel (A39V) that has been previously linked to Brugada syndrome. Using mutant CaM constructs in which the N- and/or C-lobe calcium binding sites were ablated, we were able to show that this Brugada syndrome mutation disrupts N-lobe CDI of the channel. In the course of these experiments, we discovered that all mutant CaM molecules were able to alter the kinetics of channel activation even in the absence of calcium for WT-Cav1.2, but not A39V-Cav1.2 channels. Moreover, CaM mutants differentially shifted the voltage-dependence of activation for WT and A39V-Cav1.2 channels to hyperpolarized potentials. Our data therefore suggest that structural changes in CaM that arise directly from site directed mutagenesis of calcium binding domains alter activation gating of Cav1.2 channels independently of their effects on calcium binding, and that the N-terminus of the channel contributes to this CaM dependent process.

Conclusions

Our data indicate that caution must be exercised when interpreting the effects of CaM mutants on ion channel gating.

The calmodulin-like protein CML43 functions as a salicylic-acid-inducible root-specific Ca(2+) sensor in Arabidopsis

Many signalling pathways in plants are regulated by the second messenger calcium (Ca(2+)). In the standard model, Ca(2+)-sensor proteins, such as CaM (calmodulin), detect Ca(2+) signals and subsequently regulate downstream targets to advance the signal transduction cascade. In addition to CaM, plants possess many CMLs (CaM-like proteins) that are predicted to function as Ca(2+) sensors, but which remain largely uncharacterized. In the present study, we examined the biochemical properties, subcellular localization and tissue-specific distribution of Arabidopsis CML43. Our data indicate that CML43 displays characteristics typical of Ca(2+) sensors, including high-affinity Ca(2+) binding, conformational changes upon Ca(2+) binding that expose hydrophobic regions and stabilization of structure in the presence of Mg(2+) or Ca(2+). In vivo localization analysis demonstrates that CML43 resides in cytosolic and nuclear compartments. Transgenic plants expressing a CML43:GUS (β-glucoronidase) promoter reporter gene revealed that CML43 promoter activity is restricted almost exclusively to root tips under normal growth conditions. GUS reporter activity in these transgenic plants was strongly increased when exposed to the defence compound SA (salicylic acid). Furthermore, immunoblot analysis revealed that the CML43 protein accumulates following treatment with SA. Collectively, our findings suggest that CML43 functions as a Ca(2+) sensor in root tips during both normal growth and plant immune response.

Extramembrane control of ion channel peptide assemblies, using alamethicin as an example

Ion channels allow the influx and efflux of specific ions through a plasma membrane. Many ion channels can sense, for example, the membrane potential (the voltage gaps between the inside and the outside of the membrane), specific ligands such as neurotransmitters, and mechanical tension within the membrane. They modulate cell function in response to these stimuli. Researchers have focused on developing peptide- and non-peptide-based model systems to elucidate ion-channel protein functions and to create artificial sensing systems. In this Account, we employed a typical peptide that forms ion channels,alamethicin, as a model to evaluate our methodologies for controlling the assembly states of channel-forming molecules in membranes. As alamethicin self-assembles in membranes, it prompts channel formation, but number of peptide molecules in these channels is not constant. Using planar-lipid bilayer methods, we monitored the association states of alamethicin in real time. Many ligand-gated, natural-ion channel proteins have large extramembrane domains. As these proteins interact with specific ligands, those conformational alterations in the extramembrane domains are transmitted to the transmembrane, pore-forming domains to open and close the channels. We hypothesized that if we conjugated suitable extramembrane segments to alamethicin, ligand binding to the extramembrane segments could alter the structure of the extramembrane domains and influence the association states or association numbers of alamethicin in the membranes. We could then assess those changes by using single-channel current recording. We found that we could modulate channel assembly and eventual ion flux with attached leucine-zipper extramembrane peptide segments. Using conformationally switchable leucine-zipper extramembrane segments that respond to Fe(3+), we fabricated an artificial Fe(3+)-sensitive ion channel; a decrease in the helical content of the extramembrane segment led to an increase in the channel current. When we added a calmodulin C-terminus segment, we formed a channel that was sensitive to Ca(2+). This result demonstrated that we could prepare artificial channels that were sensitive to specific ligands by adding appropriate extramembrane segments from natural protein motifs that respond to external stimuli. In conclusion, our research points to the possibility of creating tailored sensor or signal transduction systems through the conjugation of a conformationally switchable extramembrane peptide/protein segment to a suitable transmembrane peptide segment.

Elucidating the ensemble of functionally-relevant transitions in protein systems with a robotics-inspired method

Background

Many proteins tune their biological function by transitioning between different functional states, effectively acting as dynamic molecular machines. Detailed structural characterization of transition trajectories is central to understanding the relationship between protein dynamics and function. Computational approaches that build on the Molecular Dynamics framework are in principle able to model transition trajectories at great detail but also at considerable computational cost. Methods that delay consideration of dynamics and focus instead on elucidating energetically-credible conformational paths connecting two functionally-relevant structures provide a complementary approach. Effective sampling-based path planning methods originating in robotics have been recently proposed to produce conformational paths. These methods largely model short peptides or address large proteins by simplifying conformational space.

Methods

We propose a robotics-inspired method that connects two given structures of a protein by sampling conformational paths. The method focuses on small- to medium-size proteins, efficiently modeling structural deformations through the use of the molecular fragment replacement technique. In particular, the method grows a tree in conformational space rooted at the start structure, steering the tree to a goal region defined around the goal structure. We investigate various bias schemes over a progress coordinate for balance between coverage of conformational space and progress towards the goal. A geometric projection layer promotes path diversity. A reactive temperature scheme allows sampling of rare paths that cross energy barriers.

Results and conclusions

Experiments are conducted on small- to medium-size proteins of length up to 214 amino acids and with multiple known functionally-relevant states, some of which are more than 13Å apart of each-other. Analysis reveals that the method effectively obtains conformational paths connecting structural states that are significantly different. A detailed analysis on the depth and breadth of the tree suggests that a soft global bias over the progress coordinate enhances sampling and results in higher path diversity. The explicit geometric projection layer that biases the exploration away from over-sampled regions further increases coverage, often improving proximity to the goal by forcing the exploration to find new paths. The reactive temperature scheme is shown effective in increasing path diversity, particularly in difficult structural transitions with known high-energy barriers.

Calcium binding peptide motifs from calmodulin confer divalent ion selectivity to elastin-like polypeptides

Calcium-sensitive elastin-like polypeptides (CELPs) were synthesized by periodically interspersing a calcium-binding peptide sequence from calmodulin within an elastin-like polypeptide (ELP) with the goal of creating thermal and calcium responsive peptide polymers. The CELPs exhibit high sensitivity to calcium compared to monovalent cations but do not exhibit the exquisite selectivity for calcium over other divalent cations, such as magnesium, that is displayed by calmodulin. The CELPs were further used as a building block for the synthesis of calcium-sensitive nanoparticles by fusing a hydrophilic, noncalcium-sensitive ELP block with a CELP block that becomes more hydrophobic upon calcium binding. We show that addition of calcium at concentrations between 50 and 500 mM imparts sufficient amphiphilicity to the diblock polypeptide between 33 and 46 °C to trigger its self-assembly into monodisperse spherical micelles with a hydrodynamic radius of ∼50 nm.

Construction of a Ca(2+)-gated artificial channel by fusing alamethicin with a calmodulin-derived extramembrane segment

Using native chemical ligation, we constructed a Ca(2+)-gated fusion channel protein consisting of alamethicin and the C-terminal domain of calmodulin. At pH 5.4 and in the absence of Ca(2+), this fusion protein yielded a burst-like channel current with no discrete channel conductance levels. However, Ca(2+) significantly lengthened the specific channel open state and increased the mean channel current, while Mg(2+) produced no significant changes in the channel current. On the basis of 8-anilinonaphthalene-1-sulfonic acid (ANS) fluorescent measurement, Ca(2+)-stimulated gating may be related to an increased surface hydrophobicity of the extramembrane segment of the fusion protein.

Calcium-dependent conformational transition of calmodulin determined by Fourier transform infrared spectroscopy

The Ca(2+)-induced conformational changes in calmodulin (CaM) were monitored by Fourier transform infrared spectroscopy (FT-IR) at different molar ratios of Ca(2+) to CaM. The results show that these changes occur in two distinctive transitions. The first transition involves significant changes in the overall secondary structure with a small gain in solvent accessibility, and is completed after the second Ca(2+) binds to both EF-hands of its C-terminal domain. The second transition is accompanied by CaM folding into a tighter, less hydrogen-exchangeable structure, and is completed by the addition of the fourth Ca(2+) to have four Ca(2+) per molecule. Particularly, α-helices in CaM-nCa(2+)(n=0, 1, 2) are less stable than those in CaM-nCa(2+)(n=3, 4).

Use of (113)Cd NMR to probe the native metal binding sites in metalloproteins: an overview

Our laboratories have actively published in this area for several years and the objective of this chapter is to present as comprehensive an overview as possible. Following a brief review of the basic principles associated with (113)Cd NMR methods, we will present the results from a thorough literature search for (113)Cd chemical shifts from metalloproteins. The updated (113)Cd chemical shift figure in this chapter will further illustrate the excellent correlation of the (113)Cd chemical shift with the nature of the coordinating ligands (N, O, S) and coordination number/geometry, reaffirming how this method can be used not only to identify the nature of the protein ligands in uncharacterized cases but also the dynamics at the metal binding site. Specific examples will be drawn from studies on alkaline phosphatase, Ca(2+) binding proteins, and metallothioneins.In the case of Escherichia coli alkaline phosphatase, a dimeric zinc metalloenzyme where a total of six metal ions (three per monomer) are involved directly or indirectly in providing the enzyme with maximal catalytic activity and structural stability, (113)Cd NMR, in conjunction with (13)C and (31)P NMR methods, were instrumental in separating out the function of each class of metal binding sites. Perhaps most importantly, these studies revealed the chemical basis for negative cooperativity that had been reported for this enzyme under metal deficient conditions. Also noteworthy was the fact that these NMR studies preceded the availability of the X-ray crystal structure.In the case of the calcium binding proteins, we will focus on two proteins: calbindin D(9k) and calmodulin. For calbindin D(9k) and its mutants, (113)Cd NMR has been useful both to follow actual changes in the metal binding sites and the cooperativity in the metal binding. Ligand binding to calmodulin has been studied extensively with (113)Cd NMR showing that the metal binding sites are not directly involved in the ligand binding. The (113)Cd chemical shifts are, however, exquisitely sensitive to minute changes in the metal ion environment.In the case of metallothionein, we will reflect upon how (113)Cd substitution and the establishment of specific Cd to Cys residue connectivity by proton-detected heteronuclear (1)H-(113)Cd multiple-quantum coherence methods (HMQC) was essential for the initial establishment of the 3D structure of metallothioneins, a protein family deficient in the regular secondary structural elements of α-helix and β-sheet and the first native protein identified with bound Cd. The (113)Cd NMR studies also enabled the characterization of the affinity of the individual sites for (113)Cd and, in competition experiments, for other divalent metal ions: Zn, Cu, and Hg.