Fast and Accurate Quantum Mechanical Modeling of Large Molecular Systems Using Small Basis Set Hartree-Fock Methods Corrected with Atom-Centered Potentials

Insight into the activation mechanism of carbonic anhydrase(II) through 2-(2-aminoethyl)-pyridine: a promising pathway for enhanced enzymatic activity

Activation of human carbonic anhydrase II (hCA II) holds great promise for treating memory loss symptoms associated with Alzheimer's disease. Despite its importance, the activation mechanism of hCA II has been largely overlooked in favor of the well-studied inhibition mechanism. To address this unexplored realm, we use first-principles calculations to tease out the activation mechanism of hCA II using 2-(2-aminoethyl)-pyridine (2-2AEPy), a promising in vitro activator. We explored both stepwise and concerted mechanisms via both available nitrogen sites of 2-2AEPy: (i) aminoethyl group (Nα) and (ii) pyridine ring (Nβ). Our results show that a concerted mechanism via Nα holds the key to hCA II activation. The activation process of the concerted mechanism exhibits the characteristics of an exergonic reaction, wherein the transition state resembles the reactant with a notably low imaginary frequency of 452.4i cm-1 and barrier height of 5.2 kcal mol-1. Such meager transition barriers propel the activation of hCA II at in vivo temperatures. These findings initiate future research into hCA II activation mechanisms and the development of efficient activators, which may lead to promising therapeutic interventions for Alzheimer's disease.

Encapsulation in Charged Droplets Generates Distorted Host-Guest Complexes

The ability of various hydrogen-bonded resorcinarene-based capsules to bind α,ω-alkylbisDABCOnium (DnD) guests of different lengths was investigated in solution and in the gas-phase. While no host-guest interactions were detected in solution, encapsulation could be achieved in the charged droplets formed during electrospray ionisation (ESI). This included guests, which are far too long in their most stable conformation to fit inside the cavity of the capsules. A combination of three mass spectrometric techniques, namely, collision-induced dissociation, hydrogen/deuterium exchange, and ion-mobility mass spectrometry, together with computational modelling allow us to determine the binding mode of the DnD guests inside the cavity of the capsules. Significant distortions of the guest into horseshoe-like arrangements are required to optimise cation-π interactions with the host, which also adopt distorted geometries with partially open hydrogen-bonding seams when binding longer guests. Such quasi "spring-loaded" capsules can form in the charged droplets during the ESI process as there is no competition between guest encapsulation and ion pair formation with the counterions that preclude encapsulation in solution. The encapsulation complexes are sufficiently stable in the gas-phase - even when strained - because non-covalent interactions significantly strengthen in the absence of solvent.

Solvent and Substituent Size Influence on the Cyclochiral Rigidity of Aminomethylene Derivatives of Resorcin[4]arene

Resorcin[4]arenes (R[4]A) are a group of macrocyclic compounds whose peculiar feature is the presence of eight hydroxyl groups in their structure. The directional formation of intramolecular hydrogen bonds with their participation leads to the formation of a cyclochiral racemic mixture of these compounds. Their stability strongly depends on the substituent and especially the environment in which they are located. The paper discusses the cyclochiral nature of aminomethylene derivatives of R[4]A (AMD-R[4]A). Their cyclochiral rigidity in non-polar solvents has been shown. The influence of the size of the alkyl groups in the amino substituents of AMD-R[4]A on their cyclochiral nature was noted. To calculate the reaction paths for their racemization, the nudged elastic band (NEB) method was employed using the semi-empirical DFT (GFN1-xTB) approach. The calculated activation barrier energies for their racemization in chloroform, obtained through various semi-empirical quantum chemical methods (SE), Hartree-Fock (HF), and density functionals theory (DFT), show good correlation with experimental observations. Among the tested methods, the B38LYP-D4 method is highly recommended due to its fast computational speed and accuracy, which is comparable to the time-consuming double-hybrid DH-revDSD-PBEP86 approach.

Best-Practice DFT Protocols for Basic Molecular Computational Chemistry

Nowadays, many chemical investigations are supported by routine calculations of molecular structures, reaction energies, barrier heights, and spectroscopic properties. The lion's share of these quantum-chemical calculations applies density functional theory (DFT) evaluated in atomic-orbital basis sets. This work provides best-practice guidance on the numerous methodological and technical aspects of DFT calculations in three parts: Firstly, we set the stage and introduce a step-by-step decision tree to choose a computational protocol that models the experiment as closely as possible. Secondly, we present a recommendation matrix to guide the choice of functional and basis set depending on the task at hand. A particular focus is on achieving an optimal balance between accuracy, robustness, and efficiency through multi-level approaches. Finally, we discuss selected representative examples to illustrate the recommended protocols and the effect of methodological choices.

Small-Basis Set Density-Functional Theory Methods Corrected with Atom-Centered Potentials

Density functional theory (DFT) is currently the most popular method for modeling noncovalent interactions and thermochemistry. The accurate calculation of noncovalent interaction energies, reaction energies, and barrier heights requires choosing an appropriate functional and, typically, a relatively large basis set. Deficiencies of the density-functional approximation and the use of a limited basis set are the leading sources of error in the calculation of noncovalent and thermochemical properties in molecular systems. In this article, we present three new DFT methods based on the BLYP, M06-2X, and CAM-B3LYP functionals in combination with the 6-31G* basis set and corrected with atom-centered potentials (ACPs). ACPs are one-electron potentials that have the same form as effective-core potentials, except they do not replace any electrons. The ACPs developed in this work are used to generate energy corrections to the underlying DFT/basis-set method such that the errors in predicted chemical properties are minimized while maintaining the low computational cost of the parent methods. ACPs were developed for the elements H, B, C, N, O, F, Si, P, S, and Cl. The ACP parameters were determined using an extensive training set of 118655 data points, mostly of complete basis set coupled-cluster level quality. The target molecular properties for the ACP-corrected methods include noncovalent interaction energies, molecular conformational energies, reaction energies, barrier heights, and bond separation energies. The ACPs were tested first on the training set and then on a validation set of 42567 additional data points. We show that the ACP-corrected methods can predict the target molecular properties with accuracy close to complete basis set wavefunction theory methods, but at a computational cost of double-ζ DFT methods. This makes the new BLYP/6-31G*-ACP, M06-2X/6-31G*-ACP, and CAM-B3LYP/6-31G*-ACP methods uniquely suited to the calculation of noncovalent, thermochemical, and kinetic properties in large molecular systems.