Better performance of Hartree-Fock over DFT: a quantum mechanical investigation on pyridinium benzimidazolate types of zwitterions in the light of localization/delocalization issues.

CONTEXT: With the advent of fast computing facilities, combined with rapid emerges of many new and intricate quantum mechanical functionals, computations with pure Hartree-Fock (HF) theory are now-a-days regarded as trivial or obsolete, or even considered as not reliable by many researchers. Consequently, current trends in computational chemistry show extensive use of post-HF theories for smaller molecular systems and various DFT methods for organic and inorganic chemistry related problems (larger molecules/systems). In this contribution, I have tried to show that sometimes, HF might be more suitable over DFT methodologies in addressing structure-property correlations. Molecules studied here were previously synthesized by Boyd in 1966 and important experimental data were produced by Alcalde and co-workers in 1987. Comparison of computed and experimental results clearly shows that HF method was more effective in reproducing the experimental data compared to especially the DFT methodologies. Reliability of HF method was further assured from the very similar results shown by the CCSD, CASSCF, CISD and QCISD methods. Current study also indicates that the localization issue associated with HF proved to be advantageous over delocalization issue of DFT based methodologies, in correctly describing the structure-property correlation for zwitterion systems. METHODS: All computations were performed with Gaussian 09. A wide-range of quantum mechanical methodologies, HF, B3LYP, CAM-B3LYP, BMK, B3PW91, TPSSh, LC-ωPBE, M06-2X, M06-HF, ωB97xD, MP2, CASSCF, CCSD, QCISD, CISD and semi-empirical methods like, Huckel, CNDO, AM1, PM3MM and PM6, were used for investigations.

Ortho-para interconversion of nuclear states of H2O through replica transition state: prospect of quantum entanglement at homodromic Bjerrum defect site.

CONTEXT: From a nuclear spin prospective, water exists as para and ortho nuclear spin isomers (isotopomers). Spin interconversions in isolated molecules of water are forbidden, but many recent reports have shown them to happen in bulk, through dynamic proton exchanges happening between interconnected networks of a large array of water molecules. In this contribution, a possible explanation for an unexpected slow or delayed interconversion of ortho-para water in ice observed in an earlier reported experiment is provided. Using the results of quantum mechanical investigations, we have discussed the roles played by Bjerrum defects in the dynamic proton exchanges and ortho-para spin state interconversions. We guess that at the sites of the Bjerrum defects, there are possibilities of quantum entanglements of states, through pairwise interactions. Based on the perfectly correlated exchange happening via a replica transition state, we speculate that it can have significant influences on ortho-para interconversions of water. We also conjecture that the overall ortho-para interconversion is not a continuous process, rather can be imagined to be happening serendipitously, but within the boundary of the rules of quantum mechanics. METHODS: All computations were performed with Gaussian 09 program. B3LYP/6-31++G(d,p) methodology was used to compute all the stationary points. Further energy corrections were computed using CCSD(T)/aug-cc-pVTZ methodology. Intrinsic reaction coordinate (IRC) path computations were carried out for the transition states.

Planar in Brooker's mode and twisted in Reichardt's mode: defying the steric forces in biphenyl types of zwitterionic systems through metameric resonance stabilizations.

To be planar or to be twisted at the bridge junctions in biphenyls or biaryl types of molecular systems depends on two conflicting forces: (1) steric repulsions (destabilizations) and (2) conjugation assisted electron delocalizations (resonance stabilizations). This work reports an unfamiliar kind of behaviour shown by metamers of a zwitterionic biphenyl type of system, where the Reichardt's metamer was found to be in an usual twisted conformation (delicate balance of conflicting forces), but the Brooker's metamer was found to be in a fully planar conformation. Interestingly, at the ωB97xD/aug-cc-pVDZ level, energetically (ΔE) the planar Brooker's metamer was found to be 16.7 kcal mol-1 lower (22.9 kcal mol-1 lower in the CASSCF method) in energy (more stable) than the isoelectronic twisted Reichardt's metamer, and also thermodynamic ΔG values were found to be close to ΔE values for various methods (for example, 15.6 kcal mol-1 in the above case using the ωB97xD method). When the steric repulsions are in their full potentials at the ring junction site, attainment of a conformational planarity by any biaryl type of system has not been reported previously. Without reducing the steric constraints or even without inducing any attractive forces, determining what other factors were responsible for defying the steric forces is the main focus of this investigation. Using the results of quantum mechanical computations of NBO, rotational barriers, and other saddle points (metastable conformations in singlet and triplet surfaces) in the potential energy surfaces, the dominant contribution of the resonance stabilized quinonoid form to the ground state was delineated as the possible reason for this unusual behaviour.

Formation of Formamide from HCN + H2O: A Computational Study on the Roles of a Second H2O as a Catalyst, as a Spectator, and as a Reactant.

Formamide (NH2CHO), being the smallest and fundamental building block of life (with a peptide linkage), has recently been able to attract much interests, in the field of astrochemistry, astrophysics, and astrobiology. In this work, using quantum mechanical computations, reactions between HCN and H2O, leading to the formation of formamide, have been analyzed. For the first time, an alternative and competing reaction channel, which proceeds via a geminal diol intermediate, for the formation of formamide, has been proposed. In this alternative channel, an extra water molecule (second H2O) was found to be acting as a reactant, in the second step of the reaction path. Effects of second H2O molecule in the reaction paths, providing catalytic assistance to the reaction or behaving like a spectator (concept is introduced for the first time for this reaction), have also been analyzed. Usefulness of spectator behavior is highlighted for the reactions happening on the rigid water-ice surfaces, where the water-ice may not be getting involved for any catalytic assistance. In light of catalytic assistances provided by the second H2O, prominent effects in reducing the barrier heights drastically (even for the second step of the reaction, the barrier height was found to be below the reactants), through a hydrogen relay transport mechanism, were observed. In addition to the mechanism studies, interstellar feasibilities of all the reaction channels and their significances are discussed in detail.

Hemiaminal route for the formation of interstellar glycine: a computational study.

Calculations related to two simple two-step paths (path-I: [Formula: see text] path-II: [Formula: see text]) for the formation of glycine have been discussed. Calculations show that at interstellar conditions these two paths are feasible only in hot cores, not in the cold interstellar clouds (cold core formation is possible only if CH2 = NH, H2O (excess) and CO of path-II, react in a concerted manner). For the laboratory synthesis of glycine, the possibility suggested is via path-I and the reaction being carried out as controlled temperature one-pot synthesis. This study can also be extended to other α-amino acids and possibly enantiomeric excess can be expected. We think this work will not only be able to enrich our future understanding about the formation of amino acids in interstellar medium but also be able to suggest alternative paths for laboratory synthesis of amino acids using either Strecker's or Miller's ingredients. Graphical abstract Using computational calculations, two different reaction paths which go through a hemiaminal (α-hydroxyamine) intermediate have been proposed. It has been proposed that the reaction [Formula: see text] is a thermodynamically favorable reaction path in the laboratory conditions, if carried out as a controlled temperature one-pot synthesis. On the hand, it has been argued that the reaction[Formula: see text] is a feasible reaction path in the interstellar conditions, if it proceeds not via the hemiaminal route, rather in a concerted reaction path.

Reaction between NH3 (X̌1A1) and CO (X1Σ+): A Computational Insight into the Reaction Mechanism of Formamide (H2N-CHO) Formation.

Formamide (NH2CHO) is the smallest molecular unit that contains the basic peptide linkage and thus has recently attracted a great amount of interest in the field of astrochemistry. In this work using computational calculations, we have analyzed the three possible reaction paths for the reaction between CO and NH3 to form formamide in both neutral-neutral and cation-neutral reaction surfaces. All of these three paths strongly favor the path of 1,2-hydrogen migration, which was discounted by previous studies in view of the constraints from steric factor. We have also analyzed the significant role played by prereaction complexes in these three reaction paths. We have proposed that for the neutral-neutral reaction path, formation of formamide in the low temperature interstellar clouds was hypothesized to proceed via hydrogen tunneling assisted by a tunneling ready like state as prereaction complex. On the other hand, for the two cation-neutral reactions, any tunneling cannot facilitate formation of formamide in the interstellar clouds. Rather in one case as all the stationary points are below the reactants, it can facilitate the reaction, whereas in the second case the reaction is only possible if it can get some catalytic assistance.

Efficient discrimination of natural stereoisomers of chicoric acid, an HIV-1 integrase inhibitor.

Plants from the Asteraceae family are known to contain a wide spectrum of phytochemicals with various nutraceutical properties. One important phytochemical, chicoric acid (CA), is reported to exist in plants, such as Sonchus oleraceus and Bidens pilosa, as stereoisomers. These CA molecules occur either as the naturally abundant RR-chicoric acid (RR-CA), or the less abundant RS-chicoric acid (RS-CA), also known as meso-chicoric acid. To date, little is known about the biological activity of RS-CA, but there is evidence of its anti-human immunodeficiency virus (HIV) properties. In this study, a reliable analytical method was developed to distinguish between the two stereoisomers detected in S. oleraceus and B. pilosa. For structure identification and characterization of CA molecules, liquid chromatography-mass spectrometry (LC-MS) was used in combination with ultraviolet radiation (UV)-induced geometrical isomerization, molecular dynamics (MD) simulations, and density functional theory (DFT) models. Optimized structures from DFT calculations were used for docking studies against the HIV-1 integrase enzyme. Different retention times on the reverse phase chromatograms revealed that the plants produce two different CA stereoisomers: S. oleraceus produced the RR-CA isomer, while B. pilosa produced the RS-CA isomer. DFT results demonstrated the RR-CA molecule was more stable than RS-CA due to the stabilizing force of intra-molecular hydrogen bonding. Differences in the HIV-1 integrase enzyme binding modes were observed, with the RR-CA being a more potent inhibitor than the RS-CA molecule. The results highlight the significance of plant metabolite structural complexity from both chemical and biological perspectives. Furthermore, the study demonstrates that induced-formation of geometrical isomers, in combination with the predictive ability of DFT models and the resolving power of the LC-MS, can be exploited to distinguish structurally closely related compounds, such as stereoisomers.

Possible interstellar formation of glycine through a concerted mechanism: a computational study on the reaction of CH2[double bond, length as m-dash]NH, CO2 and H2.

Glycine being the simplest amino acid and also having significant astrobiological implications, has meant that intensive investigations have been carried out in the past, starting from its detection in the interstellar medium (ISM) to analysis of meteorites and cometary samples and laboratory synthesis, as well as computational studies on the possible reaction paths. In this present work quantum chemical calculations have been performed to investigate the possible interstellar formation of glycine via two different paths; (1) in a two-step process via a dihydroxy carbene intermediate and (2) through a one-step concerted mechanism, starting from reactants like CH2[double bond, length as m-dash]NH, CO, CO2, H2O and H2. For the two reactions representing the carbene route, it was observed that the formation of dihydroxy carbene from either CO + H2O or CO2 + H2 is highly endothermic with large barrier heights, whereas the subsequent step of interaction of this carbene with CH2[double bond, length as m-dash]NH to give glycine is exothermic and the barrier is below the reactants. Based on this observation it is suggested that the formation of glycine via the carbene route is a least favourable or even unfavourable path. On the other hand, the two reactions CH2[double bond, length as m-dash]NH + CO + H2O and CH2[double bond, length as m-dash]NH + CO2 + H2 representing the concerted paths were found to be favourable in leading to the formation of glycine. After an extensive study on the first concerted reaction in our previous work (Phys. Chem. Chem. Phys., 2016, 18, 375-381), in this work a detailed investigation has been carried out for the second concerted reaction, CH2[double bond, length as m-dash]NH + CO2 + H2, which can possibly lead to the interstellar formation of glycine. It was observed that this reaction proceeds through a large barrier and at the same time the transition state shows prominent hydrogen dynamics, indicating a tunnelling possibility for this reaction. Based on these observations the possible formation of glycine via this reaction in hot-cores and in cold interstellar clouds has been proposed. The cold-core possibility of this reaction is argued on the basis of the phenomenon of tunnelling assisted by a van der Waals' complex.

Expanding the applicability of electrostatic potentials to the realm of transition states.

Central to any reaction mechanism study, and sometimes a challenging job, is tracing a transition state in a reaction path. For the first time, electrostatic potentials (ESP) of the reactants were used as guiding tactics to predict whether there is a possibility of any transition state in a reaction surface. The main motive behind this strategy is to see whether the directionality nature of the transition state has something to do with the anisotropic natures of the ESP with their embedded directionalities. Strategically, some atmospherically important, but simple, reactions have been chosen for this study, which heretofore were believed to be barrierless. By carefully analysing the ESP maps of the reactants, regions of possible interactions were located. Using the bilinear interpolation of the 2D grids of the ESP surfaces, search co-ordinates were fine-tuned for a local gradient based approach for the search of a transition state. Out of the three reactions studied in this work, we were able to successfully locate transition states, for the first time, in two cases and the third one still proved to be barrierless. This gives a clear indication that though ESP maps can qualitatively predict the possibility of a transition state; it is not always true that there should definitely be a transition state, as some of the reaction surfaces may genuinely be barrierless. But, nevertheless this strategy definitely has credential to be tested for many more reactions, either new or already established, and may be applied to create the initial search co-ordinates for any well-established transition state search method. Moreover, we have observed that the analysis of the ESP maps of the reactants were very much useful in explaining the nature of interactions existing in those observed transition states and we hope the same can also be extended to any transition state in an electrostatically driven reaction potential energy surface.

Possible interstellar formation of glycine from the reaction of CH2=NH, CO and H2O: catalysis by extra water molecules through the hydrogen relay transport.

"How the fundamental life elements are created in the interstellar medium (ISM)?" is one of the intriguing questions related to the genesis of life. Using computational calculations, we have discussed the reaction of CH2=NH, CO and H2O for the formation of glycine, the simplest life element. This reaction proceeds through a concerted mechanism with reasonably large barriers for the cases with one and two water molecules as reactants. For the two water case we found that the extra water molecule exhibits some catalytic role through the hydrogen transport relay effect and the barrier height is reduced substantially compared to the case with one water molecule. These two cases can be treated as ideal cases for the hot-core formation of the interstellar glycine. With an increasing number of water molecules as the reactants, we found that when the numbers of water molecules are three or more than three, the barrier height reduced so drastically that the transition states were more stable than the reactants. Such a situation gives a clear indication that with excess water molecules as the reactants, this reaction will be feasible even under the low temperature conditions existing in the cold interstellar clouds and the exothermic nature of the reaction will be the driving force.

Reaction between HN and SN: a possible channel for the interstellar formation of N2 and SH in the cold interstellar clouds.

Using computational calculations the potential energy surface (PES) of the reaction between NH and NS has been analysed. The PES of the reaction shows the formation of two very stable species, HNSN and HNNS. Out of these two, HNNS which has the signature N-N linkage was found to be the most stable species in the PES. In view of the highly exothermic nature of the reaction surface, it has been proposed that these two species can possibly be detected in the interstellar space. For the first time it has also been shown that the reaction between the NH and NS can lead to the possible formation of N2via the isomer HNNS, and how the effect of tunnelling can make this reaction very much feasible, even under the extremely low temperature conditions prevailing in the interstellar medium. Based on the already reported results, a similar kind of behaviour for the NH + NO reaction surface has also been proposed. These dissociation reactions leading to the formation of N2 can be considered as potential secondary contributing channels while accounting for the total estimates of N2 in the interstellar medium, and thus HNNS as well as HNNO can be considered as stable reservoir molecules for interstellar N2. Besides the formation of N2, the formation of another astronomically important radical, SH in the cold interstellar clouds, has also been proposed.

Potential application of compliance constants in predicting the mass spectral fragmentation of metabolites.

RATIONALE: Metabolomics is a qualitative and quantitative measurement of the metabolite content of any biological system under a given physiological status. Due to the chemically diverse nature of these samples, metabolite identification is a difficult task, and development of alternative approaches, such as those based on mass spectrometry (MS), aimed at proper metabolite identification is required. METHODS: Compliance constants, a direct measure of mechanical bond strength, were used for the first time to predict the MS fragmentation patterns of different regional isomers of a ubiquitous plant metabolite, caffeoylquinic acid (CQA). The compliance constant of an ester bond linking caffeic acid and a quinic acid molecule in CQA was calculated using density functional theory and Wilson's G-matrix formalism to distinguish the isomers. The predicted fragmentation patterns were compared with mass spectra obtained using negative ion electrospray ionization ultra-high-performance liquid chromatography quadrupole time-of-flight mass spectrometry (UHPLC/QTOFMS). RESULTS: Our in silico results were found to be in correlation with our UHPLC/QTOFMS results, suggesting a potential application of compliance constant algorithms for the rationalization of complex mass spectrometric data. The results also show that the different configurations in stereochemistry that exist between different regional isomers contribute to the underlying energy of the surrounding bonds and the fragmentation thereof. CONCLUSIONS: The results of our pilot study suggest that computational modelling can be applied for metabolite identification during metabolomic data mining and Natural Product research in general.

Role of large thermal fluctuations and magnesium ions in t-RNA selectivity of the ribosome.

The fidelity of translation selection begins with the base pairing of codon-anticodon complex between the m-RNA and tRNAs. Binding of cognate and near-cognate tRNAs induces 30S subunit of the ribosome to wrap around the ternary complex, EF-Tu(GTP)aa-tRNA. We have proposed that large thermal fluctuations play a crucial role in the selection process. To test this conjecture, we have developed a theoretical technique to determine the probability that the ternary complex, as a result of large thermal fluctuations, forms contacts leading to stabilization of the GTPase activated state. We argue that the configurational searches for such processes are in the tail end of the probability distribution and show that the probability for this event is localized around the most likely configuration. Small variations in the repositioning of cognate relative to near-cognate complexes lead to rate enhancement of the cognate complex. The binding energies of over a dozen unique site-bound magnesium structural motifs are investigated and provide insights into the nature of interaction of divalent metal ions with the ribosome.

Role of the oxyallyl substructure in the near infrared (NIR) absorption in symmetrical dye derivatives: A computational study.

It is well-known from experimental studies that the oxyallyl-substructure-based squarylium and croconium dyes absorb in the NIR region of the spectrum. Recently, another dye has been reported (J. Am. Chem. Soc. 2003, 125, 348) which contains the same basic chromophore, but the absorption is red-shifted by at least 300 nm compared to the former dyes and is observed near 1100 nm. To analyze the reasons behind the large red shift, in this work we have carried out symmetry-adapted cluster-configuration interaction (SAC-CI) studies on some of these NIR dyes which contain the oxyallyl substructure. From this study, contrary to the earlier reports, it is seen that the donor groups do not seem to play a major role in the red-shift of the absorption. On the other hand, on the basis of the results of the high-level calculations carried out here and using qualitative molecular orbital theory, it is observed that the orbital interactions play a key role in the red shift. Finally, design principles for the oxyallyl-substructure-based NIR dyes are suggested.

Near-infrared absorption in symmetric squarylium and croconate dyes: a comparative study using symmetry-adapted cluster-configuration interaction methods.

Symmetric croconate (CR) and squarylium dyes (SQ) are well-known near-infrared (NIR) dyes and, in general, are considered to be donor-acceptor-donor type molecules. It is established in the literature that CR dyes absorb in a longer wavelength region than the corresponding SQ dyes. This has been attributed to the CR ring being a better acceptor than the SQ ring. Thus increasing the donor capacity should lead to a bathochromic shift in both SQ and CR. On the other hand, some experiments reported in the literature have revealed that increasing the conjugation in the donor part of the SQ molecule leads first to red shift, which upon a further increase of the conjugation changes to a blue shift. Hence, to understand the role of the central ring and the substitutions in the absorption of these dyes, we carried out high-level symmetry-adapted cluster-configuration interaction (SAC-CI) calculations of some substituted SQ and CR dyes and compare the absorption energy with the existing experimental data. We found that there is very good agreement. We also carried out SAC-CI calculations of some smaller model molecules, which contain the main oxyallyl substructure. We varied the geometry (angle) of the oxyallyl subgroup and the substitution in these model molecules to establish a correlation with the bathochromic shift. We found that the charge transfer is very small and does not play the key role in the red shift, but on the other hand, the perturbation of the HOMO-LUMO gap (HLG) from both the geometry and substitution seems to be responsible for this shift. We suggest as a design principle that increasing the donor capacity of the groups may not help in the red shift, but introducing groups which perturb the HLG and decrease it without changing the MO character should lead to a larger bathochromic shift.