Long-range interactions are known to play an important role in highly polar biomolecules like DNA. In molecular dynamics simulations of nucleic acids and proteins, an accurate treatment of the long-range interactions are crucial for achieving stable nanosecond trajectories. In this report, we evaluate the structural and dynamic effects on a highly charged oligonucleotide in aqueous solution from different long-range truncation methods. Two group-based truncation methods, one with a switching function and one with a force-switching function were found to fail to give accurate stable trajectories close to the crystal structure. For these group-based truncation methods, large root mean square (rms) deviations from the initial structure were obtained and severe distortions of the oligonucleotide were observed. Another group-based truncation scheme, which used an abrupt truncation at 8. 0 A or at 12.0 A was also investigated. For the short cutoff distance, the conformations deviated far away from the initial structure and were significantly distorted. However, for the longer cutoff, where all necessary electrostatic interactions were included, the trajectory was quite stable. For the particle mesh Ewald (PME) truncation method, a stable DNA simulation with a heavy atom rms deviation of 1.5 A was obtained. The atom-based truncation methods also resulted in stable trajectories, according to the rms deviation from the initial B-DNA structure, of between 1.5 and 1.7 A for the heavy atoms. In these stable simulations, the heavy atom rms deviations were approximately 0.6-1.0 A lower for the bases than for the backbone. An increase of the cutoff radius from 8 to 12 A decreased the rms deviation by approximately 0.2 A for the atom-based truncation method with a force-shifting function, but increased the computational time by a factor of 2. Increasing the cutoff from 12 to 18 A for the atom-based truncation method with a force-shifting function requires 2-3 times more computational time, but did not significantly change the rms deviation. Similar rms deviations from the initial structure were found for the atom-based method with a force-shifting function and for the PME method. The computational cost was longer for the PME method with a cutoff of 12. 0 A for the direct space nonbonded calculations than for the atom-based truncation method with a force-shifting function and a cutoff of 12.0 A. If a nonperiodic boundary, e.g., a spherical boundary, was used, a considerable speedup could be achieved. From the rms fluctuations, the terminal nucleotides and especially the cytidines were found to be more flexible than the nonterminal nucleotides. The B-DNA form of the oligonucleotide was maintained throughout the simulations and is judged to depend on the parameters of the energy function and not on the truncation method used to handle the long-range electrostatic interactions. To perform accurate and stable simulations of highly charged biological macromolecules, we recommend that the atom-based force-shift method or the PME method...

No abstract

Molecular dynamics simulations are widely used today to tackle problems in biochemistry and molecular biology. In the 25 years since the first simulation of a protein computers have become faster by many orders of magnitude, algorithms and force fields have been improved, and simulations can now be applied to very large systems, such as protein-nucleic acid complexes and multimeric proteins in aqueous solution. In this review we give a general background about molecular dynamics simulations, and then focus on some recent technical advances, with applications to biologically relevant problems.

Currently, we found that the 70-kDa p70 S6 kinase (p70S6K) directly phosphorylates tau at S262, S214, and T212 sites in vitro. By immunoprecipitation, p-p70S6K (T421/S424) showed a close association with p-tau (S262 and S396/404). Zinc-induced p70S6K activation could only upregulate translation of total S6 and tau but not global proteins in SH-SY5Y cells. The requirement of p70S6K activation was confirmed in the SH-SY5Y cells that overexpress wild-type htau40. Level of p-p70S6K (T421/S424) was only significantly correlated with p-tau at S262, S214, and T212, but not T212/S214, in AlzheimerÕs disease (AD) brains. These suggested that p70S6K might contribute to tau related pathologies in AD brains.

Supported liquid membrane extraction is a powerful analytical methodology for continuous extraction and enrichment of ionizable organic pollutants from aqueous environmental samples. Sample workup with this technique can be hampered by low enrichment factors due to incomplete trapping in the stagnant acceptor solution for weak basic compounds. The different parameters, such as acceptor pH and ionic strength of both the donor and acceptor solutions, related to this problem were studied theoretically and experimentally. A simple equation for the maximum enrichment factor was developed. Results showed that the ionic strengths of the two aqueous phases in the system have an effect on the maximum enrichment factor, which can be increased by increasing the ionic strength of the donor phase. For chloro-s-triazines, a satisfactory agreement was obtained between the experimental measurements and theoretical estimations. This study defines the practical analytical applications and limitations of the system where the analytes are partially ionized in the acceptor solution.