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

On the Quantum Chemical Nature of Lead(II) “Lone Pair”

We study the quantum chemical nature of the Lead(II) valence basins, sometimes called the lead “lone pair”. Using various chemical interpretation tools, such as molecular orbital analysis, natural bond orbitals (NBO), natural population analysis (NPA) and electron localization function (ELF) topological analysis, we study a variety of Lead(II) complexes. A careful analysis of the results shows that the optimal structures of the lead complexes are only governed by the 6s and 6p subshells, whereas no involvement of the 5d orbitals is found. Similarly, we do not find any significant contribution of the 6d. Therefore, the Pb(II) complexation with its ligand can be explained through the interaction of the 6s² electrons and the accepting 6p orbitals. We detail the potential structural and dynamical consequences of such electronic structure organization of the Pb (II) valence domain.

Exploitation of 2D Cu-MOF nanosheets as a unique electroactive material for ultrasensitive Cu(II) ion estimation in various real samples

Herein, hybrid carbon sensor has been developed with graphite sheets as a matrix, tricresyl phosphate (TCP) as a plasticizer and nanosheets of 2D Cu-MOF (metal-organic framework) as an electroactive material for the ultrasensitive Cu(II) ion detection in various real samples. Where, the present study proves the efficiency of 2D Cu-MOF as a promising sensing material for the development of Cu(II) ion selective carbon sensor. The developed 2D Cu-MOF nanosheets based sensor containing 2D Cu-MOF: TCP: graphite in the ratio of 2.67: 30.54: 66.79 (% wt/wt) displayed unique Nernstian behavior over two linearity ranges of 1.0 × 10^-11-1.0 × 10^-9 and 1.0 × 10^-5-1.0 × 10^-1 mol L^-1 with slopes of 29.5 ± 0.25 and 29.6 ± 0.13 mV decade^-1, respectively. The fabricated carbon sensor achieved a widely pH independency, fast response time and superior thermal stability with highly selective and ultrasensitive performance. Moreover, It has been efficiently applied for the Cu(II) ion potentiometric estimation in human hair, sesames seeds, two different tea infusions and tap water real samples.

Interfacial Water Many-Body Effects Drive Structural Dynamics and Allosteric Interactions in SARS-CoV-2 Main Protease Dimerization Interface

Following our previous work ( Chem. Sci. 2021, 12, 4889-4907), we study the structural dynamics of the SARS-CoV-2 Main Protease dimerization interface (apo dimer) by means of microsecond adaptive sampling molecular dynamics simulations (50 μs) using the AMOEBA polarizable force field (PFF). This interface is structured by a complex H-bond network that is stable only at physiological pH. Structural correlations analysis between its residues and the catalytic site confirms the presence of a buried allosteric site. However, noticeable differences in allosteric connectivity are observed between PFFs and non-PFFs. Interfacial polarizable water molecules are shown to appear at the heart of this discrepancy because they are connected to the global interface H-bond network and able to adapt their dipole moment (and dynamics) to their diverse local physicochemical microenvironments. The water-interface many-body interactions appear to drive the interface volume fluctuations and to therefore mediate the allosteric interactions with the catalytic cavity.

Sulfamethoxazole Removal from Drinking Water by Activated Carbon: Kinetics and Diffusion Process

Sulfamethoxazole (SMX), a pharmaceutical residue, which is persistent and mobile in soils, shows low biodegradability, and is frequently found in the different aquatic compartments, can be found at very low concentrations in water intended for human consumption. In conditions compatible with industrial practices, the kinetic reactivity and performance of tap water purification using activated carbon powder (ACP) are examined here using two extreme mass ratios of SMX to ACP: 2 µg/L and 2 mg/L of SMX for only 10 mg/L of ACP. In response to surface chemistry, ACP texture and the intrinsic properties of SMX in water at a pH of 8.1, four kinetic models, and two monosolute equilibrium models showed a total purification of the 2 µg/L of SMX, the presence of energetic heterogeneity of surface adsorption of ACP, rapid kinetics compatible with the residence times of industrial water treatment processes, and kinetics affected by intraparticle diffusion. The adsorption mechanisms proposed are physical mechanisms based mainly on π-π dispersion interactions and electrostatic interactions by SMX-/Divalent cation/ArO- and SMX-/Divalent cation/ArCOO- bridging. Adsorption in tap water, also an innovative element of this study, shows that ACP is very efficient for the purification of very slightly polluted water.

A Rare Natural Benzo[k,l]xanthene as a Turn-Off Fluorescent Sensor for Cu2+ Ion

Rapid and efficient analyses of copper ions are crucial to providing key information for Cu²⁺ in living cells because of their biological importance. In this study, we reported one new turn-off fluorescent sensor for Cu²⁺ with a benzo[k,l]xanthene core, which served as an efficient cation sensor for copper ion over a wide range of other cations (Na⁺, K⁺, Ag⁺, Hg²⁺, Cd²⁺, Co²⁺, Ni²⁺, Zn²⁺, Mg²⁺, and Fe³⁺) owing to the catechol group in the aromatic core. The sensor showed selectivity for Cu²⁺ over other ions; the logK_β for Cu²⁺ binding to compound 1 had a value of 13.265. In the presence of Cu²⁺, sensor 1 provided significant fluorescence decrement; Co²⁺, and Ni²⁺ caused a fluorescence decrement when employed at a higher concentration than Cu²⁺, while Na⁺, K⁺, Hg²⁺, Cd²⁺, Zn²⁺, and Mg²⁺ metal ions produced only minor changes in fluorescence intensity. Fluorescence experiments demonstrate that compound 1 may have an application as a fluorescent probe for detecting Cu²⁺ with a limit of detection of 0.574 µM.

Curly arrows, electron flow, and reaction mechanisms from the perspective of the bonding evolution theory

Despite the usefulness of curly arrows in chemistry, their relationship with real electron density flows is still imprecise, and even their direct connection to quantum chemistry is still controversial. The paradigmatic description - from first principles - of the mechanistic aspects of a given chemical process is based mainly on the relative energies and geometrical changes at the stationary points of the potential energy surface along the reaction pathway; however, it is not sufficient to describe chemical systems in terms of bonding aspects. Probing the electron density distribution during a chemical reaction can provide important insights, enabling us to understand and control chemical reactions. This aim has required an extension of the relationships between the concepts of traditional chemistry and those of quantum mechanics. Bonding evolution theory (BET), which combines the topological analysis of the electron localization function (ELF) and Thom's catastrophe theory (CT), provides a powerful method that offers insight into the molecular mechanism of chemical rearrangements. In agreement with the laws of physical and aspects of quantum theory, BET can be considered an appropriate tool to tackle chemical reactivity with a wide range of possible applications. In this work, BET is applied to address a long-standing problem: the ability to monitor the flow of electron density. BET analysis shows a connection between quantum mechanics and bond making/forming processes. Likewise, the present approach retrieves the classical curly arrows used to describe the rearrangements of chemical bonds and provides detailed physical grounds for this type of representation. We demonstrate this procedure using the test set of prototypical examples of thermal ring apertures, and the degenerated Cope rearrangement of semibullvalene.

Curly arrows meet electron density transfers in chemical reaction mechanisms: from electron localization function (ELF) analysis to valence-shell electron-pair repulsion (VSEPR) inspired interpretation

The displacement of the nuclei along the reaction path provides an explanatory interpretation of the electron density transfers making possible to understand chemical reactions.

Structure-Function Relationship of the Interaction between Tissue Factor and Factor VIIa

Interactions between tissue factor and factor VIIa are the primary initiators of coagulation in hemostasis and certain thrombotic diseases. Tissue factor, an integral membrane protein expressed extensively outside of the vasculature, is the regulatory protein cofactor for coagulation factor VIIa. Factor VIIa, a trypsin-like serine protease homologous with other blood coagulation proteases, is weakly active when free in solution and must bind its membrane-bound cofactor for physiologically relevant activity. Tissue factor allosterically activates factor VIIa by several mechanisms such as active site positioning, spatial stabilization, and direct interactions with the substrate. Protein-membrane interactions between tissue factor, factor VIIa, and substrates all play critical roles in modulating the activity of this enzyme complex. Additionally, divalent cations such as Ca(2+) and Mg(2+) are critical for correct protein folding, as well as protein-membrane and protein-protein interactions. The contributions of these factors toward tissue factor-factor VIIa activity are discussed in this review.

Quantifying vitamin K-dependent holoprotein compaction caused by differential γ-carboxylation using high-pressure size exclusion chromatography

This study uses high-pressure size exclusion chromatography (HPSEC) to quantify divalent metal ion (X(2+))-induced compaction found in vitamin K-dependent (VKD) proteins. Multiple X(2+) binding sites formed by the presence of up to 12 γ-carboxyglutamic acid (Gla) residues are present in plasma-derived FIX (pd-FIX) and recombinant FIX (r-FIX). Analytical ultracentrifugation (AUC) was used to calibrate the Stokes radius (R) measured by HPSEC. A compaction of pd-FIX caused by the filling of Ca(2+) and Mg(2+) binding sites resulted in a 5 to 6% decrease in radius of hydration as observed by HPSEC. The filling of Ca(2+) sites resulted in greater compaction than for Mg(2+) alone where this effect was additive or greater when both ions were present at physiological levels. Less X(2+)-induced compaction was observed in r-FIX with lower Gla content populations, which enabled the separation of biologically active r-FIX species from inactive ones by HPSEC. HPSEC was sensitive to R changes of approximately 0.01nm that enabled the detection of FIX compaction that was likely cooperative in nature between lower avidity X(2+) sites of the Gla domain and higher avidity X(2+) sites of the epidermal growth factor 1 (EGF1)-like domain.

Mysteries of metals in metalloenzymes

Natural metalloenzymes are often the most proficient catalysts in terms of their activity, selectivity, and ability to operate at mild conditions. However, metalloenzymes are occasionally surprising in their selection of catalytic metals, and in their responses to metal substitution. Indeed, from the isolated standpoint of producing the best catalyst, a chemist designing from first-principles would likely choose a different metal. For example, some enzymes employ a redox active metal where a simple Lewis acid is needed. Such are several hydrolases. In other cases, substitution of a non-native metal leads to radical improvements in reactivity. For example, histone deacetylase 8 naturally operates with Zn(2+) in the active site but becomes much more active with Fe(2+). For β-lactamases, the replacement of the native Zn(2+) with Ni(2+) was suggested to lead to higher activity as predicted computationally. There are also intriguing cases, such as Fe(2+)- and Mn(2+)-dependent ribonucleotide reductases and W(4+)- and Mo(4+)-dependent DMSO reductases, where organisms manage to circumvent the scarcity of one metal (e.g., Fe(2+)) by creating protein structures that utilize another metal (e.g., Mn(2+)) for the catalysis of the same reaction. Naturally, even though both metal forms are active, one of the metals is preferred in every-day life, and the other metal variant remains dormant until an emergency strikes in the cell. These examples lead to certain questions. When are catalytic metals selected purely for electronic or structural reasons, implying that enzymatic catalysis is optimized to its maximum? When are metal selections a manifestation of competing evolutionary pressures, where choices are dictated not just by catalytic efficiency but also by other factors in the cell? In other words, how can enzymes be improved as catalysts merely through the use of common biological building blocks available to cells? Addressing these questions is highly relevant to the enzyme design community, where the goal is to prepare maximally efficient quasi-natural enzymes for the catalysis of reactions that interest humankind. Due to competing evolutionary pressures, many natural enzymes may not have evolved to be ideal catalysts and can be improved for the isolated purpose of catalysis in vitro when the competing factors are removed. The goal of this Account is not to cover all the possible stories but rather to highlight how variable enzymatic catalysis can be. We want to bring up possible factors affecting the evolution of enzyme structure, and the large- and intermediate-scale structural and electronic effects that metals can induce in the protein, and most importantly, the opportunities for optimization of these enzymes for catalysis in vitro.

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.

DFT and docking studies of rhodostreptomycins A and B and their interactions with solvated/nonsolvated Mg²⁺ and Ca²⁺ ions

The interactions of L-aminoglucosidic stereoisomers such as rhodostreptomycins A (Rho A) and B (Rho B) with cations (Mg(2+), Ca(2+), and H(+)) were studied by a quantum mechanical method that utilized DFT with B3LYP/6-311G. Docking studies were also carried out in order to explore the surface recognition properties of L-aminoglucoside with respect to Mg(2+) and Ca(2+) ions under solvated and nonsolvated conditions. Although both of the stereoisomers possess similar physicochemical/antibiotic properties against Helicobacter pylori, the thermochemical values for these complexes showed that its high affinity for Mg(2+) cations caused the hydration of Rho B. According to the results of the calculations, for Rho A-Ca(2+)(H2O)6, ΔH = -72.21 kcal mol(-1); for Rho B-Ca(2+)(H2O)6, ΔH = -72.53 kcal mol(-1); for Rho A-Mg(2+)(H2O)6, ΔH = -72.99 kcal mol(-1) and for Rho B-Mg(2+)(H2O)6, ΔH = -95.00 kcal mol(-1), confirming that Rho B binds most strongly with hydrated Mg(2+), considering the energy associated with this binding process. This result suggests that Rho B forms a more stable complex than its isomer does with magnesium ion. Docking results show that both of these rhodostreptomycin molecules bind to solvated Ca(2+) or Mg(2+) through hydrogen bonding. Finally, Rho B is more stable than Rho A when protonation occurs.

Understanding the structure and electronic properties of Th4+-water complexes

We present a systematic quantum chemistry study of the [Th(H2O)n]4+ (n = 1 to 10) complexes to gain insight into their electronic structure and properties: the effect of the ligand distribution on the valence shells of the thorium(IV) ion is studied by means of the electron localization function (ELF) topological analysis. Particular care is given to the study of the mono-aqua complex both at its equilibrium geometry, using various tools such as energy decomposition analyses (EDA), and along its dissociation pathway. Indeed, as several electronic states cross the Th4 +-H2O0 ground state along the minimum energy path, we demonstrate that the diabatic representation implemented in MOLPRO is able to generate reference potential energy surfaces that will lead to the evaluation of diabatic dissociation curves. The calculated diabatic interaction energy curve will allow for a consistent parameterization of new generation force fields dedicated to heavy metals based on quantum chemistry.

Structural and functional studies of γ-carboxyglutamic acid domains of factor VIIa and activated Protein C: role of magnesium at physiological calcium

Crystal structures of factor (F) VIIa/soluble tissue factor (TF), obtained under high Mg(2+) (50mM Mg(2+)/5mM Ca(2+)), have three of seven Ca(2+) sites in the γ-carboxyglutamic acid (Gla) domain replaced by Mg(2+) at positions 1, 4, and 7. We now report structures under low Mg(2+) (2.5mM Mg(2+)/5mM Ca(2+)) as well as under high Ca(2+) (5mM Mg(2+)/45 mM Ca(2+)). Under low Mg(2+), four Ca(2+) and three Mg(2+) occupy the same positions as in high-Mg(2+) structures. Conversely, under low Mg(2+), reexamination of the structure of Gla domain of activated Protein C (APC) complexed with soluble endothelial Protein C receptor (sEPCR) has position 4 occupied by Ca(2+) and positions 1 and 7 by Mg(2+). Nonetheless, in direct binding experiments, Mg(2+) replaced three Ca(2+) sites in the unliganded Protein C or APC. Further, the high-Ca(2+) condition was necessary to replace Mg4 in the FVIIa/soluble TF structure. In biological studies, Mg(2+) enhanced phospholipid binding to FVIIa and APC at physiological Ca(2+). Additionally, Mg(2+) potentiated phospholipid-dependent activations of FIX and FX by FVIIa/TF and inactivation of activated factor V by APC. Since APC and FVIIa bind to sEPCR involving similar interactions, we conclude that under the low-Mg(2+) condition, sEPCR binding to APC-Gla (or FVIIa-Gla) replaces Mg4 by Ca4 with an attendant conformational change in the Gla domain ω-loop. Moreover, since phospholipid and sEPCR bind to FVIIa or APC via the ω-loop, we predict that phospholipid binding also induces the functional Ca4 conformation in this loop. Cumulatively, the data illustrate that Mg(2+) and Ca(2+) act in concert to promote coagulation and anticoagulation.

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.

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.

Structural biology of factor VIIa/tissue factor initiated coagulation

Factor VII (FVII) consists of an N-terminal gamma-carboxyglutamic acid domain followed by two epidermal growth factor-like (EGF1 and EGF2) domains and the C-terminal protease domain. Activation of FVII results in a two-chain FVIIa molecule consisting of a light chain (Gla-EGF1-EGF2 domains) and a heavy chain (protease domain) held together by a single disulfide bond. During coagulation, the complex of tissue factor (TF, a transmembrane glycoprotein) and FVIIa activates factor IX (FIX) and factor X (FX). FVIIa is structurally "zymogen-like" and when bound to TF, it is more "active enzyme-like." FIX and FX share structural homology with FVII. Three structural biology aspects of FVIIa/TF are presented in this review. One, regions in soluble TF (sTF) that interact with FVIIa as well as mapping of Ca2+, Mg2+, Na+ and Zn2+ sites in FVIIa and their functions; two, modeled interactive regions of Gla and EGF1 domains of FXa and FIXa with FVIIa/sTF; and three, incompletely formed oxyanion hole in FVIIa/sTF and its induction by substrate/inhibitor. Finally, an overview of the recognition elements in TF pathway inhibitor is provided.

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.

The reaction mechanism of type I phosphomannose isomerases: new information from inhibition and polarizable molecular mechanics studies

Type I phosphomannose isomerases (PMIs) are zinc-dependent metalloenzymes involved in the reversible isomerization of D-mannose 6-phosphate (M6P) and D-fructose 6-phosphate (F6P). 5-Phospho-D-arabinonohydroxamic acid (5PAH), an inhibitor endowed with nanomolar affinity for yeast (Type I) and Pseudomonas aeruginosa (Type II) PMIs (Roux et al., Biochemistry 2004; 43:2926-2934), strongly inhibits human (Type I) PMI (for which we report an improved expression and purification procedure), as well as Escherichia coli (Type I) PMI. Its K(i) value of 41 nM for human PMI is the lowest value ever reported for an inhibitor of PMI. 5-Phospho-D-arabinonhydrazide, a neutral analogue of the reaction intermediate 1,2-cis-enediol, is about 15 times less efficient at inhibiting both enzymes, in accord with the anionic nature of the postulated high-energy reaction intermediate. Using the polarizable molecular mechanics, sum of interactions between fragments ab initio computed (SIBFA) procedure, computed structures of the complexes between Candida albicans (Type I) PMI and the cyclic substrate β-D-mannopyranose 6-phosphate (β-M6P) and between the enzyme and the high-energy intermediate analogue inhibitor 5PAH are reported. Their analysis allows us to identify clearly the nature of each individual active site amino acid and to formulate a hypothesis for the overall mechanism of the reaction catalyzed by Type I PMIs, that is, the ring-opening and isomerization steps, respectively. Following enzyme-catalyzed ring-opening of β-M6P by zinc-coordinated water and Gln111 ligands, Lys136 is identified as the probable catalytic base involved in proton transfer between the two carbon atoms C1 and C2 of the substrate D-mannose 6-phosphate.

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.

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.