CO2 capture in amine solutions: modelling and simulations with non-empirical methods

Wanda Andreoni, Fabio Pietrucci
Iop Publishing Ltd, 2016
Article

Résumé

Absorption in aqueous amine solutions is the most advanced technology for the capture of CO2, although suffering from drawbacks that do not allow exploitation on large scale. The search for optimum solvents has been pursued with empirical methods and has also motivated a number of computational approaches over the last decade. However, a deeper level of understanding of the relevant chemical reactions in solution is required so as to contribute to this effort. We present here a brief critical overview of the most recent applications of computer simulations using ab initio methods. Comparison of their outcome shows a strong dependence on the structural models employed to represent the molecular systems in solution and on the strategy used to simulate the reactions. In particular, the results of very recent ab initio molecular dynamics augmented with metadynamics are summarized, showing the crucial role of water, which has been so far strongly underestimated both in the calculations and in the interpretation of experimental data. Indications are given for advances in computational approaches that are necessary if meant to contribute to the rational design of new solvents.

Official source

https://infoscience.epfl.ch/record/224912?ln=fr

À propos de ce résultat

Cette page est générée automatiquement et peut contenir des informations qui ne sont pas correctes, complètes, à jour ou pertinentes par rapport à votre recherche. Il en va de même pour toutes les autres pages de ce site. Veillez à vérifier les informations auprès des sources officielles de l'EPFL.

Concepts associés

Chargement

Publications associées

Chargement

Bridging Static and Dynamical Descriptions of Chemical Reactions: An ab Initio Study of CO2 Interacting with Water Molecules

Wanda Andreoni, Fabio Pietrucci

Extracting reliable thermochemical parameters from molecular dynamics simulations of chemical reactions, although based on ab initio methods, is generally hampered by difficulties in reproducing the results and controlling the statistical errors. This is a serious drawback with respect to the quantum chemical description based on potential energy surfaces. This work is an attempt to fill this gap. We apply molecular dynamics, based on density functional theory (DFT) and empowered by path metadynamics (MTD), to simulate the reaction of CO2 with (one, two, and three) water molecules in the gas phase. This study relies on a strategy that ensures a precise control of the accuracy of the reaction coordinates and of the reconstructed free-energy surface on this space, namely, on (i) fully reversible MTD simulations, (ii) a committor probability analysis for the diagnosis of the collective variables, and (iii) a cluster analysis for the characterization of the reconstructed free energy surfaces. This robust procedure permits a meaningful comparison with more traditional calculations of the potential energy surfaces that we also perform within the same DFT computational scheme. This comparison shows in particular that the reactants and products of systems with only three water molecules can no longer be understood in terms of one structure but must be described as statistical configuration ensembles. Calculations carried out with different prescriptions for the exchange-correlation functionals also allow us to establish their quantitative effect on the activation barriers for the formation and the dissociation of carbonic acid. Their decrease induced by the addition of one water molecule (catalytic effect) is found to be largely independent of the specific functional.

Amer Chemical Soc2012

Chargement

Over the past years, molecular modeling and simulation techniques have had a major impact on experimental life sciences. They are capable of providing accurate insight into microscopic mechanisms, which are usually difficult to investigate experimentally. Moreover, the integration of experimental data with molecular modeling appears to be a promising strategy to better understand complex biological processes. In this thesis, molecular dynamics simulation has been used in combination with experimental data to investigate two transmembrane proteins: (i) the bacterial chemoreceptor PhoQ and (ii) the Amyloid Precursor Protein (APP). (i) Bacterial two-component system PhoQ and bacterial membranes. Two-component systems (TCSs) are signaling complexes essential for bacterial survival and virulence. PhoQ is the histidine kinase chemoreceptor of the PhoQ-PhoP tandem machine that detects the concentration of cationic species at the inner membrane of Gram-negative bacteria. A full understanding of the PhoQ signal transduction mechanism is currently hindered by the lack of a complete atomistic structure. In this thesis project, the first structural model of the transmembrane (TM) portion of PhoQ from Escherichia coli was assembled, by using molecular simulations integrated with cross-linking disulfide scanning data. Its structural and dynamic features induce a concerted displacement of the TM helices at the periplasmic side, which modulates a rotation at the cytoplasmic end. This supports the idea that signal transduction is promoted through a combination of scissoring and rotational movements of the TM helices. Knowledge of this complex mechanism is essential in order to understand how the chemical stimuli sensed by the periplasmic sensor domain trigger, via the relay of the HAMP domain, the histidine auto-phosphorylation and kinase/phosphatase activity at the cytoplasmic end. The PhoQ sensor domain lies in close proximity to the membrane. Interaction with anionic lipids, such as phosphophatidylglycerols (PG) and cardiolipins (CL), are thought to play a key role in PhoQ activity. Present in bacterial and mitochondrial membranes, cardiolipins have a unique dimeric structure, which carries up to two charges, i.e. one per phosphate group, and under physiological conditions, can be unprotonated or singly protonated. Exhaustive models and characterization of cardiolipins are to-date scarce; therefore an ab initio parameterization of cardiolipin species for molecular simulation consistent with commonly used force fields is proposed here. Molecular dynamics (MD) simulations based on these models indicate a protonation-dependent lipid packing. A noteworthy interaction with solvating mono- and divalent cations is also observed. The proposed models will contribute to the biophysical and biochemical characterizations of bacterial and mitochondrial membranes and membrane-embedded proteins. (ii) Structural and dynamic properties of the Amyloid Precursor Protein. The Amyloid Precursor Protein (APP) is a type I membrane glycoprotein present at the neuronal synapsis. The proteolytic cleavage of its C-terminal segment produces amyloid-β (Aβ) peptides of different lengths, the deposition of which is an early indicator of Alzheimer"s disease (AD). Recently, the backbone structure of the APP transmembrane (TM) domain in detergent micelles was determined by nuclear magnetic resonance (NMR, independently by two different experimental groups). The TM conformations of these two structures are however markedly different. One is characterized by a highly kinked α-helix, whereas the other is mainly straight. Molecular dynamics simulations showed that the APP TM region is highly flexible and its secondary structure is influenced by the surrounding lipid environment. The size of the embedding detergent micelles strongly affects the conformation of the APP α-helix, with solvation being the main driving force for the development of a helical curvature. Once embedded in a membrane bilayer, APP systematically prefers a straight helical conformation. This is also confirmed when analyzing in silico the atomistic APP population observed in double electron-electron resonance (DEER) spectroscopy. In summary, the APP transmembrane domain is highly flexible due to the presence of glycine residues and can readily respond to the lipid environment, a property that might be critical for proteolytic processing by γ-secretase enzymes. The presented thesis work clearly shows how molecular simulations and their interplay with available experimental input can help advance the understanding of the mechanism of complex biological systems and processes on a molecular scale. These results, in particular, go well beyond the current understanding of the functioning of two transmembrane proteins relevant for human health. Furthermore, the computational approaches and procedures developed in these projects will hopefully promote novel integrated strategies for investigating biological systems.

EPFL2013

Chargement

Modeling non-covalent interactions in condensed phase

Veronika Juraskova

The modeling of non-covalent interactions, solvation effects, and chemical reactions in complex molecular environment is a challenging task. Current state-of-the-art approaches often rely on static computations using implicit solvent models and harmonic approximations. Owing to the rich interaction network characteristic of molecular species in the condensed phase, the neglect of the dynamic fluctuations and of the mutual competition with the solvent molecules often leads to an incomplete picture.While ab initio molecular dynamics can, in principle, capture all these effects, the associated computational cost requires compromising the accuracy of the electronic structure method and the statistical convergence of the simulation.The main objective of this thesis is to develop a robust simulation workflow based on neural network potentials and enhanced sampling techniques to describe molecules and reactions in complex environment. The validity and usefulness of the proposed workflow is illustrated on three applications: (a) the solvation environment and protonation state of a mixture of phenol, hydrogen peroxide and methanesulfonic acid relevant to catalyzed industrial phenol hydroxylation, (b) the extension of the above to model the hydroxylation reaction, and (c) the description of a chalcogen-bonded complex in tetrahydrofuran featuring multiple and heavy chemical elements as a potential route for anion transport and catalysis. The developed workflow is illustrated in detail through the first example, for which the solvation environment and protonation states are dominated by strongly competing hydrogen bonds. We demonstrate the relevance of relying upon surrogate neural network potentials, enhanced sampling techniques, and nuclear quantum effects to achieve an accurate description of the complex mixture. The second example extends the simulation pipeline, combining the neural network potentials with metadynamics to enable the modeling of the catalyzed phenol hydroxylation reaction with a special emphasis on its regioselectivity. Finally, the third example broadens the application of neural network potentials to the description of molecular systems containing up to seven chemical elements, which participate in very subtle and specific types of non-covalent interactions.Overall, this work highlights the potential of using neural network-based potentials to model complex reactive systems in condensed phase and to accurately describe competing non-covalent interactions. The developed workflow extends the range of applicability of ab initio molecular dynamics toolbox to real-life chemical systems.

EPFL2022

Chargement

Modeling non-covalent interactions in condensed phase

Veronika Juraskova

EPFL2022

EPFL2013