Accuracy of Rotational Parameters Predicted by High-Level Quantum-Chemical Calculations: Case Study of Sulfur-Containing Molecules of Astrochemical Interest

Repurposing raltegravir for reducing inflammation and treating cancer: a bioinformatics analysis

Inflammation is a defensive mechanism that safeguards the human body against detrimental stimuli. Within this intricate process, ADAM17, a zinc-dependent metalloprotease, emerges as an indispensable element, fostering the activation of diverse inflammatory and growth factors within the organism. Nonetheless, ADAM17 malfunctions can augment the rate of growth, inflammatory factors, and subsequent damage. Thus, in this study, we examined and repurposed drugs to suppress the activity of ADAM17. To this end, we employed bioinformatics techniques such as molecular docking, molecular dynamics, and pharmacokinetic studies. Five FDA-approved drugs including Raltegravir, Conivaptan, Paclitaxel, Saquinavir, and Venetoclax with the ability to impede the activity of the ADAM17 metalloenzyme were identified. Moreover, these drugs did not include strong zinc-binding functional groups when verified by the ACE functional group finder. However, further in silico analysis has indicated that Raltegravir demonstrates a commendable interaction with the active site amino acids and exhibits the most favorable pharmacokinetic properties compared to others. Considering the results of bioinformatics tools, it can be concluded that Raltegravir as an antiviral drug could be repurposed to prevent severe inflammatory response and tumorigenesis resulting from ADAM17 malfunction.

The Semiexperimental Approach at Work: Equilibrium Structure of Radical Species

The so-called semiexperimental (SE) approach is a powerful technique for obtaining highly accurate equilibrium structures for isolated systems. This Featured Article describes its extension to open-shell species, thus providing the first systematic investigation on radical equilibrium geometries to be used for benchmarking purposes. The small yet significant database obtained demonstrates that there is no reduction in accuracy when moving from closed-shell species to radicals. We also provide an extension of the applicability of the SE approach to medium-/large-sized radicals by exploiting the so-called "Lego-brick" approach, which is based on the assumption that a molecular system can be seen as formed by smaller fragments for which the SE equilibrium structure is available. In this Featured Article we show that this model can be successfully applied also to open-shell species.

Isomerism of CH 2 SO : Accurate structural, energetic, and spectroscopic characterization

A recent work [Ye et al. Mon. Not. R. Astron. Soc. 2023, 525, 1158] on the gas-phase formation of t-HC(O)SH, already detected in the interstellar medium, pointed out that the trans form of HC(S)OH is a potential candidate for astronomical observations. Prompted by these results, theCH 2 SO family of isomers has been investigated from an energetic point of view using a double-hybrid density functional in combination with a partially augmented triple-zeta basis set. This preliminary study showed that the most stable species of the family are the cis and trans forms of HC(O)SH and HC(S)OH. For their structural and spectroscopic characterization, a composite scheme based on coupled cluster (CC) calculations that incorporates up to the quadruple excitations and accounts for the extrapolation to the complete basis set limit and core correlation effects has been employed. This approach opens to the prediction of rotational constants with an accuracy of 0.1%. A hybrid scheme, based on harmonic frequencies computed using the CC singles, doubles and a perturbative treatment of triples method (CCSD(T)) in conjunction with a quadruple-zeta basis set, allowed us to obtain fundamental vibrational frequencies with a mean absolute error of about 1%.

Rotational spectra and semi-experimental structures of furonitrile and its water cluster

Rotational spectroscopy represents an invaluable tool for several applications: from the identification of new molecules in interstellar objects to the characterization of van der Waals complexes, but also for the determination of very accurate molecular structures and for conformational analyses. In this work, we used high-resolution rotational spectroscopic techniques in combination with high-level quantum-chemical calculations to address all these aspects for two isomers of cyanofuran, namely 2-furonitrile and 3-furonitrile. In particular, we have recorded and analyzed the rotational spectra of both of them from 6 to 320 GHz; rotational transitions belonging to several singly-substituted isotopologues have been identified as well. The rotational constants derived in this way have been used in conjunction with computed rotation-vibration interaction constants in order to derive a semi-experimental equilibrium structure for both isomers. Moreover, we observed the rotational spectra of four different intermolecular adducts formed by furonitrile and water, whose identification has been supported by a conformational analysis and a theoretical spectroscopic characterization. A semi-experimental determination of the intermolecular parameters has been achieved for all of them and the results have been compared with those obtained for the analogous system formed by benzonitrile and water.

Hunting for interstellar molecules: rotational spectra of reactive species

Interstellar molecules are often highly reactive species, which are unstable under terrestrial conditions, such as radicals, ions and unsaturated carbon chains. Their detection in space is usually based on the astronomical observation of their rotational fingerprints. However, laboratory investigations have to face the issue of efficiently producing these molecules and preserving them during rotational spectroscopy measurements. A general approach for producing and investigating unstable/reactive species is presented by means of selected case-study molecules. The overall strategy starts from quantum-chemical calculations that aim at obtaining accurate predictions of the missing spectroscopic information required to guide spectral analysis and assignment. Rotational spectra of these species are then recorded by exploiting the approach mentioned above, and their subsequent analysis leads to accurate spectroscopic parameters. These are then used for setting up accurate line catalogs for astronomical searches.

Design, synthesis, computational study and cytotoxic evaluation of some new quinazoline derivatives containing pyrimidine moiety

Quinazoline derivatives, as an important category of heterocyclic compounds, have received much attention for the design and development of new drugs due to their various pharmacological properties. Besides, there is a great deal of evidence showing pyrimidine analogs as anticancer agents. Thus, in the present study, for the design of new target compounds with cytotoxic activity, we focused on various quinazolinone and pyrimidine hybrids. A new series of quinazoline-pyrimidine hybrid derivatives (6a-6n) have been designed and synthesized as novel antiproliferative agents. All the synthesized compounds characterized based on their IR, NMR and Mass spectroscopic data. Antiproliferative activities of the new compounds were evaluated against three human cancer cell lines (MCF-7, A549, SW-480). The compounds were found to have appropriate potential with IC50 values ranging from 2.3 ± 5.91 to 176.5 ± 0.7 μM against the tested cell lines. Compound 6n exerted the highest antiproliferative activity with IC50 values of 5.9 ± 1.69 μM, 2.3 ± 5.91 μM and 5.65 ± 2.33 μM against A549, SW-480 and MCF-7 respectively. The results indicated that 6n could induce apoptosis in A549 cell line in a dose dependent manner and arrest in the S phase of cell cycle. Docking studies were also done to investigate the detailed binding pattern of the synthesized compounds against EGFR. Furthermore, molecular dynamic simulation and binding free energy calculation have been done to rescore initial docking pose of the synthesized compounds using ensemble-based MMGB/PBSA free energy method. According to the results, free energy calculation confirmed biological activity of compounds and also, Arg 817 and Lys 721 residues had the pivotal role in the high potency of 6n. Finally, the drug likeness and in silico ADME study were also predicted.

Describing nuclear quantum effects in vibrational properties using molecular dynamics with Wigner sampling

In this work we discuss the generally applicable Wigner sampling and introduce a new, simplified Wigner sampling method, for computationally effective modeling of molecular properties containing nuclear quantum effects and vibrational anharmonicity. For various molecular systems test calculations of (a) vibrationally averaged rotational constants, (b) vibrational IR spectra and (c) photoelectron spectra have been performed. The performance of Wigner sampling has been assessed by comparing with experimental data and with results of other theoretical models, including harmonic and VPT2 approximations. The developed simplified Wigner sampling method shows advantages in application to large and flexible molecules.

Insights into the molecular structure and infrared spectrum of the prebiotic species aminoacetonitrile

Aminoacetonitrile is an interstellar molecule with a prominent prebiotic role, already detected in the chemically-rich molecular cloud Sagittarius B2(N) and postulated to be present in the atmosphere of the largest Saturn's moon, Titan. To further support its observation in such remote environments and laboratory experiments aimed at improving our understanding of interstellar chemistry, we report a thorough spectroscopic and structural characterization of aminoacetonitrile. Equilibrium geometry, fundamental bands as well as spectroscopic and molecular parameters have been accurately computed by exploiting a composite scheme rooted in the coupled-cluster theory that accounts for the extrapolation to the complete basis set limit and core-correlation effects. In addition, a semi-experimental approach that combines ground-state rotational constants for different isotopic species and calculated vibrational corrections has been employed for the structure determination. From the experimental side, we report the analysis of the three strongest fundamental bands of aminoacetonitrile observed between 500 and 1000 cm^-1 in high-resolution infrared spectra. More generally, all computed band positions are in excellent agreement with the present and previous experiments. The only exception is the ν₁₅ band, for which we provide a revision of the experimental assignment, now in good agreement with theory.

Spectroscopic and Computational Characterization of 2-Aza-1,3-butadiene, a Molecule of Astrochemical Significance

Being N-substituted unsaturated species, azabutadienes are molecules of potential relevance in astrochemistry, ranging from the interstellar medium to Titan’s atmosphere. 2-Azabutadiene and butadiene share a similar conjugated π system, thus allowing investigation of the effects of heteroatom substitution. More interestingly, 2-azabutadiene can be used to proxy the abundance of interstellar butadiene. To enable future astronomical searches, the rotational spectrum of 2-azabutadiene has been investigated up to 330 GHz. The experimental work has been supported and guided by accurate computational characterization of the molecular structure, energetics, and spectroscopic properties of the two possible forms, trans and gauche. The trans species, more stable by about 7 kJ/mol than gauche-2-azabutadiene, has been experimentally observed, and its rotational and centrifugal distortion constants have been obtained with remarkable accuracy, while theoretical estimates of the spectroscopic parameters are reported for gauche-2-azabutadiene.

Ab Initio Study of Fine and Hyperfine Interactions in Triplet POH

Phosphorous-containing molecules have a great relevance in prebiotic chemistry in view of the fact that phosphorous is a fundamental constituent of biomolecules, such as RNA, DNA, and ATP. Its biogenic importance has led astrochemists to investigate the possibility that P-bearing species could have formed in the interstellar medium (ISM) and subsequently been delivered to early Earth by rocky bodies. However, only two P-bearing molecules have been detected so far in the ISM, with the chemistry of interstellar phosphorous remaining poorly understood. Here, in order to shed further light on P-carriers in space, we report a theoretical spectroscopic characterisation of the rotational spectrum of POH in its 3A″ ground electronic state. State-of-the-art coupled-cluster schemes have been employed to derive rotational constants, centrifugal distortion terms, and most of the fine and hyperfine interaction parameters, while the electron spin-spin dipolar coupling has been investigated using the multi-configuration self-consistent-field method. The computed spectroscopic parameters have been used to simulate the appearance of triplet POH rotational and ro-vibrational spectra in different conditions, from cold to warm environments, either in gas-phase experiments or in molecular clouds. Finally, we point out that the predicted hyperfine structures represent a key pattern for the recognition of POH in laboratory and interstellar spectra.

Exploiting the “Lego brick” approach to predict accurate molecular structures of PAHs and PANHs

Polycyclic aromatic hydrocarbons (PAHs) and polycyclic aromatic nitrogen heterocycles (PANHs) are important and ubiquitous species in space. However, their accurate structural and spectroscopic characterization is often missing. To fill this...

Extending the Applicability of the Semi-experimental Approach by Means of "Template Molecule" and "Linear Regression" Models on Top of DFT Computations

The accurate determination of equilibrium structures for isolated molecules plays a central role in the evaluation and interpretation of stereoelectronic, thermodynamic, and spectroscopic properties. For small semi-rigid systems, state-of-the-art quantum-chemical computations can rival the most sophisticated experimental results. For larger molecules, cheaper yet accurate approaches need to be defined. The double-hybrid rev-DSD-PBEP86 functional already delivers remarkable results that can be further improved by means of a "Lego brick" model. This is based on the idea that a molecular system can be seen as formed by different fragments (the "Lego bricks"), whose accurate semi-experimental (SE) equilibrium geometries are available. The template molecule (TM) approach can be used to account for the modifications occurring when going from the isolated fragment to the molecular system under investigation, with the linear regression (LR) model employed to correct the linkage between the different fragments. The resulting TM-SE_LR approach has been tested with respect to available SE equilibrium structures and rotational constants. Indeed, the latter parameters straightforwardly depend on the equilibrium geometry of the system under consideration. The main outcome of our study is the reliability, robustness, and accuracy of this novel approach. The molecular systems considered for benchmarking the TM-SE_LR scheme are those formally issued from addition/elimination reactions of nucleophilic unsaturated radicals (e.g., CN, C2H, and phenyl) to alkenes, imines, and aldehydes, whose rotational spectra have been investigated, but accurate structural determinations are not yet available.

Highly-accurate quartic force fields for the prediction of anharmonic rotational constants and fundamental vibrational frequencies

The CcCR quartic force field (QFF) methodology is capable of computing B₀ and C₀ rotational constants to within 35 MHz (0.14%) of experiment for triatomic and larger molecules with at least two heavy atoms. Additionally, the same constants for molecules with four or more atoms agree to within 20 MHz (0.12%) of experiment for the current test set. This work also supports previous claims that the same QFF methodology can produce fundamental vibrational frequencies with a deviation less than 5.7 cm^-1 from experiment. Consequently, this approach of augmenting complete basis set extrapolated energies with treatments of core electron correlation and scalar relativity produces some of the most accurate rovibrational spectroscopic data available.

In search of phosphorus in astronomical environments: The reaction between the CP radical (X2Σ+) and methanimine

Phosphorus is of particular interest in astrochemistry because it is a biogenic element together with hydrogen, carbon, nitrogen, oxygen, and sulfur. However, the chemical evolution of such element in the interstellar medium (ISM) is still far from an accurate characterization, with the chemistry of P-bearing molecules being poorly understood. To provide a contribution in this direction, we have carried out an accurate investigation of the potential energy surface for the reaction between the CP radical and methanimine (CH₂NH), two species already detected in the ISM. In analogy to similar systems, i.e., CH₂NH + X, with X = OH, CN, and CCH, this reaction can occur-from an energetic point of view-under the harsh conditions of the ISM. Furthermore, since the major products of the aforementioned reaction, namely, E- and Z-2-phosphanylidyneethan-1-imine (HN=CHCP) and N-(phosphaneylidynemethyl)methanimine (H₂C=NCP), have not been spectroscopically characterized yet, some effort has been made for filling this gap by means of accurate computational approaches.

Looking for the bricks of the life in the interstellar medium: The fascinating world of astrochemistry

The discovery in the interstellar medium of molecules showing a certain degree of complexity, and in particular those with a prebiotic character, has attracted great interest. A complex chemistry takes place in space, but the processes that lead to the production of molecular species are a matter of intense discussion, the knowledge still being at a rather primitive stage. Debate on the origins of interstellar molecules has been further stimulated by the identification of biomolecular building blocks, such as nucleobases and amino acids, in meteorites and comets. Since many of the molecules found in space play a role in the chemistry of life, the issue of their molecular genesis and evolution might be related to the profound question of the origin of life itself. Understanding the underlying chemical processes, including the production, reactions and destruction of compounds, requires the concomitant study of spectroscopy, gas-phase reactivity, and heterogeneous processes on dust-grains. The aim of this contribution is to provide a general view of a complex and multifaceted challenge, while focusing on the role played by molecular spectroscopy and quantum-chemical computations. In particular, the derivation of the molecular spectroscopic features and the investigation of gas-phase formation routes of prebiotic species in the interstellar medium are addressed from a computational point of view.

Structure, Stability, and Spectroscopic Properties of Small Acetonitrile Cation Clusters

Ionization and fragmentation pathways induced by ionizing agents are key to understanding the formation of complex molecules in astrophysical environments. Acetonitrile (CH₃CN), the simplest organic nitrile, is an important molecule present in the interstellar medium. In this work, DFT and MP2 calculations were performed in order to obtain the low energy structures of the most relevant cations formed from electron-stimulated ion desorption of CH₃CN ices. Selected reaction pathways and spectroscopic properties were also calculated. Our results indicate that the most stable acetonitrile cation structure is CH₂CNH⁺ and that hydrogenation can occur successively without isomerization steps until its complete saturation. Moreover, the stability of distinct cluster families formed from the interaction of acetonitrile with small fragments, such as CH_n⁺, C₂H_n⁺, and CH_nCNH⁺, is discussed in terms of their respective binding energies. Some of these molecular clusters are stabilized by hydrogen bonds, leading to species whose infrared features are characterized by a strong redshift of the N-H stretching mode. Finally, the rotational spectra of CH₃CN and protonated acetonitrile, CH₃CNH⁺, were simulated using distinct computational protocols based on DFT, MP2, and CCSD(T) considering centrifugal distortion, vibrational-rotational coupling, and vibrational anharmonicity corrections. By adopting an empirical scaling procedure for calculating spectroscopic parameters, we were able to estimate the rotational frequencies of CH₃CNH⁺ with an expected average error below 1 MHz for J values up to 10.

The Role of State-of-the-Art Quantum-Chemical Calculations in Astrochemistry: Formation Route and Spectroscopy of Ethanimine as a Paradigmatic Case

The gas-phase formation and spectroscopic characteristics of ethanimine have been re-investigated as a paradigmatic case illustrating the accuracy of state-of-the-art quantum-chemical (QC) methodologies in the field of astrochemistry. According to our computations, the reaction between the amidogen, NH, and ethyl, C2H5, radicals is very fast, close to the gas-kinetics limit. Although the main reaction channel under conditions typical of the interstellar medium leads to methanimine and the methyl radical, the predicted amount of the two E,Z stereoisomers of ethanimine is around 10%. State-of-the-art QC and kinetic models lead to a [E-CH3CHNH]/[Z-CH3CHNH] ratio of ca. 1.4, slightly higher than the previous computations, but still far from the value determined from astronomical observations (ca. 3). An accurate computational characterization of the molecular structure, energetics, and spectroscopic properties of the E and Z isomers of ethanimine combined with millimeter-wave measurements up to 300 GHz, allows for predicting the rotational spectrum of both isomers up to 500 GHz, thus opening the way toward new astronomical observations.

The challenge of non-covalent interactions: theory meets experiment for reconciling accuracy and interpretation

In the past decade, many gas-phase spectroscopic investigations have focused on the understanding of the nature of weak interactions in model systems. Despite the fact that non-covalent interactions play a key role in several biological and technological processes, their characterization and interpretation are still far from being satisfactory. In this connection, integrated experimental and computational investigations can play an invaluable role. Indeed, a number of different issues relevant to unraveling the properties of bulk or solvated systems can be addressed from experimental investigations on molecular complexes. Focusing on the interaction of biological model systems with solvent molecules (e.g., water), since the hydration of the biomolecules controls their structure and mechanism of action, the study of the molecular properties of hydrated systems containing a limited number of water molecules (microsolvation) is the basis for understanding the solvation process and how structure and reactivity vary from gas phase to solution. Although hydrogen bonding is probably the most widespread interaction in nature, other emerging classes, such as halogen, chalcogen and pnicogen interactions, have attracted much attention because of the role they play in different fields. Their understanding requires, first of all, the characterization of the directionality, strength, and nature of such interactions as well as a comprehensive analysis of their competition with other non-covalent bonds. In this review, it is shown how state-of-the-art quantum-chemical computations combined with rotational spectroscopy allow for fully characterizing intermolecular interactions taking place in molecular complexes from both structural and energetic points of view. The transition from bi-molecular complex to microsolvation and then to condensed phase is shortly addressed.

The challenging playground of astrochemistry: an integrated rotational spectroscopy - quantum chemistry strategy

While it is now well demonstrated that the interstellar medium (ISM) is characterized by a diverse and complex chemistry, a significant number of features in radioastronomical spectra are still unassigned and call for new laboratory efforts, which are increasingly based on integrated experimental and computational strategies. In parallel, the identification of an increasing number of molecules containing more than five atoms and at least one carbon atom (the so-called "interstellar" complex organic molecules), which can play a relevant role in the chemistry of life, raises the additional issue of how these species can be produced in the typical harsh conditions of the ISM. On these grounds, this perspective aims to present an integrated rotational spectroscopy - quantum chemistry approach for supporting radioastronomical observations and a computational strategy for contributing to the elucidation of chemical reactivity in the interstellar space.

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 C₂ symmetry. The rotational spectra of monosubstituted ¹³C and ³⁴S isotopologues have also been recorded for G'GG', leading to an accurate structural determination of this conformer. Two additional ³⁴S 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 N_G'GG'/N_G'GT/N_GGG' ≈ 5:3:2.

Interaction and Photodissociation of Electronic Excited States of HS2 in the Ultraviolet Region: A Theoretical Contribution

HS₂ molecules play an important role in the photochemical processes in combustion, atmosphere, and interstellar medium, yet our knowledge about the electronic excited states in the UV region is limited. In this study, we perform high-level ab initio calculations on electronic states of HS₂ using the internally contracted multireference configuration interaction method including the Davidson correction (icMRCI + Q) method. The vertical transition energies, oscillator strengths, electron configurations, and transitions of thirteen electronic states of HS₂ with energy up to 8 eV are calculated at the icMRCI + Q/aug-cc-pv(5 + d)Z level. Based on the calculated potential energy curves, we investigate the interaction and photodissociation mechanism of the electronic states, which should shed some light on the decomposition processes of the gas-phase HS₂ molecules in the ultraviolet region.