Structure and dynamics of calcium-activated calmodulin in solution.

Two 4-ns molecular dynamics simulations of calcium loaded calmodulin in solution have been performed, using both standard nonbonded cutoffs and Ewald summation to treat electrostatic interactions. Our simulation results are generally consistent with solution experimental studies of calmodulin structure and dynamics, including NMR, cross-linking, fluorescence and x-ray scattering. The most interesting result of the molecular dynamics simulations is the detection of large-scale structural fluctuations of calmodulin in solution. The globular N- and C-terminal domains tend to move approximately like rigid bodies, with fluctuations of interdomain distances within a 7 A range and of interdomain angles by up to 60 deg. Essential dynamics analysis indicates that the three dominant types of motion involve bending of the central helix in two perpendicular planes and a twist in which the domains rotate in opposite directions around the central helix. In the more realistic Ewald trajectory the protein backbone remains mostly within a 2-3 A root-mean-square distance from the crystal structure, the secondary structure within the domains is conserved and middle part of the central helix becomes disordered. The central helix itself exhibits limited fluctuations, with its bend angle exploring the 0-50 degrees range and the end-to-end distance falling in 39-43 A. The results of the two simulations were similar in many respects. However, the cutoff trajectory exhibited a larger deviation from the crystal, loss of several helical hydrogen bonds in the N-terminal domain and lack of structural disorder in the central helix.     

The fractal structure of homeopathic preparations

It was demonstrated (by the MFL method) in experiments the collective dynamical structural states of the liquid medium are the basis of the biological activity of homeopathic preparations (HP). The dynamical structure of water and HP determines the specific structure of the proton potentials in these media and the EM radiation of very low level and uniform spectrum. The vibration structure in medium (specific molecular dynamics, proton potential structure and spectrum of very low level eradiation) have fractal organization--spectrum of processes from high-frequency region (intramolecular vibrations) to low-frequency region (creation and destruction of clatrates) consists of consequence of corresponding regions with the repetition of frequency correlation in every region.     

Hydration of the dTn.dAn x dTn parallel triple helix: a Fourier transform infrared and gravimetric study correlated with molecular dynamics simulations.

We present a comparative analysis of the water organization around the dTn.dAn x dTn triple helix and the Watson-Crick double helix dTn.dAn respectively by means of gravimetric measurements, infrared spectroscopy and molecular dynamics simulations. The hydration per nucleotide determined by gravimetric and spectroscopic methods correlated with the molecular dynamics simulations shows that at high relative humidity (98% RH) the triple helix is less solvated than the duplex (17 +/- 2 water molecules per nucleotide instead of 21 +/-1). The experimental desorption curves are different for both structures and indicate that below 81% RH the triplex becomes more hydrated than the duplex. At this RH the FTIR spectra show the emergence of N-type sugars in the adenosine strand of the triplex. When the third strand is bound in the major groove of the Watson-Crick duplex molecular dynamics simulations show the formation of a spine of water molecules between the two thymidine strands.     

Water-mediated hydrogen-bonded network on the cytoplasmic side of the Schiff base of the L photointermediate of bacteriorhodopsin

The L intermediate in the proton-motive photocycle of bacteriorhodopsin is the starting state for the first proton transfer, from the Schiff base to Asp85, in the formation of the M intermediate. Previous FTIR studies of L have identified unique vibration bands caused by the perturbation of several polar amino acid side chains and several internal water molecules located on the cytoplasmic side of the retinylidene chromophore. In the present FTIR study we describe spectral features of the L intermediate in D(2)O in the frequency region which includes the N-D stretching vibrations of the backbone amides. We show that a broad band in the 2220-2080 cm(-1) region appears in L. By use of appropriate (15)N labeling and mutants, the lower frequency side of this band in L is assigned to the amides of Lys216 and Gly220. These amides are coupled to each other, and interact with Thr46 and Val49 in helix B and Asp96 in helix C via weakly H-bonding water molecules that exhibit O-D stretching vibrations at 2621 and 2605 cm(-1). These water molecules are part of a hydrogen-bonded network characteristic of L which includes other water molecules located closer to the chromophore that exhibit an O-D stretching vibration at 2589 cm(-1). This structure, extending from the Schiff base to the internal proton donor Asp96, stabilizes L and affects the L-to-M transition.     

Recognition complex between the HMG domain of LEF-1 and its cognate DNA studied by molecular dynamics simulations with explicit solvation

Molecular dynamics simulations of the complex formed between the HMG box of the lymphoid enhancer-binding factor (LEF-1) and its cognate DNA duplex were carried out with explicit inclusion of water. The simulation started with an NMR-based model (pdb code 2LEF) and the dynamics was pursued for 10 nanoseconds without constraints. It revealed that water intervenes in many ionic/polar interactions, establishing in particular local equilibria between direct and water-mediated hydrogen bonds, and thus increasing the entropy of the complex. Quite unexpectedly, the simulation indicated that a binding pocket for a specific water molecule may be reversibly formed at the apex of the bend induced in the DNA helix by LEF-1 binding, where a methionine side chain intercalates between two destacked adenines. We observed that the specific water molecule can temporarily replace the intercalated S-CH(3) group, acting as a sort of "extension" of the side chain. The residence time of this water molecule was about 3.5 ns. Simulations of the cognate DNA alone showed that this sequence has no intrinsic tendency to bend; therefore, the bending occurs solely as a consequence of the recognition, following the "induced-fit" mechanism.     

General dynamic properties of Abeta12-36 amyloid peptide involved in Alzheimer's disease from unfolding simulation

To study the folding/unfolding properties of a beta-amyloid peptide Abeta(12-36) of Alzheimer's disease, five molecular dynamics simulations of Abeta(12-36) in explicit water were done at 450 K starting from a structure that is stable in trifluoroethanol/water at room temperature with two alpha-helices. Due to high temperature, the initial helical structure unfolded during the simulation. The observed aspects of the unfolding were as follows. 1) One helix (helix 1) had a longer life than the other (helix 2), which correlates well with the theoretically computed Phi values. 2) Temporal prolongation of helix 1 was found before unfolding. 3) Hydrophobic cores formed frequently with rearrangement of amino-acid residues in the hydrophobic cores. The formation and rearrangement of the hydrophobic cores may be a general aspect of this peptide in the unfolded state, and the structural changes accompanied by the hydrophobic-core rearrangement may lead the peptide to the most stable structure. 4) Concerted motions (collective modes) appeared to unfold helix 1. The collective modes were similar with those observed in another simulation at 300 K. The analysis implies that the conformation moves according to the collective modes when the peptide is in the initial stage of protein unfolding and in the final stage of protein folding.     

Transmembrane helix structure, dynamics, and interactions: multi-nanosecond molecular dynamics simulations.

To probe the fundamentals of membrane/protein interactions, all-atom multi-nanosecond molecular dynamics simulations were conducted on a single transmembrane poly(32)alanine helix in a fully solvated dimyristoyphosphatidylcholine (DMPC) bilayer. The central 12 residues, which interact only with the lipid hydrocarbon chains, maintained a very stable helical structure. Helical regions extended beyond these central 12 residues, but interactions with the lipid fatty-acyl ester linkages, the lipid headgroups, and water molecules made the helix less stable in this region. The C and N termini, exposed largely to water, existed as random coils. As a whole, the helix tilted substantially, from perpendicular to the bilayer plane (0 degree) to a 30 degrees tilt. The helix experienced a bend at its middle, and the two halves of the helix at times assumed substantially different tilts. Frequent hydrogen bonding, of up to 0.7 ns in duration, occurred between peptide and lipid molecules. This resulted in correlated translational diffusion between the helix and a few lipid molecules. Because of the large variation in lipid conformation, the lipid environment of the peptide was not well defined in terms of "annular" lipids and on average consisted of 18 lipid molecules. When compared with a "neat" bilayer without peptide, no significant difference was seen in the bilayer thickness, lipid conformations or diffusion, or headgroup orientation. However, the lipid hydrocarbon chain order parameters showed a significant decrease in order, especially in those methylene groups closest to the headgroup.     

Brownian dynamics simulations of supercoiled DNA with bent sequences.

The recently presented Brownian dynamics model for superhelical DNA is extended to include local curvature of the DNA helix axis. Here we analyze the effect of a permanent bend on the structure and dynamics of an 1870-bp superhelix with delta Lk = -10. Furthermore, we define quantitative expressions for computing structural parameters such as loop positions, superhelix diameter, and plectonemic content for trajectories of superhelical DNA, and assess the convergence toward global equilibrium. The structural fluctuations in an interwound superhelix, as reflected in the change in end loop positions, seem to occur by destruction/creation of loops rather than by a sliding motion of the DNA around its contour. Their time scale is on the order of 30-100 microseconds. A permanent bend changes the structure and the internal motions of the DNA drastically. The position of the end loop is fixed at the permanent bend, and the local motions of the chain are enhanced near the loops. A displacement of the bend from the end loop to a position inside the plectonemic part of the superhelix results in the formation of a new loop and the disappearance of the old one; we estimate the time involved in this process to be about 0.5 ms.     

Spontaneous conformational changes in the E. coli GroEL subunit from all-atom molecular dynamics simulations

The Escherichia coli chaperonin GroEL is a complex of identical subunit proteins (57 kDa each) arranged in a back-to-back stacking of two heptameric rings. Its hallmarks include nested positive intra-ring and negative inter-ring cooperativity in adenosine trisphosphate (ATP) binding and the ability to mediate the folding of newly transcribed and/or denatured substrate proteins. We performed unbiased molecular dynamics simulations of the GroEL subunit protein in explicit water both with and without the nucleotide KMgATP to understand better the details of the structural transitions that enable these behaviors. Placing KMgATP in the equatorial domain binding pocket of a t state subunit, which corresponds to a low ATP-affinity state, produced a short-lived (6 ns) state that spontaneously transitioned to the high ATP-affinity r state. The important feature of this transition is a large-scale rotation of the intermediate domain's helix M to close the ATP binding pocket. Pivoting of helix M is accompanied by counterclockwise rotation and slight deformation of the apical domain, important for lowering the affinity for substrate protein. Aligning simulation conformations into model heptamer rings demonstrates that the t-->r transition in one subunit is not sterically hindered by t state neighbors, but requires breakage of Arg(197)-Glu(386) intersubunit salt bridges, which are important for inter-ring positive cooperativity. Lowest-frequency quasi-harmonic modes of vibration computed pre- and post-transition clearly show that natural vibrations facilitate the transition. Finally, we propose a novel mechanism for inter-ring cooperativity in ATP binding inspired by the observation of spontaneous insertion of the side chain of Ala(480) into the empty nucleotide pocket.     

Structural change of threonine 89 upon photoisomerization in bacteriorhodopsin as revealed by polarized FTIR spectroscopy.

The all-trans to 13-cis photoisomerization of the retinal chromophore of bacteriorhodopsin occurs selectively, efficiently, and on an ultrafast time scale. The reaction is facilitated by the surrounding protein matrix which undergoes further structural changes during the proton-transporting reaction cycle. Low-temperature polarized Fourier transform infrared difference spectra between bacteriorhodopsin and the K intermediate provide the possibility to investigate such structural changes, by probing O-H and N-H stretching vibrations [Kandori, Kinoshita, Shichida, and Maeda (1998) J. Phys. Chem. B 102, 7899-7905]. The measurements of [3-18O]threonine-labeled bacteriorhodopsin revealed that one of the D2O-sensitive bands (2506 cm(-1) in bacteriorhodopsin and 2466 cm(-1) in the K intermediate, in D2O exhibited 18(O)-induced isotope shift. The O-H stretching vibrations of the threonine side chain correspond to 3378 cm(-1) in bacteriorhodopsin and to 3317 cm(-1) in the K intermediate, indicating that hydrogen bonding becomes stronger after the photoisomerization. The O-H stretch frequency of neat secondary alcohol is 3340-3355 cm(-1). The O-H stretch bands are preserved in the T46V, T90V, T142N, T178N, and T205V mutant proteins, but diminished in T89A and T89C, and slightly shifted in T89S. Thus, the observed O-H stretching vibration originates from Thr89. This is consistent with the atomic structure of this region, and the change of the S-H stretching vibration of the T89C mutant in the K intermediate [Kandori, Kinoshita, Shichida, Maeda, Needleman, and Lanyi (1998) J. Am. Chem. Soc. 120, 5828-5829]. We conclude that all-trans to 13-cis isomerization causes shortening of the hydrogen bond between the OH group of Thr89 and a carboxyl oxygen atom of Asp85.     

Theory of Low-Frequency Vibrations in DNA Macromolecules

Abstract

To describe low-frequency dynamics of DNA macromolecules a model is developed taking into account the hydrogen bond stretching in base pairs, the backbone flexibility and intranucleoside mobility. For double-stranded DNA the normal vibrations are found and the structure of low- frequency spectrum is determined. The agreement between theory and Raman spectroscopy data for B-DNA is demonstrated. Conformational dependences of vibration spectrum during the B→A and helix→coil DNA transitions are studied. The contribution coming from low-frequency mobility to the nucleic-protein recognition processes is discussed.     

Molecular dynamics simulations of the glucocorticoid receptor DNA-binding domain in complex with DNA and free in solution.

Molecular dynamics simulations have been performed on the glucocorticoid receptor DNA binding domain (GR DBD) in aqueous solution as a dimer in complex with DNA and as a free monomer. In the simulated complex, we find a slightly increased bending of the DNA helix axis compared with the crystal structure in the spacer region of DNA between the two half-sites that are recognized by GR DBD. The bend is mainly caused by an increased number of interactions between DNA and the N-terminal extended region of the sequence specifically bound monomer. The recognition helices of GR DBD are pulled further into the DNA major groove leading to a weakening of the intrahelical hydrogen bonds in the middle of the helices. Many ordered water molecules with long residence times are found at the intermolecular interfaces of the complex. The hydrogen-bonding networks (including water bridges) on either side of the DNA major groove involve residues that are highly conserved within the family of nuclear receptors. Very similar hydrogen-bonding networks are found in the estrogen receptor (ER) DBD in complex with DNA, which suggests that this is a common feature for proper positioning of the recognition helix in ER DBD and GR DBD.     

Molecular dynamics study of onset of water gelation around the collagen triple helix

The onset of water gelation around a collagen-like triple helix peptide was studied at ambient temperature and pressure by performing Molecular Dynamics simulations. The radial distribution functions of the oxygen and hydrogen atoms of water are distorted below 4 A from the peptide. The distortion is accompanied by the breakdown of the tetrahedral coordination of the hydrogen-bonded network of water molecules. The water shell around the peptide consists of alternating regions of higher and lower density. In agreement with experiments we find that the first hydration shell is kinetically labile, with a residence time in the order of picoseconds for a water molecule. From the computed diffusion coefficient, a key measure of the collective dynamics, we estimate the average diffusion speed decreases by a factor of 1.5 close to the peptide compared to the liquid. Our results give new insight in gel formation and structure on a molecular level.     

Distortions induced in double-stranded oligonucleotides by the binding of cis- or trans-diammine-dichloroplatinum(II) to the d(GTG) sequence.

Conformational changes induced in double-stranded oligonucleotides by the binding of trans- or cis-diamminedichloro platinum(II) to the d(GTG) sequence have been characterized by means of melting temperatures, electrophoretic migrations in non-denaturing polyacrylamide gels, reactivities with the artificial nuclease Phenanthroline-copper and with chemical probes. The cis-platinum adduct behaves more as a centre of directed bend than as a hinge joint, the induced bend angle being of the order of 25-30 degrees. The double helix is locally denatured over 2 base pairs (corresponding to the platinated 5'G residue and the central T residue) and is distorted over 4-5 base pairs. The trans-platinum adduct behaves also more as a centre of directed bend than as a hinge joint, the induced bend angle being of the order of 60 degrees. The double helix is locally denatured over 4 base pairs (corresponding to the immediately 5'T residue adjacent to the adduct and to the three base residues of the adduct). Both the cis- and trans-platinum adducts decrease the thermal stability of the double helix.     

Local dynamics of DNA probed with optical absorption spectroscopy of bound ethidium bromide.

We have studied the local dynamics of calf thymus double-helical DNA by means of an "optical labeling" technique. The study has been performed by measuring the visible absorption band of the cationic dye ethidium bromide, both free in solution and bound to DNA, in the temperature interval 360-30 K and in two different solvent conditions. The temperature dependence of the absorption line shape has been analyzed within the framework of the vibronic coupling theory, to extract information on the dynamic properties of the system; comparison of the thermal behavior of the absorption band of free and DNA-bound ethidium bromide gave information on the local dynamics of the double helix in the proximity of the chromophore. For the dye free in solution, large spectral heterogeneity and coupling to a "bath" of low-frequency (soft) modes is observed; moreover, anharmonic motions become evident at suitably high temperatures. The average frequency of the soft modes and the amplitude of anharmonic motions depend upon solvent composition. For the DNA-bound dye, at low temperatures, heterogeneity is decreased, the average frequency of the soft modes is increased, and anharmonic motions are hindered. However, a new dynamic regime characterized by a large increase in anharmonic motions is observed at temperatures higher than approximately 280 K. The DNA double helix therefore appears to provide, at low temperatures, a rather rigid environment for the bound chromophore, in which conformational heterogeneity is reduced and low-frequency motions (both harmonic vibrations and anharmonic contributions) are hindered. The system becomes anharmonic at approximately 180 K; however, above approximately 280 K, anharmonicity starts to increase much more rapidly than for the dye free in solution; this can be attributed to the onset of wobbling of the dye in its intercalation site, which is likely connected with the onset of (functionally relevant) DNA motions, involving local opening/unwinding of the double helix. As shown by parallel measurements of the melting curves, these motions precede the melting of the double helix and depend upon solvent composition much more than does the melting itself.     

Matching Rules for Collective Behaviors on Complex Networks: Optimal Configurations for Vibration Frequencies of Networked Harmonic Oscillators

The structure-dynamics-function has become one of central problems in modern sciences, and it is a great challenge to unveil the organization rules for different dynamical processes on networks. In this work, we study the vibration spectra of the classical mass spring model with different masses on complex networks, and pay our attention to how the mass spatial configuration influences the second-smallest vibrational frequency () and the largest one (). For random networks, we find that becomes maximal and becomes minimal if the node degrees are point-to-point-positively correlated with the masses. In these cases, we call it point-to-point matching. Moreover, becomes minimal under the condition that the heaviest mass is placed on the lowest-degree vertex, and is maximal as long as the lightest mass is placed on the highest-degree vertex, and in both cases all other masses can be arbitrarily settled. Correspondingly, we call it single-point matching. These findings indicate that the matchings between the node dynamics (parameter) and the node position rule the global systems dynamics, and sometimes only one node is enough to control the collective behaviors of the whole system. Therefore, the matching rules might be the common organization rules for collective behaviors on networks.     

Thermodynamics of the pseudo‐knot in helix 18 of 16S ribosomal RNA

A fragment of E. coli 16S rRNA formed by nucleotides 500 to 545 is termed helix 18. Nucleotides 505‐507 and 524‐526 form a pseudo‐knot and its distortion affects ribosome function. Helix 18 isolated from the ribosome context is thus an interesting fragment to investigate the structural properties and folding of RNA with pseudo‐knots. With all‐atom molecular dynamics simulations, spectroscopic and gel electrophoresis experiments, we investigated thermodynamics of helix 18, with a focus on its pseudo‐knot. In solution studies at ambient conditions we observed dimerization of helix 18. We proposed that the loop, containing nucleotides forming the pseudo‐knot, interacts with another monomer of helix 18. The native dimer is difficult to break but introducing mutations in the pseudo‐knot indeed assured a monomeric form of helix 18. Molecular dynamics simulations at 310 K confirmed the stability of the pseudo‐knot but at elevated temperatures this pseudo‐knot was the first part of helix 18 to lose the hydrogen bond pattern. To further determine helix 18 stability, we analyzed the interactions of helix 18 with short oligomers complementary to a nucleotide stretch containing the pseudo‐knot. The formation of higher‐order structures by helix 18 impacts hybridization efficiency of peptide nucleic acid and 2'‐O methyl RNA oligomers.     

Biological functions of low-frequency vibrations (phonons). III. Helical structures and microenvironment.

Low-frequency vibrations in biomacromolecules possess significant biological functions. In this paper, the alpha-helix element is compared with a mass-distributed spring. Based on this, a set of intuitive and easily handled equations are derived for predicting the fundamental frequencies of helical structures in protein molecules. As shown in the equations, the fundamental frequency depends not only on the constituents of a helix itself but also on its microenvironment. The calculated results agree with the observations. The calculations also demonstrate that the low-frequency vibrations with wave number of approximately 30 cm-1 do not necessarily arise from motions that involve either all or very large portions of the protein molecule as previously thought; a piece of helix containing more than 10 residues and surrounded by a proper microenvironment can also generate such low-frequency motions. Furthermore , we illustrate that the low-frequency motions are closely related to the native state of a protein molecule. Upon denaturation, which is accompanied by a radical change of the relevant microenvironment, the original fundamental frequency also disappears. Consequently, this kind of special frequency termed activating low frequency can serve as a dynamic criterion in identifying whether a biomacromolecule is in its native state. The energy of a phonon excited by this kind of low-frequency vibration is of the same order of magnitude as the average enthalpy value per residue measured during conformational change in some protein molecules. Therefore, there must be some intrinsic relation between the allosteric transitions of protein molecules and their low-frequency motions.     

Molecular dynamics of spermine-DNA interactions: sequence specificity and DNA bending for a simple ligand.

We used molecular dynamics to model interactions between the physiologically important polyamine spermine and two B-DNA oligomers, the homopolymer (dG)10-(dC)10 and the heteropolymer (dGdC)5-(dGdC)5. Water and counterions were included in the simulation. Starting coordinates for spermine-DNA complexes were structures obtained by molecular mechanics modeling of spermine with the two oligomers; in these models, spermine binding induced a bend in the heteropolymer but not in the homopolymer. During approximately 40 psec of molecular dynamics simulation, spermine moves away from the floor of the major groove and interacts nospecifically with d(G)10-d(C)10. In contrast, a spermine-induced bend in the helix of (dGdC)5-(dGdC)5 is maintained throughout the simulation and spermine remains closely associated with the major groove. These results provide further evidence that the binding of spermine to nucleic acids can be sequence specific and that bending of alternating purine-pyrimidine sequences may be a physiologically important result of spermine binding.