Solvating,manipulating, damaging,and repairing DNA in a computer

This work highlights four different topics in modeling of DNA: (i) the importance of water and ions together with the structure and function of DNA; the hydration structure around the ions appears to be the determining factor in the ion coordination to DNA, as demonstrated in the results of our MD simulations; (ii) how MD simulations can be used to simulate single molecule manipulation experiments as a complement to reveal the structural dynamics of the studied biomolecules; (iii) how damaged DNA can be studied in computer simulations; and (iv) how repair of damaged DNA can be studied theoretically. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Electrostatics of phosphatidic acid monolayers: Insights from computer simulations

In this paper, we argue that many of the fascinating electrostatic effects that take place in amphiphilic systems are strongly related to the particular organization of the oxygen atoms within each individual molecule. In particular, we focus on two effects: charge inversion and dielectric overscreening. For that purpose, we present molecular dynamics simulations of phosphatidic acid (DMPA²⁻) in the presence of divalent counterions. Our results show that the many oxygens present in DMPA²⁻ cooperatively create strong binding sites for counterions, which in some cases lead to charge inversion. We also present an analysis of the role of interfacial water and relate our analysis to the phenomenon of dielectric overscreening. Several experimental implications are discussed in the conclusions.     

On the relationship betweenT 1 andT 2 for stochastic relaxation models

We consider the relaxation dynamics of two quantum levels coupled to a stochastic bath. We emphasize that even if the matrix elements of the fluctuating Hamiltonian are Gaussian, a second-order cumulant truncation is not exact. For various stochastic models, including the case of a spin-1/2 particle in a fluctuating magnetic field, we calculate 1/T ₁, the population relaxation rate, and 1/T ₂, the phase relaxation rate, up to fourth order in perturbation theory. We show that unlike the commonly accepted second-order result that 1/T ₂1/2T ₁, when fourth-order terms are included, in some instances 1/T ₂<1/2T ₁.     

Analogy between the Lorenz strange attractor and a bistable stochastic oscillator

The relation between the aperiodic solution of the Lorenz model and that of a stochastic anharmonic oscillator is explored. The stochastic oscillator is constructed by replacing (t) in the Lorenz model by a stochastic variable(t) of specified statistics. The resulting system is of course not isomorphic to the Lorenz model, but does share with it a number of statistical properties. Thus, within the confines of these measures the two systems are physically very similar.     

Short range photoinduced electron transfer in proteins: QM-MM simulations of tryptophan and flavin fluorescence quenching in proteins

Hybrid quantum mechanical-molecular mechanics (dynamics) were performed on flavin reductase (Fre) and flavodoxin reductase (Fdr), both from Escherichia coli. Each was complexed with riboflavin (Rbf) or flavin mononucleotide (FMN). During 50 ps trajectories, the relative energies of the fluorescing state (S₁) of the isoalloxazine ring and the lowest charge transfer state (CT) were assessed to aid prediction of fluorescence lifetimes that are shortened due to quenching by electron transfer from tyrosine. The simulations for the four cases display a wide range in CT–S₁ energy gap caused by the presence of phosphate, other charged and polar residues, water, and by intermolecular separation between donor and acceptor. This suggests that the Gibbs energy change (ΔG₀) and reorganization energy (λ) for the electron transfer may differ in different flavoproteins.     

Incorporating the effect of ionic strength in free energy calculations using explicit ions

The incorporation of explicit ions to mimic the effect of ionic strength or to neutralize the overall charge on a system in free energy calculations using molecular dynamics simulations is investigated. The difference in the free energy of hydration between two triosephosphate isomerase inhibitors calculated at five different ion concentrations is used as an example. We show that the free energy difference can be highly sensitive to the presence of explicit ions even in cases where the mutation itself does not involve a change in the overall charge. The effect is most significant if the molecule carries a net charge close to the site mutated. Furthermore, it is shown that the introduction of a small number of ions can lead to very severe sampling problems suggesting that in practical calculations convergence can best be achieved by incorporating either no counterions or by simulating at high ionic strength to ensure sufficient sampling of the ion distribution.     

Molecular Dynamics Simulations of Energetic Solids

A continuing objective in the area of energetic materials is to reduce sensitivity toward impact and shock. One approach is to develop a better understanding of how factors related to the crystal lattice, e.g., defects, influence the initiation and propagation of detonation. Molecular dynamics is a useful tool for this purpose. This paper presents an overview of molecular dynamics treatments of energetic solids. Some of these have simulated initiation and propagation in idealized systems; others have focused on developing a satisfactory procedure for describing molecular crystals of practical significance. Our emphasis in this discussion is on the progress that has been made along the second lines.     

Monte Carlo Simulation of the Assembly of bis‐Biotinylated DNA and Streptavidin

We present Monte Carlo simulations of the self‐assembly of bivalent bis‐biotinylated DNA molecules with the tetravalent biotin‐binding protein streptavidin (STV). By fitting the STV binding probabilities for the four possible valencies, the modelling correctly reproduces the dependencies of various network parameters experimentally observed in an earlier study. The combined results from the experimental and theoretical studies suggest that the binding probability for divalent STV formation is about 50 times larger than for the formation of trivalent and about 200 times larger than for tetravalent STV. In accordance with the experimental results, the modelling also indicates that the mixture of an equimolar ratio of DNA and STV leads to a maximum in size of the oligomeric DNA–STV clusters formed. Furthermore, we found a percolation transition in which the DNA cluster size increases rapidly with increasing DNA concentration resulting in the formation of a single supercluster at elevated concentrations. This behaviour coincides with the occurrence of an immobile band previously observed in electrophoretic experiments, indicating the formation of extremely large DNA–STV aggregate networks.     

Anion mediated polytype selectivity among the basic salts of Co(II)

Basic salts of Co(II) crystallize in the rhombohedral structure. Two different polytypes, 3R₁ and 3R₂, with distinct stacking sequences of the metal hydroxide slabs, are possible within the rhombohedral structure. These polytypes are generated by simple translation of successive layers by (2/3, 1/3, z) or (1/3, 2/3, z). The symmetry of the anion and the mode of coordination influences polytype selection. Cobalt hydroxynitrate crystallizes in the structure of the 3R₂ polytype while the hydroxytartarate, hydroxychloride and α-cobalt hydroxide crystallize in the structure of the 3R₁ polytype. Cobalt hydroxysulfate is turbostratically disordered. The turbostratic disorder is a direct consequence of the mismatch between the crystallographically defined interlayer sites generated within the crystal and the tetrahedral symmetry of the SO₄²⁻ ions.     

Estimation of the absolute internal-rotation entropy of molecules with two torsional degrees of freedom from stochastic simulations

A method of statistical estimation is applied to the problem of evaluating the absolute entropy of internal rotation in a molecule with two torsional degrees of freedom. The configurational part of the entropy is obtained as that of the joint probability density of an arbitrary form represented by a two-dimensional Fourier series, the coefficients of which are statistically estimated using a sample of the torsional angles of the molecule obtained by a stochastic simulation. The internal rotors in the molecule are assumed to be attached to a common frame, and their reduced moments of inertia are initially calculated as functions of the two torsional angles, but averaged over all the remaining internal degrees of freedom using the stochastic-simulation sample of the atomic configurations of the molecule. The torsional-angle dependence of the reduced moments of inertia can be also averaged out, and the absolute internal-rotation entropy of the molecule is obtained in a good approximation as the sum of the configurational entropy and a kinetic contribution fully determined by the averaged reduced moments of inertia. The method is illustrated using Monte Carlo simulations of isomers of stilbene and halogenated derivatives of propane. The two torsional angles in cis-stilbene are found to be much more strongly correlated than those in trans-stilbene, while the degree of the angular correlation in propane increases strongly on substitution of hydrogen atoms with chlorine.