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.

Accurate geometries for "Mountain pass" regions of the Ramachandran plot using quantum chemical calculations.

Unusual local arrangements of protein in Ramachandran space are not well represented by standard geometry tools used in either protein structure refinement using simple harmonic geometry restraints or in protein simulations using molecular mechanics force fields. In contrast, quantum chemical computations using small poly-peptide molecular models can predict accurate geometries for any well-defined backbone Ramachandran orientation. For conformations along transition regions-ϕ from -60 to 60°-a very good agreement with representative high-resolution experimental X-ray (≤1.5 Å) protein structures is obtained for both backbone C-1 -N-Cα angle and the nonbonded O-1 …C distance, while "standard geometry" leads to the "clashing" of O…C atoms and Amber FF99SB predicts distances too large by about 0.15 Å. These results confirm that quantum chemistry computations add valuable support for detailed analysis of local structural arrangements in proteins, providing improved or missing data for less understood high-energy or unusual regions.

Structural and Vibrational Properties of Iodopentafluorobenzene: A Combined Raman and Infrared Spectral and Theoretical Study.

A combined study of the vibrational spectroscopy of iodopentafluorobenzene by new Raman and Fourier transform infrared (FTIR) spectroscopies, over the spectral range 300 to 3200 cm-1 (Raman) and 50 to 3400 cm-1 (FTIR), with a state-of-the-art theoretical investigation is reported. This has enabled reliable identification of numerous fundamental, overtone, and combination band transitions in unprecedented detail. The theoretical analysis, beyond the double-harmonic approximation, is based on generalized second-order vibrational perturbation theory (GVPT2) with a hybrid coupled cluster/density functional theory (CC/DFT) approach. Anharmonic contributions to structural parameters, rotational constants, vibrational frequencies, and spectral intensities are incorporated. The procedures, of general applicability, enable rigorous comparison of theoretical methods with experimental results in vibrational spectroscopy.

The ionic states of difluoromethane: A reappraisal of the low energy photoelectron spectrum including ab initio configuration interaction computations.

A new synchrotron-based study of the photoelectron spectrum (PES) of difluoromethane is interpreted by an ab initio analysis of the ionic states, which includes Franck-Condon (FC) factors. Double differentiation of the spectrum leads to significant spectral sharpening; the vibrational structure observed is now measured with greater accuracy than in previous studies. Several electronic structure methods are used, including equation of motion coupled cluster calculations with single and double excitations (EOM-CCSD), its ionization potential variant EOM-IP-CCSD, 4th order Møller-Plesset perturbation theory (MP4SDQ) configuration interaction (CI), and complete active space self-consistent-field (CASSCF) methods. The adiabatic ionization energies (AIEs) confirm the assignments as band I, one state 12B1 (12.671 eV); band II, three states, 12B2 (14.259) < 12A1 (15.030) < 12A2 (15.478 eV); and band III, three states, 22B2 (18.055) < 22A1 (18.257) < 22B1 (18.808 eV). The three ionizations in each of the bands II and III lead to selective line broadening of the PES structure, which is attributed to vibronic overlap. The apparent lack of a vibrational structure attributable to both the 12A1 and 22A1 states in the PES arises from line broadening with the preceding states 12B2 and 22B2, respectively. Although these 2A1 states clearly overlap with their adjacent higher IE, some vibrational structure is observed on the higher IE. The effects of vibronic coupling are evident since the observed structure does not fit closely with the calculated Born-Oppenheimer FC profiles. Correlation of the lowest group of four AIEs in the PES of other members of the CH2X2 group, where X = F, Cl, Br, and I, clearly indicate these effects are more general.

Accurate spectroscopic characterization of the HOC(O)O radical: A route toward its experimental identification.

A set of accurate spectroscopic parameters for the detection of the atmospherically important HOC(O)O radical has been obtained by means of state-of-the-art ab initio computations. These include advanced coupled cluster treatments, involving both standard and explicitly correlated approaches, to correctly account for basis set incompleteness and core-valence effects. Geometric parameters for the X̃2A' and Ã2A'' states and, for the ground state only, vibrationally corrected rotational constants including quartic and sextic centrifugal distortion terms are reported. The infrared spectrum of the X̃2A' state has been simulated in the 4000-400 cm-1 wavenumber interval with an approach based on second order vibrational perturbation theory that allows accounting for anharmonic effects in both energies and intensities. Finally, the vibronic spectrum for the à ← X̃ transition has been calculated at three different temperatures in the 9000-3000 cm-1 energy range with a time-independent technique based on the Franck-Condon approximation.