A simple model for predicting the free energy of binding between anthracycline antibiotics and DNA.

A theoretical model for predicting the free energy of binding between anthracycline antibiotics and DNA was developed using the electron density functional (DFT) and molecular mechanics (MM) methods. Partial DFT-ESP charges were used in calculating the MM binding energies for complexes formed between anthracycline antibiotics and oligodeoxynucleotides. These energies were then compared with experimental binding free energies. The good correlation between the experimental and theoretical energies allowed us to propose a model for predicting the binding free energy for derivatives of anthracycline antibiotics and for quickly screening new anthracycline derivatives.

A mezoscopic model of nucleic acids. Part 1. Lagrangian and quaternion molecular dynamics.

This study presents a model for mezoscopic molecular dynamics simulations with objects of different scale and properties e.e. atoms, pseudoatoms, rigid and pseudo-elastic bodies, described by the external coordinates and internal degrees of freedom. The Lagrangian approach is used to derive equations of motion and a quaternion representation is used for the description of the dynamics of rigid and pseudo-elastic molecular elements. Stability of the LQMD algorithm was tested for a 10-base pair deoxynucleotide. The total energy, momentum and angular momentum are conserved for time-steps up to 20 fs.

A mezoscopic model of nucleic acids. Part 2. An effective potential energy function for DNA.

In this study we present an effective Potential of Mean Force (PMF) designed for Lagrangian and Quaternion Molecular Dynamics (LQMD) of DNA. The DNA model is built from pseudoatoms as well as rigid and pseudo-elastic bodies described by a limited number of selected Cartesian and internal degrees of freedom. Phosphate groups, deoxyribose rings and nucleic acid bases are represented by pseudoparticles, some of them with internal degrees of freedom. PMF is defined as the sum of effective bonded and long-range potentials. The potentials were fitted to numerical free energy surfaces. Over 50 free energy surfaces, each depending on a conformational variable (pseudobond length, angle or dihedral angle) and the pseudorotation phase of the nearest neighbour deoxribose ring, were computed. The numerical free energy surfaces were obtained from probability distributions derived from a 1.5 ns conventional, microscopic MD simulation of the d(GpC)9 double helical DNA molecule. An umbrella sampling method was used to simulate transitions between the A and B DNA forms, and PMF reproduces these transitions.

Binding of SPXK- and APXK-peptide motifs to AT-rich DNA. Experimental and theoretical studies.

The binding properties of the SPXK- and APXK-type peptides to the AT-rich DNA fragments of different length were studied by measuring the competition of peptides with Hoechst 33258 dye for DNA binding and by the gel shift assay analysis. In parallel to the experimental studies, molecular modeling techniques were used to analyze possible binding modes of the SPXZ and APXK motifs to the AT-rich DNA. The results of the competition measurements and gel shift assays suggest that serine at the i-1 position (i is proline) can be replaced by alanine without affecting the binding properties of the motif. Thus, the presence of the conserved serine in this motif in many DNA-binding proteins is probably not dictated by structural requirements. Based on the results of molecular modeling studies we propose that the binding mode of the SPXK-type motifs to the AT-rich DNA resembles closely that between the N-terminal arm of the homeodomain and DNA. This model confirms that serine in the SPXK motifs is not essential for the DNA binding. The model also indicates that if X in the motif is glutamic acid, this residue is probably protonated in the complex with DNA.

Modeling the DNA-solvent interface.

We extend the technique of using perpendicular distribution functions to salt solutions around nucleic acids. Both solute density averaged and nonaveraged reference frames are considered and compared. Using a previous simulation of DNA in salt water of over a nanosecond in duration, the aqueous distribution functions were found to be well coveraged, whereas the salt perpendicular distribution functions were less well determined. Three-dimensional density reconstructions reliably showed the prominent solvation features with transferable functions. The number of solute atom types needed for reconstructions of a given precision was determined in the context of the reference simulation data set with the goal of achieving a required level of reconstruction quality.

Lagrangian molecular dynamics using selected conformational degrees of freedom, with application to the pseudorotation dynamics of furanose rings.

Using internal conformational degrees of freedom for biopolymers as natural variables, and introducing a Lagrangian dynamics approach, one can simulate time-dependent processes over a much longer time scale than in classical Newtonian molecular dynamics (MD) techniques. Two factors contribute to this: a substantial reduction in the number of degrees of freedom and a very large increase in the size of the time step. We present the Lagrangian equations of motion for repuckering transitions in model furanose (F), ribose (R), and 2'-deoxyribose (dR) ring systems using the pseudorotation phase angle as the single dynamic variable. As in most Lagrangian analyses, the effective masses for the R and dR models are dependent on conformation, and we test the behavior of this variable mass (VM) model. Since the variation in effective mass is small, the VM model is compared with a simplified constant mass (CM) model, which is shown to be an excellent approximation. The equations of motion for the CM and VM models are integrated with the leapfrog and the iterative leapfrog algorithms, respectively. The Lagrangian dynamics approach reduces the number of degrees of freedom from about 40 to 1, and allows the use of time steps on the order of 20 fs, about an order of magnitude greater than is used in conventional MD simulations.