Backbone and side chain dynamics of mutant calmodulin-peptide complexes

The mechanism of long-range coupling of allosteric sites in calcium-saturated calmodulin (CaM) has been explored by characterizing structural and dynamics effects of mutants of calmodulin in complex with a peptide corresponding to the smooth muscle myosin light chain kinase calmodulin-binding domain (smMLCKp). Four CaM mutants were examined: D95N and D58N, located in Ca2+-binding loops; and M124L and E84K, located in the target domain-binding site of CaM. Three of these mutants have altered allosteric coupling either between Ca2+-binding sites (D58N and D95N) or between the target- and Ca2+-binding sites (E84K). The structure and dynamics of the mutant calmodulins in complex with smMLCKp were characterized using solution NMR. Analysis of chemical shift perturbations was employed to detect largely structural perturbations. 15N and 2H relaxation was employed to detect perturbations of the dynamics of the backbone and methyl-bearing side chains of calmodulin. The least median squares method was found to be robust in the detection of perturbed sites. The main chain dynamics of calmodulin are found to be largely unresponsive to the mutations. Three mutants show significantly perturbed dynamics of methyl-bearing side chains. Despite the pseudosymmetric location of Ca2+-binding loop mutations D58N and D95N, the dynamic response of CaM is asymmetric, producing long-range perturbation in D58N and almost none in D95N. The mutations located at the target domain-binding site have quite different effects. For M124L, a local perturbation of the methyl dynamics is observed, while the E84K mutation produces a long-range propagation of dynamic perturbations along the target domain-binding site.     

Dynamics and entropy of a calmodulin-peptide complex studied by NMR and molecular dynamics

All-atom, explicit water molecular dynamics simulations of calcium-loaded calmodulin complexed with a peptide corresponding to the smooth muscle myosin light chain kinase target were carried out at 295 and 346 K. Amide and side chain methyl angular generalized order parameters were calculated and analyzed in the context of the protein's structure and dynamics. The agreement between amide order parameters measured by NMR and those from the simulations was found to be good, especially at the higher temperature, indicating both better convergence for the latter and excellent transferrability of the CHARMM parameters to the higher temperature. Subtle dynamical features such as helix fraying were reproduced. A large range of order parameters for the nine calmodulin methionines was observed in the NMR, and reproduced quite well in the simulations. The major determinant of the methionine order parameter was found to be the proximity to side chains of aromatic residues. An upper bound estimate of the difference in backbone entropy between loop and helical regions was extracted from the order parameters using a model of motion in an effective potential. Although loop regions are more flexible than helical regions, it was found that the entropy loss per residue upon folding was only approximately 20% less for loops than for helices. Pairwise correlated motions, which could significantly lower entropy estimates obtained from order parameter analysis alone, were found to be largely absent.     

The role of electrostatic interactions in calmodulin-peptide complex formation

The complex between calmodulin and the calmodulin-binding portion of smMLCKp has been studied. Electrostatic interactions have been anticipated to be important in this system where a strongly negative protein binds a peptide with high positive charge. Electrostatic interactions were probed by varying the pH in the range from 4 to 11 and by charge deletions in CaM and smMLCKp. The change in net charge of CaM from approximately -5 at pH 4.5 to -15 at pH 7.5 leaves the binding constant virtually unchanged. The affinity was also unaffected by mutations in CaM and charge substitutions in the peptide. The insensitivity of the binding constant to pH may seem surprising, but it is a consequence of the high charge on both protein and peptide. At low pH it is further attenuated by a charge regulation mechanism. That is, the protein releases a number of protons when binding the positively charged peptide. We speculate that the role of electrostatic interactions is to discriminate against unbound proteins rather than to increase the affinity for any particular target protein.     

Molecular dynamics simulations of a calmodulin-peptide complex in solution

We have performed an 4-ns MD simulation of calmodulin complexed with a target peptide in explicit water, under realistic conditions of constant temperature and pressure, in the presence of a physiological concentration of counterions and using Ewald summation to avoid truncation of long-range electrostatic forces. During the simulation the system tended to perform small fluctuations around a structure similar to, but somewhat looser than the starting crystal structure. The calmodulin-peptide complex was quite rigid and did not exhibit any large amplitude domain motions such as previously seen in apo- and calcium-bound calmodulin. We analyzed the calmodulin-peptide interactions by calculating buried surface areas, CHARMM interaction energies and continuum model interaction free energies. In the trajectory, the protein surface area buried by contact with the peptide is 1373 A(2) approximately evenly divided between the calmodulin N-terminal, C-terminal and central linker regions. A majority of this buried surface, 803 A(2), comes from nonpolar residues, in contrast to the protein as a whole, for which the surface is made up of mostly polar and charged groups. Our continuum calculations indicate that the largest favorable contribution to peptide binding comes from burial of molecular surface upon complex formation. Electrostatic contributions are favorable but smaller in the trajectory structures, and actually unfavorable for binding in the crystal structure. Since nonpolar groups make up most of buried surface of the protein, our calculations suggest that the hydrophobic effect is the main driving force for binding the helical peptide to calmodulin, consistent with thermodynamic analysis of experimental data. Besides the burial of nonpolar surface area, secondary contributions to peptide binding come from burial of polar surface and electrostatic interactions. In the nonpolar interactions a crucial role is played by the nine methionines of calmodulin. In the electrostatic interactions the negatively charged protein residues and positively charged peptide residues play a dominant role.     

Temperature dependence of the internal dynamics of a calmodulin-peptide complex

The temperature dependence of the fast internal dynamics of calcium-saturated calmodulin in complex with a peptide corresponding to the calmodulin-binding domain of the smooth muscle myosin light chain kinase is examined using 15N and 2H NMR relaxation methods. NMR relaxation studies of the complex were carried out at 13 temperatures that span 288-346 K. The dynamics of the backbone and over four dozen methyl-bearing side chains, distributed throughout the calmodulin molecule, were probed. The side chains show a much more variable and often considerably larger response to temperature than the backbone. A significant variation in the temperature dependence of the amplitude of motion of individual side chains is seen. The amplitude of motion of some side chains is essentially temperature-independent while many show a simple roughly linear temperature dependence. In a few cases, angular order increases with temperature, which is interpreted as arising from interactions with neighboring residues. In addition, a number of side chains display a nonlinear temperature dependence. The significance of these and other results is illuminated by several simple interpretative models. Importantly, analysis of these models indicates that changes in generalized order parameters can be robustly related to corresponding changes in residual entropy. A simple cluster model that incorporates features of cooperative or conditional motion reproduces many of the unusual features of the experimentally observed temperature dependence and illustrates that side chain interactions result in a dynamically changing environment that significantly influences the motion of internal side chains. This model also suggests that the intrinsic entropy of interacting clusters of side chains is only modestly reduced from that of independent side chain motion. Finally, estimates of protein heat capacity support the view that the major contribution to the heat capacity of protein solutions largely arises from local bond vibrations and solvent interactions and not from torsional oscillations of side chains.     

Characterization of the backbone and side chain dynamics of the CaM-CaMKIp complex reveals microscopic contributions to protein conformational entropy

Calmodulin is a central mediator of calcium-dependent signal transduction pathways and regulates the activity of a large number of diverse targets. Calcium-dependent interactions of calmodulin with regulated proteins are of generally high affinity but of quite variable thermodynamic origins. Here we investigate the influence of the binding of the calmodulin-binding domain of calmodulin kinase I on the fast internal dynamics of calcium-saturated calmodulin. NMR relaxation was used to probe motion on the backbone (viewed through the backbone amide NH group) and the side chains (viewed through methyl groups). The distribution of the amplitudes of side chain dynamics is trimodal. The microscopic details of side chain motion are compared with those of a thermodynamically and structurally similar complex of calmodulin with the calmodulin-binding domain of the smooth muscle myosin light chain kinase. While there are no significant differences in backbone dynamics and no net change in methyl-bearing side chain dynamics, a large redistribution of the amplitude of methyl dynamics is observed between the two complexes. The variation in dynamics was largely localized to the heterogeneously dynamic target-binding interface, suggesting that differential dynamics of the binding surface plays a functional role in the high-affinity binding interactions of calmodulin. These results begin to reveal a fundamental role for residual protein entropy in molecular recognition by calmodulin.     

Monomeric variants of IL-8: effects of side chain substitutions and solution conditions upon dimer formation.

IL-8 dimers have been observed in NMR and X-ray structures of the protein. We have engineered IL-8 monomers by mutations of residues throughout the dimer interface, which introduce hindrance determinants to dimerization. These IL-8 variants are shown by NMR to have wild-type monomer folding, but by ultracentrifugation to have a range of dimerization constants from microM to mM, as compared with a dimerization constant of about 10 microM for wild-type IL-8, under physiological salt and temperature conditions. The monomeric variants of IL-8 bind the erythrocyte chemokine receptor DARC, as well as the neutrophil IL-8 receptors CXCR1 and CXCR2 with affinities similar to that of wild-type IL-8. In addition, the monomeric variants were shown to have agonist activity, with similar potency to wild-type, in both Ca(2+)-flux assays on CXCR1 and CXCR2 transfected cells, and in chemotaxis assays on neutrophils. Thus, these variants confirm that monomeric IL-8 is functionally equivalent to wild-type in vitro assays. We have also investigated the effects of various solution conditions upon IL-8 dimer formation using analytical ultracentrifugation. At salt concentrations, temperatures, and pH conditions lower than physiological, the dimerization affinity of IL-8 is greatly enhanced. This suggests that, under some conditions, IL-8 dimer formation may occur at concentrations of IL-8 considerably lower than 10 microM, with consequences in vivo that are yet to be determined.     

Conformational changes of aminoacyl-tRNA and uncharged tRNA upon complex formation with polypeptide chain elongation factor Tu

The conformation change of Thermus thermophilus tRNA(1Ile) upon complex formation with T. thermophilus elongation factor Tu (EF-Tu) was studied by analysis of the circular dichroism (CD) bands at 315 nm (due to the 2-thioribothymidine residue in the T-loop) and at 295 nm (due to the core structure of tRNA). Formation of the ternary complex of isoleucyl-tRNA(1Ile) and EF-Tu.GTP increased the intensities of these CD bands, indicating stabilization of the association between the T-loop and the D-loop and also a significant conformation change of the core region. Upon complex formation of EF-Tu.GTP and uncharged tRNA, however, the conformation of the core region is not changed, while the association of the two loops is still stabilized. On the other hand, the binding with EF-Tu.GDP does not appreciably affect the conformation of isoleucyl-tRNA or uncharged tRNA. These indicate the importance of the gamma-phosphate group of GTP and the aminoacyl group in the formation of the active complex of aminoacyl-tRNA and EF-Tu.GTP.     

An increase in side chain entropy facilitates effector binding: NMR characterization of the side chain methyl group dynamics in Cdc42Hs

Cdc42Hs is a signal transduction protein that is involved in cytoskeletal growth and organization. We describe here the methyl side chain dynamics of three forms of (2)H,(13)C,(15)N-Cdc42Hs [GDP-bound (inactive), GMPPCP-bound (active), and GMPPCP/PBD46-bound (effector-bound)] from (13)C-(1)H NMR measurements of deuterium T(1) and T(1 rho) relaxation times. A wide variation in flexibility was observed throughout the protein, with methyl axis order parameters (S(2)(axis)) ranging from 0.2 to 0.4 (highly disordered) in regions near the PBD46 binding site to 0.8--1.0 (highly ordered) in some helices. The side chain dynamics of the GDP and GMPPCP forms are similar, with methyl groups on the PBD46 binding surface experiencing significantly greater mobility (lower S(2)(axis)) than those not on the binding surface. Binding of PBD46 results in a significant increase in the disorder and a corresponding increase in entropy for the majority of methyl groups. Many of the methyl groups that experience an increase in mobility are found in residues that are not part of the PBD46 binding interface. This entropy gain represents a favorable contribution to the overall entropy of effector binding and partially offsets unfavorable entropy losses such as those that occur in the backbone.     

Critical evaluation of methods to incorporate entropy loss upon binding in high-throughput docking

Proper accounting of the positional/orientational/conformational entropy loss associated with protein-ligand binding is important to obtain reliable predictions of binding affinity. Herein, we critically examine two simplified statistical mechanics-based approaches, namely a constant penalty per rotor method, and a more rigorous method, referred to here as the partition function-based scoring (PFS) method, to account for such entropy losses in high-throughput docking calculations. Our results on the estrogen receptor beta and dihydrofolate reductase proteins demonstrate that, while the constant penalty method over-penalizes molecules for their conformational flexibility, the PFS method behaves in a more "DeltaG-like" manner by penalizing different rotors differently depending on their residual entropy in the bound state. Furthermore, in contrast to no entropic penalty or the constant penalty approximation, the PFS method does not exhibit any bias towards either rigid or flexible molecules in the hit list. Preliminary enrichment studies using a lead-like random molecular database suggest that an accurate representation of the "true" energy landscape of the protein-ligand complex is critical for reliable predictions of relative binding affinities by the PFS method.     

Direct estimation of entropy loss due to reduced translational and rotational motions upon molecular binding

The entropic cost due to the loss of translational and rotational (T-R) degree of freedom upon binding has been well recognized for several decades. Tightly bound ligands have higher entropic costs than loosely bound ligands. Quantifying the ligand's residual T-R motions after binding, however, is not an easy task. We describe an approach that uses a reduced Hessian matrix to estimate the contributions due to translational and rotational degrees of freedom to entropy change upon molecular binding. The calculations use a harmonic model for the bound state but only include the T-R degrees of freedom. This approximation significantly speeds up entropy calculations because only 6 x 6 matrices need to be treated, which makes it easier to be used in computer-aided drug design for studying many ligands. The methodological connection with other methods is discussed as well. We tested this approximation by applying it to study the binding of ATP, peptide inhibitor (PKI), and several bound water molecules to protein kinase A (PKA). These ligands span a wide range in size. The model gave reasonable estimates of the residual T-R entropy of bound ligands or water molecules. The residual T-R entropy demonstrated a wide range of values, e.g., 4 to 16 cal/K.mol for the bound water molecules of PKA.     

Side chain to side chain cyclization of an enkephalin analog results in loss of opioid receptor selectivity

The cyclic enkephalin analog H-Tyr-D-Lys-Gly-Phe-Glu-NH2 (I) and the structurally related open chain analogs H-Tyr-D-Nle-Gly-Phe-Gln-NH2 (II) and H-Tyr-D-Lys(For)-Gly-Phe-Abu-NH2 (III) were tested in mu and delta opioid receptor-representative binding assays and bioassays. Whereas both linear analogs showed a pronounced preference for mu receptors over delta receptors, the conformationally restricted cyclic peptide I was found to be unselective. This finding represents the first reported example of a peptide cyclization resulting in a loss of receptor selectivity. From this and earlier studies, it was concluded that the receptor selectivity of cyclized peptide analogs relative to that of their linear correlates may depend on the size and relative rigidity of their ring structures.     

Processive ubiquitin chain formation by the anaphase-promoting complex

Progression through mitosis requires the sequential ubiquitination of cell cycle regulators by the anaphase-promoting complex, resulting in their proteasomal degradation. Although several mechanisms contribute to APC/C regulation during mitosis, the APC/C is able to discriminate between its many substrates by exploiting differences in the processivity of ubiquitin chain assembly. Here, we discuss how the APC/C achieves processive ubiquitin chain formation to trigger the sequential degradation of cell cycle regulators during mitosis.     

Helix formation and capping energetics of arginine analogs with varying side chain length

Arginine (Arg) has been used for recognizing negatively charged biological molecules, cell penetration, and oligosaccharide mass signal enhancement. The versatility of Arg has inspired the need to develop Arg analogs and to research the structural effects of incorporating Arg analogs. Accordingly, we investigated the effect of Arg side chain length on helix formation by studying 12 Ala-based peptides containing the Arg analogs (S)-2-amino-6-guanidino-hexanoic acid (Agh), (S)-2-amino-4-guanidinobutyric acid (Agb), and (S)-2-amino-3-guanidinopropionic acid (Agp). Solid phase guanidinylation with orthogonal protection strategies was necessary to synthesize Agb- and Agp-containing peptides using Fmoc-based chemistry. The fraction helix for the peptides was determined by circular dichroism spectroscopy, and used to derive the statistical mechanical parameters and energetics for N-capping, C-capping, and helix propagation (propensity). All four Arg analogs were unfavorable for N-capping. The C-cap parameter followed the trend AgpAgh, highlighting the uniqueness of the Arg side chain length in helix formation. Molecular mechanics calculations and a survey on protein structures were consistent with the experimental results. Furthermore, calculations and survey both showed that the g- conformation for the χ1 dihedral was present for the first two residues at the N-terminus of helices, but not favored in the center or C-terminus of helices due to sterics. These results should serve as the foundation for developing Arg-related bioactive compounds and technologies.     

Fragmentation of a 70000-dalton calpastatin molecule upon its complex formation with calpain

Homogenously purified porcine calpain I (Mr 112000), a low-Ca2+-requiring form of Ca2+-dependent cysteine proteinase [EC 3.4.22.17], was coupled to Sepharose 4B gel as an active form. It was used as a ligand to calpastatin (Mr 70000), calpain-specific inhibitor protein, for an affinity chromatography. Only in the presence of Ca2+, calpastatin bound to calpain-Sepharose, but the interaction resulted in rather extensive fragmentation of a calpastatin molecule into several peptides of Mr 14000 to 70000, which still retain inhibitory activities against calpain. Fragmentation was demonstrated both by polyacrylamide gel electrophoresis in the presence of sodium dodecyl sulfate and by high-performance liquid chromatography in the presence of 6 M guanidine-HCl.     

A model for the calmodulin-peptide complex based on the troponin C crystal packing and its similarity to the NMR structure of the calmodulin-myosin light chain kinase peptide complex.

In the crystal structure of troponin C, the holo C-domain is bound in a head-to-tail fashion to the A-helix of the apo N-domain of a symmetry-related molecule. Using this interaction, we have proposed a model for the calmodulin-peptide complex. We find that the interaction of the C-domain with the A-helix is similar to that observed in the NMR structure of the calmodulin-myosin light chain kinase (MLCK) peptide complex. This similarity in binding has enabled us to make a precise sequence alignment of the target peptides in the calmodulin-binding cleft and to rationalize the amino acid sequence-dependent binding strengths of various peptides. Our model differs from that proposed by Strynadka and James (Proteins Struct. Funct. Genet. 7, 234-248, 1990) in that the peptides are rotated by 100 degrees in the calmodulin binding cleft.     

Change in entropy of the free polypeptide chain during formation of hydrogen bonds]

Formation probabilities of different hydrogen bonds between carbonyl oxygen and amide hydrogen were determined by Monte Carlo simulations using a computer model in the space of sterically allowable conformations of alanine and glycine oligopeptides, and the corresponding entropy losses for the peptide backbone, T delta S, were calculated. The model was studied at different criteria of steric interactions. Comparison with the data of other authors showed the values of T delta S to be mainly determined by overall extent and type of the state space and to be only slightly dependent on its energy profile. Both short-range and long-range steric interactions were shown to prevent hydrogen bonding, especially in alanine peptides. In the model studied, the initiation of alpha(R)-helices is associated with T delta S = 8-10 kT, and prior formation of a 3/10-turn or one three-center H-bond does not appreciably decrease this entropy barrier. Elongation of the alpha(R)-helix by one residue leads to T delta S = 3.0-3.7 kT, the helices begin to stabilize after at least three sequential H-bonds are formed. The difference in the probability of insertion of Ala and Gly into the helix is lower than it follows from comparison of their mobility. The results could be explained assuming that factors different from helical H-bonds take part in the stabilization of the helices. One may suppose upon modeling of folding that even three sequential H-bonds are unable to fix the structure of a flexible peptide loop, while the elongation of alpha(R)-helices in the supersecondary helix-loop-helix structure is favorable as long as the loop conformation remains nearly optimal.     

Screening of metal complex–amino acid side chain interactions by potentiometric titration

Reversible coordination of amino acid side chains to metal complexes is widely used in protein purification (IMAC technique), but available data on affinity and selectivity of such binding processes are limited. We use potentiometric titration of a series of metal complexes with vacant coordination sites in the presence of molecules resembling amino acid side chain functionalities to screen for new affinities. The investigation confirms documented affinities of imidazole to nickel(II) and copper(II) IDA and NTA complexes, and discovers a hitherto unknown binding of zinc(II)- and cadmium(II) cyclen complexes to imidazole.     

Vitamin K2 and geranylgeraniol, its side chain component, inhibited osteoclast formation in a different manner

We comparatively examined the mechanism by which vitamin K(2) (Menatetrenone, MK4) and its side chain component, geranylgeraniol (GGO), inhibited osteoclast formation in the co-culture system of stromal cells with spleen cells. Both MK4 and GGO inhibited osteoclast formation induced by 1alpha,25-dihydroxyvitamin D3 (1,25(OH)(2)D(3)). MK4, but not GGO, inhibited cyclooxygenase-2 (COX-2) expression and prostaglandin E(2) (PGE(2)) production in the co-culture system. To elucidate the precise mechanism of the inhibitory effect of GGO on osteoclast formation, the co-cultured cells were stimulated with PGE(2). GGO, but not MK4, inhibited osteoclast formation via suppression of the receptor activator of NF-kappaB ligand (RANKL) expression. Moreover, GGO abolished the disruption of osteoclastic actin rings induced by nitrogen-containing bisphosphonate (N-BP), whereas MK4 did not affect it at all. These data suggest that MK4 inhibited osteoclast formation independently of GGO, and that MK4, but not GGO, has no competitive action on the anti-osteoporotic effect of N-BP.