On the Binding Strength Sequence for Nucleic Acid Bases and C(60) with Density Functional and Dispersion-corrected Density Functional Theories: Whether C(60) could protect nucleic acid bases from radiation-induced damage?

Inquest for the interaction of canonical and non-canonical DNA/RNA bases with ternary based 2D Si2BN and doped Si2BN for biosensing applications

Density functional theory (DFT) is invoked to investigate the interaction between the canonical (CN) and non-canonical (NC) bases with pristine Si2BN (Si2BN) and Phosphorous-doped Si2BN (P-dop-Si2BN) sheets. Inquest for the better sensing substrate is decided through the adsorption energy calculation which reveals that doping of phosphorous atom enhances the adsorption strength of AT (-83.74 kcal/mol) AU (-82.77 kcal/mol) and GC (-96.36 kcal/mol) base pairs. The CN and NC bases have higher adsorption energy than the previous reported values which concludes that the P-dop-Si2BN sheet will be optimal substrate to sense the bases. Meanwhile, the selected CN and NC (except hypoxanthine) bases interact with sheet in parallel manner which infers the π-π interaction with Si2BN and P-dop-Si2BN sheets. The energy gap variation (ΔEg%) of the P-dop-Si2BN complexes has a noticeable change, ranging from -24.75 to -197.28% which thrust the sensitivity of the P-dop-Si2BN sheet over the detection of CN and NC bases. The natural population analysis (NPA) and electron density difference map (EDDM) confirms that charges are transferred from CN and NC bases to Si2BN and P-dop-Si2BN sheet. The optical property of the P-dop-Si2BN complexes reveals that the noticeable red and blue shift in the visible and near-infrared regions (778 nm to 1143 nm) has been observed. Therefore, the above results conclude that the P-dop-Si2BN sheet plays a potential candidate to detect the CN and NC bases which contribute to the development of biosensors and DNA/RNA sequencing devices.Communicated by Ramaswamy H. Sarma.

Graphynes and Graphdiynes for Energy Storage and Catalytic Utilization: Theoretical Insights into Recent Advances

Carbon allotropes have contributed to all aspects of people's lives throughout human history. As emerging carbon-based low-dimensional materials, graphyne family members (GYF), represented by graphdiyne, have a wide range potential applications due to their superior physical and chemical properties. In particular, graphdiyne (GDY), as the leader of the graphyne family, has been practically applied to various research fields since it was first successfully synthesized. GYF have a large surface area, both sp and sp² hybridization, and a certain band gap, which was considered to originate from the overlap of carbon 2p_z orbitals and the inhomogeneous π-bonds of carbon atoms in different hybridization forms. These properties mean GYF-based materials still have many potential applications to be developed, especially in energy storage and catalytic utilization. Since most of the GYF have yet to be synthesized and applications of successfully synthesized GYF have not been developed for a long time, theoretical results in various application fields should be shared to experimentalists to attract more intentions. In this Review, we summarized and discussed the synthesis, structural properties, and applications of GYF-based materials from the theoretical insights, hoping to provide different viewpoints and comments.

Microsolvation of lithium cation in xenon clusters: An octahedral growth pattern

The structural and energetic proprieties for the Li ⁺ Xe_n (n = 1-18) clusters are investigated using both Basin-Hopping combined with Potential Model description (BH-PM) and DFT methods. A structural transition from tetrahedral (4 coordination) form to octahedral (6 coordination) one is observed for n = 6. Above this size, all structures have an octahedral core. The cubic-face-centered arrangement for xenon atoms is detected for Li ⁺ Xe₁₄. To the best of our knowledge, the Li ⁺ Xe_n (n = 1-18) clusters are studied in the present work for the first time using the DFT theoretical approach. The M062X functional combined with aug-cc-pVDZ (for Li) and def2-TZVP (for Xe) basis sets reproduces accurately the CCSD(T) potential energy curve of Li + Xe system. Atom-Centered Density Matrix Propagation (ADMP) molecular dynamic calculations have been carried. Moreover, we investigate the larger sizes n = 31-35, 44, and 55 for the first time using the BH-PM theoretical approach. The closing of the first and second octahedron shells are proved for the n = 6 and 34 sizes, respectively. The relative stabilities of the Li ⁺ Xe_n molecules are also studied by computing the total energy, the binding energy per atoms for each size n. Then, the second energy difference between the size n and its two near neighbors allows identifying the magic number series. Our present data are analyzed, discussed and compared with the available theoretical and experimental data.

γ-Graphyne nanotubes as defect-free catalysts of the oxygen reduction reaction: a DFT investigation

γGyNTs as excellent metal-free ORR catalysts without any defects.

Synergy of sp-N and sp2-N codoping endows graphdiyne with comparable oxygen reduction reaction performance to Pt

Nitrogen doped graphdiyne (NGDY) has been reported to have comparable oxygen reduction reaction (ORR) performance to Pt-based catalysts. However, the source of this enhanced ORR performance is not clearly understood. Herein, density functional theory calculations were performed to study the detailed ORR process on NGDY. The theoretically predicted overpotential (η) of GDY materials was 0.442 V, which is comparable to that of Pt-based catalysts, suggesting that GDY is a candidate for non-expensive metal-free ORR catalyst. Our results revealed that the good ORR performance of NGDY originates from the synergy of sp-N and sp²-N, which rules out the experimental proposal that sp-N doping is the dominating factor. Our results further suggest that local positive charge is not a definite descriptor to predict the ORR performance of GDY; instead ΔG_O shows a better correlation with performance. Furthermore, it was revealed that the adsorption site is crucial for determining ORR performance, which should not be ignored to fully understand the catalytic activity of GDY-based materials.

Nitrogen Doping of Carbon Nanoelectrodes for Enhanced Control of DNA Translocation Dynamics

Controlling the dynamics of DNA translocation is a central issue in the emerging nanopore-based DNA sequencing. To address the potential of heteroatom doping of carbon nanostructures and for achieving this goal, herein, we carry out atomistic molecular dynamics simulations for single-stranded DNAs translocating between two pristine or doped carbon nanotube (CNT) electrodes. Specifically, we consider the substitutional nitrogen doping of capped CNT (capCNT) electrodes and perform two types of molecular dynamics simulations for the entrapped and translocating single-stranded DNAs. We find that the substitutional nitrogen doping of capCNTs facilitates and stabilizes the edge-on nucleobase configurations rather than the original face-on ones and slows down the DNA translocation speed by establishing hydrogen bonds between the N dopant atoms and nucleobases. Due to the enhanced interactions between DNAs and N-doped capCNTs, the duration time of nucleobases within the nanogap was extended by up to ∼300%. Given the possibility to be combined with the extrinsic light or gate voltage modulation methods, the current work demonstrates that the substitutional nitrogen doping is a promising direction for the control of DNA translocation dynamics through a nanopore or nanogap, based of carbon nanomaterials.

TDDFT Study of the Optical Excitation of Nucleic Acid Bases-C60 Complexes

The potential of C₆₀ as a nucleic acid base (NAB) optical sensor is theoretically explored. We investigate the adsorption of four NABs, namely, adenine, cytosine, guanine, and thymine, on C₆₀ in the gas phase. For the optimal NAB@C₆₀ adsorption configurations, obtained using a dispersion-corrected density functional, we calculate the vis-near-ultraviolet optical response using time-dependent density functional theory. While the isolated C₆₀ and NAB molecules do not exhibit visible optical excitation, we find that C₆₀/NAB conjugation gives rise to distinct spectral features in the visible range. These results suggest that C₆₀ conjugation can be applied for photodetection of individual NABs.

The effect of defect types on the electronic and optical properties of graphene nanoflakes physisorbed by ionic liquids

Defect engineering and non-covalent interaction strategies allow for dramatically tuning the optoelectronic properties of graphene. Using ab initio density functional theory (M06-2X/cc-pVDZ), we find that the nature of defects on the graphene nanoflakes (GNFs) and the size of defective GNF (DGNF) surfaces affect the binding energy (ΔE_b) of ionic liquids (ILs) and the UV-Vis absorption spectra of DGNFIL complexes. Further, our results indicate that increasing the size of DGNFs affects the geometrical structure of the surfaces and increases the binding energy of ILs by about 10%. Analysis based on AIM and EDA shows that the interactions between ILs and DGNFs are non-covalent in nature (dispersion energy being dominant) and associated with charge transfer between the IL and nanoflakes. A comparison between the ΔE_b values of ILs on DGNFs, GNFs, and h-BN nanoflakes (h-BNNF) shows that the presence of defects on the GNF surfaces increases the binding energy values as follows: DGNFIL > pristine GNFIL > h-BNNFIL. Our calculations indicate that increasing the size of DGNF surfaces leads to a decrease in the HOMO-LUMO energy gap (E_g) of the DGNF surfaces. Orbital energy and density of state calculations show that the E_g of DV(SW)-GNFs decreases upon IL adsorption and their Fermi energy level is shifted depending on the type of IL, thus enabling better conductivity. Reactivity descriptors generally indicate that the chemical potential (μ) and chemical hardness (η) of nanoflakes decrease upon IL adsorption, whereas the electrophilicity index (ω) increases. The UV-Vis absorption spectrum of DV-GNF and SW-GNF shows four bands in the visible spectrum which correspond to π → π* transitions with the absorption bands of SW-GNF appearing at higher wavelengths than those of DV-GNF. The most intense absorption bands in DV-GNF (λ = 348 nm) and SW-GNF (λ = 375 nm) are associated with electronic transitions HOMO-1 → LUMO+2 and HOMO → LUMO+1, respectively. In addition, these absorption bands undergo a red-shift by both increasing the size of the DV(SW)-GNF surfaces and IL adsorption. We also observe that the energy gaps and absorption spectra can be altered by varying the defect types and the type of IL adsorbate, where the defect types affect the spectral shapes of the bands and adsorbates at the first absorption peak, thus having potential application for light-emitting devices.

Interaction Mechanism of Doxorubicin and SWCNT: Protonation and Diameter Effects on the Drug Loading and Releasing

In the present work the adsorption of doxorubicin (DOX) on the surface of single-walled carbon nanotube (SWCNT) as well as its encapsulation in SWCNT, and their dependence on the protonation of NH2 group of DOX, solvent, and the diameter of armchair (n,n) SWCNT were systematically investigated using theoretical methods such as PM6-DH2 and M06-2X in the scheme of OMIOM. It was found that the two loadings, adsorption on the sidewall of CNT and the encapsulation in CNT, have distinct solvent, protonation and diameter dependences. The encapsulation is much stronger than the adsorption of DOX on the sidewall of CNT, and the former also has significantly higher solvent and protonation effects than the latter. The adsorption primarily occurs through π-π stacking and just becomes slightly stronger as the diameter of CNT increases, while besides π-π stacking the additional C-H/N-H/O-H…π and C=O…π also contribute to the encapsulation of DOX in CNT. It seems that (8,8) CNT (diameter ~ 11Å) energetically is an onset for the encapsulation since the encapsulation turns from endothermic to exothermic as the diameter is larger than approximately 11 Å, and the optimal diameter for the encapsulation is 14Å corresponding to (10,10) CNT. Thus for the thick CNT the encapsulation may also play an important role in the loading and releasing for the CNT-based drug delivery system of the DOX.

Substituent Effects in π-Stacking of Histidine on Functionalized-SWNT and Graphene

Adsorptions of histidine on the functionalized (10,0) single-walled carbon nanotube (SWNT) and graphene were investigated using density function theory methods, M05-2x and DFT-D. The results show that the binding of the histidine ring to the functionalized SWNT is weaker than that to the pristine SWNT for both singlet and triplet complexes, regardless of the electron-donating (-OH, -NH₂) or electron-withdrawing (-COOH) character and their attached sites. The present decreased binding is opposite to the well-known enhanced binding in the substituted benzene dimers. Since the atoms of the histidine are distant from the substituent atoms by over 6Å, there would be no direct interaction between histidine and the substituent as in the case of the substituted benzene systems. The decreased binding can be mainly driven by the aromaticity of the functionalized SWNT. The nucleus-independent chemical shift (NICS) index analysis for the functionalized SWNTs in deed shows that local aromaticity of SWNT is decreased because of the electron redistribution induced by functional groups, and the π-π stacking between the histidine ring and functionalized-SWNT is therefore decreased as compared to the pristine SWNT. However, the above trend does not remain for the binding between the histidine and graphene. The binding of the histidine to the functionalized graphene with -OH and -NH₂ is just slightly weaker than that to the pristine graphene, while its binding to COOH-SWNT becomes a little bit stronger.

The role of N7 protonation of guanine in determining the structure, stability and function of RNA base pairs

Ab initio computations and bioinformatics studies reveal that stabilization of some important RNA structural motifs might involve N7 protonation of guanine.

Size dependent electronic band structures of β- and γ-graphyne nanotubes

In the present paper, density functional theory calculations have been implemented by using Dmol³ to study the electronic band structures of β-graphyne nanotubes (βGyNTs) and γ-graphyne nanotubes (γGyNTs).

Insight into the interaction between DNA bases and defective graphenes: covalent or non-covalent

Although some metal clusters and molecules were found to more significantly bind to defective graphenes than to pristine graphenes, exhibiting chemisorptions on defective graphenes, the present investigation shows that the adsorption of DNA bases on mono- and di-vacant defective graphenes does not show much difference from that on pristine graphene, and is still dominantly driven by noncovalent interactions. In the present study the adsorptions of the nucleobases, adenine (A), cytosine (C), guanine, (G), and thymine (T) on pristine and defective graphenes, are fully optimized using a hybrid-meta GGA density functional theory (DFT), M06-2X/6-31G*, and the adsorption energies are then refined with both M06-2X and B97-D/6-311++G**. Graphene is modeled as nano-clusters of C₇₂H₂₄, C₇₁H₂₄, and C₇₀H₂₄ for pristine, mono- and di-vacant defective graphenes, respectively, supplemented by a few larger ones. The result shows that guanine has the maximum adsorption energy in all of the three adsorption systems; and the sequence of the adsorption strength is G>A>T>C on the pristine and di-vacant graphene and G>T>A>C on the mono-vacant graphene. In addition, the binding energies of the DNA bases with the pristine graphene are less than the corresponding ones with di-vacant defective graphene; however, they are greater than those of mono-vacant graphene with guanine and adenine, while it is dramatic that the binding energies of mono-vacant graphene with thymine and cytosine appear larger than those of pristine graphene.

Advancing risk assessment of engineered nanomaterials: application of computational approaches

Nanotechnology that develops novel materials at size of 100nm or less has become one of the most promising areas of human endeavor. Because of their intrinsic properties, nanoparticles are commonly employed in electronics, photovoltaic, catalysis, environmental and space engineering, cosmetic industry and - finally - in medicine and pharmacy. In that sense, nanotechnology creates great opportunities for the progress of modern medicine. However, recent studies have shown evident toxicity of some nanoparticles to living organisms (toxicity), and their potentially negative impact on environmental ecosystems (ecotoxicity). Lack of available data and low adequacy of experimental protocols prevent comprehensive risk assessment. The purpose of this review is to present the current state of knowledge related to the risks of the engineered nanoparticles and to assess the potential of efficient expansion and development of new approaches, which are offered by application of theoretical and computational methods, applicable for evaluation of nanomaterials.

Interaction and protection mechanism between Li@C(60) and nucleic acid bases (NABs): performance of PM6-DH2 on noncovalent interaction of NABs-Li@C60

To discuss the protection mechanism of DNA from radiation as well as assess the performance of PM6-DH2 on noncovalent interactions, the interaction of four nucleic acid bases (NABs) such as adenine (A), cytosine (C), guanine (G), and thymine (T), with Li@C(60) was extensively investigated with the-state-of-art theoretical methods describing noncovalent systems, like M06-2x, PBE-D, and PM6-DH2 methods. In the gas phase, the binding strength of NABs to Li@C(60) from M06-2x decreases in the sequence, G>C>A>T. As dispersion was explicitly included, PBE-D relatively enhances the binding of A and T and corrects the sequence to, G>A>C∼T. PM6-DH2 predicted similar binding energies to those from PBE-D within 0.5 kcal/mol and the same binding sequence, suggesting that the PM6-DH2 method is promising for nano-scale systems. In the aqueous solution, binding of NABs-Li@C(60) is considerably decreased, and the M06-2X and PM6-D methods yield a different sequence from the gas phase, G>A>T>C. The encapsulation of Li atom results in a lower IP for Li@C(60) than those of NABs, and the dominant localization of single-occupied molecular orbital on Li@C(60) moiety of the complexes NABs-Li@C(60) further indicates that an electron would be ejected from Li@C(60) upon radiation and Li@C(60) is therefore able to protect DNA bases from radiation. In addition, it was revealed that Li prefers coordination with the hexagonal ring at Li@C(60) , which clarifies the existing controversy in this respect. Finally, Yang's reduced density gradient approach clearly shows that the weak and strong noncovalent interaction regions in the complexes, NABs-Li@C(60) and (NABs-Li@C(60) )(+).