Noncovalent Interactions in Paired DNA Nucleobases Investigated by Terahertz Spectroscopy and Solid-State Density Functional Theory

A new horizon for neuroscience: terahertz biotechnology in brain research

Terahertz biotechnology has been increasingly applied in various biomedical fields and has especially shown great potential for application in brain sciences. In this article, we review the development of terahertz biotechnology and its applications in the field of neuropsychiatry. Available evidence indicates promising prospects for the use of terahertz spectroscopy and terahertz imaging techniques in the diagnosis of amyloid disease, cerebrovascular disease, glioma, psychiatric disease, traumatic brain injury, and myelin deficit. In vitro and animal experiments have also demonstrated the potential therapeutic value of terahertz technology in some neuropsychiatric diseases. Although the precise underlying mechanism of the interactions between terahertz electromagnetic waves and the biosystem is not yet fully understood, the research progress in this field shows great potential for biomedical noninvasive diagnostic and therapeutic applications. However, the biosafety of terahertz radiation requires further exploration regarding its two-sided efficacy in practical applications. This review demonstrates that terahertz biotechnology has the potential to be a promising method in the field of neuropsychiatry based on its unique advantages.

An intelligent sensing platform for detecting and identifying biochemical substances based on terahertz spectra

This paper presents the development of an intelligent sensing platform dedicated to accurately identifying terahertz (THz) spectra obtained from various biochemical substances. The platform currently has two distinct identification modes, which focus on identifying five amino acids, namely phenylalanine, methionine, lysine, leucine, and threonine, and five carbohydrates, namely aspartame, fructose, glucose, lactose monohydrate, and sucrose based on their THz spectra. The first mode, called One-dimensional THz Spectrum Identification (OTSI), combines THz time-domain spectroscopy (THz-TDS) with the proposed mini convolutional neural network (MCNN) model. THz-TDS detects biochemical substances, while the MCNN model identifies the THz spectra. The MCNN model has a simple structure and only needs to deal with the THz absorption coefficients of biochemical substances, which are less computationally intensive and easily converged. The model can achieve 99.07 % accuracy in identifying one-dimensional THz spectra of the ten biochemical substances. The second mode, THz Spectrum Image-based Identification (TSII), applies the YOLO-v5 target detection model to THz spectral image recognition. The YOLO-v5 model uses THz absorption peaks as identification features and can identify biochemical substances based on only one or several THz absorption peaks. The overall identifying accuracy of the YOLO-v5 model for ten biochemical substances is 96.20 %. We also compared the MCNN and YOLO-v5 models with other deep learning and machine learning models, which demonstrate that they have better performance. This feature broadens the platform's utility in biomolecular analysis and paves the way for further research and development in detecting and analyzing diverse biological compounds.

Study of temperature-dependent terahertz spectra of isonicotinamide in the form I using the quasi-harmonic approximation

Anharmonicity of molecular vibrational motions is closely associated with the thermal property of crystals. However, the origin of anharmonicity is still not fully understood. Low-frequency vibrations, which are usually defined in the terahertz (THz) range, show excellent sensitivity to anharmonicity. In this work, anharmonicity of isonicotinamide in the form I was investigated by using temperature-dependent terahertz time-domain spectroscopy and quasi-harmonic approximation (QHA) approach at PBE-D3 and PBE-MBD levels. Both DFT calculations suggest the variation of p-p stacking conformation dominates in the thermal expansion of unit cell. The frequency shifts of the modes in THz range obtained by QHA approach are found to be qualitatively consistent with the observation, demonstrating QHA approach is a useful tool for the interpretation of frequency shifts of modes induced by temperature.

Terahertz spectra and weak intermolecular interactions of nucleosides or nucleoside drugs

In this paper, terahertz (THz) spectra of four DNA nucleosides (Adenosine, Thymidine, Cytidine and Guanosine) and two nucleoside derivatives (Ribavirin and Entecavir, first time reported) in the solid phase were studied experimentally by Fourier Transform Infrared Spectroscopy (FTIR) in the frequency of 1-10 THz. The lattice energy, geometric structure, vibration spectrum of them were analyzed theoretically by the generalized energy-based fragmentation approach under periodic boundary conditions (denoted as PBC-GEBF) and the density functional theory (DFT). The intra- and inter-molecular weak interactions corresponding to the vibrational modes of the crystal, polymer and monomer were obtained, with the help of the potential energy distribution (PED) and reduced density gradient (RDG) methods. It was found that the sum of electronic and thermal free energies increased from the monomer to polymer, and from the polymer to crystal. For example, the inter-molecular interaction energy from the monomer to dimer of adenosine increased 6.969 kcal/mol, and that from the dimer to crystal (the periodic boundary conditions were considered) increased 666.792 kcal/mol. Therefore, only the crystal structure constrained the periodic boundary conditions could well describe the experimental results, although the former scholars chose the monomer or polymer as the initial configuration due to the limitation of computing resources and methods. In THz band, the vibrational modes were generally originated from the collective vibration (more than 99% of them were vibration, only less than 1% of them were rotation and translation) of all molecules involved, which could reflect the molecular structure and spatial distribution of different substances. In order to accurately identify the spectra, we studied the location, type and contribution of all weak interactions, and found that the strong characteristic peaks corresponding to the strong hydrogen bonds came from inter-molecular, while the weak hydrogen bonds mainly originated from intra- and inter-molecular, the out-of-plane bending made the largest contribution, accounting for more than 90%. Furthermore, taking guanine, guanosine and two guanosine derivatives (Ribavirin and Entecavir) as examples, the differences of weak interaction among them caused by different molecular configuration, arrangement and substituent position were studied, and the fundamental reason of THz spectrum change was found. This research can lay a foundation for crystal engineering, supramolecular chemistry, molecular recognition and self-assembly, protein-ligand interaction, etc.

Probing the Hydrogen-Bonding Environment of Individual Bases in DNA Duplexes with Isotope-Edited Infrared Spectroscopy

Measuring the strength of the hydrogen bonds between DNA base pairs is of vital importance for understanding how our genetic code is physically accessed and recognized in cells, particularly during replication and transcription. Therefore, it is important to develop probes for these key hydrogen bonds (H-bonds) that dictate events critical to cellular function, such as the localized melting of DNA. The vibrations of carbonyl bonds are well-known probes of their H-bonding environment, and their signals can be observed with infrared (IR) spectroscopy. Yet, pinpointing a single bond of interest in the complex IR spectrum of DNA is challenging due to the large number of carbonyl signals that overlap with each other. Here, we develop a method using isotope editing and infrared (IR) spectroscopy to isolate IR signals from the thymine (T) C2=O carbonyl. We use solvatochromatic studies to show that the TC2=O signal’s position in the IR spectrum is sensitive to the H-bonding capacity of the solvent. Our results indicate that C2=O of a single T base within DNA duplexes experiences weak H-bonding interactions. This finding is consistent with the existence of a third, noncanonical CH···O H-bond between adenine and thymine in both Watson–Crick and Hoogsteen base pairs in DNA.

Crystallographic evidence of Watson-Crick connectivity in the base pair of anionic adenine with thymine

Utilizing an ionic liquid strategy, we report crystal structures of salts of free anionic nucleobases and base pairs previously studied only computationally and in the gas phase. Reaction of tetrabutylammonium ([N₄₄₄₄]⁺) or tetrabutylphosphonium ([P₄₄₄₄]⁺) hydroxide with adenine (HAd) and thymine (HThy) led to hydrated salts of deprotonated adenine, [N₄₄₄₄][Ad]·2H₂O, and thymine, [P₄₄₄₄][Thy]·2H₂O, as well as the double salt cocrystal, [P₄₄₄₄]₂[Ad][Thy]·3H₂O·2HThy. The cocrystal includes the anionic [Ad^-(HThy)] base pair which is a stable formation in the solid state that has previously not even been suggested. It exhibits Watson-Crick connectivity as found in DNA but which is unusual for the free neutral base pairs. The stability of the observed anionic bases and their supramolecular formations and hydrates has also been examined by electronic structure calculations, contributing to more insight into how base pairs can bind when a proton is removed and highlighting mechanisms of stabilization or chemical transformation in the DNA chains.

DFT Computed Dielectric Response and THz Spectra of Organic Co-Crystals and Their Constituent Components

Terahertz (THz) spectroscopy has been put forth as a non-contact, analytical probe to characterize the intermolecular interactions of biologically active molecules, specifically as a way to understand, better develop, and use active pharmaceutical ingredients. An obstacle towards fully utilizing this technique as a probe is the need to couple features in the THz regions to specific vibrational modes and interactions. One solution is to use density functional theory (DFT) methods to assign specific vibrational modes to signals in the THz region, coupling atomistic insights to spectral features. Here, we use open source planewave DFT packages that employ ultrasoft pseudopotentials to assess the infrared (IR) response of organic compounds and complex co-crystal formulations in the solid state, with and without dispersion corrections. We compare our DFT computed lattice parameters and vibrational modes to experiment and comment on how to improve the agreement between theory and modeling to allow for THz spectroscopy to be used as an analytical probe in complex biologically relevant systems.

THz vibrational spectroscopy for RNA basepair cocrystals and oligonucleotide sequences

Fourier infrared spectrometer and X-ray diffractometer were used to detect the spectra of lab-made U:A (uracil and adenine hydrate dried at room temperature), lab-made C:G (cytosine and guanine hydrate dried at room temperature), U + A (mixture of uracil and adenine), and C + G (mixture of cytosine and guanine). The results of our testing showed that U:A did not form a eutectic, but C:G did. In order to further characterize the vibrational modes of RNA base pair crystals, the absorption spectra of 1‑methyl‑5‑bromouracil‑9‑methyladenine (abbreviated as 1M5BU:9MA) and 1‑methylcytosine‑9‑ethylguanine (abbreviated as 1MC:9EG) were calculated based on the PBC-GEBF (generalized energy-based fragmentation approach under periodic boundary conditions) method. To further study the effect of substituents on the vibrational mode of the crystal structure, the substituents of 1M5BU:9MA and 1MC:9EG were artificially removed. The results of calculation brought out that methyl and ethyl as substituents have little effect on the vibrational spectrum, but halogen atoms such as Br atom in 1M5BU:9MA has a certain influence on the spectrum. Furthermore, THz (terahertz) spectra of the RNA nucleotide sequence 5'-AUCG-3' was analyzed. In the perspective of the spectra with DNA signatures, their spectra show a great deal of similarity regardless of RNA or DNA, or the base sequence difference. This study will provide a very important information for revealing the role of RNA in the formation of biological macromolecules and its mechanism.

THz spectra and corresponding vibrational modes of DNA base pair cocrystals and polynucleotides

The generalized energy-based fragmentation (GEBF) approach has been applied to study the THz spectra and vibrational modes of base pair cocrystals under periodic boundary conditions (denoted as PBC-GEBF). Results of vibrational mode reveal that hydrogen bonds play a pivotal role in the pairing process of base crystals, where most NH and CH bonds stretch to some extent. We also found that hydrogen bonds of a self-made A:T cocrystal completely break in a transition from liquid to the solid state, while self-made C:G cocrystal is different and easier to form a cocrystal, as confirmed by X-ray diffraction (XRD) and terahertz (THz) spectra. Furthermore, we have studied DNA polynucleotides (in both A and B forms) found that the vibrational modes changed a lot during the process of their forming double strand. Despite the key role played by hydrogen bonds, the key contribution originates from collective motions of the main skeleton. A comparative study of the spectra of some stranded fragments suggests that different sequences or forms have similar spectra in THz band. They distinguish from each other mainly in the low-frequency regions, especially below 1 THz. This study would make great contributions to the molecular dynamics model based DNA long-chain structure simulation in the future study.

Intermolecular vibrational modes and H-bond interactions in crystalline urea investigated by terahertz spectroscopy and theoretical calculation

The characteristic absorption spectra of crystalline urea in 0.6-1.8 THz region have been measured by terahertz time-domain spectroscopy at room temperature experimentally. Five broad absorption peaks were observed at 0.69, 1.08, 1.27, 1.47 and 1.64 THz respectively. Moreover, density functional theory (DFT) calculation has been performed for the isolated urea molecule, and there is no infrared intensity in the region below 1.8 THz. This means that single molecule calculations are failure to predict the experimental spectra of urea crystals. To simulate these spectra, calculations on a cluster of seven urea molecules using M06-2X and B3LYP-D3 are performed, and we found that M06-2X perform better. The observed THz vibrational modes are assigned to bending and torsional modes related to the intermolecular H-bond interactions with the help of potential energy distribution (PED) method. Using the reduced-density-gradient (RDG) analysis, the positions and types of intermolecular H-bond interactions in urea crystals are visualized. Therefore, we can confirm that terahertz spectroscopy can be used as an effective means to detect intermolecular H-bond interactions in molecular crystals.

Intermolecular hydrogen bond stretching vibrations observed in terahertz spectra of crystalline vitamins

Intermolecular hydrogen bond stretching vibrations in the crystal of vitamins has several specific properties in the THz spectrum.

Ultrasensitive terahertz sensing of gold nanoparticles inside nano slot antennas

We introduce a robust control method of terahertz (THz) transmission by tuning filling factors of Au nanoparticles (AuNPs) inside nano slot antennas. AuNPs in sub-100 nm diameters were spread over the nano slot antennas, followed by sweeping them into the slots. AuNPs can be efficiently localized and inserted into nano slots where the THz fields are greatly enhanced, by a "squeegee" made of the polydimethylsiloxane (PDMS). The sweeping of the AuNPs results in further dramatic reduction of THz transmission by suppressing the fundamental resonance mode of the nano slot, as compared to a typical random dropping case. It definitely works for an accurate THz transmission control, as well as the removal of unwanted ions that occasionally confuse signal accuracy from the target signals. Our approach provides a complete reinterpretation of sample deposition for further steady demands in developing ultrasensitive terahertz (THz) molecule sensors.

Intramolecular vibrations in low-frequency normal modes of amino acids: L-alanine in the neat solid state

This paper presents a theoretical analysis of the low-frequency phonons of L-alanine by using the solid-state density functional theory at the Γ point. We are particularly interested in the intramolecular vibrations accessing low-frequency phonons via harmonic coupling with intermolecular vibrations. A new mode-analysis method is introduced to quantify the vibrational characteristics of such intramolecular vibrations. We find that the torsional motions of COO(-) are involved in low-frequency phonons, although COO(-) is conventionally assumed to undergo localized torsion. We also find the broad distributions of intramolecular vibrations relevant to important functional groups of amino acids, e.g., the COO(-) and NH3(+) torsions, in the low-frequency phonons. The latter finding is illustrated by the concept of frequency distribution of vibrations. These findings may lead to immediate implications in other amino acid systems.

Terahertz time-domain and low-frequency Raman spectroscopy of organic materials

With the ongoing proliferation of terahertz time-domain instrumentation from semiconductor physics into applied spectroscopy over the past decade, measurements at terahertz frequencies (1 THz ≡ 10(12) Hz ≡ 33 cm(-1)) have attracted a sustained growing interest, in particular the investigation of hydrogen-bonding interactions in organic materials. More recently, the availability of Raman spectrometers that are readily able to measure in the equivalent spectral region very close to the elastic scattering background has also grown significantly. This development has led to renewed efforts in performing spectroscopy at the interface between dielectric relaxation phenomena and vibrational spectroscopy. In this review, we briefly outline the underlying technology, the physical phenomena governing the light-matter interaction at terahertz frequencies, recent examples of spectroscopic studies, and the current state of the art in assigning spectral features to vibrational modes based on computational techniques.

Molecular recognition and interaction between uracil and urea in solid-state studied by terahertz time-domain spectroscopy

Using terahertz time-domain spectroscopy characterization, we observe that urea is able to recognize and interact with uracil efficiently even in the solid phase without involving water or solvents. A cocrystal configuration linked by a pair of hydrogen bonds between uracil and urea was formed. The terahertz absorption spectrum of the cocrystal shows a distinct new absorption at 0.8 THz (26.7 cm(-1)), which originates from the intermolecular hydrogen bonding. Both mechanical milling and heating can accelerate the reaction efficiently. Density functional theory was adopted to simulate the vibrational modes of the cocrystal, and the results agree well with the experimental observation. Multiple techniques, including powder X-ray diffraction, scanning electron microscopy, and Fourier transform infrared spectroscopy, were performed to investigate the reaction process, and they presented supportive evidence. This work enables in-depth understanding of recognition and interaction of urea with nucleobases and comprehension of the denaturation related to RNA. We also demonstrate that terahertz spectroscopy is an effective and alternative tool for online measurement and quality control in pharmaceutical and chemical industries.

Revealing the true crystal structure of L-phenylalanine using solid-state density functional theory

Solid-state density functional theory can be used for crystal structure determination from powder X-ray diffraction data of molecular crystals that are too large and complex for conventional refinement methods.

Noncovalent interactions between modified cytosine and guanine DNA base pair mimics investigated by terahertz spectroscopy and solid-state density functional theory

Modified cytosine and guanine nucleobases cocrystallize in a hydrogen bonding configuration similar to that observed in native DNA. The noncovalent interactions binding these base pairs in the crystalline solid were investigated using terahertz (THz) spectroscopy and solid-state density functional theory (DFT). While stronger hydrogen bonding interactions are responsible for the general molecular orientations in the crystalline state, it is the weaker dipole-dipole and dispersion forces that determine the overall packing arrangement. The inclusion of dispersion interactions in the DFT calculations was found to be necessary to accurately simulate the unit cell structure and THz vibrational spectrum. Using properly modeled intermolecular potentials, the lattice vibrational motions of the cytosine and guanine derivatives were calculated. The vibrational characters of the modes exhibited by the DNA base pair mimic in the THz region were primarily rotational motions and are indicative of the energies and the nature of vibrations that would likely be observed between similar base pairs in DNA molecules.