Thesis (Ph. D.)--University of North Texas, May, 2009. / Title from title page display. Includes bibliographical references.

Durant cette thèse, nous avons abordé l'étude de la réactivité en phase gazeuse des biomolécules. L’avènement des techniques d’ionisation douces telle que l’ionisation par éléctronébulisation, a rendu possible ces dernières années, la formation d'ions en phase gazeuse sans dégrader la biomolécule étudiée.La Dissociation Induite par Collision (CID) est un cas particulier de spectrométrie de masse en tandem, que nous avons utilisée durant ce travail. Le principe du CID est d'activer les modes rovibrationnelles d’un système moléculaire ionique par collision avec un gaz inerte, ce qui augmente la probabilité de fragmentation de l'ion. Bien qu'étant une technique très utile d'un point de vue analytique, la spectrométrie de masse en tandem ne donne pas d'informations sur les mécanismes des réactions se produisant dans la cellule de collision; afin d’obtenir ces informations, les simulations de dynamique chimiques apparaissent comme un outil satisfaisant. En effet, en utilisant la dynamique directe, nous évitons ainsi d'explorer la totalité de la surface d'énergie potentielle, qui devient compliquée lors de l’étude d’édifices moléculaires de grande taille. Etant donné que les simulations de dynamique chimiques sont limitées à de courtes échelles, de l’ordre de la dizaine de picosecondes, nous avons également employé la théorie unimoléculaire RRKM (Rice-Ramsperger-Kassel-Marcus) pour étudier la réactivité à des temps plus longs, en vue de comprendre les processus réactionnels se produisant à l’issue du processus de relaxation vibrationnelle intramoléculaire (IVR). Durant ce travail de thèse, nous avons choisi d'étudier comme système modèle de base nucléique la molécule d'uracile. Par ailleurs,nous avons aussi étudié la réactivité en phase gazeuse de sucres (cellobiose, maltose et gentiobiose), qui ont été au préalable dérivatisés afin de localiser la charge sur la molécule et ainsi simplifier l’étude théorique associée. / In the present thesis, we address the study of the reactivity of biomolecules in the gasphase.The advent of soft ionization techniques such as electrospray ionization, made possible, in the last years, the gentle formation of ions in the gas phase without breaking the molecule understudy.Collision Induced Dissociation (CID) is aparticular case of tandem mass spectrometrydynamics simulations are pointed like asatisfactory tool. Using direct dynamics weavoid exploring the whole potential energysurface, which becomes really complicatedwhen dealing with big molecules.Since chemical dynamics simulations arerestricted to the short time scale reactivity,typically ~10ps, we make use of the Rice–Ramsperger–Kassel–Marcus (RRKM)unimolecular theory to study the reactivity atUniversité Paris-SaclayEspace Technologique / Immeuble DiscoveryRoute de l’Orme aux Merisiers RD 128 / 91190 Saint-Aubin, Francethat we use in the present thesis. The aim of CIDis to activate the rovibrational modes of an ionicmolecular system by collisions with an inert gas,increasing the probability of the ion of beingfragmented.Despite being a really useful technique, tandemmass spectrometry does not give informationabout the mechanisms of the reactions takingplace in the collision cell; in order to obtain suchinformation, chemicallonger time scales to understand reaction pathsthat take place after intramolecular vibrationrelaxation (IVR).In the present thesis we have chosen to study asmodel system of nucleobase the uracil molecule.Furthermore, we also studied the gas-phase reactivity of carbohydrates (cellobiose, maltose and gentiobiose), which were preliminarily derivatized in order to simplify the charge localization, and consequently the theoretical study.

When an anode electrode potential is larger than the lowest unoccupied molecular orbital (LUMO) of the electrolyte, Li-ions and electrolyte molecules will participate in reduction reactions on the anode surface and form a solid electrolyte interface (SEI) layer. Active Li-ion consumption in the formation reactions is the main source of capacity loss (>50) and ageing in Li-ion batteries (LIBs). Due to the fast-occurring and complex nature of the electrochemical processes, conventional experimental techniques are not a feasible approach for capturing and characterizing the SEI formation phenomenon. The lack of experimental data and consequently the absence of potential parameters for crystal structures in this layer makes molecular dynamics~(MD) simulations inapplicable to it. Also, due to the multi-component multi-layer structure of the SEI, the smallest system representing an SEI layer is too large for employing the principles of quantum mechanics~(QM), that traditionally work with much smaller system sizes. Addressing this, this thesis presents a novel computational framework for coupling QM and MD calculations to simulate a system with the size limits of MD simulations independent of the experimental data. The QM evaluates sub-atomic properties such as energy barriers against diffusion and employs seven new algorithms to estimate potential parameters as the input of the MD simulations. Then MD simulations forecast SEI's properties including density, Young's Modules, Poisson's Ratio, thermal conductivity, and diffusion coefficient mechanisms. The output of the QM and MD calculations are employed to develop two macro-scale mathematical models for predicting battery ageing and battery performance, incorporating the impact of the SEI layer in addition to the cathode, anode, and separator parts. Finally, the results obtained have been validated with respect to the experimental data in different operational conditions. / Thesis / Doctor of Philosophy (PhD) / The limited lifespan of expensive batteries is the main obstacle to electrification of the transport sector, despite its necessity for addressing the current environmental issues. Li+/electrolyte reduction on the electrode surface is responsible for more than 50% of capacity loss and the consequent ageing is a complex and fast-occurring phenomenon (few ns) that cannot be easily resolved using conventional experimental and computational techniques. This thesis presents the development of some computational frameworks and demonstrates their employment to investigate this phenomenon from a multi-scale perspective, i.e., from a few electrons to an entire battery length scale, with the operating cycles ranging from a few ps to several months, employing Quantum Mechanics, Molecular Dynamics, and Macro-Scale Modeling. The frameworks have been successfully validated with respect to experimental data from the literature and have been applied successfully to highlight the parameters that impact ageing in batteries. The findings presented in this thesis can be used as the base for further research on next-gen durable batteries with liquid and solid-state electrolytes.

<p>Bridging the accuracy of ab initio (AI) QM/MM with the efficiency of semi-empirical<br> (SE) QM/MM methods has long been a goal in computational chemistry. This dissertation<br> presents four ∆-Machine learning schemes aimed at achieving this objective. Firstly, the in-<br> corporation of negative force observations into the Gaussian process regression (GPR) model,<br> resulting in GPR with derivative observations, demonstrates the remarkable capability to<br> attain high-quality potential energy surfaces, accurate Cartesian force descriptions, and reli-<br> able free energy profiles using a training set of just 80 points. Secondly, the adaptation of the<br> sparse streaming GPR algorithm showcases the potential of memory retention from previous<br> phasespace, enabling energy-only models to converge using simple descriptors while faith-<br> fully reproducing high-quality potential energy surfaces and accurate free energy profiles.<br> Thirdly, the utilization of GPR with atomic environmental vectors as input features proves<br> effective in enhancing both potential energy surface and free energy description. Further-<br> more, incorporating derivative information on solute atoms further improves the accuracy<br> of force predictions on molecular mechanical (MM) atoms, addressing discrepancies arising<br> from QM/MM interaction energies between the target and base levels of theory. Finally, a<br> comprehensive comparison of three distinct GPR schemes, namely GAP, GPR with an aver-<br> age kernel, and GPR with a system-specific sum kernel, is conducted to evaluate the impact<br> of permutational invariance and atomistic learning on the model’s quality. Additionally, this<br> dissertation introduces the adaptation of the GAP method to be compatible with the sparse<br> variational Gaussian processes scheme and the streaming sparse GPR scheme, enhancing<br> their efficiency and applicability. Through these four ∆-Machine learning schemes, this dis-<br> sertation makes significant contributions to the field of computational chemistry, advancing<br> the quest for accurate potential energy surfaces, reliable force descriptions, and informative<br> free energy profiles in QM/MM simulations.<br> </p>

Fourier transform infrared difference spectroscopy (FTIR DS) is widely used to study the structural details of electron transfer cofactors in photosynthetic protein complexes. In photosynthetic proteins quinones play an important role, functioning as a cofactor in light-driven electron transfer. In photosystem I (PS I) phylloquinone (PhQ) functions as an intermediary in electron transfer. To investigate the properties of PhQ that occupies the, so called, A1 binding site in PS I, time-resolved step-scan FTIR DS, with 5µs time resolution at 77K has been used. By replacing PhQ in the A1 binding site with specifically isotope labeled version, information on the vibrational frequencies associated specifically with the quinone in the binding site were obtained, which could be compared to the vibrational properties of quinone in solution or quinones in other protein binding sites. To further aid in assessing the origin of bands in the spectra, quantum mechanics /molecular mechanics (QM/MM) ONIOM type calculations were undertaken. ONIOM is an acronym for Our own N-layered Integrated molecular Orbital and molecular Mechanics. We find that the phytyl tail of PhQ does not play an important role in the orientation of PhQ in the A1 binding site. We also find that PhQ, in both neutral and reduced states, is strongly hydrogen bonded. To test and verify the applicability of our QM/MM approach, ONIOM calculations were also undertaken for ubiquinone and a variety of other quinones incorporated into the, so called, QA binding site in purple bacteria photosynthetic reaction centers. The calculated and experimental spectra agree well, demonstrating the utility and applicability of our ONIOM approach. Hydrogen bonding to the carbonyl groups of quinones in the QA binding site was shown to be relatively weak, and it was found that hydrogen bonding to neutral ubiquinone in purple bacterial reaction centers can be considered in purely electrostatic terms, contrary to the widely held belief that the hydrogen bonding amino acids should be treated quantum mechanically.

No description available.

We have investigated the noncovalent interactions in carbohydrate-aromatic interactions which are pivotal to the recognition of carbohydrates in proteins. We have employed quantum mechanical methods to study carbohydrate-aromatic complexes. Due to the importance of dispersion contribution to the interaction energy, we mainly use density functional theory augmented with an empirical correction for the dispersion interactions (DFT-D). We have validated this method with a limited number of high level ab initio calculations. We have also analysed the vibrational and NMR chemical shift characteristics using the DFT-D method. We have mainly studied the complexes involving β-glucose with 3-methylindole and p-hydroxytoluene, which are analogues of tryptophan and tyrosine, respectively. We find that the contribution for interaction energy mainly comes from CH/π and OH/π interactions. We find that the interaction energy of complexes involving CH/π and OH/π interactions is reflected in the associated blue and red shifts of vibrational spectrum. We also find that the interactions involving 3-methylindole are somewhat greater than those for p-hydroxytoluene. The C-H blueshifts are also in parallel with the predicted NMR proton shift. We have also tested different density functionals including both standard density functionals and newly developed M0x functionals and MP2 method for studying carbohydrate-aromatic complexes. The DFT-D method and M06 functionals of the M0x family are found to perform better, while B3LYP and BLYP functionals perform poorly. We find that the inclusion of a dispersion term to BLYP is found to perform better. The dispersion energy dominates over the interaction energy of carbohydrate-aromatic complexes. From the DFT-D calculations, we found that the complexes would be unstable without the contribution from dispersive energy. We have also studied the importance of noncovalent interactions in functionalization of nanotubes by nucleic acid bases and aromatic amino acids by using semi-empirical methods with dispersion term such asPM3-D and PM3-D*. We find that the both semi-empirical schemes give reasonable interaction energies with respect to DFT-D interaction energies. We have also used PM3-D method to study the adsorption of organic pollutants on graphene sheet and on nanotubes. We found that the semi-empirical schemes, which are faster and cheaper, are suitable to study these larger molecules involving noncovalent interactions and can be used as an alternative to DFT-D method. We have also studied the importance of dispersion interaction and the effect of steric hindrance in aggregation of functionalized anthracenes and pentacenes. We have also employed molecular dynamics simulation methods to study the aggregation of anthracene molecules in toluene solution.

Noncovalent interactions are of great importance in studies on crystal design and drug discovery. One such noncovalent interaction, halogen bonding, is present between a covalently bound halogen atom and a Lewis base. A halogen bond is a directional interaction caused by the anisotropic distribution of charge on a halogen atom X covalently bound to A, which in turn forms a positive region called σ-hole on the A–X axis. Utilization of halogen bonds in lead optimization have been rarely considered in drug discovery until recently and yet more than 50% of the drug candidates are halogenated. To date, the halogen bond has not been subjected to practical molecular mechanical-molecular dynamics (MM-MD) study, where this noncovalent interaction cannot be described by conventional force fields because they do not account for the anisotropic distribution of the charge density on the halogen atoms. This problem was solved by the author and, for the first time, an extra-point of positive charge was used to represent the σ-hole on the halogen atom. This approach is called positive extra-point (PEP) approach. Interestingly, it was found that the performance of the PEP approach in describing halogen bond was better than the semiempirical methods including the recent halogen-bond corrected PM6 (PM6-DH2X) method. The PEP approach also gave promising results in describing other noncovalent halogen interactions, such as C–X···H and C–X···π-systems. The PEP resulted in an improvement in the accuracy of the electrostatic-potential derived charges of halogen-containing molecules, giving in turn better dipole moments and solvation free energies compared to high-level quantum mechanical and experimental data.With the aid of our PEP approach, the first MM-molecular dynamics (MM-MD) study of inhibitors that form a halogen bond with a receptor was performed for tetrahalobenzotriazole inhibitors complexed to cyclin-dependent protein kinase (CDK2). When the PEP approach was used, the calculated MM-generalized Born surface area (MM-GBSA)//MM-MD binding energies for halobenzimidazole and halobenzotriazole inhibitors complexed with protein kinase CK2 were found to correlate well with the corresponding experimental data, with correlation coefficients R2 of greater than 0.90. The nature and strength of halogen bonding in halo molecule···Lewis base complexes were studied in terms of molecular mechanics using our PEP approach. The contributions of the σ-hole (i.e., positively charged extra-point) and the halogen atom to the strength of this noncovalent interaction were clarified using the atomic parameter contribution to the molecular interaction approach. The molecular mechanical results revealed that the halogen bond is electrostatic and van der Waals in nature. The strength of the halogen bond increases with increasing the magnitude of the extra-point charge. The van der Waals interaction’s contribution to the halogen bond strength is most favorable in chloro complexes, whereas the electrostatic interaction is dominant in iodo complexes.The failure of the PM6 semiempirical method in describing noncovalent halogen interactions —not only halogen bonds, but also hydrogen bonds involving halogen atoms— was reported and corrected by the introduction of a second and third generation of noncovalent halogen interactions correction. The developed correction yielded promising results for the four examined noncovalent halogen interactions, namely: C–X···O, C–X···N, C–X···π-system, and C–X···H interactions.

Des études antérieures ont mis en évidence que le transfert de solutés neutres à travers des membranes est influencé par la présence d'ions en solution. Ainsi, la connaissance des interactions multiples à l'échelle nanométrique, entre le polymère, l'eau et les solutés (ions, espèces organiques) constituent un verrou pour l'amélioration des performances des procédés membranaires. Dans cette étude une approche multi-échelle fondamentale est proposée, combinant des outils théoriques et expérimentaux, afin d'obtenir les paramètres microscopiques et macroscopiques caractérisant les interactions étudiées pour différentes compositions ioniques. Plus précisément, il s'agit de comprendre comment les ions affectent le transfert d'un soluté organique. Dans un premier temps, certaines propriétés caractérisant l'hydratation des ions sont calculées et comparées aux flux de diffusions de sucres à travers des membranes de Nanofiltration et échangeuses d'ions obtenus pour différentes compositions ioniques. Dans un deuxième temps, des systèmes constitués d'une membrane échangeuse de cations (CMX) équilibrée avec différents cations ainsi que le glucose hydraté sont modélisés en utilisant une approche combinée Mécanique Quantique/ Mécanique Moléculaire. Cette approche a permis d'étudier la solubilité du sucre dans la matrice polymère ainsi que les interactions polymère-polymère comme l'énergie de cohésion. Enfin, l'influence des ions sur les caractéristiques physiques de la membrane CMX est étudiée en utilisant diverses méthodes expérimentales comme la détermination des angles de contacts et des spectres IR ou la mesure de la température de solidification par DSC. Les propriétés vibrationnelles sont également calculées dans le cadre de la théorie de la fonctionnelle de la densité (DFT). L'ensemble de ces données sont comparées avec les grandeurs de transport afin de valider les mécanismes moléculaires proposés. Ce travail montre que la nature des contre-ions de la membrane modifie l'énergie de cohésion entre les fragments de la membrane. Ainsi, l'énergie de cohésion influe sur la diffusion des composés organiques neutres à travers les membranes. / Previous works have shown that the transfer of neutral solutes through membranes is influenced by the presence of ions in solution. In the framework of process intensification, the knowledge of the molecular mechanisms involved is of fundamental importance to increase and predict the process performances. The aim of this Thesis is to use a combined quantum/molecular computational approach and experimental methodologies to better understand how ions can affect the solute flux. In the first part of the work, some properties of ions in solution are computed and compared with sugar fluxes through membranes for nanofiltration and electrodialysis. In the following, systems composed of Cation-exchange membrane equilibrated by different counter-ion and hydrated glucose are examined by Quantum Mechanics/Molecular Mechanics. This is done mainly to investigate the sugar solubility in the polymer matrix and diffusion related interactions like polymer chain-chain cohesion energy. In the last part, contact angle, differential scanning calorimetry and Infra-Red spectra are measured to characterize the physical properties of the membrane and possible influence of the counter-ion on cation exchange membrane. This work shows that the nature of the counter-ions modifies the cohesion energy between the membrane polymer fragments. In its turn, the cohesion energy affects the diffusion of neutral organic compounds through the membranes.

Glycosaminoglycans (GAGs) play an important role in many biological processes in the extracellular matrix. In a theoretical approach, structures of monosaccharide building blocks of natural GAGs and their sulfated derivatives were optimized by a B3LYP6311ppdd//B3LYP/ 6-31+G(d) method. The dependence of the observed conformational properties on the applied methodology is described. NMR chemical shifts and proton-proton spin-spin coupling constants were calculated using the GIAO approach and analyzed in terms of the method's accuracy and sensitivity towards the influence of sulfation, O1-methylation, conformations of sugar ring, and ω dihedral angle. The net sulfation of the monosaccharides was found to be correlated with the 1H chemical shifts in the methyl group of the N-acetylated saccharides both theoretically and experimentally. The ω dihedral angle conformation populations of free monosaccharides and monosaccharide blocks within polymeric GAG molecules were calculated by a molecular dynamics approach using the GLYCAM06 force field and compared with the available NMR and quantum mechanical data. Qualitative trends for the impact of sulfation and ring conformation on the chemical shifts and proton-proton spin-spin coupling constants were obtained and discussed in terms of the potential and limitations of the computational methodology used to be complementary to NMR experiments and to assist in experimental data assignment.