Competing interactions in DNA assembly on graphene

We study the patterns that short strands of single-stranded DNA form on the top graphene surface of graphite. We find that the DNA assembles into two distinct patterns, small spherical particles and elongated networks. Known interaction models based on DNA-graphene binding, hydrophobic interactions, or models based on the purine/pyrimidine nature of the bases do not explain our observed crossover in pattern formation. We argue that the observed assembly behavior is caused by a crossover in the competition between base-base pi stacking and base-graphene pi stacking and we infer a critical crossover energy of eV. The experiments therefore provide a projective measurement of the base-base interaction strength.

A van der Waals density functional study of adenine on graphene: single-molecular adsorption and overlayer binding

The adsorption of an adenine molecule on graphene is studied using a first-principles van der Waals functional, vdW-DF (Dion et al 2004 Phys. Rev. Lett. 92 246401). The cohesive energy of an ordered adenine overlayer is also estimated. For the adsorption of a single molecule, we determine the optimal binding configuration and adsorption energy by translating and rotating the molecule. The adsorption energy for a single molecule of adenine is found to be 711 meV, which is close to the calculated adsorption energy of the similarly sized naphthalene. On the basis of the single-molecular binding configuration, we estimate the cohesive energy of a two-dimensional ordered overlayer. We find a significantly stronger binding energy for the ordered overlayer than for single-molecule adsorption.

Adsorption and desorption of DNA on graphene oxide studied by fluorescently labeled oligonucleotides

Being the newest member of the carbon materials family, graphene possesses many unique physical properties resulting is a wide range of applications. Recently, it was discovered that graphene oxide can effectively adsorb DNA, and at the same time, it can completely quench adsorbed fluorophores. These properties make it possible to prepare DNA-based optical sensors using graphene oxide. While practical analytical applications are being demonstrated, the fundamental understanding of binding between graphene oxide and DNA in solution received relatively less attention. In this work, we report that the adsorption of 12-, 18-, 24-, and 36-mer single-stranded DNA on graphene oxide is affected by several factors. For example, shorter DNAs are adsorbed more rapidly and bind more tightly to the surface of graphene. The adsorption is favored by a lower pH and a higher ionic strength. The presence of organic solvents such as ethanol can either increase or decrease adsorption depending on the ionic strength of the solution. By adding the cDNA, close to 100% desorption of the absorbed DNA on graphene can be achieved. On the other hand, if temperature is increased, only a small percentage of DNA is desorbed. Further, the adsorbed DNA can also be exchanged by free DNA in solution. These findings are important for further understanding of the interactions between DNA and graphene and for the optimization of DNA and graphene-based devices and sensors.

Fast DNA sequencing with a graphene-based nanochannel device

Devices in which a single strand of DNA is threaded through a nanopore could be used to efficiently sequence DNA. However, various issues will have to be resolved to make this approach practical, including controlling the DNA translocation rate, suppressing stochastic nucleobase motions, and resolving the signal overlap between different nucleobases. Here, we demonstrate theoretically the feasibility of DNA sequencing using a fluidic nanochannel functionalized with a graphene nanoribbon. This approach involves deciphering the changes that occur in the conductance of the nanoribbon as a result of its interactions with the nucleobases via π-π stacking. We show that as a DNA strand passes through the nanochannel, the distinct conductance characteristics of the nanoribbon (calculated using a method based on density functional theory coupled to non-equilibrium Green function theory) allow the different nucleobases to be distinguished using a data-mining technique and a two-dimensional transient autocorrelation analysis. This fast and reliable DNA sequencing device should be experimentally feasible in the near future.

Parameterization of a B3LYP specific correction for non-covalent interactions and basis set superposition error on a gigantic dataset of CCSD(T) quality non-covalent interaction energies

A vast number of non-covalent interaction energies at the counterpoise corrected CCSD(T) level have been collected from the literature to build a diverse new dataset. The whole dataset, which consists of 2027 CCSD(T) energies, includes most of the published data at this level. A large subset of the data was then used to train a novel, B3LYP specific, empirical correction scheme for non-covalent interactions and basis set superposition error (abbreviated as B3LYP-MM). Results obtained with our new correction scheme were directly compared to benchmark results obtained with B3LYP-D3(1) and M06-2X(2) (two popular density functions designed specifically to accurately model non-covalent interactions). For non-covalent complexes dominated by dispersion or dipole-dipole interactions all three tested methods give accurate results with the medium size aug-cc-pVDZ(3-6) basis set with MUE's of 0.27 (B3LYP-MM), 0.32 (B3LYP-D3) and 0.47 kcal/mol (M06-2X) (with explicit counterpoise corrections). These results validate both B3LYP-D3 and M06-2X for interactions of this type using a much larger data set than was presented in prior work. However, our new dispersion correction scheme shows some clear advantages for dispersion and dipole-dipole dominated complexes with the small LACVP* basis set, which is very popular in use due to its low associated computational cost: The MUE for B3LYP-MM with the LACVP* basis set for this subset of complexes (without explicit counterpoise corrections) is only 0.28 kcal/mol, compared to 0.65 kcal/mol for M06-2X or 1.16 kcal/mol for B3LYP-D3. Additionally, our new correction scheme also shows major improvements in accuracy for hydrogen-bonded systems and for systems involving ionic interactions, for example cation-π interactions. Compared to B3LYP-D3 and M06-2X, we also find that our new B3LYP-MM correction scheme gives results of higher or equal accuracy for a large dataset of conformer energies of di- and tripeptides, sugars, and cysteine.

Efficient and Accurate Double-Hybrid-Meta-GGA Density Functionals-Evaluation with the Extended GMTKN30 Database for General Main Group Thermochemistry, Kinetics, and Noncovalent Interactions

We present an extended and improved version of our recently published database for general main group thermochemistry, kinetics, and noncovalent interactions [J. Chem. Theory Comput. 2010, 6, 107], which is dubbed GMTKN30. Furthermore, we suggest and investigate two new double-hybrid-meta-GGA density functionals called PTPSS-D3 and PWPB95-D3. PTPSS-D3 is based on reparameterized TPSS exchange and correlation contributions; PWPB95-D3 contains reparameterized PW exchange and B95 parts. Both functionals contain fixed amounts of 50% Fock-exchange. Furthermore, they include a spin-opposite scaled perturbative contribution and are combined with our latest atom-pairwise London-dispersion correction [J. Chem. Phys. 2010, 132, 154104]. When evaluated with the help of the Laplace transformation algorithm, both methods scale as N(4) with system size. The functionals are compared with the double hybrids B2PLYP-D3, B2GPPLYP-D3, DSD-BLYP-D3, and XYG3 for GMTKN30 with a quadruple-ζ basis set. PWPB95-D3 and DSD-BLYP-D3 are the best functionals in our study and turned out to be more robust than B2PLYP-D3 and XYG3. Furthermore, PWPB95-D3 is the least basis set dependent and the best functional at the triple-ζ level. For the example of transition metal carbonyls, it is shown that, mainly due to the lower amount of Fock-exchange, PWPB95-D3 and PTPSS-D3 are better applicable than the other double hybrids. Finally, we discuss in some detail the XYG3 functional [Proc. Nat. Acad. Sci. U.S.A. 2009, 106, 4963], which makes use of B3LYP orbitals and electron densities. We show that it is basically a highly nonlocal variant of B2PLYP and that its partially good performance is mainly due to a larger effective amount of perturbative correlation compared to other double hybrids. We finally recommend the PWPB95-D3 functional in general chemistry applications.

Label-free electrical detection of DNA hybridization using carbon nanotubes and graphene

The interface between biosystems and nanomaterials is emerging for detection of various biomolecules and subtle cellular activities. In particular, the development of cost-effective and sequence-selective DNA detection is urgent for the diagnosis of genetic or pathogenic diseases. Graphene-based nanocarbon materials, such as carbon nanotubes and thin graphene layers, have been employed as biosensors because they are biocompatible, extraordinarily sensitive, and promising for large-area detection. Electrical and label-free detection of DNA can be achieved by monitoring the conductance change of devices fabricated from these carbon materials. Here, the recent advances in this research area are briefly reviewed. The key issues and perspectives of future development are also discussed.

Understanding effects of molecular adsorption at a single-wall boron nitride nanotube interface from density functional theory calculations

In this paper, we explored computationally the feasibility of modulating the bandgap in a single-wall BN nanotube (BNNT) upon noncovalent adsorption of organic molecules, combined with the application of a transverse electric field. Effects of analytes' physisorption on the surface of BNNTs regarding structural and electronic properties were delineated. Relatively large binding energies were calculated, however, with minimal perturbation of the structural framework. Electronic structure calculations indicated that the bandgap of BNNTs can be modified by weak adsorption due to the presence of adsorbate states in the gap of the host system. Furthermore, we have shown that the application of a transverse electric field can tune the bandgap by shifting adsorbate states, consistent with calculated current-voltage characteristics.

An analysis of the different behavior of DNA and RNA through the study of the mutual relationship between stacking and hydrogen bonding

The mutual relationship between stacking and hydrogen bonding and the possible influence of stacking in the different behavior of DNA and RNA base pairs have been studied through complete DFT optimization of different structures of A-U and A-T dimers (i.e., A-U/U-A and A-T/T-A), using some functionals developed by the group of Truhlar. The results obtained in this work clearly show that stacking and hydrogen bonding are deeply connected. The different behavior of DNA and RNA when replacing uracil by thymine can be interpreted through the formation of a stabilizing CH/pi interaction between the methyl group of thymine and the five-member ring of adenine.

Interaction of substituted aromatic compounds with graphene

We have modeled the adsorption of various substituted derivatives of benzene on a graphene sheet, using a first-principles density functional theory-local density approximation method. The presence of functional groups can significantly alter the overall magnitude of pi-pi interactions between the adsorbed molecules and graphene by giving rise to strong medium-range interactions involving pi-orbitals of the substituents. When the substituents can simultaneously permit the formation of hydrogen bonds between adsorbed molecules, it is possible to evaluate the relative contributions of hydrogen bonding and pi-based interactions to the overall adsorption. Adsorption of individual molecules and hydrogen-bonded aggregates reflects a hierarchical balance of the different interactions that determine the overall energy of adsorption.