CPL of Mellein and Related Natural Compounds: Analysis of the ESIPT Phenomenon.

(R)-(-)-Mellein, (3R,4R)-4-hydroxymellein and (3R,4S)-4-hydroxymellein obtained from fungi, i.e. from Diplodia globulosa, were investigated as a class of natural products presenting ESIPT (excited state intramolecular proton transfer) phenomenon, through fluorescence and CPL (circularly polarized luminescence). The study was preceded by the assessment of the absolute configuration through ECD and VCD (electronic and vibrational circular dichroism) spectroscopies in addition to NMR spectra. It is found that ESIPT takes place in these systems very rapidly, and no dual fluorescence has been observed. The experimental study is backed up by TD-DFT calculations of ECD and CPL spectra, plus MD dynamics to follow proton transfer in the excited state and careful analysis of the puckering dynamics of the lactone ring. Deprotonated forms of the three compounds were also investigated by the same chiroptical experimental and theoretical methods, showing how one can find in natural compounds not only biological activity but also biologically compatible sensing probes.

Mid-IR and CH stretching vibrational circular dichroism spectroscopy to distinguish various sources of chirality: The case of quinophaneoxazoline based ruthenium(II) complexes.

Five diastereomers of ruthenium(II) complexes based on quinolinophaneoxazoline ligands were investigated by vibrational circular dichroism (VCD) in the mid-IR and CH stretching regions. Diastereomers differ in three sources of chirality: the planar chirality of the quinolinophane moiety, the central chirality of an asymmetric carbon atom of the oxazoline ring, and the chirality of the ruthenium atom. VCD, allied to DFT calculations, has been found to be effective in disentangling the various forms of chirality. In particular, a VCD band is identified in the CH stretching region directly connected to the chirality of the metal. The analysis of the calculated VCD spectra is carried out by partitioning the complexes into fragments. The anharmonic analysis is also performed with a recently proposed reduced-dimensionality approach: such treatment is particularly important when examining spectroscopic regions highly perturbed by resonances, like the CH stretching region.

Circularly polarized luminescence of natural products lycorine and narciclasine: role of excited-state intramolecular proton-transfer and test of pH sensitivity.

Circularly polarized luminescence (CPL) is increasingly gaining interest not only for its applicative potentialities but also for providing an understanding of the excited state properties of chiral molecules. However, applications of CPL are mainly in the field of materials science: special organic molecules and polymers, metal (lanthanide) complexes, and organic dyes are actively and intensely studied. So far natural compounds have not been investigated much. We fill the gap here by measuring circular dichroism (CD) and CPL of lycorine and narciclasine, the most abundant known alkaloid and isocarbostyril from Amaryllidaceae, which exhibit a large spectrum of biological activities and are promising anticancer compounds. Dual fluorescence detection in narciclasine led us to unveil an occurring excited-state intramolecular proton transfer (ESIPT) process, this mechanism well accounts for the Stokes shift and CPL spectra observed in narciclasine. The same molecule is interesting also as a pH chiroptical switch. Both in absorption and emission, lycorine and narciclasine are also studied computationally via density functional theory (DFT) calculations further shedding light on their properties.

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.

Insights into the Structures of Bilirubin and Biliverdin from Vibrational and Electronic Circular Dichroism: History and Perspectives.

In this work we review research activities on a few of the most relevant structural aspects of bilirubin (BR) and biliverdin (BV). Special attention is paid to the exocyclic C=C bonds being in mostly Z rather than E configurations, and to the overall conformation being essentially different for BR and BV due to the presence or absence of the double C=C bond at C-10. In both cases, racemic mixtures of each compound of either M or P configuration are present in achiral solutions; however, imbalance between the two configurations may be easily achieved. In particular, results based on chiroptical spectroscopies, both electronic and vibrational circular dichroism (ECD and VCD) methods, are presented for chirally derivatized BR and BV molecules. Finally, we review deracemization experiments monitored with ECD data from our lab for BR in the presence of serum albumin and anesthetic compounds.

Anharmonic Aspects in Vibrational Circular Dichroism Spectra from 900 to 9000 cm-1 for Methyloxirane and Methylthiirane.

Vibrational circular dichroism (VCD) spectra and the corresponding IR spectra of the chiral isomers of methyloxirane and of methylthiirane have been reinvestigated, both experimentally and theoretically, with particular attention to accounting for anharmonic corrections, as calculated by the GVPT2 approach. De novo recorded VCD spectra in the near IR (NIR) range regarding CH-stretching overtone transitions, together with the corresponding NIR absorption spectra, were also considered and accounted for, both with the GVPT2 and with the local mode approaches. Comparison of the two methods has permitted us to better describe the nature of active "anharmonic" modes in the two molecules and the role of mechanical and electrical anharmonicity in determining the intensities of VCD and IR/NIR data. Finally, two nonstandard IR/NIR regions have been investigated: the first one about ≈2000 cm-1, involving mostly two-quanta bending mode transitions, the second one between 7000 and 7500 cm-1 involving three-quanta transitions containing CH-stretching overtones and HCC/HCH bending modes.

Unbiased disentanglement of conformational baths with the help of microwave spectroscopy, quantum chemistry, and artificial intelligence: The puzzling case of homocysteine.

An integrated experimental-computational strategy for the accurate characterization of the conformational landscape of flexible biomolecule building blocks is proposed. This is based on the combination of rotational spectroscopy with quantum-chemical computations guided by artificial intelligence tools. The first step of the strategy is the conformer search and relative stability evaluation performed by means of an evolutionary algorithm. In this step, last generation semiempirical methods are exploited together with hybrid and double-hybrid density functionals. Next, the barriers ruling the interconversion between the low-lying conformers are evaluated in order to unravel the possible fast relaxation paths. The relative stabilities and spectroscopic parameters of the "surviving" conformers are then refined using state-of-the-art composite schemes. The reliability of the computational procedure is further improved by the inclusion of vibrational and thermal effects. The final step of the strategy is the comparison between experiment and theory without any ad hoc adjustment, which allows an unbiased assignment of the spectroscopic features in terms of different conformers and their spectroscopic parameters. The proposed approach has been tested and validated for homocysteine, a highly flexible non-proteinogenic α-amino acid. The synergism of the integrated strategy allowed for the characterization of five conformers stabilized by bifurcated N-H2⋯O=C hydrogen bonds, together with an additional conformer involving a more conventional HN⋯H-O hydrogen bond. The stability order estimated from the experimental intensities as well as the number and type of conformers observed in the gas phase are in full agreement with the theoretical predictions. Analogously, a good match has been found for the spectroscopic parameters.

Alternative Strategy to Obtain Artificial Imine Reductase by Exploiting Vancomycin/D-Ala-D-Ala Interactions with an Iridium Metal Complex.

Based on the supramolecular interaction between vancomycin (Van), an antibiotic glycopeptide, and D-Ala-D-Ala (DADA) dipeptides, a novel class of artificial metalloenzymes was synthesized and characterized. The presence of an iridium(III) ligand at the N-terminus of DADA allowed the use of the metalloenzyme as a catalyst in the asymmetric transfer hydrogenation of cyclic imines. In particular, the type of link between DADA and the metal-chelating moiety was found to be fundamental for inducing asymmetry in the reaction outcome, as highlighted by both computational studies and catalytic results. Using the [IrCp*(m-I)Cl]Cl ⊂ Van complex in 0.1 M CH3COONa buffer at pH 5, a significant 70% (S) e.e. was obtained in the reduction of quinaldine B.

Exploring the Maze of Cycloserine Conformers in the Gas Phase Guided by Microwave Spectroscopy and Quantum Chemistry.

Cycloserine has in common with isoxazolidines the saturated five-membered ring, which is an important scaffold for drug design, exhibiting diverse biological activities. The most remarkable feature of these compounds is the presence of the N-O bond framed in a cyclic moiety. The lack of an accurate characterization of this structural feature in an isolated system calls for a state-of-the-art theoretical-experimental study. A quantum-chemical investigation of cycloserine unveiled the presence of 11 local energy minima, with only two of them being separated by significant barriers. This picture has been experimentally confirmed: two species have been unequivocally detected in the gas phase by means of laser ablation microwave spectroscopy, also disentangling the complicated hyperfine structure originating from the presence of two nitrogen atoms. A thorough characterization of cycloserine and isoxazolidine, benchmarked by the semiexperimental investigation of hydroxylamine, provided the first accurate determination of their structures and pointed out that the rev-DSD-PBEP86 functional is competitive with respect to explicitly correlated coupled-cluster computations. This outcome paves the way toward accurate studies of large flexible molecules.

A Computational Journey across Nitroxide Radicals: From Structure to Spectroscopic Properties and Beyond.

Nitroxide radicals are characterized by a long-lived open-shell electronic ground state and are strongly sensitive to the chemical environment, thus representing ideal spin probes and spin labels for paramagnetic biomolecules and materials. However, the interpretation of spectroscopic parameters in structural and dynamic terms requires the aid of accurate quantum chemical computations. In this paper we validate a computational model rooted into double-hybrid functionals and second order vibrational perturbation theory. Then, we provide reference quantum chemical results for the structures, vibrational frequencies and other spectroscopic features of a large panel of nitroxides of current biological and/or technological interest.

Chemical bonding in cuprous complexes with simple nitriles: octet rule and resonance concepts versus quantitative charge-redistribution analysis.

Chemical bonding in a set of six cuprous complexes with simple nitriles (CN-, HNC, HCN, CH3NC, and CH3CN) is investigated by means of a recently devised analysis scheme framed in density-functional theory and quantitatively singling out concurrent charge flows such as σ donation and π backdonation. The results of our analysis are comparatively assessed against qualitative models for charge redistribution based on the popular concepts of octet rule and resonance structures, and the relative importance of different charge-flow channels relating to σ donation, π back-donation, polarization, and hyperconjugation is discussed on a quantitative basis.

Toward Fully Unsupervised Anharmonic Computations Complementing Experiment for Robust and Reliable Assignment and Interpretation of IR and VCD Spectra from Mid-IR to NIR: The Case of 2,3-Butanediol and trans-1,2-Cyclohexanediol.

The infrared (IR) and vibrational circular dichroism (VCD) spectra of 2,3-butanediol and trans-1,2-cyclohexanediol from 900 to 7500 cm-1 (including mid-IR, fundamental CH and OH stretchings, and near-infrared regions) have been investigated by a combined experimental and computational strategy. The computational approach is rooted in density functional theory (DFT) computations of harmonic and leading anharmonic mechanical, electrical, and magnetic contributions, followed by a generalized second-order perturbative (GVPT2) evaluation of frequencies and intensities for all the above regions without introducing any ad hoc scaling factor. After proper characterization of large-amplitude motions, all resonances plaguing frequencies and intensities are taken into proper account. Comparison of experimental and simulated spectra allows unbiased assignment and interpretation of the most interesting features. The reliability of the GVPT2 approach for OH stretching fundamentals and overtones is confirmed by the remarkable agreement with a local mode model purposely tailored for the latter two regions. Together with the specific interest of the studied molecules, our results confirm that an unbiased assignment and interpretation of vibrational spectra for flexible medium-size molecules can be achieved by means of a nearly unsupervised reliable, robust, and user-friendly DFT/GVPT2 model.

Unsupervised search of low-lying conformers with spectroscopic accuracy: A two-step algorithm rooted into the island model evolutionary algorithm.

The fruitful interplay of high-resolution spectroscopy and quantum chemistry has a long history, especially in the field of small, semi-rigid molecules. However, in recent years, the targets of spectroscopic studies are shifting toward flexible molecules, characterized by a large number of closely spaced energy minima, all contributing to the overall spectrum. Here, artificial intelligence comes into play since it is at the basis of powerful unsupervised techniques for the exploration of soft degrees of freedom. Integration of such algorithms with a two-stage QM/QM' (Quantum Mechanical) exploration/refinement strategy driven by a user-friendly graphical interface is the topic of the present paper. We will address in particular: (i) the performances of different semi-empirical methods for the exploration step and (ii) the comparison between stochastic and meta-heuristic algorithms in achieving a cheap yet complete exploration of the conformational space for medium sized chromophores. As test cases, we choose three amino acids of increasing complexity, whose full conformer enumeration has been reached only very recently. Next, we show that systems in condensed phases can be treated at the same level and with the same efficiency when employing a polarizable continuum description of the solvent. Finally, the challenging issue represented by the vibrational circular dichroism spectra of some rhodium complexes with flexible ligands has been addressed, showing that our fully unsupervised approach leads to remarkable agreement with the experiment.

Discovering the Elusive Global Minimum in a Ternary Chiral Cluster: Rotational Spectra of Propylene Oxide Trimer.

The chirality controlled conformational landscape of the trimer of propylene oxide (PO), a prototypical chiral molecule, was investigated using rotational spectroscopy and a range of theoretical tools for conformational searches and for evaluating vibrational contributions to effective structures. Two sets of homochiral (PO)3 rotational transitions were assigned and the associated conformers identified with theoretical support. One set of heterochiral (PO)3 transitions was assigned, but no structures generated by one of the latest, advanced conformational search codes could account for them. With the aid of a Python program, the carbon atom backbone and then the heterochiral (PO)3 structure were generated using 13 C isotopic data. Excellent agreement between theoretical and experimental rotational constants and relative dipole moment components of all three conformers was achieved, especially after applying vibrational corrections to the rotational constants.

Vibrational circular dichroism under the quantum magnifying glass: from the electronic flow to the spectroscopic observable.

We present a comprehensive methodology for the analysis and interpretation of vibrational circular dichroism spectra supported by novel graphical representations. The tools rely on the vibrational transition current density (VTCD) associated with a molecular vibration, whose visualization allows exploration of the physical origin of the electronic contribution to the electric and magnetic vibrational dipole transition moments. Different ways of visualizing VTCD from 2D maps to 3D representations are reported and applied to molecular systems of growing complexity. An extension of the VTCD analysis to fully anharmonic spectra within the second-order vibrational perturbation theory (VPT2) is discussed. The analysis is applied to different types of chiral systems: the doubly deuterated oxirane (2S,3S)-oxirane-d2, taken as a reference to validate our implementation; 1,3-difluoroallene, a simple rigid system to explore the application of VTCD to anharmonic VCD spectra. Finally, the analysis of VTCD has been used to better understand the origin of the signal enhancement in peptides linked to the ferrocene groups. The merits and shortcomings of the methods are discussed, and some perspectives for future developments are offered.

Two-level stochastic search of low-energy conformers for molecular spectroscopy: implementation and validation of MM and QM models.

The search for stationary points in the molecular potential energy surfaces (PES) is a problem of increasing relevance in different fields of molecular sciences especially for large, flexible systems characterized by several large-amplitude internal motions leading to shallow minima with comparable energies and separated by small barriers. After structural biology and medicinal chemistry, also high-resolution molecular spectroscopy, which is the focus of our research activity, is nowadays shifting its attention to this kind of molecular systems. In such circumstances, accurate geometrical structures and relative stabilities of all these minima are a mandatory prerequisite for the vis-à-vis comparison between computed and experimental spectra. This task raises, in turn, the problem of the best compromise between accuracy and feasibility. In our opinion, a promising route is offered by a two-level stochastic search in which a relatively inexpensive MM or QM method is used in the initial search, followed by single point energy evaluation at a higher QM level of a relatively large number of low-energy structures in order to select a final short-list of candidates, whose geometries are fully optimized at the higher QM level. Finally, the relative stabilities and properties of the final short-list of energy minima can be computed by a state-of-the-art QM approach. This strategy defines a general two-level search/three-level evaluation approach, which can be finely tuned in terms of the accuracy of the sought results. Setup of the procedure, interface with a general purpose electronic structure code and validation of the most effective low-level methods for some representative molecular systems (three already well characterized and one new) ended up with a general, robust and user-friendly tool, which can be easily used and extended also by non-specialists to aid experimental spectroscopic studies.

Diving into chemical bonding: An immersive analysis of the electron charge rearrangement through virtual reality.

An integrated environment for the analysis of chemical bonding based on immersive virtual reality is presented. Using a multiscreen stereoscopic projection system, researchers are cast into the world of atoms and molecules, where they can visualize at a human scale the electron charge rearrangement (computed via state-of-the-art quantum-chemical methods) occurring on bond formation throughout the molecular region. Thanks to specifically designed features, such a virtual laboratory couples the immediacy of an immersive experience with a powerful, recently developed method yielding quantitative, spatially detailed pictures of the several charge flows involved in the formation of a chemical bond. By means of two case studies on organometallic complexes, we show how familiar concepts in coordination chemistry, such as donation and back-donation charge flows, can be effectively identified and quantified to predict experimental observables. © 2018 Wiley Periodicals, Inc.

Towards the SMART workflow system for computational spectroscopy.

Is it possible to convert highly specialized research in the field of computational spectroscopy into robust and user-friendly aids to experiments and industrial applications? What kind of tools should be created to increase the interactions between researchers with different backgrounds and push towards new frontiers in computational chemistry? The outstanding advances in computational spectroscopy and the wide availability of computational and analytical tools are paving the route toward the study of problems that were previously difficult or impossible to solve and enable the imagination of even more ambitious targets for fundamental and applied research. The combination of new computer- and data-centric technologies is transforming data analysis from an uncommon and retrospective practice into a proactive process of strategic decision and action. This paper starts from these premises and proposes a perspective for a new cyberinfrastructure aimed at integrating developments in theory, algorithms and software with new tools for workflow management, data mining and visualization. We make a case for this approach by means of a few examples that deal with unmanageable types of data in molecular modelling and results obtained using different unsupervised learning algorithms.

Computational simulation of vibrationally resolved spectra for spin-forbidden transitions.

In this computational study, we illustrate a method for computing phosphorescence and circularly polarized phosphorescence spectra of molecular systems, which takes into account vibronic effects including both Franck-Condon and Herzberg-Teller contributions. The singlet and triplet states involved in the phosphorescent emission are described within the harmonic approximation, and the method fully takes mode-mixing effects into account when evaluating Franck-Condon integrals. Spin-orbit couplings, which are responsible for these otherwise forbidden phenomena, are accounted for by means of a relativistic two-component time-dependent density functional theory method. The model is applied to two types of chiral systems: camphorquinone, a rigid organic system that allows for an extensive benchmark, and some members of a class of iridium complexes. The merits and shortcomings of the methods are discussed, and some perspectives for future developments are offered.

On the relation between carbonyl stretching frequencies and the donor power of chelating diphosphines in nickel dicarbonyl complexes.

The relation between spectroscopic observables and the detailed metal-ligand bonding features in chelation complexes is addressed using both experimental and state-of-the-art theoretical and computational methods. We synthesized and characterized a set of six nickel dicarbonyl complexes of general formula [Ni(CO)2(PP)], where PP is an atropoisomeric chelating diphosphine ligand. The analysis of the obtained experimental data and the basicity and oxidative potentials of the free ligands suggests a close relation between the donor ability of the chelating ligand and the carbonyl stretching frequencies observed in the complexes. We then use theory to unravel the detailed mechanisms of chelation-bond formation in terms of partial charge flows between the molecular orbitals of the fragments. By extending the promising, recently published natural orbitals for chemical valence/charge displacement (NOCV/CD) analysis scheme we provide a thorough, quantitative description of the several charge fluxes following the metal-ligand bond formation and demonstrate that the carbonyl stretching frequencies in the considered complexes selectively respond to the σ-donation charge flow from the phosphorus lone pairs of the ligands, with the frequency shift being in quantitative correlation with the extent of the ligand-to-metal charge transfer.