Bringing Machine-Learning Enhanced Quantum Chemistry and Microwave Spectroscopy to Conformational Landscape Exploration: the Paradigmatic Case of 4-Fluoro-Threonine

Engineering Amino Acid and Peptide Supramolecular Architectures through Fluorination

Fluorinated non-natural amino acids are attracting considerable research interest, especially in the biomedical field and in materials science, thanks to their ability to self-assemble into peculiar supramolecular structures. The conformational changes induced by the presence of fluorine atoms obviously affect their functions, as well as the biological activity of the deriving peptides and proteins. Here, we will briefly describe the main effects of fluorination on the aggregation behavior of such building blocks, focusing in particular on their improved tendency to form fibrils, and gels therefrom. Our aim is to underline the promising potential of fluorination as a tool to affect the self-assembly features of amino acids, both when used alone and when inserted into polypeptide sequences. The ability of fluorine to influence physical, chemical, and structural properties of these substrates offers the possibility to engineer bioinspired materials with specific and tunable functions.

Unbiased Comparison between Theoretical and Experimental Molecular Structures and Properties: Toward an Accurate Reduced-Cost Evaluation of Vibrational Contributions

The tremendous development of hardware and software is constantly increasing the role of quantum chemical (QC) computations in the assignment and interpretation of experimental results. However, an unbiased comparison between theory and experiment requires the proper account of vibrational averaging effects. In particular, high-resolution spectra in the gas phase are now available for molecules containing up to about 50 atoms, which are too large for a brute-force approach with the available QC methods of sufficient accuracy. In the present paper, we introduce hybrid approaches, which allow the accurate evaluation of vibrational averaging effects for molecules of this size beyond the harmonic approximation, with special attention being devoted to rotational constants. After the validation of new tools for relatively small molecules, the β-estradiol hormone and a prototypical molecular motor have been considered to witness the feasibility of accurate computations for large molecules.

Quantum chemistry meets high-resolution spectroscopy for characterizing the molecular bricks of life in the gas-phase

Computation of accurate geometrical structures and spectroscopic properties of large flexible molecules in the gas-phase is tackled at an affordable cost using a general exploration/exploitation strategy. The most distinctive feature of the approach is the careful selection of different quantum chemical models for energies, geometries and vibrational frequencies with the aim of maximizing the accuracy of the overall description while retaining a reasonable cost for all the steps. In particular, a composite wave-function method is used for energies, whereas a double-hybrid functional (with the addition of core-valence correlation) is employed for geometries and harmonic frequencies and a cheaper hybrid functional for anharmonic contributions. A thorough benchmark based on a wide range of prototypical molecular bricks of life shows that the proposed strategy is close to the accuracy of state-of-the-art composite wave-function methods, and is applicable to much larger systems. A freely available web-utility post-processes the geometries optimized by standard electronic structure codes paving the way toward the accurate yet not prohibitively expensive study of medium- to large-sized molecules by experimentally-oriented researchers.

Gas Phase Conformation of Trisaccharides and Core Pentasaccharide: A Three-Step Tree-Based Sampling and Quantum Mechanical Computational Approach

As an important component of N-linked glycoproteins, the core pentasaccharide is highly crucial to the potential application prospect of glycoprotein. However, the gas phase conformation study is a challenging one due to the size and complexity of the molecule, together with the necessity to rely on quantum chemistry modeling for relevant energetics and structures. In this paper, the structures of the trisaccharides and core pentasaccharides in N-linked glycans in the gas phase were constructed by a three-step tree-based (TSTB) sampling. Since single point energies of all the conformers are calculated at the temperature of zero, it is necessary to evaluate the stability at a high temperature. We calculate the Gibbs free energies using the standard thermochemistry model (T = 298.15 K). For trimannose, the energetic ordering at 298.15 K can be strongly changed compared to 0 K. Moreover, two structures of trimannose with high energies at 0 K are considered to provide a much better match of IR vibration signatures with the low Gibbs free energies. On this basis, the core pentasaccharide was constructed in three ways. The building configurations of core pentasaccharide were optimized to obtain reasonable low-energy stable conformers. Fortunately, the lowest-energy structure of core pentasaccharide is eventually the minimum at 0 K and 298.15 K. Furthermore, spectrum analysis of core pentasaccharide was carried out. Although poorly resolved, its contour from the experiment was in qualitative correspondence with the computed IR spectrum associated with its minimum free energy structure. A large number of strongly and weakly hydrogen-bonded hydroxyl and acetylamino groups contribute to a highly congested set of overlapping bands. Compared with traditional conformation generators, the TSTB sampling is employed to efficiently and comprehensively obtain preferred conformers of larger saccharides with lower energy.

Accurate structures and spectroscopic parameters of α,α-dialkylated α-amino acids in the gas-phase: a joint venture of DFT and wave-function composite methods

Accurate computations of structural, conformational and spectroscopic properties in the gas phase have been performed for two α,α-dialkylated α-amino acids, namely aminoisobutyric acid and cyclopropylglycine. Thanks to the integration of modern double hybrid functionals and wave-function methods, several low-energy structures of the title molecules could be analyzed employing standard computer resources. The computed features of all the most stable conformers of the target amino acids closely match the corresponding spectroscopic parameters issued from microwave spectroscopic studies in the gas-phase. Together with their intrinsic interest, the accuracy of the results obtained with reasonable computer times paves the way for accurate investigations of other flexible bricks of life.

DFT Meets Wave-Function Composite Methods for Characterizing Cytosine Tautomers in the Gas Phase

A general strategy for the accurate computation of structural and spectroscopic properties of biomolecule building blocks in the gas phase has been further improved and validated with a special reference to tautomeric equilibria. The main improvements concern the use of the cc-pVTZ-F12 basis set in both DFT and CCSD(T)-F12 computations, the inclusion of core-valence correlation in geometry optimizations by double hybrid functionals, and the use of the cc-pVQZ-F12 basis set for complete basis set extrapolation at the MP2-F12 level. The resulting model chemistry is applied to the challenging problem of cytosine tautomers in the gas phase. The results are in remarkable agreement with experiment concerning both rotational and vibrational spectroscopic parameters and permit their unbiased interpretation in terms of structural and thermochemical features. Together with the intrinsic interest of the studied molecule, the accuracy of the results obtained at reasonable cost without any empirical parameter suggests that the proposed composite method can be profitably employed for accurate investigations of other molecular bricks of life.

Accurate Structures and Spectroscopic Parameters of Phenylalanine and Tyrosine in the Gas Phase: A Joint Venture of DFT and Composite Wave-Function Methods

A general strategy for the accurate computation of conformational and spectroscopic properties of flexible molecules in the gas phase is applied to two representative proteinogenic amino acids with aromatic side chains, namely, phenylalanine and tyrosine. The main features of all the most stable conformers predicted by this computational strategy closely match those of the species detected in microwave and infrared experiments. Together with their intrinsic interest, the accuracy of the results obtained with reasonable computer times paves the route for accurate investigations of other flexible bricks of life.

Computational Protocol for the Identification of Candidates for Radioastronomical Detection and Its Application to the C3H3NO Family of Isomers

The C3H3NO family of isomers is relevant in astrochemistry, even though its members are still elusive in the interstellar medium. To identify the best candidate for astronomical detection within this family, we developed a new computational protocol based on the minimum-energy principle. This approach aims to identify the most stable isomer of the family and consists of three steps. The first step is an extensive investigation that characterizes the vast number of compounds having the C3H3NO chemical formula, employing density functional theory for this purpose. The second step is an energy refinement, which is used to select isomers and relies on coupled cluster theory. The last step is a structural improvement with a final energy refinement that provides improved energies and a large set of accurate spectroscopic parameters for all isomers lying within 30 kJ mol−1 above the most stable one. According to this protocol, vinylisocyanate is the most stable isomer, followed by oxazole, which is about 5 kJ mol−1 higher in energy. The other stable species are pyruvonitrile, cyanoacetaldehyde, and cyanovinylalcohol. For all of these species, new computed rotational and vibrational spectroscopic data are reported, which complement those already available in the literature or fill current gaps.