The solution structure of apocalmodulin from Saccharomyces cerevisiae implies a mechanism for its unique Ca2+ binding property

Accurately Predicting Protein pKa Values Using Nonequilibrium Alchemy

The stability, solubility, and function of a protein depend on both its net charge and the protonation states of its individual residues. pKa is a measure of the tendency for a given residue to (de)protonate at a specific pH. Although pKa values can be resolved experimentally, theory and computation provide a compelling alternative. To this end, we assess the applicability of a nonequilibrium (NEQ) alchemical free energy method to the problem of pKa prediction. On a data set of 144 residues that span 13 proteins, we report an average unsigned error of 0.77 ± 0.09, 0.69 ± 0.09, and 0.52 ± 0.04 pK for aspartate, glutamate, and lysine, respectively. This is comparable to current state-of-the-art predictors and the accuracy recently reached using free energy perturbation methods (e.g., FEP+). Moreover, we demonstrate that our open-source, pmx-based approach can accurately resolve the pKa values of coupled residues and observe a substantial performance disparity associated with the lysine partial charges in Amber14SB/Amber99SB*-ILDN, for which an underused fix already exists.

Canonical structural-binding modes in the calmodulin-target protein complexes

Intracellular calcium sensor protein calmodulin (CaM) belongs to the large EF-hand protein superfamily. CaM shows a unique and not fully understood ability to bind to multiple targets, allows them to participate in a variety of regulatory processes. The protein has two approximately symmetrical globular domains (the N- and C-lobes). Analysis of the CaM-binding sites of target proteins showed that they have two hydrophobic 'anchor' amino acids separated by 10 to 17 residues. Consequently, several CaM-binding motifs: {1-10}, {1-11}, {1-13}, {1-14}, {1-16}, {1-17}, differing by the distance between the two anchor residues along the amino acid sequence, have been identified. Despite extensive structural information on the role of target-protein amino acid residues in the formation of complexes with CaM, much less is known about the role of amino acids from CaM contributing to these interactions. In this work, a quantitative analysis of the contact surfaces of CaM and target proteins has been carried out for 35 representative three-dimensional structures. It has been shown that, in addition to the two hydrophobic terminal residues of the target fragment, the interaction also involves residues that are 4 residues earlier in the sequence (binding mode {1-5}). It has also been found that the N- and C-lobes of CaM bind the {1-5} motif located at the ends of the target in a structurally identical manner. Methionine residues at positions 51 (corresponding to 124 in the C-lobe), 71 (144), and 72 (145) of the CaM amino acid sequence are key hydrophobic residues for this interaction. They are located at the N- and C-boundaries of the even EF-hand motifs. The hydrophobic core of CaM ('Ф-quatrefoil') consists of 10 amino acids in the N-lobe (and in the C-lobe): Phe16 (Phe89), Phe19 (Phe92), Ile27 (Ile100), Thr29 (Ala102), Leu32 (Leu105), Ile52 (Ile125), Val55 (Ala128), Ile63 (Val136), Phe65 (Tyr138), and Phe68 (Phe141) and do not intersect with the target-binding methionine residues. CaM belongs to the 'dynamic' group of EF-hand proteins, in which calcium and protein ligand binding causes only global conformational changes but does not alter the conservative 'black' and 'grey' clusters described in our earlier works (PLoS One. 2014; 9(10):e109287). The membership of CaM in the 'dynamic' group is determined by the triggering and protective methionine layer: Met51 (Met124), Met71 (Met144) and Met72 (Met145). HIGHLIGHTSInterchain interactions in the unique 35 CaM complex structures were analyzed.Methionine amino acids of the N- and C-lobes of CaM form triggering and protective layers.Interactions of the target terminal residues with these methionine layers are structurally identical.CaM belonging to the 'dynamic' group is determined by the triggering and protective methionine layer.Communicated by Ramaswamy H. Sarma.

Combined Pulsed Electron Double Resonance EPR and Molecular Dynamics Investigations of Calmodulin Suggest Effects of Crowding Agents on Protein Structures

Calmodulin (CaM) is a highly dynamic Ca²⁺-binding protein that exhibits large conformational changes upon binding Ca²⁺ and target proteins. Although it is accepted that CaM exists in an equilibrium of conformational states in the absence of target protein, the physiological relevance of an elongated helical linker region in the Ca²⁺-replete form has been highly debated. In this study, we use PELDOR (pulsed electron–electron double resonance) EPR measurements of a doubly spin-labeled CaM variant to assess the conformational states of CaM in the apo-, Ca²⁺-bound, and Ca²⁺ plus target peptide-bound states. Our findings are consistent with a three-state conformational model of CaM, showing a semi-open apo-state, a highly extended Ca²⁺-replete state, and a compact target protein-bound state. Molecular dynamics simulations suggest that the presence of glycerol, and potentially other molecular crowding agents, has a profound effect on the relative stability of the different conformational states. Differing experimental conditions may explain the discrepancies in the literature regarding the observed conformational state(s) of CaM, and our PELDOR measurements show good evidence for an extended conformation of Ca²⁺-replete CaM similar to the one observed in early X-ray crystal structures.

Modeling Calcium Signaling in S. cerevisiae Highlights the Role and Regulation of the Calmodulin-Calcineurin Pathway in Response to Hypotonic Shock

Calcium homeostasis and signaling processes in Saccharomyces cerevisiae, as well as in any eukaryotic organism, depend on various transporters and channels located on both the plasma and intracellular membranes. The activity of these proteins is regulated by a number of feedback mechanisms that act through the calmodulin-calcineurin pathway. When exposed to hypotonic shock (HTS), yeast cells respond with an increased cytosolic calcium transient, which seems to be conditioned by the opening of stretch-activated channels. To better understand the role of each channel and transporter involved in the generation and recovery of the calcium transient—and of their feedback regulations—we defined and analyzed a mathematical model of the calcium signaling response to HTS in yeast cells. The model was validated by comparing the simulation outcomes with calcium concentration variations before and during the HTS response, which were observed experimentally in both wild-type and mutant strains. Our results show that calcium normally enters the cell through the High Affinity Calcium influx System and mechanosensitive channels. The increase of the plasma membrane tension, caused by HTS, boosts the opening probability of mechanosensitive channels. This event causes a sudden calcium pulse that is rapidly dissipated by the activity of the vacuolar transporter Pmc1. According to model simulations, the role of another vacuolar transporter, Vcx1, is instead marginal, unless calcineurin is inhibited or removed. Our results also suggest that the mechanosensitive channels are subject to a calcium-dependent feedback inhibition, possibly involving calmodulin. Noteworthy, the model predictions are in accordance with literature results concerning some aspects of calcium homeostasis and signaling that were not specifically addressed within the model itself, suggesting that it actually depicts all the main cellular components and interactions that constitute the HTS calcium pathway, and thus can correctly reproduce the shaping of the calcium signature by calmodulin- and calcineurin-dependent complex regulations. The model predictions also allowed to provide an interpretation of different regulatory schemes involved in calcium handling in both wild-type and mutants yeast strains. The model could be easily extended to represent different calcium signals in other eukaryotic cells.

A Free-Energy Landscape Analysis of Calmodulin Obtained from an NMR Data-Utilized Multi-Scale Divide-and-Conquer Molecular Dynamics Simulation

Calmodulin (CaM) is a multifunctional calcium-binding protein, which regulates a variety of biochemical processes. CaM acts through its conformational changes and complex formation with its target enzymes. CaM consists of two globular domains (N-lobe and C-lobe) linked by an extended linker region. Upon calcium binding, the N-lobe and C-lobe undergo local conformational changes, followed by a major conformational change of the entire CaM to wrap the target enzyme. However, the regulation mechanisms, such as allosteric interactions, which regulate the large structural changes, are still unclear. In order to investigate the series of structural changes, the free-energy landscape of CaM was obtained by multi-scale divide-and-conquer molecular dynamics (MSDC-MD). The resultant free-energy landscape (FEL) shows that the Ca²⁺ bound CaM (holo-CaM) would take an experimentally famous elongated structure, which can be formed in the early stage of structural change, by breaking the inter-domain interactions. The FEL also shows that important interactions complete the structural change from the elongated structure to the ring-like structure. In addition, the FEL might give a guiding principle to predict mutational sites in CaM. In this study, it was demonstrated that the movement process of macroscopic variables on the FEL may be diffusive to some extent, and then, the MSDC-MD is suitable to the parallel computation.

The mechanism of complex formation between calmodulin and voltage gated calcium channels revealed by molecular dynamics

Calmodulin, a ubiquitous eukaryotic calcium sensor responsible for the regulation of many fundamental cellular processes, is a highly flexible protein and exhibits an unusually wide range of conformations. Furthermore, CaM is known to interact with more than 300 cellular targets. Molecular dynamics (MD) simulation trajectories suggest that EF-hand loops show different magnitudes of flexibility. Therefore, the four EF-hand motifs have different affinities for Ca2+ ions, which enables CaM to function on wide range of Ca2+ ion concentrations. EF-hand loops are 2-3 times more flexible in apo CaM whereas least flexible in Ca2+/CaM-IQ motif complexes. We report a unique intermediate conformation of Ca2+/CaM while transitioning from extended to compact form. We also report the complex formation process between Ca2+/CaM and IQ CaM-binding motifs. Our results showed how IQ motif recognise its binding site on the CaM and how CaM transforms from extended to compact form upon binding to IQ motif.

A Change in Domain Cooperativity Drives the Function of Calnuc

With the increasing incidence of neurodegenerative disorders, there is an urgent need to understand the protein folding process. Examining the folding process of multidomain proteins remains a prime challenge, as their complex conformational dynamics make them highly susceptible to misfolding and/or aggregation. The presence of multiple domains in a protein can lead to interaction between the partially folded domains, thereby driving misfolding and/or aggregation. Calnuc is one such multidomain protein for which Ca²⁺ binding plays a pivotal role in governing its structural dynamics and stability and, presumably, in directing its interactions with other proteins. We demonstrate differential structural dynamics between the Ca²⁺-free and Ca²⁺-bound forms of calnuc. In the absence of Ca²⁺, full-length calnuc displays equilibrium structural transitions with four intermediate states, reporting a sum of the behavioral properties of its individual domains. Fragment-based studies illustrate the sequential events of structure adoption proceeding in the following order: EF domain followed by the NT and LZ domains in the apo state. On the other hand, Ca²⁺ binding increases domain cooperativity and enables the protein to fold as a single unit. Single-tryptophan mutant proteins, designed in a domain-dependent manner, confirm an increase in the number of interdomain interactions in the Ca²⁺-bound form as compared to the Ca²⁺-free state of the protein, thereby providing insight into its folding process. The attenuated domain crosstalk in apo-calnuc is likely to influence and regulate its physiologically important intermolecular interactions.

Ca2+-Dependent Switch of Calmodulin Interaction Mode with Tandem IQ Motifs in the Scaffolding Protein IQGAP1

IQ domain GTPase-activating scaffolding protein 1 (IQGAP1) mediates cytoskeleton, cell migration, proliferation, and apoptosis events. Calmodulin (CaM) modulates IQGAP1 functions by binding to its four tandem IQ motifs. Exactly how CaM binds the IQ motifs and which functions of IQGAP1 CaM regulates and how are fundamental mechanistic questions. We combine experimental pull-down assays, mutational data, and molecular dynamics simulations to understand the IQ-CaM complexes with and without Ca2+ at the atomic level. Apo-CaM favors the IQ3 and IQ4 motifs but not the IQ1 and IQ2 motifs that lack two hydrophobic residues for interactions with apo-CaM's hydrophobic pocket. Ca2+-CaM binds all four IQ motifs, with both N- and C-lobes tightly wrapped around each motif. Ca2+ promotes IQ-CaM interactions and increases the amount of IQGAP1-loaded CaM for IQGAP1-mediated signaling. Collectively, we describe IQ-CaM binding in atomistic detail and feature the emergence of Ca2+ as a key modulator of the CaM-IQGAP1 interactions.

Inhibition of bacterial ubiquitin ligases by SidJ-calmodulin catalysed glutamylation

The family of bacterial SidE enzymes catalyzes phosphoribosyl-linked (PR) serine ubiquitination and promotes infectivity of Legionella pneumophilia, a pathogenic bacterium causing Legionnaires’ disease^1,2,3. SidEs share the genetic locus with the Legionella effector SidJ that spatiotemporally opposes their toxicity in yeast and mammalian cells, through an unknown mechanism^4–6. Deletion of SidJ leads to a significant defect in the growth of Legionella in both its natural host amoeba and in murine macrophages^4,5. Here, we demonstrate that SidJ is a glutamylase that modifies the catalytic glutamate in the mono-ADPribosyl transferase (mART) domain of SdeA thus blocking its ubiquitin (Ub) ligase activity. SidJ glutamylation activity requires interaction with Calmodulin (CaM), a eukaryotic specific co-factor, and can be regulated by intracellular changes in Ca²⁺ concentrations. The cryo-EM structure of SidJ/human apo-CaM complex revealed the architecture of this unique heterodimeric glutamylase. In infected cells, we show that SidJ mediates glutamylation of SidEs on the surface of Legionella-containing vacuoles (LCVs). Using quantitative proteomics, we also uncovered multiple host proteins as putative targets of SidJ-mediated glutamylation. Collectively, this study reveals the mechanism of SidE ligases inhibition by a SidJ/CaM glutamylase and opens new avenues for studying protein glutamylation, an understudied protein modification in higher eukaryotes.

Effect of Ca2+ on the promiscuous target-protein binding of calmodulin

Calmodulin (CaM) is a calcium sensing protein that regulates the function of a large number of proteins, thus playing a crucial part in many cell signaling pathways. CaM has the ability to bind more than 300 different target peptides in a Ca2+-dependent manner, mainly through the exposure of hydrophobic residues. How CaM can bind a large number of targets while retaining some selectivity is a fascinating open question. Here, we explore the mechanism of CaM selective promiscuity for selected target proteins. Analyzing enhanced sampling molecular dynamics simulations of Ca2+-bound and Ca2+-free CaM via spectral clustering has allowed us to identify distinct conformational states, characterized by interhelical angles, secondary structure determinants and the solvent exposure of specific residues. We searched for indicators of conformational selection by mapping solvent exposure of residues in these conformational states to contacts in structures of CaM/target peptide complexes. We thereby identified CaM states involved in various binding classes arranged along a depth binding gradient. Binding Ca2+ modifies the accessible hydrophobic surface of the two lobes and allows for deeper binding. Apo CaM indeed shows shallow binding involving predominantly polar and charged residues. Furthermore, binding to the C-terminal lobe of CaM appears selective and involves specific conformational states that can facilitate deep binding to target proteins, while binding to the N-terminal lobe appears to happen through a more flexible mechanism. Thus the long-ranged electrostatic interactions of the charged residues of the N-terminal lobe of CaM may initiate binding, while the short-ranged interactions of hydrophobic residues in the C-terminal lobe of CaM may account for selectivity. This work furthers our understanding of the mechanism of CaM binding and selectivity to different target proteins and paves the way towards a comprehensive model of CaM selectivity.

Inference of Calmodulin's Ca2+-Dependent Free Energy Landscapes via Gaussian Mixture Model Validation

A free energy landscape estimation method based on the well-known Gaussian mixture model (GMM) is used to compare the efficiencies of thermally enhanced sampling methods with respect to regular molecular dynamics. The simulations are carried out on two binding states of calmodulin, and the free energy estimation method is compared with other estimators using a toy model. We show that GMM with cross-validation provides a robust estimate that is not subject to overfitting. The continuous nature of Gaussians provides better estimates on sparse data than canonical histogramming. We find that diffusion properties determine the sampling method effectiveness, such that diffusion-dominated apo calmodulin is most efficiently sampled by regular molecular dynamics, while holo calmodulin, with its rugged free energy landscape, is better sampled by enhanced sampling methods.

Conformational landscapes of membrane proteins delineated by enhanced sampling molecular dynamics simulations

The expansion of computational power, better parameterization of force fields, and the development of novel algorithms to enhance the sampling of the free energy landscapes of proteins have allowed molecular dynamics (MD) simulations to become an indispensable tool to understand the function of biomolecules. The temporal and spatial resolution of MD simulations allows for the study of a vast number of processes of interest. Here, we review the computational efforts to uncover the conformational free energy landscapes of a subset of membrane proteins: ion channels, transporters and G-protein coupled receptors. We focus on the various enhanced sampling techniques used to study these questions, how the conclusions come together to build a coherent picture, and the relationship between simulation outcomes and experimental observables.

Residue-residue interactions regulating the Ca2+-induced EF-hand conformation changes in calmodulin

Calmodulin (CaM) is a Ca2+-binding messenger protein having four Ca2+-binding motifs named 'EF-hand'; the EF-hand motifs undergo a conformation change induced by Ca2+-binding. In order to study how Ca2+-binding induces the conformation change of EF-hand motifs and which residues are involved in the reaction, two 1μ second long MD simulations were independently performed from the apo- and holo-CaM and their structures and interactions were compared. The Ca2+-binding weakens the helix-helix interaction in all EF-hand, however, the holo-CaM MD adopted the close-like form. The correlation coefficients obtained from the two MDs show the residues comprising interactions being involved in their close-open conformation changes; most of these residues are hydrophobic amino acids but some of them are hydrophilic (T34, H107, N111 and Q143). The hydrophilic residues are expected to lock the EF-hands by their side-chains and main-chain carbonyl oxygen of another hydrophobic residue. Furthermore, the interaction pattern of EF-hand3 and 4 are similar to each other. On the other hand, the interaction pattern of EF-hand2 is different from others; its polar residues are expected to play an important role in regulating the EF-hand2 conformation.

Size-dependent impact of CNTs on dynamic properties of calmodulin

There are growing concerns about the biosafety of nanomaterials such as carbon nanotubes (CNTs) as their applications become more widespread. We report here a theoretical and experimental study of the binding of various sizes of CNTs [CNT (4,4), (5,5), (6,6) and (7,7)] to calmodulin (CaM) protein and, in particular, their impact on the Ca(2+)-dependent dynamic properties of CaM. Our simulations show that all the CNTs can plug into the hydrophobic binding pocket of Ca(2+)-bound CaM with binding affinities comparable with the native substrate M13 peptide. Even though CNT (4,4) shows a similar behavior to the M13 peptide in its dissociation from Ca(2+)-free CaM, wider CNTs still bind firmly to CaM, indicating a potential failure of Ca(2+) regulation. Such a size-dependent impact of CNTs on the dynamic properties of CaM is a result of the excessively strong hydrophobic interactions between the wider CNTs and CaM. These simulation results were confirmed by circular dichroism spectroscopy, which showed that the secondary structures of CaM become insensitive to Ca(2+) concentrations after the addition of CNTs. Our findings indicate that the cytotoxicity of nanoparticles to proteins arises not only from the inhibition of static protein structures (binding pockets), but also from impacts on their dynamic properties.

Mapping conformational dynamics of proteins using torsional dynamics simulations

All-atom molecular dynamics simulations are widely used to study the flexibility of protein conformations. However, enhanced sampling techniques are required for simulating protein dynamics that occur on the millisecond timescale. In this work, we show that torsional molecular dynamics simulations enhance protein conformational sampling by performing conformational search in the low-frequency torsional degrees of freedom. In this article, we use our recently developed torsional-dynamics method called Generalized Newton-Euler Inverse Mass Operator (GNEIMO) to study the conformational dynamics of four proteins. We investigate the use of the GNEIMO method in simulations of the conformationally flexible proteins fasciculin and calmodulin, as well as the less flexible crambin and bovine pancreatic trypsin inhibitor. For the latter two proteins, the GNEIMO simulations with an implicit-solvent model reproduced the average protein structural fluctuations and sample conformations similar to those from Cartesian simulations with explicit solvent. The application of GNEIMO with replica exchange to the study of fasciculin conformational dynamics produced sampling of two of this protein's experimentally established conformational substates. Conformational transition of calmodulin from the Ca(2+)-bound to the Ca(2+)-free conformation occurred readily with GNEIMO simulations. Moreover, the GNEIMO method generated an ensemble of conformations that satisfy about half of both short- and long-range interresidue distances obtained from NMR structures of holo to apo transitions in calmodulin. Although unconstrained all-atom Cartesian simulations have failed to sample transitions between the substates of fasciculin and calmodulin, GNEIMO simulations show the transitions in both systems. The relatively short simulation times required to capture these long-timescale conformational dynamics indicate that GNEIMO is a promising molecular-dynamics technique for studying domain motion in proteins.

The many faces of calmodulin in cell proliferation, programmed cell death, autophagy, and cancer

Calmodulin (CaM) is a ubiquitous Ca(2+) receptor protein mediating a large number of signaling processes in all eukaryotic cells. CaM plays a central role in regulating a myriad of cellular functions via interaction with multiple target proteins. This review focuses on the action of CaM and CaM-dependent signaling systems in the control of vertebrate cell proliferation, programmed cell death and autophagy. The significance of CaM and interconnected CaM-regulated systems for the physiology of cancer cells including tumor stem cells, and processes required for tumor progression such as growth, tumor-associated angiogenesis and metastasis are highlighted. Furthermore, the potential targeting of CaM-dependent signaling processes for therapeutic use is discussed.

Surface dynamics in allosteric regulation of protein-protein interactions: modulation of calmodulin functions by Ca2+

Knowledge of the structural basis of protein-protein interactions (PPI) is of fundamental importance for understanding the organization and functioning of biological networks and advancing the design of therapeutics which target PPI. Allosteric modulators play an important role in regulating such interactions by binding at site(s) orthogonal to the complex interface and altering the protein's propensity for complex formation. In this work, we apply an approach recently developed by us for analyzing protein surfaces based on steered molecular dynamics simulation (SMD) to the study of the dynamic properties of functionally distinct conformations of a model protein, calmodulin (CaM), whose ability to interact with target proteins is regulated by the presence of the allosteric modulator Ca(2+). Calmodulin is a regulatory protein that acts as an intracellular Ca(2+) sensor to control a wide variety of cellular processes. We demonstrate that SMD analysis is capable of pinpointing CaM surfaces implicated in the recognition of both the allosteric modulator Ca(2+) and target proteins. Our analysis of changes in the dynamic properties of the CaM backbone elicited by Ca(2+) binding yielded new insights into the molecular mechanism of allosteric regulation of CaM-target interactions.

The Neuronal Calcium Sensor protein Acrocalcin: a potential target of calmodulin regulation during development in the coral Acropora millepora

To understand the calcium-mediated signalling pathways underlying settlement and metamorphosis in the Scleractinian coral Acropora millepora, a predicted protein set derived from larval cDNAs was scanned for the presence of EF-hand domains (Pfam Id: PF00036). This approach led to the identification of a canonical calmodulin (AmCaM) protein and an uncharacterised member of the Neuronal Calcium Sensor (NCS) family of proteins known here as Acrocalcin (AmAC). While AmCaM transcripts were present throughout development, AmAC transcripts were not detected prior to gastrulation, after which relatively constant mRNA levels were detected until metamorphosis and settlement. The AmAC protein contains an internal CaM-binding site and was shown to interact in vitro with AmCaM. These results are consistent with the idea that AmAC is a target of AmCaM in vivo, suggesting that this interaction may regulate calcium-dependent processes during the development of Acropora millepora.

Specific conformation and Ca(2+)-binding mode of yeast calmodulin: insight into evolutionary development

The vertebrate calmodulin is configured with two structurally independent globular lobes in N- and C-terminus, and a flexible central linker. Distinctly, two lobes of calmodulin from Saccharomyces cerevisiae (yCaM) interact and influence the Ca(2+)-binding profile of each other. We explored this further using the mutant proteins with eliminated Ca(2+)-binding ability in one of the lobes and found that the Ca(2+)-bound N-lobe associates with the Ca(2+)-free C-lobe to gain the Ca(2+) affinity of a wild-type level. Next, analysing series of C-terminal residue truncation mutant, we found that the truncation of C-terminal three residues induce the hyper Ca(2+) affinity. These residues are also important for the general structural behaviour of calmodulin, such as Ca(2+)-induced slow mobility shift in polyacrylamide gel electrophoresis and for the ability to activate Cmk1p (yeast calmodulin kinase). These suggest: (i) when Ca(2+) occupies only N-lobe, two lobes interact and form the stable intermediate leading to a proper level of Ca(2+) affinity; (ii) the C-terminal three residues are required to prohibit abnormal stabilization of the intermediate promoting abnormally high Ca(2+) affinity and for recognition of target enzymes. A model for Ca(2+) and target bindings of yCaM is proposed. Evolutional aspect concerning the biological significance of this model was discussed.

Similarities between protein folding and granular jamming

Grains and glasses, widely different materials, arrest their motions upon decreasing temperature and external load, respectively, in common ways, leading to a universal jamming phase diagram conjecture. However, unified theories are lacking, mainly because of the disparate nature of the particle interactions. Here we demonstrate that folded proteins exhibit signatures common to both glassiness and jamming by using temperature- and force-unfolding molecular dynamics simulations. Upon folding, proteins develop a peak in the interatomic force distributions that falls on a universal curve with experimentally measured forces on jammed grains and droplets. Dynamical signatures are found as a dramatic slowdown of stress relaxation upon folding. Together with granular similarities, folding is tied not just to the jamming transition, but a more nuanced picture of anisotropy, preparation protocol and internal interactions emerges. Results have implications for designing stable polymers and can open avenues to link protein folding to jamming theory.

The dynamics of Ca2+ ions within the solvation shell of calbindin D9k

The encounter of a Ca²⁺ ion with a protein and its subsequent binding to specific binding sites is an intricate process that cannot be fully elucidated from experimental observations. We have applied Molecular Dynamics to study this process with atomistic details, using Calbindin D9k (CaB) as a model protein. The simulations show that in most of the time the Ca²⁺ ion spends within the Debye radius of CaB, it is being detained at the 1st and 2nd solvation shells. While being detained near the protein, the diffusion coefficient of the ion is significantly reduced. However, due to the relatively long period of detainment, the ion can scan an appreciable surface of the protein. The enhanced propagation of the ion on the surface has a functional role: significantly increasing the ability of the ion to scan the protein's surface before being dispersed to the bulk. The contribution of this mechanism to Ca²⁺ binding becomes significant at low ion concentrations, where the intervals between successive encounters with the protein are getting longer. The efficiency of the surface diffusion is affected by the distribution of charges on the protein's surface. Comparison of the Ca²⁺ binding dynamics in CaB and its E60D mutant reveals that in the wild type (WT) protein the carboxylate of E60 function as a preferred landing-site for the Ca²⁺ arriving from the bulk, followed by delivering it to the final binding site. Replacement of the glutamate by aspartate significantly reduced the ability to transfer Ca²⁺ ions from D60 to the final binding site, explaining the observed decrement in the affinity of the mutated protein to Ca²⁺.

Modulation of calmodulin plasticity by the effect of macromolecular crowding

In vitro biochemical reactions are most often studied in dilute solution, a poor mimic of the intracellular space of eukaryotic cells, which are crowded with mobile and immobile macromolecules. Such crowded conditions exert volume exclusion and other entropic forces that have the potential to impact chemical equilibria and reaction rates. In this article, we used the well-characterized and ubiquitous molecule calmodulin (CaM) and a combination of theoretical and experimental approaches to address how crowding impacts CaM's conformational plasticity. CaM is a dumbbell-shaped molecule that contains four EF hands (two in the N-lobe and two in the C-lobe) that each could bind Ca(2+), leading to stabilization of certain substates that favor interactions with other target proteins. Using coarse-grained molecular simulations, we explored the distribution of CaM conformations in the presence of crowding agents. These predictions, in which crowding effects enhance the population of compact structures, were then confirmed in experimental measurements using fluorescence resonance energy transfer techniques of donor- and acceptor-labeled CaM under normal and crowded conditions. Using protein reconstruction methods, we further explored the folding-energy landscape and examined the structural characteristics of CaM at free-energy basins. We discovered that crowding stabilizes several different compact conformations, which reflects the inherent plasticity in CaM's structure. From these results, we suggest that the EF hands in the C-lobe are flexible and can be thought of as a switch, while those in the N-lobe are stiff, analogous to a rheostat. New combinatorial signaling properties may arise from the product of the differential plasticity of the two distinct lobes of CaM in the presence of crowding. We discuss the implications of these results for modulating CaM's ability to bind Ca(2+) and target proteins.

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.

Structural basis for the interaction of the myosin light chain Mlc1p with the myosin V Myo2p IQ motifs

Calmodulin, regulatory, and essential myosin light chain are evolutionary conserved proteins that, by binding to IQ motifs of target proteins, regulate essential intracellular processes among which are efficiency of secretory vesicles release at synapsis, intracellular signaling, and regulation of cell division. The yeast Saccharomyces cerevisiae calmodulin Cmd1 and the essential myosin light chain Mlc1p share the ability to interact with the class V myosin Myo2p and Myo4 and the class II myosin Myo1p. These myosins are required for vesicle, organelle, and mRNA transport, spindle orientation, and cytokinesis. We have used the budding yeast model system to study how calmodulin and essential myosin light chain selectively regulate class V myosin function. NMR structural analysis of uncomplexed Mlc1p and interaction studies with the first three IQ motifs of Myo2p show that the structural similarities between Mlc1p and the other members of the EF-hand superfamily of calmodulin-like proteins are mainly restricted to the C-lobe of these proteins. The N-lobe of Mlc1p presents a significantly compact and stable structure that is maintained both in the free and complexed states. The Mlc1p N-lobe interacts with the IQ motif in a manner that is regulated both by the IQ motifs sequence as well as by light chain structural features. These characteristic allows a distinctive interaction of Mlc1p with the first IQ motif of Myo2p when compared with calmodulin. This finding gives us a novel view of how calmodulin and essential light chain, through a differential binding to IQ1 of class V myosin motor, regulate this activity during vegetative growth and cytokinesis.

Prediction of volatile anesthetic binding sites in proteins

Computational methods designed to predict and visualize ligand protein binding interactions were used to characterize volatile anesthetic (VA) binding sites and unoccupied pockets within the known structures of VAs bound to serum albumin, luciferase, and apoferritin. We found that both the number of protein atoms and methyl hydrogen, which are within approximately 8 A of a potential ligand binding site, are significantly greater in protein pockets where VAs bind. This computational approach was applied to structures of calmodulin (CaM), which have not been determined in complex with a VA. It predicted that VAs bind to [Ca(2+)](4)-CaM, but not to apo-CaM, which we confirmed with isothermal titration calorimetry. The VA binding sites predicted for the structures of [Ca(2+)](4)-CaM are located in hydrophobic pockets that form when the Ca(2+) binding sites in CaM are saturated. The binding of VAs to these hydrophobic pockets is supported by evidence that halothane predominantly makes contact with aliphatic resonances in [Ca(2+)](4)-CaM (nuclear Overhauser effect) and increases the Ca(2+) affinity of CaM (fluorescence spectroscopy). Our computational analysis and experiments indicate that binding of VA to proteins is consistent with the hydrophobic effect and the Meyer-Overton rule.

A molecular dynamics study of the effect of Ca2+ removal on calmodulin structure

Calmodulin is a small (148 residues), ubiquitous, highly-conserved Ca(2+) binding protein serving as a modulator of many calcium-dependent processes. In this study, we followed, by means of molecular dynamics, the structural stability of the protein when one of its four bound Ca(2+) ions is removed, and compared it to a simulation of the fully Ca(2+) bound protein. We found that the removal of a single Ca(2+) ion from the N-lobe of the protein, which has a lower affinity for the ion, is sufficient to initiate a considerable structural rearrangement. Although the overall structure of the fully 4 Ca(2+) bound protein remained intact in the extended conformation, the Ca(2+)-removed protein changed its conformation into a compact state. The observation that the 3 Ca(2+) loaded protein assumes a compacted solution state is in accord with experimental observation that the NSCP protein, which binds only three Ca(2+) ions, is natively in a compact state. Examination of the folding dynamics reveals a cooperation between the C-lobe, N-lobe, and the interdomain helix that enable the conformation change. The forces driving this conformational change are discussed.

Al3+ interaction sites of calmodulin and the Al3+ effect on target binding of calmodulin

The interaction between calmodulin (CaM) and Al(3+) was studied by spectroscopic methods. Heteronuclear two-dimensional NMR data indicated that peaks related to the both lobes and middle of the central helix of CaM are largely affected by Al(3+). But chemical shift perturbation suggested that overall conformation of Ca(2+)-loaded CaM is not changed by Al(3+) binding. It is thought that Al(3+) interaction to the middle of the central helix is a key for the property of CaM's target recognition. If the structure and/or flexibility of the central helix are/is changed by Al(3+), target affinity to CaM must be influenced by Al(3+). Thus, we performed surface plasmon resonance experiments to observe the effect of Al(3+) on the target recognition by CaM. The data clearly indicated that target affinity to CaM is reduced by addition of Al(3+). All the results presented here support a hypothesis that Al(3+) may affect on the Ca(2+) signaling pathway in cells.

Sites on Calmodulin That Interact with the C-terminal Tail of Cav1.2 Channel

Two fragments of the C-terminal tail of the alpha(1) subunit (CT1, amino acids 1538-1692 and CT2, amino acids 1596-1692) of human cardiac L-type calcium channel (Ca(V)1.2) have been expressed, refolded, and purified. A single Ca(2+)-calmodulin binds to each fragment, and this interaction with Ca(2+)-calmodulin is required for proper folding of the fragment. Ca(2+)-calmodulin, bound to these fragments, is in a more extended conformation than calmodulin bound to a synthetic peptide representing the IQ motif, suggesting that either the conformation of the IQ sequence is different in the context of the longer fragment, or other sequences within CT2 contribute to the binding of calmodulin. NMR amide chemical shift perturbation mapping shows the backbone conformation of calmodulin is nearly identical when bound to CT1 and CT2, suggesting that amino acids 1538-1595 do not contribute to or alter calmodulin binding to amino acids 1596-1692 of Ca(V)1.2. The interaction with CT2 produces the greatest changes in the backbone amides of hydrophobic residues in the N-lobe and hydrophilic residues in the C-lobe of calmodulin and has a greater effect on residues located in Ca(2+) binding loops I and II in the N-lobe relative to loops III and IV in the C-lobe. In conclusion, Ca(2+)-calmodulin assumes a novel conformation when part of a complex with the C-terminal tail of the Ca(V)1.2 alpha(1) subunit that is not duplicated by synthetic peptides corresponding to the putative binding motifs.

Calmodulin's flexibility allows for promiscuity in its interactions with target proteins and peptides

The small bilobal calcium regulatory protein calmodulin (CaM) activates numerous target enzymes in response to transient changes in intracellular calcium concentrations. Binding of calcium to the two helix-loop-helix calcium-binding motifs in each of the globular domains induces conformational changes that expose a methionine-rich hydrophobic patch on the surface of each domain of the protein, which it uses to bind to peptide sequences in its target enzymes. Although these CaM-binding domains typically have little sequence identity, the positions of several bulky hydrophobic residues are often conserved, allowing for classification of CaM-binding domains into recognition motifs, such as the 1-14 and 1-10 motifs. For calcium-independent binding of CaM, a third motif known as the IQ motif is also common. Many CaM-peptide complexes have globular conformations, where CaM's central linker connecting the two domains unwinds, allowing the protein to wrap around a single predominantly alpha-helical target peptide sequence. However, novel structures have recently been reported where the conformation of CaM is highly dissimilar to these globular complexes, in some instances with less than a full compliment of bound calcium ions, as well as novel stoichiometries. Furthermore, many divergent CaM isoforms from yeast and plant species have been discovered with unique calcium-binding and enzymatic activation characteristics compared to the single CaM isoform found in mammals.