Metropolis Monte Carlo calculations of DNA structure using internal coordinates and NMR distance restraints: an...

A new method, a restrained Monte Carlo (rMC) calculation, is demonstrated for generating high-resolution structures of DNA oligonucleotides in solution from interproton distance restraints and bounds derived from complete relaxation matrix analysis of two-dimensional nuclear Overhauser effect (NOE) spectral peak intensities. As in the case of restrained molecular dynamics (rMD) refinement of structures, the experimental distance restraints and bounds are incorporated as a pseudo-energy term (or penalty function) into the mathematical expression for the molecular energy. However, the use of generalized helical parameters, rather than Cartesian coordinates, to define DNA conformation increases efficiency by decreasing by an order of magnitude the number of parameters needed to describe a conformation and by simplifying the potential energy profile. The Metropolis Monte Carlo method is employed to simulate an annealing process. The rMC method was applied to experimental 2D NOE data from the octamer duplex d(GTATAATG).d(CATTATAC). Using starting structures from different locations in conformational space (e.g. A-DNA and B-DNA), the rMC calculations readily converged, with a root-mean-square deviation (RMSD) of < 0.3 A between structures generated using different protocols and starting structures. Theoretical 2D NOE peak intensities were calculated for the rMC-generated structures using the complete relaxation matrix program CORMA, enabling a comparison with experimental intensities via residual indices. Simulation of the vicinal proton coupling constants was carried out for the structures generated, enabling a comparison with the experimental deoxyribose ring coupling constants, which were not utilized in the structure determination in the case of the rMC simulations. Agreement with experimental 2D NOE and scalar coupling data was good in all cases. The rMC structures are quite similar to that refined by a traditional restrained MD approach (RMSD < 0.5 A) despite the different force fields used and despite the fact that MD refinement was conducted with additional restraints imposed on the endocyclic torsion angles of deoxyriboses. The computational time required for the rMC and rMD calculations is about the same. A comparison of structural parameters is made and some limitations of both methods are discussed with regard to the average nature of the experimental restraints used in the refinement.