Vitifolignans, 3,4-dibenzyltetrahydrofuran lignans from Anemone vitifolia Buch.-Ham

Anemone vitifolia is a small herb found in Asia that is used to treat a range of diseases in Chinese traditional medicine. GNPS-based molecular networking of an Anemone vitifolia specimen revealed the presence of a network containing numerous ions indicating the presence of lignans, several of which suggested that there might be previously undescribed compounds in the extract. Fractionation of the organic extract yielded five undescribed lignans, the vitifolignans, together with one known. The structures were identified based on extensive spectroscopic data analysis (NMR, HR-ESI-MS, and UV), coupling constant calculation and comparison with reported data. Their absolute configurations were determined by comparison of experimental ECD spectra with calculated spectra. Compounds 4/5 showed weak inhibition of LPS-induced NO production in mouse mononuclear macrophages.

chARpack: The Chemistry Augmented Reality Package

Off-loading visualization and interaction into virtual reality (VR) using head-mounted displays (HMDs) has gained considerable popularity in simulation sciences, particularly in chemical modeling. Because of its unique way of soft immersion, augmented reality (AR) HMD technology has even more potential to be integrated into the everyday workflow of computational chemists. In this work, we present our environment to explore the prospects of AR in chemistry and general molecular sciences: The chemistry in Augmented Reality package (chARpack). Besides providing an extensible framework, our software focuses on a seamless transition between a 3D stereoscopic view with true 3D interactions and the traditional desktop PC setup to provide users with the best setup for all tasks in their workflow. Using feedback from domain experts, we discuss our design requirements for this kind of hybrid working environment (AR + PC), regarding input, features, degree of immersion, and collaboration.

Chelate Complexes of 3d Transition Metal Ions─A Challenge for Electronic-Structure Methods?

Different electronic-structure methods were assessed for their ability to predict two important properties of the industrially relevant chelating agent nitrilotriacetic acid (NTA): its selectivity with respect to six different first-row transition metal ions and the spin-state energetics of its complex with Fe(III). The investigated methods encompassed density functional theory (DFT), the random phase approximation (RPA), coupled cluster (CC) theory, and the auxiliary-field quantum Monte Carlo (AFQMC) method, as well as the complete active space self-consistent field (CASSCF) method and the respective on-top methods: second-order N-electron valence state perturbation theory (NEVPT2) and multiconfiguration pair-density functional theory (MC-PDFT). Different strategies for selecting active spaces were explored, and the density matrix renormalization group (DMRG) approach was used to solve the largest active spaces. Despite somewhat ambiguous multi-reference diagnostics, most methods gave relatively good agreement with experimental data for the chemical reactions connected to the selectivity, which only involved transition-metal complexes in their high-spin state. CC methods yielded the highest accuracy followed by range-separated DFT and AFQMC. We discussed in detail that even higher accuracies can be obtained with NEVPT2, under the prerequisite that consistent active spaces along the entire chemical reaction can be selected, which was not the case for reactions involving Fe(III). A bigger challenge for electronic-structure methods was the prediction of the spin-state energetics, which additionally involved lower spin states that exhibited larger multi-reference diagnostics. Conceptually different, typically accurate methods ranging from CC theory via DMRG-NEVPT2 in combination with large active spaces to AFQMC agreed well that the high-spin state is energetically significantly favored over the other spin states. This was in contrast to most DFT functionals and RPA which yielded a smaller stabilization and some common DFT functionals and MC-PDFT even predicting the low-spin state to be energetically most favorable.

Conjugated Polymer/Recombinant Escherichia coli Biohybrid Systems for Photobiocatalytic Hydrogen Production

Biohybrid photocatalysts are composite materials that combine the efficient light-absorbing properties of synthetic materials with the highly evolved metabolic pathways and self-repair mechanisms of biological systems. Here, we show the potential of conjugated polymers as photosensitizers in biohybrid systems by combining a series of polymer nanoparticles with engineered Escherichia coli cells. Under simulated solar light irradiation, the biohybrid system consisting of fluorene/dibenzo [b,d]thiophene sulfone copolymer (LP41) and recombinant E. coli (i.e., a LP41/HydA BL21 biohybrid) shows a sacrificial hydrogen evolution rate of 3.442 mmol g-1 h-1 (normalized to polymer amount). It is over 30 times higher than the polymer photocatalyst alone (0.105 mmol g-1 h-1), while no detectable hydrogen was generated from the E. coli cells alone, demonstrating the strong synergy between the polymer nanoparticles and bacterial cells. The differences in the physical interactions between synthetic materials and microorganisms, as well as redox energy level alignment, elucidate the trends in photochemical activity. Our results suggest that organic semiconductors may offer advantages, such as solution processability, low toxicity, and more tunable surface interactions with the biological components over inorganic materials.

How to Decarbonize N-Heterocyclic Carbenes (NHCs): The simple Alane Adducts (NHC) ⋅ AlR3 (R=H, Me, Et)

The reaction of the amine-stabilized alane (NMe3) ⋅ AlH3 1 with the backbone-saturated N-heterocyclic carbene (NHC) SIDipp (SIDipp=1,3-bis-{2,6-di-iso-propyl-phenyl}-imidazolidin-2-ylidene) at 0 °C yielded the NHC alane adduct (SIDipp) ⋅ AlH3 2. Reaction at elevated temperatures or prolonged reaction at room temperature gave the product of a ring expansion reaction (RER) of the NHC, (NMe3) ⋅ AlH(RER-SIDippH2) 3 ⋅ (NMe3). Subsequent reaction of the latter with sterically less hindered NHCs (IMeMe {=1,3,4,5-tetramethyl-imidazolin-2-ylidene}, IiPrMe {=1,3-di-iso-propyl-4,5-dimethyl-imidazolin-2-ylidene}, and IiPr {=1,3-di-iso-propyl-imidazolin-2-ylidene}) afforded the NHC-stabilized RER-products (NHC) ⋅ AlH(RER-SIDippH2) 3 ⋅ (NHC) (NHC=IMeMe, IiPrMe, IiPr), while no reaction was observed with the sterically more demanding NHCs IDipp (=1,3-bis-{2,6-di-iso-propyl-phenyl}-imidazolin-2-ylidene), SIDipp and ItBu (=1,3-di-tert-butyl-imidazolin-2-ylidene). The compounds 3 ⋅ (NHC) were also obtained starting from (SIDipp) ⋅ AlH3 2 and NHC at room temperature. Heating solutions of (SIDipp) ⋅ AlH3 2 without additional base to 95 °C resulted in decarbonization of the NHC and substitution of the carbene carbon atom with aluminum hydride under loss of ethene. Subsequent dimerization afforded cis-[AlH{μ-N(Dipp)CH2CH2N(Dipp)}]2 4_dimer. Heating solutions of the NHC-ligated aluminum alkyls (SIDipp) ⋅ AlR3 2R (R=Me, Et) to 145 °C instead led to complete scission of the NHC backbone with evolution of ethene and isolation of the dialkylaluminium(III) amidinates {DippNC(R)NDipp}AlR2 5R (R=Me, Et).

Mechanistic Investigation on the Regioselectivity of Electrochemical Co(II)-Catalyzed [2 + 2 + 2] Cycloaddition of Terminal Acetylenes

In this paper, the regioselectivity of electrochemical Co(II)-catalyzed [2 + 2 + 2] cycloaddition of terminal alkynes was investigated using density functional theory. We explored in detail the energy profiles for both 1,2,4- and 1,3,5-regioselectivity pathways and revealed the origin of the regioselectivity. Two kinds of conformational isomers derived from the different coordination modes of alkynes with cobaltacyclopentadiene have been found, which were formed through electrochemically mediated redox processes. The regioselectivity of the reaction depends on the two coordination modes. When the Co(II) center attacks α-C of the third alkyne, while β2-C in cyclopentadiene bonds to β-C of the alkyne, the reaction favors the formation of 1,2,4-products. In contrast, when the Co(II) center connects to β-C of the alkyne, it forms only the 1,3,5-products via [4 + 2] cycloaddition because of the steric repulsion between the bulky ligand on Co(II) and the phenyl group in the alkyne.

A Vision for the Future of Multiscale Modeling

The rise of modern computer science enabled physical chemistry to make enormous progresses in understanding and harnessing natural and artificial phenomena. Nevertheless, despite the advances achieved over past decades, computational resources are still insufficient to thoroughly simulate extended systems from first principles. Indeed, countless biological, catalytic and photophysical processes require ab initio treatments to be properly described, but the breadth of length and time scales involved makes it practically unfeasible. A way to address these issues is to couple theories and algorithms working at different scales by dividing the system into domains treated at different levels of approximation, ranging from quantum mechanics to classical molecular dynamics, even including continuum electrodynamics. This approach is known as multiscale modeling and its use over the past 60 years has led to remarkable results. Considering the rapid advances in theory, algorithm design, and computing power, we believe multiscale modeling will massively grow into a dominant research methodology in the forthcoming years. Hereby we describe the main approaches developed within its realm, highlighting their achievements and current drawbacks, eventually proposing a plausible direction for future developments considering also the emergence of new computational techniques such as machine learning and quantum computing. We then discuss how advanced multiscale modeling methods could be exploited to address critical scientific challenges, focusing on the simulation of complex light-harvesting processes, such as natural photosynthesis. While doing so, we suggest a cutting-edge computational paradigm consisting in performing simultaneous multiscale calculations on a system allowing the various domains, treated with appropriate accuracy, to move and extend while they properly interact with each other. Although this vision is very ambitious, we believe the quick development of computer science will lead to both massive improvements and widespread use of these techniques, resulting in enormous progresses in physical chemistry and, eventually, in our society.

Beyond Explored Functionals: A Computational Journey of Two-Photon Absorption

We present a thorough investigation into the efficacy of 19 density functional theory (DFT) functionals, relative to RI-CC2 results, for computing two-photon absorption (2PA) cross sections (σ2PA) and key dipole moments (|μ00|, |μ11|, |Δμ|, |μ01|) for a series of coumarin dyes in the gas-phase. The functionals include different categories, including local density approximation (LDA), generalized gradient approximation (GGA), hybrid-GGA (H-GGA), range-separated hybrid-GGA (RSH-GGA), meta-GGA (M-GGA), and hybrid M-GGA (HM-GGA), with 14 of them being subjected to analysis for the first time with respect to predicting σ2PA values. Analysis reveals that functionals integrating both short-range (SR) and long-range (LR) corrections, particularly those within the RSH-GGA and HM-GGA classes, outperform the others. Furthermore, the range-separation approach was found more impactful compared to the varying percentages of Hartree-Fock exchange (HF Ex) within different functionals. The functionals traditionally recommended for 2PA do not appear among the top 9 in our study, which is particularly interesting, as these top-performing functionals have not been previously investigated in this context. This list is dominated by M11, QTP variants, ωB97X, ωB97X-V, and M06-2X, surpassing the performance of other functionals, including the commonly used CAM-B3LYP.

Accurate Machine Learning for Predicting the Viscosities of Deep Eutectic Solvents

Deep eutectic solvents (DESs) are emerging as environmentally friendly designer solvents for mass transport and heat transfer processes in industrial applications; however, the lack of accurate tools to predict and thus control their viscosities under both a range of environmental factors and formulations hinders their general application. While DESs may serve as designer solvents, with nearly unlimited combinations, this unfortunately makes it experimentally infeasible to comprehensively measure the viscosities of all DESs of potential industrial interest. To assist in the design of DESs, we have developed several new machine learning (ML) models that accurately and rapidly predict the viscosities of a diverse group of DESs at different temperatures and molar ratios using, to date, one of the most comprehensive data sets containing the properties of over 670 DESs over a wide range of temperatures (278.15-385.25 K). Three ML models, including support vector regression (SVR), feed forward neural networks (FFNNs), and categorical boosting (CatBoost), were developed to predict DES viscosity as a function of temperature and molar ratio and contrasted with multilinear and two-factor polynomial regression baselines. Quantum chemistry-based, COSMO-RS-derived sigma profile (σ-profile) features were used as inputs for the ML models. The CatBoost model is excellent at externally predicting DES viscosity, as indicated by high R2 (0.99) and low root-mean-square-error (RMSE) and average absolute relative deviations (AARD) (5.22%) values for the testing data sets, and 98% of the data points lie within the 15% of AARD deviations. Furthermore, SHapley additive explanation (SHAP) analysis was employed to interpret the ML results and rationalize the viscosity predictions. The result is an ML approach that accurately predicts viscosity and will aid in accelerating the design of appropriate DESs for industrial applications.

Computational Investigations of the Detailed Mechanism of Reverse Intersystem Crossing in Inverted Singlet-Triplet Gap Molecules

Inverted singlet-triplet gap (INVEST) materials have promising photophysical properties for optoelectronic applications due to an inversion of their lowest singlet (S1) and triplet (T1) excited states. This results in an exothermic reverse intersystem crossing (rISC) process that potentially enhances triplet harvesting, compared to thermally activated delayed fluorescence (TADF) emitters with endothermic rISCs. However, the processes and phenomena that facilitate conversion between excited states for INVEST materials are underexplored. We investigate the complex potential energy surfaces (PESs) of the excited states of three heavily studied azaphenalene INVEST compounds, namely, cyclazine, pentazine, and heptazine using two state-of-the-art computational methodologies, namely, RMS-CASPT2 and SCS-ADC(2) methods. Our findings suggest that ISC and rISC processes take place directly between the S1 and T1 electronic states in all three compounds through a minimum-energy crossing point (MECP) with an activation energy barrier between 0.11 to 0.58 eV above the S1 state for ISC and between 0.06 and 0.36 eV above the T1 state for rISC. We predict that higher-lying triplet states are not populated, since the crossing point structures to these states are not energetically accessible. Furthermore, the conical intersection (CI) between the ground and S1 states is high in energy for all compounds (between 0.4 to 2.0 eV) which makes nonradiative decay back to the ground state a relatively slow process. We demonstrate that the spin-orbit coupling (SOC) driving the S1-T1 conversion is enhanced by vibronic coupling with higher-lying singlet and triplet states possessing vibrational modes of proper symmetry. We also rationalize that the experimentally observed anti-Kasha emission of cyclazine is due to the energetically inaccessible CI between the bright S2 and the dark S1 states, hindering internal conversion. Finally, we show that SCS-ADC(2) is able to qualitatively reproduce excited state features, but consistently overpredict relative energies of excited state structural minima compared to RMS-CASPT2. The identification of these excited state features elaborates design rules for new INVEST emitters with improved emission quantum yields.

Fluoride Abstraction Induced by Tris(pentafluoroethyl)difluorophosphorane: A Convenient Way to Synthesize Cationic N-Heterocyclic Carbene- and Cyclic (Alkyl)(amino)carbene-Ligated Copper Alkyne and Arene Complexes

We herein report the convenient synthesis of different N-heterocyclic carbene (NHC)- and cyclic (alkyl)(amino)carbene (cAAC)-ligated copper cations using the weakly coordinating tris(pentafluoroethyl)trifluorophosphate counterion (FAP anion, [(C2F5)3PF3]-). The reaction of the fluorido complexes [(carbene)CuF] (carbene = NHC, cAACMe) 2a-2f and the tris(pentafluoroethyl)difluorophosphorane (C2F5)3PF2 in the presence of alkynes or arenes led to fluoride transfer from Cu to the phosphorane with formation of the cationic transition metal complexes [(carbene)Cu(L)]+ and the weakly coordinating counteranion [(C2F5)3PF3]- (FAP). Using this method, the complexes [(IDipp)Cu(L)]+FAP- (IDipp = 1,3-bis(2,6-di-iso-propylphenyl)-imidazolin-2-ylidene; L = PhC≡CPh, 4d; PhC≡CMe, 5d), [(cAACMe)Cu(L)]+FAP- (cAACMe = 1-(2,6-di-iso-propylphenyl)-3,3,5,5-tetramethyl-pyrrolidin-2-ylidene; L = PhC≡CPh, 4f; PhC≡CMe, 5f), [(SIDipp)Cu(C6Me6)]+FAP- (6e), (SIDipp = 1,3-bis(2,6-di-iso-propylphenyl)-imidazolidine-2-ylidene), and [(cAACMe)Cu(C6Me6)]+FAP- (6f) have been synthesized and characterized. The complexes [(IDipp)Cu(C6Me6)]+FAP- (6d) and [(cAACMe)Cu(C6Me6)]+FAP- (6f) have been used as catalysts for the copper(I)-catalyzed cycloaddition of benzyl azide to terminal alkynes.

Interplay of halogen bonding and solvation in protein-ligand binding

Halogen bonding is increasingly utilized in efforts to achieve high affinity and selectivity of molecules designed to bind proteins, making it paramount to understand the relationship between structure, dynamics, and thermodynamic driving forces. We present a detailed analysis addressing this problem using a series of protein-ligand complexes involving single halogen substitutions - F, Cl, Br, and I - and nearly identical structures. Isothermal titration calorimetry reveals an increasingly favorable binding enthalpy from F to I that correlates with the halogen size and σ-hole electropositive character, but is partially counteracted by unfavorable entropy, which is constant from F to Cl and Br, but worse for I. Consequently, the binding free energy is roughly equal for Cl, Br, and I. QM and solvation-free-energy calculations reflect an intricate balance between halogen bonding, hydrogen bonds, and solvation. These advances have the potential to aid future drug design initiatives involving halogenated compounds.

Massively scalable workflows for quantum chemistry: BigChem and ChemCloud

Electronic structure theory, i.e., quantum chemistry, is the fundamental building block for many problems in computational chemistry. We present a new distributed computing framework (BigChem), which allows for an efficient solution of many quantum chemistry problems in parallel. BigChem is designed to be easily composable and leverages industry-standard middleware (e.g., Celery, RabbitMQ, and Redis) for distributed approaches to large scale problems. BigChem can harness any collection of worker nodes, including ones on cloud providers (such as AWS or Azure), local clusters, or supercomputer centers (and any mixture of these). BigChem builds upon MolSSI packages, such as QCEngine to standardize the operation of numerous computational chemistry programs, demonstrated here with Psi4, xtb, geomeTRIC, and TeraChem. BigChem delivers full utilization of compute resources at scale, offers a programable canvas for designing sophisticated quantum chemistry workflows, and is fault tolerant to node failures and network disruptions. We demonstrate linear scalability of BigChem running computational chemistry workloads on up to 125 GPUs. Finally, we present ChemCloud, a web API to BigChem and successor to TeraChem Cloud. ChemCloud delivers scalable and secure access to BigChem over the Internet.

Excited-state antiaromaticity relief drives facile photoprotonation of carbons in aminobiphenyls

A combined computational and experimental study reveals that ortho-, meta- and para-aminobiphenyl isomers undergo distinctly different photochemical reactions involving proton transfer. Deuterium exchange experiments show that the ortho-isomer undergoes a facile photoprotonation at a carbon atom via excited-state intramolecular proton transfer (ESIPT). The meta-isomer undergoes water-assisted excited-state proton transfer (ESPT) and a photoredox reaction via proton-coupled electron transfer (PCET). The para-isomer undergoes a water-assisted ESPT reaction. All three reactions take place in the singlet excited-state, except for the photoredox process of the meta-isomer, which involves a triplet excited-state. Computations illustrate the important role of excited-state antiaromaticity relief in these photoreactions.

Benchmark Study on the Calculation of 207Pb NMR Chemical Shifts

A benchmark set for the computation of 207Pb nuclear magnetic resonance (NMR) chemical shifts is presented. The PbS50 set includes conformer ensembles of 50 lead-containing molecular compounds and their experimentally measured 207Pb NMR chemical shifts. Various bonding motifs at the Pb center with up to seven bonding partners are included. Six different solvents were used in the measurements. The respective shifts lie in the range between +10745 and -5030 ppm. Several calculation settings are assessed by evaluating computed 207Pb NMR shifts for the use with different density functional approximations (DFAs), relativistic approaches, treatment of the conformational space, and levels for geometry optimization. Relativistic effects were included explicitly with the zeroth order regular approximation (ZORA), for which only the spin-orbit variant was able to yield reliable results. In total, seven GGAs and three hybrid DFAs were tested. Hybrid DFAs significantly outperform GGAs. The most accurate DFAs are mPW1PW with a mean absolute deviation (MAD) of 429 ppm and PBE0 with an MAD of 446 ppm. Conformational influences are small as most compounds are rigid, but more flexible structures still benefit from Boltzmann averaging. Including explicit relativistic treatments such as SO-ZORA in the geometry optimization does not show any significant improvement over the use of effective core potentials (ECPs).

QM/MM study of the catalytic reaction of aphid myrosinase

Brevicoryne brassicae, an aphid species, exclusively consumes plants from the Brassicaceae family and employs a sophisticated defense mechanism involving a myrosinase enzyme that breaks down glucosinolates obtained from its host plants. In this work, we employed combined quantum mechanical and molecular mechanical (QM/MM) calculations and molecular dynamics (MD) simulations to study the catalytic reaction of aphid myrosinase. A proper QM region to study the myrosinase reaction should contain the whole substrate, models of Gln-19, His-122, Asp-124, Asn-166, Glu-167, Lys-173, Tyr-180, Val-228, Tyr-309, Tyr-346, Ile-347, Glu-374, Glu-423, Trp-424, and a water molecule. The calculations show that Asp-124 and Glu-423 must be charged, His-122 must be protonated on NE2, and Glu-167 must be protonated on OE2. Our model reproduces the anomeric retaining characteristic of myrosinase and indicates that the deglycosylation reaction is the rate-determining step of the reaction. Based on the calculations, we propose a reaction mechanism for aphid myrosinase-mediated hydrolysis of glucosinolates with an overall barrier of 15.2 kcal/mol. According to the results, removing a proton from Arg-312 or altering it to valine or methionine increases glycosylation barriers but decreases the deglycosylation barrier.

A Theoretical Perspective on the Photochemistry of Boron-Nitrogen Lewis Adducts

Boron-Nitrogen (B-N) Lewis adducts form a versatile family of compounds with numerous applications in functional molecules. Despite the growing interest in this family of compounds for optoelectronic applications, little is currently known about their photophysics and photochemistry. Even the electronic absorption spectrum of ammonia borane, the textbook example of a B-N Lewis adduct, is unavailable. Given the versatility of the light-induced processes exhibited by these molecules, we propose in this work a detailed theoretical study of the photochemistry and photophysics of simple B-N Lewis adducts. We used advanced techniques in computational photochemistry to identify and characterize the possible photochemical pathways followed by ammonia borane and extended this knowledge to the substituted B-N Lewis adducts pyridine-borane and pyridine-boric acid. The photochemistry observed for this series of molecules allows us to extract qualitative rules to rationalize the light-induced behavior of more complex B-N-containing molecules.

Machine learning-based correction for spin-orbit coupling effects in NMR chemical shift calculations

As one of the most powerful analytical methods for molecular and solid-state structure elucidation, NMR spectroscopy is an integral part of chemical laboratories associated with a great research interest in its computational simulation. Particularly when heavy atoms are present, a relativistic treatment is essential in the calculations as these influence also the nearby light atoms. In this work, we present a Δ-machine learning method that approximates the contribution to 13C and 1H NMR chemical shifts that stems from spin-orbit (SO) coupling effects. It is built on computed reference data at the spin-orbit zeroth-order regular approximation (ZORA) DFT level for a set of 6388 structures with 38 740 13C and 64 436 1H NMR chemical shifts. The scope of the methods covers the 17 most important heavy p-block elements that exhibit heavy atom on the light atom (HALA) effects to covalently bound carbon or hydrogen atoms. Evaluated on the test data set, the approach is able to recover roughly 85% of the SO contribution for 13C and 70% for 1H from a scalar-relativistic PBE0/ZORA-def2-TZVP calculation at virtually no extra computational costs. Moreover, the method is transferable to other baseline DFT methods even without retraining the model and performs well for realistic organotin and -lead compounds. Finally, we show that using a combination of the new approach with our previous Δ-ML method for correlation contributions to NMR chemical shifts, the mean absolute NMR shift deviations from non-relativistic DFT calculations to experimental values can be halved.

Generalized perturbative singles corrections to the random phase approximation method: Impact on noncovalent interaction energies of closed- and open-shell dimers

The post-Kohn-Sham (KS) random phase approximation (RPA) method may provide a poor description of interaction energies of weakly bonded molecules due to inherent density errors in approximate KS functionals. To overcome these errors, we develop a generalized formalism to incorporate perturbative singles (pS) corrections to the RPA method using orbital rotations as a perturbation parameter. The pS schemes differ in the choice of orbital-rotation gradient and Hessian. We propose a pS scheme termed RPA singles (RPAS)[Hartree-Fock (HF)] that uses the RPA orbital-rotation gradient and time-dependent HF Hessian. This correction reduces the errors in noncovalent interaction energies of closed- and open-shell dimers. For the open-shell dimers, the RPAS(HF) method leads to a consistent error reduction by 50% or more compared to the RPA method for the cases of hydrogen-bonding, metal-solvent, carbene-solvent, and dispersion interactions. We also find that the pS corrections are more important in error reduction compared to higher-order exchange corrections to the RPA method. Overall, for open shells, the RPAS(HF)-corrected RPA method provides chemical accuracy for noncovalent interactions and is more reliable than other perturbative schemes and dispersion-corrected density functional approximations, highlighting its importance as a reliable beyond-RPA correction.

Computer-Aided Molecular Design of Ionic Liquids as Advanced Process Media: A Review from Fundamentals to Applications

The unique physicochemical properties, flexible structural tunability, and giant chemical space of ionic liquids (ILs) provide them a great opportunity to match different target properties to work as advanced process media. The crux of the matter is how to efficiently and reliably tailor suitable ILs toward a specific application. In this regard, the computer-aided molecular design (CAMD) approach has been widely adapted to cover this family of high-profile chemicals, that is, to perform computer-aided IL design (CAILD). This review discusses the past developments that have contributed to the state-of-the-art of CAILD and provides a perspective about how future works could pursue the acceleration of the practical application of ILs. In a broad context of CAILD, key aspects related to the forward structure-property modeling and reverse molecular design of ILs are overviewed. For the former forward task, diverse IL molecular representations, modeling algorithms, as well as representative models on physical properties, thermodynamic properties, among others of ILs are introduced. For the latter reverse task, representative works formulating different molecular design scenarios are summarized. Beyond the substantial progress made, some future perspectives to move CAILD a step forward are finally provided.

Hexaguanidino-Triptycenes and Triphenylenes: Electronic Coupling in Molecules Containing Three Redox-Active o-Diguanidinobenzene Units Connected either Directly or Interacting Through Homoconjugation

Novel redox-active hexaguanidine molecules with multiple redox states were synthesized by connecting three o-diguanidinobenzene units. In 2,3,6,7,14,15-hexaguanidino-triptycenes, the three redox-active o-diguanidinobenzene units are connected through C-C bonds to the sp3 -hybridized bridgehead C atoms, and in 2,3,6,7,10,11-hexaguanidino-triphenylenes they are directly connected. The connectivity difference leads to different electronic coupling between the three redox-active o-diguanidinobenzene units, with homoconjugation being present in the triptycene, but not in the triphenylene compounds. Motivated by the appearance of an intense low-energy electronic transition, we especially analysed the effect of homoconjugation on the electronic structure and charge delocalization in the dicationic redox state of the triptycene derivatives. Then, several trinuclear high-spin cobalt (and copper) complexes were synthesized with the triphenylene and triptycene ligands, and the magnetic coupling and redox properties analysed. By choice of the coligands (hexafluoroacetylacetonate, trifluoroacetylacetonate and acetylacetonate), oxidation could be switched between metal- and ligand-centered redox events, leading to drastic changes in the magnetic or optical properties, especially as a consequence of homoconjugation in the triptycene derivatives.

Merocyanines form bacteriorhodopsins with strongly bathochromic absorption maxima

There is a need to shift the absorbance of biomolecules to the optical transparency window of tissue for applications in optogenetics and photo-pharmacology. There are a few strategies to achieve the so-called red shift of the absorption maxima. Herein, a series of 11 merocyanine dyes were synthesized and employed as chromophores in place of retinal in bacteriorhodopsin (bR) to achieve a bathochromic shift of the absorption maxima relative to bR's [Formula: see text] of 568 nm. Assembly with the apoprotein bacterioopsin (bO) led to stable, covalently bound chromoproteins with strongly bathochromic absorbance bands, except for three compounds. Maximal red shifts were observed for molecules 9, 2, and 8 in bR where the [Formula: see text] was 766, 755, and 736 nm, respectively. While these three merocyanines have different end groups, they share a similar structural feature, namely, a methyl group which is located at the retinal equivalent position 13 of the polyene chain. The absorption and fluorescence data are also presented for the retinal derivatives in their aldehyde, Schiff base (SB), and protonated SB (PSB) forms in solution. According to their hemicyanine character, the PSBs and their analogue bRs exhibited fluorescence quantum yields (Φf) several orders of magnitude greater than native bR (Φf 0.02 to 0.18 versus 1.5 × 10-5 in bR) while also exhibiting much smaller Stokes shifts than bR (400 to 1000 cm-1 versus 4030 cm-1 in bR). The experimental results are complemented by quantum chemical calculations where excellent agreement between the experimental [Formula: see text] and the calculated [Formula: see text] was achieved with the second-order algebraic-diagrammatic construction [ADC(2)] method. In addition, quantum mechanics/molecular mechanics (QM/MM) calculations were employed to shed light on the origin of the bathochromic shift of merocyanine 2 in bR compared with native bR.

On the description of conical intersections between excited electronic states with LR-TDDFT and ADC(2)

Conical intersections constitute the conceptual bedrock of our working understanding of ultrafast, nonadiabatic processes within photochemistry (and photophysics). Accurate calculation of potential energy surfaces within the vicinity of conical intersections, however, still poses a serious challenge to many popular electronic structure methods. Multiple works have reported on the deficiency of methods like linear-response time-dependent density functional theory within the adiabatic approximation (AA LR-TDDFT) or algebraic diagrammatic construction to second-order [ADC(2)]-approaches often used in excited-state molecular dynamics simulations-to describe conical intersections between the ground and excited electronic states. In the present study, we focus our attention on conical intersections between excited electronic states and probe the ability of AA LR-TDDFT and ADC(2) to describe their topology and topography, using protonated formaldimine and pyrazine as two exemplar molecules. We also take the opportunity to revisit the performance of these methods in describing conical intersections involving the ground electronic state in protonated formaldimine-highlighting in particular how the intersection ring exhibited by AA LR-TDDFT can be perceived either as a (near-to-linear) seam of intersection or two interpenetrating cones, depending on the magnitude of molecular distortions within the branching space.

Protonation of Homocitrate and the E1 State of Fe-Nitrogenase Studied by QM/MM Calculations

Nitrogenase is the only enzyme that can cleave the strong triple bond in N2, making nitrogen available for biological life. There are three isozymes of nitrogenase, differing in the composition of the active site, viz., Mo, V, and Fe-nitrogenase. Recently, the first crystal structure of Fe-nitrogenase was presented. We have performed the first combined quantum mechanical and molecular mechanical (QM/MM) study of Fe-nitrogenase. We show with QM/MM and quantum-refinement calculations that the homocitrate ligand is most likely protonated on the alcohol oxygen in the resting E0 state. The most stable broken-symmetry (BS) states are the same as for Mo-nitrogenase, i.e., the three Noodleman BS7-type states (with a surplus of β spin on the eighth Fe ion), which maximize the number of nearby antiferromagnetically coupled Fe-Fe pairs. For the E1 state, we find that protonation of the S2B μ2 belt sulfide ion is most favorable, 14-117 kJ/mol more stable than structures with a Fe-bound hydride ion (the best has a hydride ion on the Fe2 ion) calculated with four different density-functional theory methods. This is similar to what was found for Mo-nitrogenase, but it does not explain the recent EPR observation that the E1 state of Fe-nitrogenase should contain a photolyzable hydride ion. For the E1 state, many BS states are close in energy, and the preferred BS state differs depending on the position of the extra proton and which density functional is used.

NHC-Stabilized Dialanes(4) of Al2 Mes4

The synthesis and characterization of novel N-heterocyclic carbene (NHC) stabilized dialanes Al2 Mes4 as well as first investigations concerning the reactivity of these compounds are reported. The synthesis of these compounds proceeds via the mesityl-substituted alanes (NHC)⋅AlHMes2 (NHC=IMeMe {=1,3,4,5-tetramethyl-imidazolin-2-ylidene}, IiPrMe {=1,3-di-iso-propyl-4,5-dimethylimidazolin-2-ylidene}) and iodo-alanes (NHC)⋅AlIMes2 (NHC=IMeMe , IiPrMe ). Metallic reduction of (NHC)⋅AlIMes2 afforded the new NHC-stabilized dialanes (NHC)2 ⋅Al2 Mes4 (NHC=IMeMe , IiPrMe ). The NHC-ligated dialanes are thermally robust and storable synthons for the dialane Al2 Mes4 . First reactivity studies on (IMeMe )2 ⋅Al2 Mes4 towards small molecules confirm this, as this compound shows controlled and selective reactions with several substrates. Reaction with CuCl leads to oxidation of the dialane and formation of (IMeMe )⋅AlClMes2 , reactions with pyridine N-oxide and t Bu-N=C=S, respectively, gave the chalcogenide-bridged dimers {(IMeMe )⋅AlMes2 }2 -μ-E (E=O, S), and reaction with acetylene afforded the dimetallaacetylide {(IMeMe )⋅AlMes2 }2 -μ-(C≡C).

ConfSolv: Prediction of Solute Conformer-Free Energies across a Range of Solvents

Predicting Gibbs free energy of solution is key to understanding the solvent effects on thermodynamics and reaction rates for kinetic modeling. Accurately computing solution free energies requires the enumeration and evaluation of relevant solute conformers in solution. However, even after generation of relevant conformers, determining their free energy of solution requires an expensive workflow consisting of several ab initio computational chemistry calculations. To help address this challenge, we generate a large data set of solution free energies for nearly 44,000 solutes with almost 9 million conformers calculated in 41 different solvents using density functional theory and COSMO-RS and quantify the impact of solute conformers on the solution free energy. We then train a message passing neural network to predict the relative solution free energies of a set of solute conformers, enabling the identification of a small subset of thermodynamically relevant conformers. The model offers substantial computational time savings with predictions usually substantially within 1 kcal/mol of the free energy of the solution calculated by using computational chemical methods.

A Chiral Titanocene Complex as Regiodivergent Photoredox Catalyst: Synthetic Scope and Mechanism of Catalyst Generation

We describe a combined synthetic, spectroscopic, and computational study of a chiral titanocene complex as a regiodivergent photoredox catalyst (PRC). With Kagan's complex catCl2 either monoprotected 1,3-diols or 1,4-diols can be obtained in high selectivity from a common epoxide substrate in a regiodivergent epoxide opening depending on which enantiomer of the catalyst is employed. Due to the catalyst-controlled regioselectivity of ring opening and the broader substrate scope, the PRC with catCl2 is also a highly attractive branching point for diversity-oriented synthesis. The photochemical processes of cat(NCS)2, a suitable model for catCl2, were probed by time-correlated single-photon counting. The photoexcited complex displays a thermally activated delayed fluorescence as a result of a singlet-triplet equilibration, S1 ⇄ T1, via intersystem crossing and recrossing. Its triplet state is quenched by electron transfer to the T1 state. Computational and cyclic voltammetry studies highlight the importance of our sulfonamide additive. By bonding to sulfonamide additives, chloride abstraction from [catCl2]- is facilitated, and catalyst deactivation by coordination of the sulfonamide group is circumvented.

TURBOMOLE: Today and Tomorrow

TURBOMOLE is a highly optimized software suite for large-scale quantum-chemical and materials science simulations of molecules, clusters, extended systems, and periodic solids. TURBOMOLE uses Gaussian basis sets and has been designed with robust and fast quantum-chemical applications in mind, ranging from homogeneous and heterogeneous catalysis to inorganic and organic chemistry and various types of spectroscopy, light-matter interactions, and biochemistry. This Perspective briefly surveys TURBOMOLE's functionality and highlights recent developments that have taken place between 2020 and 2023, comprising new electronic structure methods for molecules and solids, previously unavailable molecular properties, embedding, and molecular dynamics approaches. Select features under development are reviewed to illustrate the continuous growth of the program suite, including nuclear electronic orbital methods, Hartree-Fock-based adiabatic connection models, simplified time-dependent density functional theory, relativistic effects and magnetic properties, and multiscale modeling of optical properties.

Titanium-Catalyzed Polymerization of a Lewis Base-Stabilized Phosphinoborane

The reaction of the Lewis base-stabilized phosphinoborane monomer tBuHPBH2 NMe3 (2 a) with catalytic amounts of bis(η5 :η1 -adamantylidenepentafulvene)titanium (1) provides a convenient new route to the polyphosphinoborane [tBuPH-BH2 ]n (3 a). This method offers access to high molar mass materials under mild conditions and with short reaction times (20 °C, 1 h in toluene). It represents an unprecedented example of a transition metal-mediated polymerization of a Lewis base-stabilized Group 13/15 compound. Preliminary studies of the substrate scope and a potential mechanism are reported.

Novel hybrid QSPR-GPR approach for modeling of carbon dioxide capture using deep eutectic solvents

In recent years, deep eutectic solvents (DESs) have garnered considerable attention for their potential in carbon capture and utilization processes. Predicting the carbon dioxide (CO2) solubility in DES is crucial for optimizing these solvent systems and advancing their application in sustainable technologies. In this study, we presented an evolving hybrid Quantitative Structure-Property Relationship and Gaussian Process Regression (QSPR-GPR) model that enables accurate predictions of CO2 solubility in various DESs. The QSPR-GPR model combined the strengths of both approaches, leveraging molecular descriptors and structural features of DES components to establish a robust and adaptable predictive framework. Through a systematic evolution process, we iteratively refined the model, enhancing its performance and generalization capacity. By incorporating experimental CO2 solubility data in varied DES compositions and temperatures, we trained the model to capture the intricate solubility behaviour precisely. The analytical capability of the evolving hybrid model was validated against an extensive dataset of experimental CO2 solubility values, demonstrating its superiority over individual QSPR and GPR models. The model achieves high accuracy, capturing the complex interactions between CO2 and DES components under varying thermodynamic conditions. The versatility of the evolving hybrid model was highlighted by its ability to accommodate new experimental data and adapt to different DES compositions and temperatures. The proposed QSPR-GPR model presented a powerful tool for predicting CO2 solubility in DES, providing valuable insights for designing and optimizing solvent systems in carbon capture technologies. The model's remarkable performance enhances our understanding of CO2 solubility mechanisms and contributes to sustainable solutions for mitigating greenhouse gas emissions. As research in DESs progresses, the evolving hybrid QSPR-GPR model offers a versatile and accurate means for predicting CO2 solubility, supporting advancements in carbon capture and utilization processes towards a greener and more sustainable future.

Simulation workflows to predict the circular dichroism and circularly polarized luminescence of chiral materials

Chiral materials are attracting considerable interest in various fields in view of their unique properties and optical activity. Indeed, the peculiar features of chiral materials to absorb and emit circularly polarized light enable their use in an extensive range of applications. Motivated by the interest in boosting the development of chiral materials characterized by enhanced chiroptical properties such as circular dichroism (CD) and circular polarized luminescence (CPL), we herein illustrate in this tutorial how theoretical simulations can be used for the predictions and interpretations of chiroptical data and for the identification of chiral geometries. We are focusing on computational frameworks that can be used to investigate the theoretical aspects of chiral materials' photophysical and conformational characteristics. We will then illustrate ab initio methods based on density functional theory (DFT) and its time-dependent extension (TD-DFT) to simulate CD and CPL signals, and we will exemplify a variety of enhanced sampling techniques useful for an adequate sampling of the configurational space for chiral systems.

Regularized and Opposite Spin-Scaled Functionals from Møller-Plesset Adiabatic Connection─Higher Accuracy at Lower Cost

Noncovalent interactions (NCIs) play a crucial role in biology, chemistry, material science, and everything in between. To improve pure quantum-chemical simulations of NCIs, we propose a methodology for constructing approximate correlation energies by combining an interpolation along the Møller-Plesset adiabatic connection (MP AC) with a regularization and spin-scaling strategy applied to MP2 correlation energies. This combination yields cosκos-SPL2, which exhibits superior accuracy for NCIs compared to any of the individual strategies. With the N4 formal scaling, cosκos-SPL2 is competitive or often outperforms more expensive dispersion-corrected double hybrids for NCIs. The accuracy of cosκos-SPL2 particularly shines for anionic halogen bonded complexes, where it surpasses standard dispersion-corrected DFT by a factor of 3 to 5.

N-Heterocyclic carbene and cyclic (alkyl)(amino)carbene ligated half-sandwich complexes of chromium(II) and chromium(I)

The synthesis and characterization of a series of Cr(II) N-Heterocyclic Carbene (NHC) complexes of the type [{Cr(NHC)Cl(μ-Cl)}2] and [(Cyp)Cr(NHC)X] (Cyp = η5-C5H5, cyclopentadienyl; η5-C5Me5, pentamethylcyclopentadienyl; X = Cl, η3-C3H5; NHC = IMeMe, IiPrMe, IMes, IDipp) as well as the cyclic (alkyl)(amino)carbene cAACMe ligated complexes [(η5-C5H5)Cr(cAACMe)X] (X = Cl, NPh2), [(η5-C9H7)Cr(cAACMe)Cl] (C9H7 = Ind, indenyl) and [(η5-C13H9)Cr(cAACMe)Cl] (C13H9 = Fl, fluorenyl) are reported. The reduction of [(η5-C5Me5)Cr(IMeMe)Cl] with KC8 in the presence of CO afforded the NHC ligated Cr(I) metallo-radical [(η5-C5Me5)Cr(IMeMe)(CO)2]. Quantum chemical calculations performed on [(η5-C5Me5)Cr(IMeMe)(CO)2] confirm for this complex a predominantly chromium centered radical.

Disequilibrating azobenzenes by visible-light sensitization under confinement

Photoisomerization of azobenzenes from their stable E isomer to the metastable Z state is the basis of numerous applications of these molecules. However, this reaction typically requires ultraviolet light, which limits applicability. In this study, we introduce disequilibration by sensitization under confinement (DESC), a supramolecular approach to induce the E-to-Z isomerization by using light of a desired color, including red. DESC relies on a combination of a macrocyclic host and a photosensitizer, which act together to selectively bind and sensitize E-azobenzenes for isomerization. The Z isomer lacks strong affinity for and is expelled from the host, which can then convert additional E-azobenzenes to the Z state. In this way, the host-photosensitizer complex converts photon energy into chemical energy in the form of out-of-equilibrium photostationary states, including ones that cannot be accessed through direct photoexcitation.

Luminescence and Stability of 1,4,5-Triaryl-1,2,3-Triazoles

The fluorescence, phosphorescence, and photochemical properties of di- and triaryl-substituted-1,2,3-triazoles are reported in this work. The ease of synthesis of regioisomers of substituted triazoles enables a systematic study on the correlation between regiochemistry and excited state properties, which include the solvent dependence of fluorescence, energy gap between singlet and triplet emitters, and propensity to photon-triggered transformations. The triazoles that carry electron (e)-donor and e-acceptor aryl substituents show high fluorescence quantum yields in weakly polar solvents and exhibit solvent-dependent fluorescence. The luminescence properties of these compounds in glass matrices at 77 K are characterized. The thermal and photo-stability, two parameters that are crucial to their potential utilities in optical devices, of these compounds are determined. The position of the e-donor substituent has a significant impact on the fluorescence emission energy and solvent sensitivity, singlet-triplet energy gap, and photochemical reactivity and stability. The experimental observations on the structural correlation with the photophysical and photochemical properties are explained by quantum chemical calculations. This study provides a rationale on the placement of substituent on a donor-acceptor type fluorophore to maneuver a range of photo-related properties.

Deciphering the Influence of Ground-State Distributions on the Calculation of Photolysis Observables

Nonadiabatic molecular dynamics offers a powerful tool for studying the photochemistry of molecular systems. Key to any nonadiabatic molecular dynamics simulation is the definition of its initial conditions (ICs), ideally representing the initial molecular quantum state of the system of interest. In this work, we provide a detailed analysis of how ICs may influence the calculation of experimental observables by focusing on the photochemistry of methylhydroperoxide (MHP), the simplest and most abundant organic peroxide in our atmosphere. We investigate the outcome of trajectory surface hopping simulations for distinct sets of ICs sampled from different approximate quantum distributions, namely harmonic Wigner functions and ab initio molecular dynamics using a quantum thermostat (QT). Calculating photoabsorption cross-sections, quantum yields, and translational kinetic energy maps from the results of these simulations reveals the significant effect of the ICs, in particular when low-frequency (∼ a few hundred cm-1) normal modes are connected to the photophysics of the molecule. Overall, our results indicate that sampling ICs from ab initio molecular dynamics using a QT is preferable for flexible molecules with photoactive low-frequency modes. From a photochemical perspective, our nonadiabatic dynamics simulations offer an explanation for a low-energy tail observed at high excitation energy in the translational kinetic energy map of MHP.

NEXMD v2.0 Software Package for Nonadiabatic Excited State Molecular Dynamics Simulations

We present NEXMD version 2.0, the second release of the NEXMD (Nonadiabatic EXcited-state Molecular Dynamics) software package. Across a variety of new features, NEXMD v2.0 incorporates new implementations of two hybrid quantum-classical dynamics methods, namely, Ehrenfest dynamics (EHR) and the Ab-Initio Multiple Cloning sampling technique for Multiconfigurational Ehrenfest quantum dynamics (MCE-AIMC or simply AIMC), which are alternative options to the previously implemented trajectory surface hopping (TSH) method. To illustrate these methodologies, we outline a direct comparison of these three hybrid quantum-classical dynamics methods as implemented in the same NEXMD framework, discussing their weaknesses and strengths, using the modeled photodynamics of a polyphenylene ethylene dendrimer building block as a representative example. We also describe the expanded normal-mode analysis and constraints for both the ground and excited states, newly implemented in the NEXMD v2.0 framework, which allow for a deeper analysis of the main vibrational motions involved in vibronic dynamics. Overall, NEXMD v2.0 expands the range of applications of NEXMD to a larger variety of multichromophore organic molecules and photophysical processes involving quantum coherences and persistent couplings between electronic excited states and nuclear velocity.

Assessment of DFT functionals for a minimal nitrogenase [Fe(SH)4H]- model employing state-of-the-art ab initio methods

We have designed a [Fe(SH)4H]- model with the fifth proton binding either to Fe or S. We show that the energy difference between these two isomers (∆E) is hard to estimate with quantum-mechanical (QM) methods. For example, different density functional theory (DFT) methods give ∆E estimates that vary by almost 140 kJ/mol, mainly depending on the amount of exact Hartree-Fock included (0%-54%). The model is so small that it can be treated by many high-level QM methods, including coupled-cluster (CC) and multiconfigurational perturbation theory approaches. With extrapolated CC series (up to fully connected coupled-cluster calculations with singles, doubles, and triples) and semistochastic heat-bath configuration interaction methods, we obtain results that seem to be converged to full configuration interaction results within 5 kJ/mol. Our best result for ∆E is 101 kJ/mol. With this reference, we show that M06 and B3LYP-D3 give the best results among 35 DFT methods tested for this system. Brueckner doubles coupled cluster with perturbaitve triples seems to be the most accurate coupled-cluster approach with approximate triples. CCSD(T) with Kohn-Sham orbitals gives results within 4-11 kJ/mol of the extrapolated CC results, depending on the DFT method. Single-reference CC calculations seem to be reasonably accurate (giving an error of ∼5 kJ/mol compared to multireference methods), even if the D1 diagnostic is quite high (0.25) for one of the two isomers.

Two local minima for structures of [4Fe-4S] clusters obtained with density functional theory methods

[4Fe-4S] clusters are essential cofactors in many proteins involved in biological redox-active processes. Density functional theory (DFT) methods are widely used to study these clusters. Previous investigations have indicated that there exist two local minima for these clusters in proteins. We perform a detailed study of these minima in five proteins and two oxidation states, using combined quantum mechanical and molecular mechanical (QM/MM) methods. We show that one local minimum (L state) has longer Fe-Fe distances than the other (S state), and that the L state is more stable for all cases studied. We also show that some DFT methods may only obtain the L state, while others may obtain both states. Our work provides new insights into the structural diversity and stability of [4Fe-4S] clusters in proteins, and highlights the importance of reliable DFT methods and geometry optimization. We recommend r2SCAN for optimizing [4Fe-4S] clusters in proteins, which gives the most accurate structures for the five proteins studied.

Computation of CCSD(T)-Quality NMR Chemical Shifts via Δ-Machine Learning from DFT

NMR spectroscopy undoubtedly plays a central role in determining molecular structures across different chemical disciplines, and the accurate computational prediction of NMR parameters is highly desirable. In this work, a new Δ-machine learning approach is presented to correct DFT-computed NMR chemical shifts using input features from the calculation and in addition highly accurate reference data at the CCSD(T)/pcSseg-2 level of theory with a basis set extrapolation scheme. The model is trained on a data set containing 1000 optimized and geometrically distorted structures of small organic molecules comprising most elements of the first three periods and containing data for 7090 ¹H and 4230 ¹³C NMR chemical shifts. Applied to the PBE0/pcSseg-2 method, the mean absolute deviation (MAD) on the internal NMR shift test set is reduced by 81% for ¹H and 92% for ¹³C at virtually no additional computational cost. For 12 different DFT functional and basis set combinations, the MAD of the ML-corrected NMR shifts ranges from 0.021 to 0.039 ppm (¹H) and from 0.38 to 1.07 ppm (¹³C). Importantly, the new method consistently outperforms the simple and widely used linear regression correction technique. This behavior is reproduced on three different external benchmark sets, confirming the generality and robustness of the correction scheme, which can easily be applied in DFT-based spectral simulations.

Latent Carbene in Diaminomethylation of Benzenes: Mechanism and Practical Application

Silylformamidine 1 exists in equilibrium with its carbenic form 1' due to an easy migration of the silyl group. The reaction of 1 with variously substituted fluorobenzenes proceeds as an insertion of the nucleophilic carbene 1' into the most acidic C-H bond upon mixing the reagents and does not require any catalyst. According to DFT calculations, the classical interpretation of the insertion reaction proceeding via a three-membered transition state structure requires high activation energy. Instead, low activation barriers are predicted for a transfer of the most acidic proton in the aromatic substrate to the carbene carbon. As the next step, a barrierless rearrangement of the formed ion pair toward the product completes the process. The reactivity of substituted benzenes in the reaction with silylformamidine can be roughly assessed by calculated pK_a (DMSO) values for the C-H hydrogens. Benzene derivatives having pK_a approx. less than 31 can undergo C-H insertion. The reaction provides aminals as the first products, which can easily be transformed into the corresponding aldehydes via acidic hydrolysis. As silylformamidine 1 is tolerant to many functional groups, the reaction can be applied to numerous benzene derivatives, making it a reliable strategy for application in organic synthesis.

Catalytic Reaction Mechanism of Glyoxalase II: A Quantum Mechanics/Molecular Mechanics Study

Methylglyoxal (MG) is a reactive and toxic compound produced in carbohydrate, lipid, and amino acid metabolism. The glyoxalase system is the main detoxifying route for MG and consists of two enzymes, glyoxalase I (GlxI) and glyoxalase II (GlxII). GlxI catalyzes the formation of S-d-lactoylglutathione from hemithioacetal, and GlxII converts this intermediate to d-lactate. A relationship between the glyoxalase system and some diseases like diabetes has been shown, and inhibiting enzymes of this system may be an effective means of controlling certain diseases. A detailed understanding of the reaction mechanism of an enzyme is essential to the rational design of competitive inhibitors. In this work, we use quantum mechanics/molecular mechanics (QM/MM) calculations and energy refinement utilizing the big-QM and QM/MM thermodynamic cycle perturbation methods to propose a mechanism for the GlxII reaction that starts with a nucleophilic attack of the bridging OH^- group on the substrate. The coordination of the substrate to the Zn ions places its electrophilic center close to the hydroxide group, enabling the reaction to proceed. Our estimated reaction energies are in excellent agreement with experimental data, thus demonstrating the reliability of our approach and the proposed mechanism. Additionally, we examined alternative protonation states of Asp-29, Asp-58, Asp-134, and the bridging hydroxide ion in the catalytic process. However, these give less favorable reactions, a poorer reproduction of the crystal structure geometry of the active site, and higher root-mean-squared deviations of the active site residues in molecular dynamics simulations.

Evaluating the Transition State Stabilization/Destabilization Effects of the Electric Fields from Scaffold Residues by a QM/MM Approach

The protein scaffolds of enzymes not only provide structural support for the catalytic center but also exert preorganized electric fields for electrostatic catalysis. In recent years, uniform oriented external electric fields (OEEFs) have been widely applied to enzymatic reactions to mimic the electrostatic effects of the environment. However, the electric fields exerted by individual residues in proteins may be quite heterogeneous across the active site, with varying directions and strengths at different positions of the active site. Here, we propose a QM/MM-based approach to evaluate the effects of the electric fields exerted by individual residues in the protein scaffold. In particular, the heterogeneity of the residue electric fields and the effect of the native protein environment can be properly accounted for by this QM/MM approach. A case study of the O-O heterolysis reaction in the catalytic cycle of TyrH shows that (1) for scaffold residues that are relatively far from the active site, the heterogeneity of the residue electric field in the active site is not very significant and the electrostatic stabilization/destabilization due to each residue can be well approximated with the interaction energy between a uniform electric field and the QM region dipole; (2) for scaffold residues near the active site, the residue electric fields can be highly heterogeneous along the breaking O-O bond. In such a case, approximating the residue electric fields as uniform fields may misrepresent the overall electrostatic effect of the residue. The present QM/MM approach can be applied to evaluate the residues' electrostatic impact on enzymatic reactions, which also can be useful in computational optimization of electric fields to boost the enzyme catalysis.

Unveiling the tartrazine binding mode with ds-DNA by UV-visible spectroscopy, electrochemical, and QM/MM methods

Here, we studied the interaction between the food colorant tartrazine (TZ) and double stranded DNA (dsDNA), using spectroscopic, electrochemical, and computational methods such as QM/MM combined with TD-DFT. Despite the UV-vis spectroscopy is widely used to study the interaction between molecules, for the case of TZ there are discrepancies in the analyses presented in the literature available, presenting both hyperchromic and hypochromic effects and consequently different rationalizations for their results. Herein we propose the combination of UV-vis experiments with the design of high-level computational models capable of reproducing the experimental behavior to finally define the proper binding mode at the molecular scale together with the rationalization of the experimental optical response due to the complex formation. To complement the UV-vis experiments, we propose the use of electrochemical measurements, to support the results obtained through UV-vis spectroscopy, as it has been successfully used for the determination of interaction modes between small molecules and biomolecules in any condition. Our UV-vis spectroscopy experiments showed only a hypochromic effect of the absorption spectra of TZ after interaction with DNA, indicative of TZ being deeply buried in the DNA structure. The effect of ionic strength in the experimental procedures led to the dissociation of TZ, thus indicating that the interaction mode was groove binding. On the other hand, the electrochemical studies showed an irreversible reduction peak of TZ, which after the interaction with DNA exhibited a positive shift in potential that can be attributed to groove binding. The binding constant for TZ-DNA was calculated as 4.45x10⁴M^-1 (UV-vis) and 5.75x10⁴M^-1 (electrochemistry), in line with other groove binder azo dyes. Finally, through the QM/MM calculations we found that the minor-groove binding mode interacting in zones rich in adenine and thymine was the model best suited to reproduce the experimental UV-vis response.

Histidine oxidation in lytic polysaccharide monooxygenase

The lytic polysaccharide monooxygenases (LPMOs) comprise a super-family of copper enzymes that boost the depolymerisation of polysaccharides by oxidatively disrupting the glycosidic bonds connecting the sugar units. Industrial use of LPMOs for cellulose depolymerisation has already begun but is still far from reaching its full potential. One issue is that the LPMOs self-oxidise and thereby deactivate. The mechanism of this self-oxidation is unknown, but histidine residues coordinating to the copper atom are the most susceptible. An unusual methyl modification of the NE2 atom in one of the coordinating histidine residues has been proposed to have a protective role. Furthermore, substrate binding is also known to reduce oxidative damage. We here for the first time investigate the mechanism of histidine oxidation with combined quantum and molecular mechanical (QM/MM) calculations, with outset in intermediates previously shown to form from a reaction with peroxide and a reduced LPMO. We show that an intermediate with a [Cu-O]⁺ moiety is sufficiently potent to oxidise the nearest C-H bond on both histidine residues, but methylation of the NE2 atom of His-1 increases the reaction barrier of this reaction. The substrate further increases the activation barrier. We also investigate a [Cu-OH]²⁺ intermediate with a deprotonated tyrosine radical. This intermediate was previously proposed to have a protective role, and we also find it to have higher barriers than the corresponding a [Cu-O]⁺ intermediate.

Reducing Exact Two-Component Theory for NMR Couplings to a One-Component Approach: Efficiency and Accuracy

The self-consistent and complex spin-orbit exact two-component (X2C) formalism for NMR spin-spin coupling constants [ J. Chem. Theory Comput. 17, 2021, 3874-3994] is reduced to a scalar one-component ansatz. This way, the first-order response term can be partitioned into the Fermi-contact (FC) and spin-dipole (SD) interactions as well as the paramagnetic spin-orbit (PSO) contribution. The FC+SD terms are real and symmetric, while the PSO term is purely imaginary and antisymmetric. The relativistic one-component approach is combined with a modern density functional treatment up to local hybrid functionals including the response of the current density. Computational demands are reduced by factors of 8-24 as shown for a large tin compound consisting of 137 atoms. Limitations of the current ansatz are critically assessed for Sn, Pb, Pd, and Pt compounds, i.e. the one-component treatment is not sufficient for tin compounds featuring a few heavy halogen atoms.

Insights into molecular cluster materials with adamantane-like core structures by considering dimer interactions

A class of adamantane-like molecular materials attracts attention because they exhibit an extreme non-linear optical response and emit a broad white-light spectrum after illumination with a continuous-wave infrared laser source. According to recent studies, not only the nature of the cluster molecules, but also the macroscopic structure of the materials determines their non-linear optical properties. Here we present a systematic study of cluster dimers of the compounds AdR₄ and [(RT)₄ S₆ ] (T = Si, Ge, Sn) with R = methyl, phenyl or 1-naphthyl to gain fundamental knowledge about the interactions in the materials. For all compounds, a similar type of dimer structures with a staggered arrangement of substituents was determined as the energetically most favorable configuration. The binding energy between the dimers, determined by including London dispersion interactions, increases with the size of the core and the substituents. The cluster interactions can be classified as substituent-substituent-dominated (small cores, large substituents) or core-core-dominated (large cores, small substituents). Among various possible dimer conformers, those with small core-core distances are energetically preferred. Trimer and tetramer clusters display similar trends regarding the minimal core-core distances and binding energies. The much lower energy barrier determined for the rotation of substituents as compared to the rotation of the cluster dimers past each other indicates that the rotation of substituents more easily leads to different conformers in the material. Thus, understanding the interaction of the cluster dimers allows an initial assessment of the interactions in the materials.

Ultrafast structural changes direct the first molecular events of vision

Vision is initiated by the rhodopsin family of light-sensitive G protein-coupled receptors (GPCRs)1. A photon is absorbed by the 11-cis retinal chromophore of rhodopsin, which isomerizes within 200 femtoseconds to the all-trans conformation2, thereby initiating the cellular signal transduction processes that ultimately lead to vision. However, the intramolecular mechanism by which the photoactivated retinal induces the activation events inside rhodopsin remains experimentally unclear. Here we use ultrafast time-resolved crystallography at room temperature3 to determine how an isomerized twisted all-trans retinal stores the photon energy that is required to initiate the protein conformational changes associated with the formation of the G protein-binding signalling state. The distorted retinal at a 1-ps time delay after photoactivation has pulled away from half of its numerous interactions with its binding pocket, and the excess of the photon energy is released through an anisotropic protein breathing motion in the direction of the extracellular space. Notably, the very early structural motions in the protein side chains of rhodopsin appear in regions that are involved in later stages of the conserved class A GPCR activation mechanism. Our study sheds light on the earliest stages of vision in vertebrates and points to fundamental aspects of the molecular mechanisms of agonist-mediated GPCR activation.

Interaction of gentamicin and gentamicin-AOT with poly-(lactide-co-glycolate) in a drug delivery system - density functional theory calculations and molecular dynamics simulation

Gentamicin is used to treat brucellosis, an infectious disease caused by the Brucella species but the drug faces several issues such as low efficacy, instability, low solubility, and toxicity. It also has a very short half-life, therefore, requiring frequent dosing. Consequently, several other antibiotics are also being used for the treatment of brucellosis as a single dose as well as in combination with other antibiotics but none of these therapies are satisfactory. Nanoparticles in particular polymer-based ones utilizing polymers that are biodegradable and biocompatible for instance PLGA are a method of choice to overcome such drug delivery issues and enable potential targeted delivery. The current study focuses on the evaluation of the structural and dynamical properties of a drug-polymer system consisting of gentamicin drug and PLGA polymer nanoparticles in the water representing a targeted drug delivery system for the treatment of brucellosis. For this purpose, all-atom molecular dynamics simulations were carried out on the drug-polymer systems in the absence and presence of the surfactant bis(2-Ethylhexyl) sulfosuccinate (AOT) to determine the structural and dynamical properties as well as the effect of the surfactant on these properties. We also investigated systems in which the polymer constituents were in the form of monomeric units toward decoupling the primary interactions of the monomer units and polymer effects. The simulation results explain the nature of the interactions between the drug and the polymer as well as transport properties in terms of drug diffusion coefficients, which characterize the molecular behavior of gentamicin-polymer nanoparticles for use in brucellosis.

Cooperativity-Driven Reactivity of a Dinuclear Copper Dimethylglyoxime Complex

In this report, we present the dinuclear copper(II) dimethylglyoxime (H₂ dmg) complex [Cu₂ (H₂ dmg)(Hdmg)(dmg)]⁺ (1), which, in contrast to its mononuclear analogue [Cu(Hdmg)₂ ] (2), is subject to a cooperativity-driven hydrolysis. The combined Lewis acidity of both copper centers increases the electrophilicity of the carbon atom in the bridging μ₂ -O-N=C-group of H₂ dmg and thus, facilitates the nucleophilic attack of H₂ O. This hydrolysis yields butane-2,3-dione monoxime (3) and NH₂ OH that, depending on the solvent, is then either oxidized or reduced. In ethanol, NH₂ OH is reduced to NH₄ ⁺ , yielding acetaldehyde as the oxidation product. In contrast, in CH₃ CN, NH₂ OH is oxidized by Cu^II to form N₂ O and [Cu(CH₃ CN)₄ ]⁺ . Herein are presented the combined synthetic, theoretical, spectroscopic and spectrometric methods that indicate and establish the reaction pathway of this solvent-dependent reaction.