Theoretical analysis of photosensitization of DNA by thionine

CONTEXT

In this work, we are the first to perform a theoretical analysis of photoinduced charge transfer in the intercalation complex of thionine (TH) with double-stranded DNA, which was observed in experiments. Efficient DNA binding and long-wave absorption maximum make TH an attractive photosensitizer. d(CpG)2 tetranucleotide was used as a minimal model DNA fragment. Intercalation of TH between pairs of nucleobases causes the transfer of a small negative charge (0.24 e) from the tetranucleotide to the dye. S0 → S1 photoexcitation of their complex using visible light leads to the transfer in the same direction of a significant negative charge (0.9 e). This electronic transition has a HOMO → LUMO electronic configuration, with HOMO localized on one of the two phosphate groups of the tetranucleotide, and LUMO on TH; the latter has the same shape as the LUMO of free dye. In the complex, TH, by its amino groups, forms two intermolecular H-bonds: with the deoxyribose oxygen atom of one d(CpG)2 strand and with the non-bridging oxygen atom of the phosphate group of the other strand. In this case, the H-bond TH with the phosphate group is stronger than with the sugar, but the charge transfer is carried out from another phosphate group through the sugar to the dye. Thus, charge transfer occurs along the longer of the two paths. However, the path of charge transfer depends on the parameters of the excitation since higher electronic transitions also include the second phosphate group, i.e., a short way is also used.

METHODS

For the calculations of the excitation of the complex, TD-DFT was used in combination with a set of ten functionals (CAM-B3LYP + D3BJ, ωB97XD, LC-ωHPBE, M052X, M062X, M06HF, M08HX, M11, MN15, and SOGGA11X), which have proven themselves well in modeling the excitation of dimers of aromatic molecules. Of these, LC-ωHPBE, which gave the best agreement with the experiment, was selected for the final calculations. It was used in combination with the 6-31 + + G(d,p) basis set and the IEFPCM solvent model. The photoinduced charge redistribution was quantitatively estimated using natural population analysis, and visually by building the frontier molecular and natural transition orbitals.

Nonproductive Binding Modes as a Prominent Feature of Aβ40 Fiber Elongation: Insights from Molecular Dynamics Simulation

The formation of amyloid fibers has been implicated in a number of neurodegenerative diseases. The growth of amyloid fibers is strongly thermodynamically favorable, but kinetic traps exist where the incoming monomer binds in an incompatible conformation that blocks further elongation. Unfortunately, this process is difficult to follow experimentally at the atomic level. It is also too complex to simulate in full detail and to date has been explored either through coarse-grained simulations, which may miss many important interactions, or full atomic simulations, in which the incoming peptide is constrained to be near the ideal fiber geometry. Here we use an alternate approach starting from a docked complex in which the monomer is from an experimental NMR structure of one of the major conformations in the unbound ensemble, a largely unstructured peptide with the central hydrophobic region in a 3₁₀ helix. A 1000 ns full atomic simulation in explicit solvent shows the formation of a metastable intermediate by sequential, concerted movements of both the fiber and the monomer. A Markov state model shows that the unfolded monomer is trapped at the end of the fiber in a set of interconverting antiparallel β-hairpin conformations. The simulation here may serve as a model for the binding of other non-β-sheet conformations to amyloid fibers.

Biomolecular Simulations under Realistic Macroscopic Salt Conditions

Biomolecular simulations are typically performed in an aqueous environment where the number of ions remains fixed for the duration of the simulation, generally with either a minimally neutralizing ion environment or a number of salt pairs intended to match the macroscopic salt concentration. In contrast, real biomolecules experience local ion environments where the salt concentration is dynamic and may differ from bulk. The degree of salt concentration variability and average deviation from the macroscopic concentration remains, as yet, unknown. Here, we describe the theory and implementation of a Monte Carlo osmostat that can be added to explicit solvent molecular dynamics or Monte Carlo simulations to sample from a semigrand canonical ensemble in which the number of salt pairs fluctuates dynamically during the simulation. The osmostat reproduces the correct equilibrium statistics for a simulation volume that can exchange ions with a large reservoir at a defined macroscopic salt concentration. To achieve useful Monte Carlo acceptance rates, the method makes use of nonequilibrium candidate Monte Carlo (NCMC) moves in which monovalent ions and water molecules are alchemically transmuted using short nonequilibrium trajectories, with a modified Metropolis-Hastings criterion ensuring correct equilibrium statistics for an ( Δμ, N, p, T) ensemble to achieve a ∼1046× boost in acceptance rates. We demonstrate how typical protein (DHFR and the tyrosine kinase Src) and nucleic acid (Drew-Dickerson B-DNA dodecamer) systems exhibit salt concentration distributions that significantly differ from fixed-salt bulk simulations and display fluctuations that are on the same order of magnitude as the average.

Molecular dynamics simulation of the opposite-base preference and interactions in the active site of formamidopyrimidine-DNA glycosylase

Background

Formamidopyrimidine-DNA glycosylase (Fpg) removes abundant pre-mutagenic 8-oxoguanine (oxoG) bases from DNA through nucleophilic attack of its N-terminal proline at C1′ of the damaged nucleotide. Since oxoG efficiently pairs with both C and A, Fpg must excise oxoG from pairs with C but not with A, otherwise a mutation occurs. The crystal structures of several Fpg–DNA complexes have been solved, yet no structure with A opposite the lesion is available.

Results

Here we use molecular dynamic simulation to model interactions in the pre-catalytic complex of Lactococcus lactis Fpg with DNA containing oxoG opposite C or A, the latter in either syn or anti conformation. The catalytic dyad, Pro1–Glu2, was modeled in all four possible protonation states. Only one transition was observed in the experimental reaction rate pH dependence plots, and Glu2 kept the same set of interactions regardless of its protonation state, suggesting that it does not limit the reaction rate. The adenine base opposite oxoG was highly distorting for the adjacent nucleotides: in the more stable syn models it formed non-canonical bonds with out-of-register nucleotides in both the damaged and the complementary strand, whereas in the anti models the adenine either formed non-canonical bonds or was expelled into the major groove. The side chains of Arg109 and Phe111 that Fpg inserts into DNA to maintain its kinked conformation tended to withdraw from their positions if A was opposite to the lesion. The region showing the largest differences in the dynamics between oxoG:C and oxoG:A substrates was unexpectedly remote from the active site, located near the linker joining the two domains of Fpg. This region was also highly conserved among 124 analyzed Fpg sequences. Three sites trapping water molecules through multiple bonds were identified on the protein–DNA interface, apparently helping to maintain enzyme-induced DNA distortion and participating in oxoG recognition.

Conclusion

Overall, the discrimination against A opposite to the lesion seems to be due to incorrect DNA distortion around the lesion-containing base pair and, possibly, to gross movement of protein domains connected by the linker.

Electronic supplementary material

The online version of this article (doi:10.1186/s12900-017-0075-y) contains supplementary material, which is available to authorized users.

Design and characterization of symmetric nucleic acids via molecular dynamics simulations

Asymmetry (5'→3') associated with each strand of the deoxyribonucleic acid (DNA) is inherent in the sugar-phosphate backbone connectivity and is essential for replication and transcription. We note that this asymmetry is due to one single chemical bond (C_3' to C_2' ) in each nucleotide unit, and the absence of this bond results in directionally symmetric nucleic acids. We also discovered that creation of an extra chemical bond (C_5' to C_2' ) can lead to a symmetric backbone. Keeping their potential synthetic and therapeutic interest in mind, we designed a few novel symmetric nucleic acids. We investigated their conformational stability and flexibility via detailed all atom explicit solvent 100-ns long molecular dynamics simulations and compared the resulting structures with that of regular B-DNA. Quite interestingly, some of the symmetric nucleic acids retain the overall double helical structure indicating their potential for integration in physiological DNA without causing major structural perturbations.

8-Oxoguanine Affects DNA Backbone Conformation in the EcoRI Recognition Site and Inhibits Its Cleavage by the Enzyme

8-oxoguanine is one of the most abundant and impactful oxidative DNA lesions. However, the reasons underlying its effects, especially those not directly explained by the altered base pairing ability, are poorly understood. We report the effect of the lesion on the action of EcoRI, a widely used restriction endonuclease. Introduction of 8-oxoguanine inside, or adjacent to, the GAATTC recognition site embedded within the Drew—Dickerson dodecamer sequence notably reduced the EcoRI activity. Solution NMR revealed that 8-oxoguanine in the DNA duplex causes substantial alterations in the sugar—phosphate backbone conformation, inducing a BI→BII transition. Moreover, molecular dynamics of the complex suggested that 8-oxoguanine, although does not disrupt the sequence-specific contacts formed by the enzyme with DNA, shifts the distribution of BI/BII backbone conformers. Based on our data, we propose that the disruption of enzymatic cleavage can be linked with the altered backbone conformation and dynamics in the free oxidized DNA substrate and, possibly, at the protein—DNA interface.

Nucleic acids: theory and computer simulation, Y2K

Molecular dynamics simulations on DNA and RNA that include solvent are now being performed under realistic environmental conditions of water activity and salt. Improvements to force-fields and treatments of long-range interactions have significantly increased the reliability of simulations. New studies of sequence effects, axis bending, solvation and conformational transitions have appeared.

Molecular dynamic simulations of environment and sequence dependent DNA conformations: the development of the BMS nucleic acid force field and comparison with experimental results

Molecular dynamic (MD) simulations using the BMS nucleic acid force field produce environment and sequence dependent DNA conformations that closely mimic experimentally derived structures. The parameters were initially developed to reproduce the potential energy surface, as defined by quantum mechanics, for a set of small molecules that can be used as the building blocks for nucleic acid macromolecules (dimethyl phosphate, cyclopentane, tetrahydrofuran, etc.). Then the dihedral parameters were fine tuned using a series of condensed phase MD simulations of DNA and RNA (in zero added salt, 4M NaCl, and 75% ethanol solutions). In the tuning process the free energy surface for each dihedral was derived from the MD ensemble and fitted to the conformational distributions and populations observed in 87 A- and B-DNA x-ray and 17 B-DNA NMR structures. Over 41 nanoseconds of MD simulations are presented which demonstrate that the force field is capable of producing stable trajectories, in the correct environments, of A-DNA, double stranded A-form RNA, B-DNA, Z-DNA, and a netropsin-DNA complex that closely reproduce the experimentally determined and/or canonical DNA conformations. Frequently the MD averaged structure is closer to the experimentally determined structure than to the canonical DNA conformation. MD simulations of A- to B- and B- to A-DNA transitions are also shown. A-DNA simulations in a low salt environment cleanly convert into the B-DNA conformation and converge into the RMS space sampled by a low salt simulation of the same sequence starting from B-DNA. In MD simulations using the BMS force field the B-form of d(GGGCCC)2 in a 75% ethanol solution converts into the A-form. Using the same methodology, parameters, and conditions the A-form of d(AAATTT)2 correctly converts into the B-DNA conformation. These studies demonstrate that the force field is capable of reproducing both environment and sequence dependent DNA structures. The 41 nanoseconds (nsec) of MD simulations presented in this paper paint a global picture which suggests that the DNA structures observed in low salt solutions are largely due to the favorable internal energy brought about by the nearly uniform screening of the DNA electrostatics. While the conformations sampled in high salt or mixed solvent environments occur from selective and asymmetric screening of the phosphate groups and DNA grooves, respectively, brought about by sequence induced ion and solvent packing.

Structural equilibrium of DNA represented with different force fields

We have recently indicated preliminary evidence of different equilibrium average structures with the CHARMM and AMBER force fields in explicit solvent molecular dynamics simulations on the DNA duplex d(C5T5) . d(A5G5) (Feig, M. and B.M. Pettitt, 1997, Experiment vs. Force Fields: DNA conformation from molecular dynamics simulations. J. Phys. Chem. B. (101:7361-7363). This paper presents a detailed comparison of DNA structure and dynamics for both force fields from extended simulation times of 10 ns each. Average structures display an A-DNA base geometry with the CHARMM force field and a base geometry that is intermediate between A- and B-DNA with the AMBER force field. The backbone assumes B form on both strands with the AMBER force field, while the CHARMM force field produces heterogeneous structures with the purine strand in A form and the pyrimidine strand in dynamical equilibrium between A and B conformations. The results compare well with experimental data for the cytosine/guanine part but fail to fully reproduce an overall B conformation in the thymine/adenine tract expected from crystallographic data, particularly with the CHARMM force field. Fluctuations between A and B conformations are observed on the nanosecond time scale in both simulations, particularly with the AMBER force field. Different dynamical behavior during the first 4 ns indicates that convergence times of several nanoseconds are necessary to fully establish a dynamical equilibrium in all structural quantities on the time scale of the simulations presented here.