Energetics and dynamics of the non-natural fluorescent 4AP:DAP base pair

Abstract

The fluorescent non-natural 4-aminophthalimide (4AP) base, when paired to the complementary 2,4-diaminopyrimidine (DAP) nucleobase, is accommodated in a B-DNA duplex being efficiently recognized and incorporated by DNA polymerases. To complement the experimental studies and rationalize the impact of the above non-natural bases on the structure, stability and dynamics of nucleic acid structures, we performed quantum mechanics (QM) calculations along with classical molecular dynamics (MD) simulations. QM calculations were initially focused on the geometry and energetics of the 4AP:DAP non-natural pair and of H-bonded base pairs between 4AP and all the natural bases in their classical Watson-Crick geometries. The QM calculations indicate that the 4AP:DAP pair, despite the fact that it can form 3 H-bonds in a classic Watson-Crick geometry, has a stability comparable to the A:T pair. Then, we extended the study to reverse Watson-Crick geometries, characteristic of parallel strands. MD simulations were carried out on two 13-mer DNA duplexes, featuring a central 4AP:DAP or A:T pair, respectively. No major structural deformation of the duplex was observed during the MD simulation. Snapshots from the MD simulations were subjected to QM calculations to investigate the 4AP:DAP interaction energy when embedded into a duplex structure, and to investigate the impact of the two non-natural bases on the stacking interactions with adjacent bases in the DNA duplex. We found a slight increase in stacking interactions involving the 4AP:DAP pair, counterbalanced by a moderate decrease in H-bonding interactions of the 4AP:DAP and of the adjacent base pairs in the duplex. The results of our study are in agreement with experimental data and complement them by providing an insight into which factors contribute positively and which factors contribute negatively to the structural compatibility of the fluorescent 4AP:DAP pair with a B-DNA structure.

Immunoinformatics Aided Design and In-Vivo Validation of a Cross-Reactive Peptide Based Multi-Epitope Vaccine Targeting Multiple Serotypes of Dengue Virus

“…Molecular dynamics (MD) simulations have been proven to be very useful and effective method in order to study the stability and analysis of biological systems (60). In fact, our research group has successfully utilized and employed MD simulations to check the stability of different protein complexes (61)(62)(63) and nucleic acid systems (64)(65)(66)(67). In this study, we used the GROMACS (45) software for molecular dynamics simulation to understand the structural properties and interaction between TLR-5 and the predicted multi-epitope vaccine.…”

The allosteric interplay between S‐nitrosylation and glycine binding controls the activity of human serine racemase

“…To perform MD simulations, the force field parameters of N'-pyridoxyl-lysine-5'-monophosphate (KPLP) were derived and implemented into the Amber-ILDN force field libraries, following previous protocols [55][56][57][58] and described below. The force constants were assigned taking into account the analogous values in standard amino acid fragments as included in the Amber-ILDN force field library.…”

Rational Design of a User-Friendly Aptamer/Peptide-Based Device for the Detection of Staphylococcus aureus

“…Counter-ions were added to neutralize the whole system. Temperature and pressure were controlled with the Berendsen algorithm [63], following previous protocols [64][65][66][67]. The Ewald method was applied to model electrostatic interactions.…”