The impact of G-quadruplex dynamics on inter-tetrad electronic couplings: a hybrid computational study.

The G-quadruplex is a fascinating nucleic acid motif with implications in biology, medicine, and nanotechnologies. G-quadruplexes can form in the telomeres at the edges of chromosomes and in other guanine-rich regions of the genome. They can also be engineered for exploitation as biological materials for nanodevices. Their higher stiffness and higher charge transfer rates make them better candidates in nanodevices than duplex DNA. For the development of molecular nanowires, it is important to optimize electron transport along the wire axis. One powerful basis to do so is by manipulating the structure, based on known effects that structural changes have on electron transport. Here, we investigate such effects, by a combination of classical simulations of the structure and dynamics and quantum calculations of electronic couplings. We find that this structure-function relationship is complex. A single helix shape parameter alone does not embody such complexity, but rather a combination of distances and angles between stacked bases influences charge transfer efficiency. By analyzing linear combinations of shape descriptors for different topologies, we identify the structural features that most affect charge transfer efficiency. We discuss the transferability of the proposed model and the limiting effects of inherent flexibility.

Landscape of DNA binding signatures of myocyte enhancer factor-2B reveals a unique interplay of base and shape readout.

Myocyte enhancer factor-2B (MEF2B) has the unique capability of binding to its DNA target sites with a degenerate motif, while still functioning as a gene-specific transcriptional regulator. Identifying its DNA targets is crucial given regulatory roles exerted by members of the MEF2 family and MEF2B's involvement in B-cell lymphoma. Analyzing structural data and SELEX-seq experimental results, we deduced the DNA sequence and shape determinants of MEF2B target sites on a high-throughput basis in vitro for wild-type and mutant proteins. Quantitative modeling of MEF2B binding affinities and computational simulations exposed the DNA readout mechanisms of MEF2B. The resulting binding signature of MEF2B revealed distinct intricacies of DNA recognition compared to other transcription factors. MEF2B uses base readout at its half-sites combined with shape readout at the center of its degenerate motif, where A-tract polarity dictates nuances of binding. The predominant role of shape readout at the center of the core motif, with most contacts formed in the minor groove, differs from previously observed protein-DNA readout modes. MEF2B, therefore, represents a unique protein for studies of the role of DNA shape in achieving binding specificity. MEF2B-DNA recognition mechanisms are likely representative for other members of the MEF2 family.

The RNA Polymerase α Subunit Recognizes the DNA Shape of the Upstream Promoter Element.

We demonstrate here that the α subunit C-terminal domain of Escherichia coli RNA polymerase (αCTD) recognizes the upstream promoter (UP) DNA element via its characteristic minor groove shape and electrostatic potential. In two compositionally distinct crystallized assemblies, a pair of αCTD subunits bind in tandem to the UP element consensus A-tract that is 6 bp in length (A6-tract), each with their arginine 265 guanidinium group inserted into the minor groove. The A6-tract minor groove is significantly narrowed in these crystal structures, as well as in computationally predicted structures of free and bound DNA duplexes derived by Monte Carlo and molecular dynamics simulations, respectively. The negative electrostatic potential of free A6-tract DNA is substantially enhanced compared to that of generic DNA. Shortening the A-tract by 1 bp is shown to "knock out" binding of the second αCTD through widening of the minor groove. Furthermore, in computationally derived structures with arginine 265 mutated to alanine in either αCTD, either with or without the "knockout" DNA mutation, contact with the DNA is perturbed, highlighting the importance of arginine 265 in achieving αCTD-DNA binding. These results demonstrate that the importance of the DNA shape in sequence-dependent recognition of DNA by RNA polymerase is comparable to that of certain transcription factors.

Vacancies in graphene: an application of adiabatic quantum optimization.

Quantum annealers have grown in complexity to the point that quantum computations involving a few thousand qubits are now possible. In this paper, with the intentions to show the feasibility of quantum annealing to tackle problems of physical relevance, we used a simple model, compatible with the capability of current quantum annealers, to study the relative stability of graphene vacancy defects. By mapping the crucial interactions that dominate carbon-vacancy interchange onto a quadratic unconstrained binary optimization problem, our approach exploits the ground state as well as the excited states found by the quantum annealer to extract all the possible arrangements of multiple defects on the graphene sheet together with their relative formation energies. This approach reproduces known results and provides a stepping stone towards applications of quantum annealing to problems of physical-chemical interest.

Adsorption Mechanisms of Nucleobases on the Hydrated Au(111) Surface.

The solution environment is of fundamental importance in the adsorption of molecules on surfaces, a process that is strongly affected by the capability of the adsorbate to disrupt the hydration layer above the surface. Here we disclose how the presence of interface water influences the adsorption mechanism of DNA nucleobases on a gold surface. By means of metadynamics simulations, we describe the distinctive features of a complex free-energy landscape for each base, which manifests activation barriers for the adsorption process. We characterize the different pathways that allow each nucleobase to overcome the barriers and be adsorbed on the surface, discussing how they influence the kinetics of adsorption of single-stranded DNA oligomers with homogeneous sequences. Our findings offer a rationale as to why experimental data on the adsorption of single-stranded homo-oligonucleotides do not straightforwardly follow the thermodynamics affinity rank.

Conformations of Human Telomeric G-Quadruplex Studied Using a Nucleotide-Independent Nitroxide Label.

Guanine-rich oligonucleotides can form a unique G-quadruplex (GQ) structure with stacking units of four guanine bases organized in a plane through Hoogsteen bonding. GQ structures have been detected in vivo and shown to exert their roles in maintaining genome integrity and regulating gene expression. Understanding GQ conformation is important for understanding its inherent biological role and for devising strategies to control and manipulate functions based on targeting GQ. Although a number of biophysical methods have been used to investigate structure and dynamics of GQs, our understanding is far from complete. As such, this work explores the use of the site-directed spin labeling technique, complemented by molecular dynamics simulations, for investigating GQ conformations. A nucleotide-independent nitroxide label (R5), which has been previously applied for probing conformations of noncoding RNA and DNA duplexes, is attached to multiple sites in a 22-nucleotide DNA strand derived from the human telomeric sequence (hTel-22) that is known to form GQ. The R5 labels are shown to minimally impact GQ folding, and inter-R5 distances measured using double electron-electron resonance spectroscopy are shown to adequately distinguish the different topological conformations of hTel-22 and report variations in their occupancies in response to changes of the environment variables such as salt, crowding agent, and small molecule ligand. The work demonstrates that the R5 label is able to probe GQ conformation and establishes the base for using R5 to study more complex sequences, such as those that may potentially form multimeric GQs in long telomeric repeats.

Control of DNA minor groove width and Fis protein binding by the purine 2-amino group.

The width of the DNA minor groove varies with sequence and can be a major determinant of DNA shape recognition by proteins. For example, the minor groove within the center of the Fis-DNA complex narrows to about half the mean minor groove width of canonical B-form DNA to fit onto the protein surface. G/C base pairs within this segment, which is not contacted by the Fis protein, reduce binding affinities up to 2000-fold over A/T-rich sequences. We show here through multiple X-ray structures and binding properties of Fis-DNA complexes containing base analogs that the 2-amino group on guanine is the primary molecular determinant controlling minor groove widths. Molecular dynamics simulations of free-DNA targets with canonical and modified bases further demonstrate that sequence-dependent narrowing of minor groove widths is modulated almost entirely by the presence of purine 2-amino groups. We also provide evidence that protein-mediated phosphate neutralization facilitates minor groove compression and is particularly important for binding to non-optimally shaped DNA duplexes.

DNAshape: a method for the high-throughput prediction of DNA structural features on a genomic scale.

We present a method and web server for predicting DNA structural features in a high-throughput (HT) manner for massive sequence data. This approach provides the framework for the integration of DNA sequence and shape analyses in genome-wide studies. The HT methodology uses a sliding-window approach to mine DNA structural information obtained from Monte Carlo simulations. It requires only nucleotide sequence as input and instantly predicts multiple structural features of DNA (minor groove width, roll, propeller twist and helix twist). The results of rigorous validations of the HT predictions based on DNA structures solved by X-ray crystallography and NMR spectroscopy, hydroxyl radical cleavage data, statistical analysis and cross-validation, and molecular dynamics simulations provide strong confidence in this approach. The DNAshape web server is freely available at http://rohslab.cmb.usc.edu/DNAshape/.

Unraveling the mechanism of tip-enhanced molecular energy transfer.

Electronic Energy Transfer (EET) between chromophores is fundamental in many natural light-harvesting complexes, serving as a critical step for solar energy funneling in photosynthetic plants and bacteria. The complicated role of the environment in mediating this process in natural architectures has been addressed by recent scanning tunneling microscope experiments involving EET between two molecules supported on a solid substrate. These measurements demonstrated that EET in such conditions has peculiar features, such as a steep dependence on the donor-acceptor distance, reminiscent of a short-range mechanism more than of a Förster-like process. By using state of the art hybrid ab initio/electromagnetic modeling, here we provide a comprehensive theoretical analysis of tip-enhanced EET. In particular, we show that this process can be understood as a complex interplay of electromagnetic-based molecular plasmonic processes, whose result may effectively mimic short range effects. Therefore, the established identification of an exponential decay with Dexter-like effects does not hold for tip-enhanced EET, and accurate electromagnetic modeling is needed to identify the EET mechanism.

A Perspective on Sustainable Computational Chemistry Software Development and Integration.

The power of quantum chemistry to predict the ground and excited state properties of complex chemical systems has driven the development of computational quantum chemistry software, integrating advances in theory, applied mathematics, and computer science. The emergence of new computational paradigms associated with exascale technologies also poses significant challenges that require a flexible forward strategy to take full advantage of existing and forthcoming computational resources. In this context, the sustainability and interoperability of computational chemistry software development are among the most pressing issues. In this perspective, we discuss software infrastructure needs and investments with an eye to fully utilize exascale resources and provide unique computational tools for next-generation science problems and scientific discoveries.

Electrical behaviour of heterobimetallic [MM'(EtCS2)4] (MM'=NiPd, NiPt, PdPt) and MM'X-chain polymers [PtM(EtCS2)4I] (M=Ni, Pd).

Herein, we report the isolation of new heterobimetallic complexes [Ni(0.6)Pd(1.4)(EtCS(2))(4)] (1), [NiPt(EtCS(2))(4)] (2) and [Pd(0.4)Pt(1.6)(EtCS(2))(4)] (3), which were constructed by using transmetallation procedures. Subsequent oxidation with iodine furnished the MM'X monodimensional chains [Ni(0.6)Pt(1.4)(EtCS(2))(4)I] (4) and [Ni(0.1)Pd(0.3)Pt(1.6)(EtCS(2))(4)I] (5). The physical properties of these systems were investigated and the chain structures 4 and 5 were found to be reminiscent of the parent [Pt(2)(EtCS(2))(4)I] species. However, they were more sensitively dependent on the localised nature of the charge on the Ni ion, which caused spontaneous breaking of the conduction bands.

Combined effects of metal complexation and size expansion in the electronic structure of DNA base pairs.

Novel DNA derivatives have been recently investigated in the pursuit of modified DNA duplexes to tune the electronic structure of DNA-based assemblies for nanotechnology applications. Size-expanded DNAs (e.g., xDNA) and metalated DNAs (M-DNA) may enhance stacking interactions and induce metallic conductivity, respectively. Here we explore possible ways of tailoring the DNA electronic structure by combining the aromatic size expansion with the metal-doping. We select the salient structures from our recent study on natural DNA pairs complexed with transition metal ions and consider the equivalent model configurations for xDNA pairs. We present the results of density functional theory electronic structure calculations of the metalated expanded base-pairs with various localized basis sets and exchange-correlation functionals. Implicit solvent and coordination water molecules are also included. Our results indicate that the effect of base expansion is largest in Ag-xGC complexes, while Cu-xGC complexes are the most promising candidates for nanowires with enhanced electron transfer and also for on-purpose modification of the DNA double-helix for signal detection.

LCAO Electronic Structure of Nucleic Acid Bases and Other Heterocycles and Transfer Integrals in B-DNA, Including Structural Variability.

To describe the molecular electronic structure of nucleic acid bases and other heterocycles, we employ the Linear Combination of Atomic Orbitals (LCAO) method, considering the molecular wave function as a linear combination of all valence orbitals, i.e., 2s, 2px, 2py, 2pz orbitals for C, N, and O atoms and 1s orbital for H atoms. Regarding the diagonal matrix elements (also known as on-site energies), we introduce a novel parameterization. For the non-diagonal matrix elements referring to neighboring atoms, we employ the Slater-Koster two-center interaction transfer integrals. We use Harrison-type expressions with factors slightly modified relative to the original. We compare our LCAO predictions for the ionization and excitation energies of heterocycles with those obtained from Ionization Potential Equation of Motion Coupled Cluster with Singles and Doubles (IP-EOMCCSD)/aug-cc-pVDZ level of theory and Completely Normalized Equation of Motion Coupled Cluster with Singles, Doubles, and non-iterative Triples (CR-EOMCCSD(T))/aug-cc-pVDZ level of theory, respectively, (vertical values), as well as with available experimental data. Similarly, we calculate the transfer integrals between subsequent base pairs, to be used for a Tight-Binding (TB) wire model description of charge transfer and transport along ideal or deformed B-DNA. Taking into account all valence orbitals, we are in the position to treat deflection from the planar geometry, e.g., DNA structural variability, a task impossible for the plane Hückel approach (i.e., using only 2pz orbitals). We show the effects of structural deformations utilizing a 20mer evolved by Molecular Dynamics.

Effects of Structural Dynamics on Charge Carrier Transfer in B-DNA: A Combined MD and RT-TDDFT Study.

Hole transfer along the axis of duplex DNA has been the focus of physical chemistry research for decades, with implications in diverse fields, from nanotechnology to cell oxidative damage. Computational approaches are particularly amenable for this problem, to complement experimental data for interpretation of transfer mechanisms. To be predictive, computational results need to account for the inherent mobility of biological molecules during the time frame of experimental measurements. Here, we address the structural variability of B-DNA and its effects on hole transfer in a combined molecular dynamics (MD) and real-time time-dependent density functional theory (RT-TDDFT) study. Our results show that quantities that characterize the charge transfer process, such as the time-dependent dipole moment and hole population at a specific site, are sensitive to structural changes that occur on the nanosecond time scale. We extend the range of physical properties for which such a correlation has been observed, further establishing the fact that quantitative computational data on charge transfer properties should include statistical averages. Furthermore, we use the RT-TDDFT results to assess an efficient tight-binding method suitable for high-throughput predictions. We demonstrate that charge transfer, although affected by structural variability, on average, remains strong in AA and GG dimers.

Hydroxyl-rich beta-sheet adhesion to the gold surface in water by first-principle simulations.

Proteins able to recognize inorganic surfaces are of paramount importance for living organisms. Mimicking nature, surface-recognizing proteins and peptides have also been man-made by combinatorial biochemistry. However, to date the recognition mechanisms remain elusive, and the underlying physicochemical principles are still unknown. Selectivity of gold-binding peptides (cysteine-free and rich in hydroxyl amino acids) is particularly puzzling, since the most relevant gold surface, Au(111), is known to be chemically inert and atomically flat. Using atomistic first-principle simulations we show that weak chemical interactions of dative-bond character confer to a prototype secondary structure (an antiparallel beta-sheet made of hydroxyl amino acids) and its hydration layer the capability of discriminating among gold surface sites. Our results highlight the unexpected role of hydration water in this process, suggesting that hydrophilic amino acids and their hydration shell cooperate to contribute to protein-gold surface recognition.

Protein-surface interactions: challenging experiments and computations.

Protein-surface interactions are fundamental in natural processes, and have great potential for applications ranging from nanotechnology to medicine. A recent workshop highlighted the current achievements and the main challenges in the field.

Protein-Mutation-Induced Conformational Changes of the DNA and Nuclease Domain in CRISPR/Cas9 Systems by Molecular Dynamics Simulations.

Class 2 CRISPR (clustered regularly interspaced short palindromic repeats) systems offer a unique protocol for genome editing in eukaryotic cells. The nuclease activity of Cas9 has been harnessed to perform precise genome editing by creating double-strand breaks. However, the nuclease activity of Cas9 can be triggered when there is imperfect complementarity between the RNA guide sequence and an off-target genomic site, which is a major limitation of the CRISPR technique for practical applications. Hence, understanding the binding mechanisms in CRISPR/Cas9 for predicting ways to increase cleavage specificity is a timely research target. One way to understand and tune the binding strength is to study wild-type and mutant Cas9, in complex with a guide RNA and a target DNA. We have performed classical all-atom MD simulations over a cumulative time scale of 13.5 µs of CRISPR/Cas9 ternary complexes with the wild-type Cas9 from Streptococcus pyogenes and three of its mutants: K855A, H982A, and the combination K855A+H982A, selected from the outcome of experimental work. Our results reveal significant structural impact of the mutations, with implications for specificity. We find that the "unwound" part of the nontarget DNA strand exhibits enhanced flexibility in complexes with Cas9 mutants and tries to move away from the HNH/RuvC interface, where it is otherwise stabilized by electrostatic couplings in the wild-type complex. Our findings refine an electrostatic model by which cleavage specificity can be optimized through protein mutations.

Porphyrin intercalation in G4-DNA quadruplexes by molecular dynamics simulations.

We investigated the intercalation of a porphyrin derivative (TMPyP) between guanine tetrads (G4-tetrads, G4t's) in 4-stranded G4-DNA oligomers by classical molecular dynamics simulations. Contrary to experimental evidence on very short oligomers that contain stabilizing cations, we find that TMPyP can stack with the G4-tetrads in the absence of interplane cations. A high TMPyP/G4t stoichiometric ratio of 1/2 induces strong deformations of the G4-quadruplexes. A lower ratio of 1/8 is better compatible with the helical conformation. When a TMPyP is accommodated between two tetrads, the stacking distance between the intercalated molecule and a G4-tetrad is approximately 4.3-4.7 A. We find that the possibility of regular TMPyP intercalation depends on the length of the quadruplex, on the stoichiometric ratio and on the edge termination motif.

First principles effective electronic couplings for hole transfer in natural and size-expanded DNA.

Hole transfer processes between base pairs in natural DNA and size-expanded DNA (xDNA) are studied and compared, by means of an accurate first principles evaluation of the effective electronic couplings (also known as transfer integrals), in order to assess the effect of the base augmentation on the efficiency of charge transport through double-stranded DNA. According to our results, the size expansion increases the average electronic coupling, and thus the CT rate, with potential implications in molecular biology and in the implementation of molecular nanoelectronics. Our analysis shows that the effect of the nucleobase expansion on the charge-transfer (CT) rate is sensitive to the sequence of base pairs. Furthermore, we find that conformational variability is an important factor for the modulation of the CT rate. From a theoretical point of view, this work offers a contribution to the CT chemistry in pi-stacked arrays. Indeed, we compare our methodology against other standard computational frameworks that have been adopted to tackle the problem of CT in DNA, and unravel basic principles that should be accounted for in selecting an appropriate theoretical level.

Fingerprints of bonding motifs in DNA duplexes of adenine and thymine revealed from circular dichroism: synchrotron radiation experiments and TDDFT calculations.

Synchrotron radiation circular dichroism (SRCD) spectra were recorded for a family of 12 DNA duplexes that all contain nine adenines (A) and nine thymines (T) in each strand but in different combinations. The total number of AT Watson-Crick (WC) base pairs is constant (18), but the number of cross-strand (CS) hydrogen bonds between A and T varies between 0 and 16, the maximum possible. Eleven of the duplexes have one or more A tracts, and one duplex has T tracts. The signals due to hybridization were found from subtraction of spectra of single strands from spectra of the duplexes. The residual spectrum of the T-tract duplex T(9)A(9):A(9)T(9) (5'-3':3'-5') significantly differs from that of the A-tract duplex A(9)T(9):T(9)A(9), but only below 210 nm, which suggests that the signal in this region depends on the superhelicity of the duplex. A principal component analysis of all residual spectra reveals that spectra of A-tract duplexes can be obtained to a good approximation as a linear combination of just two basis spectra. The first component is assigned to the spectrum of 18 WC and 8 CS pairs, whereas the second component is that of 8 CS pairs. This interpretation is supported by separate experiments on duplexes of varying lengths but with similar arrangements of the A and T's and by experiments on two other duplex families of 14 and 30 base pairs. The best correlation is obtained by the assumption that cross-strand interactions occur as long as there are two adenine neighbors in a strand. Our data indicate that a circular dichroism spectrum of a duplex containing only A and T can simply be inferred from the number of WC base pairs and the number of CS interactions, and we provide reference spectra for these two interactions. Finally, time dependent density functional theory calculations of the circular dichroism spectra for an isolated WC base pair and two different CS base pairs (between adenine N-6 amine and thymine O-4 or between adenine C-2-H and thymine O-2) were performed to provide some additional support for the interpretation of the experimental spectra. We find large differences between the two calculated CS spectra. However, there is a reasonable qualitative agreement between the calculated WC and the C-2-H...O-2 CS spectra and those deduced from the experimental data.