Pierre Mignon

Influence of the π–π interaction on the hydrogen bonding capacity of stacked DNA/RNA bases

The interplay between aromatic stacking and hydrogen bonding in nucleobases has been investigated via high-level quantum chemical calculations. The experimentally observed stacking arrangement between consecutive bases in DNA and RNA/DNA double helices is shown to enhance their hydrogen bonding ability as opposed to gas phase optimized complexes. This phenomenon results from more repulsive electrostatic interactions as is demonstrated in a model system of cytosine stacked offset-parallel with substituted benzenes. Therefore, the H-bonding capacity of the N3 and O2 atoms of cytosine increases linearly with the electrostatic repulsion between the stacked rings. The local hardness, a density functional theory-based reactivity descriptor, appears to be a key index associated with the molecular electrostatic potential (MEP) minima around H-bond accepting atoms, and is inversely proportional to the electrostatic interaction between stacked molecules. Finally, the MEP minima on surfaces around the bases in experimental structures of DNA and RNA–DNA double helices show that their hydrogen bonding capacity increases when taking more neighboring (intra-strand) stacking partners into account.

Cooperative effects at water–crystalline silica interfaces strengthen surface silanol hydrogen bonding. An ab initio molecular dynamics study

Silica and silica based materials are widely used in chemistry and materials science due to their importance in many technological fields. The properties of these materials, which are crucial for their applications, are mainly determined by the presence of hydrogen bonding between surface silanols. Here, we present ab initio molecular dynamics simulations (AIMD) on different surfaces derived from the crystallographic α-quartz (100) and the α-cristobalite (001) and (101) faces, both free and at the interface with liquid water. The focus was on studying whether water adsorption can disrupt the H-bond pattern at the pristine free silica surface and how deep the perturbation due to the contact with the surface affects the structure of the water multilayer. Results highlight that the water phase is over structured at the interface with silica, as compared to water bulk. Furthermore, an apparent counterintuitive behavior has been observed for quartz (100) and cristobalite (001) surfaces: the interaction with water does not cleave the pre-existent H-bonds between the surface silanol groups. On the contrary, in several cases, it is observed that SiOH···OHSi H-bonds are even strengthened, as the result of a mutual cooperative H-donor/H-acceptor enhancement between silanols and water molecules, which may alter the adsorption capability of these silica surfaces.

Understanding the Reactivity and Basicity of Zeolites: A Periodic DFT Study of the Disproportionation of N2O4 on Alkali‐Cation‐Exchanged Zeolite Y

The disproportionation of N(2)O(4) into NO(3)(-) and NO(+) on Y zeolites has been studied through periodic DFT calculations to unravel 1) the role of metal cations and the framework oxygen atoms and 2) the relationship between the NO(+) stretching frequency and the basicity of zeolites. We have considered three situations: adsorption on site II cations with and without a cation at site III and adsorption on a site III cation. We observed that cations at sites II and III cooperate to stabilize N(2)O(4) and that the presence of a cation at site III is necessary to allow the disproportionation reaction. The strength of the stabilization is due to the number of stabilizing interactions increasing with the size of the cation and to the Lewis acidity of the alkali cations, which increases as the size of the cations decreases. In the product, NO(3)(-) interacts mainly with the cations and NO(+) with the basic oxygen atoms of the tetrahedral aluminium through its nitrogen atom. As the cation size increases, the NO(3)(-)...cation interaction increases. As a result, the negative charge of the framework is less well screened by the larger cations and the interaction between NO(+) and the basic oxygen atoms becomes stronger. NO(+) appears to be a good probe of zeolite basicity, in agreement with experimental observations.

Theoretical study of the adsorption of DNA bases on the acidic external surface of montmorillonite.

In the present study, DFT periodic plane wave calculations, at the PBE-D level of theory, were carried out to investigate the interaction of DNA nucleobases with acidic montmorillonite. The surface model was considered in its octahedral (Osub) and tetrahedral (Tsub) substituted forms, known to have different acidic properties. The adsorption of adenine, guanine and cytosine was considered in both orthogonal and coplanar orientations with the surface, interacting with the proton via a given heteroatom. In almost all considered cases, adsorption involved the spontaneous proton transfer to the nucleobase, with a more pronounced character in the Osub structures. The binding energy is about 10 kcal mol(-1) larger for Osub than for Tsub complexes mainly due to the larger acidity in Osub surfaces and due to the better stabilization by H-bond contacts between the negatively charged surface and the protonated base. The binding energy of coplanar orientations of the base is observed to be as large as the orthogonal ones due to a balance between electrostatic and dispersion contributions. Finally the binding of guanine and adenine on the acidic surface amounts to 50 kcal mol(-1) while that of cytosine rises to 44 kcal mol(-1).

DFT study of the fragmentation mechanism of uracil RNA base

The fragmentation process of the uracil RNA base has been investigated via DFT calculations in order to assign fragments to the ionisation mass spectrum obtained after dissociation induced by collision experiments. The analysis of the electronic distribution and geometry parameters of the cation allows selection of several bonds that may be cleaved and lead to the formation of various fragments. Differences are observed in the electronic behaviour of the bond breaking as well as the energy required for the cleavage. It is reported that N(3)-C(4) and N(1)-C(2) bonds are more easily cleaved than the C(5)-C(6) bond, since the corresponding energy barriers amount to ΔG = +1.627, +1.710, +5.459 eV, respectively, which makes the C(5)-C(6) bond cleavage almost prohibited. Among all possible formed fragments, the formation of the OCN(+) fragment for the peak at m/z = 42 Da is excluded because of an intermediate that was not observed experimentally and too a large free energy barrier. Based on the required free energy, it is observed that two fragment derivatives: C(2)H(4)N(+) and C(2)H(2)O˙(+) may be formed, with a small preference for C(2)H(4)N(+). This latter product is not formed through a retro Diels Alder reaction in contrast to C(2)H(2)O˙(+). The following sequence is proposed for the peak at 42 Da: C(2)H(4)N(+) (from N(1)-C(2), C(4)-C(5) cleavages) > C(2)H(2)O˙(+) (from N(3)-C(4), N(1)-C(2) and C(5)-C(6) cleavages) > C(2)H(4)N(+) (from N(1)-C(2), N(3)-C(4) and C(4)-C(5)) > C(2)H(2)O˙(+) (from C(5)-C(6), N(1)-C(2) and N(3)-C(4) cleavages) > NCO(+) (from N(1)-C(2), C(4)-C(5) and N(3)-C(4) cleavages). Finally the peak at 28 Da is assigned to CNH(2)(+) derivatives that can be formed through two different paths, the easiest one requiring 5.4 eV.

Dual descriptor and molecular electrostatic potential: complementary tools for the study of the coordination chemistry of ambiphilic ligands

In this paper, we show that the ambiphilic properties of some organic ligands in organometallic complexes may be retrieved readily from simple calculations in the framework of conceptual density functional theory (C-DFT): namely, the dual descriptor (DD) and the molecular electrostatic potential (MEP) of the ligands afford a rather straightforward interpretation of experimental trends such as the bonding geometry and the electronic properties of complexes in terms of σ-, π- and back-bonding. The studied ligands were chosen to be representative of the wide variety organometallic chemistry offers, ranging from neutral to charged systems and from diatomic to polyatomic molecules. The present approach is general since all relevant parameters are retrieved from the electron density, obtained either from a DFT or post-Hartree-Fock calculation. It is believed to be helpful for organometallic chemists, since it allows a deep understanding and may be used as a predictive tool of the coordinating properties of ligands.

Quantitative structure–activity relationship to predict acute fish toxicity of organic solvents

REACH regulation requires ecotoxicological data to characterize industrial chemicals. To limit in vivo testing, Quantitative Structure-Activity Relationships (QSARs) are advocated to predict toxicity of a molecule. In this context, the topic of this work was to develop a reliable QSAR explaining the experimental acute toxicity of organic solvents for fish trophic level. Toxicity was expressed as log(LC50), the concentration in mmol.L(-1) producing the 50% death of fish. The 141 chemically heterogeneous solvents of the dataset were described by physico-chemical descriptors and quantum theoretical parameters calculated via Density Functional Theory. The best subsets of solvent descriptors for LC50 prediction were chosen both through the Kubinyi function associated with Enhanced Replacement Method and a stepwise forward multiple linear regressions. The 4-parameters selected in the model were the octanol-water partition coefficient, LUMO energy, dielectric constant and surface tension. The predictive power and robustness of the QSAR developed were assessed by internal and external validations. Several techniques for training sets selection were evaluated: a random selection, a LC50-based selection, a balanced selection in terms of toxic and non-toxic solvents, a solvent profile-based selection with a space filling technique and a D-optimality onions-based selection. A comparison with fish LC50 predicted by ECOSAR model validated for neutral organics confirmed the interest of the QSAR developed for the prediction of organic solvent aquatic toxicity regardless of the mechanism of toxic action involved.

Study of Prepolymerization Complex Formation in the Synthesis of Steroid-Based Molecularly Imprinted Polymers

The preparation of steroid-based molecularly imprinted polymers (MIPs) based upon noncovalent interaction is particularly suited for selective capture of steroid hormones in biological and environmental samples. The success of this method lies in the optimization of the interaction between steroids (template) and methacrylic acid (functional monomer) in the prepolymerization mixture. NMR techniques coupled with DFT calculations were used to evaluate the capacity of the methacrylic acid to bind a steroid for future applications. The androstane derivative steroids considered in the present study have two functional groups at C(3) and C(17), which may interact with the methacrylic acid. These can be hydroxyl or ketone groups. Experimental results show that the steroids can be divided in three groups corresponding to the ketone type at C(3), the H-bond strength increasing with the number of double bonds. DFT calculations are in very good agreement with experimental results, showing increasing binding energy from no bonds, a single bond, and two double bond steroids. For steroids holding a hydroxyl group the binding energy obtained in the solvent model is comparable to the binding energy of single bond ketone steroids.

Oxo iron(IV) as an oxidative active intermediate of p-chlorophenol in the Fenton reaction: a DFT study

Debate continues over which active species plays the role of oxidative agent during the Fenton reaction-the HO˙ radical or oxo iron [Fe(IV)O](2+). In this context, the present study investigates the oxidation of p-chlorophenol by [Fe(IV)O(H(2)O)(5)](2+) using DFT calculations, within gas-phase and micro-solvated models, in order to explore the possible role of oxo iron as a reactant. The results show that the chlorine atom substitution of p-chlorophenol by oxo iron is a highly stabilising step (ΔH = -83 kcal mol(-1)) with a free energy barrier of 5.8 kcal mol(-1) in the micro-solvated model. This illustrates the high oxidising power of the [Fe(IV)O(H(2)O)(5)](2+) complex. On the other hand, the breaking of the Fe-O bond, leading to the formation of hydroquinone, is observed to be the rate-determining step of the reaction. The rather large free energy barrier corresponding to this bond cleavage amounts to 10.2 and 9.3 kcal mol(-1) in the gas-phase and micro-solvated models, respectively. Elsewhere, the lifetime of the HO˙ radical has previously been shown to be extremely small. These facts, combined with observations of oxo iron under certain experimental conditions, suggest that oxo iron is a highly plausible oxidative species of the reaction. In addition, a trigonal bipyramidal iron complex, coordinated either by hydroxyl groups and/or by water molecules, has been found in all described mechanisms. This structure appears to be a stable intermediate; and to our knowledge, it has not been characterised by previous studies.

Hydrogen-Induced Adsorption of Carbon Monoxide on the Gold Dimer Cation: A Joint Experimental and DFT Investigation

It is demonstrated, using tandem mass spectrometry and radio frequency ion trap, that the adsorption of a H atom on the gold dimer cation, Au2H+, prevents its dissociation and allows for adsorption of CO. Reaction kinetics are measured by employing a radio frequency ion trap, where Au2+ and CO interact for a given reaction time. The effect of a hydrogen atom is evaluated by comparing reaction rate constants measured for Au2+ and Au2H+. The theoretical results for the adsorption of CO molecules and their reaction characteristics with Au2+ and Au2H+ are found to agree with the experimental findings. The joint investigations provide insights into hydrogen atom adsorption effects and consequent reaction mechanisms.