CHARMM force field parameters for simulation of reactive intermediates in native and thio-substituted ribozymes.

Force field parameters specifically optimized for residues important in the study of RNA catalysis are derived from density-functional calculations, in a fashion consistent with the CHARMM27 all-atom empirical force field. Parameters are presented for residues that model reactive RNA intermediates and transition state analogs, thio-substituted phosphates and phosphoranes, and bound Mg(2+) and di-metal bridge complexes. Target data was generated via density-functional calculations at the B3LYP/6-311++G(3df,2p)// B3LYP/6-31++G(d,p) level. Partial atomic charges were initially derived from CHelpG electrostatic potential fitting and subsequently adjusted to be consistent with the CHARMM27 charges. Lennard-Jones parameters were determined to reproduce interaction energies with water molecules. Bond, angle, and torsion parameters were derived from the density-functional calculations and renormalized to maintain compatibility with the existing CHARMM27 parameters for standard residues. The extension of the CHARMM27 force field parameters for the nonstandard biological residues presented here will have considerable use in simulations of ribozymes, including the study of freeze-trapped catalytic intermediates, metal ion binding and occupation, and thio effects.

QCRNA 1.0: a database of quantum calculations for RNA catalysis.

This work outlines a new on-line database of quantum calculations for RNA catalysis (QCRNA) available via the worldwide web at http://theory.chem.umn.edu/QCRNA. The database contains high-level density functional calculations for a large range of molecules, complexes and chemical mechanisms important to phosphoryl transfer reactions and RNA catalysis. Calculations are performed using a strict, consistent protocol such that a wealth of cross-comparisons can be made to elucidate meaningful trends in biological phosphate reactivity. Currently, around 2000 molecules have been collected in varying charge states in the gas phase and in solution. Solvation was treated with both the PCM and COSMO continuum solvation models. The data can be used to study important trends in reactivity of biological phosphates, or used as benchmark data for the design of new semiempirical quantum models for hybrid quantum mechanical/molecular mechanical simulations.

The contribution of phosphate-phosphate repulsions to the free energy of DNA bending.

DNA bending is important for the packaging of genetic material, regulation of gene expression and interaction of nucleic acids with proteins. Consequently, it is of considerable interest to quantify the energetic factors that must be overcome to induce bending of DNA, such as base stacking and phosphate-phosphate repulsions. In the present work, the electrostatic contribution of phosphate-phosphate repulsions to the free energy of bending DNA is examined for 71 bp linear and bent-form model structures. The bent DNA model was based on the crystallographic structure of a full turn of DNA in a nucleosome core particle. A Green's function approach based on a linear-scaling smooth conductor-like screening model was applied to ascertain the contribution of individual phosphate-phosphate repulsions and overall electrostatic stabilization in aqueous solution. The effect of charge neutralization by site-bound ions was considered using Monte Carlo simulation to characterize the distribution of ion occupations and contribution of phosphate repulsions to the free energy of bending as a function of counterion load. The calculations predict that the phosphate-phosphate repulsions account for approximately 30% of the total free energy required to bend DNA from canonical linear B-form into the conformation found in the nucleosome core particle.

Structure and binding of Mg(II) ions and di-metal bridge complexes with biological phosphates and phosphoranes.

Divalent Mg(2+) ions often serve as cofactors in enzyme or ribozyme-catalyzed phosphoryl transfer reactions. In this work, the interaction of Mg(2+) ions and di-metal bridge complexes with phosphates, phosphoranes, and other biological ligands relevant to RNA catalysis are characterized with density functional methods. The effect of bulk solvent is treated with two continuum solvation methods (PCM and COSMO) for comparison. The relative binding affinity for different biological ligands to Mg(2+) are quantified in different protonation states. The structure and stability of the single-metal and di-metal complexes are characterized, and the changes in phosphate and phosphorane geometry induced by metal ion binding are discussed. Di-metal bridge complexes are a ubiquitous motif and the key factors governing their electrostatic stabilization are outlined. The results presented here provide quantitative characterization of metal ion binding to ligands of importance to RNA catalysis, and lay the groundwork for design of new generation quantum models that can be applied to the full biological enzymatic systems.