Comparison among Four Different Ways to Condense the Fukui Function

Molecular interactions from the density functional theory for chemical reactivity: Interaction chemical potential, hardness, and reactivity principles

In the first paper of this series, the authors derived an expression for the interaction energy between two reagents in terms of the chemical reactivity indicators that can be derived from density functional perturbation theory. While negative interaction energies can explain reactivity, reactivity is often more simply explained using the “|dμ| big is good” rule or the maximum hardness principle. Expressions for the change in chemical potential (μ) and hardness when two reagents interact are derived. A partial justification for the maximum hardness principle is that the terms that appear in the interaction energy expression often reappear in the expression for the interaction hardness, but with opposite sign.

Transference Number Reinforced-Based Gel Copolymer Electrolyte for Dendrite-Free Lithium Metal Batteries

The progress of electric vehicles is highly inhibited by the limited energy density and growth of dendrite Li in current power batteries. Breakthroughs and improvements in electrolyte chemistry are highlighted to directly address the above issues, namely, the development of electrolytes with a high lithium-ion transference number (tLi+), enabling one to effectively restrict the concentration polarization during repetitious cycling. Herein, we propose a novel ether-based copolymer-based gel polymer electrolyte (ECP-based GPE) by in situ copolymerization as an intriguing strategy to achieve a high tLi+ of ∼0.64. Molecular dynamics simulations and finite element method analyses illustrate the enhanced Li+ diffusion process (DLi+, ∼1.76 × 10-10 m2 s-1) in ECP-based GPE with a homogeneous electric potential accommodated around the lithium metal anode. Therefore, such a high-tLi+-based electrolyte renders a high reversibility of dendrite-free lithium plating/stripping at a high areal capacity (5 mA cm-2/5 mA h cm-2) in an Li||Li symmetric cell and facilitates superior cycling performances (over 1000 cycles) at a high rate (5 C) with a capacity retention of ∼91.1% in Li||LiFePO4 batteries, promoting the practical application of solid-state lithium metal batteries.

Anisotropic growth of ZnO nanoparticles driven by the structure of amine surfactants: the role of surface dynamics in nanocrystal growth

Herein, we elucidate the key role of amine surfactants in the controlled anisotropic growth of ZnO nanoparticles that is achieved under mild conditions by organometallic hydrolysis. The structuring influence of alkyl substituents on the nitrogen atom of amines is jointly analyzed theoretically by DFT modeling, and experimentally by multinuclear NMR (¹H, ¹³C and ¹⁷O) spectroscopy. We demonstrate that in initial steps leading to the growth of colloidal ZnO particles, the nature of molecular species that are involved in the solution strongly depends on the structure of the amine surfactant. By using tertiary, secondary or primary amines, no or weak adducts between the amine and zinc, or stable adducts, or adduct oligomers were identified, respectively. Afterwards, following the course of the reaction, the dynamic behavior of the amines on the grown ZnO nanocrystal surfaces is also strongly correlated with their structure. We identified that in the presence of tertiary, secondary or primary amines, no significant [Zn⋯N] adsorption, or surface adsorption with notable surface mobility, or a very strong adsorption is achieved, respectively. The last case, primary amines, significantly involves the structuring of a hydrogen bonding network. Therefore, such surface dynamic behavior has a predominant role in driving the nanocrystal growth, and orienting the ZnO material final morphology. By forming hydrogen bonds at the nanoparticle surface during the growth process, primary amines specifically lead to the formation of nanorods. Conversely, isotropic nanoparticles and aggregates are obtained when secondary and tertiary amines are used, respectively. These findings shed light on the role of weak surface interactions, herein H-bonding, that rule the growth of nano-objects and are as such crucial to identify, study, and control for achieving progress in nanoscience.

Kick-Fukui: A Fukui Function-Guided Method for Molecular Structure Prediction

Here, we introduce a hybrid method, named Kick-Fukui, to explore the potential energy surface (PES) of clusters and molecules using the Coulombic integral between the Fukui functions in the first screening of the best individuals. In the process, small stable molecules or clusters whose combination has the stoichiometry of the explored species are used as assembly units. First, a small set of candidates has been selected from a large and stochastically generated (Kick) population according to the maximum value of the Coulombic integral between the Fukui functions of both fragments. Subsequently, these few candidates are optimized using a gradient method and density functional theory (DFT) calculations. The performance of the program has been evaluated to explore the PES of various systems, including atomic and molecular clusters. In most cases studied, the global minimum (GM) has been identified with a low computational cost. The strategy does not allow to identify the GM of some silicon clusters; however, it predicts local minima very close in energy to the GM that could be used as the initial population of evolutionary algorithms.

New insights in chemical reactivity from quantum chemical topology

Based on the quantum chemical topology of the modified electron localization function ELF_x , an efficient and robust mechanistic methodology designed to identify the favorable reaction pathway between two reactants is proposed. We first recall and reshape how the supermolecular interaction energy can be evaluated from only three distinct terms, namely the intermolecular coulomb energy, the intermolecular exchange-correlation energy and the intramolecular energies of reactants. Thereafter, we show that the reactivity between the reactants is driven by the first-order variation in the coulomb intermolecular energy defined in terms of the response to changes in the number of electrons. Illustrative examples with the formation of the dative bond B-N involved in the BH₃ NH₃ molecule and the typical formation of the hydrogen bond in the canonical water dimer are presented. For these selected systems, our approach unveils a noticeable mimicking of E_dual onto the DFT intermolecular interaction energy surface calculated between the both reactants. An automated reaction-path algorithm aimed to determine the most favorable relative orientations when the two molecules approach each other is also outlined.

Structure-activity study of furyl aryloxazole fluorescent probes for the detection of singlet oxygen

In this study, we report the synthesis and the photochemical behavior of a series of new "click-on" fluorescent probes designed to detect singlet oxygen. They include a highly fluorescent chemical structure, an aryloxazole ring, linked to a furan moiety operating as singlet oxygen trap. Their activity depends on both the structure of the aryloxazole fluorophore and the electron-donating and electron-accepting properties of the substituents attached to the C-5 of the furan ring. All probes are selectively oxidized by singlet oxygen to give a single fluorescent product in methanol and produce negligible amounts of singlet oxygen themselves by self-sensitization. The most promising dyad, (E)-2-(2-(5-methylfuran-2-yl)vinyl)naphtho[1,2-d]oxazole, FN-6, shows outstanding reactivity and sensitivity: it traps singlet oxygen with a rate constant (5,8 ± 0.1) x 1(07) M-1 s-1 and its fluorescence increases by a factor of 500 upon reaction. Analysis of the dyads reactivity in terms of linear free energy relationships using the modified Swain and Lupton parameter F and the Fukui condensed function for the electrophilic attack, suggests that cycloaddition of singlet oxygen to the furan ring is partially concerted and possibly involves an exciplex with a "more open" structure than could be expected for a concerted cycloaddition.

Synthesis, spectroscopic characterization (FT-IR, FT-Raman, and NMR), quantum chemical studies and molecular docking of 3-(1-(phenylamino)ethylidene)-chroman-2,4-dione

The experimental and theoretical investigations of structure of the 3-(1-(phenylamino)ethylidene)-chroman-2,4-dione were performed. X-ray structure analysis and spectroscopic methods (FTIR and FT-Raman, ¹H and ¹³C NMR), along with the density functional theory calculations (B3LYP functional with empirical dispersion corrections D3BJ in combination with the 6-311 + G(d,p) basis set), were used in order to characterize the molecular structure and spectroscopic behavior of the investigated coumarin derivative. Molecular docking analysis was carried out to identify the potency of inhibition of the title molecule against human's Ubiquinol-Cytochrome C Reductase Binding Protein (UQCRB) and Methylenetetrahydrofolate reductase (MTHFR). The inhibition activity was obtained for ten conformations of ligand inside the proteins.

Predicting the Metabolic Sites by Flavin-Containing Monooxygenase on Drug Molecules Using SVM Classification on Computed Quantum Mechanics and Circular Fingerprints Molecular Descriptors

As an important enzyme in Phase I drug metabolism, the flavin-containing monooxygenase (FMO) also metabolizes some xenobiotics with soft nucleophiles. The site of metabolism (SOM) on a molecule is the site where the metabolic reaction is exerted by an enzyme. Accurate prediction of SOMs on drug molecules will assist the search for drug leads during the optimization process. Here, some quantum mechanics features such as the condensed Fukui function and attributes from circular fingerprints (called Molprint2D) are computed and classified using the support vector machine (SVM) for predicting some potential SOMs on a series of drugs that can be metabolized by FMO enzymes. The condensed Fukui function fA- representing the nucleophilicity of central atom A and the attributes from circular fingerprints accounting the influence of neighbors on the central atom. The total number of FMO substrates and non-substrates collected in the study is 85 and they are equally divided into the training and test sets with each carrying roughly the same number of potential SOMs. However, only N-oxidation and S-oxidation features were considered in the prediction since the available C-oxidation data was scarce. In the training process, the LibSVM package of WEKA package and the option of 10-fold cross validation are employed. The prediction performance on the test set evaluated by accuracy, Matthews correlation coefficient and area under ROC curve computed are 0.829, 0.659, and 0.877 respectively. This work reveals that the SVM model built can accurately predict the potential SOMs for drug molecules that are metabolizable by the FMO enzymes.

Understanding the Highly Varying pKa of Arylamines. A Perspective from the Average Local Ionization Condensed-to-Atom Framework

The highly varying experimental pKa values for 36 arylamines spanning 7 orders of magnitude is carefully examined. Within this framework, a valence condensed-to-atom model for the average ionization energy is introduced and tested. The theoretical approach is connected to orbital Fukui functions directly mapped into semilocal or regional site-specific responses. It is revealed that the average local ionization energies associated with the amino nitrogen atom is linearly correlated to the basicity of the substituted arylamines, properly reproducing the experimental ordering of basicity. The condensed-to-atom descriptor exhibits a high predictive power, providing a new direct reactivity evaluation of significant value.

Condensed Fukui function predicts innate C–H radical functionalization sites on multi-nitrogen containing fused arenes

The condensed Fukui function (values in blue and red) predicts innate radical C–H functionalization sites (dark green circles) observed in experiments.

An information-theoretic resolution of the ambiguity in the local hardness

A definition of the local hardness, suitable for application in the local hard/soft acid/base principle, is derived by applying information theory.

Computational nanochemistry report on the oxicams--conceptual DFT indices and chemical reactivity

A density functional theory study of eight oxicams was carried out in order to determine their global and local reactivities. These types of reactivities were measured by means of global and local reactivity descriptors coming from the conceptual density functional theory. Net electrophilicity as a global reactivity descriptor and local hypersoftness as a local reactivity descriptor were the used tools to distinguish reactivity and selectivity among these oxicams. Globally, isoxicam presents the highest electron donating capacity; meanwhile, the highest electron accepting capacity is exhibited by droxicam. Locally, two oxicams present neither nucleophilic nor electrophilic relevant reactivity in their peripheral pyridine ring, droxicam and tenoxicam, so that their more reactive zones are found on the respective fused rings. Oxicams have been divided into two subgroups in order to facilitate the local analysis of reactivity. One group is characterized because their most important condensed values for local hypersoftnes are well-separated: 4-meloxicam, lornoxicam, meloxicam, and normeloxicam. Meanwhile, the opposite situation is found in droxicam, isoxicam, piroxicam, and tenoxicam. As a whole, the nucleophilic characteristic noticeably predominates in these eight oxicams instead of an electrophilic behavior, thus meaning a greater tendency to donate electrons rather than withdrawing them; a consequence of this behavior implies a favorable interaction with a hypothetical receptor bearing one or more electron acceptor functional groups rather than electron donor functional groups; this would imply a maximization of this interaction from the covalent point of view.

Global and local chemical reactivities of mutagen X and simple derivatives

Registered by the World Health Organization (WHO), 3-chloro-4-(dichloromethyl)-5-hydroxy-2(5H)-furanone (MX) is one of the strongest bacterial mutagens ever tested, as highlighted by the Ames Salmonella typhimurium TA100 assay. We provide new insights concerning this mutagenic activity on the basis of global and local theoretically defined electrophilicity indices. Our results further support the idea that mutagenicity of MX and its analogues is related more closely to one-electron transfer processes from the electron-rich biological environment than to adduct formation processes. We also stress that, although the Z-open tautomers are intrinsically more electrophilic than furanone ring analogues, the observed mutagenic activity is significantly correlated only to the electrophilicity response of the ring forms. In that context, we also emphasize that it is electrophilicity at the C α in the α-β unsaturated carbonyl moiety that exhibits a strong correlation with the observed mutagenic activity.

Why is quercetin a better antioxidant than taxifolin? Theoretical study of mechanisms involving activated forms

The stronger antioxidant capacity of the flavonoid quercetin (Q) compared with taxifolin (dihydroquercetin, T) has been the subject of previous experimental and theoretical studies. Theoretical work has focused on the analysis of hydrogen bond dissociation energies (BDE) of the OH phenolic groups, but consider mechanisms that only involve the transfer of one hydrogen atom. In the present work we consider other mechanisms involving a second hydrogen transfer in reactions with free radicals. The relative stability of the radicals formed after the first hydrogen transfer reaction is considered in discussing the antioxidant activity of Q and T. In terms of global and local theoretical reactivity descriptors, we propose that the radical arising from Q should be more persistent in the environment and with the capability to react with a second radical by hydrogen transfer, proton transfer and electron transfer mechanisms. These mechanisms could be responsible of the stronger antioxidant capacity of Q.

In pursuit of negative Fukui functions: examples where the highest occupied molecular orbital fails to dominate the chemical reactivity

In our quest to explore molecules with chemically significant regions where the Fukui function is negative, we explored reactions where the frontier orbital that indicates the sites for electrophilic attack is not the highest occupied molecular orbital. The highest occupied molecular orbital (HOMO) controls the location of the regions where the Fukui function is negative, supporting the postulate that negative values of the Fukui function are associated with orbital relaxation effects and nodal surfaces of the frontier orbitals. Significant negative values for the condensed Fukui function, however, were not observed.

Singlet oxygen reactions with flavonoids. A theoretical-experimental study

Detection of singlet oxygen emission, λ(max) = 1270 nm, following laser excitation and steady-state methods were employed to measure the total reaction rate constant, k(T), and the reactive reaction rate constant, k(r), for the reaction between singlet oxygen and several flavonoids. Values of k(T) determined in deuterated water, ranging from 2.4×10(7) M(-1) s(-1) to 13.4×10(7) M(-1) s(-1), for rutin and morin, respectively, and the values measured for k(r), ranging from 2.8×10(5) M(-1) s(-1) to 65.7×10(5) M(-1) s(-1) for kaempferol and morin, respectively, being epicatechin and catechin chemically unreactive. These results indicate that all the studied flavonoids are good quenchers of singlet oxygen and could be valuable antioxidants in systems under oxidative stress, in particular if a flavonoid-rich diet was previously consumed. Analysis of the dependence of rate constant values with molecular structure in terms of global descriptors and condensed Fukui functions, resulting from electronic structure calculations, supports the formation of a charge transfer exciplex in all studied reactions. The fraction of exciplex giving reaction products evolves through a hydroperoxide and/or an endoperoxide intermediate produced by singlet oxygen attack on the double bond of the ring C of the flavonoid.

Revisiting caffeate's capabilities as a complexation agent to silver cation in mining processes by means of the dual descriptor--a conceptual DFT approach

Caffeic acid (C(9)H(8)O(4)) and its conjugate base C(9)H(7)O(4) (-) (anionic form-known as caffeate) were analyzed computationally through the use of quantum chemistry to assess their intrinsic global and local reactivity using the tools of conceptual density functional theory. The anionic form was found to be better at coordinating the silver cation than caffeic acid thus suggesting the use of caffeate as a complexation agent. The complexation capability of caffeate was compared with that of some of the most common ligand agents used to coordinate silver cations. Local reactivity descriptors allowed identification of the preferred sites on caffeate for silver cation coordination thus generating a plausible silver complex. All silver complexes were analyzed thermodynamically considering interaction energies in both gas and aqueous phases; the complexation free energy in aqueous phase was also determined. These results suggest that more attention be paid to the caffeate anion and its derivatives because this work has shed new light on the behavior of this anion in the recovery of silver cations that could be exploited in silver mining processes in a environmentally friendly way.

Correlation between electron localization and metal ion mutagenicity in DNA synthesis from QM/MM calculations

DNA polymerases require two divalent metal ions in the active site for catalysis. Mg(2+) has been confirmed to be the most probable cation utilized by most polymerases in vivo. Other metal ions are either potent mutagens or inhibitors. We used structural and topological analyses based on ab initio QM/MM calculations to study human DNA polymerase λ (Polλ) with different metals in the active site. Our results indicate a slightly longer O3'-Pα distance (∼3.6 Å) for most inhibitor cations compared to the natural and mutagenic metals (∼3.3-3.4 Å). Optimization with a larger basis set for the previously reported transition state (TS) structures (Cisneros et al., DNA Repair, 2008, 7, 1824.) gives barriers of 17.4 kcal mol(-1) and 15.1 kcal mol(-1) for the Mg(2+) and Mn(2+) catalyzed reactions respectively. Relying on the key relation between the topological signature of a metal cation and its selectivity within biological systems (de Courcy et al., J. Chem. Theor. Comput., 2010, 6, 1048.) we have performed electron localization function (ELF) topological analyses. These analyses show that all inhibitor and mutagenic metals considered, except Na(+), present a "split" of the outer-shell density of the metal. This "splitting" is not observed for the non-mutagenic Mg(2+) metal. Population and multipole analyses on the ELF basins reveal that the electronic dipolar and quadrupolar polarization is significantly different with Mg(2+) compared to all other cations. Our results shed light at the atomic level on the subtle differences between Mg(2+), mutagenic, and inhibitor metals in DNA polymerases. These results provide a correlation between the electronic distribution of the cations in the active site and the possible consequences on DNA synthesis.

Understanding selectivity of hard and soft metal cations within biological systems using the subvalence concept. I. Application to blood coagulation: direct cation-protein electronic effects vs. indirect interactions through water networks

Following a previous study by de Courcy et al. ((2009) Interdiscip. Sci. Comput. Life Sci. 1, 55-60), we demonstrate in this contribution, using quantum chemistry, that metal cations exhibit a specific topological signature in the electron localization of their density interacting with ligands according to its "soft" or "hard" character. Introducing the concept of metal cation subvalence, we show that a metal cation can split its outer-shell density (the so-called subvalent domains or basins) according to it capability to form a partly covalent bond involving charge transfer. Such behaviour is investigated by means of several quantum chemical interpretative methods encompasing the topological analysis of the Electron Localization Function (ELF) and Bader's Quantum Theory of Atoms in Molecules (QTAIM) and two energy decomposition analyses (EDA), namely the Restricted Variational Space (RVS) and Constrained Space Orbital Variations (CSOV) approaches. Further rationalization is performed by computing ELF and QTAIM local properties such as electrostatic distributed moments and local chemical descriptors such as condensed Fukui Functions and dual descriptors. These reactivity indexes are computed within the ELF topological analysis in addition to QTAIM offering access to non atomic reactivity local index, for example on lone pairs. We apply this "subvalence" concept to study the cation selectivity in enzymes involved in blood coagulation (GLA domains of three coagulation factors). We show that the calcium ions are clearly able to form partially covalent charge transfer networks between the subdomain of the metal ion and the carboxylate oxygen lone pairs whereas magnesium does not have such ability. Our analysis also explains the different role of two groups (high affinity and low affinity cation binding sites) present in GLA domains. If the presence of Ca(II) is mandatory in the central "high affinity" region to conserve a proper folding and a charge transfer network, external sites are better stabilised by Mg(II), rather than Ca(II), in agreement with experiment. The central role of discrete water molecules is also discussed in order to understand the stabilities of the observed X-rays structures of the Gla domain. Indeed, the presence of explicit water molecules generating indirect cation-protein interactions through water networks is shown to be able to reverse the observed electronic selectivity occuring when cations directly interact with the Gla domain without the need of water.

Addition of Hetero Allenyl Copper Reagents to Aldehydes: Scope and Behavior

Cyanocuprates derived from propargylic amines or ethers react with aldehydes to give regioselectively the corresponding anti-homopropargylic alcohols with a high level of diastereoselectivity. Such selectivity could be obtained independently of the nature of the heteroatom (amine or ethers) or the acetylenic substituents. Excellent selectivities can be reached regardless of the aldehydes used, remarkably also with vinylic or acetylenic ones. A reactivity scale for the cuprates bearing different acetylenic substituents was established to be: SiMe3 >Ph> Et. The rate of the addition reaction to aldehydes was also found to be slowed down in the presence of HMPA underlining the crucial role of the Li counterion. DFT calculations have shown that the relationship between the rate and the acetylenic substituent is not connected to a possible metallotropic equilibrium but to the stability of the reactive allenic species compared to the less-stable propargylic isomer.