Uncovering the mechanism of selective stabilization of high-energy diastereoisomers via inclusion.

Supramolecular systems may be used to stabilize otherwise unstable isomers to find alternative synthetic pathways. It has been reported that cucurbit[8]uril can stabilize trans-I and trans-II CuII cyclam, whereas trans-III is the only non-substituted trans CuII cyclam diastereoisomer found outside of the host molecule experimentally. Quantum chemistry methods can provide valuable insight into the intermolecular interactions involved in these inclusion complexes. All five possible trans diastereoisomers of CuII cyclam were studied within the host molecule to calculate the interaction energy and free energy of association for each complex. The relative free energies of the five free cyclams confirm that trans-I and trans-II are the most energetically accessible diastereoisomers from the initial trans-III starting point. Energy decomposition analysis was used to identify the attractive and repulsive interactions between cyclam and cucurbit[8]uril and showed that trans-II encounters repulsive forces almost three times greater than trans-I, which may explain the 7:3 ratio of trans-I to trans-II within cucurbit[8]uril that occurs experimentally. Optimized complex geometries with trans-III, IV, and V show that the cyclams protrude out of cucurbit[8]uril, whereas trans-I and trans-II become more encapsulated and elongate the host, suggesting that the position of the cyclam is extremely important when forming non-covalent interactions. Our results agree with the experimental findings and provide greater insight into why the most stable isolated cyclam diastereoisomer, trans-III, does not form a complex. Supplementary Information: The online version contains supplementary material available at 10.1007/s00214-023-03077-7.

A shared-weight neural network architecture for predicting molecular properties.

Quantum chemical methods scale poorly with increasing molecular size and machine learning models have emerged as a promising, computationally-efficient alternative. We present a shared-weight neural network architecture based on modified atom-centered symmetry functions (ACSFs) and show that it performs similarly to the more computationally expensive per-element neural networks of previous work with ACSFs. The model achieves chemically accurate predictions, with a mean absolute error as low as 0.63 kcal mol-1 on energy predictions in the QM9 data set. Additionally, we show that it can reliably predict atomic forces.

Exploring electron pair behaviour in chemical bonds using the extracule density.

We explore explicit electron pair behaviour within the chemical bond (and lone pairs) by calculating the probability distribution for the center-of-mass (extracule) of an electron pair described by single localized orbitals. Using Edmiston-Ruedenberg localized orbitals in a series of 61 chemical systems, we demonstrate the utility of the extracule density as an interpretive tool in chemistry. By accessing localized regions of chemical space we simplify the interpretation of the extracule density and afford a quantum mechanical interpretation of "chemically intuitive" features of electronic structure. Specifically, we describe the localized effects on chemical bonds due to changes in electronegativities of bonded neighbours, bond strain, and non-covalent interactions. We show that the extracule density offers unique insight into electronic structure and allows one to readily quantify the effects of changing the chemical environment.

Assessment and Application of Density Functional Theory for the Prediction of Structure and Reactivity of Vanadium Complexes.

We assess the performance of six density functionals, each paired with one of five basis sets (a total of 30 model chemistries) for the prediction of geometrical parameters in the coordination sphere of nine vanadium complexes (for a total of 270 structural analyses). We find that results are generally consistent over the range of functionals tested and that none fail drastically. For bond lengths, the model chemistry PBE0/QZVP performed the best overall (having a MAD of only 0.02 Å from experiment) yet PBE0/6-31G* provides nearly identical results. For bond angles, PBE0 also performed best overall and, when combined with the 6-31G* basis, produces one of the smallest error distributions of any model chemistry tested. We subsequently applied the PBE0/6-31G* model chemistry to understanding the mechanism of action of a [BIMPY]VCl3 catalyst in the polymerization of styrene (Sty) and vinyl acetate (VAc). Our results indicate that the [BIMPY]VCl3 catalyst operates through a unique, two-step reaction pathway: dehalogenation to form a reactive V(II) intermediate (a highly favorable process) followed by a potentially reversible OMRP to control the polymerization of vinyl acetate. Control over vinyl acetate is facilitated by both the higher reactivity of the radical species and the participation of the ester group in the trapping step. In both the Sty and VAc cases we predict relatively poor control of the polymerization with the vanadium catalyst, which is in good agreement with our experimental results.

Predicting bond strength from a single Hartree-Fock ground state using the localized pair model.

We present an application of the recently introduced Localized Pair Model (LPM) [Z. A. Zielinksi and J. K. Pearson, Comput. Theor. Chem., 2013, 1003, 7990] to characterize and quantify properties of the chemical bond in a series of substituted benzoic acid molecules. By computing interelectronic distribution functions for doubly-occupied Edmiston-Ruedenberg localized molecular orbitals (LMOs), we show that chemically intuitive electron pairs may be uniquely classified and bond strength may be predicted with remarkable accuracy. Specifically, the HF/u6-311G(d,p) level (where u denotes a complete uncontraction of the basis set) is used to generate the relevant LMOs and their respective interelectronic distribution functions can be linearly correlated to the well-known Hammett σp or σm parameters with near-unity correlation coefficients.

Machine Learned Composite Methods for Electronic Structure Theory.

Because of the prohibitive scaling of ab initio techniques for modeling chemical species with high accuracy, they are not generally tractable for large systems. It is therefore of considerable interest to develop high-accuracy computational models with low computational cost that can afford predictions of electronic structure and properties of macromolecular species. Composite methods, as first introduced by Pople [Pople, J. A.; Head-Gordon, M.; Fox, D. J.; Raghavachari, K.; Curtiss, L. A. J. Chem. Phys.1989, 90, 5622.], are an intuitive solution to this problem as they seek to systematically increase accuracy in model chemistries by taking advantage of favorable error cancellation among reasonably low-cost models. By linearly combining a series of carefully chosen model chemistries, the result of a prohibitive-scaling correlated model chemistry with a large basis set may be approximated with relatively good fidelity. However, the full extent to which the choice of low-cost models dictates the predictive accuracy of composite methods is not known, and a full exploration of all model chemistries would be advantageous for the design and validation of a generalizable composite method for widespread application. Here, we show that remarkable accuracy can be generally achieved with composite methods that are more judiciously constructed, leading to increased accuracy with significantly reduced computational cost. By designing a systematic procedure for the automated generation and assessment of over 10 billion unique composite methods, we have extensively explored the space of modern model chemistries to elucidate important design principles in the construction of reliable composite procedures. We anticipate our work to be the starting point in the pursuit of creative approaches to modeling large chemical systems with high accuracy by using novel combinatorial modeling.

A simultaneous probability density for the intracule and extracule coordinates.

We introduce the intex density X(R,u), which combines both the intracular and extracular coordinates to yield a simultaneous probability density for the position of the center-of-mass radius (R) and relative separation (u) of electron pairs. One of the principle applications of the intex density is to investigate the origin of the recently observed secondary Coulomb hole. The Hartree-Fock (HF) intex densities for the helium atom and heliumlike ions are symmetric functions that may be used to prove the isomorphism 2I(2R)=E(R), where I(u) is the intracule density and E(R) is the extracule density. This is not true of the densities that we have constructed from explicitly correlated wave functions. The difference between these asymmetric functions and their symmetric HF counterparts produces a topologically rich intex correlation hole. From the intex hole distributions (X(exact)(R,u)-X(HF)(R,u)), we conclude that the probability of observing an electron pair with a very large interelectronic separation increases with the inclusion of correlation only when their center-of-mass radius is close to half of their separation.

Isolation, biomimetic synthesis, and cytotoxic activity of bis(pseudopterane) amines.

LC-MS/MS-based screening of the dichloromethane extract of the gorgonian coral Pseudopterogorgia acerosa led to the isolation of a novel bis(pseudopterane) amine (1). The structural assignment of 1 was achieved by 1D and 2D NMR and mass spectrometry analysis. A biomimetic synthesis of 1 and the known symmetrical diterpene 2 from pseudopterolide (3) is described in this report. Bis(pseudopterane) amine showed selective growth inhibition activity against cancer cell lines with IC(50) values of 4.2 microM (HCT116) and 42 microM (HeLa).

Intracule functional models. IV. Basis set effects.

We have calculated position and dot intracules for a series of atomic and molecular systems, starting from an unrestricted Hartree-Fock wave function, expanded using the STO-3G, 6-31G, 6-311G, 6-311++G, 6-311++G(d,p), 6-311++G(3d,3p), and 6-311++G(3df,3pd) basis sets as well as the nonpolarized part of Dunning's cc-pV5Z basis. We find that the basis set effects on the intracules are small and that correlation energies from the dot intracule ansatz are remarkably insensitive to the basis set quality. Mean absolute errors in correlation energies across the G1 data set agree to within 2 mE(h) for all basis sets tested.

Effect of substituents on the GPx-like activity of ebselen: steric versus electronic.

The direct oxidation of ebselen and several derivatives by hydrogen peroxide is investigated using the B3LYP/6-31G(d,p) method to elucidate the effects of substituents on GPx-like activity. While previous studies have attributed the differences in GPx activity of substituted ebselen compounds to the electronic nature of the substituents, the influence of functional groups is poorly understood. The effects of various solvents are incorporated by employing the CPCM method. It is shown that a substituent in the ortho position to the selenium atom sterically hinders attack of a nucleophile at selenium and thus increases the barrier to reaction. The observed increase in GPx-like activity of an ebselen derivative with an ortho substituent is explained by the fact that the steric hindrance prevents thiol exchange reactions.

Understanding the electronic structure, reactivity, and hydrogen bonding for a 1,2-diphosphonium dication.

The experimental charge density for hexamethyldiphosphonium ditriflate has been determined from low-temperature high-resolution X-ray diffraction data. These results have been compared with theoretically calculated values for the isolated gas-phase compound. Analysis of the topological and atomic basin properties has provided insight into the exact nature of the P-P bond in both the crystalline and the gas-phase structures. The rho(b)(r) and nabla2rho(b)(r) values highlight the covalent nature of the P-P bond, while the atomic charges indicate a localization of the positive charges on the two phosphorus atoms. This seems to indicate that a covalent bond is formed despite a strong electrostatic repulsion between these two heteroatoms. The topological properties and electrostatic potentials have also been shown to provide significant insight into the chemical reactivity of the title compound. A topological analysis of P2Me4, P2Me5(+), and P2Me6(+2) species has provided information about the progression of the P-P bond in the synthesis of the title compound. An investigation of the different hydrogen-bonding networks present in the crystalline and gas-phase structures, along with their affect on the electronic structure of the title compound has also been investigated. This has all led to significant new insight into the electronic structure, reactivity, and weak hydrogen bonding in prototypical 1,2-diphosphonium dications.

Novel Bonding Mode in Phosphine Haloboranes.

We have predicted, using a wide range of theoretical models, the potential energy surfaces of dative bond stretching in some phosphine haloboranes and closely associated analogues. It is shown that these dative complexes demonstrate unusual bond stretching potentials that are characterized by having multiple inflection points and are not able to be fit to any traditional Morse or Lennard-Jones-type curve. Specifically, in the case of Cl3B-PH3, this effect is so pronounced that the surface actually exhibits two distinct minima. To the best of our knowledge, this is the first example of such a unique bonding phenomenon reported for these species and is explained by the competition between the energetic cost of the required pyramidalization of the Lewis acid to form a dative bond and the stabilization from the favorable attraction between the Lewis acid and base. When the cost of pyramidalization of the Lewis acid is high relative to the strength of the interaction between the acid and base, the potential well associated with dative bonding is significantly weakened and the result is a relatively flat potential energy surface that is susceptible to the unusual characteristics described herein.

Isolation and structure elucidation of satosporin A and B: new polyketides from Kitasatospora griseola.

Satosporins A and B, two novel glucosylated polyketides, were isolated from the actinomycete Kitasatospora griseola MF730-N6. The polyketides possess an unprecedented tricyclic ring system that was fully characterized using a combination of spectroscopic analyses and computational calculations. Satosporin A was quantitatively converted into its aglycon homologue, satosporin C, using a ß-glucosidase. The determination of the absolute stereochemistry was achieved using solution TDDFT/ECD calculations and chemical derivatization methods.

Density functional theory study of the reaction mechanism and energetics of the reduction of hydrogen peroxide by ebselen, ebselen diselenide, and ebselen selenol.

Density functional theory calculations at the B3LYP/6-311++G(3df,3pd)//B3LYP/6-31G(d,p) level have been performed to elucidate the mechanism and reaction energetics for the reduction of hydrogen peroxide by ebselen, ebselen diselenide, ebselen selenol, and their sulfur analogues. The effects of solvation have been included with the CPCM model, and in the case of the selenol anion reaction, diffuse functions were used on heavy atoms for the geometry optimizations and thermochemical calculations. The topology of the electron density in each system was investigated using the quantum theory of atoms in molecules, and a detailed interpretation of the electronic charge and population data as well as the atomic energies is presented. Reaction free energy barriers for the oxidation of ebselen, ebselen diselenide, and ebselen selenol are 36.8, 38.4, and 32.5 kcal/mol, respectively, in good qualitative agreement with experiment. It is demonstrated that the oxidized selenium atom is significantly destabilized in all cases and that the exothermicity of the reactions is attributed to the peroxide oxygen atoms via reduction. The lower barrier to oxidation exhibited by the selenol is largely due to entropic effects in the reactant complex.

Modeling the reduction of hydrogen peroxide by glutathione peroxidase mimics.

Theoretical calculations have been performed on three model reactions representing the reduction of hydrogen peroxide by ebselen, ebselen selenol, and ebselen diselenide. The reaction surfaces have been investigated at the B3PW91/6-311G(2df,p) level, and single-point energies were calculated using the 6-311++G(3df,3pd) basis set. Solvent effects were included implicitly with the conductor-like polarizable continuum model and in one case with explicit inclusion of three water molecules. Mechanistic information is gained from investigating the critical points using the quantum theory of atoms in molecules. The barriers for the reduction of hydrogen peroxide with the ebselen, ebselen selenol, and ebselen diselenide models are 56.7, 53.4, and 35.3 kcal/mol, respectively, suggesting that ebselen diselenide may be the most active antioxidant in the ebselen GPx redox pathway. Results are also compared to that of the sulfur analogues of the model compounds.

An evaluation of various computational methods for the treatment of organoselenium compounds.

A reliable computational method for the prediction of organoselenium geometries and bond dissociation energies (BDEs) has been determined on the basis of the performance of density functional theory (DFT: B3LYP and B3PW91) and ab initio molecular orbital procedures (Hartree-Fock (HF)) in conjunction with various Pople basis sets including (but not limited to) the 6-31G(d), 6-31G(d,p), 6-311G(d), 6-311G(d,p), 6-311G(2df,p), and 6-311G(3df,3pd) sets. Predicted geometries and BDEs are compared with available experimental data and quadratic configuration interaction including single and double substitutions (QCISD) results. The B3PW91/6-311G(2df,p) level of theory is recommended for the prediction of the geometries and energetics of organoselenium compounds.