Computational study of the interaction of proflavine with d(ATATATATAT)2 and d(GCGCGCGCGC)2

Ultrafast spectroscopy study of DNA photophysics after proflavine intercalation

Proflavine (PF), an acridine DNA intercalating agent, has been widespread applied as an anti-microbial and topical antiseptic agent due to its ability to suppress DNA replication. On the other hand, various studies show that PF intercalation to DNA can increase photogenotoxicity and has potential chances to induce carcinomas of skin appendages. However, the effects of PF intercalation on the photophysical and photochemical properties of DNA have not been sufficiently explored. In this study, the excited state dynamics of the PF intercalated d(GC)9 • d(GC)9 and d(AT)9 • d(AT)9 DNA duplex are investigated in an aqueous buffer solution. Under 267 nm excitation, we observed ultrafast charge transfer (CT) between PF and d(GC)9 • d(GC)9 duplex, generating a CT state with an order of magnitude longer lifetime compared to that of the intrinsic excited state reported for the d(GC)9 • d(GC)9 duplex. In contrast, no excited state interaction was detected between PF and d(AT)9 • d(AT)9. Nevertheless, a localized triplet state with a lifetime over 5 µs was identified in the PF-d(AT)9 • d(AT)9 duplex.

Interaction of osmium(ii) redox probes with DNA: insights from theory

In the course of developing ultrasensitive and quantitative electrochemical point-of-care analytical tools for genetic detection of infectious diseases, osmium(ii) metallointercalators were revealed to be suitable and efficient redox probes to monitor the in vitro DNA amplification [Defever etal, Anal. Chem., 2011, 83, 1815-1821]. In this work, we thus propose a complete computational protocol in order to evaluate the affinity between Os(ii) complexes with double-stranded DNA. This protocol is based on molecular dynamics, with the parametrization of the GAFF force field for the Os(ii) complexes presenting an octahedral environment with polypyridine ligands, and QM/QM' calculations to evaluate the binding energy. For three Os(ii) probes and different binding sites, molecular dynamics simulations and interaction energies calculated at the QM/QM' level are successively discussed and compared to experimental data in order to identify the most stable binding sites. The computational protocol we propose should then be used to design more efficient Os(ii) metallointercalators.

Structure and dynamics of proflavine association around DNA

Proflavine is a small molecule that intercalates into DNA and, thereby, acts as an anticancer agent. Intercalation of proflavine is shown to be a two-step process in which the first step is believed to be the formation of a pre-intercalative outside bound state. Experimental studies so far have been unable to capture the nature of the outside bound state. However, the sub-millisecond timescale observed in fluorescence kinetic experiments is often attributed to the binding of proflavine outside of DNA. Here, we have performed molecular dynamics simulations with multiple proflavine molecules to study the structure and dynamics of the formation of the outside bound state of DNA at different ion concentrations. We observed that the timescale of the outside bound state formation is, at least, five orders of magnitude faster (in nanoseconds) than the experimentally reported timescale (sub-milliseconds) attributed to binding outside DNA. Moreover, we also observed the stacked arrangement of proflavine all around DNA, which is different from the experimentally predicted stacking arrangement perpendicular to the helical axis of DNA in the close vicinity of the phosphate groups. This study, therefore, provides insight into the molecular structure and dynamics of the pre-intercalative outside bound state and will help in understanding the overall intercalation mechanism.

B-DNA characteristics are preserved in double stranded d(A)3·d(T)3 and d(G)3·d(C)3 mini-helixes: conclusions from DFT/M06-2X study

We report the results of the first comprehensive DFT study on the d(A)3·d(T)3 and d(G)3·d(C)3 nucleic acid duplexes. The ability of mini-helixes to preserve the conformation of B-DNA in the gas phase and under the influence of such factors as: solvent, uncompensated charge, and counter-ions was evaluated using M06-2X functional with 6-31G(d,p) basis set. The accuracy of the models was ascertained based on their ability to reproduce key structural features of natural B-DNA. Analysis of the helicity suggests that the helical conformations adopt geometrical parameters which are close to those of the B-DNA form. The torsion angles fall somewhere between the values observed for BI/BII conformational classes. The comparative analysis of parameters of isolated Watson-Crick base pairs versus B-DNA-like conformations indicates the same tendency of base-pair polarization and hydration. Specifically, effects of polarization of nucleobases in continuum type dielectric medium mimicking water are stronger than those caused by the presence of backbone. Polar environment as well as the presence of counterions stabilizes duplexes, facilitating helix formation. Substantial conformational changes of nucleotides upon duplex formation decrease the binding energy. In spite of structural and energetic changes, the placement of a mini-helix into the gas phase does not lead to significant disruption of the structure. On the contrary, the duplex preserves its helicity and the strands remain bound.

Toward a rationale for the PTC124 (Ataluren) promoted readthrough of premature stop codons: a computational approach and GFP-reporter cell-based assay

The presence in the mRNA of premature stop codons (PTCs) results in protein truncation responsible for several inherited (genetic) diseases. A well-known example of these diseases is cystic fibrosis (CF), where approximately 10% (worldwide) of patients have nonsense mutations in the CF transmembrane regulator (CFTR) gene. PTC124 (3-(5-(2-fluorophenyl)-1,2,4-oxadiazol-3-yl)-benzoic acid), also known as Ataluren, is a small molecule that has been suggested to allow PTC readthrough even though its target has yet to be identified. In the lack of a general consensus about its mechanism of action, we experimentally tested the ability of PTC124 to promote the readthrough of premature termination codons by using a new reporter. The reporter vector was based on a plasmid harboring the H2B histone coding sequence fused in frame with the green fluorescent protein (GFP) cDNA, and a TGA stop codon was introduced in the H2B-GFP gene by site-directed mutagenesis. Additionally, an unprecedented computational study on the putative supramolecular interaction between PTC124 and an 11-codon (33-nucleotides) sequence corresponding to a CFTR mRNA fragment containing a central UGA nonsense mutation showed a specific interaction between PTC124 and the UGA codon. Altogether, the H2B-GFP-opal based assay and the molecular dynamics (MD) simulation support the hypothesis that PTC124 is able to promote the specific readthrough of internal TGA premature stop codons.

Proposed model for in vitro interaction between fenitrothion and DNA, by using competitive fluorescence, (31)P NMR, (1)H NMR, FT-IR, CD and molecular modeling

In this work we proposed a model for in vitro interaction of fenitrothion (FEN) with calf thymus-DNA by combination of multispectroscopic and two dimensional molecular modeling (ONIOM) methods. The circular dichroism results showed that FEN changes the conformation of B-DNA and caused some changes to C-DNA form. The FT-IR results confirmed a partial intercalation between FEN and edges of all base pairs. The competitive fluorescence, using methylene blue as fluorescence probe, in the presence of increasing amounts of FEN, revealed that FEN is able to release the non-intercalated methylene blue from the DNA. The weak chemical shift and peak broadening of (1)H NMR spectrum of FEN in the presence of DNA confirmed a non-intercalation mode. The (31)P NMR showed that FEN interacts more with DNA via its -NO2 moiety. The ONIOM, based on the hybridization of QM/MM (DFT, 6.31++G (d,p)/UFF) methodology, was also performed by Gaussian 2003 package. The results revealed that the interaction is base sequence dependent, and FEN interacts more with AT base sequences.

A three-layer ONIOM model for the outside binding of cationic porphyrins and nucleotide pair DNA

In this work we investigated the outside binding mode between a cationic porphyrin and a nucleotide pair of DNA, adenine-thymine and guanine-cytosine, in a supramolecular assembly. We used two structural models (semi-extended, extended) that differ in the size of porphyrin, two kinds of theoretical methods: a three layer ONIOM (B3LYP/6-31G(d)/PM3/UFF), and DFT B3LYP/6-31G(d,p), and three cationic porphyrins. ONIOM method was first tested on the semi-extended model that was calculated in four environments: gas phase, solution phase using an explicit solvent model (H(2)O), in the presence of a sodium cation (Na(+)) and in both (H(2)O + Na(+)). From interaction energy results, we found that the affinity of the cationic substituent by the adenine nucleotide is favored upon the thymine nucleotide. The extended model that considers the whole porphyrin was applied in the gas phase to the four nucleotides. All the cationic porphyrins showed affinity by the nucleotides in the order adenine > guanine > thymine > cytosine. The interaction energy values for outside binding showed a strong porphyrin-nucleotide interaction (≈-90 kcal mol(-1)), that slightly varies between the nucleotides suggesting that this kind of cationic porphyrin has a little selectivity for some of them. We also found that the effect of the nature of the cationic substituent (chain length) in the porphyrin on the outside binding is small (≈2-13 kcal mol(-1)). Coherence between the results showed that ONIOM is a useful tool to get a reasonable molecular geometry to be used as a starting point in calculations of density functional theory.

Electrolytes in biomolecular systems studied with the 3D-RISM/RISM theory

We reviewed our recent studies on the molecular recognition and stability of biomolecules in aqueous solutions, which have been carried out based on the statistical mechanics of molecular liquids, or the 3D-RISM/RISM theory. A special stress is put on roles of electrolytes in determining the stability of biomolecules.

The experimental and theoretical QM/MM study of interaction of chloridazon herbicide with ds-DNA

We report a multispectroscopic, voltammetric and theoretical hybrid of QM/MM study of the interaction between double-stranded DNA containing both adenine-thymine and guanine-cytosine alternating sequences and chloridazon (CHL) herbicide. The electrochemical behavior of CHL was studied by cyclic voltammetry on HMDE, and the interaction of ds-DNA with CHL was investigated by both cathodic differential pulse voltammetry (CDPV) at a hanging mercury drop electrode (HMDE) and anodic differential pulse voltammetry (ADPV) at a glassy carbon electrode (GCE). The constant bonding of CHL-DNA complex that was obtained by UV/vis, CDPV and ADPV was 2.1×10(4), 5.1×10(4) and 2.6×10(4), respectively. The competition fluorescence studies revealed that the CHL quenches the fluorescence of DNA-ethidium bromide complex significantly and the apparent Stern-Volmer quenching constant has been estimated to be 1.71×10(4). Thermal denaturation study of DNA with CHL revealed the ΔTm of 8.0±0.2°C. Thermodynamic parameters, i.e., enthalpy (ΔH), entropy (ΔS°), and Gibbs free energy (ΔG) were 98.45 kJ mol(-1), 406.3 J mol(-1) and -22.627 kJ mol(-1), respectively. The ONIOM, based on the hybridization of QM/MM (DFT, 6.31++G(d,p)/UFF) methodology, was also performed using Gaussian 2003 package. The results revealed that the interaction is base sequence dependent, and the CHL has more interaction with ds-DNA via the GC base sequence. The results revealed that CHL may have an interaction with ds-DNA via the intercalation mode.

Developing a computational model that accurately reproduces the structural features of a dinucleoside monophosphate unit within B-DNA

The ability of a dinucleoside monophosphate to mimic the conformation of B-DNA was tested using a combination of different phosphate models (anionic, neutral, counterion), environments (gas, water), and density functionals (B3LYP, MPWB1K, M06-2X) with the 6-31G(d,p) basis set. Three sequences (5'-GX(Py)-3', where X(Py) = T, U or (Br)U) were considered, which vary in the (natural or modified) 3' pyrimidine nucleobase (X(Py)). These bases were selected due to their presence in natural DNA, structural similarity to T and/or applications in radical-initiated anti-tumour therapies. The accuracy of each of the 54 model, method and sequence combinations was judged based on the ability to reproduce key structural features of natural B-DNA, including the stacked base-base orientation and important backbone torsion angles. B3LYP yields distorted or tilted relative base-base orientations, while many computational challenges were encountered for MPWB1K. Despite wide use in computational studies of DNA, the neutral (protonated) phosphate model could not consistently predict a stacked arrangement for all sequences. In contrast, stacked base-base arrangements were obtained for all sequences with M06-2X in conjunction with both the anionic and (sodium) counterion phosphate models. However, comparison of calculated and experimental backbone conformations reveals the charge-neutralized counterion phosphate model best mimics B-DNA. Structures optimized with implicit solvent (water) are comparable to gas-phase structures, which suggests similar results should be obtained in an environment of intermediate polarity. We recommend the use of either gas-phase or water M06-2X optimizations with the counterion phosphate model to study the structure and/or reactivity of natural or modified DNA with a dinucleoside monophosphate.

Synthesis, characterization and cytotoxicity studies of palladium(II)-proflavine complexes

An investigation of the reaction of Pd(II) complexes with proflavine (3,6-diaminoacridine) resulted in the isolation of the compounds [Pd(terpy)(proflavine)](NO(3))(HSO(4))*3H(2)O, 1, (terpy = 2,2':6',2″-terpyridine), [Pd(en)(proflavineH))](NO(3))(SO(4)), 2, (en = ethylenediamine), and [Pd(proflavineH)Cl(2)](SO(4))(0.5)*H(2)O, 3. They have been isolated and characterized by NMR, IR, and electro-spray ionization mass spectrometry techniques and by elemental analyses. The proflavine was bonded to the Pd(II) through the endocyclic nitrogen in 1, but through the proflavine NH(2) in 2. Compound 3 appeared to be polymeric in the solid state with a 1:1 mole ratio of Pd(II):proflavine. Upon solution of 3 in DMSO, two unique species were formed. In one species the Pd(II) was bonded to two proflavines through the endocyclic nitrogen (1:2 mole ratio) and in the other species, a Pd(II) was bonded to each NH(2) group of a single proflavine (2:1 mole ratio). Molecular modeling of the equilibrium geometry by Spartan 8 produced structures which were consistent with the experimental data on the solutions of the three compounds. In vitro cytotoxicity testing against two breast cancer cell lines and one ovarian cancer cell line showed that compounds 1 and 3 had significant activity.