Pushing the Frontiers of First-Principles Based Computer Simulations of Chemical and Biological Systems

Negar Ashari Astani, Elizabeth Claire Brunk, Pablo Campomanes Ramos, Basile Curchod, Polydefkis Diamantis, Andrey Laktionov, Marilisa Neri, Giulia Palermo, Thomas James Penfold, Ursula Röthlisberger, Ivano Tavernelli, Stefano Vanni
2011
Journal paper

Abstract

The Laboratory of Computational Chemistry and Biochemistry is active in the development and application of first-principles based simulations of complex chemical and biochemical phenomena. Here, we review some of our recent efforts in extending these methods to larger systems, longer time scales and increased accuracies. Their versatility is illustrated with a diverse range of applications, ranging from the determination of the gas phase structure of the cyclic decapeptide gramicidin S, to the study of G protein coupled receptors, the interaction of transition metal based anti-cancer agents with protein targets, the mechanism of action of DNA repair enzymes, the role of metal ions in neurodegenerative diseases and the computational design of dye-sensitized solar cells. Many of these projects are done in collaboration with experimental groups from the Institute of Chemical Sciences and Engineering (ISIC) at the EPFL.

Official source

https://infoscience.epfl.ch/record/170025?ln=en

About this result

This page is automatically generated and may contain information that is not correct, complete, up-to-date, or relevant to your search query. The same applies to every other page on this website. Please make sure to verify the information with EPFL's official sources.

Related concepts

Related publications

Radicals play an important role in many areas of chemistry, such as atmospheric, aquatic, polymer, and biological, to name a few. Radicals are often highly reactive species which are short-lived and therefore harder to study by experimental techniques. In principle, computational quantum chemistry enables the study of arbitrary chemical species, including elusive radicals. However, the approximate quantum chemical methods which are used for production simulations need to be validated for the type of problem in question. First, it is necessary to validate that highly accurate reference methods commonly applied to closed-shell systems are equally applicable to radicals. I show that the “gold standard” for closed-shell systems (coupled cluster with singles, doubles, and perturbative triples excitations (CCSD(T))) can produce molecular geometries and vibrational frequencies suitable for (sub-)kJ/mol thermochemistry for radical species as well. However, additional diagnosis is necessary, especially for vibrational frequency calculations. Coupled cluster methods are then used to compute benchmark binding energies between small radicals and water or hydrogen fluoride as a proxy to the aqueous condensed phase. I find that many common density functional theory (DFT) methods cannot reproduce these binding energies to within chemical accuracy (1 kcal/mol). However, some functionals, for example MPW1K and BHandHLYP, perform better than others. If more accurate ab initio methods are computationally too demanding for the system of interest, these DFT methods should be used for (micro-)solvated radicals. Complexes of radicals with polar solvent are found to be more strongly bound than expected. This can not be explained by dispersive and electrostatic interactions alone. Using natural bond orbitals, I show that additional donor-acceptor interactions may be responsible for the strong intermolecular binding observed. One of the most fundamental processes involving aqueous organic radicals is the one-electron oxidation of closed-shell organics. I study the vertical ionization of neutral organics. That is, the ionization is faster than the reorganization of the solvent (and solute), and the conformation is assumed to be unchanged during the ionization process. I find that the difference between the gas-phase and aqueous-phase vertical ionization energy, ∆VIE = VIEaq − VIEgas, ranges from ≈ 0 to -1.0 eV for small aromatic compounds. The differences in ∆VIE between different compounds seem to be caused by the solvent structures of the ≈ 70 water molecules closest to the solute. In natural surface waters, oxidized organics arise through the reaction of photochemically created radicals with organic molecules already present in the water column. One example for such a reaction is the oxidation of sulfonamide antibiotics. I show that upon oxidation of sulfadiazine and related compounds, an aniline radical cation is formed. The reactivity of this oxidized species is strikingly different from its reduced counterpart: the barrier for a nucleophilic attack on the aniline ring is drastically lowered, and a fast rearrangement reaction can take place. I demonstrate that current quantum chemical protocols are suitable for mechanistic studies of solvated radicals, and I successfully apply them to such problems. To achieve chemical accuracy for quantum computational protocols, further developments are necessary. To validate novel approaches, the benchmark data presented in this thesis may be used.

EPFL2013

Multi-objective optimisation (MOO) of super-structured process designs are expensive in CPU- time because of the high number of potential configurations and operation conditions to be calculated. Thus single process units are generally represented by simple models like equilibrium based (chemical or phase equilibrium) or specific short cut models. In the development of new processes, kinetic effects or mass transport limitations in certain process units may play an important role, especially in multiphase chemical reactors. Therefore, it is desirable to represent such process unit by experimentally derived RBMs (rate based models) in the process flowsheet simulators used for the extensive MOO. This increases the trust engineers have in the results and allows enriching the process simulations with newest experimental findings. As most RBMs are iteratively solved, a direct incorporation would cause higher CPU-time that penalises the use of MOO. A global surrogate model (SUMO) of a RBM was successfully generated to allow its incorporation into a process design & optimisation (PDO) tool which makes use of an evolutionary MOO. The methodology was applied to a fluidised bed methanation reactor in the process chain from wood to Synthetic Natural Gas (SNG). Two types of surrogate model, an ordinary Kriging interpolation and an artificial neural network, were generated and compared to its underlying rate based model and the chemical equilibrium model. The analysis showed that kinetic limitations have significant influence on the result already for standard bulk gas chemical components. From experimental experience it can be stated that a significant amount of chemical compounds, which are near to complete conversion according the thermodynamic equilibrium, are measured in the gas phase after the reactor. Thus including RBMs in PDO improves the quality of results. This approach allows a significant improvement of information exchange between process design & optimisation workflow and experimental development of PUs in early stages of process development.

2013

Modeling, predicting and mining metabolism at atom-level resolution

Jasmin Maria Hafner

Living organisms can catalyze many thousands of biochemical reactions that they use to convert energy and matter, which provides them with the essentials for life. The sum of these chemical reactions happening in an organism is called metabo-lism. Understanding metabolism is crucial to elucidate the fundamental principles of biology, and further to enable us to redesign it for the sustainable, biosynthetic pro-duction of bio-based fuels, commodity chemicals and medicines. Mathematical mod-els are essential to organize and understand the complexity of metabolism. They usu-ally represent metabolism as a network of reactions, but they tend to neglect the exact molecular structure of the metabolites. To redesign metabolic reactions, how-ever, a mechanistic understanding of metabolic reactions and their catalysts, proteins called enzymes, is essential. In this work, a mathematical description of enzymatic reaction mechanism, called generalized reaction rules, is applied to computationally simulate and predict meta-bolic processes at the level of atoms. Each reaction rule describes the catalytic activi-ty of an enzyme, or a group of enzymes, at the mechanistic level by encoding the re-arrangement of atoms in the reaction. The reaction rules are called âgeneralizedâ, because they mimic the ability of a single enzyme to catalyze multiple reactions by acting on a range of substrates. Using these reaction rules, we first developed a computational representation of me-tabolism that allowed tracking single atoms throughout complex metabolic reaction networks. The principle of atom-tracking was then used to develop a graph-theory based method to represent and analyze metabolic networks, and to reliably identify of metabolic pathways for the biosynthesis of chemicals. Next, we applied the gener-alized reaction rules to predict all possible novel, hypothetical reactions from known biological compounds, and we stored the five million generated novel reactions them in a database called ATLAS. Finally, the developed tools and resources were applied to specific engineering and research problems, such as the biosynthetic pathway design for the biofuel bisabolene and the plastic precursor 1,4-butanediol. We further pre-dicted a biosynthesis route for the pharmaceutical tetrahydropalmatine and engi-neered a yeast strain to produce it. Finally, we show that our tools can be used to mine available genome sequences to find organisms that can degrade xenobiotics. Our findings suggest that the atom-level representation of metabolism can greatly contribute to its understanding, exploration and prediction. Given the complexity of atom-level modeling of metabolic processes, we propose metrics that can approxi-mate the atom-level information to conserve the information at the level of big, hy-pothetical metabolic networks like ATLAS. This database plus the developed pathway search techniques form a valuable resource for scientist to help characterizing un-known biosynthesis pathways towards secondary metabolites, and for metabolic engi-neers for designing novel bioproduction pathways for chemicals. Hopefully, these considerations will contribute to a better understanding of metabolism, advance the exploration of the bioproduction of drugs and other valuable molecules, and acceler-ate metabolic engineering efforts to realize the switch from a petroleum-based chem-ical industry towards a more sustainable, bio-based production of societyâs chemical needs.

EPFL2020

EPFL2013

2013

Modeling, predicting and mining metabolism at atom-level resolution

Jasmin Maria Hafner

EPFL2020