Ante la evidente necesidad de comprender los frecuentes cambios negativos que experimenta el planeta Tierra en los suelos, los ríos y la atmosfera y con la idea de intentar mitigar el cambio climático, aparece la química computacional como una rama de la química que nos permite estudiar y predecir diversas propiedadesfisicoquímicas de sistemas moleculares complejos de interés ambiental las cuales son, algunas veces, imposibles de estudiar de forma experimental. El presente artículo muestra un panorama molecular del comportamiento y características fisicoquímicas de algunos contaminantes atmosféricos estudiados a través de la química computacional / The urgent need to understand the frequent negative changes that the planet Earth suffer on its soil, rivers and atmosphere and the need to mitigate climate change makes computational chemistry stand out as a branch of chemistry that allows thestudy and prediction of different physicochemical properties of several complex molecular systems of environmental significance which are sometimes impossible to study in a laboratory. The present article shows a molecular overview of the behavior and physicochemical characteristics of some atmospheric pollutants studied through computational chemistry.

El Premio Nobel de Química de 2013 ha sido otorgado a A. Warshel, M. Levitt y M. Karplus debido al desarrollo de métodos híbridos de cálculo para química computacional. En este artículo se presentará una breve introducción del uso de los métodos de química computacional. Se describirá cómo se desarrollaron, y por qué, los métodos híbridos de cálculo, conocidos como QM/MM (Quantum Mechanics/Molecular Mechanics) para el estudio de sistemas macromoleculares, sobre todo para el caso de su aplicación en enzimas y bioquímica. Finalmente, se comentarán los alcances y expectativas futuras para estos métodos, desarrollados en los años 70. / The 2013 Nobel Prize in chemistry was awarded to A.Warshel, M.Levitt and M.Karplus for their contribution to the development of hybrid methods for computational chemistry. In this article a brief introduction about computational chemistry methods is presented. This paper will show the order in which the QM/MM (Quantum Mechanics/Molecular Mechanics) methods were developed for the study of macromolecular systems and specially their application in enzymes and biochemistry. Finally, the reach and future prospects of these methods originally developed by A. Warshel, M.Levitt and M. Karplus in the seventies will be discussed.

No description available.

Abstract The present thesis investigates the electronic structure of selected atoms and molecules in vapor phase. Electron spectroscopy is applied for studying the electronic transitions following excitation and ionization with electron and photon bombardment. The work focuses on the photoionization and Auger decay of selected noble gasses, and on the photoionization and Auger decay of core ionized or resonant excited alkali halide molecules. The experimental results are compared with theoretical predictions.

This thesis is concerned with the simulation of various electrochemical experiments and its application to electroanalysis. Chapter 1 outlines the fundamental principles of electrochemistry which are of importance for the reading of this thesis. Chapter 2 then outlines the methods used in the numerical simulation of electrochemical experiments. Chapters 3 and 4 are concerned with the electrochemistry of nanoparticles, and how this is affected by the presence of near wall hindered diffusion. In Chapter 3, a computational model to simulate anodic particle coulometry of nanoparticles in the presence of hindered diffusion is developed, and the effect of this hindered diffusion investigated. The model is then applied to simulate experimental data. Chapter 4 looks at the effect of hindered diffusion on the adsorption of nanoparticles on electrode surfaces, and investigates the effects of this adsorption on electrochemical experiments with nanoparticles generally. Chapters 5, 6 and 7 are concerned with band electrodes in isolation, in a pair and in an array respectively. In Chapter 5, a model to simulate double potential step chronoamperometry at an individual band electrode is developed, and used to successfully simulate experimental data. Chapter 6 looks at dual band electrodes used in generator-collector mode, and how this can be used to simultaneously measure the concentration of two species in solution. Chapter 7 looks at interdigitated arrays of band and ring electrodes in generator-collector mode, and develops a model to simulate cyclic voltammetry in both cases, as well as investigating under what conditions interdigitated ring arrays may be modelled as interdigitated band arrays. Chapter 8 develops a model to simulate chronoamperometry and cyclic voltammetry at porous electrodes, and investigates the consequences for electroanalysis of having a porous layer. Finally, Chapter 9 investigates the Marcus-Hush theory of electron transfer kinetics, and looks at the effect of the kinetically limited currents resulting from this theory to the equivalence relation between microdisc electrodes and sphere-on-a-surface electrodes.

Ce travail de thèse vise à étudier l'endommagement de l'ADN, de la formation de lésions à leur réparation par des méthodes de modélisation moléculaire. Plusieurs projets ont pris forme dans ce contexte, lesquels peuvent être classés en trois grandes catégories. D'un côté, nous nous sommes intéressés la formation de lésions induites par des agents mutagènes. Nous avons étudié les mécanismes de formation de la 8-oxo-7,8-dihydroguanine (8oxoG), mais aussi le caractère de photosensibilisateur endogène de la pyrimidine 6-4 pyrimidone (6-4PP), et la photosensibilisation de l'ADN par deux anti-inflammatoires : le kétoprofène et l'ibuprofène. D'un autre côté, les propriétés mécaniques de l'ADN endommagé ont été simulées. La structure de lésions complexes est d'une importance capitale pour comprendre la manière dont elles sont réparées. Malheureusement, seulement peu de structures RMN et cristallographiques sont disponibles à ce jour. Pour pallier à ce manque et obtenir des informations sur leur dynamique, nous avons étudié un panel de lésions complexes : les clusters de sites abasiques, les pontages inter-brins, et la photolésion 6-4PP. De même, nous nous sommes penchés sur les modes d'interaction de certaines polyamines avec l'ADN, ces molécules étant connues pour interagir avec la double hélice. Enfin, latroisième partie de cette thèse concerne les interactions ADN-enzyme de réparation. En perspective avec l'étude de clusters d'abasiques, nous avons étudié le comportement dynamique du même système, cette fois-ci en interaction avec l'endonucléase APE1. Nous nous sommes également penchés sur les interactions entre la glycosylase Fpg avec un oligonucléotide contenant un tandem de lésions 8-oxoG d'un côté, etun cluster de lésions 8-oxoG - site abasique de l'autre. Ces multiples projets ont permis l'accumulation de nouvelles connaissances à propos des lésions complexes de l'ADN, et ont également apporté un appuicomputationnel aux expérimentations, qui peuvent se révéler très délicates dans ce domaine. Nos résultats ouvrent de larges perspectives dans le domaine de la pharmacologie, la cosmétique et plus généralementla compréhension du vivant / This thesis work is focused on the theoretical modelling of DNA damages, from formation to repair. Several projects have been led in this framework, which can be sorted into three different parts. One on hand, we studied complex DNA reactivity. It included a study about 8-oxo-7,8-dihydroguanine (8oxoG) mechanisms of formation, a project concerning the UV-induced pyrimidine 6-4 pyrimidone (6-4PP) endogenous photosensitizer features, and an other one about DNA photosensitizationby nonsteroidal anti-inflammatory drugs (ie ketoprofen and ibuprofen). On the other hand, we investigated mechanical properties of damaged DNA. The structural signature of a DNA lesion is of major importance for their repair, unfortunately only few NMR and X-ray structures of such systems are available. In order to gain insights into their dynamical structure, we investigated a series of complex damages : clustered abasic sites, interstrand cross-links, and the 6-4PP photolesion. Likewise, we studied the interaction modes DNA with several polyamines, which are well known to interact with the double helix, but also with the perspective to model DNA-protein cross-linking. The third part concerned the study of DNA interactions with repair enzymes. In line with the structural study about clustered abasic sites, we investigated the dynamics of the same system, but this time interacting with the APE1 endonuclease. We also studied interactions between the Fpg glycosylase with an oligonucleotides containing tandem 8-oxoG on one hand and 8-oxoG - abasic site as multiply damaged sites. Thus, we shed new lights on damaged DNA reactivity, structure and repair, which provides perspectives for biomedicine and life's mechanisms understanding as we begin to describe nucleosomal DNA

This work explores the use of computational methods to aid in the design of Metal Organic Frameworks (MOFs) for use as CO2 scrubbers in carbon capture and storage applications. One of the main challenges in this field is in identifying important MOF design characteristics which optimize the complex interactions governing surface adsorption. We approach this in a high-throughput manner, determining properties important to CO2 adsorption from generating and sampling a large materials search space. The utilization of MOFs as potential carbon scrubbing agents is a recent phenomenon, as such, many of the computational tools necessary to perform high-throughput screening of MOFs and subsequent analysis are either underdeveloped or non-existent. A large portion of this work therefore involved the development of novel tools designed specifically for this task. The chapters in this work are contiguous with the goal of designing MOFs for CO¬2 capture, and somewhat chronological in order and complexity, meaning as time and expertise progressed, more advanced tools were developed and utilized for the purposes of computational MOF discovery. Initial work towards MOF design involved the detailed analysis of two experimental structures; CALF-15 and CALF-16 using classical molecular dynamics, grand canonical Monte Carlo simulations, and DFT to determine the structural features which promote CO2 adsorption. An unprecedented level of agreement was found between theory and experiment, as we are able to capture, with simulation, the X-ray resolved binding sites of CO2 in the confined pores of CALF-15. Molecular simulation was then used to provide a detailed breakdown of the energy contributions from nearby functional groups in both CALF-15 and CALF-16. A large database of hypothetical MOF structures is constructed for the purposes of screening for CO2 adsorption. The database contains 1.3 million hypothetical structures, generated with an algorithm which snaps together rigid molecular building blocks extracted from existing MOF crystal structures. The algorithm for constructing the hypothetical MOFs and the building blocks themselves were all developed in-house to form the resulting database. The topological, chemical, and physical features of these MOFs are compared to recently developed materials databases to demonstrate the larger structural and chemical space sampled by our database. In order to rapidly and accurately describe the electrostatic interactions of CO2 in the hypothetical database of MOFs, parameters were developed for use with the charge equilibration method. This method assigns partial charges on the framework atoms based on a set of parameters assigned to each atom type. An evolutionary algorithm was used to optimize the charge equilibration parameters on a set of 543 hypothetical MOFs such that the partial charges generated would reproduce each MOFs DFT-derived electrostatic potential. Validation of these parameters was performed by comparing the CO2 adsorption from the charge equilibration method vs DFT-derived charges on a separate set of 693 MOFs. Our parameter set were found to reproduce DFT-derived CO2 adsorption extremely well using only a fraction of the time, making this method ideal for rapid and accurate high-throughput MOF screening. A database of 325,000 MOFs was then screened for CO2 capture and storage applications. From this study we identify important binding pockets for CO2 in MOFs using a binding site analysis tool. This tool uses a pattern recognition method to compare the 3-D configurations of thousands of pore structures surrounding strong CO2 adsorption sites, and present common features found amongst them. For the purposes of developing larger databases which sample a more diverse materials space, a novel MOF construction tool is devloped which builds MOFs based on abstract graphs. The graph theoretical foundations of this method are discussed and several examples of MOF construction are presented to demonstrate its use. Notably, not only can it build existing MOFs with complicated geometries, but it can sample a wide range of unique structures not yet discovered by experimental means.

A series of para-substituted phenylethylamine analogues bound to the enzyme aromatic amine dehydrogenase have been simulated using quantum mechanical electronic structure calculations and molecular mechanical molecular dynamics simulations. Trends have been verified connecting bond dissociation energy (and thus driving force) to observed rate constants and activation enthalpy. Trends have been identified in connecting statistics drawn from molecular dynamics simulations and the temperature dependence of the kinetic isotope effect, notably that as the temperature dependence of the kinetic isotope effect increases the flexibility of the promoting vibration decreases. This is explained as being more effected by thermal energy put into the system, and therefore more affected by temperature.

Accumulation and aggregation of the 42-residue amyloid-[beta] (A[beta]) protein fragment, which originates from the cleavage of amyloid precursor protein by beta and gamma secretase, correlates with the pathology of Alzheimer's disease (AD). Possible therapies for AD include peptides based on the A[beta] sequence, and recently identified small molecular weight compounds designed to mimic these, that interfere with the aggregation of A[beta] and prevent its toxic effects on neuronal cells in culture. Here, we use molecular dynamics simulations to compare the mode of interaction of an active (LPFFD) and inactive (LHFFD) [beta]-sheet breaker peptide with an A[beta] fibril structure from solid state NMR studies. We found that LHFFD had a weaker interaction with the fibril than the active peptide, LPFFD, from geometric and energetic considerations, as estimated by the MM/PBSA approach. Cluster analysis and computational alanine scanning identified important ligand-fibril contacts, including a possible difference in the effect of histidine on ligand-fibril [pi]-stacking interactions, and the role of the proline residue establishing contacts that compete with those essential for maintenance of the inter-monomer [beta]-sheet structure of the fibril. Our results show that molecular dynamics simulations can be a useful way to classify the stability of docking sites. These mechanistic insights into the ability of LPFFD to reverse aggregation of toxic A[beta] will guide the redesign of lead compounds, and aid in developing realistic therapies for AD and other diseases of protein aggregation. We have also performed long explicit solvent MD simulations of unliganded amyloid fibril in three putative protonation states, in order to better understand the energetic and mechanical features of the fibril receptor. Over 100 ns MD simulations, the trajectories where fibril has Glu11 and Glu22 side-chains protonated exhibit the least deviation from the initial solid state NMR structures. Free energy calculations on these rajectories suggest that the weakest fibril interface lies in the lateral rather than transverse direction and that there is little dependence on whether the lateral interface is situated at the edge or middle of the fibril. This agrees with recent reported steered molecular dynamics calculations. Secondly, in an effort to improve the ability of atomistic simulation techniques to directly resolve protein tertiary structure from primary amino acid sequence, we explore the use of a molecular dynamics technique based on swarm intelligence, called SWARM-MD, to identify the native states of two peptides, polyalanine and AEK17, as well as Trp-cage miniprotein. We find that the presence of cooperative swarm interactions significantly enhanced the efficiency of molecular dynamics simulations in predicting native conformation. However, it also is evident that the presence of outlying simulation replicas can adversely impact correctly folded replica structures. By slowly removing the swarm potential after folding simulations, the negative effect of the swarm potential can be alleviated and better agreement with experiment obtained.

Enzymes direct chemical reactions with precision and speed, making life as we know it possible. How they do this is still not completely understood, but the relatively recent discovery of subpicosecond protein motion coupled to the reaction coordinate has provided a crucial piece of the puzzle. This type of motion is called a rate-promoting vibration (RPV) and has been seen in a number of different enzymatic systems. It typically involves a compression of the active site of the enzyme which lowers the barrier for the reaction to occur. In this work we present a number of studies that probe these motions in two dehydrogenase enzymes, yeast alcohol dehydrogenase (YADH) and homologs of lactate dehydrogenase (LDH). The goal of the study on the reaction of YADH was to probe the role of the protein in proton tunneling in the enzyme, which was suggested to occur from experimental kinetic isotope effect studies. We did this using transition path sampling (TPS), which perturbatively generates ensembles of reactive trajectories to observe transitions between stable states, such as chemical reactions. By applying a quantum method that can account for proton tunneling, centroid molecular dynamics, and generating reactive trajectory ensembles with and without the method, we were able to observe the change in barrier to proton transfer upon application of the tunneling method. We found that there was little change in the barrier, showing that classical over-the-barrier transfer is dominant over tunneling in the proton transfer in YADH. We also applied the knowledge of RPVs to identify a new class of allosteric molecules, which modulate enzymatic reaction not by changing a binding affinity, but by disrupting the reactive motion of enzymes. We showed, through design of a novel allosteric effector for human heart LDH, applying TPS to a system with and without the small molecule bound, and analysis of the reaction coordinate of the reactive trajectory ensemble, that the molecule was able to disrupt the motion of the protein such that it was no longer coupled to the reaction. We also examined the subpicosecond motions of two other LDHs, from Plasmodium falciparum and Cryptosporidium parvum, which evolved separately from previously studied LDHs. We found, using TPS and reaction coordinate identification, that while the LDH from C. parvum had similar dynamics to the earlier LDHs, the LDH from P. falciparum had a earlier transition-state associated with proton transfer, not hydride transfer. This is likely due to this LDH having a larger active site pocket, increasing the amount of motion necessary for proton transfer, and, thus, the barrier to proton transfer. More work is necessary in this system to determine whether the protein is coupled with the search for the reactive conformation for proton transfer. Protein motion coupled to the particle transfer in dehydrogenases plays an important role in their reactions and there is still much work to be done to understand the extent and role of RPVs.