Advancing beyond charge analysis using the electronic localization function: Chemically intuitive distribution of electrostatic moments

A comprehensive electron wavefunction analysis toolbox for chemists, Multiwfn

Analysis of electron wavefunction is a key component of quantum chemistry investigations and is indispensable for the practical research of many chemical problems. After more than ten years of active development, the wavefunction analysis program Multiwfn has accumulated very rich functions, and its application scope has covered numerous aspects of theoretical chemical research, including charge distribution, chemical bond, electron localization and delocalization, aromaticity, intramolecular and intermolecular interactions, electronic excitation, and response property. This article systematically introduces the features and functions of the latest version of Multiwfn and provides many representative examples. Through this article, readers will be able to fully understand the characteristics and recognize the unique value of Multiwfn. The source code and precompiled executable files of Multiwfn, as well as the manual containing a detailed introduction to theoretical backgrounds and very rich tutorials, can all be downloaded for free from the Multiwfn website (http://sobereva.com/multiwfn).

Exploring the effects of mono-bromination on hole-electron transport and distribution in dibenzofuran and dibenzothiophene isomers: a first-principles study

CONTEXT

This study delves into hole-electron transport and distribution properties inherent in mono-brominated dibenzofuran (DBF) and dibenzothiophene (DBT) isomers. As determined by frontier molecular orbitals, all brominated structures have narrower bandgaps than their primary structures. The TD-DFT calculation showed that 2BDBT had the highest absorption wavelength of all molecules at 315.35 nm. Notably, the study unveils remarkably low electron and hole reorganization energies due to bromine substitution in DBF and DBT molecules. Specifically, the 4BDBF has the lowest hole reorganization energy of all DBF configurations, 0.229 eV. In addition, 3BDBF has 0.226 eV less electron reorganization energy than all other molecules. Compared to DBT, 3BDBT has the lowest electron reorganization energy of 0.254 eV. Overall, this research sheds significant light on the fundamental electronic and hole transport characteristics of bromine-substituted DBF and DBT isomers, highlighting their promising role in polymer design as donors/acceptors for advanced organic electronic applications.

METHODS

Molecular structures were optimized using Density Functional Theory (DFT) B3LYP/6-311 +  + G (d, p) level of theory, and the study further elucidates these molecules' energy levels and absorption spectra through Time-Dependent Density Functional Theory TD-DFT; these calculations were performed using Gaussian 09W software package. The key parameters such as reorganization energies, Electron Localization Function map, Laplacian Bond Order, and NCI-RDG were meticulously examined for the molecules with the results of DFT calculations were analyzed and displayed by utilizing the software packages VMD 1.9.4 and Multiwfn 3.8, aiming to comprehend their charge transport and distribution properties.

Distal Mutations in the β-Clamp of DNA Polymerase III* Disrupt DNA Orientation and Affect Exonuclease Activity

DNA polymerases are responsible for the replication and repair of DNA found in all DNA-based organisms. DNA Polymerase III is the main replicative polymerase of E. coli and is composed of over 10 proteins. A subset of these proteins (Pol III*) includes the polymerase (α), exonuclease (ϵ), clamp (β), and accessory protein (θ). Mutations of residues in, or around the active site of the catalytic subunits (α and ϵ), can have a significant impact on catalysis. However, the effects of distal mutations in noncatalytic subunits on the activity of catalytic subunits are less well-characterized. Here, we investigate the effects of two Pol III* variants, β-L82E/L82'E and β-L82D/L82'D, on the proofreading reaction catalyzed by ϵ. MD simulations reveal major changes in the dynamics of Pol III*, which extend throughout the complex. These changes are mostly induced by a shift in the position of the DNA substrate inside the β-clamp, although no major structural changes are observed in the protein complex. Quantum mechanics/molecular mechanics (QM/MM) calculations indicate that the β-L82D/L82'D variant has reduced catalytic proficiency due to highly endoergic reaction energies resulting from structural changes in the active site and differences in the electric field at the active site arising from the protein and substrate. Conversely, the β-L82E/L82'E variant is predicted to maintain proofreading activity, exhibiting a similar reaction barrier for nucleotide excision compared with the WT system. However, significant differences in the reaction mechanism are obtained due to the changes induced by the mutations on the β-clamp.

Combining Evolutionary Conservation and Quantum Topological Analyses To Determine Quantum Mechanics Subsystems for Biomolecular Quantum Mechanics/Molecular Mechanics Simulations

Selection of residues and other molecular fragments for inclusion in the quantum mechanics (QM) region for QM/molecular mechanics (MM) simulations is an important step for these calculations. Here, we present an approach that combines protein sequence/structure evolution and electron localization function (ELF) analyses. The combination of these two analyses allows the determination of whether a residue needs to be included in the QM subsystem or can be represented by the MM environment. We have applied this approach on two systems previously investigated by QM/MM simulations, 4-oxalocrotonate tautomerase (4OT) and ten-eleven translocation-2 (TET2), that provide examples where fragments may or may not need to be included in the QM subsystem. Subsequently, we present the use of this approach to determine the appropriate QM subsystem to calculate the minimum energy path (MEP) for the reaction catalyzed by human DNA polymerase λ (Polλ) with a third cation in the active site. Our results suggest that the combination of protein evolutionary and ELF analyses provides insights into residue/molecular fragment selection for QM/MM simulations.

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.

Polarizable ab initio QM/MM Study of the Reaction Mechanism of N-tert-Butyloxycarbonylation of Aniline in [EMIm][BF₄]

N-t e r t-butoxycarbonylation of amines in solution (water, organic solvents, or ionic liquids) is a common reaction for the preparation of drug molecules. To understand the reaction mechanism and the role of the solvent, quantum mechanical/molecular mechanical simulations using a polarizable multipolar force field with long⁻range electrostatic corrections were used to optimize the minimum energy paths (MEPs) associated with various possible reaction mechanisms employing the nudged elastic band (NEB) and the quadratic string method (QSM). The calculated reaction energies and energy barriers were compared with the corresponding gas-phase and dichloromethane results. Complementary Electron Localization Function (ELF)/NCI analyses provide insights on the critical structures along the MEP. The calculated results suggest the most likely path involves a sequential mechanism with the rate⁻limiting step corresponding to the nucleophilic attack of the aniline, followed by proton transfer and the release of CO 2 without the direct involvement of imidazolium cations as catalysts.

Quantum Chemical Topology of the Electron Localization Function in the Field of Attosecond Electron Dynamics

We report original analyses of attosecond electron dynamics of molecules subject to collisions by high energy charged particles based on Real-Time Time-Dependent-Density-Functional-Theory simulations coupled to Topological Analyses of the Electron Localization Function (TA-TD-ELF). We investigate irradiation of water and guanine. TA-TD-ELF enables qualitative and quantitative characterizations of bond breaking and formation, of charge migration within topological basins, or of electron attachment to the colliding particle. Whereas the Lewis-VSEPR structure of gas phase water is blown out within a few attoseconds after collision, that of guanine is far more robust and reconstitutes rapidly after impact even though the molecule remains electronically excited. This difference is accounted by the presence of the electron bath surrounding the impact point which enables energy relaxation within the molecule. Our approach should stimulate future studies to unravel the early steps following irradiation of various types of systems (isolated molecules, biomolecules, nanoclusters, solids, etc.) and is also readily applicable to irradiation by photons of various energies.

The bonding picture in hypervalent XF3 (X = Cl, Br, I, At) fluorides revisited with quantum chemical topology

Hypervalent XF₃ (X = Cl, Br, I, At) fluorides exhibit T-shaped C_2V equilibrium structures with the heavier of them, AtF₃ , also revealing an almost isoenergetic planar D_3h structure. Factors explaining this behavior based on simple "chemical intuition" are currently missing. In this work, we combine non-relativistic (ClF₃ ), scalar-relativistic and two-component (X = Br - At) density functional theory calculations, and bonding analyses based on the electron localization function and the quantum theory of atoms in molecules. Typical signatures of charge-shift bonding have been identified at the bent T-shaped structures of ClF₃ and BrF₃ , while the bonds of the other structures exhibit a dominant ionic character. With the aim of explaining the D_3h structure of AtF₃ , we extend the multipole expansion analysis to the framework of two-component single-reference calculations. This methodological advance enables us to rationalize the relative stability of the T-shaped C_2v and the planar D_3h structures: the Coulomb repulsions between the two lone-pairs of the central atom and between each lone-pair and each fluorine ligand are found significantly larger at the D_3h structures than at the C_2v ones for X = Cl - I, but not with X = At. This comes with the increasing stabilization, along the XF₃ series, of the planar D_3h structure with respect to the global T-shaped C_2v minima. Hence, we show that the careful use of principles that are at the heart of the valence shell electron pair repulsion model provides reasonable justifications for stable planar D_3h structures in AX₃ E₂ systems. © 2017 Wiley Periodicals, Inc.

Dimerization of quercetin, Diels-Alder vs. radical-coupling approach: a joint thermodynamics, kinetics, and topological study

Quercetin is a prototypical antioxidant and prominent member of flavonoids, a large group of natural polyphenols. The oxidation of quercetin may lead to its dimerization, which is a paradigm of the more general polyphenol oligomerization. There exist two opposing mechanisms to describe the dimerization process, namely radical-coupling or Diels-Alder reactions. This work presents a comprehensive rationalization of this dimerization process, acquired from density functional theory (DFT) calculations. It is found that the two-step radical-coupling pathway is thermodynamically and kinetically preferred over the Diels-Alder reaction. This is in agreement with the experimental results showing the formation of only one isomer, whereas the Diels-Alder mechanism would yield two isomers. The evolution in bonding, occurring during these two processes, is investigated using the atoms in molecules (AIM) and electron localization function (ELF) topological approaches. It is shown that some electron density is accumulated between the fragments in the transition state of the radical-coupling reaction, but not in the transition state of the Diels-Alder process. Graphical Abstract Quantum chemistry calculations of the dimerization process of quercetin show that a radical coupling approach is preferred to a Diels-Alder type reaction, in agreement with experimental results. Analysis of the bonding evolution highlights the reaction mechanism.

Non-adiabatic dynamics simulation exploration of the wavelength-dependent photoinduced relaxation mechanism of trans-N-1-methyl-2-(tolylazo) imidazole in the gas phase

A comprehensive picture of the photoinduced non-adiabatic relaxation dynamics of trans-N-1-methyl-2-(tolylazo) imidazole (trans-MTAI) in different electronic excited states has been revealed using the on-the-fly surface hopping method at the ab initio CASSCF level.

QTAIM Analysis in the Context of Quasirelativistic Quantum Calculations

Computational chemistry currently lacks ad hoc tools for probing the nature of chemical bonds in heavy and superheavy-atom systems where the consideration of spin-orbit coupling (SOC) effects is mandatory. We report an implementation of the Quantum Theory of Atoms-In-Molecules in the framework of two-component relativistic calculations. Used in conjunction with the topological analysis of the Electron Localization Function, we show for astatine (At) species that SOC significantly lowers At electronegativity and boosts its propensity to make charge-shift bonds. Relativistic spin-dependent effects are furthermore able to change some bonds from mainly covalent to charge-shift type. The implication of the disclosed features regarding the rationalization of the labeling protocols used in nuclear medicine for (211)At radioisotope nicely illustrates the potential of the introduced methodology for investigating the chemistry of (super)heavy elements.

Revisiting H2O Nucleation around Au+ and Hg2+: The Peculiar "Pseudo-Soft" Character of the Gold Cation

In this contribution, we propose a deeper understanding of the electronic effects affecting the nucleation of water around the Au+ and Hg2+ metal cations using quantum chemistry. To do so, and in order to go beyond usual energetical studies, we make extensive use of state of the art quantum interpretative techniques combining ELF/NCI/QTAIM/EDA computations to capture all ranges of interactions stabilizing the well characterized microhydrated structures. The Electron Localization Function (ELF) topological analysis reveals the peculiar role of the Au+ outer-shell core electrons (subvalence) that appear already spatially preorganized once the addition of the first water molecule occurs. Thus, despite the addition of other water molecules, the electronic structure of Au(H2O)+ appears frozen due to relativistic effects leading to a maximal acceptation of only two waters in gold's first hydration shell. As the values of the QTAIM (Quantum Theory of Atoms in Molecules) cations's charge is discussed, the Non Covalent Interactions (NCI) analysis showed that Au+ appears still able to interact through longer range van der Waals interaction with the third or fourth hydration shell water molecules. As these types of interaction are not characteristic of either a hard or soft metal cation, we introduced the concept of a "pseudo-soft" cation to define Au+ behavior. Then, extending the study, we performed the same computations replacing Au+ with Hg2+, an isoelectronic cation. If Hg2+ behaves like Au+ for small water clusters, a topological, geometrical, and energetical transition appears when the number of water molecules increases. Regarding the HSAB theory, this transition is characteristic of a shift of Hg2+ from a pseudosoft form to a soft ion and appears to be due to a competition between the relativistic and correlation effects. Indeed, if relativistic effects are predominant, then mercury will behave like gold and have a similar subvalence/geometry; otherwise when correlation effects are predominant, Hg2+ behaves like a soft cation.

Organohelium compounds: structures, stabilities and chemical bonding analyses

This paper deals with the possibility of forming short and relatively strong carbon-helium bonds in small typical organic molecules through substitution of one or several H atoms by He(+). A structural and energetics study (based on high-level calculations) of this unusual bonding, as well as a topological characterization of the resulting cations, is undertaken. Stable species generally requires substitution of about half of the hydrogen atoms for formation. Under these conditions, the number of such species appears to be potentially unlimited. "True" C-He bonds exhibit equilibrium distances ranging from 1.327 (C2H2He2(2+)) to 1.129 Å (He2CO(2+)). The energies of neutral He releasing range from approximately 5 kcal mol(-1) [He2CO(2+), (Z)-C2H2He2(2+)] to 25 kcal mol(-1) (C2HHe3(3+)), but remain most frequently around 10 kcal mol(-1). However, most of He(+)-substituted hydrocarbons are metastable with respect to C-C cleavage, except derivatives of ethene. Atoms in molecules (AIM) and electron localization function (ELF) topological descriptors classify the C-He bond as a weak charge-shift interaction [S. Shaik, D. Danovich, B. Silvi, D. L. Lauvergnat, P. C. Hiberty, Chem. Eur. J. 2005, 11, 6358-6371] in agreement with a recent publication by Rzepa [S. H. Rzepa, Nat. Chem. 2010, 2, 390-393]. He2CO(2+) is the only investigated compound that presents a C-He bonding ELF basin, which indicates a non-negligible covalent contribution to the bond. Other modifications in the electronic structure, such as the breaking of the triple bond in ethyne derivatives or the loss of aromaticity in C6H3He3(3+), are also nicely revealed by the ELF topology.

Pseudobond parameters for QM/MM studies involving nucleosides, nucleotides, and their analogs

In biological systems involving nucleosides, nucleotides, or their respective analogs, the ribose sugar moiety is the most common reaction site, for example, during DNA replication and repair. However, nucleic bases, which comprise a sizable portion of nucleotide molecules, are usually unreactive during such processes. In quantum mechanical∕molecular simulations of nucleic acid reactivity, it may therefore be advantageous to describe specific ribosyl or ribosyl phosphate groups quantum mechanically and their respective nucleic bases with a molecular mechanics potential function. Here, we have extended the pseudobond approach to enable quantum mechanical∕molecular mechanical simulations involving nucleotides, nucleosides, and their analogs in which the interface between the two subsystems is located between the sugar and the base, namely, the C(sp(3))-N(sp(2)) bond. The pseudobond parameters were optimized on a training set of 10 molecules representing several nucleotide and nucleoside bases and analogs, and they were then tested on a larger test set of 20 diverse molecules. Particular emphasis was placed on providing accurate geometries and electrostatic properties, including electrostatic potential, natural bond orbital (NBO) and atoms in molecules (AIM) charges and AIM first moments. We also tested the optimized parameters on five nucleotide and nucleoside analogues of pharmaceutical relevance and a small polypeptide (triglycine). Accuracy was maintained for these systems, which highlights the generality and transferability of the pseudobond approach.

Ab initio QM/MM calculations show an intersystem crossing in the hydrogen abstraction step in dealkylation catalyzed by AlkB

AlkB is a bacterial enzyme that catalyzes the dealkylation of alkylated DNA bases. The rate-limiting step is known to be the abstraction of an H atom from the alkyl group on the damaged base by a Fe(IV)-oxo species in the active site. We have used hybrid ab initio quantum mechanical/molecular mechanical methods to study this step in AlkB. Instead of forming an Fe(III)-oxyl radical from Fe(IV)-oxo near the C-H activation transition state, the reactant is found to be an Fe(III)-oxyl with an intermediate-spin Fe (S = 3/2) ferromagnetically coupled to the oxyl radical, which we explore in detail using molecular orbital and quantum topological analyses. The minimum energy pathway remains on the quintet surface, but there is a transition between (IS)Fe(III)-oxyl and the state with a high-spin Fe (S = 5/2) antiferromagnetically coupled to the oxyl radical. These findings provide clarity for the evolution of the well-known π and σ channels on the quintet surface in the enzyme environment. Additionally, an energy decomposition analysis reveals nine catalytically important residues for the C-H activation step, some of which are conserved in two human homologues. These conserved residues are proposed as targets for experimental mutagenesis studies.

Toward a Deeper Understanding of Enzyme Reactions Using the Coupled ELF/NCI Analysis: Application to DNA Repair Enzymes

The combined Electron Localization Funtion (ELF)/ Noncovalent Interaction (NCI) topological analysis (Gillet et al. J. Chem. Theory Comput.2012, 8, 3993) has been extended to enzymatic reaction paths. We applied ELF/NCI to the reactions of DNA polymerase λ and the ε subunit of DNA polymerase III. ELF/NCI is shown to provide insights on the interactions during the evolution of enzymatic reactions including predicting the location of TS from structures located earlier along the reaction coordinate, differential metal coordination, and on barrier differences with two different cations.

Coupling quantum interpretative techniques: another look at chemical mechanisms in organic reactions

A cross ELF-NCI analysis is tested over prototypical organic reactions. The synergetic use of ELF and NCI enables the understanding of reaction mechanisms since each method can respectively identify regions of strong and weak electron pairing. Chemically intuitive results are recovered and enriched by the identification of new features. Non covalent interactions are found to foresee the evolution of the reaction from the initial steps. Within NCI, no topological catastrophe is observed as changes are continuous to such an extent that future reaction steps can be predicted from the evolution of the initial NCI critical points. Indeed, strong convergences through the reaction paths between ELF and NCI critical points enable to identify key interactions at the origin of the bond formation. VMD scripts enabling the automatic generation of movies depicting the cross NCI/ELF analysis along a reaction path (or following a Born-Oppenheimer molecular dynamics trajectory) are provided as S.I.

Introducing the ELF Topological Analysis in the Field of Quasirelativistic Quantum Calculations

We present an original formulation of the electron localization function (ELF) in the field of relativistic two-component DFT calculations. Using I2 and At2 species as a test set, we show that the ELF analysis is suitable to evaluate the spin-orbit effects on the electronic structure. Beyond these examples, this approach opens up new opportunities for the bonding analysis of large molecular systems involving heavy and superheavy elements.

Evolution of amide stacking in larger γ-peptides: triamide H-bonded cycles

The single-conformation spectroscopy of two model γ-peptides has been studied under jet-cooled conditions in the gas phase. The methyl-capped triamides, Ac-γ(2)-hPhe-γ(2)-hAla-NHMe and Ac-γ(2)-hAla-γ(2)-hPhe-NHMe, were probed by resonant two-photon ionization (R2PI) and resonant ion-dip infrared (RIDIR) spectroscopies. Four conformers of Ac-γ(2)-hPhe-γ(2)-hAla-NHMe and three of Ac-γ(2)-hAla-γ(2)-hPhe-NHMe were observed and spectroscopically interrogated. On the basis of comparison with the predictions of density functional theory calculations employing a dispersion-corrected functional (ωB97X-D/6-311++G(d,p)), all seven conformers have been assigned to particular conformational families. The preference for formation of nine-membered rings (C9) observed in a previous study [James, W. H., III et al., J. Am. Chem. Soc. 2009, 131, 14243] of the smaller analog, Ac-γ(2)-hPhe-NHMe, carries over to these triamides, with four of the seven conformers forming C9/C9 sequential double-ring structures, and one conformer a C9/C14 bifurcated double ring. The remaining two conformers form C7/C7/C14 H-bonded cycles involving all three amide NH groups, unprecedented in other peptides and peptidomimetics. The amide groups in these structures form a H-bonded triangle with the two trimethylene bridges forming loops above and below the molecule's midsection. The structure is a natural extension of amide stacking, with the two terminal amides blocked from forming the amide tristack by formation of the C14 H-bond. Pair interaction energy decomposition analysis based on the fragment molecular orbital method (FMO-PIEDA) is used to determine the nonbonded contributions to the stabilization of these conformers. Natural bond orbital (NBO) analysis identifies amide stacking with a pair of n → π* interactions between the nitrogen lone pairs and π* orbitals on the carbonyl of the opposing amide groups.

New insights in quantum chemical topology studies using numerical grid-based analyses

New insights in Quantum Chemical Topology of one-electron density functions have been proposed here by using a recent grid-based algorithm (Tang et al., J Phys Condens Matter 2009, 21, 084204), initially designed for the decomposition of the electron density. Beyond the charge analysis, we show that this algorithm is suitable for different scalar functions showing a more complex topology, that is, the Laplacian of the electron density, the electron localization function (ELF), and the molecular electrostatic potential (MEP). This algorithm makes use of a robust methodology enabling to numerically assign the data points of three-dimensional grids to basin volumes, and it has the advantage of requiring only the values of the scalar function without details on the wave function used to build the grid. Our implementation is briefly outlined (program named TopChem), its capabilities are examined, and technical aspects in terms of CPU requirement and accuracy of the results are discussed. Illustrative examples for individual molecules and crystalline solids obtained with gaussian and plane-wave-based density functional theory calculations are presented. Special attention was given to the MEP because its topological analysis is complex and scarce.

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.

Study of the docking of competitive inhibitors at a model of tyrosinase active site: insights from joint broken-symmetry/Spin-Flip DFT computations and ELF topological analysis

Following our previous study (Piquemal et al., New J. Chem., 2003, 27, 909), we present here a DFT study of the inhibition of the Tyrosinase enzyme. Broken-symmetry DFT computations are supplemented with Spin-Flip TD-DFT calculations, which, for the first time, are applied to such a dicopper enzyme. The chosen biomimetic model encompasses a dioxygen molecule, two Cu(II) cations, and six imidazole rings. The docking energy of a natural substrate, namely phenolate, together with those of several inhibitor and non-inhibitor compounds, are reported and show the ability of the model to rank the most potent inhibitors in agreement with experimental data. With respect to broken-symmetry calculations, the Spin-Flip TD-DFT approach reinforces the possibility for theory to point out potent inhibitors: the need for the deprotonation of the substrates, natural or inhibitors, is now clearly established. Moreover, Electron Localization Function (ELF) topological analysis computations are used to deeply track the particular electronic distribution of the Cu-O-Cu three-center bonds involved in the enzymatic Cu(2)O(2) metallic core (Piquemal and Pilmé, J. Mol. Struct.: Theochem, 2006, 77, 764). It is shown that such bonds exhibit very resilient out-of-plane density expansions that play a key role in docking interactions: their 3D-orientation could be the topological electronic signature of oxygen activation within such systems.

Polarizable molecular dynamics simulation of Zn(II) in water using the AMOEBA force field

The hydration free energy, structure, and dynamics of the zinc divalent cation are studied using a polarizable force field in molecular dynamics simulations. Parameters for the Zn(2+) are derived from gas-phase ab initio calculation of Zn(2+)-water dimer. The Thole-based dipole polarization is adjusted based on the Constrained Space Orbital Variations (CSOV) calculation while the Symmetry Adapted Perturbation Theory (SAPT) approach is also discussed. The vdW parameters of Zn(2+) have been obtained by comparing the AMOEBA Zn(2+)-water dimerization energy with results from several theory levels and basis sets over a range of distances. Molecular dynamics simulations of Zn(2+) solvation in bulk water are subsequently performed with the polarizable force field. The calculated first-shell water coordination number, water residence time and free energy of hydration are consistent with experimental and previous theoretical values. The study is supplemented with extensive Reduced Variational Space (RVS) and Electron Localization Function (ELF) computations in order to unravel the nature of the bonding in Zn(2+)(H(2)O)(n) (n=1,6) complexes and to analyze the charge transfer contribution to the complexes. Results show that the importance of charge transfer decreases as the size of Zn-water cluster grows due to anticooperativity and to changes in the nature of the metal-ligand bonds. Induction could be dominated by polarization when the system approaches condensed-phase and the covelant effects are eliminated from the Zn(II)-water interaction. To construct an "effective" classical polarizable potential for Zn(2+) in bulk water, one should therefore avoid over-fitting to the ab initio charge transfer energy of Zn(2+)-water dimer. Indeed, in order to avoid overestimation of condensed-phase many-body effects, which is crucial to the transferability of polarizable molecular dynamics, charge transfer should not be included within the classical polarization contribution and should preferably be either incorporated in to the pairwise van der Waals contribution or treated explicitly.

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.