Femtosecond direct observation of charge transfer between bases in DNA

Specificity of DNA damage formation induced by femtosecond near-infrared laser filamentation in water

We investigated the deoxyribonucleic acid (DNA) damage induced by laser filamentation, which was generated by focusing femtosecond near-infrared Ti:Sapphire laser light in water at several repetition rates ranging from 1000 Hz to 10 Hz. Using plasmid DNA (pUC19), the single-strand break, double-strand break, nucleobase lesions, and the fragmented DNA were analyzed and quantified by agarose gel electrophoresis. Additionally, the H2O2 concentration after irradiation was determined. We observed that (1) the DNA damage per laser shot and (2) the enzyme-sensitive base lesions per total DNA damage decreased as the laser repetition rate increased. Furthermore, (3) extraordinarily short DNA fragments were likely to be produced, compared with those produced using X-rays, and (4) most OH radicals could be eliminated by recombination to generate H2O2, preventing them from damaging the DNA. The Monte-Carlo simulation of the strand break formation implies that the observed dependency of strand break efficiency on the laser repetition rate is mainly due to diffusion of DNA molecules. These findings quantitatively and qualitatively revealed that an intense laser pulse induces a specific DNA damage profile that is not induced by X-rays, a sparsely ionizing radiation source.

Electronic structure, absorption spectra and oxidation dynamics in polyynes and dicyanopolyynes

The advent of femtosecond to attosecond experimental tools has made now possible to study such ultrafast carrier dynamics, e.g., the spatial and temporal charge density evolution, after an initial oxidation or reduction in molecules, candidates for atomic wires like polyynes and dicyanopolyynes. Here, we study the electronic structure and hole transfer in symmetric molecules containing carbon, nitrogen and hydrogen, the first members in the series of polyynic carbynes and dicyanopolyynes, using methods based on density functional theory (DFT): constrained DFT (CDFT), time-dependent DFT (TDDFT) and real-time TDDFT (RT-TDDFT), with Löwdin population analysis, comparing many levels of theory and obtaining convergence of the results. For the same purposes, we develop a tight binding (TB) variant using all valence orbitals of all atoms. This TB variant is applied here in linear molecules, but it is also adequate for electronic structure, charge transfer and charge transport of non-linear molecules and clusters of molecules. We calculate the electronic structure, the time-dependent dipole moment and the probabilities of finding the hole at each site, their mean over time values, the mean transfer rates from the oxidation site to other sites and the frequency content (using charge as well as dipole moment oscillations). We take into account zero-point motion. The initial conditions for RT-TDDFT are obtained by CDFT. For TB, we explore different initial conditions: we place the hole at a particular orbital or distribute it among a number of orbitals; it is also possible to include phase differences between orbitals. Finally, we compare with available experimental data.

2-Aminopurine-based quencher-free DNA tweezers with fluorescence properties well tuned by surrounding bases

Reversible structural changes in DNA nanomachines have great potential in the field of bioanalysis. Here, we demonstrate an assembly strategy for quencher-free and tunable DNA tweezers based on 2-aminopurine (2-AP), avoiding the tedious fluorescence labelling step. The conformational state of the tweezers could be controlled by specific oligonucleotides (fuel or anti-fuel). Taking advantage of the local environmental sensitivity of 2-AP, the structural changes of the tweezers were easily tracked, and multiple cyclic switching of the tweezers between the open and closed states was achieved. In addition, the influence of oligonucleotide structure on the fluorescence properties of 2-AP was deeply explored. We figured out that the fluorescence of 2-AP was highly quenched by the base-stacking of natural bases in DNA oligonucleotides. Moreover, by comprehensively regulating the type of bases surrounding the inserted 2-AP site, a sensitive fluorescence response towards dynamic change can be obtained. This principle of quencher-free nanodevices based on 2-AP provides a convenient method for monitoring the structural changes of DNA nanomachines.

Geometry-structure models for liquid crystal interfaces, drops and membranes: wrinkling, shape selection and dissipative shape evolution

We review our recent contributions to anisotropic soft matter models for liquid crystal interfaces, drops and membranes, emphasizing validations with experimental and biological data, and with related theory and simulation literature. The presentation aims to illustrate and characterize the rich output and future opportunities of using a methodology based on the liquid crystal-membrane shape equation applied to static and dynamic pattern formation phenomena. The geometry of static and kinetic shapes is usually described with dimensional curvatures that co-mingle shape and curvedness. In this review, we systematically show how the application of a novel decoupled shape-curvedness framework to practical and ubiquitous soft matter phenomena, such as the shape of drops and tactoids and bending of evolving membranes, leads to deeper quantitative insights than when using traditional dimensional mean and Gaussian curvatures. The review focuses only on (1) statics of wrinkling and shape selection in liquid crystal interfaces and membranes; (2) kinetics and dissipative dynamics of shape evolution in membranes; and (3) computational methods for shape selection and shape evolution; due to various limitations other important topics are excluded. Finally, the outlook follows a similar structure. The main results include: (1) single and multiple wavelength corrugations in liquid crystal interfaces appear naturally in the presence of surface splay and bend orientation distortions with scaling laws governed by ratios of anchoring-to-isotropic tension energy; adding membrane elasticity to liquid crystal anchoring generates multiple scales wrinkling as in tulips; drops of liquid crystals encapsulates in membranes can adopt, according to the ratios of anchoring/tension/bending, families of shapes as multilobal, tactoidal, and serrated as observed in biological cells. (2) Mapping the liquid crystal director to a membrane unit normal. The dissipative shape evolution model with irreversible thermodynamics for flows dominated by bending rates, yields new insights. The model explains the kinetic stability of cylinders, while spheres and saddles are attractors. The model also adds to the evolving understanding of outer hair cells in the inner ear. (3) Computational soft matter geometry includes solving shape equations, trajectories on energy and orientation landscapes, and shape-curvedness evolutions on entropy production landscape with efficient numerical methods and adaptive approaches.

Theoretical Studies for the Discovery of Potential Sucrase-Isomaltase Inhibitors from Maize Silk Phytochemicals: An Approach to Treatment of Type 2 Diabetes

A theoretical analysis of the potential inhibition of human sucrase-isomaltase (SI) by flavonoids was carried out with the aim of identifying potential candidates for an alternative treatment of type 2 diabetes. Two compounds from maize silks, maysin and luteolin, were selected to be studied with the structure-based density functional theory (DFT), molecular docking (MDock), and molecular dynamics (MD) approaches. The docking score and MD simulations suggested that the compounds maysin and luteolin presented higher binding affinities in N-terminal sucrase-isomaltase (NtSI) than in C-terminal sucrase-isomaltase (CtSI). The reactivity parameters, such as chemical hardness (η) and chemical potential (µ), of the ligands, as well as of the active site amino acids of the NtSI, were calculated by the meta-GGA M06 functional in combination with the 6-31G(d) basis set. The lower value of chemical hardness calculated for the maysin molecule indicated that this might interact more easily with the active site of NtSI, in comparison with the values of the acarbose and luteolin structures. Additionally, a possible oxidative process was proposed through the quantum chemical calculations of the electronic charge transfer values (∆N) between the active site amino acids of the NtSI and the ligands. In addition, maysin displayed a higher ability to generate more oxidative damage in the NtSI active site. Our results suggest that maysin and luteolin can be used to develop novel α-glucosidase inhibitors via NtSI inhibition.

Screening for Novel Fluorescent Nucleobase Analogues Using Computational and Experimental Methods: 2-Amino-6-chloro-8-vinylpurine (2A6Cl8VP) as a Case Study

Novel fluorescent nucleic acid base analogues (FBAs) with improved optical properties are needed in a variety of biological applications. 2-Amino-6-chloro-8-vinylpurine (2A6Cl8VP) is structural analogue of two existing highly fluorescent FBAs, 2-aminopurine (2AP) and 8-vinyladenine (8VA), and can therefore be expected to have similar base pairing as well as better optical properties compared to its counterparts. In order to determine the absorption and fluorescence properties of 2A6Cl8VP, as a first step, we used TD-DFT calculations and the polarizable continuum model for simulating the solvents and computationally predicted absorption and fluorescence maxima. To test the computational predictions, we also synthesized 2A6Cl8VP and measured its UV/vis absorbance, fluorescence emission, and fluorescence lifetime. The computationally predicted absorbance and fluorescence maxima of 2A6Cl8VP are in reasonable agreement to the experimental values and are significantly redshifted compared to 2AP and 8VA, allowing for its specific excitation. The fluorescence quantum yield of 2A6Cl8VP, however, is significantly lower than those of 2AP and 8VA. Overall, 2A6Cl8VP is a novel fluorescent nucleobase analogue, which can be useful in studying structural, biophysical, and biochemical applications.

Reaction mechanism - explored with the unified reaction valley approach

One of the ultimate goals of chemistry is to understand and manipulate chemical reactions, which implies the ability to monitor the reaction and its underlying mechanism at an atomic scale. In this article, we introduce the Unified Reaction Valley Approach (URVA) as a tool for elucidating reaction mechanisms, complementing existing computational procedures. URVA combines the concept of the potential energy surface with vibrational spectroscopy and describes a chemical reaction via the reaction path and the surrounding reaction valley traced out by the reacting species on the potential energy surface on their way from the entrance to the exit channel, where the products are located. The key feature of URVA is the focus on the curving of the reaction path. Moving along the reaction path, any electronic structure change of the reacting species is registered by a change in the normal vibrational modes spanning the reaction valley and their coupling with the path, which recovers the curvature of the reaction path. This leads to a unique curvature profile for each chemical reaction, with curvature minima reflecting minimal change and curvature maxima indicating the location of important chemical events such as bond breaking/formation, charge polarization and transfer, rehybridization, etc. A decomposition of the path curvature into internal coordinate components or other coordinates of relevance for the reaction under consideration, provides comprehensive insight into the origin of the chemical changes taking place. After giving an overview of current experimental and computational efforts to gain insight into the mechanism of a chemical reaction and presenting the theoretical background of URVA, we illustrate how URVA works for three diverse processes, (i) [1,3] hydrogen transfer reactions; (ii) α-keto-amino inhibitor for SARS-CoV-2 Mpro; (iii) Rh-catalyzed cyanation. We hope that this article will inspire our computational colleagues to add URVA to their repertoire and will serve as an incubator for new reaction mechanisms to be studied in collaboration with our experimental experts in the field.

Programmed Control of Fluorescence Blinking Patterns based on Electron Transfer in DNA

Fluorescence imaging uses changes in the fluorescence intensity and emission wavelength to analyze multiple targets simultaneously. To increase the number of targets that can be identified simultaneously, fluorescence blinking can be used as an additional parameter. To understand and eventually control blinking, we used DNA as a platform to elucidate the processes of electron transfer (ET) leading to blinking, down to the rate constants. With a fixed ET distance, various blinking patterns were observed depending on the DNA sequence between the donor and acceptor units of the DNA platform. The blinking pattern was successfully described with a combination of ET rate constants. Therefore, molecules with various blinking patterns can be developed by tuning ET. It is expected that the number of targets that can be analyzed simultaneously will increase by the power of the number of blinking patterns.

Base-Stacking Heterogeneity in RNA Resolved by Fluorescence-Detected Circular Dichroism Spectroscopy

RNA plays a critical role in many biological processes, and the structures it adopts are intimately linked to those functions. Among many factors that contribute to RNA folding, van der Waals interactions between adjacent nucleobases stabilize structures in which the bases are stacked on top of one another. Here, we utilize fluorescence-detected circular dichroism spectroscopy (FDCD) to investigate base-stacking heterogeneity in RNA labeled with the fluorescent adenine analogue 2-aminopurine (2-AP). Comparison of standard (transmission-detected) CD and FDCD spectra reveals that in dinucleotides, 2-AP fluorescence is emitted almost exclusively by unstacked molecules. In a trinucleotide, some fluorescence is emitted by a population of stacked and highly quenched molecules, but more than half originates from a minor ∼10% population of unstacked molecules. The combination of FDCD and standard CD measurements reveals the prevalence of stacked and unstacked conformational subpopulations as well as their relative fluorescence quantum yields.

Nonstatistical Photoinduced Processes in Gaseous Organic Molecules

Processes that proceed in femtoseconds are usually referred to as being ultrafast, and they are investigated in experiments that involve laser pulses with femtosecond duration in so-called pump probe schemes, where a light pulse triggers a molecular process and a second light pulse interrogates the temporal evolution of the molecular population. The focus of this review is on the reactivity patterns that arise when energy is not equally distributed on all the available degrees of freedom as a consequence of the very short time scale in play and on how the localization of internal energy in a specific mode can be thought of as directing a process toward (or away from) a certain outcome. The nonstatistical aspects are illustrated with examples from photophysics and photochemistry for a range of organic molecules. The processes are initiated by a variety of nuclear motions that are all governed by the energy gradients in the Franck-Condon region. Essentially, the molecules will start to adapt to the new electronic environment on the excited state to eventually reach the equilibrium structure. It is this structural change that is enabling an ultrafast electronic transition in cases where the nuclear motion leads to a transition point with significant coupling between to electronic states and to ultrafast reaction if there is a coupling to a reactive mode at the transition point between the involved states. With the knowledge of the relation between electronic excitation and equilibrium structure, it is possible to predict how the nuclei move after excitation and often whether an ultrafast (and inherently nonstatistical) electronic transition or even a bond breakage will take place. In addition to the understanding of how nonstatistical photoinduced processes proceed from a given excited state, it has been found that randomization of the energy does not even always take place when the molecule takes part in processes that are normally considered statistical, such as for example nonradiative transitions between excited states. This means that energy can be localized in a specific degree of freedom on a state other than the one that is initially prepared. This is a finding that could kickoff the ultimate dream in applied photochemistry; namely light excitation that leads to the rupture of a specific bond.

The mystery of sub-picosecond charge transfer following irradiation of hydrated uridine monophosphate

The early mechanisms by which ionizing rays damage biological structures by so-called direct effects are largely elusive. In a recent picosecond pulse radiolysis study of concentrated uridine monophosphate solutions [J. Ma, S. A. Denisov, J.-L. Marignier, P. Pernot, A. Adhikary, S. Seki and M. Mostafavi, J. Phys. Chem. Lett., 2018, 9, 5105], unexpected results were found regarding the oxidation of the nucleobase. The signature of the oxidized nucleobase could not be detected 5 ps after the electron pulse, but only the oxidized phosphate, raising intriguing questions about the identity of charge-transfer mechanisms that could explain the absence of U⁺. We address here this question by means of advanced first-principles atomistic simulations of solvated uridine monophosphate, combining Density Functional Theory (DFT) with polarizable embedding schemes. We contrast three very distinct mechanisms of charge transfer covering the atto-, femto- and pico-second timescales. We first investigate the ionization mechanism and subsequent hole/charge migrations on a timescale of attoseconds to a few femtoseconds under the frozen nuclei approximation. We then consider a nuclear-driven phosphate-to-oxidized-nucleobase electron transfer, showing that it is an uncompetitive reaction channel on the sub-picosecond timescale, despite its high exothermicity and significant electronic coupling. Finally, we show that non-adiabatic charge transfer is enabled by femtosecond nuclear relaxation after ionization. We show that electronic decoherence and the electronic coupling strength are the key parameters that determine the hopping probabilities. Our results provide important insight into the interplay between electronics and nuclear motions in the early stages of the multiscale responses of biological matter subjected to ionizing radiation.

Transition between Random and Periodic Electron Currents on a DNA Chain

By resorting to a model inspired to the standard Davydov and Holstein-Fröhlich models, in the present paper we study the motion of an electron along a chain of heavy particles modeling a sequence of nucleotides proper to a DNA fragment. Starting with a model Hamiltonian written in second quantization, we use the Time Dependent Variational Principle to work out the dynamical equations of the system. It can be found that, under the action of an external source of energy transferred to the electron, and according to the excitation site, the electron current can display either a broad frequency spectrum or a sharply peaked frequency spectrum. This sequence-dependent charge transfer phenomenology is suggestive of a potentially rich variety of electrodynamic interactions of DNA molecules under the action of electron excitation. This could imply the activation of interactions between DNA and transcription factors, or between DNA and external electromagnetic fields.

Dependence of Fluorescence Quenching of CY3 Oligonucleotide Conjugates on the Oxidation Potential of the Stacking Base Pair

To understand the complex fluorescence properties of astraphloxin (CY3)-labelled oligonucleotides, it is necessary to take into account the redox properties of the nucleobases. In oligonucleotide hybrids, we observed a dependence of the fluorescence intensity on the oxidation potential of the neighbouring base pair. For the series I < A < G < 8-oxoG, the extent of fluorescence quenching follows the trend of decreasing oxidation potentials. In a series of 7 nt hybrids, stacking interactions of CY3 with perfect match and mismatch base pairs were found to stabilise the hybrid by 7–8 kJ/mol. The fluorescence measurements can be explained by complex formation resulting in fluorescence quenching that prevails over the steric effect of a reduced excited state trans-cis isomerisation, which was expected to increase the fluorescence efficiency of the dye when stacking to a base pair. This can be explained by the fact that, in a double strand, base pairing and stacking cause a dramatic change in the oxidation potential of the nucleobases. In single-molecule fluorescence measurements, the oxidation of G to 8-oxoG was observed as a result of photoinduced electron transfer and subsequent chemical reactions. Our results demonstrate that covalently linked CY3 is a potent oxidant towards dsDNA. Sulfonated derivatives should be used instead.

Hole Transfer in Open Carbynes

We investigate hole transfer in open carbynes, i.e., carbon atomic nanowires, using Real-Time Time-Dependent Density Functional Theory (RT-TDDFT). The nanowire is made of N carbon atoms. We use the functional B3LYP and the basis sets 3-21G, 6-31G*, cc-pVDZ, cc-pVTZ, cc-pVQZ. We also utilize a few Tight-Binding (TB) wire models, a very simple model with all sites equivalent and transfer integrals given by the Harrison ppπ expression (TBI) as well as a model with modified initial and final sites (TBImod) to take into account the presence of one or two or three hydrogen atoms at the edge sites. To achieve similar site occupations in cumulenes with those obtained by converged RT-TDDFT, TBImod is sufficient. However, to achieve similar frequency content of charge and dipole moment oscillations and similar coherent transfer rates, the TBImod transfer integrals have to be multiplied by a factor of four (TBImodt4times). An explanation for this is given. Full geometry optimization at the B3LYP/6-31G* level of theory shows that in cumulenes bond length alternation (BLA) is not strictly zero and is not constant, although it is symmetrical relative to the molecule center. BLA in cumulenic cases is much smaller than in polyynic cases, so, although not strictly, the separation to cumulenes and polyynes, approximately, holds. Vibrational analysis confirms that for N even all cumulenes with coplanar methylene end groups are stable, for N odd all cumulenes with perpendicular methylene end groups are stable, and the number of hydrogen atoms at the end groups is clearly seen in all cumulenic and polyynic cases. We calculate and discuss the Density Functional Theory (DFT) ground state energy of neutral molecules, the CDFT (Constrained DFT) "ground state energy" of molecules with a hole at one end group, energy spectra, density of states, energy gap, charge and dipole moment oscillations, mean over time probabilities to find the hole at each site, coherent transfer rates, and frequency content, in general. We also compare RT-TDDFT with TB results.

About electron transfer over long distances with tunable sub/supersonic velocities

Provided in this paper is a theory of long-range electron transfer with near sound (supersonic or subsonic) velocity along one-dimensional crystal lattices. The theory represents the development of an earlier work by introducing Marcus formulation. To illustrate its application to a realistic case, the theory is used to offer an explanation of two puzzling observations made by Donovan and Wilson in transient photoconduction experiments with non-dopable perfectly crystalline polydiacetylene crystals in the presence of an electric field: transport velocity value close to sound velocity being independent of field for four orders of magnitude of field (10² V/m-10⁶ V/m) and, in the low field values, an ultra-high mobility greater than 20 m²/V s. We also study factors eventually leading to lowering of the transport velocity.

Dynamic and static quenching of 2-aminopurine fluorescence by the natural DNA nucleotides in solution

2-aminopurine (2AP) is a responsive fluorescent base analogue that is used widely as a probe of the local molecular environment in DNA. The ability of 2AP to report changes in local conformation and base-stacking interactions arises from the efficient quenching of its fluorescence by the natural DNA bases. However, the mechanism of this inter-base quenching remains imperfectly understood. Two previous studies of the collisional quenching of 2AP by the natural bases, in different buffer solutions, showed that dynamic quenching efficiency depends on the identity of the natural base, but disagreed on the relative quenching efficiencies of the bases. We report a comprehensive investigation of inter-base quenching of 2AP by the natural nucleoside monophosphates (NMPs), replicating the buffer conditions used in the previous studies. Using time-resolved fluorescence measurements to distinguish between dynamic and static quenching, we find that the dynamic quenching rate constants of the different bases show a consistent trend across both buffers, and this is in line with a charge-transfer mechanism. Time-resolved measurements also provide insight into static quenching, revealing formation of 2AP-NMP ground-state complexes in which 2AP displays a very short fluorescence lifetime, comparable to that seen in oligonucleotides. In these complexes, the dependence of the rate of quenching on the partner base also supports a charge-transfer mechanism.

Polyethylene glycols affect electron transfer rate in phenosafranin-DNA complex

Long distance electron transfer (ET) between small ligands and DNA is a much studied phenomenon and is principally believed to occur through electron (or hole) hopping. Several studies have been carried out in aqueous environments while in real biological milieu the DNA molecules experience a more dense and heterogeneous environment containing otherwise indifferent molecular crowders. It is therefore expected that the ET could get modified in the presence of crowding agent and to investigate that we have made elaborate studies on steady state and time-resolved (picosecond (ps) and femtosecond (fs)-resolved) emission properties of a phenosafranine (PSF) intercalated to calf thymus (CT) DNA in the presence of ethylene glycol (EG) and polyethylene glycols (PEG) of different chain lengths (PEG 200, 400 and 1000). The emission of PSF gets considerably quenched when intercalated to DNA; the quenching is released when PEGs are added into it. The structural integrity of the CT DNA has been established using circular dichroism spectroscopy. CD measurements have evidenced only marginal changes in the DNA structure upon the addition of PEGs. ps-Resolved fluorescence measurements show significant decrease in the contribution of the DNA induced quenched time-constant of PSF upon the addition of PEGs, however, fs-resolved measurements show less noticeable changes in the time constants. Our study shows that the electron hopping rate through the guanine base in DNA core remains unaffected whereas the 'through space' electron transfer process does get affected in the presence of molecular crowders.

Electron Transfer in DNA at Electrified Interfaces

The ability of the DNA double helix to transport electrons underlies many life-centered biological processes and bio-electronic applications. However, there is little consensus on how efficiently the base pair π-stacks of DNA mediate electron transport. This minireview scrutinizes the current state-of-the-art knowledge on electron transfer (ET) properties of DNA and its long-range ability to transfer (mediate) electrical signals at electrified interfaces, without being oxidized or reduced. Complex changes an electric field induces in the DNA structure and its electronic properties govern the efficiency of DNA-mediated ET at electrodes and allow addressing the existing phenomenological riddles, while recently discovered rectifying properties of DNA contribute both to our understanding of DNA's ET in living systems and to advances in molecular bioelectronics.

Quasi-Periodic and Fractal Polymers: Energy Structure and Carrier Transfer

We study the energy structure and the coherent transfer of an extra electron or hole along aperiodic polymers made of N monomers, with fixed boundaries, using B-DNA as our prototype system. We use a Tight-Binding wire model, where a site is a monomer (e.g., in DNA, a base pair). We consider quasi-periodic (Fibonacci, Thue-Morse, Double-Period, Rudin-Shapiro) and fractal (Cantor Set, Asymmetric Cantor Set) polymers made of the same monomer (I polymers) or made of different monomers (D polymers). For all types of such polymers, we calculate the highest occupied molecular orbital (HOMO) eigenspectrum and the lowest unoccupied molecular orbital (LUMO) eigenspectrum, the HOMO-LUMO gap and the density of states. We examine the mean over time probability to find the carrier at each monomer, the frequency content of carrier transfer (Fourier spectra, weighted mean frequency of each monomer, total weighted mean frequency of the polymer), and the pure mean transfer rate k. Our results reveal that there is a correspondence between the degree of structural complexity and the transfer properties. I polymers are more favorable for charge transfer than D polymers. We compare k ( N ) of quasi-periodic and fractal sequences with that of periodic sequences (including homopolymers) as well as with randomly shuffled sequences. Finally, we discuss aspects of experimental results on charge transfer rates in DNA with respect to our coherent pure mean transfer rates.

Remote control of charge transport and chiral induction along a DNA-metallohelicate

Herein we present a new strategy to achieve chiral induction and redox switching along the backbone of metallohelicate architectures, wherein a DNA duplex directs the handedness and charge transport properties of a metal-organic assembly more than 60 bonds away (a distance of >10 nm). The quantitative and site-specific binding of copper(i) ions to DNA-templated coordination sites imparts enhanced thermodynamic stability to the assembly, while the DNA duplex transfers its natural right-handed helicity to the proximal and distal metal centers of the helicates. When copper(ii) ions are employed instead of copper(i) ions, spontaneous DNA-mediated reduction occurs, which we propose is followed by a slower change in coordination environment (from pentacoordinate CuII to tetrahedral CuI) to generate copper(i) helicates. We demonstrate that the reduction of the adjacent and distal bis-phenanthroline sites is dependent on their proximity to DNA guanine bases (which act as the electron source). The kinetics of helical charge transport can thus be tuned based on guanine-CuII separation, resulting in a sequence- and distance-dependent redox switch that transfers electronic information from DNA to multiple linearly-arranged metal centers.

Theoretical modeling of DNA electron hole transport through polypyrimidine sequences: a QM/MM study

The phenomenon of DNA hole transport (HT) has attracted of scientists for several decades, mainly due to its potential application in molecular electronics. As electron holes mostly localize on purine bases in DNA, the majority of scientific effort has been invested into chemically modifying the structures of adenine and guanine in order to increase their HT-mediating properties. In this work we examine an alternative, never yet explored, way of affecting the HT efficiency by forcing electron holes to localize on pyrimidine bases and move between them. Using an enhanced and revised version of our previously developed QM/MM model, we perform simulations of HT through polyadenine, polycytosine, polyguanine, and polythymine stacks according to a multistep hopping mechanism. From these simulations, kinetic parameters for HT are obtained. The results indicate a particularly high efficiency of cytosine→cytosine hopping, which is about ten times higher than the G → G hopping. We also discuss possible improvement of cytosine HT by modifying the oxidoreductive properties of complementary guanine residues.

Genotoxic and cytotoxic effects of 50 Hz 1 mT electromagnetic field on larval rainbow trout (Oncorhynchus mykiss), Baltic clam (Limecola balthica) and common ragworm (Hediste diversicolor)

The aim of this research was to assess genotoxicity and cytotoxicity responses in aquatic animals exposed to 50 Hz 1 m T electromagnetic field (EMF). Rainbow trout (Oncorhynchus mykiss) at early stages of development were exposed to EMF for 40 days, whereas marine benthic invertebrates - the common ragworm Hediste diversicolor and the Baltic clam Limecola balthica - for 12 days. To define genotoxicity and cytotoxicity responses in selected animals, assays of nuclear abnormalities in peripheral blood erythrocytes of O. mykiss, coelomocytes of H. diversicolor and gill cells of L. balthica were performed. Induction of formation of micronuclei (MN), nuclear buds (NB), nuclear buds on filament cells (NBf) and cells with blebbed nuclei (BL) were assessed as genotoxicity endpoints, and 8-shaped nuclei, fragmented (Fr), apoptotic (Ap) and binucleated (BN) cells as cytotoxicity endpoints. Exposure to EMF affected all studied species but with varying degrees. The strongest responses to EMF treatment were elicited in L. balthica, in which six out of the total eight analyzed geno- and cytotoxicity endpoints were significantly elevated. Significantly induced frequencies of MN were detected in O. mykiss and H. diversicolor cells, NBf and BL only in gill cells of L. balthica, and NB in analyzed tissues of all the test species. As cytotoxicity endpoints, a significant elevation in frequencies of cells with 8-shaped nuclei was found in O. mykiss and L. balthica, while Ap and BN was observed only in L. balthica. EMF exposure did not induce any significant cytotoxic activity in H. diversicolor coelomocytes. The present study is the first to reveal the genotoxic and cytotoxic activity of 1 m T EMF in aquatic animals, and, consequently, the first one to report the adverse effect of this factor on common marine invertebrates and early life stages of fish.

Why Are DNA and Protein Electron Transfer So Different?

The corpus of electron transfer (ET) theory provides considerable power to describe the kinetics and dynamics of electron flow at the nanoscale. How is it, then, that nucleic acid (NA) ET continues to surprise, while protein-mediated ET is relatively free of mechanistic bombshells? I suggest that this difference originates in the distinct electronic energy landscapes for the two classes of reactions. In proteins, the donor/acceptor-to-bridge energy gap is typically several-fold larger than in NAs. NA ET can access tunneling, hopping, and resonant transport among the bases, and fluctuations can enable switching among mechanisms; protein ET is restricted to tunneling among redox active cofactors and, under strongly oxidizing conditions, a few privileged amino acid side chains. This review aims to provide conceptual unity to DNA and protein ET reaction mechanisms. The establishment of a unified mechanistic framework enabled the successful design of NA experiments that switch electronic coherence effects on and off for ET processes on a length scale of multiple nanometers and promises to provide inroads to directing and detecting charge flow in soft-wet matter.

Stability of RNA duplexes containing inosine·cytosine pairs

Inosine is found naturally in the anticodon loop of tRNA, is a product of adenosine deaminases that act on RNA, and can be used in oligonucleotide probes or to investigate the role of the exocyclic amino group of guanosine. Although the thermodynamics of I·U pairs in RNA have been systematically studied [Wright, D. J., Rice, J. L., Yanker, D. M., and Znosko, B. M. (2007) Biochemistry 46, 4625-4634], the thermodynamics of I·C pairs in RNA have not. Here, we have performed optical melting experiments on a series of RNA duplexes containing I·C pairs and compared their thermodynamics to the same duplexes containing A·C and G-C pairs. Nearest neighbor parameters for single I·C pairs adjacent to Watson-Crick pairs were derived. The derived nearest neighbor parameters are compared to those previously predicted blindly through a reweighting of energy-function collection with conformational ensemble sampling in Rosetta [Chou, F.-C., Kladwang, W., Kappel, K., and Das, R. (2016) Proc. Natl. Acad. Sci. U.S.A. 113, 8430-8435]. Scientists can use these nearest neighbor parameters to calculate the stability of ADAR products and to calculate the stability of an RNA duplex in which G-to-I substitution was used to determine the role of the exocyclic amino group of G.

Interaction of Tamoxifen Analogs With the Pocket Site of Some Hormone Receptors. A Molecular Docking and Density Functional Theory Study

In this paper, the antiestrogenic properties of Tamoxifen analogs have been investigated and a theoretical report of its analogs interaction with the pocket site of some hormone receptors are presented. Analogs were generated by modification of the hydrophilic functional group of Tamoxifen by hydroxyl, amide, carboxyl, and sulfhydryl functional groups, in an attempt to improve their activity and selectivity. The analogs exhibit a negative binding energy in the estrogen and progesterone receptors, which indicates a spontaneous interaction between the analogs and the pocket site in the hormone receptors. The values of the molecular polar surface area indicate that the analogs have good permeability and are strong electrophiles. The couplings showed electrostatic interactions such as hydrogen bond and π-π interactions. According with the Lipinsky Rule of Five, the four analogs presented a good biodistribution, permeability, and pharmacological action on the hormone receptors. The analysis of the charge transfer suggests a limited enhanced oxidative damage in the estrogen receptor that not takes place with the progesterone receptor.

Fluorescent excimers and exciplexes of the purine base derivative 8-phenylethynyl-guanine in DNA hairpins

The ground- and excited-state electronic interactions between the nucleobase analog 8-(4'-phenylethynyl)deoxyguanosine, EG, with natural nucleobases and 7-deazaguanine, as well as between adjacent EG base analogs, have been characterized using a combination of steady-state spectroscopy and time-resolved fluorescence, absorption, and stimulated Raman spectroscopies. The properties of the nucleoside EG-H2 are only weakly perturbed upon incorporation into synthetic DNA hairpins in which thymine, cytosine or adenine are the bases flanking EG. Incorporation of the nucleoside to be adjacent to guanine or deazaguanine results in the formation of short-lived (40-80 ps) exciplexes, the charge transfer character of which increases as the oxidation potential of the donor decreases. Hairpins possessing two or three adjacent EG base analogs display exciton-coupled circular dichroism in the ground state and form long-lived fluorescent excited states upon electronic excitation. Incorporation of EG into the helical scaffold of the DNA hairpins places it adjacent to its neighboring nucleobases or a second EG, thus providing the close proximity required for the formation of exciplex or excimer intermediates upon geometric relaxation of the short-lived EG excited state. The three time-resolved spectroscopic methods employed permit both the characterization of the several intermediates and the kinetics of their formation and decay.

Physico-chemical profiles of the wobble ↔ Watson-Crick G*·2AP(w) ↔ G·2AP(WC) and A·2AP(w) ↔ A*·2AP(WC) tautomerisations: a QM/QTAIM comprehensive survey

This study is intended to clarify in detail the tautomeric transformations of the wobble (w) G*·2AP(w) and A·2AP(w) nucleobase mispairs involving 2-aminopurine (2AP) into the Watson-Crick (WC) G·2AP(WC) and A*·2AP(WC) base mispairs (asterisks denote mutagenic tautomers of the DNA bases), respectively, by quantum-mechanical methods and Bader's Quantum Theory of Atoms in Molecules. Our previously reported methodology has been used, which allows the evolution of the physico-chemical parameters to be tracked along the entire internal reaction coordinate (IRC), not exclusively in the stationary states of these reactions. These biologically important G*·2AP(w) ↔ G·2AP(WC) and A·2AP(w) ↔ A*·2AP(WC) w ↔ WC tautomerisations, which are involved in mutagenic tautomerically-conformational pathways, determine the origin of the transitions and transversions induced by 2AP. In addition, it is established that they proceed through planar, highly stable, zwitterionic transition states and they exhibit similar physico-chemical profiles and stages of sequential intrapair proton transfer, followed by spatial rearrangement of the nucleobases relative to each other within the base pairs. These w ↔ WC tautomerisations occur non-dissociatively and are accompanied by a significant alteration in geometry (from wobble to Watson-Crick and vice versa) and redistribution of the specific intermolecular interactions, which can be divided into 10 patterns including AHB H-bonds and loosened A-H-B covalent bridges along the IRC of tautomerisation. Based on the redistribution of the geometrical and electron-topological parameters of the intrapair hydrogen bonds, exactly 9 key points have been allocated to characterize the evolution of these reactions.

HnRNP A1 Alters the Structure of a Conserved Enterovirus IRES Domain to Stimulate Viral Translation

Enteroviruses use a type I Internal Ribosome Entry Site (IRES) structure to facilitate protein synthesis and promote genome replication. Type I IRES elements require auxiliary host proteins to organize RNA structure for 40S ribosomal subunit assembly. Heterogeneous nuclear ribonucleoprotein A1 stimulates enterovirus 71 (EV71) translation in part through specific interactions with its stem loop II (SLII) IRES domain. Here, we determined a conjoined NMR-small angle x-ray scattering structure of the EV71 SLII domain and a mutant that significantly attenuates viral replication by abrogating hnRNP A1 interactions. Native SLII adopts a locally compact structure wherein stacking interactions in a conserved 5'-AUAGC-3' bulge preorganize the adjacent helices at nearly orthogonal orientations. Mutating the bulge sequence to 5'-ACCCC-3' ablates base stacking in the loop and globally reorients the SLII structure. Biophysical titrations reveal that the 5'-AUAGC-3' bulge undergoes a conformational change to assemble a functional hnRNP A1-RNA complex. Importantly, IRES mutations that delete the bulge impair viral translation and completely inhibit replication. Thus, this work provides key details into how an EV71 IRES structure adapts to hijack a cellular protein, and it suggests that the SLII domain is a potential target for antiviral therapy.

Alternative DNA structure formation in the mutagenic human c-MYC promoter

Mutation 'hotspot' regions in the genome are susceptible to genetic instability, implicating them in diseases. These hotspots are not random and often co-localize with DNA sequences potentially capable of adopting alternative DNA structures (non-B DNA, e.g. H-DNA and G4-DNA), which have been identified as endogenous sources of genomic instability. There are regions that contain overlapping sequences that may form more than one non-B DNA structure. The extent to which one structure impacts the formation/stability of another, within the sequence, is not fully understood. To address this issue, we investigated the folding preferences of oligonucleotides from a chromosomal breakpoint hotspot in the human c-MYC oncogene containing both potential G4-forming and H-DNA-forming elements. We characterized the structures formed in the presence of G4-DNA-stabilizing K+ ions or H-DNA-stabilizing Mg2+ ions using multiple techniques. We found that under conditions favorable for H-DNA formation, a stable intramolecular triplex DNA structure predominated; whereas, under K+-rich, G4-DNA-forming conditions, a plurality of unfolded and folded species were present. Thus, within a limited region containing sequences with the potential to adopt multiple structures, only one structure predominates under a given condition. The predominance of H-DNA implicates this structure in the instability associated with the human c-MYC oncogene.

On the origin of multiexponential fluorescence decays from 2-aminopurine-labeled dinucleotides

The fluorescent probe 2-aminopurine (2Ap) has been used for decades to study local conformational fluctuations in DNA. Steady-state and time-resolved measurements of 2Ap fluorescence have been used to predict specific conformational states through suitable modeling of the quenching of the fluorescence of a 2Ap residue incorporated site-specifically into a DNA strand. The success of this approach has been limited by a lack of understanding of the precise factors responsible for the complex, multiexponential decays observed experimentally. In this study, dinucleotides composed of 2Ap and adenine were studied by the time-correlated single-photon counting technique to investigate the causes of heterogeneous emission kinetics. Contrary to previous reports, we argue that emission from 2Ap that is stacked with a neighboring base contributes negligibly to the emission signals recorded more than 50 ps after excitation, which are instead dominated by emission from unstacked 2Ap. We find that the decay kinetics can be modeled using a continuous lifetime distribution, which arises from the inherent distance dependence of electron transfer rates without the need to postulate a small number of discrete states with decay times derived from multiexponential fits. These results offer a new perspective on the quenching of 2Ap fluorescence and expand the information that can be obtained from experiments.

Polarization-controlled optimal scatter suppression in transient absorption spectroscopy

Ultrafast transient absorption spectroscopy is a powerful technique to study fast photo-induced processes, such as electron, proton and energy transfer, isomerization and molecular dynamics, in a diverse range of samples, including solid state materials and proteins. Many such experiments suffer from signal distortion by scattered excitation light, in particular close to the excitation (pump) frequency. Scattered light can be effectively suppressed by a polarizer oriented perpendicular to the excitation polarization and positioned behind the sample in the optical path of the probe beam. However, this introduces anisotropic polarization contributions into the recorded signal. We present an approach based on setting specific polarizations of the pump and probe pulses, combined with a polarizer behind the sample. Together, this controls the signal-to-scatter ratio (SSR), while maintaining isotropic signal. We present SSR for the full range of polarizations and analytically derive the optimal configuration at angles of 40.5° between probe and pump and of 66.9° between polarizer and pump polarizations. This improves SSR by (or compared to polarizer parallel to probe). The calculations are validated by transient absorption experiments on the common fluorescent dye Rhodamine B. This approach provides a simple method to considerably improve the SSR in transient absorption spectroscopy.

Acridone in a biological nanocavity: detailed spectroscopic and docking analyses of probing both the tryptophan residues of bovine serum albumin

AD initially gets hooked to Trp 212 housed in domain IIA, inducing conformational changes in the protein and paving the way for the ligand to reach Trp 134 located in domain IB.

The Importance of Short- and Long-Range Exchange on Various Excited State Properties of DNA Monomers, Stacked Complexes, and Watson-Crick Pairs

We present a detailed analysis of several time-dependent DFT (TD-DFT) methods, including conventional hybrid functionals and two types of nonempirically tuned range-separated functionals, for predicting a diverse set of electronic excitations in DNA nucleobase monomers and dimers. This large and extensive set of excitations comprises a total of 50 different transitions (for each tested DFT functional) that includes several n → π and π → π* valence excitations, long-range charge-transfer excitations, and extended Rydberg transitions (complete with benchmark calculations from high-level EOM-CCSD(T) methods). The presence of localized valence excitations as well as extreme long-range charge-transfer excitations in these systems poses a serious challenge for TD-DFT methods that allows us to assess the importance of both short- and long-range exchange contributions for simultaneously predicting all of these various transitions. In particular, we find that functionals that do not have both short- and full long-range exchange components are unable to predict the different types of nucleobase excitations with the same accuracy. Most importantly, the current study highlights the importance of both short-range exchange and a nonempirically tuned contribution of long-range exchange for accurately predicting the diverse excitations in these challenging nucleobase systems.

Monitoring guanine photo-oxidation by enantiomerically resolved Ru(II) dipyridophenazine complexes using inosine-substituted oligonucleotides

The intercalating [Ru(TAP)2(dppz)](2+) complex can photo-oxidise guanine in DNA, although in mixed-sequence DNA it can be difficult to understand the precise mechanism due to uncertainties in where and how the complex is bound. Replacement of guanine with the less oxidisable inosine (I) base can be used to understand the mechanism of electron transfer (ET). Here the ET has been compared for both Λ- and Δ-enantiomers of [Ru(TAP)2(dppz)](2+) in a set of sequences where guanines in the readily oxidisable GG step in {TCGGCGCCGA}2 have been replaced with I. The ET has been monitored using picosecond and nanosecond transient absorption and picosecond time-resolved IR spectroscopy. In both cases inosine replacement leads to a diminished yield, but the trends are strikingly different for Λ- and Δ-complexes.

Nucleotides sequestered at different subsite loci within DNA-binding pockets of two OB-fold single-stranded DNA-binding proteins are unstacked to different extents

The gene 5 protein (g5p) encoded by the Ff strains of Escherichia coli bacteriophages is a dimeric single-stranded DNA-binding protein (SSB) that consists of two identical OB-fold (oligonucleotide/oligosaccharide-binding) motifs. Ultrafast time-resolved fluorescence measurements were carried out to investigate the effect of g5p binding on the conformation of 2-aminopurine (2AP) labels positioned between adenines or cytosines in the 16-nucleotide antiparallel tails of DNA hairpins. The measurements revealed significant changes in the conformational heterogeneity of the 2AP labels caused by g5p binding. The extent of the changes was dependent on sub-binding-site location, but generally resulted in base unstacking. When bound by g5p, the unstacked 2AP population increased from ∼ 22% to 59-67% in C-2AP-C segments and from 39% to 77% in an A-2AP-A segment. The OB-fold RPA70A domain of the human replication protein A also caused a significant amount of base unstacking at various locations within the DNA binding site as evidenced by steady-state fluorescence titration measurements using 2AP-labeled 5-mer DNAs. These solution studies support the concept that base unstacking at most of a protein's multiple sub-binding-site loci may be a feature that allows non-sequence specific OB-fold proteins to bind to single-stranded DNAs (ssDNAs) with minimal preference for particular sequences.

Power Frequency Magnetic Fields Affect the p38 MAPK-Mediated Regulation of NB69 Cell Proliferation Implication of Free Radicals

The proliferative response of the neuroblastoma line NB69 to a 100 µT, 50 Hz magnetic field (MF) has been shown mediated by activation of the MAPK-ERK1/2 pathway. This work investigates the MF effect on the cell cycle of NB69, the participation of p38 and c-Jun N-terminal (JNK) kinases in the field-induced proliferative response and the potential involvement of reactive oxygen species (ROS) in the activation of the MAPK-ERK1/2 and -p38 signaling pathways. NB69 cultures were exposed to the 100 µT MF, either intermittently for 24, 42 or 63 h, or continuously for periods of 15 to 120 min, in the presence or absence of p38 or JNK inhibitors: SB203580 and SP600125, respectively. Antioxidant N-acetylcysteine (NAC) was used as ROS scavenger. Field exposure induced transient activation of p38, JNK and ERK1/2. The MF proliferative effect, which was mediated by changes in the cell cycle, was blocked by the p38 inhibitor, but not by the JNK inhibitor. NAC blocked the field effects on cell proliferation and p38 activation, but not those on ERK1/2 activation. The MF-induced proliferative effects are exerted through sequential upregulation of MAPK-p38 and -ERK1/2 activation, and they are likely mediated by a ROS-dependent activation of p38.

Quantitative sampling of conformational heterogeneity of a DNA hairpin using molecular dynamics simulations and ultrafast fluorescence spectroscopy

Molecular dynamics (MD) simulations and time resolved fluorescence (TRF) spectroscopy were combined to quantitatively describe the conformational landscape of the DNA primary binding sequence (PBS) of the HIV-1 genome, a short hairpin targeted by retroviral nucleocapsid proteins implicated in the viral reverse transcription. Three 2-aminopurine (2AP) labeled PBS constructs were studied. For each variant, the complete distribution of fluorescence lifetimes covering 5 orders of magnitude in timescale was measured and the populations of conformers experimentally observed to undergo static quenching were quantified. A binary quantification permitted the comparison of populations from experimental lifetime amplitudes to populations of aromatically stacked 2AP conformers obtained from simulation. Both populations agreed well, supporting the general assumption that quenching of 2AP fluorescence results from pi-stacking interactions with neighboring nucleobases and demonstrating the success of the proposed methodology for the combined analysis of TRF and MD data. Cluster analysis of the latter further identified predominant conformations that were consistent with the fluorescence decay times and amplitudes, providing a structure-based rationalization for the wide range of fluorescence lifetimes. Finally, the simulations provided evidence of local structural perturbations induced by 2AP. The approach presented is a general tool to investigate fine structural heterogeneity in nucleic acid and nucleoprotein assemblies.

Whether the amino–imino tautomerism of 2-aminopurine is involved into its mutagenicity? Results of a thorough QM investigation

2AP* mutagenic tautomer is able to induce only one incorporation error – transversion – by pairing through the H-bonds into the G·2AP* mispair.

Site-Selective Monitoring of the Interaction of the SRA Domain of UHRF1 with Target DNA Sequences Labeled with 2-Aminopurine

UHRF1 plays a central role in the maintenance and transmission of epigenetic modifications by recruiting DNMT1 to hemimethylated CpG sites via its SET and RING-associated (SRA) domain, ensuring error-free duplication of methylation profiles. To characterize SRA-induced changes in the conformation and dynamics of a target 12 bp DNA duplex as a function of the methylation status, we labeled duplexes by the environment-sensitive probe 2-aminopurine (2-Ap) at various positions near or far from the central CpG recognition site containing either a nonmodified cytosine (NM duplex), a methylated cytosine (HM duplex), or methylated cytosines on both strands (BM duplex). Steady-state and time-resolved fluorescence indicated that binding of SRA induced modest conformational and dynamical changes in NM, HM, and BM duplexes, with only slight destabilization of base pairs, restriction of global duplex flexibility, and diminution of local nucleobase mobility. Moreover, significant restriction of the local motion of residues flanking the methylcytosine in the HM duplex suggested that these residues are more rigidly bound to SRA, in line with a slightly higher affinity of the HM duplex as compared to that of the NM or BM duplex. Our results are consistent with a "reader" role, in which the SRA domain scans DNA sequences for hemimethylated CpG sites without perturbation of the structure of contacted nucleotides.

Mode of DNA binding with γ-butyrolactone receptor protein CprB from Streptomyces coelicolor revealed by site-specific fluorescence dynamics

BACKGROUND

The γ-butyrolactone (GBL) binding transcription factors in Streptomyces species are known for their involvement in quorum sensing where they control the expression of various genes initiating secondary metabolic pathways. The structurally characterized member of this family CprB from Streptomyces coelicolor had earlier been demonstrated to bind a multitude of sequences containing a specific binding signature. Though structural breakthrough has been obtained for its complex with a consensus DNA sequence there is, however a dearth of information regarding the overall and site specific dynamics of protein-DNA interaction.

METHODS

To delineate the effect of CprB on the bound DNA, changes in motional dynamics of the fluorescent probe 2-aminopurine were monitored at three conserved base positions (5th, 12th and 23rd) for two DNA sequences: the consensus and the biologically relevant cognate element, on complex formation.

RESULTS

The changes in lifetime and generalized order parameter revealed a similarity in the binding pattern of the protein to both sequences with greater dynamic restriction at the end positions, 5th and 23rd, as compared to the middle 12th position. Also differences within this pattern demonstrated the influence of even small changes in sequence on protein interactions.

CONCLUSIONS

Here the study of motional dynamics was instrumental in establishing a structural footprint for the cognate DNA sequence and explaining the dynamics for the consensus DNA from structural correspondence.

GENERAL SIGNIFICANCE

Motional dynamics can be a powerful tool to efficiently study the mode of DNA binding to proteins that interact differentially with a plethora of DNA sequences, even in the absence of structural breakthrough.