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

Ligand-induced protein transition state stabilization switches the binding pathway from conformational selection to induced fit

Protein-ligand complex formation is fundamental to biological function. A central question is whether proteins spontaneously adopt binding-competent conformations to which ligands bind conformational selection (CS) or whether ligands induce the binding-competent conformation induced fit (IF). Here, we resolve the CS and IF binding pathways by characterizing protein conformational dynamics over a wide range of ligand concentrations using NMR relaxation dispersion. We determined the relative flux through the two pathways using a four-state binding model that includes both CS and IF. Experiments conducted without ligand show that galectin-3 exchanges between the ground-state conformation and a high-energy conformation similar to the ligand-bound conformation, demonstrating that CS is a plausible pathway. Near-identical crystal structures of the apo and ligand-bound states suggest that the high-energy conformation in solution corresponds to the apo crystal structure. Stepwise additions of the ligand lactose induce progressive changes in the relaxation dispersions that we fit collectively to the four-state model, yielding all microscopic rate constants and binding affinities. The ligand affinity is higher for the bound-like conformation than for the ground state, as expected for CS. Nonetheless, the IF pathway contributes greater than 70% of the total flux even at low ligand concentrations. The higher flux through the IF pathway is explained by considerably higher rates of exchange between the two protein conformations in the ligand-associated state. Thus, the ligand acts to decrease the activation barrier between protein conformations in a manner reciprocal to enzymatic transition-state stabilization of reactions involving ligand transformation.

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

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

Time-resolved solid state NMR of biomolecular processes with millisecond time resolution

We review recent efforts to develop and apply an experimental approach to the structural characterization of transient intermediate states in biomolecular processes that involve large changes in molecular conformation or assembly state. This approach depends on solid state nuclear magnetic resonance (ssNMR) measurements that are performed at very low temperatures, typically 25-30 K, with signal enhancements from dynamic nuclear polarization (DNP). This approach also involves novel technology for initiating the process of interest, either by rapid mixing of two solutions or by a rapid inverse temperature jump, and for rapid freezing to trap intermediate states. Initiation by rapid mixing or an inverse temperature jump can be accomplished in approximately-one millisecond. Freezing can be accomplished in approximately 100 microseconds. Thus, millisecond time resolution can be achieved. Recent applications to the process by which the biologically essential calcium sensor protein calmodulin forms a complex with one of its target proteins and the process by which the bee venom peptide melittin converts from an unstructured monomeric state to a helical, tetrameric state after a rapid change in pH or temperature are described briefly. Future applications of millisecond time-resolved ssNMR are also discussed briefly.

Role of water-bridged interactions in metal ion coupled protein allostery

Allosteric communication between distant parts of proteins controls many cellular functions, in which metal ions are widely utilized as effectors to trigger the allosteric cascade. Due to the involvement of strong coordination interactions, the energy landscape dictating the metal ion binding is intrinsically rugged. How metal ions achieve fast binding by overcoming the landscape ruggedness and thereby efficiently mediate protein allostery is elusive. By performing molecular dynamics simulations for the Ca²⁺ binding mediated allostery of the calmodulin (CaM) domains, each containing two Ca²⁺ binding helix-loop-helix motifs (EF-hands), we revealed the key role of water-bridged interactions in Ca²⁺ binding and protein allostery. The bridging water molecules between Ca²⁺ and binding residue reduces the ruggedness of ligand exchange landscape by acting as a lubricant, facilitating the Ca²⁺ coupled protein allostery. Calcium-induced rotation of the helices in the EF-hands, with the hydrophobic core serving as the pivot, leads to exposure of hydrophobic sites for target binding. Intriguingly, despite being structurally similar, the response of the two symmetrically arranged EF-hands upon Ca²⁺ binding is asymmetric. Breakage of symmetry is needed for efficient allosteric communication between the EF-hands. The key roles that water molecules play in driving allosteric transitions are likely to be general in other metal ion mediated protein allostery.

Natural proteins often utilize allostery in executing a variety of functions. Metal ions are typical cofactors to trigger the allosteric cascade. In this work, using the Ca²⁺ sensor protein calmodulin as the model system, we revealed crucial roles of water-bridged interactions in the metal ion coupled protein allostery. The coordination of the Ca²⁺ to the binding site involves an intermediate in which the water molecule bridges the Ca²⁺ and the liganding residue. The bridging water reduces the free energy barrier height of ligand exchange, therefore facilitating the ligand exchange and allosteric coupling by acting as a lubricant. We also showed that the response of the two symmetrically arranged EF-hand motifs of CaM domains upon Ca²⁺ binding is asymmetric, which is directly attributed to the differing dehydration process of the Ca²⁺ ions and is needed for efficient allosteric communication.

1H R1ρ relaxation dispersion experiments in aromatic side chains

Aromatic side chains are attractive probes of protein dynamic, since they are often key residues in enzyme active sites and protein binding sites. Dynamic processes on microsecond to millisecond timescales can be studied by relaxation dispersion experiments that attenuate conformational exchange contributions to the transverse relaxation rate by varying the refocusing frequency of applied radio-frequency fields implemented as either CPMG pulse trains or continuous spin-lock periods. Here we present an aromatic ¹H R_1ρ relaxation dispersion experiment enabling studies of two to three times faster exchange processes than achievable by existing experiments for aromatic side chains. We show that site-specific isotope labeling schemes generating isolated ¹H-¹³C spin pairs with vicinal ²H-¹²C moieties are necessary to avoid anomalous relaxation dispersion profiles caused by Hartmann-Hahn matching due to the ³J_HH couplings and limited chemical shift differences among ¹H spins in phenylalanine, tyrosine and the six-ring moiety of tryptophan. This labeling pattern is sufficient in that remote protons do not cause additional complications. We validated the approach by measuring ring-flip kinetics in the small protein GB1. The determined rate constants, k_flip, agree well with previous results from ¹³C R_1ρ relaxation dispersion experiments, and yield ¹H chemical shift differences between the two sides of the ring in good agreement with values measured under slow-exchange conditions. The aromatic¹H R_1ρ relaxation dispersion experiment in combination with the site-selective ¹H-¹³C/²H-¹²C labeling scheme enable measurement of exchange rates up to k_ex = 2k_flip = 80,000 s^-1, and serve as a useful complement to previously developed ¹³C-based methods.

Coarse-Grained Modeling and Molecular Dynamics Simulations of Ca2+-Calmodulin

Calmodulin (CaM) is a calcium-binding protein that transduces signals to downstream proteins through target binding upon calcium binding in a time-dependent manner. Understanding the target binding process that tunes CaM’s affinity for the calcium ions (Ca²⁺), or vice versa, may provide insight into how Ca²⁺-CaM selects its target binding proteins. However, modeling of Ca²⁺-CaM in molecular simulations is challenging because of the gross structural changes in its central linker regions while the two lobes are relatively rigid due to tight binding of the Ca²⁺ to the calcium-binding loops where the loop forms a pentagonal bipyramidal coordination geometry with Ca²⁺. This feature that underlies the reciprocal relation between Ca²⁺ binding and target binding of CaM, however, has yet to be considered in the structural modeling. Here, we presented a coarse-grained model based on the Associative memory, Water mediated, Structure, and Energy Model (AWSEM) protein force field, to investigate the salient features of CaM. Particularly, we optimized the force field of CaM and that of Ca²⁺ ions by using its coordination chemistry in the calcium-binding loops to match with experimental observations. We presented a “community model” of CaM that is capable of sampling various conformations of CaM, incorporating various calcium-binding states, and carrying the memory of binding with various targets, which sets the foundation of the reciprocal relation of target binding and Ca²⁺ binding in future studies.

Millisecond Time-Resolved Solid-State NMR Reveals a Two-Stage Molecular Mechanism for Formation of Complexes between Calmodulin and a Target Peptide from Myosin Light Chain Kinase

Calmodulin (CaM) mediates a wide range of biological responses to changes in intracellular Ca²⁺ concentrations through its calcium-dependent binding affinities to numerous target proteins. Binding of two Ca²⁺ ions to each of the two four-helix-bundle domains of CaM results in major conformational changes that create a potential binding site for the CaM binding domain of a target protein, which also undergoes major conformational changes to form the complex with CaM. Details of the molecular mechanism of complex formation are not well established, despite numerous structural, spectroscopic, thermodynamic, and kinetic studies. Here, we report a study of the process by which the 26-residue peptide M13, which represents the CaM binding domain of skeletal muscle myosin light chain kinase, forms a complex with CaM in the presence of excess Ca²⁺ on the millisecond time scale. Our experiments use a combination of selective ¹³C labeling of CaM and M13, rapid mixing of CaM solutions with M13/Ca²⁺ solutions, rapid freeze-quenching of the mixed solutions, and low-temperature solid state nuclear magnetic resonance (ssNMR) enhanced by dynamic nuclear polarization. From measurements of the dependence of 2D ¹³C-¹³C ssNMR spectra on the time between mixing and freezing, we find that the N-terminal portion of M13 converts from a conformationally disordered state to an α-helix and develops contacts with the C-terminal domain of CaM in about 2 ms. The C-terminal portion of M13 becomes α-helical and develops contacts with the N-terminal domain of CaM more slowly, in about 8 ms. The level of structural order in the CaM/M13/Ca²⁺ complexes, indicated by ¹³C ssNMR line widths, continues to increase beyond 27 ms.

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.

Resolved Structural States of Calmodulin in Regulation of Skeletal Muscle Calcium Release

Calmodulin (CaM) is proposed to modulate activity of the skeletal muscle sarcoplasmic reticulum (SR) calcium release channel (ryanodine receptor, RyR1 isoform) via a mechanism dependent on the conformation of RyR1-bound CaM. However, the correlation between CaM structure and functional regulation of RyR in physiologically relevant conditions is largely unknown. Here, we have used time-resolved fluorescence resonance energy transfer (TR-FRET) to study structural changes in CaM that may play a role in the regulation of RyR1. We covalently labeled each lobe of CaM (N and C) with fluorescent probes and used intramolecular TR-FRET to assess interlobe distances when CaM is bound to RyR1 in SR membranes, purified RyR1, or a peptide corresponding to the CaM-binding domain of RyR (RyRp). TR-FRET resolved an equilibrium between two distinct structural states (conformations) of CaM, each characterized by an interlobe distance and Gaussian distribution width (disorder). In isolated CaM, at low Ca²⁺, the two conformations of CaM are resolved, centered at 5 nm (closed) and 7 nm (open). At high Ca²⁺, the equilibrium shifts to favor the open conformation. In the presence of RyRp at high Ca²⁺, the closed conformation shifts to a more compact conformation and is the major component. When CaM is bound to full-length RyR1, either purified or in SR membranes, strikingly different results were obtained: 1) the two conformations are resolved and more ordered, 2) the open state is the major component, and 3) Ca²⁺ stabilized the closed conformation by a factor of two. We conclude that the Ca²⁺-dependent structural distribution of CaM bound to RyR1 is distinct from that of CaM bound to RyRp. We propose that the function of RyR1 is tuned to the Ca²⁺-dependent structural dynamics of bound CaM.

Site-selective 1H/2H labeling enables artifact-free 1H CPMG relaxation dispersion experiments in aromatic side chains

Aromatic side chains are often key residues in enzyme active sites and protein binding sites, making them attractive probes of protein dynamics on the millisecond timescale. Such dynamic processes can be studied by aromatic ¹³C or ¹H CPMG relaxation dispersion experiments. Aromatic ¹H CPMG relaxation dispersion experiments in phenylalanine, tyrosine and the six-ring moiety of tryptophan, however, are affected by ³J ¹H-¹H couplings which are causing anomalous relaxation dispersion profiles. Here we show that this problem can be addressed by site-selective ¹H/²H labeling of the aromatic side chains and that artifact-free relaxation dispersion profiles can be acquired. The method has been further validated by measuring folding-unfolding kinetics of the small protein GB1. The determined rate constants and populations agree well with previous results from ¹³C CPMG relaxation dispersion experiments. Furthermore, the CPMG-derived chemical shift differences between the folded and unfolded states are in excellent agreement with those obtained directly from the spectra. In summary, site-selective ¹H/²H labeling enables artifact-free aromatic ¹H CPMG relaxation dispersion experiments in phenylalanine and the six-ring moiety of tryptophan, thereby extending the available methods for studying millisecond dynamics in aromatic protein side chains.

Cation and peptide binding properties of CML7, a calmodulin-like protein from Arabidopsis thaliana

Plants contain a large family of so-called calmodulin-like proteins (CMLs) which differ from canonical calmodulin in that they show greater variability in sequence, length, and number of EF-hand domains. The presence of this extended CML family has raised questions regarding the role of these proteins: are they functionally redundant or do they play specific functions in physiological plant processes? To answer these questions, comprehensive biochemical and structural information on CML proteins is fundamental. Among the 50 CMLs from Arabidopsis thaliana, herein we described the ability of CML7 to bind metal ions focusing on the Ca²⁺ and Mg²⁺ sensing properties, as well as on metal-induced conformational changes. Circular dichroism and nuclear magnetic resonance (NMR) studies indicated that both Ca²⁺ and Mg²⁺ stabilize CML7, as reflected in conformational rearrangements in secondary and tertiary structure and in increases in thermal stability of the protein. However, the conformational changes that binding induces differ between the two metal ions, and only Ca²⁺ binding controls a structural transition that leads to hydrophobic exposure, as suggested by 8-anilino-1-naphthalenesulfonic acid fluorescence. Isothermal titration calorimetry data coupled with NMR experiments revealed the presence of two high affinity Ca²⁺-binding sites in the C-lobe of CML7 and two weaker sites in the N-lobe. The paired nature of these CML7 EF-hands enables them to bind Ca²⁺ with positive cooperativity within each globular domain. Our results clearly place CML7 in the category of Ca²⁺ sensors. Along with this, the protein can bind to a model target peptide (melittin) in a Ca²⁺-dependent manner.

Consequences of Hydrophobic Nanotube Binding on the Functional Dynamics of Signaling Protein Calmodulin

The wide applications of nanomaterials in industry and our daily life have raised growing concerns on their toxicity to human body. Increasing evidence links the cytotoxicity of nanoparticles to the disruption of cellular signaling pathways. Here, we report a computational study on the mechanisms of the cytotoxicity of carbon nanotubes (CNTs) by investigating the direct impacts of CNTs on the functional motions of calmodulin (CaM), which is one of the most important signaling proteins in a cell, and its signaling function relies on the Ca2+ binding-coupled conformational switching. Computational simulations with a coarse-grained model showed that binding of CNTs modifies the conformational equilibrium of CaM and induces the closed-to-open conformational transition, leading to the loss of its Ca2+-sensing ability. In addition, the binding of CNTs drastically increases the calcium affinity of CaM, which may disrupt the Ca2+ homeostasis in a cell. These results suggest that the binding of hydrophobic nanotubes not only inhibits the signaling function of CaM as a calcium sensor but also renders CaM to toxic species through sequestering Ca2+ from other competing calcium-binding proteins, suggesting a new physical mechanism of the cytotoxicity of nanoparticles.

Conformational exchange of aromatic side chains by 1H CPMG relaxation dispersion

Aromatic side chains are attractive probes of protein dynamics on the millisecond time scale, because they are often key residues in enzyme active sites and protein binding sites. Further they allow to study specific processes, like histidine tautomerization and ring flips. Till now such processes have been studied by aromatic ¹³C CPMG relaxation dispersion experiments. Here we investigate the possibility of aromatic ¹H CPMG relaxation dispersion experiments as a complementary method. Artifact-free dispersions are possible on uniformly ¹H and ¹³C labeled samples for histidine δ2 and ε1, as well as for tryptophan δ1. The method has been validated by measuring fast folding-unfolding kinetics of the small protein CspB under native conditions. The determined rate constants and populations agree well with previous results from ¹³C CPMG relaxation dispersion experiments. The CPMG-derived chemical shift differences between the folded and unfolded states are in good agreement with those obtained directly from the spectra. In contrast, the ¹H relaxation dispersion profiles in phenylalanine, tyrosine and the six-ring moiety of tryptophan, display anomalous behavior caused by ³J ¹H-¹H couplings and, if present, strong ¹³C-¹³C couplings. Therefore they require site-selective ¹H/²H and, in case of strong couplings, ¹³C/¹²C labeling. In summary, aromatic ¹H CPMG relaxation dispersion experiments work on certain positions (His δ2, His ε1 and Trp δ1) in uniformly labeled samples, while other positions require site-selective isotope labeling.

Binding and backbone dynamics of protein under topological constraint: calmodulin as a model system

Herein we present the effect of artificially imposed topological constraint on calmodulin (CaM) backbone dynamics and its molecular recognition behavior. While backbone dynamics of CaM remain largely unperturbed, the thermodynamic profile of CaM binding to the smooth-muscle myosin light-chain kinase (smMLCK) peptide is modulated significantly.

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.

Towards Understanding Plant Calcium Signaling through Calmodulin-Like Proteins: A Biochemical and Structural Perspective

Ca2+ ions play a key role in a wide variety of environmental responses and developmental processes in plants, and several protein families with Ca2+-binding domains have evolved to meet these needs, including calmodulin (CaM) and calmodulin-like proteins (CMLs). These proteins have no catalytic activity, but rather act as sensor relays that regulate downstream targets. While CaM is well-studied, CMLs remain poorly characterized at both the structural and functional levels, even if they are the largest class of Ca2+ sensors in plants. The major structural theme in CMLs consists of EF-hands, and variations in these domains are predicted to significantly contribute to the functional versatility of CMLs. Herein, we focus on recent advances in understanding the features of CMLs from biochemical and structural points of view. The analysis of the metal binding and structural properties of CMLs can provide valuable insight into how such a vast array of CML proteins can coexist, with no apparent functional redundancy, and how these proteins contribute to cellular signaling while maintaining properties that are distinct from CaM and other Ca2+ sensors. An overview of the principal techniques used to study the biochemical properties of these interesting Ca2+ sensors is also presented.

Conformational landscape mapping the difference between N-lobes and C-lobes of calmodulin

Calmodulin is a calcium binding protein that consists of four EF-hand domains. The two EF-lobes of calmodulin, called the N-lobe and the C-lobe, arose from duplication and fusion of a precursor EF-hand. The amino acid sequences and the structures of the N-lobe and of the C-lobe are quite similar to each other. The N-lobe and the C-lobe, however, have subtle differences in structure and function. We analyzed the helix positions of calmodulin lobes by the alignment with the pseudo-two fold axis of the EF-lobe. We made a map of conformational landscape of helix positions. The four states of the EF-lobe appeared on two lines in the landscape; these two lines show the trajectory of opening and closing of the EF-lobe. For the N-lobe of calmodulin, the calcium bound form and the apo-forms are on the lower line. The two apo-forms of the C-lobe of calmodulin, with target and without target, are on the upper line. The calcium bound form of the C-lobe is on the lower line. The rearrangement of helix interaction between two the EF-hands is necessary for calcium binding in the C-lobe. The hydrophobic packing in the apo-form of the N-lobe is similar to the packing of the N- and C-lobes of the calcium bound form. However, the packing of C-lobe side chains in the apo-form is different from these other three structures. Our detailed analysis should serve as an example that can be applied to other proteins that undergo changes in conformation upon binding effectors.

On the Ca2+ binding and conformational change in EF-hand domains: Experimental evidence of Ca2+-saturated intermediates of N-domain of calmodulin

Double mutation of Q41L and K75I in the N-domain of calmodulin (N-Cam) stabilizes the closed form of N-Cam such that binding of Ca²⁺ in solution no longer triggers a conformational change to the open form, and its Ca²⁺ binding affinity decreases dramatically. To further investigate the solvation effects on the structure, Ca²⁺ binding affinity and conformational dynamics of this N-Cam double mutant in the Ca²⁺ saturated state, we solved its X-ray structure. Surprisingly, the structure revealed an open conformation of the domain which contradicts its closed conformation in solution. Here we provide evidence that crystallization conditions were responsible for this Ca²⁺-saturated domain open conformation in the crystal. Importantly, we demonstrate that the presence of the crystallization co-precipitant and alcohols were able to induce a progressive opening of the closed form of this domain, in Ca²⁺ saturated state, in solution. However, in the Ca²⁺ depleted state, addition of alcohols was unable to induce any opening of this domain in solution. In addition, in the Ca²⁺ saturated state, the molecular dynamics simulations show that while N-Cam can populate the open and closed conformation, the N-Cam double mutant exclusively populates the closed conformation. Our results provide experimental evidence of intermediate conformations of Ca²⁺-saturated N-Cam in solution. We propose that conformational change of Ca²⁺ sensor EF-hand domains depends on solvation energetics, Ca²⁺ binding to promote the full open form, Ca²⁺ depleted state conformational dynamics, and the chemical properties of the molecules nearby key residues such as those at positions 41 and 75 in N-Cam.

Site-selective 13C labeling of proteins using erythrose

NMR-spectroscopy enables unique experimental studies on protein dynamics at atomic resolution. In order to obtain a full atom view on protein dynamics, and to study specific local processes like ring-flips, proton-transfer, or tautomerization, one has to perform studies on amino-acid side chains. A key requirement for these studies is site-selective labeling with ¹³C and/or ¹H, which is achieved in the most general way by using site-selectively ¹³C-enriched glucose (1- and 2-¹³C) as the carbon source in bacterial expression systems. Using this strategy, multiple sites in side chains, including aromatics, become site-selectively labeled and suitable for relaxation studies. Here we systematically investigate the use of site-selectively ¹³C-enriched erythrose (1-, 2-, 3- and 4-¹³C) as a suitable precursor for ¹³C labeled aromatic side chains. We quantify ¹³C incorporation in nearly all sites in all 20 amino acids and compare the results to glucose based labeling. In general the erythrose approach results in more selective labeling. While there is only a minor gain for phenylalanine and tyrosine side-chains, the ¹³C incorporation level for tryptophan is at least doubled. Additionally, the Phe ζ and Trp η2 positions become labeled. In the aliphatic side chains, labeling using erythrose yields isolated ¹³C labels for certain positions, like Ile β and His β, making these sites suitable for dynamics studies. Using erythrose instead of glucose as a source for site-selective ¹³C labeling enables unique or superior labeling for certain positions and is thereby expanding the toolbox for customized isotope labeling of amino-acid side-chains.

Comparing allosteric transitions in the domains of calmodulin through coarse-grained simulations

Calmodulin (CaM) is a ubiquitous Ca(2+)-binding protein consisting of two structurally similar domains with distinct stabilities, binding affinities, and flexibilities. We present coarse grained simulations that suggest that the mechanism for the domain's allosteric transitions between the open and closed conformations depends on subtle differences in the folded state topology of the two domains. Throughout a wide temperature range, the simulated transition mechanism of the N-terminal domain (nCaM) follows a two-state transition mechanism while domain opening in the C-terminal domain (cCaM) involves unfolding and refolding of the tertiary structure. The appearance of the unfolded intermediate occurs at a higher temperature in nCaM than it does in cCaM consistent with nCaM's higher thermal stability. Under approximate physiological conditions, the simulated unfolded state population of cCaM accounts for 10% of the population with nearly all of the sampled transitions (approximately 95%) unfolding and refolding during the conformational change. Transient unfolding significantly slows the domain opening and closing rates of cCaM, which can potentially influence its Ca(2+)-binding mechanism.

PEP-19 modulates calcium binding to calmodulin by electrostatic steering

PEP-19 is a small protein that increases the rates of Ca²⁺ binding to the C-domain of calmodulin (CaM) by an unknown mechanism. Although an IQ motif promotes binding to CaM, an acidic sequence in PEP-19 is required to modulate Ca²⁺ binding and to sensitize HeLa cells to ATP-induced Ca²⁺ release. Here, we report the NMR solution structure of a complex between PEP-19 and the C-domain of apo CaM. The acidic sequence of PEP-19 associates between helices E and F of CaM via hydrophobic interactions. This allows the acidic side chains in PEP-19 to extend toward the solvent and form a negatively charged surface that resembles a catcher's mitt near Ca²⁺ binding loop III of CaM. The topology and gradients of negative electrostatic surface potential support a mechanism by which PEP-19 increases the rate of Ca²⁺ binding to the C-domain of CaM by ‘catching' and electrostatically steering Ca²⁺ to site III.

The protein PEP-19 increases the rates of calcium binding to calmodulin. Here, the authors report the structure of PEP-19 bound to the C-terminal domain of calmodulin, and are able to propose a mechanism for the observed increased calcium association rate.

Spindle function in Xenopus oocytes involves possible nanodomain calcium signaling

Injection of dibromo-BAPTA caused immediate collapse of meiotic spindles in frog oocytes. In contrast, EGTA had no effect on the spindle or polar body emission. The disruption of spindle integrity by the fast but not slow calcium chelators suggests that meiotic spindle function in the oocytes involves nanodomain calcium signaling.

Intracellular calcium transients are a universal phenomenon at fertilization and are required for egg activation, but the exact role of Ca²⁺ in second-polar-body emission remains unknown. On the other hand, similar calcium transients have not been demonstrated during oocyte maturation, and yet, manipulating intracellular calcium levels interferes with first-polar-body emission in mice and frogs. To determine the precise role of calcium signaling in polar body formation, we used live-cell imaging coupled with temporally precise intracellular calcium buffering. We found that BAPTA-based calcium chelators cause immediate depolymerization of spindle microtubules in meiosis I and meiosis II. Surprisingly, EGTA at similar or higher intracellular concentrations had no effect on spindle function or polar body emission. Using two calcium probes containing permutated GFP and the calcium sensor calmodulin (Lck-GCaMP3 and GCaMP3), we demonstrated enrichment of the probes at the spindle but failed to detect calcium increase during oocyte maturation at the spindle or elsewhere. Finally, endogenous calmodulin was found to colocalize with spindle microtubules throughout all stages of meiosis. Our results—most important, the different sensitivities of the spindle to BAPTA and EGTA—suggest that meiotic spindle function in frog oocytes requires highly localized, or nanodomain, calcium signaling.

Structure and calcium-binding studies of calmodulin-like domain of human non-muscle α-actinin-1

The activity of several cytosolic proteins critically depends on the concentration of calcium ions. One important intracellular calcium-sensing protein is α-actinin-1, the major actin crosslinking protein in focal adhesions and stress fibers. The actin crosslinking activity of α-actinin-1 has been proposed to be negatively regulated by calcium, but the underlying molecular mechanisms are poorly understood. To address this, we determined the first high-resolution NMR structure of its functional calmodulin-like domain (CaMD) in calcium-bound and calcium-free form. These structures reveal that in the absence of calcium, CaMD displays a conformationally flexible ensemble that undergoes a structural change upon calcium binding, leading to limited rotation of the N- and C-terminal lobes around the connecting linker and consequent stabilization of the calcium-loaded structure. Mutagenesis experiments, coupled with mass-spectrometry and isothermal calorimetry data designed to validate the calcium binding stoichiometry and binding site, showed that human non-muscle α-actinin-1 binds a single calcium ion within the N-terminal lobe. Finally, based on our structural data and analogy with other α-actinins, we provide a structural model of regulation of the actin crosslinking activity of α-actinin-1 where calcium induced structural stabilisation causes fastening of the juxtaposed actin binding domain, leading to impaired capacity to crosslink actin.

Rigid Residue Scan Simulations Systematically Reveal Residue Entropic Roles in Protein Allostery

Intra-protein information is transmitted over distances via allosteric processes. This ubiquitous protein process allows for protein function changes due to ligand binding events. Understanding protein allostery is essential to understanding protein functions. In this study, allostery in the second PDZ domain (PDZ2) in the human PTP1E protein is examined as model system to advance a recently developed rigid residue scan method combining with configurational entropy calculation and principal component analysis. The contributions from individual residues to whole-protein dynamics and allostery were systematically assessed via rigid body simulations of both unbound and ligand-bound states of the protein. The entropic contributions of individual residues to whole-protein dynamics were evaluated based on covariance-based correlation analysis of all simulations. The changes of overall protein entropy when individual residues being held rigid support that the rigidity/flexibility equilibrium in protein structure is governed by the La Châtelier's principle of chemical equilibrium. Key residues of PDZ2 allostery were identified with good agreement with NMR studies of the same protein bound to the same peptide. On the other hand, the change of entropic contribution from each residue upon perturbation revealed intrinsic differences among all the residues. The quasi-harmonic and principal component analyses of simulations without rigid residue perturbation showed a coherent allosteric mode from unbound and bound states, respectively. The projection of simulations with rigid residue perturbation onto coherent allosteric modes demonstrated the intrinsic shifting of ensemble distributions supporting the population-shift theory of protein allostery. Overall, the study presented here provides a robust and systematic approach to estimate the contribution of individual residue internal motion to overall protein dynamics and allostery.

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.

Probing conformational and functional substates of calmodulin by high pressure FTIR spectroscopy: influence of Ca2+ binding and the hypervariable region of K-Ras4B

Using pressure perturbation, conformational substates of CaM could be uncovered that conceivably facilitate target recognition by exposing the required binding surfaces.

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.

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

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

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.

Impact of methionine oxidation on calmodulin structural dynamics

We have used electron paramagnetic resonance (EPR) to examine the structural impact of oxidizing specific methionine (M) side chains in calmodulin (CaM). It has been shown that oxidation of either M109 or M124 in CaM diminishes CaM regulation of the muscle calcium release channel, the ryanodine receptor (RyR), and that mutation of M to Q (glutamine) in either case produces functional effects identical to those of oxidation. Here we have used site-directed spin labeling and double electron-electron resonance (DEER), a pulsed EPR technique that measures distances between spin labels, to characterize the structural changes resulting from these mutations. Spin labels were attached to a pair of introduced cysteine residues, one in the C-lobe (T117C) and one in the N-lobe (T34C) of CaM, and DEER was used to determine the distribution of interspin distances. Ca binding induced a large increase in the mean distance, in concert with previous X-ray crystallography and NMR data, showing a closed structure in the absence of Ca and an open structure in the presence of Ca. DEER revealed additional information about CaM's structural heterogeneity in solution: in both the presence and absence of Ca, CaM populates both structural states, one with probes separated by ∼4nm (closed) and another at ∼6nm (open). Ca shifts the structural equilibrium constant toward the open state by a factor of 13. DEER reveals the distribution of interprobe distances, showing that each of these states is itself partially disordered, with the width of each population ranging from 1 to 3nm. Both mutations (M109Q and M124Q) decrease the effect of Ca on the structure of CaM, primarily by decreasing the closed-to-open equilibrium constant in the presence of Ca. We propose that Met oxidation alters CaM's functional interaction with its target proteins by perturbing this Ca-dependent structural shift.

Mechanisms of regulation of olfactory transduction and adaptation in the olfactory cilium

Olfactory adaptation is a fundamental process for the functioning of the olfactory system, but the underlying mechanisms regulating its occurrence in intact olfactory sensory neurons (OSNs) are not fully understood. In this work, we have combined stochastic computational modeling and a systematic pharmacological study of different signaling pathways to investigate their impact during short-term adaptation (STA). We used odorant stimulation and electroolfactogram (EOG) recordings of the olfactory epithelium treated with pharmacological blockers to study the molecular mechanisms regulating the occurrence of adaptation in OSNs. EOG responses to paired-pulses of odorants showed that inhibition of phosphodiesterases (PDEs) and phosphatases enhanced the levels of STA in the olfactory epithelium, and this effect was mimicked by blocking vesicle exocytosis and reduced by blocking cyclic adenosine monophosphate (cAMP)-dependent protein kinase (PKA) and vesicle endocytosis. These results suggest that G-coupled receptors (GPCRs) cycling is involved with the occurrence of STA. To gain insights on the dynamical aspects of this process, we developed a stochastic computational model. The model consists of the olfactory transduction currents mediated by the cyclic nucleotide gated (CNG) channels and calcium ion (Ca²⁺)-activated chloride (CAC) channels, and the dynamics of their respective ligands, cAMP and Ca²⁺, and it simulates the EOG results obtained under different experimental conditions through changes in the amplitude and duration of cAMP and Ca²⁺ response, two second messengers implicated with STA occurrence. The model reproduced the experimental data for each pharmacological treatment and provided a mechanistic explanation for the action of GPCR cycling in the levels of second messengers modulating the levels of STA. All together, these experimental and theoretical results indicate the existence of a mechanism of regulation of STA by signaling pathways that control GPCR cycling and tune the levels of second messengers in OSNs, and not only by CNG channel desensitization as previously thought.

Mapping central α-helix linker mediated conformational transition pathway of calmodulin via simple computational approach

The effects of intrinsic structural flexibility of calmodulin protein on the mechanism of its allosteric conformational transition are investigated in this article. Using a novel in silico approach, the conformational transition pathways of intact calmodulin as well as the isolated N- and C- terminal domains are identified and energetically characterized. It is observed that the central α-helix linker amplifies the structural flexibility of intact Ca(2+)-free calmodulin, which might facilitate the transition of the two domains. As a result, the global conformational transition of Ca(2+)-free calmodulin is initiated by the barrierless transition of two domains and proceeds through the barrier associated unwinding and bending of the central α-helix linker. The binding of Ca(2+) cations to calmodulin further increases the structural flexibility of the C-terminal domain and results in a downhill transition pathway of which all regions transit in a concerted manner. On the other hand, the separation of the N- and C-terminal domains from calmodulin protein loses the mediating function of central α-helix linker, leading to more difficult conformational transitions of both domains. The present study provides novel insights into the correlation of the integrity of protein, the structural flexibility, and the mechanism of conformational transition of proteinlike calmodulin.

Energy landscape views for interplays among folding, binding, and allostery of calmodulin domains

Ligand binding modulates the energy landscape of proteins, thus altering their folding and allosteric conformational dynamics. To investigate such interplay, calmodulin has been a model protein. Despite much attention, fully resolved mechanisms of calmodulin folding/binding have not been elucidated. Here, by constructing a computational model that can integrate folding, binding, and allosteric motions, we studied in-depth folding of isolated calmodulin domains coupled with binding of two calcium ions and associated allosteric conformational changes. First, mechanically pulled simulations revealed coexistence of three distinct conformational states: the unfolded, the closed, and the open states, which is in accord with and augments structural understanding of recent single-molecule experiments. Second, near the denaturation temperature, we found the same three conformational states as well as three distinct binding states: zero, one, and two calcium ion bound states, leading to as many as nine states. Third, in terms of the nine-state representation, we found multiroute folding/binding pathways and shifts in their probabilities with the calcium concentration. At a lower calcium concentration, "combined spontaneous folding and induced fit" occurs, whereas at a higher concentration, "binding-induced folding" dominates. Even without calcium binding, we observed that the folding pathway of calmodulin domains can be modulated by the presence of metastable states. Finally, full-length calmodulin also exhibited an intriguing coupling between two domains when applying tension.

Distinguishing unfolding and functional conformational transitions of calmodulin using ultraviolet resonance Raman spectroscopy

Calmodulin (CaM) is a ubiquitous moderator protein for calcium signaling in all eukaryotic cells. This small calcium-binding protein exhibits a broad range of structural transitions, including domain opening and folding-unfolding, that allow it to recognize a wide variety of binding partners in vivo. While the static structures of CaM associated with its various binding activities are fairly well-known, it has been challenging to examine the dynamics of transition between these structures in real-time, due to a lack of suitable spectroscopic probes of CaM structure. In this article, we examine the potential of ultraviolet resonance Raman (UVRR) spectroscopy for clarifying the nature of structural transitions in CaM. We find that the UVRR spectral change (with 229 nm excitation) due to thermal unfolding of CaM is qualitatively different from that associated with opening of the C-terminal domain in response to Ca(2+) binding. This spectral difference is entirely due to differences in tertiary contacts at the interdomain tyrosine residue Tyr138, toward which other spectroscopic methods are not sensitive. We conclude that UVRR is ideally suited to identifying the different types of structural transitions in CaM and other proteins with conformation-sensitive tyrosine residues, opening a path to time-resolved studies of CaM dynamics using Raman spectroscopy.

Population shift of binding pocket size and dynamic correlation analysis shed new light on the anticooperative mechanism of PII protein

PII protein is one of the largest families of signal transduction proteins in archaea, bacteria, and plants, controlling key processes of nitrogen assimilation. An intriguing characteristic for many PII proteins is that the three ligand binding sites exhibit anticooperative allosteric regulation. In this work, PII protein from Synechococcus elongatus, a model for cyanobacteria and plant PII proteins, is utilized to reveal the anticooperative mechanism upon binding of 2-oxoglutarate (2-OG). To this end, a method is proposed to define the binding pocket size by identifying residues that contribute greatly to the binding of 2-OG. It is found that the anticooperativity is realized through population shift of the binding pocket size in an asymmetric manner. Furthermore, a new algorithm based on the dynamic correlation analysis is developed and utilized to discover residues that mediate the anticooperative process with high probability. It is surprising to find that the T-loop, which is believed to be responsible for mediating the binding of PII with its target proteins, also takes part in the intersubunit signal transduction process. Experimental results of PII variants further confirmed the influence of T-loop on the anticooperative regulation, especially on binding of the third 2-OG. These discoveries extend our understanding of the PII T-loop from being essential in versatile binding of target protein to signal-mediating in the anticooperative allosteric regulation.

Neurogranin alters the structure and calcium binding properties of calmodulin

Neurogranin (Ng) is a member of the IQ motif class of calmodulin (CaM)-binding proteins, and interactions with CaM are its only known biological function. In this report we demonstrate that the binding affinity of Ng for CaM is weakened by Ca(2+) but to a lesser extent (2-3-fold) than that previously suggested from qualitative observations. We also show that Ng induced a >10-fold decrease in the affinity of Ca(2+) binding to the C-terminal domain of CaM with an associated increase in the Ca(2+) dissociation rate. We also discovered a modest, but potentially important, increase in the cooperativity in Ca(2+) binding to the C-lobe of CaM in the presence of Ng, thus sharpening the threshold for the C-domain to become Ca(2+)-saturated. Domain mapping using synthetic peptides indicated that the IQ motif of Ng is a poor mimetic of the intact protein and that the acidic sequence just N-terminal to the IQ motif plays an important role in reproducing Ng-mediated decreases in the Ca(2+) binding affinity of CaM. Using NMR, full-length Ng was shown to make contacts largely with residues in the C-domain of CaM, although contacts were also detected in residues in the N-terminal domain. Together, our results can be consolidated into a model where Ng contacts residues in the N- and C-lobes of both apo- and Ca(2+)-bound CaM and that although Ca(2+) binding weakens Ng interactions with CaM, the most dramatic biochemical effect is the impact of Ng on Ca(2+) binding to the C-terminal lobe of CaM.

Protein conformational exchange measured by 1H R1ρ relaxation dispersion of methyl groups

Activated dynamics plays a central role in protein function, where transitions between distinct conformations often underlie the switching between active and inactive states. The characteristic time scales of these transitions typically fall in the microsecond to millisecond range, which is amenable to investigations by NMR relaxation dispersion experiments. Processes at the faster end of this range are more challenging to study, because higher RF field strengths are required to achieve refocusing of the exchanging magnetization. Here we describe a rotating-frame relaxation dispersion experiment for (1)H spins in methyl (13)CHD2 groups, which improves the characterization of fast exchange processes. The influence of (1)H-(1)H rotating-frame nuclear Overhauser effects (ROE) is shown to be negligible, based on a comparison of R 1ρ relaxation data acquired with tilt angles of 90° and 35°, in which the ROE is maximal and minimal, respectively, and on samples containing different (1)H densities surrounding the monitored methyl groups. The method was applied to ubiquitin and the apo form of calmodulin. We find that ubiquitin does not exhibit any (1)H relaxation dispersion of its methyl groups at 10 or 25 °C. By contrast, calmodulin shows significant conformational exchange of the methionine methyl groups in its C-terminal domain, as previously demonstrated by (1)H and (13)C CPMG experiments. The present R 1ρ experiment extends the relaxation dispersion profile towards higher refocusing frequencies, which improves the definition of the exchange correlation time, compared to previous results.

Binding free-energy calculation of an ion-peptide complex by constrained dynamics

Binding free energy is the most important physical parameter that describes the binding affinity of a receptor-ligand complex. Conventionally, it was obtained based on the thermodynamic cycle or alchemical reaction. These strategies have been widely used, but they would be problematic if the receptors and/or ligands have large conformational changes during the binding processes. In this paper, we present a way to calculate the binding free energy: constrained dynamics along a fragmental and high-dimensional transition path. This method directly considers unbound states in the simulation. The application to the calmodulin loop-calcium complexes shows that it is practical and the calculated relative binding affinities are in good agreement with experimental results.

Specific 12CβD(2)12CγD(2)S13CεHD(2) isotopomer labeling of methionine to characterize protein dynamics by 1H and 13C NMR relaxation dispersion

Protein dynamics on the micro- to millisecond time scale is increasingly found to be critical for biological function, as demonstrated by numerous NMR relaxation dispersion studies. Methyl groups are excellent probes of protein interactions and dynamics because of their favorable NMR relaxation properties, which lead to sharp signals in the ¹H and ¹³C NMR spectra. Out of the six different methyl-bearing amino acid residue types in proteins, methionine plays a special role because of its extensive side-chain flexibility and the high polarizability of the sulfur atom. Methionine is over-represented in many protein–protein recognition sites, making the methyl group of this residue type an important probe of the relationships among dynamics, interactions, and biological function. Here we present a straightforward method to label methionine residues with specific ¹³CHD₂ methyl isotopomers against a deuterated background. The resulting protein samples yield NMR spectra with improved sensitivity due to the essentially 100% population of the desired ¹³CHD₂ methyl isotopomer, which is ideal for ¹H and ¹³C spin relaxation experiments to investigate protein dynamics in general and conformational exchange in particular. We demonstrate the approach by measuring ¹H and ¹³C CPMG relaxation dispersion for the nine methionines in calcium-free calmodulin (apo-CaM). The results show that the C-terminal domain, but not the N-terminal domain, of apo-CaM undergoes fast exchange between the ground state and a high-energy state. Since target proteins are known to bind specifically to the C-terminal domain of apo-CaM, we speculate that the high-energy state might be involved in target binding through conformational selection.

Conformational exchange of aromatic side chains characterized by L-optimized TROSY-selected ¹³C CPMG relaxation dispersion

Protein dynamics on the millisecond time scale commonly reflect conformational transitions between distinct functional states. NMR relaxation dispersion experiments have provided important insights into biologically relevant dynamics with site-specific resolution, primarily targeting the protein backbone and methyl-bearing side chains. Aromatic side chains represent attractive probes of protein dynamics because they are over-represented in protein binding interfaces, play critical roles in enzyme catalysis, and form an important part of the core. Here we introduce a method to characterize millisecond conformational exchange of aromatic side chains in selectively (13)C labeled proteins by means of longitudinal- and transverse-relaxation optimized CPMG relaxation dispersion. By monitoring (13)C relaxation in a spin-state selective manner, significant sensitivity enhancement can be achieved in terms of both signal intensity and the relative exchange contribution to transverse relaxation. Further signal enhancement results from optimizing the longitudinal relaxation recovery of the covalently attached (1)H spins. We validated the L-TROSY-CPMG experiment by measuring fast folding-unfolding kinetics of the small protein CspB under native conditions. The determined unfolding rate matches perfectly with previous results from stopped-flow kinetics. The CPMG-derived chemical shift differences between the folded and unfolded states are in excellent agreement with those obtained by urea-dependent chemical shift analysis. The present method enables characterization of conformational exchange involving aromatic side chains and should serve as a valuable complement to methods developed for other types of protein side chains.

Calcium-dependent folding of single calmodulin molecules

Calmodulin is the primary calcium binding protein in living cells. Its function and structure depend strongly on calcium concentration. We used single molecule force spectroscopy by optical tweezers to study the folding of calmodulin in the physiologically relevant range. We find that full-length calmodulin switches from a rich and complex folding behavior at high calcium to a simple folding pathway at apo conditions. Using truncation mutants, we studied the individual domains separately. Folding and stability of the individual domains differ significantly at low calcium concentrations. With increasing calcium, the folding rate constants increase while unfolding rate constants decrease. The complete kinetic as well as energetic behavior of both domains could be modeled using a calcium-dependent three-pathway model. We find that the dominant folding pathway at high calcium concentrations proceeds via a transition state capable of binding one calcium ion. The folding of calmodulin seems to be designed to occur fast robustly over a large range of calcium concentrations and hence energetic stabilities.

Barium ions selectively activate BK channels via the Ca2+-bowl site

Activation of Ca(2+)-dependent BK channels is increased via binding of micromolar Ca(2+) to two distinct high-affinity sites per BK α-subunit. One site, termed the Ca(2+) bowl, is embedded within the second RCK domain (RCK2; regulator of conductance for potassium) of each α-subunit, while oxygen-containing residues in the first RCK domain (RCK1) have been linked to a separate Ca(2+) ligation site. Although both sites are activated by Ca(2+) and Sr(2+), Cd(2+) selectively favors activation via the RCK1 site. Divalent cations of larger ionic radius than Sr(2+) are thought to be ineffective at activating BK channels. Here we show that Ba(2+), better known as a blocker of K(+) channels, activates BK channels and that this effect arises exclusively from binding at the Ca(2+)-bowl site. Compared with previous estimates for Ca(2+) bowl-mediated activation by Ca(2+), the affinity of Ba(2+) to the Ca(2+) bowl is reduced about fivefold, and coupling of binding to activation is reduced from ∼3.6 for Ca(2+) to about ∼2.8 for Ba(2+). These results support the idea that ionic radius is an important determinant of selectivity differences among different divalent cations observed for each Ca(2+)-binding site.

Conformational flexibility and the mechanisms of allosteric transitions in topologically similar proteins

Conformational flexibility plays a central role in allosteric transition of proteins. In this paper, we extend the analysis of our previous study [S. Tripathi and J. J. Portman, Proc. Natl. Acad. Sci. U.S.A. 106, 2104 (2009)] to investigate how relatively minor structural changes of the meta-stable states can significantly influence the conformational flexibility and allosteric transition mechanism. We use the allosteric transitions of the domains of calmodulin as an example system to highlight the relationship between the transition mechanism and the inter-residue contacts present in the meta-stable states. In particular, we focus on the origin of transient local unfolding (cracking), a mechanism that can lower free energy barriers of allosteric transitions, in terms of the inter-residue contacts of the meta-stable states and the pattern of local strain that develops during the transition. We find that the magnitude of the local strain in the protein is not the sole factor determining whether a region will ultimately crack during the transition. These results emphasize that the residue interactions found exclusively in one of the two meta-stable states is the key in understanding the mechanism of allosteric conformational change.

Transient, sparsely populated compact states of apo and calcium-loaded calmodulin probed by paramagnetic relaxation enhancement: interplay of conformational selection and induced fit

Calmodulin (CaM) is the universal calcium sensor in eukaryotes, regulating the function of numerous proteins. Crystallography and NMR show that free CaM-4Ca(2+) exists in an extended conformation with significant interdomain separation, but clamps down upon target peptides to form a highly compact structure. NMR has revealed substantial interdomain motions in CaM-4Ca(2+), enabled by a flexible linker. In one instance, CaM-4Ca(2+) has been crystallized in a compact configuration; however, no direct evidence for transient interdomain contacts has been observed in solution, and little is known about how large-scale interdomain motions contribute to biological function. Here, we use paramagnetic relaxation enhancement (PRE) to characterize transient compact states of free CaM that are too sparsely populated to observe by traditional NMR methods. We show that unbound CaM samples a range of compact structures, populated at 5-10%, and that Ca(2+) dramatically alters the distribution of these configurations in favor of states resembling the peptide-bound structure. In the absence of Ca(2+), the target peptide binds only to the C-terminal domain, and the distribution of compact states is similar with and without peptide. These data suggest an alternative pathway of CaM action in which CaM remains associated with its kinase targets even in the resting state. Only CaM-4Ca(2+), however, shows an innate propensity to form the physiologically active compact structures, suggesting that Ca(2+) activates CaM not only through local structural changes within each domain but also through more global remodeling of interdomain interactions. Thus, these findings illustrate the subtle interplay between conformational selection and induced fit.

A new concept to reveal protein dynamics based on energy dissipation

Protein dynamics is essential for its function, especially for intramolecular signal transduction. In this work we propose a new concept, energy dissipation model, to systematically reveal protein dynamics upon effector binding and energy perturbation. The concept is applied to better understand the intramolecular signal transduction during allostery of enzymes. The E. coli allosteric enzyme, aspartokinase III, is used as a model system and special molecular dynamics simulations are designed and carried out. Computational results indicate that the number of residues affected by external energy perturbation (i.e. caused by a ligand binding) during the energy dissipation process shows a sigmoid pattern. Using the two-state Boltzmann equation, we define two parameters, the half response time and the dissipation rate constant, which can be used to well characterize the energy dissipation process. For the allostery of aspartokinase III, the residue response time indicates that besides the ACT2 signal transduction pathway, there is another pathway between the regulatory site and the catalytic site, which is suggested to be the β15-αK loop of ACT1. We further introduce the term “protein dynamical modules” based on the residue response time. Different from the protein structural modules which merely provide information about the structural stability of proteins, protein dynamical modules could reveal protein characteristics from the perspective of dynamics. Finally, the energy dissipation model is applied to investigate E. coli aspartokinase III mutations to better understand the desensitization of product feedback inhibition via allostery. In conclusion, the new concept proposed in this paper gives a novel holistic view of protein dynamics, a key question in biology with high impacts for both biotechnology and biomedicine.