A QM:MM model for the interaction of DNA nucleotides with carbon nanotubes

The molecular origin of the electrostatic gating of single-molecule field-effect biosensors investigated by molecular dynamics simulations

Field-effect biosensors (bioFETs) offer a novel way to measure the kinetics of biomolecular events such as protein function and DNA hybridization at the single-molecule level on a wide range of time scales. These devices generate an electrical current whose fluctuations are correlated to the kinetics of the biomolecule under study. BioFETs are indeed highly sensitive to changes in the electrostatic potential (ESP) generated by the biomolecule. Here, using all-atom solvent explicit molecular dynamics simulations, we further investigate the molecular origin of the variation of this ESP for two prototypical cases of proteins or nucleic acids attached to a carbon nanotube bioFET: the function of the lysozyme protein and the hybridization of a 10-nt DNA sequence, as previously done experimentally. Our results show that the ESP changes significantly on the surface of the carbon nanotube as the state of these two biomolecules changes. More precisely, the ESP distributions calculated for these molecular states explain well the magnitude of the conductance fluctuations measured experimentally. The dependence of the ESP with salt concentration is found to agree with the reduced conductance fluctuations observed experimentally for the lysozyme, but to differ for the case of DNA, suggesting that other mechanisms might be at play in this case. Furthermore, we show that the carbon nanotube does not impact significantly the structural stability of the lysozyme, corroborating that the kinetic rates measured using bioFETs are similar to those measured by other techniques. For DNA, we find that the structural ensemble of the single-stranded DNA is significantly impacted by the presence of the nanotube, which, combined with the ESP analysis, suggests a stronger DNA-device interplay. Overall, our simulations strengthen the comprehension of the inner working of field-effect biosensors used for single-molecule kinetics measurements on proteins and nucleic acids.

Molecular insights on the dynamic stability of peptide nucleic acid functionalized carbon and boron nitride nanotubes

The strategic approaches to the design of self-assembled hybrids of biomolecular systems at the nanoscale such as deoxyribonucleic acid (DNA) with single-wall carbon nanotubes (CNTs) and their structural analog, boron nitride nanotubes (BNNTs), rely on understanding how biomolecules recognize and mediate intermolecular interactions with the nanomaterial's surface. In this paper, we consider peptide nucleic acid (PNA), which is a synthetic analog of DNA, and investigate its interaction with a zigzag CNT and BNNT of similar diameter. The results based on the molecular dynamics method find that PNA provides definitive contrasts in the adsorption on the tubular surface in aqueous solution: it prefers to wrap along the circumferential direction on a (11,0) CNT, whereas it binds along the axial direction adopting an extended configuration on a (11,0) BNNT. Moreover, gas-phase Monte Carlo simulations show a dependence of the nanotube diameter on the calculated adsorption energy, with BNNTs exhibiting higher adsorption energy compared to CNTs, and the largest-diameter (25,0) tubular configuration facilitates encapsulation of PNA rather than PNA being adsorbed on its sidewall. The results are expected to be of relevance in the design of novel PNA-based archetypal hybrid materials for nanoscale applications in health-related areas including biosensing.

Nucleotide conjugated (ZnO)3 cluster: Interaction and optical characteristics using TDDFT

Binding of four DNA nucleotide units with (ZnO)₃ cluster in an aqueous phase has been investigated using density functional theory (DFT) and time dependent-density functional theory (TDDFT) method and the stability order for (ZnO)₃-nucleobases/sugar/phosphate systems is predicted as phosphate > C > A > S > T ∼ G. The order of binding energy for (ZnO)₃-nucleotide hybrid systems is observed to be (ZnO)₃ + nuc-C ˃ (ZnO)₃ + nuc-A ˃ (ZnO)₃ + nuc-G ˃ (ZnO)₃ + nuc-T. The binding of nucleotide units with the cluster has been explained on the basis of molecular electrostatic potential (MEP) plots, hydrogen bonding, glycosidic torsion angles, density of state (DOS) plots. The photophysical properties of (ZnO)₃-nucleotide complexes have been studied using TDDFT approach. Among all (ZnO)₃-nucleotide complexes, the absorption spectra of (ZnO)₃ + nuc-A and (ZnO)₃ + nuc-C complexes are seen to undergo red shift with respect to their bare nucleotide units that would be useful in the optical sensing of the respective nucleotides of DNA. It is interesting to note that binding of the nucleotide unit with the cluster makes it fluorescent, the study reports the fluorescence activity of (ZnO)₃ + nuc-T complex.

Electronic properties of carbon nanotubes complexed with a DNA nucleotide

Electronic properties of carbon nanotubes (CNTs) play an important role in their interactions with nano-structured materials. In this work, interactions of adenosine monophosphate (AMP), a DNA nucleotide, with metallic and semi-conducting CNTs are studied using the density functional tight binding (DFTB) method. The electronic structure of semi-conducting CNTs was found to be changed as they turned to metallic CNTs in a vacuum upon interaction with the nucleotide while metallic CNTs remain metallic. Specifically, the band gap of semi-conducting CNTs was decreased by 0.79 eV on average while nearly no change was found in the metallic tubes. However, our investigations showed that the presence of explicit water molecules prevents the metallicity change and only small changes in the CNT band gap occur. According to our charge analysis, the average negative charge accumulated on CNTs upon interaction with the AMP was determined to be 0.77 e in a vacuum while it was 0.03 e in solution. Therefore, it is essential to include explicit water molecules in simulating complexes formed by DNA nucleotides and CNTs which were ignored in several past studies performed using quantum mechanical approaches.

Computational investigation of the α2β1 integrin–collagen triple helix complex interaction

In this paper, quantum biochemistry methods have been used to describe important protein–protein interactions for the complex integrin–collagen.

Quantum mechanical investigation of G-quartet systems of DNA

Minima of the electric field and positions of K⁺ and Na⁺ (zero of the x-coordinate is the center of the cavity).

Evaluating dispersion forces for optimization of van der Waals complexes using a non-empirical functional

Modelling dispersion interactions with traditional density functional theory (DFT) is a challenge that has been extensively addressed in the past decade. The exchange-dipole moment (XDM), among others, is a non-empirical add-on dispersion correction model in DFT. The functional PW86+PBE+XDM for exchange, correlation and dispersion, respectively, compromises an accurate functional for thermochemistry and for van der Waals (vdW) complexes at equilibrium and non-equilibrium geometries. To use this functional in optimizing vdW complexes, rather than computing single point energies, it is necessary to evaluate accurate forces. The purpose of this paper is to validate that, along the potential energy surface, the distance at which the energy is minimum is commensurate with the distance at which the forces vanish to zero. This test was validated for 10 rare gas diatomic molecules using various integration grids and different convergence criteria. It was found that the use of either convergence criterion, 10^-6 or 10^-8, in Gaussian09, does not affect the accuracy of computed optimal distances and binding energies. An ultra-fine grid needs to be used when computing accurate energies using generalized gradient approximation functionals.This article is part of the themed issue 'Multiscale modelling at the physics-chemistry-biology interface'.

Electrostatics of DNA nucleotide-carbon nanotube hybrids evaluated from QM:MM simulations

Biomolecule-functionalized carbon nanotubes (CNTs) have been studied vastly in recent years due to their potential applications for instance in cancer detection, purification and separation of CNTs, and nanoelectronics. Studying the electrostatic potential generated by a biomolecule-CNT hybrid is important in predicting its interactions with the surrounding environment such as charged particles and surfaces. In this paper, we performed atomistic simulations using a QM:MM approach to evaluate the electrostatic potential and charge transfer for a hybrid structure formed by a DNA nucleotide and a CNT in solution. Four types of DNA nucleotides and two CNTs with chiralities of (4,4) and (7,0) were considered. The types of nucleotides and CNTs were both found to play important roles in the electrostatic potential and charge transfer of the hybrid. At the same distance from the CNT axis, the electrostatic potential for the nucleotide-(4,4) CNT hybrids was found to be stronger compared with that for the nucleotide-(7,0) CNT hybrids. Higher electric charge was also shown to be transferred from the DNA nucleotides to the (7,0) CNT compared with the (4,4) CNT. These results correlate with the previous finding that the nucleotides bound more tightly to the (7,0) CNT compared with the (4,4) CNT.