Bending of the calmodulin central helix: a theoretical study.

The crystal structure of calcium-calmodulin (CaM) reveals a protein with a typical dumbbell structure. Various spectroscopic studies have suggested that the central linker region of CaM, which is alpha-helical in the crystal structure, is flexible in solution. In particular, NMR studies have indicated the presence of a flexible backbone between residues Lys 77 and Asp 80. This flexibility is related directly to the function of the protein because it enables the N- and C-terminal domains of the protein to move toward each other and bind to the CaM-binding domain of a target protein. We have investigated the flexibility of the CaM central helix by a variety of computational techniques: molecular dynamics (MD) simulations, normal mode analysis (NMA), and essential dynamics (ED) analysis. Our MD results reproduce the experimentally determined location of the bend in a simulation of only the CaM central helix, indicating that the bending point is an intrinsic property of the alpha-helix, for which the remainder of the protein is not important. Interestingly, the modes found by the ED analysis of the MD trajectory are very similar to the lowest frequency modes from the NM analysis and to modes found by an ED analysis of different structures in a set of NMR structures. Electrostatic interactions involving residues Arg 74 and Asp 80 seem to be important for these bending motions and unfolding, which is in line with pH-dependent NMR and CD studies.     

Regulation of mammalian nitric oxide synthases by electrostatic interactions in the linker region of calmodulin

Calmodulin (CaM), the ubiquitous Ca(2+)-sensing protein, consists of two globular domains separated by a flexible central linker that properly orients CaM's globular domains to bind and regulate various intracellular proteins, including the nitric oxide synthase (NOS) enzymes. In the present study we determined that the charge and length of the central linker of CaM has an effect on the binding and activation of the NOS isozymes by using a variety of charge CaM mutants (T79D, S81D, T79D/S81D, S101D and E84R/E87K) and CaM mutants with residues removed (Delta84, Delta83-84, and Delta81-84). Our kinetic and spectropolarimetry results demonstrate that the NOS enzymes are not adversely affected by the CaM mutants with the exceptions of S101D, E84R/E87K and the deletion of residue 84. Electrostatic interactions in the central linker between residues 82-87 in combination with hydrophobic interactions in the globular domains of CaM are important for its tight association to inducible NOS.     

Secondary structure and side-chain 1H and 13C resonance assignments of calmodulin in solution by heteronuclear multidimensional NMR spectroscopy

Heteronuclear 2D and 3D NMR experiments were carried out on recombinant Drosophila calmodulin (CaM), a protein of 148 residues and with molecular mass of 16.7 kDa, that is uniformly labeled with 15N and 13C to a level of greater than 95%. Nearly complete 1H and 13C side-chain assignments for all amino acid residues are obtained by using the 3D HCCH-COSY and HCCH-TOCSY experiments that rely on large heteronuclear one-bond scalar couplings to transfer magnetization and establish through-bond connectivities. The secondary structure of this protein in solution has been elucidated by a qualitative interpretation of nuclear Overhauser effects, hydrogen exchange data, and 3JHNH alpha coupling constants. A clear correlation between the 13C alpha chemical shift and secondary structure is found. The secondary structure in the two globular domains of Drosophila CaM in solution is essentially identical with that of the X-ray crystal structure of mammalian CaM [Babu, Y., Bugg, C. E., & Cook, W.J. (1988) J. Mol. Biol. 204, 191-204], which consists of two pairs of a "helix-loop-helix" motif in each globular domain. The existence of a short antiparallel beta-sheet between the two loops in each domain has been confirmed. The eight alpha-helix segments identified from the NMR data are located at Glu-6 to Phe-19, Thr-29 to Ser-38, Glu-45 to Glu-54, Phe-65 to Lys-77, Glu-82 to Asp-93, Ala-102 to Asn-111, Asp-118 to Glu-127, and Tyr-138 to Thr-146. Although the crystal structure has a long "central helix" from Phe-65 to Phe-92 that connects the two globular domains, NMR data indicate that residues Asp-78 to Ser-81 of this central helix adopt a nonhelical conformation with considerable flexibility.     

NMR derived solution structure of an EF-hand calcium-binding protein from Entamoeba Histolytica.

We present the three-dimensional (3D) solution structure of a calcium-binding protein from Entamoeba histolytica (EhCaBP), an etiologic agent of amoebiasis affecting millions worldwide. EhCaBP is a 14.7 kDa (134 residues) monomeric protein thought to play a role in the pathogenesis of amoebiasis. The 3D structure of Ca(2+)-bound EhCaBP has been derived using multidimensional nuclear magnetic resonance (NMR) spectroscopic techniques. The study reveals the presence of two globular domains connected by a flexible linker region spanning 8 amino acid residues. Each domain consists of a pair of helix-loop-helix motifs similar to the canonical EF-hand motif of calcium-binding proteins. EhCaBP binds to four Ca(2+) with high affinity (two in each domain), and it is structurally related to calmodulin (CaM) and troponin C (TnC) despite its low sequence homology ( approximately 29%) with these proteins. NMR-derived structures of EhCaBP converge within each domain with low RMSDs and angular order-parameters for backbone torsion angles close to 1.0. However, the presence of a highly flexible central linker region results in an ill-defined orientation of the two domains relative to one other. These findings are supported by backbone (15)N relaxation rate measurements and deuterium exchange studies, which reveal low structural order parameters for residues in the central linker region. Earlier, biochemical studies showed that EhCaBP is involved in a novel signal transduction mechanism, distinct from CaM. A possible reason for such a functional diversity is revealed by a detailed comparison of the 3D structure of EhCaBP with that of CaM and TnC. The studies indicate a more open C-terminal domain for EhCaBP with larger water exposed total hydrophobic surface area as compared to CaM and TnC. Further dissimilarities between the structures include the presence of two Gly residues (G63 and G67) in the central linker region of EhCaBP, which seem to impart it a greater flexibility compared to CaM and TnC and also play crucial role in its biological function. Thus, unlike in CaM and TnC, wherein the length and/or composition of the central linker have been found to be crucial for their function, in EhCaBP, both flexibility as well as amino acid composition is required for the function of the protein.     

Solution structure of the two-iron rubredoxin of Pseudomonas oleovorans determined by NMR spectroscopy and solution X-ray scattering and interactions with rubredoxin reductase

Here we provide insights into the molecular structure of the two-iron 19-kDa rubredoxin (AlkG) of Pseudomonas oleovorans using solution-state nuclear magnetic resonance (NMR) and small-angle X-ray scattering studies. Sequence alignment and biochemical studies have suggested that AlkG comprises two rubredoxin folds connected by a linker region of approximately 70 amino acid residues. The C-terminal domain (C-Rb) of this unusual rubredoxin, together with approximately 35 amino acid residues of the predicted linker region, was expressed in Escherichia coli, purified in the one-iron form and the structure of the cadmium-substituted form determined at high-resolution by NMR spectroscopy. The structure shows that the C-Rb domain is similar in fold to the conventional one-iron rubredoxins from other organisms, whereas the linker region does not have any discernible structure. This tandem "flexible-folded" structure of the polypeptide chain derived for the C-Rb protein was confirmed using solution X-ray scattering methods. X-ray scattering studies of AlkG indicated that the 70-amino acid residue linker forms a structured, yet mobile, polypeptide segment connecting the globular N- and C-terminal domains. The X-ray scattering studies also showed that the N-terminal domain (N-Rb) has a molecular conformation similar to that of C-Rb. The restored molecular shape indicates that the folded N-Rb and C-Rb domains of AlkG are noticeably separated, suggesting some domain movement on complex formation with rubredoxin reductase to allow interdomain electron transfer between the metal centers in AlkG. This study demonstrates the advantage of combining X-ray scattering and NMR methods in structural studies of dynamic, multidomain proteins that are not suited to crystallographic analysis. The study forms a structural foundation for functional studies of the interaction and electron-transfer reactions of AlkG with rubredoxin reductase, also reported herein.     

Molecular dynamics simulation of proteins under high pressure: Structure,function and thermodynamics

BackgroundMolecular dynamics (MD) simulation is well-recognized as a powerful tool to investigate protein structure, function, and thermodynamics. MD simulation is also used to investigate high pressure effects on proteins. For conducting better MD simulation under high pressure, the main issues to be addressed are: (i) protein force fields and water models were originally developed to reproduce experimental properties obtained at ambient pressure; and (ii) the timescale to observe the pressure effect is often much longer than that of conventional MD simulations.Scope of reviewFirst, we describe recent developments in MD simulation methodologies for studying the high-pressure structure and dynamics of protein molecules. These developments include force fields for proteins and water molecules, and enhanced simulation techniques. Then, we summarize recent studies of MD simulations of proteins in water under high pressure.Major conclusionsRecent MD simulations of proteins in solution under pressure have reproduced various phenomena identified by experiments using high pressure, such as hydration, water penetration, conformational change, helix stabilization, and molecular stiffening.General significanceMD simulations demonstrate differences in the properties of proteins and water molecules between ambient and high-pressure conditions. Comparing the results obtained by MD calculations with those obtained experimentally could reveal the mechanism by which biological molecular machines work well in collaboration with water molecules.     

Conformational study of cyclolinopeptide A. A distance geometry and molecular dynamics approach

The conformation of cyclolinopeptide A, c(Pro-Pro-Phe-Phe-Leu-Ile-Ile-Leu-Val), a naturally occurring peptide with remarkable cytoprotective activity, has been investigated by means of distance geometry calculations and molecular dynamics simulations. The starting points for all the calculations were an X-ray structure and other structures obtained from distance geometry calculations based on NMR data. Restrained and unrestrained molecular dynamics simulations are reported in vacuo and in CCl4. Structural and dynamic properties are investigated and compared with those experimentally determined. The conformation obtained from the MD simulations which best reproduces the NMR parameters is at the same time one of the most stable ones and is also fairly similar to the crystal structure. An explanation for the occurrence of multiple conformations in solution at room temperature is given.     

Solution X-ray scattering reveals a novel structure of calmodulin complexed with a binding domain peptide from the HIV-1 matrix protein p17

The solution structures of complexes between calcium-saturated calmodulin (Ca (2+)/CaM) and a CaM-binding domain of the HIV-1 matrix protein p17 have been determined by small-angle X-ray scattering with use of synchrotron radiation as an intense and stable X-ray source. We used three synthetic peptides of residues 11-28, 26-47, and 11-47 of p17 to demonstrate the diversity of CaM-binding conformation. Ca (2+)/CaM complexed with residues 11-28 of p17 adopts a dumbbell-like structure at a molar ratio of 1:2, suggesting that the two peptides bind each lobe of CaM, respectively. Ca (2+)/CaM complexed with residues 26-47 of p17 at a molar ratio of 1:1 adopts a globular structure similar to the NMR structure of Ca (2+)/CaM bound to M13, which adopted a compact globular structure. In contrast to these complexes, Ca (2+)/CaM binds directly with both CaM-binding sites of residues 11-47 of p17 at a molar ratio of 1:1, which induces a novel structure different from known structures previously reported between Ca (2+)/CaM and peptide. A tertiary structural model of the novel structure was constructed using the biopolymer module of Insight II 2000 on the basis of the scattering data. The two domains of CaM remain essentially unchanged upon complexation. The hinge motions, however, occur in a highly flexible linker of CaM, in which the electrostatic residues 74Arg, 78Asp, and 82Glu interact with N-terminal electrostatic residues of the peptide (residues 12Glu, 15Arg, and 18Lys). The acidic residues in the N-terminal domain of CaM interact with basic residues in a central part of the peptide, thereby enabling the central part to change the conformations, while an acidic residue in the C-terminal domain interacts with two basic residues in the two helical sites of the peptide. The overall structure of the complex adopts an extended structure with the radius of gyration of 20.5 A and the interdomain distance of 34.2 A. Thus, the complex is principally stabilized by electrostatic interactions. The hydrophobic patches of Ca (2+)/CaM are not responsible for the binding with the hydrophobic residues in the peptide, suggesting that CaM plays a role to sequester the myristic acid moiety of p17.     

Molecular dynamics simulation of hen egg white lysozyme: a test of the GROMOS96 force field against nuclear magnetic resonance data

Biomolecular force fields for use in molecular dynamics (MD) simulations of proteins, DNA, or membranes are generally parametrized against ab initio quantum-chemical and experimental data for small molecules. The application of a force field in a simulation of a biomolecular system, such as a protein in solution, may then serve as a test of the quality and transferability of the force field. Here, we compare various properties obtained from two MD simulations of the protein hen egg white lysozyme (HEWL) in aqueous solution using the latest version, GROMOS96, of the GROMOS force field and an earlier version, GROMOS87+, with data derived from nuclear magnetic resonance (NMR) experiments: NOE atom-atom distance bounds, (3)J(HNalpha)-coupling constants, and backbone and side-chain order parameters. The convergence of these quantities over a 2-ns period is considered, and converged values are compared to experimental ones. The GROMOS96 simulation shows better agreement with the NMR data and also with the X-ray crystal structure of HEWL than the GROMOS87+ simulation, which was based on an earlier version of the GROMOS force field.     

Interdomain orientation of cardiac Troponin C characterized by paramagnetic relaxation enhancement NMR reveals a compact state

Cardiac troponin C (cTnC) is the calcium binding subunit of the troponin complex that triggers the thin filament response to calcium influx into the sarcomere. cTnC consists of two globular EF-hand domains (termed the N- and C-domains) connected by a flexible linker. While the conformation of each domain of cTnC has been thoroughly characterized through NMR studies involving either the isolated N-domain (N-cTnC) or C-domain (C-cTnC), little attention has been paid to the range of interdomain orientations possible in full-length cTnC that arises as a consequence of the flexibility of the domain linker. Flexibility in the domain linker of cTnC is essential for effective regulatory function of troponin. We have therefore utilized paramagnetic relaxation enhancement (PRE) NMR to assess the interdomain orientation of cTnC. Ensemble fitting of our interdomain PRE measurements reveals that isolated cTnC has considerable interdomain flexibility and preferentially adopts a bent conformation in solution, with a defined range of relative domain orientations.     

Temperature‐accelerated molecular dynamics gives insights into globular conformations sampled in the free state of the AC catalytic domain

The catalytic domain of the adenyl cyclase (AC) toxin from Bordetella pertussis is activated by interaction with calmodulin (CaM), resulting in cAMP overproduction in the infected cell. In the X‐ray crystallographic structure of the complex between AC and the C terminal lobe of CaM, the toxin displays a markedly elongated shape. As for the structure of the isolated protein, experimental results support the hypothesis that more globular conformations are sampled, but information at atomic resolution is still lacking. Here, we use temperature‐accelerated molecular dynamics (TAMD) simulations to generate putative all‐atom models of globular conformations sampled by CaM‐free AC. As collective variables, we use centers of mass coordinates of groups of residues selected from the analysis of standard molecular dynamics (MD) simulations. Results show that TAMD allows extended conformational sampling and generates AC conformations that are more globular than in the complexed state. These structures are then refined via energy minimization and further unrestrained MD simulations to optimize inter‐domain packing interactions, thus resulting in the identification of a set of hydrogen bonds present in the globular conformations. Proteins 2014; 82:2483–2496. © 2014 Wiley Periodicals, Inc.     

A structural comparison of ‘real’ and ‘model’ calmodulin clarified allosteric interactions regulating domain motion

Calmodulin (CaM) is a multifunctional calcium-binding protein, which regulates various biochemical processes. CaM acts via structural changes and complex forming with its target enzymes. CaM has two globular domains (N-lobe and C-lobe) connected by a long linker region. Upon calcium binding, the N-lobe and C-lobe undergo local conformational changes, after that, entire CaM wraps the target enzyme through a large conformational change. However, the regulation mechanism, such as allosteric interactions regulating the conformational changes, is still unclear. In order to clarify the allosteric interactions, in this study, experimentally obtained ‘real’ structures are compared to ‘model’ structures lacking the allosteric interactions. As the allosteric interactions would be absent in calcium-free CaM (apo-CaM), allostery-eliminated calcium-bound CaM (holo-CaM) models were constructed by combining the apo-CaM’s linker and the holo-CaM’s N- and C-lobe. Before the comparison, the ‘real’ and ‘model’ structures were clustered and cluster–cluster relationship was determined by a principal component analysis. The structures were compared based on the relationship, then, a distance map and a contact probability analysis clarified that the inter-domain motion is regulated by several groups of inter-domain contacting residue pairs. The analyses suggested that these residues cause inter-domain translation and rotation, and as a consequence, the motion encourage structural diversity. The resultant diversity would contribute to the functional versatility of CaM.     

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

We have determined the solution structure of calmodulin (CaM) from yeast (Saccharomyces cerevisiae) (yCaM) in the apo state by using NMR spectroscopy. yCaM is 60% identical in its amino acid sequence with other CaMs, and exhibits its unique biological features. yCaM consists of two similar globular domains (N- and C-domain) containing three Ca(2+)-binding motifs, EF-hands, in accordance with the observed 3 mol of Ca(2+) binding. In the solution structure of yCaM, the conformation of the N-domain conforms well to the one of the expressed N-terminal half-domains of yCaM [Ishida, H., et al. (2000) Biochemistry 39, 13660-13668]. The conformation of the C-domain basically consists of a pair of helix-loop-helix motifs, though a segment corresponding to the forth Ca(2+)-binding site of CaM deviates in its primary structure from a typical EF-hand motif and loses the ability to bind Ca(2+). Thus, the resulting conformation of each domain is essentially identical to the corresponding domain of CaM in the apo state. A flexible linker connects the two domains as observed for CaM. Any evidence for the previously reported interdomain interaction in yCaM was not observed in the solution structure of the apo state. Hence, the interdomain interaction possibly occurs in the course of Ca(2+) binding and generates a cooperative Ca(2+) binding among all three sites. Preliminary studies on a mutant protein of yCaM, E104Q, revealed that the Ca(2+)-bound N-domain interacts with the apo C-domain and induces a large conformational change in the C-domain.     

Conformational dynamics of yeast calmodulin in the Ca(2+)-bound state probed using NMR relaxation dispersion

Most calmodulin (CaM) in apo and Ca(2+)-bound states show a dumb-bell-like structure, involving the N- and C-terminal domains, connected with a flexible linker. However, Ca(2+)-bound yeast calmodulin (yCaM) takes on a unique globular structure; the target-binding site of this protein is autoinhibited. We applied NMR relaxation dispersion experiments to yCaM in the Ca(2+)-bound state. The amide (15)N and (1)H(N) relaxation dispersion profiles indicated the presence of conformational dynamics for specific residues at the interface between the N- and C-terminal domains. We conclude that these conformational dynamics were derived from the mobility of the C-terminal domain.     

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

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

A link between sequence conservation and domain motion within the AAA+ family

The AAA+ family of proteins play fundamental roles in all three kingdoms of life. It is thought that they act as molecular chaperones in aiding the assembly or disassembly of proteins or protein complexes. Recent structural studies on a number of AAA+ family proteins have revealed that they share similar structural elements. These structures provide a possible link between nucleotide binding/hydrolysis and the conformational changes which are then amplified to generate mechanical forces for their specific functions. However, from these individual studies it is far from clear whether AAA+ proteins in general share properties in terms of nucleotide induced conformational changes. In this study, we analyze sequence conservation within the AAA+ family and identify two subfamilies, each with a distinct conserved linker sequence that may transfer conformational changes upon ATP binding/release to movements between subdomains and attached domains. To investigate the relation of these linker sequences to conformational changes, molecular dynamics (MD) simulations on X-ray structures of AAA+ proteins from each subfamily have been performed. These simulations show differences in both the N-linker peptide, subdomain motion, and cooperativity between elements of quaternary structure. Extrapolation of subdomain movements from one MD simulation enables us to produce a structure in close agreement with cryo-EM experiments.     

Detecting protein kinase recognition modes of calmodulin by residual dipolar couplings in solution NMR

Calmodulin-regulated serine/threonine kinases (CaM kinases) play crucial roles in Ca2+-dependent signaling transduction pathways in eukaryotes. Despite having a similar overall molecular architecture of catalytic and regulatory domains, CaM kinases employ different binding modes for Ca2+/CaM recruitment which is required for their activation. Here we present a residual dipolar coupling (RDC)-based NMR approach to characterizing the molecular recognition of CaM with five different CaM kinases. Our analyses indicate that CaM kinase I and likely IV use the same CaM binding mode as myosin light chain kinase (1-14 motif), distinct from those of CaM kinase II (1-10 motif) and CaM kinase kinase (1-16- motif). This NMR approach provides an efficient experimental guide for homology modeling and structural characterization of CaM-target complexes.     

The solution structure of cytochrome c 6 from the green alga Monoraphidium braunii

 The three-dimensional structure in solution of the reduced form of cytochrome c ₆ from the green alga Monoraphidium braunii has been solved through NMR data. Cytochrome c ₆ acts as a small mono-heme electron carrier protein between the two membrane-embedded complexes cytochrome f and photosystem I. The structure was determined using 1278 relevant interproton NOEs out of 1776 assigned NOEs with distance geometry (DG) calculations which included 36 stereospecific assignments and 20 experimentally found angle constraints. The family of structures obtained from the DG calculations was subjected to energy minimization and molecular dynamics simulation using previously defined force field parameters for the heme and its ligands. In all stages of the calculations, the solution structure is well defined and similar to the now available X-ray structure. Received: 18 January 1996 / Accepted: 19 April 1996     

Conformational transitions of flanking purines in HIV-1 RNA dimerization initiation site kissing complexes studied by CHARMM explicit solvent molecular dynamics

Dimerization of HIV-1 genomic RNA is initiated by kissing loop interactions at the Dimerization Initiation Site (DIS). Dynamics of purines that flank the 5' ends of the loop-loop helix in HIV-1 DIS kissing complex were explored using explicit solvent molecular dynamics (MD) simulations with the CHARMM force field. Multiple MD simulations (200 ns in total) of X-ray structures for HIV-1 DIS Subtypes A, B, and F revealed conformational variability of flanking purines. In particular, the flanking purines, which in the starting X-ray structures are bulged-out and stack in pairs, formed a consecutive stack of four bulged-out adenines at the beginning of several simulations. This conformation is seen in the crystal structure of DIS Subtype F with no interference from crystal packing, and was frequently reported in our preceding MD studies performed with the AMBER force field. However, as CHARMM simulations progressed, the four continuously stacked adenines showed conformational transitions from the bulged-out into the bulged-in geometries. Although such an arrangement has not been seen in any X-ray structure, it has been suggested by a recent NMR investigation. In CHARMM simulations, in the longer time scale, the flanking purines display the tendency to move to bulged-in conformations. This is in contrast with the AMBER simulations, which indicate a modest prevalence for bulged-out flanking base positions in line with the X-ray data. The simulations also suggest that the intermolecular stacking between purines from the opposite hairpins can additionally stabilize the kissing complex.