Accurate Solution Structures of Proteins from X-ray Data and a Minimal Set of NMR Data: Calmodulin−Peptide Complexes As Examples

Exploring regions of conformational space occupied by two-domain proteins

The presence of heterogeneity in the interdomain arrangement of several biomolecules is required for their function. Here we present a method to obtain crucial clues to distinguish between different kinds of protein conformational distributions based on experimental NMR data. The method explores subregions of the conformational space and provides both upper and lower bounds of probability for the system to be in each subregion.

The many structural faces of calmodulin: a multitasking molecular jackknife

Calmodulin (CaM) is a highly conserved protein and a crucial calcium sensor in eukaryotes. CaM is a regulator of hundreds of diverse target proteins. A wealth of studies has been carried out on the structure of CaM, both in the unliganded form and in complexes with target proteins and peptides. The outcome of these studies points toward a high propensity to attain various conformational states, depending on the binding partner. The purpose of this review is to provide examples of different conformations of CaM trapped in the crystal state. In addition, comparisons are made to corresponding studies in solution. The different CaM conformations in crystal structures are also compared based on the positions of the metal ions bound to their EF hands, in terms of distances, angles, and pseudo-torsion angles. Possible caveats and artifacts in CaM crystal structures are discussed, as well as the possibilities of trapping biologically relevant CaM conformations in the crystal state.

Simultaneous use of solution NMR and X-ray data in REFMAC5 for joint refinement/detection of structural differences

The program REFMAC5 from CCP4 was modified to allow the simultaneous use of X-ray crystallographic data and paramagnetic NMR data (pseudocontact shifts and self-orientation residual dipolar couplings) and/or diamagnetic residual dipolar couplings. Incorporation of these long-range NMR restraints in REFMAC5 can reveal differences between solid-state and solution conformations of molecules or, in their absence, can be used together with X-ray crystallographic data for structural refinement. Since NMR and X-ray data are complementary, when a single structure is consistent with both sets of data and still maintains reasonably `ideal' geometries, the reliability of the derived atomic model is expected to increase. The program was tested on five different proteins: the catalytic domain of matrix metalloproteinase 1, GB3, ubiquitin, free calmodulin and calmodulin complexed with a peptide. In some cases the joint refinement produced a single model consistent with both sets of observations, while in other cases it indicated, outside the experimental uncertainty, the presence of different protein conformations in solution and in the solid state.

SedNMR: a web tool for optimizing sedimentation of macromolecular solutes for SSNMR

We have proposed solid state NMR (SSNMR) of sedimented solutes as a novel approach to sample preparation for biomolecular SSNMR without crystallization or other sample manipulations. The biomolecules are confined by high gravity--obtained by centrifugal forces either directly in a SSNMR rotor or in a ultracentrifugal device--into a hydrated non-crystalline solid suitable for SSNMR investigations. When gravity is removed, the sample reverts to solution and can be treated as any solution NMR sample. We here describe a simple web tool to calculate the relevant parameters for the success of the experiment.

Interdomain dynamics explored by paramagnetic NMR

An ensemble-based approach is presented to explore the conformational space sampled by a multidomain protein showing moderate interdomain dynamics in terms of translational and rotational motions. The strategy was applied on a complex of calmodulin (CaM) with the IQ-recognition motif from the voltage-gated calcium channel Ca(v)1.2 (IQ), which adopts three different interdomain orientations in the crystal. The N60D mutant of calmodulin was used to collect pseudocontact shifts and paramagnetically induced residual dipolar couplings for six different lanthanide ions. Then, starting from the crystal structure, pools of conformations were generated by free MD. We found the three crystal conformations in solution, but four additional MD-derived conformations had to be included into the ensemble to fulfill all the paramagnetic data and cross-validate optimally against unused paramagnetic data. Alternative approaches led to similar ensembles. Our "ensemble" approach is a simple and efficient tool to probe and describe the interdomain dynamics and represents a general method that can be used to provide a proper ensemble description of multidomain proteins.

Modulating alignment of membrane proteins in liquid-crystalline and oriented gel media by changing the size and charge of phospholipid bicelles

We demonstrate that alignment of a structured peptide or small protein solubilized in mixed phospholipid:detergent micelles or bicelles, when embedded in a compressed gel or liquid crystalline medium, can be altered by either changing the phospholipid aggregate shape, charge, or both together. For the hemagglutinin fusion peptide solubilized in bicelles, we show that bicelle shape and charge do not change its helical hairpin structure but impact its alignment relative to the alignment medium, both in charged compressed acrylamide gel and in liquid crystalline d(GpG). The method can be used to generate sets of residual dipolar couplings that correspond to orthogonal alignment tensors, and holds promise for high-resolution structural refinement and dynamic mapping of membrane proteins.

On the use of ultracentrifugal devices for sedimented solute NMR

We have recently proposed sedimented solute NMR (SedNMR) as a solid-state method to access biomolecules without the need of crystallization or other sample manipulation. The drawback of SedNMR is that samples are intrinsically diluted and this is detrimental for the signal intensity. Ultracentrifugal devices can be used to increase the amount of sample inside the rotor, overcoming the intrinsic sensitivity limitation of the method. We designed two different devices and we here report the directions for using such devices and the relevant equations for determining the parameters for sedimentation.

Interactions of calmodulin with death-associated protein kinase peptides: experimental and modeling studies

We have studied the interactions between calmodulin (CaM) and three target peptides from the death-associated protein kinase (DAPK) protein family using both experimental and modeling methods, aimed at determining the details of the underlying biological regulation mechanisms. Experimentally, calorimetric binding free energies were determined for the complexes of CaM with peptides representing the DAPK2 wild-type and S308D mutant, as well as DAPK1. The observed affinity of CaM was very similar for all three studied peptides. The DAPK2 and DAPK1 peptides differ significantly in sequence and total charge, while the DAPK2 S308D mutant is designed to model the effects of DAPK2 Ser308 phosphorylation. The crystal structure of the CaM-DAPK2 S308D mutant peptide is also reported. The structures of CaM-DAPK peptide complexes present a mode of CaM-kinase interaction, in which bulky hydrophobic residues at positions 10 and 14 are both bound to the same hydrophobic cleft. To explain the microscopic effects underlying these interactions, we performed free energy calculations based on the approximate MM-PBSA approach. For these highly charged systems, standard MM-PBSA calculations did not yield satisfactory results. We proposed a rational modification of the approach which led to reasonable predictions of binding free energies. All three complexes are strongly stabilized by two effects: electrostatic interactions and buried surface area. The strong favorable interactions are to a large part compensated by unfavorable entropic terms, in which vibrational entropy is the largest contributor. The electrostatic component of the binding free energy followed the trend of the overall peptide charge, with strongest interactions for DAPK1 and weakest for the DAPK2 mutant. The electrostatics was dominated by interactions of the positively charged residues of the peptide with the negatively charged residues of CaM. The nonpolar binding free energy was comparable for all three peptides, the largest contribution coming from the Trp305. About two-thirds of the buried surface area corresponds to nonpolar residues, showing that hydrophobic interactions play an important role in these CaM-peptide complexes. The simulation results agree with the experimental data in predicting a small effect of the S308D mutation on CaM interactions with DAPK2, suggesting that this mutation is not a good model for the S308 phosphorylation.

Recognition pliability is coupled to structural heterogeneity: a calmodulin intrinsically disordered binding region complex

Protein interactions within regulatory networks should adapt in a spatiotemporal-dependent dynamic environment, in order to process and respond to diverse and versatile cellular signals. However, the principles governing recognition pliability in protein complexes are not well understood. We have investigated a region of the intrinsically disordered protein myelin basic protein (MBP(145-165)) that interacts with calmodulin, but that also promiscuously binds other biomolecules (membranes, modifying enzymes). To characterize this interaction, we implemented an NMR spectroscopic approach that calculates, for each conformation of the complex, the maximum occurrence based on recorded pseudocontact shifts and residual dipolar couplings. We found that the MBP(145-165)-calmodulin interaction is characterized by structural heterogeneity. Quantitative comparative analysis indicated that distinct conformational landscapes of structural heterogeneity are sampled for different calmodulin-target complexes. Such structural heterogeneity in protein complexes could potentially explain the way that transient and promiscuous protein interactions are optimized and tuned in complex regulatory networks.

Paramagnetic relaxation enhancement for the characterization of the conformational heterogeneity in two-domain proteins

Multidomain proteins are often composed of rigid domains that can reorient in solution more or less freely. Calmodulin (CaM) is a two domain protein which can experience a large degree of conformational freedom thanks to a mobile linker connecting the N-terminal and C-terminal domains of the protein. The maximum occurrences (MOs) of the possible protein conformations have been analyzed using the paramagnetic relaxation enhancements (PREs) induced by a gadolinium(III) ion together with the paramagnetic pseudocontact shift and residual dipolar coupling restraints measured in the presence of terbium(III), thulium(III) or dysprosium(III) ions. The results suggest that the PREs provide complementary information useful for improving the description of the conformational heterogeneity of the protein. The data, acquired at 298 K and at 278 K, suggest that compact conformations are disfavoured by decreasing the temperature.

Convenient method for resolving degeneracies due to symmetry of the magnetic susceptibility tensor and its application to pseudo contact shift-based protein-protein complex structure determination

Pseudo contact shifts (PCSs) induced by paramagnetic lanthanide ions fixed in a protein frame provide long-range distance and angular information, and are valuable for the structure determination of protein-protein and protein-ligand complexes. We have been developing a lanthanide-binding peptide tag (hereafter LBT) anchored at two points via a peptide bond and a disulfide bond to the target proteins. However, the magnetic susceptibility tensor displays symmetry, which can cause multiple degenerated solutions in a structure calculation based solely on PCSs. Here we show a convenient method for resolving this degeneracy by changing the spacer length between the LBT and target protein. We applied this approach to PCS-based rigid body docking between the FKBP12-rapamycin complex and the mTOR FRB domain, and demonstrated that degeneracy could be resolved using the PCS restraints obtained from two-point anchored LBT with two different spacer lengths. The present strategy will markedly increase the usefulness of two-point anchored LBT for protein complex structure determination.

Structure and dynamics of calmodulin (CaM) bound to nitric oxide synthase peptides: effects of a phosphomimetic CaM mutation

Nitric oxide synthase (NOS) plays a major role in a number of key physiological and pathological processes. Knowledge of how this is regulated is important. The small acidic calcium binding protein, calmodulin (CaM), is required to fully activate the enzyme. The exact mechanism of how CaM activates NOS is not fully understood. Studies have shown CaM to act like a switch that causes a conformational change in NOS to allow for the transfer of an electron between the reductase and oxygenase domains through a process that is thought to be highly dynamic. To investigate the dynamic properties of CaM-NOS interactions, we determined the solution structure of CaM bound to the inducible NOS (iNOS) and endothelial NOS (eNOS) CaM binding region peptides. In addition, we investigated the effect of CaM phosphorylation. Tyrosine 99 (Y99) of CaM is reported to be phosphorylated in vivo. We have produced a phosphomimetic Y99E CaM to investigate the structural and functional effects that the phosphorylation of this residue may have on nitric oxide production. All three mammalian NOS isoforms were included in the investigation. Our results show that a phosphomimetic Y99E CaM significantly reduces the maximal synthase activity of eNOS by 40% while having little effect on nNOS or iNOS activity. A comparative nuclear magnetic resonance study between phosphomimetic Y99E CaM and wild-type CaM bound to the eNOS CaM binding region peptide was performed. This investigation provides important insights into how the increased electronegativity of a phosphorylated CaM protein affects the binding, dynamics, and activation of the NOS enzymes.

Solid-state NMR crystallography through paramagnetic restraints

Pseudocontact shifts (PCSs) measured by solid-state NMR spectroscopy (SS-NMR) on microcrystalline powders of a paramagnetic metalloprotein permit NMR crystallography. Along with other restraints for SS-NMR experiments, the protein molecular structure as well as the correct crystal packing are obtained.

Maximum occurrence analysis of protein conformations for different distributions of paramagnetic metal ions within flexible two-domain proteins

Multidomain proteins are composed of rigid domains connected by (flexible) linkers. Therefore, the domains may experience a large degree of reciprocal reorientation. Pseudocontact shifts and residual dipolar couplings arising from one or more paramagnetic metals successively placed in a single metal binding site in the protein can be used as restraints to assess the degree of mobility of the different domains. They can be used to determine the maximum occurrence (MO) of each possible protein conformation, i.e. the maximum weight that such conformations can have independently of the real structural ensemble, in agreement with the provided restraints. In the case of two-domain proteins, the metal ions can be placed all in the same domain, or distributed between the two domains. It has been demonstrated that the quantity of independent information for the characterization of the system is larger when all metals are bound in the same domain. At the same time, it has been shown that there are practical advantages in placing the metals in different domains. Here, it is shown that distributing the metals between the domains provides a tool for defining a coefficient of compatibility among the restraints obtained from different metals, without a significant decrease of the capability of the MO values to discriminate among conformations with different weights.

Protein structure determination from pseudocontact shifts using ROSETTA

Paramagnetic metal ions generate pseudocontact shifts (PCSs) in nuclear magnetic resonance spectra that are manifested as easily measurable changes in chemical shifts. Metals can be incorporated into proteins through metal binding tags, and PCS data constitute powerful long-range restraints on the positions of nuclear spins relative to the coordinate system of the magnetic susceptibility anisotropy tensor (Δχ-tensor) of the metal ion. We show that three-dimensional structures of proteins can reliably be determined using PCS data from a single metal binding site combined with backbone chemical shifts. The program PCS-ROSETTA automatically determines the Δχ-tensor and metal position from the PCS data during the structure calculations, without any prior knowledge of the protein structure. The program can determine structures accurately for proteins of up to 150 residues, offering a powerful new approach to protein structure determination that relies exclusively on readily measurable backbone chemical shifts and easily discriminates between correctly and incorrectly folded conformations.

Temperature and solvent dependence of the dynamical landscape of tau protein conformations

We report the variation with temperature of the ensemble distribution of conformations spanned by the tau protein in its dynamical states measured by small-angle X-ray scattering (SAXS) using synchrotron radiation. The SAXS data show a clear temperature variation of the distribution of occupied protein conformations from 293 to 318 K. More conformations with a smaller radius of gyration are occupied at higher temperature. The protein-solvent interactions are shown by computer simulation to be essential for controlling the dynamics of protein conformations, providing evidence for the key role of water solvent in the protein dynamics, as proposed by Giorgio Careri.

Narrowing the conformational space sampled by two-domain proteins with paramagnetic probes in both domains

Calmodulin is a two-domain protein which in solution can adopt a variety of conformations upon reorientation of its domains. The maximum occurrence (MO) of a set of calmodulin conformations that are representative of the overall conformational space possibly sampled by the protein, has been calculated from the paramagnetism-based restraints. These restraints were measured after inclusion of a lanthanide binding tag in the C-terminal domain to supplement the data obtained by substitution of three paramagnetic lanthanide ions to the calcium ion in the second calcium binding loop of the N-terminal domain. The analysis shows that the availability of paramagnetic restraints arising from metal ions placed on both domains, reduces the MO of the conformations to different extents, thereby helping to identify those conformations that can be mostly sampled by the protein.

Molecular mechanisms of calmodulin action on TRPV5 and modulation by parathyroid hormone

The epithelial Ca(2+) channel transient receptor potential vanilloid 5 (TRPV5) constitutes the apical entry gate for active Ca(2+) reabsorption in the kidney. Ca(2+) influx through TRPV5 induces rapid channel inactivation, preventing excessive Ca(2+) influx. This inactivation is mediated by the last ∼30 residues of the carboxy (C) terminus of the channel. Since the Ca(2+)-sensing protein calmodulin has been implicated in Ca(2+)-dependent regulation of several TRP channels, the potential role of calmodulin in TRPV5 function was investigated. High-resolution nuclear magnetic resonance (NMR) spectroscopy revealed a Ca(2+)-dependent interaction between calmodulin and a C-terminal fragment of TRPV5 (residues 696 to 729) in which one calmodulin binds two TRPV5 C termini. The TRPV5 residues involved in calmodulin binding were mutated to study the functional consequence of releasing calmodulin from the C terminus. The point mutants TRPV5-W702A and TRPV5-R706E, lacking calmodulin binding, displayed a strongly diminished Ca(2+)-dependent inactivation compared to wild-type TRPV5, as demonstrated by patch clamp analysis. Finally, parathyroid hormone (PTH) induced protein kinase A (PKA)-dependent phosphorylation of residue T709, which diminished calmodulin binding to TRPV5 and thereby enhanced channel open probability. The TRPV5-W702A mutant exhibited a significantly increased channel open probability and was not further stimulated by PTH. Thus, calmodulin negatively modulates TRPV5 activity, which is reversed by PTH-mediated channel phosphorylation.

Structure of the dimeric autoinhibited conformation of DAPK2, a pro-apoptotic protein kinase

The death-associated protein kinase (DAPK) family has been characterized as a group of pro-apoptotic serine/threonine kinases that share specific structural features in their catalytic kinase domain. Two of the DAPK family members, DAPK1 and DAPK2, are calmodulin-dependent protein kinases that are regulated by oligomerization, calmodulin binding, and autophosphorylation. In this study, we have determined the crystal and solution structures of murine DAPK2 in the presence of the autoinhibitory domain, with and without bound nucleotides in the active site. The crystal structure shows dimers of DAPK2 in a conformation that is not permissible for protein substrate binding. Two different conformations were seen in the active site upon the introduction of nucleotide ligands. The monomeric and dimeric forms of DAPK2 were further analyzed for solution structure, and the results indicate that the dimers of DAPK2 are indeed formed through the association of two apposed catalytic domains, as seen in the crystal structure. The structures can be further used to build a model for DAPK2 autophosphorylation and to compare with closely related kinases, of which especially DAPK1 is an actively studied drug target. Our structures also provide a model for both homodimerization and heterodimerization of the catalytic domain between members of the DAPK family. The fingerprint of the DAPK family, the basic loop, plays a central role in the dimerization of the kinase domain.

Conformational space of flexible biological macromolecules from average data

The concept of maximum occurrence (MO), i.e., the maximum percent of time that flexible proteins can spend in any given conformation, is introduced, and a rigorous method is developed to extensively sample the conformational space and to construct MO maps from experimental data. The method is tested in a case study, the flexible two-domain protein calmodulin (CaM), using SAXS and NMR data (i.e., pseudocontact shifts and self-orientation residual dipolar couplings arising from the presence of paramagnetic lanthanide ions), revealing that the "closed" and "fully extended" conformations trapped in the crystalline forms of CaM have MOs of only 5 and 15%, respectively. Compact conformations in general have small MOs, whereas some extended conformations have MO as high as 35%, strongly suggesting these conformations to be most abundant in solution. The method is universally applicable as it requires only standard SAXS data and specific NMR data on lanthanide derivatives of the protein (using native metal sites or lanthanide tagging). The computer program is publicly available using the grid computing infrastructure through the authors' Web portal.

Validation of a lanthanide tag for the analysis of protein dynamics by paramagnetic NMR spectroscopy

Paramagnetic lanthanide tags potentially can enhance the effects of microsecond to millisecond dynamics in proteins on NMR signals and provide structural information on lowly populated states encoded in the pseudocontact shifts. We have investigated the microsecond to millisecond mobility of a two-point attached lanthanide tag, CLaNP-5, using paramagnetic (1)H CPMG relaxation dispersion methods. CLaNP-5 loaded with Lu(3+), Yb(3+), or Tm(3+) was attached to three sites on the surface of two proteins, pseudoazurin and cytochrome c. The paramagnetic center causes large relaxation dispersion effects for two attachment sites, suggesting that local dynamics of the protein at the attachment site causes mobility of the paramagnetic center. At one site the relaxation dispersions are small and limited to the immediate environment of the tag. It is concluded that paramagnetic relaxation dispersion could represent a sensitive method to probe protein dynamics. However, the selection of a rigid attachment site is of critical importance.

Solution structure of the 128 kDa enzyme I dimer from Escherichia coli and its 146 kDa complex with HPr using residual dipolar couplings and small- and wide-angle X-ray scattering

The solution structures of free Enzyme I (EI, ∼128 kDa, 575 × 2 residues), the first enzyme in the bacterial phosphotransferase system, and its complex with HPr (∼146 kDa) have been solved using novel methodology that makes use of prior structural knowledge (namely, the structures of the dimeric EIC domain and the isolated EIN domain both free and complexed to HPr), combined with residual dipolar coupling (RDC), small- (SAXS) and wide- (WAXS) angle X-ray scattering and small-angle neutron scattering (SANS) data. The calculational strategy employs conjoined rigid body/torsion/Cartesian simulated annealing, and incorporates improvements in calculating and refining against SAXS/WAXS data that take into account complex molecular shapes in the description of the solvent layer resulting in a better representation of the SAXS/WAXS data. The RDC data orient the symmetrically related EIN domains relative to the C(2) symmetry axis of the EIC dimer, while translational, shape, and size information is provided by SAXS/WAXS. The resulting structures are independently validated by SANS. Comparison of the structures of the free EI and the EI-HPr complex with that of the crystal structure of a trapped phosphorylated EI intermediate reveals large (∼70-90°) hinge body rotations of the two subdomains comprising the EIN domain, as well as of the EIN domain relative to the dimeric EIC domain. These large-scale interdomain motions shed light on the structural transitions that accompany the catalytic cycle of EI.

Electron-nuclear interactions as probes of domain motion in proteins

Long range interactions between nuclear spins and paramagnetic ions can serve as a sensitive monitor of internal motion of various parts of proteins, including functional loops and separate domains. In the case of interdomain motion, the interactions between the ion and NMR-observable nuclei are modulated in direction and magnitude mainly by a combination of overall and interdomain motions. The effects on observable parameters such as paramagnetic relaxation enhancement (PRE) and pseudocontact shift (PCS) can, in principle, be used to characterize motion. These parameters are frequently used for the purpose of structural refinements. However, their use to probe actual domain motions is less common and is lacking a proper theoretical treatment from a motional perspective. In this work, a suitable spin Hamiltonian is incorporated in a two body diffusion model to produce the time correlation function for the nuclear spin-paramagnetic ion interactions. Simulated observables for nuclei in different positions with respect to the paramagnetic ion are produced. Based on these simulations, it demonstrated that both the PRE and the PCS can be very sensitive probes of domain motion. Results for different nuclei within the protein sense different aspects of the motions. Some are more sensitive to the amplitude of the internal motion, others are more sensitive to overall diffusion rates, allowing separation of these contributions. Experimentally, the interaction strength can also be tuned by substitution of different paramagnetic ions or by varying magnetic field strength (in the case of lanthanides) to allow the use of more detailed diffusion models without reducing the reliability of data fitting.

PCS-based structure determination of protein-protein complexes

A simple and fast nuclear magnetic resonance method for docking proteins using pseudo-contact shift (PCS) and (1)H(N)/(15)N chemical shift perturbation is presented. PCS is induced by a paramagnetic lanthanide ion that is attached to a target protein using a lanthanide binding peptide tag anchored at two points. PCS provides long-range (approximately 40 A) distance and angular restraints between the lanthanide ion and the observed nuclei, while the (1)H(N)/(15)N chemical shift perturbation data provide loose contact-surface information. The usefulness of this method was demonstrated through the structure determination of the p62 PB1-PB1 complex, which forms a front-to-back 20 kDa homo-oligomer. As p62 PB1 does not intrinsically bind metal ions, the lanthanide binding peptide tag was attached to one subunit of the dimer at two anchoring points. Each monomer was treated as a rigid body and was docked based on the backbone PCS and backbone chemical shift perturbation data. Unlike NOE-based structural determination, this method only requires resonance assignments of the backbone (1)H(N)/(15)N signals and the PCS data obtained from several sets of two-dimensional (15)N-heteronuclear single quantum coherence spectra, thus facilitating rapid structure determination of the protein-protein complex.

DOTA-M8: An extremely rigid, high-affinity lanthanide chelating tag for PCS NMR spectroscopy

A new lanthanide chelating tag (M8) for paramagnetic labeling of biomolecules is presented, which is based on an eight-fold, stereoselectively methyl-substituted DOTA that can be covalently linked to the host molecule by a single disulfide bond. The steric overcrowding of the DOTA scaffold leads to an extremely rigid, kinetically and chemically inert lanthanide chelator. Its steric bulk restricts the motion of the tag relative to the host molecule. These properties result in very large pseudocontact shifts (>5 ppm) and residual dipolar couplings (>20 Hz) for Dy-M8 linked to ubiquitin, which are unprecedented for a small, single-point-attachment tag. Such large pseudocontact shifts should be well detectable even for larger proteins and distances beyond approximately 50 A. Due to its exceptionally high stability and lanthanide affinity M8 can be used under extreme chemical or physical conditions, such as those applied for protein denaturation, or when it is undesirable that buffer or protein react with excess lanthanide ions.

Complex formation between calmodulin and a peptide from the intracellular loop of the gap junction protein connexin43: Molecular conformation and energetics of binding

Gap junctions are formed by a family of transmembrane proteins, connexins. Connexin43 is a widely studied member of the family, being ubiquitously expressed in a variety of tissues and a target of a large number of disease mutations. The intracellular loop of connexin43 has been shown to include a calmodulin binding domain, but detailed 3-dimensional data on the structure of the complex are not available. In this study, we used a synthetic peptide from this domain to reveal the conformation of the calmodulin-peptide complex by small angle X-ray scattering. Upon peptide binding, calmodulin lost its dumbbell shape, adopting a more globular conformation. We also studied the energetics of the interaction using calorimetry and computational methods. All our data indicate that calmodulin binds to the peptide from cx43 in the classical 'collapsed' conformation.