Computational NMR Study of Benzothienoquinoline Heterohelicenes.

Early NMR studies of several heterohelicenes containing an annular nitrogen atom and a thiophene ring in their structure suggested the possibility of the lengthening of the carbon-carbon bonds in the interior of the helical turn of the molecule based on the progressive upfield shift of 13C resonances toward the center of the helical turn. We now report a comprehensive analysis of the optimized geometry and a comparison of the calculated vs. observed 1H and 13C NMR chemical shifts of nineteen representative benzothienoquinoline heterohelicenes. As was initially hypothesized on the basis of the progressive upfield shift of carbon resonances toward the center of the interior helical turn, the present computational study has demonstrated that carbon-carbon bonds indeed have more sp3 character and are longer than normal sp2 bonds to accommodate the helical twist of the molecule, as expected.

Stereochemical and Computational NMR Survey of 1,2,3-Triazoles: in Search of the Original Tauto-Conformers.

The conformational analysis of nine functionalized 1,2,3-triazoles was carried out by the correlation of calculated and experimental high-level nuclear magnetic resonance (NMR) chemical shifts. In solution, the studied triazoles are in exchange dynamic equilibrium caused by their prototropic tautomerism of the NH-proton. The experimentally unresolved NMR signals were assigned for most of the compounds. A more thorough survey was conducted for 4-t-butyl-1,2,3-triazole-5-carbaldehyde oxime. The analysis performed within the framework of the DP4+ formalism completely confirmed the hypothesis of the predominance of the 2H-tautomer. Thus, the methodology for estimating stereochemical structures in the absence of some experimental data allowed the most stable conformations for dynamic systems with different tautomeric ratios to be reliably identified.

On the Efficiency of the Density Functional Theory (DFT)-Based Computational Protocol for 1H and 13C Nuclear Magnetic Resonance (NMR) Chemical Shifts of Natural Products: Studying the Accuracy of the pecS-n (n = 1, 2) Basis Sets.

The basis set issue has always been one of the most important factors of accuracy in the quantum chemical calculations of NMR chemical shifts. In a previous paper, we developed new pecS-n (n = 1, 2) basis sets purposed for the calculations of the NMR chemical shifts of the nuclei of the most popular NMR-active isotopes of 1-2 row elements and successfully approbated these on the DFT calculations of chemical shifts in a limited series of small molecules. In this paper, we demonstrate the performance of the pecS-n (n = 1, 2) basis sets on the calculations of as much as 713 1H and 767 13C chemical shifts of 23 biologically active natural products with complicated stereochemical structures, carried out using the GIAO-DFT(PBE0) approach. We also proposed new alternative contraction schemes for our basis sets characterized by less contraction depth of the p-shell. New contraction coefficients have been optimized with the property-energy consistent (PEC) method. The accuracies of the pecS-n (n = 1, 2) basis sets of both the original and newly contracted forms were assessed on massive benchmark calculations of proton and carbon chemical shifts of a vast variety of natural products. It was found that less contracted pecS-n (n = 1, 2) basis sets provide no noticeable improvement in accuracy. These calculations represent the most austere test of our basis sets as applied to routine calculations of the NMR chemical shifts of real-life compounds.

Stereochemical Study of the Super Large Tetrakis Alkaloid Alasmontamine A by Means of an Advanced Computational NMR.

1H and 13C NMR chemical shifts of the tetrakis monoterpene indole alkaloid alasmontamine A, with a molecular formula of C84H91N8O12, have been calculated within the DFT framework. Six minimum energy conformers of this alkaloid were identified, and three key configurations that contribute to its NMR shielding constants were established. Several ambiguities in the reported assignment of the NMR chemical shifts of alasmontamine A have been resolved.

Reactions of α-Functionally Substituted Enals with Terminal Alkynes: Unexpected Assembly of 2-Amino-2-Cyclopentenones.

The reactions of α-chalcogen, α-halo, or α-amino functionally substituted enals with terminal alkynes leading either to corresponding propargyl alcohols (for O-, S-, Cl-, and Br-bearing substrates) or unexpected 2-amino-2-cyclopentenones (for aminoenals) are described. The key feature of these reactions is the rearrangement of adducts bearing amino groups on silica gel that triggered further cyclization to five-membered carbocycles.

Computational studies of [1]benzothieno[2,3-c]naphtho[1,2-f]quinoline.

Early NMR studies of several heterohelicenes containing an annular nitrogen atom and a thiophene ring in their structure suggested the possibility of the lengthening of the carbon-carbon bonds in the interior of the helical turn of the molecule based on the progressive more shielded nature of 13 C resonances toward the center of the helical turn. Computational chemistry capabilities when those NMR studies were performed were primitive in comparison to what is now possible. We now report the optimized geometry and a comparison of the calculated versus observed 1 H and 13 C NMR chemical shift assignments for [1]benzothieno[2,3-c]naphtho[1,2-f]quinoline that confirms these suspicions.

The Development of Pharmacophore Models for the Search of New Natural Inhibitors of SARS-CoV-2 Spike RBD-ACE2 Binding Interface.

To date, some succeeding variants of SARS-CoV-2 have become more contagious. This virus is known to enter human cells by binding the receptor-binding domain (RBD) of spike protein with the angiotensin-converting enzyme 2 (ACE2), the latter being a membrane protein that regulates the renin-angiotensin system. Since the host cell receptor plays a critical role in viral entry, inhibition of the RBD-ACE2 complex is a promising strategy for preventing COVID-19 infection. In the present communication, we propose and utilize an approach based on the generation of a complex of pharmacophore models and subsequent Induced Fit Docking (IFD) to identify potential inhibitors of the main binding sites of the Omicron SARS-CoV-2 RBD(S1)-ACE2 complex (PDB ID: 7T9L) among a number of natural products of various types and origins. Several natural compounds have been found to provide a high affinity for the receptor of interest. It is expected that the present results will stimulate further research aimed at the development of specialized drugs against this virus.

Four-Component Relativistic Calculations of NMR Shielding Constants of the Transition Metal Complexes-Part 2: Nitrogen-Coordinated Complexes of Cobalt.

Both four-component relativistic and nonrelativistic computations within the GIAO-DFT(PBE0) formalism have been carried out for 15N and 59Co NMR shielding constants and chemical shifts of a number of the nitrogen-coordinated complexes of cobalt. It was found that the total values of the calculated nitrogen chemical shifts of considered cobalt complexes span over a range of more than 580 ppm, varying from -452 to +136 ppm. At that, the relativistic corrections to nitrogen shielding constants and chemical shifts were demonstrated to be substantial, changing accordingly from ca. -19 to +74 ppm and from -68 to +25 ppm. Solvent effects on 15N shielding constants and chemical shifts were shown to have contributions no less important than the relativistic effects, namely from -35 to +63 ppm and from -74 to +23 ppm, respectively. Cobalt shielding constants and chemical shifts were found to vary in the ranges of, accordingly, -20,157 to -11,373 ppm and from +3781 to +13,811. The relativistic effects are of major importance in the cobalt shielding constants, resulting in about 4% for the shielding-type contributions, while solvent corrections to cobalt shielding constants appeared to be of less significance, providing corrections of about 1.4% to the gas phase values.

Digitization of the electron shell via the localized orbital locator formalism: trends in the size and electronegativity changes of atoms across the periodic table.

The parameters of the (3, -3) critical point in the topology of the localized orbital locator inside the electron shell reveal patterns that make it possible to recognize trends in the size, electronegativity, and ionization energy of atoms in the first four periods of the periodic table.

Combined Computational NMR and Molecular Docking Scrutiny of Potential Natural SARS-CoV-2 Mpro Inhibitors.

In continuation of the search for potential drugs that inhibit the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), in this work, a combined approach based on the modeling of NMR chemical shifts and molecular docking is suggested to identify the possible suppressors of the main protease of this virus among a number of natural products of diverse nature. Primarily, with the aid of an artificial neural network, the problem of the reliable determination of the stereochemical structure of a number of studied compounds was solved. Complementary to the main goal of this study, theoretical modeling of NMR spectral parameters made it feasible to perform a number of signal reassignments together with introducing some missing NMR data. Finally, molecular docking formalism was applied to the analysis of several natural products that could be chosen as prospective candidates for the role of potential inhibitors of the main protease. The results of this study are believed to assist in further research aimed at the development of specific drugs based on the natural products against COVID-19.

Computational NMR of natural products: On the way to super large molecules exemplified with alasmontamine A.

1 H and 13 C nuclear magnetic resonance (NMR) chemical shifts of a tetrakis monoterpene indole alkaloid alasmontamine A with a molecular formula of C84 H91 N8 O12 have been calculated at the PBE0/pcSseg-2//pcseg-2 level of theory on M06-2X/aug-cc-pVDZ geometry. In the course of the preliminary conformational search, six true minimum energy conformers were identified that can contribute to the actual conformation of this huge alkaloid. Calculated chemical shifts generally demonstrated a good agreement with available experimental data characterized with a corrected mean absolute error of 0.10 ppm for the range of about 7 ppm for protons and 1.1 ppm for the range of about 160 ppm for carbons.

Four-component relativistic calculations of NMR shielding constants of the transition metal complexes. Part 1: Pentaammines of cobalt, rhodium, and iridium.

The nonrelativistic and four-component fully relativistic calculations of 1 H, 15 N, 59 Co, 103 Rh, and 193 Ir shielding constants of pentaammineaquacomplexes of cobalt(III), rhodium(III), and iridium(III) were carried out at the density functional theory (DFT) level of theory. The noticeable deshielding relativistic corrections were observed for nitrogen shielding constants (chemical shifts), whereas those corrections were found to be negligible for protons. For the transition metals cobalt, rhodium, and iridium, relativistic corrections to their nuclear magnetic resonance (NMR) shielding constants were found to be rather small for cobalt and rhodium (some 5-10%), whereas they are essentially larger for iridium (up to 70%).

Simple and Versatile Scheme for the Stereochemical Identification of Natural Products and Diverse Organic Compounds with Multiple Asymmetric Centers.

A simple and versatile scheme of stereochemical identification of the stereochemically rich natural products and organic compounds with multiple asymmetric centers is proposed based on a detailed parsing of calculated 1H and 13C NMR chemical shifts in combination with their DP4+ analysis, as exemplified for three natural products: sungucine, physalin D, and anabsinthin. Performed benchmark calculations of the considered diastereomers provided very good agreement with their known experimental stereochemical structures.

Localized orbital locator as a descriptor for quantification and digital presentation of lone pairs: benchmark calculations of 4-substituted pyridines.

The parameters of the (3,-3) critical point in the localized orbital locator topology near a heteroatom have been found to reflect the changes in the size, density and electron energy of the lone pair and correlate with the donor ability of the lone pair carrying heteroatom.

Computational 1 H and 13 C NMR of the trimeric monoterpenoid indole alkaloid strychnohexamine: Selected spectral updates.

Very large trimeric indole alkaloid strychnohexamine, with empirical formula C59 H60 N6 O (66 second-row atoms and 60 protons), has been subjected to the state-of-the-art computation of the 1 H and 13 C nuclear magnetic resonance (NMR) chemical shifts of its configurational isomers at each of the 14 asymmetric centers. Several spectral reassignments and corrections of 1 H and 13 C NMR spectra of this alkaloid were suggested based on the PBE0/pcSseg-2//pcseg-2 calculation of its NMR chemical shifts. Thus, all pairs of diastereotopic protons were assigned together with four aromatic carbon resonances of C-9 and C-11, C-9″, and C-11″. In addition, the unassigned chemical shifts of carbon C-23″ and proton at C-3' in, accordingly, 13 C and 1 H NMR spectra were predicted.

Computational 1H and 13C NMR of strychnobaillonine: On the way to larger molecules calculated at lower computational costs.

The 1H and 13C NMR chemical shifts of strychnobaillonine, a very large dimeric indole alkaloid, consisting of as many as 46 nonhydrogen atoms, were calculated with using the established earlier the most effective computational protocol, PBE0/pcSseg-2//pcseg-2. A very good result was achieved at this level, characterized by the root mean square deviation of only 0.14 ppm for protons and 2.4 ppm for carbons, which enabled the verification of the configurations of its all 13 asymmetrical centers. Essential deviations of the calculated and experimental 1H NMR spectrum of strychnobaillonine were established in several cases, which enabled the performance of some additional NMR assignments and reassignments of the originally proposed structure.

A New Basis Set for the Calculation of 13C NMR Chemical Shifts within a Non-empirical Correlated Framework.

A new basis set of triple-ζ quality for carbon, 3z-S, is developed and tested at the DFT, MP2, and CCSD(T) levels, taking into account solvent and vibrational corrections for a number of molecules ranging from the smallest fluoromethane, CH3F, to the largest 5,10,15,20-tetraphenylporphyrin, C44H30N4. The proposed highly economical 3z-S basis set has been proven to provide very good accuracy in all examinations, comparable to that of the NMR-oriented Jensen's pcS-2 basis set, which is about 50% larger than 3z-S.

The 1 H and 13 C NMR chemical shifts of Strychnos alkaloids revisited at the DFT level.

The density functional theory calculation of 1 H and 13 C NMR chemical shifts in a series of ten 10 classically known Strychnos alkaloids with a strychnine skeleton was performed at the PBE0/pcSseg-2//pcseg-2 level. It was found that calculated 1 H and 13 C NMR chemical shifts provided a markedly good correlation with experiment characterized by a mean absolute error of 0.08 ppm in the range of 7 ppm for protons and 1.67 ppm in the range of 150 ppm for carbons, so that a mean absolute percentage error was as small as ~1% in both cases.

DFT computational schemes for 1 H and 13 C NMR chemical shifts of natural products, exemplified by strychnine.

A number of computational schemes based on different Density Functional Theory (DFT) functionals in combination with a number of basis sets were tested in the calculation of 1 H and 13 C NMR chemical shifts of strychnine, as a typical representative of the vitally important natural products, and used as a challenging benchmark and a rigorous test for such calculations. It was found that the most accurate computational scheme, as compared with experiment, was PBE0/pcSseg-4//pcseg-3 characterized by a mean absolute error of 0.07 ppm for the range of about 7 ppm for 1 H NMR chemical shifts and that of only 1.13 ppm for 13 C NMR chemical shifts spread over the range of about 150 ppm. For more practical purposes, including investigation of larger molecules from this series, a much more economical computational scheme, PBE0/pcSseg-2//pcseg-2, characterized by almost the same accuracy and much less computational demand, was recommended.