Search

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.

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.

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.

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.

Interest in building dedicated Quantum Information Science and Engineering (QISE) education programs has greatly expanded in recent years. These programs are inherently convergent, complex, often resource intensive and likely require collaboration with a broad variety of stakeholders. In order to address this combination of challenges, we have captured ideas from many members in the community. This manuscript not only addresses policy makers and funding agencies (both public and private and from the regional to the international level) but also contains needs identified by industry leaders and discusses the difficulties inherent in creating an inclusive QISE curriculum. We report on the status of eighteen post-secondary education programs in QISE and provide guidance for building new programs. Lastly, we encourage the development of a comprehensive strategic plan for quantum education and workforce development as a means to make the most of the ongoing substantial investments being made in QISE., 18 pages, 2 figures, 1 table

Quantum annealers have grown in complexity to the point that quantum computations involving few thousands of qubits are now possible. In this paper, \textcolor{black}{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 \textcolor{black}{the ground state as well 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.

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.

The RNA-guided CRISPR–Cas9 nuclease has revolutionized genome engineering, yet its mechanism for DNA target selection is not fully understood. A crucial step in Cas9 target recognition involves unwinding of the DNA duplex to form a three-stranded R-loop structure. Work reported here demonstrates direct detection of Cas9-mediated DNA unwinding by a combination of site-directed spin labeling and molecular dynamics simulations. The results support a model in which the unwound nontarget strand is stabilized by a positively charged patch located between the two nuclease domains of Cas9 and reveal uneven increases in flexibility along the unwound nontarget strand upon scissions of the DNA backbone. This work establishes the synergistic combination of spin-labeling and molecular dynamics to directly monitor Cas9-mediated DNA conformational changes and yields information on the target DNA in different stages of Cas9 function, thus advancing mechanistic understanding of CRISPR–Cas9 and aiding future technological development.

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.

The recent availability of quantum annealers has fueled a new area of information technology where such devices are applied to address practically motivated and computationally difficult problems with hardware that exploits quantum mechanical phenomena. D-Wave annealers are promising platforms to solve these problems in the form of quadratic unconstrained binary optimization. Here we provide a formulation of the Chinese postman problem that can be used as a tool for probing the local connectivity of graphs and networks. We treat the problem classically with a tabu algorithm and simulated annealing, and using a D-Wave device. The efficiency of quantum annealing with respect to the simulated annealing has been demonstrated using the optimal time to solution metric. We systematically analyze computational parameters associated with the specific hardware. Our results clarify how the interplay between the embedding due to limited connectivity of the Chimera graph, the definition of logical qubits, and the role of spin-reversal controls the probability of reaching the expected solution.

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.

Proteins in the proximity of inorganic surfaces and nanoparticles may undergo profound adjustments that trigger biomedically relevant processes, such as protein fibrillation. The mechanisms that govern protein-surface interactions at the molecular level are still poorly understood. In this work, we investigate the adsorption onto a gold surface, in water, of an amyloid-β (Aβ) peptide, which is the amyloidogenic peptide involved in Alzheimer's disease. The entire adsorption process, from the peptide in bulk water to its conformational relaxation on the surface, is explored by large-scale atomistic molecular dynamics (MD) simulations. We start by providing a description of the conformational ensemble of Aβ in solution by a 22 μs temperature replica exchange MD simulation, which is consistent with previous results. Then, we obtain a statistical description of how the peptide approaches the gold surface by multiple MD simulations, identifying the preferential gold-binding sites and giving a kinetic picture of the association process. Finally, relaxation of the Aβ conformations at the gold/water interface is performed by a 19 μs Hamiltonian-temperature replica exchange MD simulation. We find that the conformational ensemble of Aβ is strongly perturbed by the presence of the surface. In particular, at the gold/water interface the population of the conformers akin to amyloid fibrils is significantly enriched, suggesting that this extended contact geometry may promote fibrillation.

Genome editing allows scientists to change an organism's DNA. One promising genome editing protocol, already validated in living organisms, is based on clustered regularly interspaced short palindromic repeats (CRISPR)/Cas protein-nucleic acid complexes. When the CRISPR/Cas approach was first demonstrated in 2012, its advantages with respect to previously available techniques, such as zinc-finger nucleases (ZFNs) and transcription activator-like effector nucleases (TALENs), immediately got attention and the method has seen a surge of experimental and computational investigations since then. However, the molecular mechanisms involved in target DNA recognition and cleavage are still not completely resolved and need further attention. The large size and complex nature of CRISPR/Cas9 complexes has been a challenge for computational studies, but some seed results exist and are illuminating on the cleavage activity. In this short review, we present recent progress in studying CRISPR/Cas9 systems by molecular dynamics simulations with coarse-grained and atomistic descriptions, including enhanced sampling., Comment: 24 pages, 10 figures, few figures are from previously published work for which appropriate permissions have been taken

Recent experiments suggest that DNA-mediated charge transport might enable signaling between the [4Fe4S] clusters in the C-terminal domains of human DNA primase and polymerase α, as well as the signaling between other replication and repair high-potential [4Fe4S] proteins. Our theoretical study demonstrates that the redox signaling cannot be accomplished exclusively by DNA-mediated charge transport because part of the charge transfer chain has an unfavorable free energy profile. We show that hole or excess electron transfer between a [4Fe4S] cluster and a nucleic acid duplex through a protein medium can occur within microseconds in one direction, while it is kinetically hindered in the opposite direction. We present a set of signaling mechanisms that may occur with the assistance of oxidants or reductants, using the allowed charge transfer processes. These mechanisms would enable the coordinated action of [4Fe4S] proteins on DNA, engaging the [4Fe4S] oxidation state dependence of the protein-DNA binding affinity.

The impact of nanotechnologies in biomedicine and biotechnology is becoming more and more evident. It imposes practical challenges, for instance, raising specific issues on the biocompatibility of nanostructures. Nanoparticles are characterized by a high surface-to-volume ratio, which makes them reactive to foreign species. Thus, when proteins or peptides approach an inorganic nanoparticle, as well as a flat surface, they are likely to interact with the substrate to some extent. This interaction is crucial for applications in drug delivery, imaging, diagnostics, implants, and other medical devices. Specifically, gold nanoparticles are highly versatile and particularly appealing. It is widely accepted that the surfaces of nanoparticles adsorb proteins either transiently in the soft corona layer or permanently in the hard corona layer. As a consequence, the protein structure and/or function may undergo profound adjustments or remain conserved. Detailing the interaction of different inorganic substrates with proteins and peptides at the atomic level, and designing ways to control the interaction, is the key for biomedical applications of nanoparticles, both from a fundamental viewpoint and for practical implementations. In the last decade, we have addressed protein–nanoparticle interactions, focusing on interfaces of gold surfaces and nanoparticles with amyloidogenic peptides and protein models. We have developed classical force fields, performed advanced molecular dynamics simulations, and compared computational outcomes with data from nuclear magnetic resonance experiments. Protein–gold complexes with differently coated gold nanoparticles have been modeled to explore the effects of charge and size on the protein structure. Our work unravels that a complex interplay between surface properties and characteristics of the biological adsorbate determines whether peptide conformation is influenced and whether protein aggregation is accelerated or inhibited by the presence of the substrate. General guidelines to cope with amyloidogenic proteins could be inferred: these can be essentially summarized with the necessity of balancing the hydrophobic and electrostatic interactions that the amyloidogenic proteins establish with the coating moieties.

Organic dyes are a viable alternative to silicon for energy conversion. Using simulations from first-principles, we show that chemical manipulation is a powerful tool for tuning the optical absorption spectra of a special class of dyes in a way that is convenient for exploitation in dye-sensitized solar cells. Specifically, we have carried out density functional theory calculations on three Ru(II) polypyridyl complexes with electron-donor ancillary ligands. These complexes were recently developed to study how different electron-donor ancillary ligands affect the photophysical and electrochemical properties of these dyes for light harvesting and photon-to-electron conversion efficiency. We found that the electron-donor ancillary ligands significantly enhance the light harvesting in the visible and the near-infrared regions relative to the reference dye N3. Furthermore, we detected a decrease in the ionization potential, which improves the energy alignment with the redox potentials of the electrolyte. These findings demonstrated that better organic materials for energy applications were developed. © 2014 American Chemical Society.

Herein, we report the isolation of new heterobimetallic complexes [Ni0.6Pd1.4ACHTUNGTRENUNG(EtCS2)4] (1), [NiPtACHTUNGTRENUNG(EtCS2)4] (2) and [Pd0.4Pt1.6ACHTUNGTRENUNG(EtCS2)4] (3), which were constructed by using transmetallation procedures. Subsequent oxidation with iodine furnished the MM'X monodimensional chains [Ni0.6Pt1.4ACHTUNGTRENUNG(EtCS2)4I] (4) and [Ni0.1Pd0.3Pt1.6ACHTUNGTRENUNG(EtCS2)4I] (5). The physical properties of these systems were investigated and the chain structures 4 and 5 were found to be reminiscent of the parent [Pt2ACHTUNGTRENUNG(EtCS2)4I] 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.

We present a combined investigation of the dynamics and optics of triplex DNA, based on classical molecular dynamics and time-dependent density functional theory. Our approach is devised to include the effects of conformational fluctuations on the electronic structure and optical excitations of the system. We find that the structural flexibility has a strong role in the determination of the optical signals. Our results allow us to unravel the peculiar fingerprints of Watson-Crick and Hoogsteen H-bonding in the optical absorption spectra. We find a specific optical absorption feature that is due to the simultaneous presence of the two H-bonding patterns in C(+)GC triplets. While this peculiar triplet signal is wiped out in some structures that are representative of the finite-temperature dynamics, it can be recovered in an average view, so that it is a pristine result of this work.

The interaction between proteins and the surfaces of inorganic materials is of paramount importance in natural systems. In recent years, the scope of polypeptide–surface interactions has been enlar...

The optical response of hybrid metal-semiconductor nanoparticles exhibits different behaviors due to the proximity between the disparate materials. For some hybrid systems, such as CdS-Au matchstick-shaped hybrids, the particles essentially retain the optical properties of their original components, with minor changes. Other systems, such as CdSe-Au dumbbell-shaped nanoparticles, exhibit significant change in the optical properties due to strong coupling between the two materials. Here, we study the absorption of these hybrids by comparing experimental results with simulations using the discrete dipole approximation method (DDA) employing dielectric functions of the bare components as inputs. For CdS-Au nanoparticles, the DDA simulation provides insights on the gold tip shape and its interface with the semiconductor, information that is difficult to acquire by experimental means alone. Furthermore, the qualitative agreement between DDA simulations and experimental data for CdS-Au implies that most effects influencing the absorption of this hybrid system are well described by local dielectric functions obtained separately for bare gold and CdS nanoparticles. For dumbbell shaped CdSe-Au, we find a shortcoming of the electrodynamic model, as it does not predict the "washing out" of the optical features of the semiconductor and the metal observed experimentally. The difference between experiment and theory is ascribed to strong interaction of the metal and semiconductor excitations, which spectrally overlap in the CdSe case. The present study exemplifies the employment of theoretical approaches used to describe the optical properties of semiconductors and metal nanoparticles, to achieve better understanding of the behavior of metal-semiconductor hybrid nanoparticles. © 2011 American Chemical Society.

G4-DNA molecules were recently reported as candidates for molecular electronics because of their higher stiffness and polarizability with respect to double-stranded DNA. Short G4-DNA structures have been traditionally studied in the past because of their possible presence in the human genome. The electronic structure of these molecules, unraveled in terms of the electron energy levels, constitutes crucial information in both technology and biology contexts. The discrete electron energy-level spectrum of isolated long single G4-DNA molecules, deposited on a gold substrate, is resolved by using cryogenic scanning tunneling spectroscopy in ultra-high-vacuum conditions. The patterns of the current−voltage curves manifest a considerable, although unexpected, variability, which required the development and implementation of a data-clustering methodology in order to resolve the reproducible curve patterns. On the basis of the sorted curves, we present a statistical analysis of the electronic levels and of the wi...

We propose a new method to extract the complete size-dependent complex dielectric function e(E) of semiconductor nanoparticles (NPs) from spectroscopic data. Typical experimental absorption spectra are available in a limited energy range near the absorption edge. We show that this limitation can be overcome by matching the NP dielectric function with the bulk dielectric function at high energies beyond the experimental range. Thus, using the Kramers−Kronig relations that link the real and imaginary parts, it is possible to extract the complete e(E) of the NP. This is achieved through an iterative procedure that systematically improves the trial e(E) until the absorption spectrum of the NP is accurately reproduced in the experimental range. Here we describe the methodology and we obtain the e(E) of CdSe and CdS NPs of selected sizes from published data sets.

By means of plane-wave pseudopotential periodic-supercell density functional theory calculations with a gradient-corrected exchange-correlation functional, we investigated the formation and the electronic structure of thin Au overlayers on CdSe(0001) and CdSe(000 (1) over bar) surfaces. We explored several possible Au/CdSe interfaces, including nonstoichiometric cases in which the very interface layer is mixed, namely, contains atoms of both the metal and the semiconductor. The relative formation energies of the computed model structures indicate that the formation of a very thin Au layer on CdSe surfaces can be epitaxial in the very early deposition stages but only in rather Au-rich conditions. The analysis of the band structures, densities of states, and wave functions for the low-energy interfaces reveals that hybridization occurs between the metal and the semiconductor electron states. This hybridization is confined at the very interface and is not expected to have significant consequences on the plasmonic and excitonic excitations that are appealing for nanotechnology applications of metal-semiconductor nanoparticles.

The electronic properties of several metal-modified Watson-Crick guanine-cytosine base pairs are investigated by rneans of first-principle density functional theory calculations. Focus is placed on a new structure recently proposed as a plausible model for building an antiparallel duplex with Zn-guanine-cytosine pairs, but we also inspect several other conformations and the incorporation of Ag and Cu ions. We analyze the effects induced by the incorporation of one metal cation per base pair by comparing the structures and the electronic properties of the metalated pairs to those of the natural guanine-cytosine pair, particularly for what concerns the modifications of energy levels and charge density distributions of the frontier orbitals. Our results reveal the establishment of covalent bonding between the metal cation and the nucleobases, identified in the presence of hybrid metal-guanine and metal-cytosine orbitals. Attachment of the cation can occur either at the N1 or the N7 site of guanine and is compatible with altering or not altering the H-bond pattern of the natural pair. Cu(II) strongly contributes to the hybridization of the orbitals around the band gap, whereas Ag(l) and Zn(II) give hybrid states farther from the band gap. Most metalated pairs have smaller band gaps than the natural guanine-cytosine pair. The band gap shrinking along with the metal-base coupling suggests interesting consequences for electron transfer through DNA double helices.

Controlling the polarity of polar semiconductors on nonpolar substrates offers a wealth of device concepts in the form of heteropolar junctions. A key to realize such structures is an appropriate buffer-layer design that, in the past, has been developed by empiricism. GaN or ZnO on sapphire are prominent examples for that. Understanding the basic processes that mediate polarity, however, is still an unsolved problem. In this work, we study the structure of buffer layers for group-III nitrides on sapphire by transmission electron microscopy as an example. We show that it is the conversion of the sapphire surface into a rhombohedral aluminum-oxynitride layer that converts the initial N-polar surface to Al polarity. With the various AlxOyNz phases of the pseudobinary Al2O3-AlN system and their tolerance against intrinsic defects, typical for oxides, a smooth transition between the octahedrally coordinated Al in the sapphire and the tetrahedrally coordinated Al in AlN becomes feasible. Based on these results, we discuss the consequences for achieving either polarity and shed light on widely applied concepts in the field of group-III nitrides like nitridation and low-temperature buffer layers.

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.

The electronic structure of periodic quadruple helix guanine wires, which mimic G4-DNA molecules, was studied as a function of the stacking distance between consecutive planes, by means of first principles density functional theory calculations. We show that, whereas for the native DNA interplane separation of 3.4 angstrom the HOMO- and LUMO-derived bands are poorly dispersive, the bandwidths can be significantly increased when compressive strain is applied alone, the helical axis. Our findings indicate that efficient band conduction for both holes and electrons can be supported by such wires for stacking distances below 2.6 angstrom, which imply a huge axial deformation with respect to double and quadruple helices in solutions and in crystals.

The adsorption of cysteine on the (111) surface of gold has been studied by means of periodic supercell density-functional theory calculations. A number of different adsorption modes are examined, including adsorption through the thiol group in either thiolate or disulfide form, and adsorption through both the thiol and amino functional groups. We find that at intermediate coverage densities the latter mode of adsorption is favored, followed by thiolate adsorption at the bridge (slightly displace toward fcc) site. The N-Au and S-Au bond strengths in the amino-thiolate adsorption are estimated to be of the order of 6 and 47 kcal/mol, respectively. The electronic structure of the different systems is analyzed, with focus on the total and projected density of states, as well as on the detailed character of the electronic states at the interface. States near the Fermi energy are found to have a metal-molecule antibonding character, whereas metal-molecule bonding states mostly occur near the lower edge of the Au-d band.

The RNA-guided CRISPR–Cas9 nuclease has revolutionized genome engineering, yet its mechanism for DNA target selection is not fully understood. A crucial step in Cas9 target recognition involves unwinding of the DNA duplex to form a three-stranded R-loop structure. Work reported here demonstrates direct detection of Cas9-mediated DNA unwinding by a combination of site-directed spin labeling and molecular dynamics simulations. The results support a model in which the unwound nontarget strand is stabilized by a positively charged patch located between the two nuclease domains of Cas9 and reveal uneven increases in flexibility along the unwound nontarget strand upon scissions of the DNA backbone. This work establishes the synergistic combination of spin-labeling and molecular dynamics to directly monitor Cas9-mediated DNA conformational changes and yields information on the target DNA in different stages of Cas9 function, thus advancing mechanistic understanding of CRISPR–Cas9 and aiding future technological ...

We present the results of a first principles investigation of the equilibrium properties of c-plane $\ensuremath{\alpha}\ensuremath{-}{\mathrm{Al}}_{2}{\mathrm{O}}_{3}$ surfaces. The stable structure for the $1\ifmmode\times\else\texttimes\fi{}1$ clean surface is Al terminated with a stoichiometric composition, while other terminations are unstable independent of surface preparation conditions. We discuss the implications of our results in the frame of possible extended reconstructions. For 1 monolayer of H coverage, we find that the preferred structure has OH dimers both perpendicular and nearly parallel to the surface. H-terminated surfaces may form in suitable preparation conditions. We discuss our results in terms of water adsorption and atomic layer epitaxy of $\ensuremath{\alpha}\ensuremath{-}{\mathrm{Al}}_{2}{\mathrm{O}}_{3}(0001).$

The interaction of DNA molecules with hard substrates is of paramount importance both for the study of DNA itself and for the variety of possible technological applications. Interaction with inorganic surfaces strongly modifies the helical shape of DNA. Hence, an accurate understanding of DNA structure and function at interfaces is a fundamental question with enormous impact in science and society. This work sets the fundamentals for the simulation of entire DNA oligomers on gold surfaces in dry and wet conditions. Thanks to the new GolDNA-AMBER force field, which was derived from first principles and includes dispersion interactions and polarization effects, we simulated self-assembled guanine and adenine monolayers on Au(111) in vacuo and the adsorption of all nucleobases on the same substrate in aqueous conditions. The periodic monolayers obtained from classical simulations match very well those from first principle calculations and experiments, assessing the robustness of the force field and motivating the application to more complex systems for which quantum calculations are not affordable and experiments are elusive. The energetics of nucleobases on Au(111) in solution reveal fundamental physicochemical effects: we find that the adsorption paradigm shifts from purely enthalpic to dominantly entropic by changing the environment and aggregation phase. © 2014 American Chemical Society.

We use a stepwise pulling protocol in molecular dynamics simulations to identify how a G-quadruplex selects and conducts Na(+), K(+), and NH4(+) ions. By estimating the minimum free-energy changes of the ions along the central channel via Jarzynski's equality, we find that the G-quadruplex selectively binds the ionic species in the following order: K(+)Na(+)NH4(+). This order implies that K(+) optimally fits the channel. However, the features of the free-energy profiles indicate that the channel conducts Na(+) best. These findings are in fair agreement with experiments on G-quadruplexes and reveal a profoundly different behavior from the prototype potassium-ion channel KcsA, which selects and conducts the same ionic species. We further show that the channel can also conduct a single file of water molecules and deform to leak water molecules. We propose a range for the conductance of the G-quadruplex.

We present the results of plane-wave pseudopotential periodic density functional theory (DFT) calculations on the geometries, energetics, and electronic structure of small molecules on Au(111). The chosen molecules---benzene, ammonia and cytosine---are representative of different adsorption regimes and interaction strengths. The chosen substrate is a prototype noble-metal surface that is widely employed as a support for organic materials. We assess the relevance of van der Waals effects in the adsorption process and the accuracy of different first-principle density functionals that have been recently developed to embody such effects. We find that there is no unique functional that is optimal for any system. In particular, our results reveal that functionals designed to reduce the short-term repulsion between the adsorbate and the substrate usually overestimate the adsorption strength and may even predict the wrong adsorption orientation. We show that an accurate description of the substrate does not ensure an accurate evaluation of the adsorption energetics, while the electronic structure is less sensitive to the specific choice. We propose the best choice for DFT calculations of DNA bases on Au(111) and similar systems in which both short-range and long-range interactions exist.

EF-hand calcium sensors respond structurally to changes in intracellular Ca2+ concentration, triggering diverse cellular responses and resulting in broad interactomes. Despite impressive advances in decoding their structure-function relationships, the folding mechanism of neuronal calcium sensors is still elusive. We used single-molecule optical tweezers to study the folding mechanism of the human neuronal calcium sensor 1 (NCS1). Two intermediate structures induced by Ca2+ binding to the EF-hands were observed during refolding. The complete folding of the C domain is obligatory for the folding of the N domain, showing striking interdomain dependence. Molecular dynamics results reveal the atomistic details of the unfolding process and rationalize the different domain stabilities during mechanical unfolding. Through constant-force experiments and hidden Markov model analysis, the free energy landscape of the protein was reconstructed. Our results emphasize that NCS1 has evolved a remarkable complex interdomain cooperativity and a fundamentally different folding mechanism compared to structurally related proteins. EF-hand calcium sensors respond structurally to changes in intracellular Ca concentration, triggering diverse cellular responses and resulting in broad interactomes. Despite impressive advances in decoding their structure-function relationships, the folding mechanism of neuronal calcium sensors is still elusive. We used single-molecule optical tweezers to study the folding mechanism of the human neuronal calcium sensor 1 (NCS1). Two intermediate structures induced by Ca binding to the EF-hands were observed during refolding. The complete folding of the C domain is obligatory for the folding of the N domain, showing striking interdomain dependence. Molecular dynamics results reveal the atomistic details of the unfolding process and rationalize the different domain stabilities during mechanical unfolding. Through constant-force experiments and hidden Markov model analysis, the free energy landscape of the protein was reconstructed. Our results emphasize that NCS1 has evolved a remarkable complex interdomain cooperativity and a fundamentally different folding mechanism compared to structurally related proteins.

We address the issue of whether chemical alterations of nucleobases are an effective tool to modulate charge transfer through DNA molecules. Our investigation uses a multilevel computational approach based on classical molecular dynamics and quantum chemistry. We find yet another piece of evidence that structural fluctuations are a key factor to determine the electronic structure of double-stranded DNA. We argue that the electronic structure and charge transfer ability of flexible polymers is the result of a complex intertwining of various structural, dynamical and chemical factors. Chemical intuition may be used to design molecular wires, but this is not the sole component in the complex charge transfer mechanism through DNA.

The interaction of molecules with inorganic substrates is a crucial issue forapplications in molecular electronics. It influences important factors suchas the immobilization efficiency and the charge injection through the inter-face. Moreover, mechanical aspects connected to the unfolding of biologicalmolecules are important.We hereby present recent efforts in our group to tackle these problems,based on density functional theory calculations. In particular, we discussour results on the adsorption of cytosine on Au(111) and on the interactionof guanine, in its natural and size-expanded forms, with small Au clusters.We find that cytosine binds to the Au(111) surface with a mechanism thatinvolves charge sharing, intermediate between chemisorption and physisorp-tion. The investigation of small complexes between guanine and gold clustersreveals the formation of hydrogen bonds: these configurations with unusualbonds are relevant at the corners of nanoparticles, while they can probably beneglected when DNA binds on flat extended metal surfaces.

The adsorption process of small organic molecules that represent reactive groups in amino acids (H2O, H2S, NH3, and HCOOH) on the nonpolar ZnS(10 (1) over bar0) surface was investigated by van der Waals corrected density functional theory calculations. At the accomplished interfaces, the oxygen, sulfur and nitrogen atoms of the adsorbates point toward the zinc atoms of the substrate, realizing electronic hybridization of their p lone pairs with the s and d bands of Zn. This electronic hybridization that involves surface cations is accompanied by H-bond formation that involves surface anions: this concerted mechanism enhances the interface strength and stability. On the basis of our results, we distinguish two classes of adsorption modes: molecular adsorption pertains to H2O, NH3, and HCOOH independently of the coverage and to H2S at low coverage, while concurrent adsorption/dissociation pertains to H2S at saturation coverage as a compromise between steric repulsion and H-bond-like interactions. Our results shed light on the passivation and modification of ZnS substrates (quantum dots and flat surfaces) in the prospect of technological and biomedical applications.