The Application of Computational Chemistry and Chemometrics to Developing a Method for Online Monitoring of Polymer Degradation in the Manufacture of Bioresorbable Medical Implants

Bioresorbable polymers such as PLA have an important role to play in the development of temporary implantable medical devices with significant benefits over traditional therapies. However, development of new devices is hindered by high manufacturing costs associated with difficulties in processing the material. A major problem is the lack of insight on material degradation during processing. In this work, a method of quantifying degradation of PLA using IR spectroscopy coupled with computational chemistry and chemometric modeling is examined. It is shown that the method can predict the quantity of degradation products in solid-state samples with reasonably good accuracy, indicating the potential to adapt the method to developing an on-line sensor for monitoring PLA degradation in real-time during processing.

Single and multiple addition to C60. A computational chemistry study

Des del seu descobriment, a la molècula C60 se li coneixen una varietat de derivats segons el tipus de funcionalització amb propietats fisicoquímiques específiques de gran interès científic. Una sel·lecció de derivats corresponents a addicions simple o múltiple al C60 s'ha considerat en aquest treball d'investigació. L'estudi a nivell de química computacional de diversos tipus d'addició al C60 s'han portat a terme per tal de poder donar resposta a aspectes que experimentalment no s'entenen o són poc clars. Els sistemes estudiats en referència a l'addició simple al C60 han estat en primer lloc els monoiminoful·lerens, C60NR, (de les dues vies proposades per la seva síntesi, anàlisis cinètic i termodinàmic han ajudat a explicar els mecanismes de formació i justificar l'addició a enllaços tipus [5,6]), i en segon lloc els metanoful·lerens i els hidroful·lerens substituits, C60CHR i C60HR, (raons geomètriques, electròniques, energètiques i magnètiques justifiquen el diferent caràcter àcid ente ambdós derivats tenint en compte una sèrie de substituents R amb diferent caràcter electrònic donor/acceptor). Els fluoroful·lerens, C60Fn, i els epoxid ful·lerens, C60On, (anàlisi sistemàtic dels seus patrons d'addició en base a poder justificar la força que els governa han aportat dades complementàries a les poques que existeixen experimentalment al respecte).

Computational chemistry for elucidating protein function: energetics and dynamics of myoglobin-ligand systems

State-of-the-art computational methodologies are used to investigate the energetics and dynamics of photodissociated CO and NO in myoglobin (Mb···CO and Mb···NO). This includes the combination of molecular dynamics, ab initio MD, free energy sampling, and effective dynamics methods to compare the results with studies using X-ray crystallography and ultrafast spectroscopy metho ds. It is shown that modern simulation techniques along with careful description of the intermolecular interactions can give quantitative agreement with experiments on complex molecular systems. Based on this agreement predictions for as yet uncharacterized species can be made.

Cyclic Voltammetry and Computational Chemistry Studies on the Evaluation of the Redox Behavior of Parabens and other Analogues

Parabens are antimicrobial preservatives widely used in pharmaceutical, cosmetic and food industries. The alkyl chain connected to the ester group defines some important physicochemical characteristics of these compounds, including the partition coefficient and redox properties. The voltammetric and computational analyses were carried out in order to evaluate the redox behavior of these compounds and other phenolic analogues. A strong correlation between chemical substituents inductive effects of parabens with redox potentials was observed. Using cyclic voltammetry and glassy carbon working electrode, only one irreversible anodic peak was observed around 0.8 V for methylparaben (MP), ethylparaben (EP), propylparaben (PP), butylparaben (BP), benzylparaben (BzP) and p-substituted phenolic analogues. The electrodonating inductive effect of alkyl groups was demonstrated by the anodic oxidation potential shift to lower values as the carbon number increases and, therefore the parabens (and other phenolic analogues) oxidation processes to the quinonoidic forms showed great dependence on the substituent pattern.

Cyclic voltammetry and computational chemistry studies on the evaluation of the redox behavior of parabens and other analogues

Parabens are antimicrobial preservatives widely used in pharmaceutical, cosmetic and food industries. The alkyl chain connected to the ester group defines some important physicochemical characteristics of these compounds, including the partition coefficient and redox properties. The voltammetric and computational analyses were carried out in order to evaluate the redox behavior of these compounds and other phenolic analogues. A strong correlation between chemical substituents inductive effects of parabens with redox potentials was observed. Using cyclic voltammetry and glassy carbon working electrode, only one irreversible anodic peak was observed around 0.8 V for methylparaben (MP), ethylparaben (EP), propylparaben (PP), butylparaben (BP), benzylparaben (BzP) and p-substituted phenolic analogues. The electrodonating inductive effect of alkyl groups was demonstrated by the anodic oxidation potential shift to lower values as the carbon number increases and, therefore the parabens (and other phenolic analogues) oxidation processes to the quinonoidic forms showed great dependence on the substituent pattern.

Computational chemistry studies of UV induced processes in human skin

This thesis presents and uses the techniques of computational chemistry to explore two different processes induced in human skin by ultraviolet light. The first is the transformation of urocanic acid into a immunosuppressing agent, and the other is the enzymatic action of the 8-oxoguanine glycosylase enzyme. The photochemistry of urocanic acid is investigated by time-dependent density functional theory. Vertical absorption spectra of the molecule in different forms and environments is assigned and candidate states for the photochemistry at different wavelengths are identified. Molecular dynamics simulations of urocanic acid in gas phase and aqueous solution reveals considerable flexibility under experimental conditions, particularly for for the cis isomer where competition between intra- and inter-molecular interactions increases flexibility. A model to explain the observed gas phase photochemistry of urocanic acid is developed and it is shown that a reinterpretation in terms of a mixture between isomers significantly enhances the agreement between theory and experiment , and resolves several peculiarities in the spectrum. A model for the photochemistry in the aqueous phase of urocanic acid is then developed, in which two excited states governs the efficiency of photoisomerization. The point of entrance into a conical intersection seam is shown to explain the wavelength dependence of photoisomerization quantum yield. Finally some mechanistic aspects of the DNA repair enzyme 8-oxoguanine glycosylase is investigated with density functional theory. It is found that the critical amino acid of the active site can provide catalytic power in several different manners, and that a recent proposal involving a SN1 type of mechanism seems the most efficient one.

Investigation of a binding model for anti-cancer cyclohexadienones by computational chemistry and crystallography

DUE TO COPYRIGHT RESTRICTIONS ONLY AVAILABLE FOR CONSULTATION AT ASTON UNIVERSITY LIBRARY AND INFORMATION SERVICES WITH PRIOR ARRANGEMENT

I. Quantal effects in biochemical cooperativity and a proposed mechanism for the differentiation of calcium signaling in synaptic plasticity. II. Evolutionary algorithms for the optimization of methods in computational chemistry

In Part 1 of this thesis, we propose that biochemical cooperativity is a fundamentally non-ideal process. We show quantal effects underlying biochemical cooperativity and highlight apparent ergodic breaking at small volumes. The apparent ergodic breaking manifests itself in a divergence of deterministic and stochastic models. We further predict that this divergence of deterministic and stochastic results is a failure of the deterministic methods rather than an issue of stochastic simulations.

Ergodic breaking at small volumes may allow these molecular complexes to function as switches to a greater degree than has previously been shown. We propose that this ergodic breaking is a phenomenon that the synapse might exploit to differentiate Ca$^{2+}$ signaling that would lead to either the strengthening or weakening of a synapse. Techniques such as lattice-based statistics and rule-based modeling are tools that allow us to directly confront this non-ideality. A natural next step to understanding the chemical physics that underlies these processes is to consider \textit{in silico} specifically atomistic simulation methods that might augment our modeling efforts.

In the second part of this thesis, we use evolutionary algorithms to optimize \textit{in silico} methods that might be used to describe biochemical processes at the subcellular and molecular levels. While we have applied evolutionary algorithms to several methods, this thesis will focus on the optimization of charge equilibration methods. Accurate charges are essential to understanding the electrostatic interactions that are involved in ligand binding, as frequently discussed in the first part of this thesis.

Modern Computational Chemistry Methods for Prediction of Ground- and Excited-State Properties in Open-Shell Systems

The main goal of the research presented in this work is to provide some important insights about computational modeling of open-shell species. Such projects are: the investigation of the size-extensivity error in Equation-of-Motion Coupled Cluster methods, the analysis of the Long-Range corrected scheme in predicting UV-Vis spectra of Cu(II) complexes with the 4-imidazole acetate and its ethylated derivative, and the exploration of the importance of choosing a proper basis set for the description of systems such as the lithium monoxide anion. The most significant findings of this research are: (i) The contribution of the left operator to the size-extensivity error of the CR-EOMCC(2,3) approach, (ii) The cause of d-d shifts when varying the range-separation parameter and the amount of the exact exchange arising from the imbalanced treatment of localized vs. delocalized orbitals via the "tuned" CAM-B3LYP* functional, (iii) The proper acidity trend of the first-row hydrides and their lithiated analogs that may be reversed if the basis sets are not correctly selected.

Studies of new inorganic species using relativistic quantum chemistry

In the present work the methods of relativistic quantum chemistry have been applied to a number of small systems containing heavy elements, for which relativistic effects are important. First, a thorough introduction of the methods used is presented. This includes some of the general methods of computational chemistry and a special section dealing with how to include the effects of relativity in quantum chemical calculations. Second, after this introduction the results obtained are presented. Investigations on high-valent mercury compounds are presented and new ways to synthesise such compounds are proposed. The methods described were applied to certain systems containing short Pt-Tl contacts. It was possible to explain the interesting bonding situation in these compounds. One of the most common actinide compounds, uranium hexafluoride was investigated and a new picture of the bonding was presented. Furthermore the rareness of uranium-cyanide compounds was discussed. In a foray into the chemistry of gold, well known for its strong relativistic effects, investigations on different gold systems were performed. Analogies between Au$^+$ and platinum on one hand and oxygen on the other were found. New systems with multiple bonds to gold were proposed to experimentalists. One of the proposed systems was spectroscopically observed shortly afterwards. A very interesting molecule, which was theoretically predicted a few years ago is WAu$_{12}$. Some of its properties were calculated and the bonding situation was discussed. In a further study on gold compounds it was possible to explain the substitution pattern in bis[phosphane-gold(I)] thiocyanate complexes. This is of some help to experimentalists as the systems could not be crystallised and the structure was therefore unknown. Finally, computations on one of the heaviest elements in the periodic table were performed. Calculation on compounds containing element 110, darmstadtium, showed that it behaves similarly as its lighter homologue platinum. The extreme importance of relativistic effects for these systems was also shown.

Computational investigations on covalent dimerization/oligomerization of polyacenes: Is it relevant to soot formation?

We have postulated a novel pathway that could assist in the nucleation of soot particles through covalent dimerization and oligomerizations of a variety of PAHs. DFT calculations were performed with the objective of obtaining the relative thermal stabilities and formation probabilities of oligomeric species that exploit the facile dimerization that is known to occur in linear oligoacenes. We propose that the presence of small stretches of linear oligoacence (tetracene or longer) in extended PAH, either embedded or tethered, would be adequate for enabling the formation of such dimeric and oligomeric adducts; these could then serve as nuclei for the growth of soot particles. Our studies also reveal the importance of p-stacking interactions between extended aromatic frameworks in governing the relative stabilities of the oligomeric species that are formed. (c) 2012 Wiley Periodicals, Inc.

Computational analysis of quinoxalines N,N-Dioxide and ITS Derivates

A quinoxalina e seus derivativos são uma importante classe de compostos heterocíclicos, onde os elementos N, S e O substituem átomos de carbono no anel. A fórmula molecular da quinoxalina é C8H6N2, formada por dois anéis aromáticos, benzeno e pirazina. É rara em estado natural, mas a sua síntese é de fácil execução. Modificações na estrutura da quinoxalina proporcionam uma grande variedade de compostos e actividades, tais como actividades antimicrobiana, antiparasitária, antidiabética, antiproliferativa, anti-inflamatória, anticancerígena, antiglaucoma, antidepressiva apresentando antagonismo do receptor AMPA. Estes compostos também são importantes no campo industrial devido, por exemplo, ao seu poder na inibição da corrosão do metal. A química computacional, ramo natural da química teórica é um método bem desenvolvido, utilizado para representar estruturas moleculares, simulando o seu comportamento com as equações da física quântica e clássica. Existe no mercado uma grande variedade de ferramentas informaticas utilizadas na química computacional, que permitem o cálculo de energias, geometrias, frequências vibracionais, estados de transição, vias de reação, estados excitados e uma variedade de propriedades baseadas em várias funções de onda não correlacionadas e correlacionadas. Nesta medida, a sua aplicação ao estudo das quinoxalinas é importante para a determinação das suas características químicas, permitindo uma análise mais completa, em menos tempo, e com menos custos.

A joint computational and experimental study of a novel dioxomolybdenum(VI) complex bearing chiral N,N-dimethyllactamide ligand

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Computational Study of the Hydration of Sulfuric Acid Dimers: Implications for Acid Dissociation and Aerosol Formation

We have investigated the thermodynamics of sulfuric acid dimer hydration using ab initio quantum mechanical methods. For (H2SO4)2(H2O)n where n = 0−6, we employed high-level ab initio calculations to locate the most stable minima for each cluster size. The results presented herein yield a detailed understanding of the first deprotonation of sulfuric acid as a function of temperature for a system consisting of two sulfuric acid molecules and up to six waters. At 0 K, a cluster of two sulfuric acid molecules and one water remains undissociated. Addition of a second water begins the deprotonation of the first sulfuric acid leading to the di-ionic species (the bisulfate anion HSO4−, the hydronium cation H3O+, an undissociated sulfuric acid molecule, and a water). Upon the addition of a third water molecule, the second sulfuric acid molecule begins to dissociate. For the (H2SO4)2(H2O)3 cluster, the di-ionic cluster is a few kcal mol−1 more stable than the neutral cluster, which is just slightly more stable than the tetra-ionic cluster (two bisulfate anions, two hydronium cations, and one water). With four water molecules, the tetra-ionic cluster, (HSO4−)2(H3O+)2(H2O)2, becomes as favorable as the di-ionic cluster H2SO4(HSO4−)(H3O+)(H2O)3 at 0 K. Increasing the temperature favors the undissociated clusters, and at room temperature we predict that the di-ionic species is slightly more favorable than the neutral cluster once three waters have been added to the cluster. The tetra-ionic species competes with the di-ionic species once five waters have been added to the cluster. The thermodynamics of stepwise hydration of sulfuric acid dimer is similar to that of the monomer; it is favorable up to n = 4−5 at 298 K. A much more thermodynamically favorable pathway forming sulfuric acid dimer hydrates is through the combination of sulfuric acid monomer hydrates, but the low concentration of sulfuric acid relative to water vapor at ambient conditions limits that process.

Computational study of ligand binding in lipid transfer proteins: Structures, interfaces, and free energies of protein-lipid complexes

Plant nonspecific lipid transfer proteins (nsLTPs) bind a wide variety of lipids, which allows them to perform disparate functions. Recent reports on their multifunctionality in plant growth processes have posed new questions on the versatile binding abilities of these proteins. The lack of binding specificity has been customarily explained in qualitative terms on the basis of a supposed structural flexibility and nonspecificity of hydrophobic protein-ligand interactions. We present here a computational study of protein-ligand complexes formed between five nsLTPs and seven lipids bound in two different ways in every receptor protein. After optimizing geometries inmolecular dynamics calculations, we computed Poisson- Boltzmann electrostatic potentials, solvation energies, properties of the protein-ligand interfaces, and estimates of binding free energies of the resulting complexes. Our results provide the first quantitative information on the ligand abilities of nsLTPs, shed new light into protein-lipid interactions, and reveal new features which supplement commonly held assumptions on their lack of binding specificity.