Unwinding the helical linker of calcium-loaded calmodulin: a molecular dynamics study.

The fold of calmodulin (CaM) consists of two globular domains connected by a helical segment (the linker), whose conformational properties play a crucial role for the protein's molecular recognition processes. Here we investigate the structural properties of the linker by performing a 11.5 ns molecular dynamics (MD) simulation of calcium-loaded human CaM in aqueous solution. The calculations are based on the AMBER force field. The calculated S2 order parameters are in good accord with NMR data: The structure of the linker in our simulations is much more flexible than that emerging from the Homo sapiens X-ray structure, consistently with the helix unwinding observed experimentally in solution. This process occurs spontaneously in a nanosecond timescale, as observed also in a very recent simulation based on the GROMOS force field. A detailed description of the mechanism that determines the linker unwinding is provided, in which electrostatic contacts between the two globular domains play a critical role. The orientation of the domains emerging from our MD calculations is consistent both with former X-ray scattering data and a recent NMR work. Based on our findings, a rationale for the experimentally measured entropy cost associated to binding to the protein's cellular partners is also given.

Calmodulin tagging provides a general method of using lanthanide induced magnetic field orientation to observe residual dipolar couplings in proteins in solution.

A general method is presented for magnetic field alignment of proteins in solution. By tagging a target protein with calmodulin saturated with paramagnetic lanthanide ions it is possible to measure substantial residual dipolar couplings (RDC) whilst minimising the effects of pseudocontact shifts on the target protein. A construct was made consisting of a calmodulin-binding peptide (M13 from sk-MLCK) attached to a target protein, dihydrofolate reductase in this case. The engineered protein binds tightly to calmodulin saturated with terbium, a paramagnetic lanthanide ion. By using only a short linker region between the M13 and the target protein, some of the magnetic field alignment induced in the CaM(Tb3+)4 is effectively transmitted to the target protein (DHFR). 1H-15N HSQC IPAP experiments on the tagged complex containing 15N-labelled DHFR-M13 protein and unlabelled CaM(Tb3+)4 allow one to measure RDC contributions in the aligned complex. RDC values in the range +4.0 to -7.4 Hz were measured at 600 MHz. Comparisons of 1H-15N HSQC spectra of 15N-DHFR-M13 alone and its complexes with CaM(Ca2+)4 and CaM(Tb3+)4 indicated that (i) the structure of the target protein is not affected by the complex formation and (ii) the spectra of the target protein are not seriously perturbed by pseudocontact shifts. The use of a relatively large tagging group (CaM) allows us to use a lanthanide ion with a very high magnetic susceptibility anisotropy (such as Tb3+) to give large alignments while maintaining relatively long distances from the target protein nuclei (and hence giving only small pseudocontact shift contributions).

Calcium-induced refolding of the calmodulin V136G mutant studied by NMR spectroscopy: evidence for interaction between the two globular domains.

The Ca(2+) titration of the (15)N-labeled mutant V136G calmodulin has been monitored using (1)H-(15)N HSQC NMR spectra. Up to a [Ca(2+)]/[CaM] ratio of 2, the Ca(2+) ions bind predominantly to sites I and II on the N-domain in contrast with the behavior of the wild-type calmodulin where the C-terminal domain has the higher affinity for Ca(2+). Surprisingly, the Ca(2+)-binding affinity for the N-domain in the mutant calmodulin is greater than that for the N-domain in the wild-type protein. The mutated C-domain is observed as a mixture of unfolded, partially folded (site III occupied), and native-like folded (sites III and IV occupied) conformations, with relative populations dependent on the [Ca(2+)]/[CaM] ratio. The occupancy of site III independently of site IV in this mutant shows that the cooperativity of Ca(2+) binding in the C-domain is mediated by the integrity of the domain structure. Several NH signals from residues in the Ca(2+)-bound N-domain appear as two signals during the Ca(2+) titration indicating separate species in slow exchange, and it can be deduced that these result from the presence and absence of interdomain interactions in the mutant. It is proposed that an unfolded part of the mutated C-domain interacts with sites on the N-domain that normally bind to target proteins. This would also account for the increase in the Ca(2+) affinity for the N-domain in the mutant compared with the wild-type calmodulin. The results therefore show the wide-ranging effects of a point mutation in a single Ca(2+)-binding site, providing details of the involvement of individual residues in the calcium-induced folding reactions.

NMR approaches for monitoring domain orientations in calcium-binding proteins in solution using partial replacement of Ca2+ by Tb3+.

This work shows that the partial replacement of diamagnetic Ca2+ by paramagnetic Tb3+ in Ca2+/calmodulin systems in solution allows the measurement of interdomain NMR pseudocontact shifts and leads to magnetic alignment of the molecule such that significant residual dipolar couplings can be measured. Both these parameters can be used to provide structural information. Species in which Tb3+ ions are bound to only one domain of calmodulin (the N-domain) and Ca2+ ions to the other (the C-domain) provide convenient systems for measuring these parameters. The nuclei in the C-domain experience the local magnetic field induced by the paramagnetic Tb3+ ions bound to the other domain at distances of over 40 A from the Tb3+ ion, shifting the resonances for these nuclei. In addition, the Tb3+ ions bound to the N-domain of calmodulin greatly enhance the magnetic susceptibility anisotropy of the molecule so that a certain degree of alignment is produced due to interaction with the external magnetic field. In this way, dipolar couplings between nuclear spins are not averaged to zero due to solution molecular tumbling and yield dipolar coupling contributions to, for example, the one-bond 15N-1H splittings of up to 17 Hz in magnitude. The degree of alignment of the C-domain will also depend on the degree of orientational freedom of this domain with respect to the N-domain containing the Tb3+ ions. Pseudocontact shifts for NH groups and 1H-15N residual dipolar couplings for the directly bonded atoms have been measured for calmodulin itself, where the domains have orientational freedom, and for the complex of calmodulin with a target peptide from skeletal muscle myosin light chain kinase, where the domains have fixed orientations with respect to each other. The simultaneous measurements of these parameters for systems with domains in fixed orientations show great potential for the determination of the relative orientation of the domains.

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

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

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

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

Conformational studies of the beta-subunit of the high affinity IgE receptor: circular dichroism and molecular modelling.

The receptor with high affinity for immunoglobulin E (Fc epsilon RI) on mast cells and basophils plays an important role in mediating many of the pathophysiological phenomena associated with allergy. Fc epsilon RI is a tetrameric complex, alpha beta gamma2, of non-covalently attached subunits: one IgE-binding alpha-subunit with the binding site in the extracellular part of the chain, one beta-subunit and a dimer of disulphide linked gamma-subunits. In the present work, prediction of the three-dimensional structure of the four membrane-spanning segments of the beta-subunit has been achieved using rules of helix-helix packing arrangements and molecular dynamics calculations. It yielded a four-helix bundle with specific Van der Waals interactions between the helices. This four-helix bundle was used as a framework upon which to calculate the conformation of the beta-subunit excluding the C and N terminal cytoplasmic tails, but including the three chains that connect the four helices in the bundle. Separately, these synthetic 11, 17 and 29 residue bridge peptides were examined by circular dichroism (CD) spectroscopy and a degree of alpha-helical content in these bridge peptides was found. Additional molecular modelling of the bridge peptides indicate the central residues of these as the location of the helical moieties. Finally, in the model proposed for the beta-subunit, for each pair of consecutive transmembrane (TM) helices and its bridge peptide, a helix-loop-helix-loop-helix motif was found.

Spectroscopy and modelling of the cytoplasmic domain of the gamma-subunit of the high affinity immunoglobulin E receptor.

The high affinity receptor for IgE, Fc epsilon RI, is responsible for immediate hypersensitivity reactions. In rodents Fc epsilon RI is a tetrameric complex, alpha beta gamma 2 of non-covalently attached subunits: one IgE-binding alpha subunit with the binding site in the extracellular part of the chain, one beta-subunit and a dimer of disulphide linked gamma-subunits. Although there is an increasing evidence that the gamma-subunit chains are important signalling proteins that appear to function through a common Tyr-Leu-Tyr-Leu amino acid motif present in their cytoplasmic tails, which link the ligand binding specificity of their associated chains to signal transduction pathways, many questions related to conformation and function of this subunit remain to be answered. In the present work, the 36-residue cytoplasmic domain of the gamma-subunit has been synthesized and conformational studies by the combined use of Fourier transform infrared (FTIR), circular dichroism (CD) and nuclear magnetic resonance (NMR) have been performed. Based on the constraints found by these methods, conformational models of the cytoplasmic tail of the gamma-subunit are proposed and discussed.