Toward the identification of cyano-astroCOMs via vibrational features: benzonitrile as a test case.

The James Webb Space Telescope (JWST) opened a new era for the identification of molecular systems in the interstellar medium (ISM) by vibrational features. One group of molecules of increasing interest is cyano-derivatives of aromatic organic molecules, which have already been identified in the ISM on the basis of the analysis of rotational signatures, and so, are plausible candidates for the detection by the JWST. Benzonitrile considered in this work represents a suitable example for the validation of a computational strategy, which can be further applied for different, larger, and not-yet observed molecules. For this purpose, anharmonic simulations of infrared (IR) spectra have been compared with recent FTIR experimental studies. The anharmonic computations using the generalized second-order vibrational perturbation theory (GVPT2) in conjunction with a hybrid force field combining the harmonic part of revDSD-PBEP86-D3/jun-cc-pVTZ with anharmonic corrections from B3LYP-D3/SNSD show very good agreement with those in the experiment, with a mean error of 11 c m - 1 for all fundamental transitions overall and only 2 c m - 1 for the C ≡ N stretching fundamental at 4.49 µ m . The inclusion of overtones up to three-quanta transitions also allowed the prediction of spectra in the near-infrared region, which shows distinct features due to C ≡ N overtones at the 2.26 µ m and 1.52 µ m . The remarkable accuracy of the GVPT2 results opens a pathway for the reliable prediction of spectra for a broader range of cyano-astroCOMs.

AQuaRef: Machine learning accelerated quantum refinement of protein structures.

Cryo-EM and X-ray crystallography provide crucial experimental data for obtaining atomic-detail models of biomacromolecules. Refining these models relies on library-based stereochemical restraints, which, in addition to being limited to known chemical entities, do not include meaningful noncovalent interactions relying solely on nonbonded repulsions. Quantum mechanical (QM) calculations could alleviate these issues but are too expensive for large molecules. We present a novel AI-enabled Quantum Refinement (AQuaRef) based on AIMNet2 neural network potential mimicking QM at substantially lower computational costs. By refining 41 cryo-EM and 30 X-ray structures, we show that this approach yields atomic models with superior geometric quality compared to standard techniques, while maintaining an equal or better fit to experimental data.

Investigating the stability of aromatic carboxylic acids in hydrated magnesium sulfate under UV irradiation to assist detection of organics on Mars.

The Scanning Habitable Environments with Raman and Luminescence for Organics and Chemicals (SHERLOC) instrument onboard the Mars 2020 Perseverance rover detected so far some of the most intense fluorescence signals in association with sulfates analyzing abraded patches of rocks at Jezero crater, Mars. To assess the plausibility of an organic origin of these signals, it is key to understand if organics can survive exposure to ambient Martian UV after exposure by the Perseverance abrasion tool and prior to analysis by SHERLOC. In this work, we investigated the stability of organo-sulfate assemblages under Martian-like UV irradiation and we observed that the spectroscopic features of phthalic and mellitic acid embedded into hydrated magnesium sulfate do not change for UV exposures corresponding to at least 48 Martian sols and, thus, should still be detectable in fluorescence when the SHERLOC analysis takes place, thanks to the photoprotective properties of magnesium sulfate. In addition, different photoproduct bands diagnostic of the parent carboxylic acid molecules could be observed. The photoprotective behavior of hydrated magnesium sulfate corroborates the hypothesis that sulfates might have played a key role in the preservation of organics on Mars, and that the fluorescence signals detected by SHERLOC in association with sulfates could potentially arise from organic compounds.

Harmonic and anharmonic vibrational computations for biomolecular building blocks: Benchmarking DFT and basis sets by theoretical and experimental IR spectrum of glycine conformers.

Advanced vibrational spectroscopic experiments have reached a level of sophistication that can only be matched by numerical simulations in order to provide an unequivocal analysis, a crucial step to understand the structure-function relationship of biomolecules. While density functional theory (DFT) has become the standard method when targeting medium-size or larger systems, the problem of its reliability and accuracy are well-known and have been abundantly documented. To establish a reliable computational protocol, especially when accuracy is critical, a tailored benchmark is usually required. This is generally done over a short list of known candidates, with the basis set often fixed a priori. In this work, we present a systematic study of the performance of DFT-based hybrid and double-hybrid functionals in the prediction of vibrational energies and infrared intensities at the harmonic level and beyond, considering anharmonic effects through vibrational perturbation theory at the second order. The study is performed for the six-lowest energy glycine conformers, utilizing available "state-of-the-art" accurate theoretical and experimental data as reference. Focusing on the most intense fundamental vibrations in the mid-infrared range of glycine conformers, the role of the basis sets is also investigated considering the balance between computational cost and accuracy. Targeting larger systems, a broad range of hybrid schemes with different computational costs is also tested.

Undecanoic Acid and L-Phenylalanine in Vermiculite: Detection, Characterization, and UV Degradation Studies for Biosignature Identification on Mars.

Solar radiation that arrives on the surface of Mars interacts with organic molecules present in the soil. The radiation can degrade or transform the organic matter and make the search for biosignatures on the planet's surface difficult. Therefore, samples to be analyzed by instruments on board Mars probes for molecular content should be selectively chosen to have the highest organic preservation content. To support the identification of organic molecules on Mars, the behavior under UV irradiation of two organic compounds, undecanoic acid and L-phenylalanine, in the presence of vermiculite and two chloride salts, NaCl and MgCl, was studied. The degradation of the molecule's bands was monitored through IR spectroscopy. Our results show that, while vermiculite acts as a photoprotective mineral with L-phenylalanine, it catalyzes the photodegradation of undecanoic acid molecules. On the other hand, both chloride salts studied decreased the degradation of both organic species acting as photoprotectors. While these results do not allow us to conclude on the preservation capabilities of vermiculite, they show that places where chloride salts are present could be good candidates for in situ analytic experiments on Mars due to their organic preservation capacity under UV radiation.

Isomerization-induced fluorescence enhancement of two new viologen derivatives: mechanism insight and DFT calculations.

The dark-colored viologen radical cations are unstable in air and easily fade, thus greatly limiting their applications. If a suitable substituent is introduced into the structure, it will have the dual function of chromism and luminescence, which will broaden its application field. Here, Vio1·2Cl and Vio2·2Br were synthesized by introducing aromatic acetophenone and naphthophenone substituents into the viologen structure. The keto group (-CH2CO-) on the substituents is prone to isomerize into the enol structure (-CH[double bond, length as m-dash]COH-) in organic solvents, especially in DMSO, resulting in a larger conjugated system to stabilize the molecular structure and enhance fluorescence. The time-dependent fluorescence spectrum shows obvious keto-to-enol isomerization-induced fluorescence enhancement. The quantum yield also increased significantly (T = 1 day, ΦVio1 = 25.81%, ΦVio2 = 41.44%; T = 7 days, ΦVio1 = 31.48%, and ΦVio2 = 54.40%) in DMSO. The NMR and ESI-MS data at different times further confirmed that the fluorescence enhancement was caused by isomerization, and no other fluorescent impurities were produced in solution. DFT calculations show that the enol form is almost coplanar throughout the molecular structure, which is conducive to stabilizing the structure and enhancing fluorescence. The fluorescence emission peaks of the keto and enol structures of Vio12+ and Vio22+ were at 416-417 nm and 563-582 nm, respectively. The fluorescence relative oscillator strength of Vio12+ and Vio22+ enol structures is significantly higher than that of keto structures (f value changes from 1.53 to 2.63 for Vio12+ and from 1.62 to 2.81 for Vio22+), indicating stronger fluorescence emission of the enol structure. The calculated results are in good agreement with the experimental results. Vio1·2Cl and Vio2·2Br are the first examples of isomerization-induced fluorescence enhancement of viologen derivatives, which shows strong solvatofluorochromism under UV light, making up for the disadvantage that it is easy for a viologen radical to fade in air, and providing a new strategy for designing and synthesizing viologen materials with strong fluorescence.

A radical approach to radicals.

Nitroxide radicals are characterized by a long-lived spin-unpaired electronic ground state and are strongly sensitive to their chemical surroundings. Combined with electron paramagnetic resonance spectroscopy, these electronic features have led to the widespread application of nitroxide derivatives as spin labels for use in studying protein structure and dynamics. Site-directed spin labelling requires the incorporation of nitroxides into the protein structure, leading to a new protein-ligand molecular model. However, in protein crystallographic refinement nitroxides are highly unusual molecules with an atypical chemical composition. Because macromolecular crystallography is almost entirely agnostic to chemical radicals, their structural information is generally less accurate or even erroneous. In this work, proteins that contain an example of a radical compound (Chemical Component Dictionary ID MTN) from the nitroxide family were re-refined by defining its ideal structural parameters based on quantum-chemical calculations. The refinement results show that this procedure improves the MTN ligand geometries, while at the same time retaining higher agreement with experimental data.

Structural and Energetic Properties of Amino Acids and Peptides Benchmarked by Accurate Theoretical and Experimental Data.

Structural, energetic, and spectroscopic data derived in this work aim at the setup of an "experimentally validated" database for amino acids and polypeptides conformers. First, the "cheap" composite scheme (ChS, CCSD(T)/(CBS+CV)MP2) is tested for evaluation of conformational energies of all eight stable conformers of glycine, by comparing to the more accurate CCSD(T)/CBS+CV computations (Phys. Chem. Chem. Phys. 2013, 15, 10094-10111 and J Mol. Model. 2020, 26, 129). The recently proposed jun-ChS (J. Chem. Theory and Comput. 2020, 16, 988-1006), employing the jun-cc-pVnZ basis set family for CCSD(T) computations and CBS extrapolation, yields conformational energies accurate to 0.2 kJ·mol-1, at reduced computational cost with respect to aug-ChS employing aug-cc-pVnZ basis sets. The jun-ChS composite scheme is further applied to derive conformational energies for three dipeptide analogues Ac-Gly-NH2, Ac-Ala-NH2, and Gly-Gly. Finally, dipeptide conformational energies and semiexperimental equilibrium rotational constants along with the CCSD(T)/(CBS+CV)MP2 structural parameters (J. Phys. Chem. Lett. 2014, 5, 534-540) stand as the reference for benchmarking of selected density functional methodologies. The double-hybrid functionals B2-PLYP-D3(BJ) and DSD-PBEP86, perform best for structural and energetic characterization of all dipeptide analogues. From hybrid functionals CAM-B3LYP-D3(BJ) and ωB97X-D3(BJ) represent promising methods applicable for larger peptide-based systems for which computations with double-hybrid functionals are not feasible.

Structural and Vibrational Properties of Amino Acids from Composite Schemes and Double-Hybrid DFT: Hydrogen Bonding in Serine as a Test Case.

The structures, relative stabilities, and vibrational wavenumbers of the two most stable conformers of serine, stabilized by the O-H···N, O-H···O═C and N-H···O-H intramolecular hydrogen bonds, have been evaluated by means of state-of-the-art composite schemes based on coupled-cluster (CC) theory. The so-called "cheap" composite approach (CCSD(T)/(CBS+CV)MP2) allowed determination of accurate equilibrium structures and harmonic vibrational wavenumbers, also pointing out significant corrections beyond the CCSD(T)/cc-pVTZ level. These accurate results stand as a reference for benchmarking selected hybrid and double-hybrid, dispersion-corrected DFT functionals. B2PLYP-D3 and DSDPBEP86 in conjunction with a triple-ζ basis set have been confirmed as effective methodologies for structural and spectroscopic studies of medium-sized flexible biomolecules, also showing intramolecular hydrogen bonding. These best performing double-hybrid functionals have been employed to simulate IR spectra by means of vibrational perturbation theory, also considering hybrid CC/DFT schemes. The best overall agreement with experiment, with mean absolute error of 8 cm-1, has been obtained by combining CCSD(T)/(CBS+CV)MP2 harmonic wavenumbers with B2PLYP-D3/maug-cc-pVTZ anharmonic corrections. Finally, a composite scheme entirely based on CCSD(T) calculations (CCSD(T)/CBS+CV) has been employed for energetics, further confirming that serine II is the most stable conformer, also when zero-point vibrational energy corrections are included.

Real-space quantum-based refinement for cryo-EM: Q|R#3.

Electron cryo-microscopy (cryo-EM) is rapidly becoming a major competitor to X-ray crystallography, especially for large structures that are difficult or impossible to crystallize. While recent spectacular technological improvements have led to significantly higher resolution three-dimensional reconstructions, the average quality of cryo-EM maps is still at the low-resolution end of the range compared with crystallography. A long-standing challenge for atomic model refinement has been the production of stereochemically meaningful models for this resolution regime. Here, it is demonstrated that including accurate model geometry restraints derived from ab initio quantum-chemical calculations (HF-D3/6-31G) can improve the refinement of an example structure (chain A of PDB entry 3j63). The robustness of the procedure is tested for additional structures with up to 7000 atoms (PDB entry 3a5x and chain C of PDB entry 5fn5) using the less expensive semi-empirical (GFN1-xTB) model. The necessary algorithms enabling real-space quantum refinement have been implemented in the latest version of qr.refine and are described here.

Toward accurate prediction of amino acid derivatives structure and energetics from DFT: glycine conformers and their interconversions.

This work provides the accurate reference data for structural and energetic properties relevant for computational studies of amino acids and polypeptides. Glycine due to its small size allows for detailed theoretical explorations of its whole conformational space. The reference energies are computed at the CCSD(T)/CBS+CV level on the best estimated geometries of all local minima and the transition states. For the minima, we complete the set of reference structures reported in Phys. Chem. Chem. Phys., 2013, 15, 10094, by determining the CCSD(T)/(CBS+CV)MP2 quality geometries for the two highest energy conformers VIIp/tcc and VIIIn/gtc. These data stand as the reference to asses reliability and accuracy of less expensive computational models, with particular focus on dispersion-corrected hybrid and double-hybrid DFT approaches, considering several basis sets of double- and triple-ζ quality. Based on results for minima the B2LYP-D3(BJ)/aug-cc-pVTZ level is set as the reference for transition states geometries. Considering accuracy of both single point energies and structural parameters B2PLYP(-D3BJ) and DSDPBEP86 in conjunction with aug-cc-pVTZ basis set can be recommended for reference studies of amino-acids and small poly-peptides, while B3LYP(-D3(BJ)) shown to be the most robust from considered hybrid DFT approaches. Graphical abstract TOC: Benchmarking DFT for amino acid analogues: highly accurate characterization of energies and structures from full glycine conformational space.

Including crystallographic symmetry in quantum-based refinement: Q|R#2.

Three-dimensional structure models refined using low-resolution data from crystallographic or electron cryo-microscopy experiments can benefit from high-quality restraints derived from quantum-chemical methods. However, nonperiodic atom-centered quantum-chemistry codes do not inherently account for nearest-neighbor interactions of crystallographic symmetry-related copies in a satisfactory way. Here, these nearest-neighbor effects have been included in the model by expanding to a super-cell and then truncating the super-cell to only include residues from neighboring cells that are interacting with the asymmetric unit. In this way, the fragmentation approach can adequately and efficiently include nearest-neighbor effects. It has previously been shown that a moderately sized X-ray structure can be treated using quantum methods if a fragmentation approach is applied. In this study, a target protein (PDB entry 4gif) was partitioned into a number of large fragments. The use of large fragments (typically hundreds of atoms) is tractable when a GPU-based package such as TeraChem is employed or cheaper (semi-empirical) methods are used. The QM calculations were run at the HF-D3/6-31G level. The models refined using a recently developed semi-empirical method (GFN2-xTB) were compared and contrasted. To validate the refinement procedure for a non-P1 structure, a standard set of crystallographic metrics were used. The robustness of the implementation is shown by refining 13 additional protein models across multiple space groups and a summary of the refinement metrics is presented.

Conformational Equilibria and Molecular Structures of Model Sulfur-Sulfur Bridge Systems: Diisopropyl Disulfide.

The conformations and molecular structures of diisopropyl disulfide have been studied by high-resolution microwave spectroscopy and quantum chemical calculations. Three conformers, G'GG', G'GT, and GGG', have been observed in the jet expansion. The global minimum, G'GG', adopts a configuration with the G' orientation of H-C-S-S and S-S-C-H and the G orientation of C-S-S-C showing the C2 symmetry. The rotational spectra of monosubstituted 13C and 34S isotopologues have also been recorded for G'GG', leading to an accurate structural determination of this conformer. Two additional 34S isotopologues have also been measured for G'GT. The relative energies of three observed conformers calculated at the MP2/6-311++(d,p) level of theory are within 2 kJ mol-1, while the relative intensity measurements suggested their population ratio to be NG'GG'/NG'GT/NGGG' ≈ 5:3:2.

Absorption-emission symmetry breaking and the different origins of vibrational structures of the 1Qy and 1Qx electronic transitions of pheophytin a.

The vibrational structure of the optical absorption and fluorescence spectra of the two lowest-energy singlet electronic states (Qy and Qx) of pheophytin a were carefully studied by combining low-resolution and high-resolution spectroscopy with quantum chemical analysis and spectral modeling. Large asymmetry was revealed between the vibrational structures of the Qy absorption and fluorescence spectra, integrally characterized by the total Huang-Rhys factor and reorganization energy in absorption of Svib A = 0.43 ± 0.06, λA = 395 cm-1 and in emission of Svib E = 0.35 ± 0.06, λE = 317 cm-1. Time-dependent density-functional theory using the CAM-B3LYP, ωB97XD, and MN15 functionals could predict and interpret this asymmetry, with the exception of one vibrational mode per model, which was badly misrepresented in predicted absorption spectra; for CAM-B3LYP and ωB97XD, this mode was a Kekulé-type mode depicting aromaticity. Other computational methods were also considered but performed very poorly. The Qx absorption spectrum is broad and could not be interpreted in terms of a single set of Huang-Rhys factors depicting Franck-Condon allowed absorption, with Herzberg-Teller contributions to the intensity being critical. For it, CAM-B3LYP calculations predict that Svib A (for modes >100 cm-1) = 0.87 and λA = 780 cm-1, with effective x and y polarized Herzberg-Teller reorganization energies of 460 cm-1 and 210 cm-1, respectively, delivering 15% y-polarized intensity. However, no method was found to quantitatively determine the observed y-polarized contribution, with contributions of up to 50% being feasible.

The role of dispersion and anharmonic corrections in conformational analysis of flexible molecules: the allyl group rotamerization of matrix isolated safrole.

Conformational changes of the monomeric safrole (5-(2-propenyl)-1,3-benzodioxole) isolated in low temperature xenon matrices were induced thermally or using narrow-band UV radiation. The rotation of the allyl group taking place in the studied matrices was followed by FTIR spectroscopy. Safrole represents a challenging example of a flexible molecule highlighting the importance of dispersion interactions and anharmonic effects in the structural, spectroscopic and energetic analysis. Structures of the safrole conformers, their energetics and infrared spectra have been calculated using various computational methods ranging from density functional theory (DFT) to coupled cluster (CC). The best theoretical results were obtained by integrating CCSD(T) energies including complete basis set extrapolation and core-valence corrections with B2PLYP-D3 equilibrium structures and hybrid B2PLYP-D3/B3LYP-D3 anharmonic computations for IR spectra and thermodynamics.

Theoretical studies of atmospheric molecular complexes interacting with NIR to UV light.

The interaction of weakly bonded complexes of atmospheric constituents with the electromagnetic spectrum available in Earth's atmosphere can induce direct excitation to electronic excited states as well as the excitation of higher vibrational states (overtones) of the electronic ground state. A better understanding of these phenomena requires improved theoretical support by including the anharmonic and vibro-electronic effects on both the band positions and transition intensities. In this work, generalized second-order vibrational perturbation and time-independent Franck-Condon and Herzberg-Teller computations are exploited together with a density functional theory (DFT)/coupled cluster (CC) scheme and its extension to the excited electronic states. Structural and spectroscopic properties are calculated for isolated formaldehyde and its complexes with H2O, CO, SO2 and H2O2, focusing on how small molecules may affect the interactions with NIR to UV irradiation.

Hydrogen detachment driven by a repulsive 1πσ* state - an electron localization function study of 3-amino-1,2,4-triazole.

Electron localization function analysis reveals the details of a charge induced hydrogen detachment mechanism of 3-amino-1,2,4-triazole, identified recently to be responsible for phototautomerization of the molecule. In this process vertical excitation to the 1πσ* state is followed by the barrier-less migration of a H atom along the N-H bond toward the conical intersection with the S0 ground state. The most striking feature revealed for the 1πσ* state is partial ejection of σ* electrons outside the molecule, even beyond the NH group, at the Franck-Condon point. Further gradual spatial localization of the electron around the proton moving along the N-H stretching coordinate gives a plausible explanation for the repulsive character of the 1πσ* potential energy surface with the proton wading through the region of space where some negative charge is accumulated ('a virtual acceptor'), dragging some electron density. This mechanism resembles the one postulated for the hydrogen transfer from a donor molecule (D-H) to an acceptor one (A) in a class of vertically excited molecules with a preexisting inter- or intramolecular D-HA motif, even though the acceptor molecule is absent. The present analysis demonstrates also that the bond evolution and changes in the electron density along the excited state reaction path can be effectively studied with the use of an electron localization function.