Molecular acidity: A quantitative conceptual density functional theory description

Information-theoretic quantities as effective descriptors of electrophilicity and nucleophilicity in density functional theory

CONTEXT

Electrophilicity and nucleophilicity are two vastly important chemical concepts gauging the capability of atoms in molecules to accept and donate the maximal number of electrons. In our earlier studies, we proposed to simultaneously quantify them using the Kullback-Leibler divergence from the information-theoretic approach in density functional theory. However, several issues with this scheme remain to be clarified such as its general validity, predictability, and relationship with other information-theoretic quantities. In this work, we revisit the matter with bigger datasets and deeper theoretical insights. Five information-theoretic quantities including Kullback-Leibler divergence, Hirshfeld charge, Ghost-Berkowitz-Parr entropy, and second and third orders of relative Onicescu information energy are found to be reliable and robust descriptors of electrophilicity and nucleophilicity propensities. Employing these five descriptors, we design a list of new compounds and predict their electrophilicity and nucleophilicity scales. This work should markedly improve our confidence and capability in applying information-theoretic quantities to evaluate electrophilicity and nucleophilicity propensities and henceforth pave the route for more applications of these quantities from information-theoretic approach in density functional theory in the future.

METHODS

All structures were fully optimized at the M06-2X/6-311 + G(d) level of DFT functional using the Gaussian 16 package (version C01) with integration grids and tight self-consistent-field convergence. The solvent effect was taken into account by using the implicit solvent model (CPCM) in the CH2Cl2 solvent, and all 3D contour surfaces of Fukui function, local temperature, and ITA (information-theoretic approach) quantities were generated by GaussView. The Multiwfn 3.8 program was used to calculate the ITA indexes and atomic charges.

Harvesting Chemical Understanding with Machine Learning and Quantum Computers

It is tenable to argue that nobody can predict the future with certainty, yet one can learn from the past and make informed projections for the years ahead. In this Perspective, we overview the status of how theory and computation can be exploited to obtain chemical understanding from wave function theory and density functional theory, and then outlook the likely impact of machine learning (ML) and quantum computers (QC) to appreciate traditional chemical concepts in decades to come. It is maintained that the development and maturation of ML and QC methods in theoretical and computational chemistry represent two paradigm shifts about how the Schrödinger equation can be solved. New chemical understanding can be harnessed in these two new paradigms by making respective use of ML features and QC qubits. Before that happens, however, we still have hurdles to face and obstacles to overcome in both ML and QC arenas. Possible pathways to tackle these challenges are proposed. We anticipate that hierarchical modeling, in contrast to multiscale modeling, will emerge and thrive, becoming the workhorse of in silico simulations in the next few decades.

Some Recent Advances in Density-Based Reactivity Theory

Establishing a chemical reactivity theory in density functional theory (DFT) language has been our intense research interest in the past two decades, exemplified by the determination of steric effect and stereoselectivity, evaluation of electrophilicity and nucleophilicity, identification of strong and weak interactions, and formulation of cooperativity, frustration, and principle of chirality hierarchy. In this Featured Article, we first overview the four density-based frameworks in DFT to appreciate chemical understanding, including conceptual DFT, use of density associated quantities, information-theoretic approach, and orbital-free DFT, and then present a few recent advances of these frameworks as well as new applications from our studies. To that end, we will introduce the relationship among these frameworks, determining the entire spectrum of interactions with Pauli energy derivatives, performing topological analyses with information-theoretic quantities, and extending the density-based frameworks to excited states. Applications to examine physiochemical properties in external electric fields and to evaluate polarizability for proteins and crystals are discussed. A few possible directions for future development are followed, with the special emphasis on its merger with machine learning.

Directionality and additivity effects of molecular acidity and aromaticity for substituted benzoic acids under external electric fields

Our recent study [M. Li et al.Phys. Chem. Chem. Phys., 2023, 25, 2595-2605] unveiled that the impact of an external electric field on molecular acidity and aromaticity for benzoic acid is directional, which can be understood using changes in frontier orbitals and partial charges. However, it is unclear if the effect will disappear when substituting groups are present and whether new patterns of changes will show up. In this work, as a continuation of our efforts to appreciate the impact of external electric fields on physiochemical properties, we find that the directionality effect is still in place for substituted benzoic acid derivatives and that there exists the additivity effect with respect to the number of substituent groups, regardless of the direction of the applied field and the type of substituting groups. We confirm the findings using electron-donating and electron-accepting groups with the electric field applied either parallelly or perpendicularly to the carboxyl group along the benzene ring. The directionality and additivity effects uncovered from this work should enrich the body of our knowledge about the impact of external electric fields on physiochemical properties and could be applicable to other systems and properties as well.

Topological analysis of information-theoretic quantities in density functional theory

We have witnessed considerable research interest in the recent literature about the development and applications of quantities from the information-theoretic approach (ITA) in density functional theory. These ITA quantities are explicit density functionals, whose local distributions in real space are continuous and well-behaved. In this work, we further develop ITA by systematically analyzing the topological behavior of its four representative quantities, Shannon entropy, two forms of Fisher information, and relative Shannon entropy (also called information gain or Kullback-Leibler divergence). Our results from their topological analyses for 103 molecular systems provide new insights into bonding interactions and physiochemical properties, such as electrophilicity, nucleophilicity, acidity, and aromaticity. We also compare our results with those from the electron density, electron localization function, localized orbital locator, and Laplacian functions. Our results offer a new methodological approach and practical tool for applications that are especially promising for elucidating chemical bonding and reactivity propensity.

Discovery of potent Plasmodium falciparum protein kinase 6 (PfPK6) inhibitors with a type II inhibitor pharmacophore

Malaria is a devastating disease that causes significant global morbidity and mortality. The rise of drug resistance against artemisinin-based combination therapy demonstrates the necessity to develop alternative antimalarials with novel mechanisms of action. We report the discovery of Ki8751 as an inhibitor of essential kinase PfPK6. 79 derivatives were designed, synthesized and evaluated for PfPK6 inhibition and antiplasmodial activity. Using group efficiency analyses, we established the importance of key groups on the scaffold consistent with a type II inhibitor pharmacophore. We highlight modifications on the tail group that contribute to antiplasmodial activity, cumulating in the discovery of compound 67, a PfPK6 inhibitor (IC50 = 13 nM) active against the P. falciparum blood stage (EC50 = 160 nM), and compound 79, a PfPK6 inhibitor (IC50 < 5 nM) with dual-stage antiplasmodial activity against P. falciparum blood stage (EC50 = 39 nM) and against P. berghei liver stage (EC50 = 220 nM).

Impacts of external fields on aromaticity and acidity of benzoic acid: a density functional theory, conceptual density functional theory and information-theoretic approach study

The impact of external fields on the molecular structure and reactivity properties has been of considerable interest in the recent literature. Benzoic acid as one of the most widely used compounds in medicinal and materials sciences is known for its dual propensity in aromaticity and acidity. In this work, we systematically investigate the impact of a uniform external electric field on these properties. We apply density functional theory, conceptual density functional theory, and an information-theoretic approach to appreciate the change pattern of aromaticity and acidity properties in external fields with different strengths. Our results show that they possess different change patterns under external fields, which can be satisfactorily rationalized by variations in reactivity descriptors and partial charges. The surprising yet novel results from this study should enrich the body of our knowledge about the impact of external fields for different kinds of electronic properties and provide guidance and foundation for future studies of this phenomenon in other molecular systems.

Development and Applications of the Density-Based Theory of Chemical Reactivity

Density functional theory, which is well-recognized for its accuracy and efficiency, has become the workhorse for modeling the electronic structure of molecules and extended materials in recent decades. Nevertheless, establishing a density-based conceptual framework to appreciate bonding, stability, function, reactivity, and other physicochemical properties is still an unaccomplished task. In this Perspective, we at first provide an overview of the four pathways currently available in the literature to tackle the matter, including orbital-free density functional theory, conceptual density functional theory, direct use of density-associated quantities, and the information-theoretic approach. Then, we highlight several recent advances of employing these approaches to realize new understandings for chemical concepts such as covalent bonding, noncovalent interactions, cooperation, frustration, homochirality, chirality hierarchy, electrophilicity, nucleophilicity, regioselectivity, and stereoselectivity. Finally, we provide a few possibilities for the future development of this relatively uncharted territory. Opportunities are abundant, and they are all ours for the taking.

Machine learning methods for pKa prediction of small molecules: Advances and challenges

The acid-base dissociation constant (pK_a) is a fundamental property influencing many ADMET properties of small molecules. However, rapid and accurate pK_a prediction remains a great challenge. In this review, we outline the current advances in machine-learning-based QSAR models for pK_a prediction, including descriptor-based and graph-based approaches, and summarize their pros and cons. Moreover, we highlight the current challenges and future directions regarding experimental data, crucial factors influencing pK_a and in silico prediction tools. We hope that this review can provide a practical guidance for the follow-up studies.

Synthesis, in vitro cholinesterase inhibition, molecular docking, DFT and ADME studies of novel 1,3,4-oxadiazole 2-thiol derivatives

A sequence of 1,3,4-oxadiazole 2-thiol derivatives bearing various alkyl or aryl moieties was designed, synthesized, and characterized by modern spectroscopic methods to yield 17 compounds ( 6a - 6q ) which were screened for acetylcholinesterase (AChE) and butyrylcholinesterase (BChE) enzymes in search of 'lead' compounds for the treatment of Alzheimer disease (AD). The compounds 6q, 6p, 6k, 6o, and 6l showed inhibitory capability against AChE and BChE, with IC 50 values ranging from 11.730.49 to 27.360.29 µM for AChE and 21.830.39 to 39.430.44 µM for BChE, inhibiting both enzymes within a limited range. The SAR ascertained that the substitution of the aromatic moiety had a profound effect on the AChE and BChE inhibitory potential as compared to the aliphatic substitutions which were supported by the molecular docking studies. In silico ADME studies reinforced the drug-likeness of most of the synthesized molecules. These results were additionally supplemented by the molecular orbital analysis (HOMO-LUMO) and electrostatic potential maps got from DFT calculations. ESP maps expose that on all structures, there are two potential binding sites conquered by the most positive and most negative districts.

The mechanism of direct reductive amination of aldehyde and amine with formic acid catalyzed by boron trifluoride complexes: insights from a DFT study

A volcano diagram of BF3 catalytic species and their activities was proposed for the DRA of aldehyde and amine with formic acid.

Quantifications and Applications of Relative Fisher Information in Density Functional Theory

Though density functional theory is widely accepted as one of the most successful developments in theoretical chemistry in the past few decades, the knowledge of how to apply this new electronic structure theory, to help us better understand chemical processes and transformations, is still an unaccomplished task. The information-theoretic approach is emerging as a viable option for that purpose in the recent literature, providing new insights about steric effect, cooperativity, electrophilicity, nucleophilicity, stereoselectivity, homochirality, etc. In this work, based on the result from a recent paper by one of us [ J. Chem. Phys, 2019, 151, 141103], we present two quantifications of the relative Fisher information and discuss their physiochemical properties and possible applications. To that end, their analytical properties have been elucidated. They have also been applied to six categories of systems to illustrate their applicability. A better descriptor to quantify the single bond rotation barrier has been obtained. The relative Fisher information can also simultaneously determine electrophilicity and nucleophilicity, and effectively describe helical structures with different homochiral and heterochiral propensities. As integral parts of the information-theoretic approach, these newly introduced quantities will provide us with more analytical tools toward the long-term goal of crafting a chemical reactivity theory in the density-based language.

Using Bases as Initiators to Isomerize Allylic Alcohols: Insights from Density Functional Theory Studies

Allylic alcohols, as common and readily available building blocks, could be converted into many widely used carbonyl compounds through isomerization reactions. However, these processes often involve expensive transition metal (TM) complexes as the catalyst. What is the bottleneck in the mechanism when no TM is used? In this study, density functional theory (DFT) was employed to explore the mechanistic patterns of allylic alcohols catalyzed using bases, such as KOH, NaOH, LiOH, tBuOK, tBuONa, tBuOLi, 1,5,7-triazabicyclo[4.4.0]dec-5-ene, 1,3,4,6,7,8-hexahydro-1-methyl-2H-pyrimido[1,2-a]pyrimidine, and 1,8-diazabicyclo[5.4.0]undec-7-ene. Our results show that bases containing metal cations follow the metal cation-assisted (MCA) mechanism, whereas organic bases without metal cations follow the ion pair-assisted (IPA) mechanism. The catalytic efficiency of bases containing metal cations is higher than that of bases without metal cations, indicating that metal cations play an important role in the reaction. Additionally, the modulation of substituents R₁ and R₂ in the substrate reveals that electron-withdrawing groups are favorable for C-H bond cleavage, and electron-donating groups are favorable for hydrogen transfer. To better understand these patterns, we applied the DFT and information-theoretic approach (ITA) to examine the impact of bases and substrate substituents on the reactivity of allylic alcohol isomerization. This work should provide a much-needed theoretical guidance to design better non-TM catalysts for the isomerization of allylic alcohols and their derivatives.

Synthesis, characterization, reaction mechanism prediction and biological study of mono, bis and tetrakis pyrazole derivatives against Fusarium oxysporum f. sp. Albedinis with conceptual DFT and ligand-protein docking studies

Twelve heterocyclic compounds were prepared using the condensation of hydroxymethanol pyrazole derivatives with different primary aminesas example 2-aminothiazole and 1-aminobenzotriazole to have a diverse productin good yield up to 97%. Those ligands were tested against Fusarium oxysporum f. sp. Albedinis fungi (BAYOUD Disease) with IC₅₀ = 25.6-33.2 µg/ml. After experiments, theoretical investigations were done as DFT study to know the ligands molecular reactivity and the-ligandprotein- docking study to know the possible binding between the prepared ligands with two biological targets: FGB1 (Fusarium oxysporum Guanine nucleotide-binding protein beta subunitprimary amino acid sequence) and Fophy (Fusarium oxysporum phytase domain enzyme). Of all the obtained results, the experimental ones were well correlated with the theoretical with the most common thing between those compounds is (N^δ--N^δ+) which is the antifungal pharmacophore as proposed pincers for Foa inhibition. From docking studies over FGB1 and Fophy, the ligand 9 has the best binding energy of -6.4872 kcal/mol in FGB1 active site and -5.5282 kcal/mol in Fophy active site, but better correlation with Fophy than FGB1 which is followed by PLIF graph to get that Arg116, Arg120 and Lys336 are the vital amino acids of fophy protein based the study over the chosen active site.

A theoretical study of the hydroboration of α,β-unsaturated carbonyl compounds catalyzed by a metal-free complex and subsequent C–C coupling with acetonitrile

Herein, the density functional theory (DFT) method was employed to investigate the reaction mechanism of the selective hydroboration of α,β-unsaturated carbonyl compounds catalyzed by the metal-free complex 1,3,2-diazaphospholene (DAP).

The Effect of Intramolecular Hydrogen Bond Type on the Gas-Phase Deprotonation of ortho-Substituted Benzenesulfonic Acids. A Density Functional Theory Study

Structural factors have been identified that determine the gas-phase acidity of ortho-substituted benzenesulfonic acid, 2-XC6H4-SO3H, (X = -SO3H, -COOH, -NO2, -SO2F, -C≡N, -NH2, -CH3, -OCH3, -N(CH3)2, -OH). The DFT/B3LYP/cc-pVTZ method was used to perform conformational analysis and study the structural features of the molecular and deprotonated forms of these compounds. It has been shown that many of the conformers may contain anintramolecular hydrogen bond (IHB) between the sulfonic group and the substituent, and the sulfonic group can be an IHB donor or an acceptor. The Gibbs energies of gas-phase deprotonation ΔrG0298 (kJ mol-1) were calculated for all compounds. It has been set that in ortho-substituted benzenesulfonic acids, the formation of various types of IHB is possible, having a significant effect on the ΔrG0298 values of gas-phase deprotonation. If the -SO3H group is the IHB donor, then an ion without an IHB is formed upon deprotonation, and the deprotonation energy increases. If this group is an IHB acceptor, then a significant decrease in ΔrG0298 of gas-phase deprotonation is observed due to an increase in IHB strength and the A- anion additional stabilization. A proton donor ability comparative characteristic of the -SO3H group in the studied ortho-substituted benzenesulfonic acids is given, and the ΔrG0298 energies are compared with the corresponding values of ortho-substituted benzoic acids.

Unifying Conceptual Density Functional and Valence Bond Theory: The Hardness-Softness Conundrum Associated with Protonation Reactions and Uncovering Complementary Reactivity Modes

In this study, we address the long-standing issue-arising prominently from conceptual density functional theory (CDFT)-of the relative importance of electrostatic, i.e., "hard-hard", versus spin-pairing, i.e., "soft-soft", interactions in determining regiochemical preferences. We do so from a valence bond (VB) perspective and demonstrate that VB theory readily enables a clear-cut resolution of both of these contributions to the bond formation/breaking process. Our calculations indicate that appropriate local reactivity descriptors can be used to gauge the magnitude of both interactions individually, e.g., Fukui functions or HOMO/LUMO orbitals for the spin-pairing/(frontier) orbital interactions and molecular electrostatic potentials (and/or partial charges) for the electrostatic interactions. In contrast to previous reports, we find that protonation reactions cannot generally be classified as either charge- or frontier orbital-controlled; instead, our results indicate that these two bonding contributions generally interplay in more subtle patterns, only giving the impression of a clear-cut dichotomy. Finally, we demonstrate that important covalent, i.e., spin pairing, reactivity modes can be missed when only a single spin-pairing/orbital interaction descriptor is considered. This study constitutes an important step in the unification of CDFT and VB theory.

Extraction Behavior of Acidic Phosphorus-Containing Compounds to Some Metal Ions: A Combination Research of Experimental and Theoretical Investigations

To provide feasible methods for the extraction of valuable metals from spent batteries or low-grade primary ores, the extraction behavior of some representative acidic phosphorus-containing compounds (APCCs) as extractants is evaluated from the perspective of experimental and theoretical investigations in this work. Aqueous solutions containing five metal ions, Ca(II), Co(II), Mg(II), Mn(II), and Ni(II), were made to simulate leaching liquids, and the extraction of these metals was investigated. A simplified calculated model was used to evaluate the interaction between each extractant and metal ions. The calculation results agree well with the experimental tests in trend. This work not only provides potential extractants for the extraction of valuable metals from spent batteries or low-grade primary ores but also demonstrates the practicability of the simplified calculation model.

Is It Possible To Determine Oxidation States for Atoms in Molecules Using Density-Based Quantities? An Information-Theoretic Approach and Conceptual Density Functional Theory Study

The oxidation state, also called oxidation number, of atoms in molecules is a fundamental chemical concept. It is defined as the charge of an atom in a molecule after the ionic approximation of its heteronuclear bonds is applied. Even though for simple molecules the assignment of oxidation states is straightforward, redundancy and ambiguity do exist for others. In this work, we present a density-based framework to determine the oxidation state using the quantities from the information-theoretic approach. As a proof of concept, we examined six elements for a total of 49 molecules. Strong linear correlations were obtained with Shannon entropy, Ghosh-Berkowitz-Parr entropy, information gain, relative Rényi entropy of orders 2 and 3, and Hirshfeld charge. We also discovered that the crystal radius of elements plays the key role in rationalizing the system dependent nature of these strong linear correlations. The validity and effectiveness of our results were demonstrated by the examples of molecules containing elements with two or more oxidation states. Our results should be applicable to more complicated systems in assigning different oxidation states.

Molecular acidity: An accurate description with information-theoretic approach in density functional reactivity theory

Molecular acidity is one of the important physiochemical properties of a molecular system, yet its accurate calculation and prediction are still an unresolved problem in the literature. In this work, we propose to make use of the quantities from the information-theoretic (IT) approach in density functional reactivity theory and provide an accurate description of molecular acidity from a completely new perspective. To illustrate our point, five different categories of acidic series, singly and doubly substituted benzoic acids, singly substituted benzenesulfinic acids, benzeneseleninic acids, phenols, and alkyl carboxylic acids, have been thoroughly examined. We show that using IT quantities such as Shannon entropy, Fisher information, Ghosh-Berkowitz-Parr entropy, information gain, Onicescu information energy, and relative Rényi entropy, one is able to simultaneously predict experimental pKa values of these different categories of compounds. Because of the universality of the quantities employed in this work, which are all density dependent, our approach should be general and be applicable to other systems as well. © 2017 Wiley Periodicals, Inc.

On the physical interpretation of the nuclear molecular orbital energy

Recently, several groups have extended and implemented molecular orbital (MO) schemes to simultaneously obtain wave functions for electrons and selected nuclei. Many of these schemes employ an extended Hartree-Fock approach as a first step to find approximate electron-nuclear wave functions and energies. Numerous studies conducted with these extended MO methodologies have explored various effects of quantum nuclei on physical and chemical properties. However, to the best of our knowledge no physical interpretation has been assigned to the nuclear molecular orbital energy (NMOE) resulting after solving extended Hartree-Fock equations. This study confirms that the NMOE is directly related to the molecular electrostatic potential at the position of the nucleus.

Competition and selectivity in supramolecular synthesis: structural landscape around 1-(pyridylmethyl)-2,2′-biimidazoles

Three isomeric forms of 1-(pyridylmethyl)-2,2′-biimidazole, A1–A3, have been synthesized and subjected to systematic co-crystallizations with selected hydrogen- and halogen-bond donors in order to explore the impact of electrostatics and geometry on the resulting supramolecular architectures. The solid-state supramolecular behavior of A1–A3 is largely consistent in halogen-bonded co-crystals. Only two types of primary interactions, the N–H⋯N/N⋯H–N homomeric hydrogen-bond interactions responsible for the pairing of biimidazole moieties and the I⋯N(pyridine) halogen bonds responsible for the co-crystal formation and structure extension, are present in these systems. The co-crystallizations with hydrogen-bond donors (carboxylic acids), however, lead to multiple possible structural outcomes because of the presence of the biimidazole–acid N–H⋯OC/N⋯H–O heterosynthon that can compete with biimidazole–biimidazole N–H⋯N/N⋯H–N homosynthon. In addition, the somewhat unpredictable nature of proton transfer makes the hydrogen-bonded co-crystals structurally less consistent than their halogen-bonded counterparts.

Aromaticity and antiaromaticity of substituted fulvene derivatives: perspectives from the information-theoretic approach in density functional reactivity theory

Strong correlations among aromaticity descriptors and information-theoretic quantities are unveiled, providing novel insights about aromaticity and antiaromaticity from different perspectives.

Electronegativity and redox reactions

Using the maximum hardness principle, we show that the oxidation potential of a molecule increases as its electronegativity increases and also increases as its electronegativity in its oxidized state increases.

Scaling properties of information-theoretic quantities in density functional reactivity theory

A number of strong linear correlations between information-theoretic quantities and electron populations for atoms, molecules, and atoms-in-molecules have been disclosed.

Density functional reactivity theory study of SN2 reactions from the information-theoretic perspective

Strong linear correlations were unveiled between barrier heights of bimolecular nucleophilic substitution (S_N2) reactions and quantities from the information-theoretic approach.

MOLECULAR ACIDITY OF BUILDING BLOCKS OF BIOLOGICAL SYSTEMS: A DENSITY FUNCTIONAL REACTIVITY THEORY STUDY

An accurate prediction of the molecular acidity by employing ab initio or density functional approaches for typical molecular systems is still challenging. Recently, we proposed to utilize two quantum descriptors, molecular electrostatic potential (MEP) and the sum of valence natural atomic orbital (NAO) energies on the nucleus of both the acidic atom and leaving proton, to quantitatively evaluate the pKa values. This new approach has been validated by a number of organic and inorganic systems and justified within the framework of density functional reactivity theory (DFRT). In this work, we apply the approach to building blocks of biological systems, namely, 20 natural α-amino acids and 5 DNA/RNA bases, together with a few other biologically relevant species. Our results show that there exists a strong linear correlation between MEP on the nucleus of the N atom and the sum of N 2p NAO energies, with the correlation coefficient R2 = 0.99. Also, we observe that both MEP on the nitrogen nucleus and the sum of N 2p NAO energies correlate well with experimental pKa values, with the correlation coefficient equal to 0.91. Using this established model, we predicted the trend of pKa changes of amino acids in proteins with different dielectric constants. We also applied the model to predict pKa values for dipeptides. Implications of these linear relationships to understand functions and reactivity of biological systems are discussed as well.

A theoretical study on the gas-phase protonation of pyridine and phosphinine derivatives

In this paper, we study the protonation of pyridine and phosphinine derivatives. In particular, the geometries, the amount of charge transfer, and the nature of the created N-H and P-H bonds are discussed, underlying the fundamental differences between the phosphorus and the nitrogen atoms as proton acceptors. Conceptual density functional theory and Bader's quantum theory of atoms-in-molecules are notably used to rationalize these trends and to predict the overall energies of these prototype gas-phase acid-base reactions.

Protophilicity index and protofelicity equalization principle: new measures of Brønsted-Lowry-Lewis acid-base interactions

The simultaneous contributions of proton and electron transfer to the Brønsted-Lowry and Lewis acid-base properties of a set of p-substituted phenols are reported in this work. As a result of the analysis, a novel protophilicity index considered as the second-order energy change of a Brønsted-Lowry base as it is saturated with protons, a combined Brønsted-Lowry-Lewis acidity index (with a corresponding basicity index), and a protofelicity equalization principle (a parallel of the electronegativity equalization principle) are presented.

Molecular electrostatic potential at the atomic sites in the effective core potential approximation

Considering calculations of the molecular electrostatic potential at the atomic sites (MEP@AS) in the presence of effective core potentials (ECP), we found that the consequent use of the definition of MEP@AS based on the energy derivative with respect to nuclear charge leads to a formula that differs by one term from the result of simple application of Coulomb's law. We have developed a general method to analytically treat derivatives of ECP with respect to nuclear charge. Benchmarking calculations performed on a set of simple molecules show that our formula leads to a systematic decrease in the error connected with the introduction of ECP when compared to all-electron results. Because of a straightforward implementation and relatively low costs of the developed procedure we suggest to use it by default.

Lewis molecular acidity of ionic liquids from empirical energy-density models

Two complementary models of Lewis molecular acidity are introduced and tested in a wide series of 45 room temperature ionic liquids (RTIL). They are defined in the context of the conceptual density functional theory. The first one, which we tentatively call the excess electronic chemical potential, assesses the electron accepting power of the RTIL by relating the H-bond donor acidity with the charge transfer associated to the acidic H-atom migration at the cation of the RTIL considered as a HB-donor species. This global index accounts for the molecular acidity of the cation moiety of the ionic liquid that takes into account the perturbation of the anionic partner. The second index is defined in terms of the local charge capacity modeled through the maximum electronic charge that the cation, in its valence state, may accept from an unspecified environment. Each model is compared with the experimental HB-donor acidity parameter of the Kamlet Taft model. The best comparison is obtained for a combination of both the excess electronic chemical potential and the local charge capacity. As expected, the correlations with the Kamlet Taft α parameter do not lead to a universal model of HB-donor acidity. Reduced correlations for limited series of structurally related RTIL are obtained instead. Finally, we illustrate the reliability and usefulness of the proposed model of RTIL molecular acidity to explain the cation-dependent solvent effects on the reactivity trends for cycloaddition, Kemp elimination, and Menschutkin reactions, for which experimental rate coefficients are available from literature.

Empirical prediction of protein pKa values with residue mutation

A fast, empirical method, Mut-pKa, is presented for predicting the pKa values of ionizable residues in proteins based on mutation. The method compares the effect of mutating each residue that may act as a hydrogen bond donor or acceptor for the ionizable residue. The energetic effect of each type of mutation, along with a desolvation measure and the overall background charge, is fit against pKa data for histidine and carboxyl residues. A total of 214 residues from 35 different proteins were used in the dataset. Using 11 parameters for each type of ionizable residue, a root mean squared error (RMSE) of 0.78 and 1.12 pH units were obtained for carboxyl and histidines residues, respectively, using leave one out cross validation (LOOCV). The results were particularly promising for buried residues, which had RMSE values of 0.99 and 1.13 for carboxyl and histidine residues, respectively. A number of desolvation measures were tested. The simplest measure, the number of atoms surrounding the residue, was found to work best. The effect of using dynamics was also studied using short molecular dynamics runs, followed by minimization of the structures. Mut-pKa has significantly fewer parameters than, but similar performance to, other empirical methods. Because of this and the LOOCV results, we believe the model is robust and that overfitting is not a problem.

Lee Pedersen’s work in theoretical and computational chemistry and biochemistry

Nature at the lab level in biology and chemistry can be described by the application of quantum mechanics. In many cases, a reasonable approximation to quantum mechanics is classical mechanics realized through Newton’s equations of motion. Dr. Pedersen began his career using quantum mechanics to describe the properties of small molecular complexes that could serve as models for biochemical systems. To describe large molecular systems required a drop-back to classical means and this led surprisingly to a major improvement in the classical treatment of electrostatics for all molecules, not just biological molecules. Recent work has involved the application of quantum mechanics for the putative active sites of enzymes to gain greater insight into the key steps in enzyme catalysis.

Steric, quantum, and electrostatic effects on S(N)2 reaction barriers in gas phase

Biomolecular nucleophilic substitution reactions, S(N)2, are fundamental and commonplace in chemistry. It is the well-documented experimental finding in the literature that vicinal substitution with bulkier groups near the reaction center significantly slows the reaction due to steric hindrance, but theoretical understanding in the quantitative manner about factors dictating the S(N)2 reaction barrier height is still controversial. In this work, employing the new quantification approach that we recently proposed for the steric effect from the density functional theory framework, we investigate the relative contribution of three independent effects-steric, electrostatic, and quantum-to the S(N)2 barrier heights in gas phase for substituted methyl halide systems, R(1)R(2)R(3)CX, reacting with the fluorine anion, where R(1), R(2), and R(3) denote substituting groups and X = F or Cl. We found that in accordance with the experimental finding, for these systems, the steric effect dominates the transition state barrier, contributing positively to barrier heights, but this contribution is largely compensated by the negative, stabilizing contribution from the quantum effect due to the exchange-correlation interactions. Moreover, we find that it is the component from the electrostatic effect that is linearly correlated with the S(N)2 barrier height for the systems investigated in the present study. In addition, we compared our approach with the conventional method of energy decomposition in density functional theory as well as examined the steric effect from the wave function theory for these systems via natural bond orbital analysis.