Biogeochemical fingerprinting of magnetotactic bacterial magnetite.

Biominerals are important archives of the presence of life and environmental processes in the geological record. However, ascribing a clear biogenic nature to minerals with nanometer-sized dimensions has proven challenging. Identifying hallmark features of biologically controlled mineralization is particularly important for the case of magnetite crystals, resembling those produced by magnetotactic bacteria (MTB), which have been used as evidence of early prokaryotic life on Earth and in meteorites. We show here that magnetite produced by MTB displays a clear coupled C-N signal that is absent in abiogenic and/or biomimetic (protein-mediated) nanometer-sized magnetite. We attribute the presence of this signal to intracrystalline organic components associated with proteins involved in magnetosome formation by MTB. These results demonstrate that we can assign a biogenic origin to nanometer-sized magnetite crystals, and potentially other biominerals of similar dimensions, using unique geochemical signatures directly measured at the nanoscale. This finding is significant for searching for the earliest presence of life in the Earth's geological record and prokaryotic life on other planets.

When Tautomers Matter: UV-Vis Absorption Spectra of Hypoxanthine in Aqueous Solution from Fully Atomistic Simulations.

The UV-Vis spectrum of the solvated purine derivative Hypoxanthine (HYX) is investigated using the Quantum Mechanics/Fluctuating Charges (QM/FQ) multiscale approach combined with a sampling of configurations through atomistic Molecular Dynamics (MD) simulations. Keto 1H7H and 1H9H tautomeric forms of HYX are the most stable in aqueous solution and form different stable complexes with the surrounding water molecules, ultimately affecting the electronic absorption spectra. The final simulated spectrum resulting from the combination of the individual spectra of tautomers agrees very well with most of the characteristics in the measured spectrum. The importance of considering the effect of the solute tautomers and, in parallel, the contribution of the different solvent arrangements around the solute when modeling spectral properties, is highlighted. In addition, the high quality of the computed spectra leads to suggesting an alternative way for acquiring tautomeric populations from combined computational/experimental spectra.

Encapsulation of charged halogens by the 512 water cage.

This study focuses on the encapsulation of the entire series of halides by the 512 cage of twenty water molecules and on the characterization of water to water and water to anion interactions. State-of-the-art computations are used to determine equilibrium geometries, energy related quantities, and thermal stability towards dissociation and to dissect the nature and strength of intermolecular interactions holding the clusters as stable units. Two types of structures are revealed: heavily deformed cages for F- indicating a preference for microsolvation, and slightly deformed cages for the remaining anions indicating a preference for encapsulation. The primary variable dictating the properties of the clusters is the charge density of the central halide, with the most severe effects observed for the F- case. For the remaining halides, the anion may be safely viewed as a sort of "big electron" with little local disruptive power, enough to affect the network of non-covalent hydrogen bonds in the cage, but not enough to break it. Gibbs energies for dissociation either into cavity and halide or into water molecules and halide suggest that, in a similar way as to methane clathrate, a more weakly bonded complex that has been detected in the gas phase, all halide containing clathrate-like structures should be amenable to experimental detection in the gas phase at moderate temperature and pressure conditions.

Interactions and reactivity in crystalline intermediates of mechanochemical cyclorhodation reactions.

There is experimental evidence that solid mixtures of the rhodium dimer [Cp*RhCl2]2 and benzo[h] quinoline (BHQ) produce two different polymorphic molecular cocrystals called 4α and 4ß under ball milling conditions. The addition of NaOAc to the mixture leads to the formation of the rhodacycle [Cp*Rh-(BHQ)Cl], where the central Rh atom retains its tetracoordinate character. Isolate 4ß reacts with NaOAc leading to the same rhodacycle while isolate 4α does not under the same conditions. We show that the puzzling difference in reactivity between the two cocrystals can be traced back to fundamental aspects of the intermolecular interactions between the BHQ and [Cp*RhCl2]2 fragments in the crystalline environment. To support this view, we report a number of descriptors of the nature and strength of chemical bonds and intermolecular interactions in the extended solids and in a cluster model. We calculate formal quantum mechanical descriptors based on electronic structure, electron density, and binding and interaction energies including an energy decomposition analysis. Without exception, all descriptors point to 4ß being a transient structure higher in energy than 4α with larger local and global electrophilic and nucleophilic powers, a more favorable spatial and energetic distribution of the frontier orbitals, and a more fragile crystal structure.

The Role of Hydrogen Bonding in the Raman Spectral Signals of Caffeine in Aqueous Solution.

The identification and quantification of caffeine is a common need in the food and pharmaceutical industries and lately also in the field of environmental science. For that purpose, Raman spectroscopy has been used as an analytical technique, but the interpretation of the spectra requires reliable and accurate computational protocols, especially as regards the Resonance Raman (RR) variant. Herein, caffeine solutions are sampled using Molecular Dynamics simulations. Upon quantification of the strength of the non-covalent intermolecular interactions such as hydrogen bonding between caffeine and water, UV-Vis, Raman, and RR spectra are computed. The results provide general insights into the hydrogen bonding role in mediating the Raman spectral signals of caffeine in aqueous solution. Also, by analyzing the dependence of RR enhancement on the absorption spectrum of caffeine, it is proposed that the sensitivity of the RR technique could be exploited at excitation wavelengths moderately far from 266 nm, yet achieving very low detection limits in the quantification caffeine content.

UV-Resonance Raman Spectra of Systems in Complex Environments: A Multiscale Modeling Applied to Doxorubicin Intercalated into DNA.

UV-Resonance Raman (RR) spectroscopy is a valuable tool to study the binding of drugs to biomolecular receptors. The extraction of information at the molecular level from experimental RR spectra is made much easier and more complete thanks to the use of computational approaches, specifically tuned to deal with the complexity of the supramolecular system. In this paper, we propose a protocol to simulate RR spectra of complex systems at different levels of sophistication, by exploiting a quantum mechanics/molecular mechanics (QM/MM) approach. The approach is challenged to investigate RR spectra of a widely used chemotherapy drug, doxorubicin (DOX) intercalated into a DNA double strand. The computed results show good agreement with experimental data, thus confirming the reliability of the computational protocol.

Computational Insights into the Adsorption of Ligands on Gold Nanosurfaces.

We study the adsorption process of model peptides, nucleobases, and selected standard ligands on gold through the development of a computational protocol based on fully atomistic classical molecular dynamics (MD) simulations combined with umbrella sampling techniques. The specific features of the interface components, namely, the molecule, the metallic substrate, and the solvent, are taken into account through different combinations of force fields (FFs), which are found to strongly affect the results, especially changing absolute and relative adsorption free energies and trends. Overall, noncovalent interactions drive the process along the adsorption pathways. Our findings also show that a suitable choice of the FF combinations can shed light on the affinity, position, orientation, and dynamic fluctuations of the target molecule with respect to the surface. The proposed protocol may help the understanding of the adsorption process at the microscopic level and may drive the in-silico design of biosensors for detection purposes.

Testing the Influence of Laser Pulse Energy and Rate in the Atom Probe Tomography Analysis of Minerals.

The use of atom probe tomography (APT) for mineral analysis is contributing to fundamental studies in Earth Sciences. Meanwhile, the need for standardization of this technique is becoming evident. Pending the use of mineral standards, the optimization of analysis parameters is needed to facilitate the study of different mineral groups in terms of data collection and quality. The laser pulse rate and energy are variables that highly affect the atom evaporation process occurring during APT analysis, and their testing is important to forecast mineral behavior and obtain the best possible data. In this study, five minerals representative of major groups (albite, As-pyrite, barite, olivine, and monazite) were analyzed over a range of laser pulse energies (10-50 pJ) and rates (100-250 kHz) to assess output parameter quality and evaluate compositional estimate stoichiometry. Among the studied minerals, As-pyrite, with the higher thermal conductivity and lower band gap, was the most affected by the laser pulse variation. Chemical composition estimates equal or close to the general chemical formula were achieved for monazite and As-pyrite. The analysis of multihit events has proved to be the best strategy to verify the efficacy of the evaporation process and to evaluate the best laser pulse setting for minerals.

The Role of Spike Protein Mutations in the Infectious Power of SARS-COV-2 Variants: A Molecular Interaction Perspective.

Specific S477N, N501Y, K417N, K417T, E484K mutations in the receptor binding domain (RBD) of the spike protein in the wild type SARS-COV-2 virus have resulted, among others, in the following variants: B.1.160 (20A or EU2, first reported in continental Europe), B1.1.7 (α or 20I501Y.V1, first reported in the United Kingdom), B.1.351 (ß or 20H/501Y.V2, first reported in South Africa), B.1.1.28.1 (γ or P.1 or 20J/501Y.V3, first reported in Brazil), and B.1.1.28.2 (ζ, or P.2 or 20B/S484K, also first reported in Brazil). From the analysis of a set of bonding descriptors firmly rooted in the formalism of quantum mechanics, including Natural Bond Orbitals (NBO), Quantum Theory of Atoms In Molecules (QTAIM) and highly correlated energies within the Domain Based Local Pair Natural Orbital Coupled Cluster Method (DLPNO-CCSD(T)), and from a set of computed electronic spectral patterns with environmental effects, we show that the new variants improve their ability to recognize available sites to either hydrogen bond or to form salt bridges with residues in the ACE2 receptor of the host cells. This results in significantly improved initial virus⋅⋅⋅cell molecular recognition and attachment at the microscopic level, which trigger the infectious cycle.

Assessing the quality of QM/MM approaches to describe vacuo-to-water solvatochromic shifts.

The performance of different quantum mechanics/molecular mechanics embedding models to compute vacuo-to-water solvatochromic shifts is investigated. In particular, both nonpolarizable and polarizable approaches are analyzed and computed results are compared to reference experimental data. We show that none of the approaches outperform the others and that errors strongly depend on the nature of the molecular transition to be described. Thus, we prove that the best choice of embedding model highly depends on the molecular system and that the use of a specific approach as a black box can lead to significant errors and, sometimes, totally wrong predictions.

Dissecting Bonding Interactions in Cysteine Dimers.

Neutral (n) and zwitterionic (z) forms of cysteine monomers are combined in this work to extensively explore the potential energy surfaces for the formation of cysteine dimers in aqueous environments represented by a continuum. A simulated annealing search followed by optimization and characterization of the candidate structures afforded a total of 746 structurally different dimers held together via 80 different types of intermolecular contacts in 2894 individual non-covalent interactions as concluded from Natural Bond Orbitals (NBO), Quantum Theory of Atoms in Molecules (QTAIM) and Non-Covalent Interactions (NCI) analyses. This large pool of interaction possibilities includes the traditional primary hydrogen bonds and salt bridges which actually dictate the structures of the dimers, as well as the less common secondary hydrogen bonds, exotic X⋯Y (X = C, N, O, S) contacts, and H⋯H dihydrogen bonds. These interactions are not homogeneous but have rather complex distributions of strengths, interfragment distances and overall stabilities. Judging by their Gibbs bonding energies, most of the structures located here are suitable for experimental detection at room conditions.

Ring Vibrations to Sense Anionic Ibuprofen in Aqueous Solution as Revealed by Resonance Raman.

We unravel the potentialities of resonance Raman spectroscopy to detect ibuprofen in diluted aqueous solutions. In particular, we exploit a fully polarizable quantum mechanics/molecular mechanics (QM/MM) methodology based on fluctuating charges coupled to molecular dynamics (MD) in order to take into account the dynamical aspects of the solvation phenomenon. Our findings, which are discussed in light of a natural bond orbital (NBO) analysis, reveal that a selective enhancement of the Raman signal due to the normal mode associated with the C-C stretching in the ring, νC=C, can be achieved by properly tuning the incident wavelength, thus facilitating the recognition of ibuprofen in water samples.

Binding of SARS-CoV-2 to Cell Receptors: A Tale of Molecular Evolution.

The magnified infectious power of the SARS-CoV-2 virus compared to its precursor SARS-CoV is intimately linked to an enhanced ability in the mutated virus to find available hydrogen-bond sites in the host cells. This characteristic is acquired during virus evolution because of the selective pressure exerted at the molecular level. We pinpoint the specific residue (in the virus) to residue (in the cell) contacts during the initial recognition and binding and show that the virus⋅⋅⋅cell interaction is mainly due to an extensive network of hydrogen bonds and to a large surface of noncovalent interactions. In addition to the formal quantum characterization of bonding interactions, computation of absorption spectra for the specific virus⋅⋅⋅cell interacting residues yields significant shifts of Δλmax =47 and 66 nm in the wavelength for maximum absorption in the complex with respect to the isolated host and virus, respectively.

A polarizable three-layer frozen density embedding/molecular mechanics approach.

We present a novel multilayer polarizable embedding approach in which the system is divided into three portions, two of which are treated using density functional theory and their interaction is based on frozen density embedding (FDE) theory, and both also mutually interact with a polarizable classical layer described using an atomistic model based on fluctuating charges (FQ). The efficacy of the model is demonstrated by extending the formalism to linear response properties and applying it to the simulation of the excitation energies of organic molecules in aqueous solution, where the solute and the first solvation shell are treated using FDE, while the rest of the solvent is modeled using FQ charges.

Atom Probe Tomography Analysis of Mica.

Laser-assisted atom probe tomography (APT) is a relatively new, powerful technique for sub-nanometric mineral and biomineral analysis. However, the laser-assisted APT analysis of highly anisotropic and chemically diverse minerals, such as phyllosilicates, may prove especially challenging due to the complex interaction between the crystal structure and the laser pulse upon applying a high electric field. Micas are a representative group of nonswelling clay minerals of relevance to a number of scientific and technological fields. In this study, a Mg-rich biotite was analyzed by APT to generate preliminary data on nonisotropic minerals and to investigate the effect of the crystallographic orientation on mica chemical composition and structure estimation. The difference in results obtained for specimens extracted from the (001) and (hk0) mica surfaces indicate the importance of both experimental parameters and the crystallography. Anisotropy of mica has a strong influence on the physicochemical properties of the mineral during field evaporation and the interpretation of APT data. The promising results obtained in the present study open the way to future innovative APT applications on mica and clay minerals and contribute to the general discussion on the challenges for the analysis of geomaterials by atom probe tomography.

Molecular spectroscopy of aqueous solutions: a theoretical perspective.

Computational spectroscopy is an invaluable tool to both accurately reproduce the spectra of molecular systems and provide a rationalization for the underlying physics. However, the inherent difficulty to accurately model systems in aqueous solutions, owing to water's high polarity and ability to form hydrogen bonds, has severely hampered the development of the field. In this tutorial review we present a technique developed and tested in recent years based on a fully atomistic and polarizable classical modeling of water coupled with a quantum mechanical description of the solute. Thanks to its unparalleled accuracy and versatility, this method can change the perspective of computational and experimental chemists alike.

Theory and algorithms for chiroptical properties and spectroscopies of aqueous systems.

Chiroptical properties and spectroscopies are valuable tools to study chiral molecules and assign absolute configurations. The spectra that result from chiroptical measurements may be very rich and complex, and hide much of their information content. For this reason, the interplay between experiments and calculations is especially useful, provided that all relevant physico-chemical interactions that are present in the experimental sample are accurately modelled. The inherent difficulty associated to the calculation of chiral signals of systems in aqueous solutions requires the development of specific tools, able to account for the peculiarities of water-solute interactions, and especially its ability to form hydrogen bonds. In this perspective we discuss a multiscale approach, which we have developed and challenged to model the most used chiroptical techniques.

Absorption spectra of xanthines in aqueous solution: a computational study.

We present a detailed computational analysis of the UV/Vis spectra of caffeine, paraxanthine and theophylline in aqueous solution. A hierarchy of solvation approaches for modeling the aqueous environment have been tested, ranging from the continuum model to the non-polarizable and polarizable quantum mechanical (QM)/molecular mechanics (MM) models, with and without the explicit inclusion of water molecules in the QM portion. The computed results are directly compared with the experimental data, thus highlighting the role of electrostatic, polarization and hydrogen boding solute-solvent interactions.

Simulating Absorption Spectra of Flavonoids in Aqueous Solution: A Polarizable QM/MM Study.

We present a detailed computational study of the UV/Vis spectra of four relevant flavonoids in aqueous solution, namely luteolin, kaempferol, quercetin, and myricetin. The absorption spectra are simulated by exploiting a fully polarizable quantum mechanical (QM)/molecular mechanics (MM) model, based on the fluctuating charge (FQ) force field. Such a model is coupled with configurational sampling obtained by performing classical molecular dynamics (MD) simulations. The calculated QM/FQ spectra are compared with the experiments. We show that an accurate reproduction of the UV/Vis spectra of the selected flavonoids can be obtained by appropriately taking into account the role of configurational sampling, polarization, and hydrogen bonding interactions.

Interplay between conformational and solvent effects in UV-visible absorption spectra: curcumin tautomers as a case study.

We present a combined theoretical and experimental study on the UV-vis spectra of enol-keto (EK) and keto-keto (KK) tautomeric forms of curcumin dissolved in aqueous solution. Solvent effects have been investigated by resorting to the implicit polarizable continuum model (QM/PCM) and non-polarizable and fully polarizable QM/MM approaches, the latter based on the fluctuating charges (FQ) force-field. In particular, all methods are challenged to rationalize the contribution of conformational, electrostatic and polarization effects in the calculation of the vertical excitation spectra of curcumin tautomers. The obtained results highlight that for both tautomers specific solute-solvent hydrogen-bond interactions play a minor role with respect to conformational and electrostatic effects.