Molecular dynamics simulation techniques have been used to investigate the effect of 2,2,2-trifluoroethanol (TFE) as a cosolvent on the stability of three different secondary structure-forming peptides: the α-helix from Melittin, the three-stranded β-sheet peptide Betanova, and the β-hairpin 41–56 from the B1 domain of protein G. The peptides were studied in pure water and 30% (vol/vol) TFE/water mixtures at 300 K. The simulations suggest that the stabilizing effect of TFE is induced by the preferential aggregation of TFE molecules around the peptides. This coating displaces water, thereby removing alternative hydrogen-bonding partners and providing a low dielectric environment that favors the formation of intrapeptide hydrogen bonds. Because TFE interacts only weakly with nonpolar residues, hydrophobic interactions within the peptides are not disrupted. As a consequence, TFE promotes stability rather than inducing denaturation.

Bromine–oxygen radial distribution functions [g(r)] have been calculated by means of molecular dynamics simulations for aqueous solutions of rubidium bromide, 2‐bromopropane and bromoethane. X‐ray absorption spectra at the bromine K edge have been recorded for these solutions. The water contribution to the extended x‐ray absorption fine structure spectra has been calculated starting from the g_Br,O(r) distribution function. Fits of the x‐ray absorption spectra have been performed directly on the raw experimental data, allowing the reliability of the g(r) distribution functions to be verified. The agreement between theoretical and experimental spectra is satisfactory. A procedure to improve model g(r) functions on the basis of the short‐range structural information provided by extended x‐ray absorption fine structure data is proposed.

Comparative explicit solvent molecular dynamics (MD) simulations have been performed on a complete nucleosome core particle with and without N‐terminal histone tails for more than 20 ns. Main purpose of the simulations was to study the dynamics of mobile elements such as histone N‐terminal tails and how packing and DNA‐bending influences the fine structure and dynamics of DNA. Except for the tails, histone and DNA molecules stayed on average close to the crystallographic start structure supporting the quality of the current force field approach. Despite the packing strain, no increase of transitions to noncanonical nucleic acid backbone conformations compared to regular B‐DNA was observed. The pattern of kinks and bends along the DNA remained close to the experiment overall. In addition to the local dynamics, the simulations allowed the analysis of the superhelical mobility indicating a limited relative mobility of DNA segments separated by one superhelical turn (mean relative displacement of approximately ±0.2 nm, mainly along the superhelical axis). An even higher rigidity was found for relative motions (distance fluctuations) of segments separated by half a superhelical turn (approximately ±0.1 nm). The N‐terminal tails underwent dramatic conformational rearrangements on the nanosecond time scale toward partially and transiently wrapped states around the DNA. Many of the histone tail changes corresponded to coupled association and folding events from fully solvent‐exposed states toward complexes with the major and minor grooves of DNA. The simulations indicate that the rapid conformational changes of the tails can modulate the DNA accessibility within a few nanoseconds. © 2007 Wiley Periodicals, Inc. Biopolymers 85: 407–421, 2007.This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com

In this paper, a method of simulating the docking of small flexible ligands to flexible receptors in water is reported. The method is based on molecular dynamics simulations and is an extension of an algorithm previously reported by Di Nola et al. (Di Nola et al., Proteins 1994;19:174–182). The method allows a fast exploration of the receptor surface, using a high temperature of the center of mass translational motion, while the ligand internal motions, the solvent, and the receptor are simulated at room temperature. In addition, the method allows a fast center of mass motion of the ligand, even in solution. The dampening effect of the solvent can be overcome by applying different weights to the interactions between system subsets (solvent, receptor, and ligand). Specific ligand–receptor distances have been used to compare the results of the simulations with the crystal structure. The method is applied, as a test system, to the docking of the phosphocholine to the immunoglobulin McPC603. The results show the similarity of structure between the complex in solution and in the crystal. Proteins 1999;35:153–162.