Potential of mean force calculations of the stacking-unstacking process in single-stranded deoxyribodinucleoside monophosphates

Internal mobility of the oligonucleotide duplexes d(TCGCG)2 and d(CGCGCG)2 in aqueous solution from molecular dynamics simulations

In this paper we present longitudinal relaxation times, order parameters and effective correlation times for the base and sugar carbons in both strands of the oglionucleotide duplexes d(TCGCG)2 and d(CGCGCG)2, as calculated from 400 ps molecular dynamics trajectories in aqueous solution. The model-free approach (Lipari and Szabo, 1982) was used to determine the amplitudes and time scales of the internal motion. Comparisons were made with NMR relaxation measurements (Borer et al., 1994). The order parameters could acceptably be reproduced, and the effective correlation times were found to be lower than the experimental estimates. Reasonable T1 relaxation times were obtained in comparison with experiment for the nonterminal nucleosides. The T1 relaxation times were found to depend mainly on the order parameters and overall rotational correlation time.

Calculations of nucleic acid conformations

The present computational power and sophistication of theoretical approaches to nucleic acid structural investigation are sufficient for the realization of static and dynamic models that correlate accurately with current crystallographic, NMR and solution-probing structural data, and consequently are able to provide valuable insights and predictions for a variety of nucleic acid conformational families. In molecular dynamics simulations, the year 1995 was marked by the foray of fast Ewald methods, an accomplishment resulting from several years' work in the search for an adequate treatment of the electrostatic long-range forces so primordial in nucleic acid behavior. In very large systems, and particularly in the RNA-folding field, techniques originating from artificial intelligence research, like constraint satisfaction programming or genetic algorithms, have established their utility and potential.

Calculation of translational friction and intrinsic viscosity. II. Application to globular proteins

The translational friction coefficients and intrinsic viscosities of four proteins (ribonuclease A, lysozyme, myoglobin, and chymotrypsinogen A) are calculated using atomic-level structural details. Inclusion of a 0.9-A-thick hydration shell allows calculated results for both hydrodynamic properties of each protein to reproduce experimental data. The use of detailed protein structures is made possible by relating translational friction and intrinsic viscosity to capacitance and polarizability, which can be calculated easily. The 0.9-A hydration shell corresponds to a hydration level of 0.3-0.4 g water/g protein. Hydration levels within this narrow range are also found by a number of other techniques such as nuclear magnetic resonance spectroscopy, infrared spectroscopy, calorimetry, and computer simulation. The use of detailed protein structures in predicting hydrodynamic properties thus allows hydrodynamic measurement to join the other techniques in leading to a unified picture of protein hydration. In contrast, earlier interpretations of hydrodynamic data based on modeling proteins as ellipsoids gave hydration levels that varied widely from protein to protein and thus challenged the existence of a unified picture of protein hydration.