Unravelling molecular interactions in uracil clusters by XPS measurements assisted by ab initio and tight-binding simulations

Rapid determination of uracil in biological fluids at mercury thin film electrode for early detection of potential 5-fluorouracil toxicity due to dihydropyrimidine dehydrogenase deficiency

Determination of plasma uracil was reported as a method for evaluation of Dihydropyrimidine dehydrogenase (DPD) activity that is highly demanded to ensure the safe administration of 5-fluorouracil (5-FU)-based therapies to cancer patients. This work reports the development of a simple electroanalytical method based on adsorptive stripping square wave voltammetry (AdSWV) at mercury film-coated glassy carbon electrode (MF/GCE) for the highly sensitive determination of uracil in biological fluids that can be used for diagnosis of decreased DPD activity. Due to the formation of the HgII-Uracil complex at the electrode surface, the accuracy of the measurement was not affected by the complicated matrices in biological fluids including human serum, plasma, and urine. The high sensitivity of the developed method results in a low limit of detection (≈1.3 nM) in human plasma samples, falling below the practical cut-off level of 15 ng mL-1 (≈0.14 μM). This threshold concentration is crucial for predicting 5-FU toxicity, as reported in buffer, and ≤1.15% in biological samples), and accuracy (recovery percentage close to 100%).

Core Binding Energy Calculations: A Scalable Approach with the Quantum Embedding-Based Equation-of-Motion Coupled-Cluster Method

We investigated the use of density matrix embedding theory to facilitate the computation of core ionization energies (IPs) of large molecules at the equation-of-motion coupled-cluster singles doubles with perturbative triples (EOM-CCSD*) level in combination with the core-valence separation (CVS) approximation. The unembedded IP-CVS-EOM-CCSD* method with a triple-ζ basis set produced ionization energies within 1 eV of experiment with a standard deviation of ∼0.2 eV for the core65 data set. The embedded variant contributed very little systematic error relative to the unembedded method, with a mean unsigned error of 0.07 eV and a standard deviation of ∼0.1 eV, in exchange for accelerating the calculations by many orders of magnitude. By employing embedded EOM-CC methods, we computed the core ionization energies of the uracil hexamer, doped fullerene, and chlorophyll molecule, utilizing up to ∼4000 basis functions within 1 eV from experimental values. Such calculations are not currently possible with the unembedded EOM-CC method.

A study of the valence photoelectron spectrum of uracil and mixed water-uracil clusters

The valence ionization of uracil and mixed water-uracil clusters has been studied experimentally and by ab initio calculations. In both measurements, the spectrum onset shows a red shift with respect to the uracil molecule, with the mixed cluster characterized by peculiar features unexplained by the sum of independent contributions of the water or uracil aggregation. To interpret and assign all the contributions, we performed a series of multi-level calculations, starting from an exploration of several cluster structures using automated conformer-search algorithms based on a tight-binding approach. Ionization energies have been assessed on smaller clusters via a comparison between accurate wavefunction-based approaches and cost-effective DFT-based simulations, the latter of which were applied to clusters up to 12 uracil and 36 water molecules. The results confirm that (i) the bottom-up approach based on a multilevel method [Mattioli et al. Phys. Chem. Chem. Phys. 23, 1859 (2021)] to the structure of neutral clusters of unknown experimental composition converges to precise structure-property relationships and (ii) the coexistence of pure and mixed clusters in the water-uracil samples. A natural bond orbital (NBO) analysis performed on a subset of clusters highlighted the special role of H-bonds in the formation of the aggregates. The NBO analysis yields second-order perturbative energy between the H-bond donor and acceptor orbitals correlated with the calculated ionization energies. This sheds light on the role of the oxygen lone-pairs of the uracil CO group in the formation of strong H-bonds, with a stronger directionality in mixed clusters, giving a quantitative explanation for the formation of core-shell structures.

A systematic study of the valence electronic structure of cyclo(Gly-Phe), cyclo(Trp-Tyr) and cyclo(Trp-Trp) dipeptides in the gas phase

The electronic energy levels of cyclo(glycine-phenylalanine), cyclo(tryptophan-tyrosine) and cyclo(tryptophan-tryptophan) dipeptides are investigated with a joint experimental and theoretical approach. Experimentally, valence photoelectron spectra in the gas phase are measured using VUV radiation. Theoretically, we first obtain low-energy conformers through an automated conformer-rotamer ensemble sampling scheme based on tight-binding simulations. Then, different first principles computational schemes are considered to simulate the spectra: Hartree-Fock (HF), density functional theory (DFT) within the B3LYP approximation, the quasi-particle GW correction, and the quantum-chemistry CCSD method. Theory allows assignment of the main features of the spectra. A discussion on the role of electronic correlation is provided, by comparing computationally cheaper DFT scheme (and GW) results with the accurate CCSD method.

Water-biomolecule clusters studied by photoemission spectroscopy and multilevel atomistic simulations: hydration or solvation?

The properties of mixed water-uracil nanoaggregates have been probed by core electron-photoemission measurements to investigate supramolecular assembly in the gas phase driven by weak interactions. The interpretation of the measurements has been assisted by multilevel atomistic simulations, based on semi-empirical tight-binding and DFT-based methods. Our protocol established a positive-feedback loop between experimental and computational techniques, which has enabled a sound and detailed atomistic description of such complex heterogeneous molecular aggregates. Among biomolecules, uracil offers interesting and generalized skeletal features; its structure encompasses an alternation of hydrophilic H-bond donor and acceptor sites and hydrophobic moieties, typical in biomolecular systems, that induces a supramolecular core-shell-like organization of the mixed clusters with a water core and an uracil shell. This structure is far from typical models of both solid-state hydration, with water molecules in defined positions, or liquid solvation, where disconnected uracil molecules are completely surrounded by water.

Revisiting stacking interactions in tetrathiafulvalene and selected derivatives using tight-binding quantum chemical calculations and local coupled-cluster method

The engineering of supramolecular architectures needs accurate descriptions of the intermolecular interactions in crystal structures. Tetrathiafulvalene (TTF) is an effective building block used in the construction of promising functional materials. The parallel packing of the neutral TTF-TTF system was studied previously using the high-level quantum chemical method, advancing it as a valuable model system. The recently developed tight-binding quantum chemical method GFN2-xTB and local coupled-cluster method DLPNO-CCSD(T) were used to investigate the stacking interactions of TTF and selected derivatives deposited in the Cambridge Structural Database. Using the interaction energy of the TTF-TTF dimer calculated at the CCSD(T)/CBS level as the reference, the accuracies of the two methods are investigated. The energy decomposition analysis within the DLPNO-CCSD(T) framework reveals the importance of dispersion interaction in the TTF-related stacking systems. The dispersion interaction density plot vividly shows the magnitude and distribution of the dispersion interaction, providing a revealing insight into the stacking interactions in crystal structures. The results show that the GFN2-xTB and DLPNO-CCSD(T) methods could achieve accuracy at an affordable computational cost, which would be valuable in understanding the nature of parallel stacking in supramolecular systems.

Ultraefficiently Calming Cytokine Storm Using Ti3C2T x MXene

During the global outbreak of COVID-19 pandemic, "cytokine storm" conditions are regarded as the fatal step resulting in most mortality. Hemoperfusion is widely used to remove cytokines from the blood of severely ill patients to prevent uncontrolled inflammation induced by a cytokine storm. This article discoveres, for the first time, that 2D Ti3C2T x MXene sheet demonstrates an ultrahigh removal capability for typical cytokine interleukin-6. In particular, MXene shows a 13.4 times higher removal efficiency over traditional activated carbon absorbents. Molecular-level investigations reveal that MXene exhibits a strong chemisorption mechanism for immobilizing cytokine interleukin-6 molecules, which is different from activated carbon absorbents. MXene sheet also demonstrates excellent blood compatibility without any deleterious side influence on the composition of human blood. This work can open a new avenue to use MXene sheets as an ultraefficient hemoperfusion absorbent to eliminate the cytokine storm syndrome in treatment of severe COVID-19 patients.