Structure and dynamics of calmodulin in solution.

To characterize the dynamic behavior of calmodulin in solution, we have carried out molecular dynamics (MD) simulations of the Ca2+-loaded structure. The crystal structure of calmodulin was placed in a solvent sphere of radius 44 A, and 6 Cl- and 22 Na+ ions were included to neutralize the system and to model a 150 mM salt concentration. The total number of atoms was 32,867. During the 3-ns simulation, the structure exhibits large conformational changes on the nanosecond time scale. The central alpha-helix, which has been shown to unwind locally upon binding of calmodulin to target proteins, bends and unwinds near residue Arg74. We interpret this result as a preparative step in the more extensive structural transition observed in the "flexible linker" region 74-82 of the central helix upon complex formation. The major structural change is a reorientation of the two Ca2+-binding domains with respect to each other and a rearrangement of alpha-helices in the N-terminus domain that makes the hydrophobic target peptide binding site more accessible. This structural rearrangement brings the domains to a more favorable position for target binding, poised to achieve the orientation observed in the complex of calmodulin with myosin light-chain kinase. Analysis of solvent structure reveals an inhomogeneity in the mobility of water in the vicinity of the protein, which is attributable to the hydrophobic effect exerted by calmodulin's binding sites for target peptides.     

Molecular dynamics simulations of calcium-free calmodulin in solution

A 4-ns molecular dynamics simulation of calcium-free calmodulin in solution has been performed, using Ewald summation to treat electrostatic interactions. Our simulation results were mostly consistent with solution experimental studies, including NMR, fluorescence and x-ray scattering. The secondary structures within the N- and C-terminal domains were conserved in the simulation, with trajectory structures similar to the NMR-derived model structure 1CFD. However, the relative orientations of the domains, for which there are no NMR restraints, differed in details between the simulation and the 1CFD model. The most interesting information provided by the simulations is that the dynamics of calcium-free calmodulin in solution is dominated by slow rigid body reorientations of the domains. The interdomain distance fluctuated between 29 and 39 A, and interdomain orientation angle, defined as the pseudo-dihedral formed by the four calcium binding sites, varied between -2 degrees and 108 degrees. Similarly, the domain linker region also exhibited significant fluctuations, with its length varying in the 34-45 A range and its bend angle in the 10-100 degrees range. The simulations are in accord with fluorescence results suggesting that calcium-free calmodulin is more compact and more flexible than the calcium activated form. Surprisingly, quite similar solvent accessibilities of the hydrophobic patches were seen in the calcium-free trajectory described in this work and previously generated calcium-loaded calmodulin simulations. Thus, our simulations suggest a reexamination of the standard model of the structural change of calmodulin upon calcium binding, involving exposure of the hydrophobic patches to solvent.     

Spin labeling studies of wheat germ calmodulin in solution.

Electron paramagnetic resonance was used to investigate the physical state of plant calmodulin in solution. Wheat germ calmodulin contains a single cysteine residue (Cys-27) on the first of four calcium binding loops. In this study the nitroxide spin label 2,2,6,6-tetramethyl-4-maleimidopiperidine-1-oxyl (MAL-6) was covalently attached to Cys-27 to produce a Ca(2+)-sensitive, biologically-active, labeled protein. The rotational correlation time of the spin label, a measure of its rotational mobility and reflective of the physical state of this region of the protein, was calculated under various conditions. Relative to control, changes in the physical state of the protein reflected by increased motion of the spin label were observed at high pH, low ionic strength and upon addition of Ca2+. These results extend knowledge of the structure of the protein, previously known from solid state and biochemical studies, to calmodulin in solution.     

Structure and backbone dynamics of Apo- and holo-cellular retinol-binding protein in solution

Retinoid-binding proteins play an important role in regulating transport, storage, and metabolism of vitamin A and its derivatives. The solution structure and backbone dynamics of rat cellular retinol-binding protein type I (CRBP) in the apo- and holo-form have been determined and compared using multidimensional high resolution NMR spectroscopy. The global fold of the protein is consistent with the common motif described for members of the intracellular lipid-binding protein family. The most relevant difference between the NMR structure ensembles of apo- and holoCRBP is the higher backbone disorder, in the ligand-free form, of some segments that frame the putative entrance to the ligand-binding site. These comprise alpha-helix II, the subsequent linker to beta-strand B, the hairpin turn between beta-strands C and D, and the betaE-betaF turn. The internal backbone dynamics, obtained from 15N relaxation data (T1, T2, and heteronuclear nuclear Overhauser effect) at two different fields, indicate several regions with significantly higher backbone mobility in the apoprotein, including the betaC-betaD and betaE-betaF turns. Although apoCRBP contains a binding cavity more shielded than that of any other retinoid carrier, conformational flexibility in the portal region may assist retinol uptake. The stiffening of the backbone in the holoprotein guarantees the stability of the complex during retinol transport and suggests that targeted retinol release requires a transiently open state that is likely to be promoted by the acceptor or the local environment.     

Structural dynamics of calmodulin and troponin C

We present the results of computational simulation studies of the structures of calmodulin (CAM) and troponin C (TNC). Possible differences between the structures of these molecules in the crystal and in solution were suggested by results from some recent experimental studies, which implied that their conformations in solution may be more compacted than the characteristic dumbbell shape observed in the crystal. The molecular dynamics simulations were carried out with the CHARMM system of programs, and the environment was modeled with a distance-dependent dielectric permittivity and discrete water molecules surrounding the proteins at starting positions identified in the crystals of CAM and TNC. Methods of macromolecular structure analysis, including linear distance plots, distance matrices and a matrix representation of hydrogen bonding, were used to analyze the nature, the extent and the source of structural differences between the computed structures of the molecules and their conformations in the crystal. Following the longest simulation, in which intradomain structure was conserved, the crystallographically observed dumbbell structure of the molecule changed due to a kinking or bending in the region of the central tether helix connecting the two Ca(2+)-binding domains which moved into close proximity. The resulting structure correlates with experimental observations of complexes between CAM and peptides such as melittin and mastoparan. Analysis of the corresponding pair distance distribution functions in comparison to experimental results suggests the dynamic existence of a non-negligible fraction of the compacted structure in aqueous solutions of CAM. In this more nearly globular shape, CAM reveals to the environment two interior pockets that contain a number of hydrophobic residues, in agreement with NMR data suggesting involvement of such residues in the binding of inhibitors and proteins to CAM.     

Calcium-induced changes in calmodulin structural dynamics and thermodynamics

The thermodynamics of the interaction between Ca(2+) and calmodulin (CaM) was examined using isothermal titration calorimetry (ITC). The chemical denaturation of calmodulin was monitored spectroscopically to determine the stability of Ca(2+)-free (apo) and Ca(2+)-loaded (holo) CaMs. We explored the conformational and structural dynamics of CaM using amide hydrogen-deuterium (H-D) exchange coupled with Fourier transform infrared (FT-IR) spectroscopy. The results of H-D exchange and FT-IR suggest that CaM activation by Ca(2+) binding involves significant conformational changes. The results have also revealed that while the overall conformation of holo-CaM is more stable than that of the apo-CaM, some part of its α-helix structures, most likely the EF-hand domain region, has more solvent exposure, thus, has a faster H-D exchange rate than that of the apo-CaM. The ITC method provides a new strategy for obtaining site-specific Ca(2+) binding properties and a better estimation of the cooperativity and conformational change contributions of coupled EF-hand proteins.     

Structure, interaction, dynamics and solvent effects on the DNA-EcoRI complex in aqueous solution from molecular dynamics simulation.

A 0.7-ns molecular dynamics simulation of the DNA-EcoRI complex in a 7.0-A solvent shell indicated a stable behavior of the system. No significant evaporation or smearing of the solvent's outer boundary occurred. The structure and the intermolecular interactions were found to be well maintained during the simulation. The interaction pattern in the simulation was found to be very similar to that in the crystal structure. Most of the specific interactions between the DNA and the protein were found to be enhanced in the simulation compared to that in the crystal structure as a result of improved interaction geometry. The nonspecific interactions were found to be stronger than the specific ones. The specific interactions between the N7 atoms of Gua(4) or Ade(5) or Ade(6) and the protein were found to be present over almost the entire time of the simulation, whereas hydrogen bonds involving the amino groups of the Ade(5) and Ade(6) with the protein were found to be relatively weaker, with lower probability and shorter lifetime. The time evolution of the root mean square deviations of the DNA and the protein were highly correlated even at the later part of the simulation, showing the tight binding between them. Several long-lived water bridges were found between the DNA backbone atoms and the protein and also between the two protein monomers, which increased the overall stability of the complex. The two protein monomers were found to interact strongly with each other. The energy of the DNA kink deformation was estimated as approximately 31 kcal/mol.     

Structure and dynamics of double helices in solution: modes of DNA bending

The long range structure of DNA restriction fragments has been analysed by electro-optical measurements. The overall rotation time constants observed in a low salt buffer with monovalent ions is shown to decrease upon addition of Mg2+ or spermine. Since the circular dichroism and also the limiting value of the linear dichroism remain almost constant under these conditions, the effect is attributed to a change of the long range structure. According to a weakly bending rod model, the persistence length decreases from about 600 A in the absence of Mg2+ or spermine to about 350 A in the presence of these ions. The persistence length measured in the presence of Mg2+ is almost independent of temperature in the range of 10 to 40 degrees C. The nature of DNA bending is analysed by measurements of bending amplitudes and time constants from dichroism decay curves. The observed absence of changes in the bending amplitudes upon addition of Mg2+ or spermine, even though addition induces changes of the persistence length by a factor of 2, is hardly consistent with simple thermal bending. The combined results, including the remarkably small temperature dependence of persistence length and bending amplitude, can be explained by the existence of two bending effects: inherent curvature of DNA dominates at low temperature, whereas thermal bending prevails at high temperature. Analysis of bending amplitudes from dichroism decay curves according to an arc model provides an approximate measure for the degree of bending in restriction fragments. The model is consistent with the observed chain length dependence of bending amplitudes and provides an approximate curvature corresponding to a radius of about 400 A. Thus the curvature observed in restriction fragments is similar to that observed for high molecular DNA condensed into toroids by addition of ions like spermine. Particularly strong bending of DNA is induced by [Co(NH3)6]3+, indicated by an apparent persistence length of 200 A and an increased bending amplitude together with a reduced limit value of the linear dichroism. This effect is attributed to the high charge density of this ion and potential site binding.     

Structure, dynamics and hydration of the nogalamycin-d(ATGCAT)2Complex determined by NMR and molecular dynamics simulations in solution.

The structure of the 1:1 nogalamycin:d(ATGCAT)2 complex has been determined in solution from high-resolution NMR data and restrained molecular dynamics (rMD) simulations using an explicit solvation model. The antibiotic intercalates at the 5'-TpG step with the nogalose lying along the minor groove towards the centre of the duplex. Many drug-DNA nuclear Overhauser enhancements (NOEs) in the minor groove are indicative of hydrophobic interactions over the TGCA sequence. Steric occlusion prevents a second nogalamycin molecule from binding at the symmetry-related 5'-CpA site, leading to the conclusion that the observed binding orientation in this complex is the preferred orientation free of the complication of end-effects (drug molecules occupy terminal intercalation sites in all X-ray structures) or steric interactions between drug molecules (other NMR structures have two drug molecules bound in close proximity), as previously suggested. Fluctuations in key structural parameters such as rise, helical twist, slide, shift, buckle and sugar pucker have been examined from an analysis of the final 500 ps of a 1 ns rMD simulation, and reveal that many sequence-dependent structural features previously identified by comparison of different X-ray structures lie within the range of dynamic fluctuations observed in the MD simulations. Water density calculations on MD simulation data reveal a time-averaged pattern of hydration in both the major and minor groove, in good agreement with the extensive hydration observed in two related X-ray structures in which nogalamycin is bound at terminal 5'-TpG sites. However, the pattern of hydration determined from the sign and magnitude of NOE and ROE cross-peaks to water identified in 2D NOESY and ROESY experiments identifies only a few "bound" water molecules with long residence times. These solvate the charged bicycloaminoglucose sugar ring, suggesting an important role for water molecules in mediating drug-DNA electrostatic interactions within the major groove. The high density of water molecules found in the minor groove in X-ray structures and MD simulations is found to be associated with only weakly bound solvent in solution.     

Structure and internal dynamics of the bovine pancreatic trypsin inhibitor in aqueous solution from long-time molecular dynamics simulations

Structural and dynamic properties of bovine pancreatic trypsin inhibitor (BPTI) in aqueous solution are investigated using two molecular dynamics (MD) simulations: one of 1.4 ns length and one of 0.8 ns length in which atom-atom distance bounds derived from NMR spectroscopy are included in the potential energy function to make the trajectory satisfy these experimental data more closely. The simulated properties of BPTI are compared with crystal and solution structures of BPTI, and found to be in agreement with the available experimental data. The best agreement with experiment was obtained when atom-atom distance restraints were applied in a time-averaged manner in the simulation. The polypeptide segments found to be most flexible in the MD simulations coincide closely with those showing differences between the crystal and solution structures of BPTI. © 1995 Wiley-Liss, Inc.     

The molecular dynamics of hyaluronates in solution.

The dynamic properties of hyaluronate solutions are discussed with relevance to some problems in sensory physiology (mechanoelectrical transduction), renal physiology, interstitial fluid regulation, and especially to the causes of open-angle glaucoma. With respect to the last problem: as recent biochemical evidence indicates that the hyaloid membrane does not exist, it now seems worthwhile to consider the increase in intraocular pressure present in the eye with glaucoma to be due--at least in the open-angle case--to a change in the specific gravity and hydrophilic nature of the hyaluronic acid in the vitreous body in particular, as well as in the trabecular meshwork. Densimetric experimental evidence indicates that the hyaluronate system could, indeed, produce the pressure changes seen in glaucoma, if intraocular pH changed but slightly. A hypothesis concerning the effect of acetazol amide on intraocular pressure is also presented.     

Structure of crambin in solution, crystal and in the trajectories of molecular dynamics simulations

The mechanisms of the three-dimensional crambin structure alterations in the crystalline environments and in the trajectories of the molecular dynamics simulations in the vacuum and crystal surroundings have been analyzed. In the crystalline state and in the solution the partial regrouping of remote intramolecular packing contacts, involved in the formation and stabilization of the tertiary structure of the crambin molecule, occurs in NMR structures. In the crystalline state it is initiated by the formation of the intermolecular contacts, the conformational influence of its appearance is distributed over the structure. The changes of the conformations and positions of the residues of the loop segments, where the intermolecular contacts of the crystal surroundings are preferably concentrated, are most observable. Under the influence of these contacts the principal change of the regular secondary structure of crambin is taking place: extension of the two-strand β structure to the three-strand structure with the participation of the single last residue N46 of the C-terminal loop. In comparison with the C-terminal loop the more profound changes are observed in the conformation and the atomic positions of the backbone atoms and in the solvent accessibility of the residues of the interhelical loop. In the solution of the ensemble of the 8 NMR structures relative accessibility to the solvent differs more noticeably also in the region of the loop segments and rather markedly in the interhelical loop. In the crambin cryogenic crystal structures the positions of the atoms of the backbone and/or side chain of 14–18 of 46 residues are discretely disordered. The disorganizations of at least 8 of 14 residues occur directly in the regions of the intermolecular contacts and another 5 residues are disordered indirectly through the intramolecular contacts with the residues of the intermolecular contacts. Upon the molecular dynamics simulation in the vacuum surrounding as in the solution of the crystalline structure of crambin the essential changes of the backbone conformation are caused by the intermolecular contacts absence, but partly masked by the structure changes owing to the nonpolar H atoms absence on the simulated structure. The intermolecular contact absence is partly manifested upon the molecular dynamics simulation of the crambin crystal with one protein molecule. Compared to the crystal structure the lengths of the interpeptide hydrogen bonds and other interresidue contacts in an average solution NMR structure are somewhat shorter and accordingly the energy of the interpeptide hydrogen bonds is better. This length shortening can occur at the stage of the refinement of the NMR structures of the crambin and other proteins by its energy minimizations in the vacuum surroundings and not exist in the solution protein structures.     

Structure and dynamics of a membrane protein in micelles from three solution NMR experiments

Three solution NMR experiments on a uniformly ¹⁵N labeled membrane protein in micelles provide sufficient information to describe the structure, topology, and dynamics of its helices, as well as additional information that characterizes the principal features of residues in terminal and inter-helical loop regions. The backbone amide resonances are assigned with an HMQC-NOESY experiment and the backbone dynamics are characterized by a ¹H-¹⁵N heteronuclear NOE experiment, which clearly distinguishes between the structured helical residues and the more mobile residues in the terminal and interhelical loop regions of the protein. The structure and topology of the helices are described by Dipolar waves and PISA wheels derived from experimental measurements of residual dipolar couplings (RDCs) and residual chemical shift anisotropies (RCSAs). The results show that the membrane-bound form of Pf1 coat protein has a 20-residue trans-membrane hydrophobic helix with an orientation that differs by about 90° from that of an 8-residue amphipathic helix. This combination of three-experiments that yields Dipolar waves and PISA wheels has the potential to contribute to high-throughput structural characterizations of membrane proteins.     

Construction and molecular dynamics simulation of calmodulin in the extended and in a bent conformation.

Analysis of sequence similarity and comparison of the three-dimensional (3D) structures of troponin C and calmodulin have revealed a sequence in the central helix of calmodulin with a high probability for bending. The three amino acids known to form a bend in the N-terminal portion of troponin C are also found in the central helix of calmodulin. The modelling of a bent calmodulin structure, using the dihedral angles of the three residues in the bend of troponin C as a 3D template, results in a conformation of calmodulin where the N- and C-terminal domains are able to form contacts. Dynamics simulations starting from the X-ray structure of calmodulin and from the modelled bent calmodulin were carried out to compare flexibility and correlated movements of Ca2+ in the binding loops. Both conformations of calmodulin remained stable during the period of simulation. In the simulation of calmodulin in the extended form, the motions of the Ca2+ atoms in the two domains (Ca2+1 and Ca2+2 in one domain, and Ca2+3 and Ca2+4 in the other) are correlated. In the simulation of the bent form, an additional correlation between the Ca atoms in the two different domains is observed. The results are compatible with the occurrence of a bent conformation of calmodulin in the presence of targets, and with increased Ca2+ affinity and cooperativity of the Ca(2+)-binding loops in the calmodulin-peptide complexes.