skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


  • NSF Public Access
  • Search Results
  • Computational Investigation of the Thermochemistry of the CO2 Capture Reaction by Ethylamine, Propylamine, and Butylamine in Aqueous Solution Considering the Full Conformational Space via Boltzmann Statistics
Title: Computational Investigation of the Thermochemistry of the CO2 Capture Reaction by Ethylamine, Propylamine, and Butylamine in Aqueous Solution Considering the Full Conformational Space via Boltzmann Statistics
Rubisco is the enzyme responsible for CO2 fixation in nature, and it is activated by CO2 addition to the amine group of its lysine 201 side chain. We are designing rubisco-based biomimetic systems for reversible CO2 capture from ambient air. The oligopeptide biomimetic capture systems are employed in aqueous solution. To provide a solid foundation for the experimental solution-phase studies of the CO2 capture reaction, we report here the results of computational studies of the thermodynamics of CO2 capture by small alkylamines in aqueous solution. We studied CO2 addition to methyl-, ethyl-, propyl-, and butylamine with the consideration of the full conformational space for the amine and the corresponding carbamic acids and with the application of an accurate solvation model for the potential energy surface analyses. The reaction energies of the carbamylation reactions were determined based on just the most stable structures (MSS) and based on the ensemble energies computed with the Boltzmann distribution (BD), and it is found that ΔGBD ≈ ΔGMSS. The effect of the proper accounting for the molecular translational entropies in solution with the Wertz approach are much more significant, and the free energy of the capture reactions ΔWABD is more negative by 2.9 kcal/mol. Further accounting for volume effects in solution results in our best estimates for the reaction energies of the carbamylation reactions of ΔWABD = −5.4 kcal/mol. The overall difference is ΔGBD – ΔWABD = 2.4 kcal/mol for butylamine carbamylation. The full conformational space analyses inform about the conformational isomerizations of carbamic acids, and we determined the relevant rotational profiles and their transition-state structures. Our detailed studies emphasize that, more generally, solution-phase reaction energies should be evaluated with the Helmholtz free energy and can be affected substantially by solution effects on translational entropies.  more » « less
Award ID(s):
2153206
PAR ID:
10511047
Author(s) / Creator(s):
; ;
Publisher / Repository:
ACS Publications
Date Published:
Journal Name:
The Journal of Physical Chemistry A
Volume:
125
Issue:
44
ISSN:
1089-5639
Page Range / eLocation ID:
9578 to 9593
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; (, ChemPhysChem)
    Abstract We have been interested in the development of rubisco‐based biomimetic systems for reversible CO2capture from air. Our design of the chemical CO2capture and release (CCR) system is informed by the understanding of the binding of the activator CO2(ACO2) in rubisco (ribulose‐1,5‐bisphosphate carboxylase/oxygenase). The active site consists of the tetrapeptide sequence Lys‐Asp‐Asp‐Glu (or KDDE) and the Lys sidechain amine is responsible for the CO2capture reaction. We are studying the structural chemistry and the thermodynamics of CO2capture based on the tetrapeptide CH3CO−KDDE−NH2(“KDDE”) in aqueous solution to develop rubisco mimetic CCR systems. Here, we report the results of1H NMR and13C NMR analyses of CO2capture by butylamine and by KDDE. The carbamylation of butylamine was studied to develop the NMR method and with the protocol established, we were able to quantify the oligopeptide carbamylation at much lower concentration. We performed a pH profile in the multi equilibrium system and measured amine species and carbamic acid/carbamate species by the integration of1H NMR signals as a function of pH in the range 8≤pH≤11. The determination of ΔG1(R) for the reaction R−NH2+CO2R−NH−COOH requires the solution of a multi‐equilibrium equation system, which accounts for the dissociation constantsK2andK3controlling carbonate and bicarbonate concentrations, the acid dissociation constantK4of the conjugated acid of the amine, and the acid dissociation constantK5of the alkylcarbamic acid. We show how the multi‐equilibrium equation system can be solved with the measurements of the daughter/parent ratioX, the knowledge of the pH values, and the initial concentrations [HCO3−]0and [R‐NH2]0. For the reaction energies of the carbamylations of butylamine and KDDE, our best values are ΔG1(Bu)=−1.57 kcal/mol and ΔG1(KDDE)=−1.17 kcal/mol. Both CO2capture reactions are modestly exergonic and thereby ensure reversibility in an energy‐efficient manner. These results validate the hypothesis that KDDE‐type oligopeptides may serve as reversible CCR systems in aqueous solution and guide designs for their improvement. 
    more » « less
  2. ; ; (, The Journal of Organic Chemistry)
    We are developing energy-efficient and reversible carbon capture and release (CCR) systems that mimic the Lys201 carbamylation reaction in the active site of ribulose-1,5-bisphosphate carboxylase-oxygenase (RuBisCO). The multiequilibria scenario ammonium ion Xa ⇌ amine Xb ⇌ carbamic acid Xc ⇌ carbamate Xd requires the presence of both free amine and CO2 for carbamylation and is affected by the pKa(Xa). Two fluorination strategies aimed at ammonium ion pKa depression and low pH carbamylation were analyzed with (2,2,2-trifluoroethyl)butylamine 2b and 2,2-difluoropropylamine 3b and compared to butylamine 1b. The determination of K1 and ΔG1 of the carbamylation reactions requires the solution of multiequilibria systems of equations based on initial conditions, 1H NMR measurements of carbamylation yields over a wide pH range, and knowledge of K2– K5 values. K2 and K3 describe carbonic acid acidity, and ammonium ion acidities K4 were measured experimentally. We calibrated carbamic acid acidities K5 based on the measured value K6 of aminocarbamic acid using isodesmic reactions. The proton exchange reactions were evaluated with ab initio computations at the APFD/6-311+G* level in combination with continuum solvation models and explicit solvation. The utilities of 1–3 will be discussed as they pertain to the development of fluorine-modified RuBisCO-mimetic reversible CCR systems. 
    more » « less
  3. ; (, The journal of physical chemistry A)
    Atmospheric CO2 concentrations have been growing steadily with no sign of moderation. Rubisco is the ubiquitous plant enzyme responsible for carbon fixation in biomass. Expanding luscious forests on Earth would increase CO2 absorption, but extreme drought, deforestation, desertification, and urbanization encroaching on agricultural land demonstrate the opposite trend. Direct air capture (DAC) from ambient air is the only meaningful way to reduce atmospheric CO2. We are designing Rubisco-inspired reversible CO2 DAC systems based on the lysine carbamylation associated with Rubisco activation. Here, we report the results of computational studies of the thermodynamics of CO2 capture by small alkylamines in aqueous solution to learn about the carbamylation of the amino group in lysine's sidechain. We determined the reaction energies of the carbamylation reactions based on the traditional approach of considering just the most stable structures and based on the ensemble energies computed with the Boltzmann distribution. View the article. 
    more » « less
  4. ; ; ; ; ; ; ; (, Angewandte Chemie International Edition)
    Abstract Reaction of {LiC6H2−2,4,6‐Cyp3⋅Et2O}2(Cyp=cyclopentyl) (1) of the new dispersion energy donor (DED) ligand, 2,4,6‐triscyclopentylphenyl with SnCl2afforded a mixture of the distannene {Sn(C6H2−2,4,6‐Cyp3)2}2(2), and the cyclotristannane {Sn(C6H2−2,4,6‐Cyp3)2}3(3).2is favored in solution at higher temperature (345 K or above) whereas3is preferred near 298 K. Van't Hoff analysis revealed the3to2conversion has a ΔH=33.36 kcal mol−1and ΔS=0.102 kcal mol−1 K−1, which gives a ΔG300 K=+2.86 kcal mol−1, showing that the conversion of3to2is an endergonic process. Computational studies show that DED stabilization in3is −28.5 kcal mol−1per {Sn(C6H2−2,4,6‐Cyp3)2unit, which exceeds the DED energy in2of −16.3 kcal mol−1per unit. The data clearly show that dispersion interactions are the main arbiter of the3to2equilibrium. Both2and3possess large dispersion stabilization energies which suppress monomer dissociation (supported by EDA results). 
    more » « less
  5. ; ; ; ; (, Physical Chemistry Chemical Physics)
    Experimentally conducted reactions between CO 2 and various substrates ( i.e. , ethylenediamine (EDA), ethanolamine (ETA), ethylene glycol (EG), mercaptoethanol (ME), and ethylene dithiol (EDT)) are considered in a computational study. The reactions were previously conducted under harsh conditions utilizing toxic metal catalysts. We computationally utilize Brønsted acidic ionic liquid (IL) [Et 2 NH 2 ]HSO 4 as a catalyst aiming to investigate and propose ‘greener’ pathways for future experimental studies. Computations show that EDA is the best to fixate CO 2 among the tested substrates: the nucleophilic EDA attack on CO 2 is calculated to have a very small energy barrier to overcome (TS1EDA, Δ G ‡ = 1.4 kcal mol −1 ) and form I1EDA (carbamic acid adduct). The formed intermediate is converted to cyclic urea (PEDA, imidazolidin-2-one) via ring closure and dehydration of the concerted transition state (TS2EDA, Δ G ‡ = 32.8 kcal mol −1 ). Solvation model analysis demonstrates that nonpolar solvents (hexane, THF) are better for fixing CO 2 with EDA. Attaching electron-donating and -withdrawing groups to EDA does not reduce the energy barriers. Modifying the IL via changing the anion part (HSO 4 − ) central S atom with 6 A and 5 A group elements (Se, P, and As) shows that a Se-based IL can be utilized for the same purpose. Molecular dynamics (MD) simulations reveal that the IL ion pairs can hold substrates and CO 2 molecules via noncovalent interactions to ease nucleophilic attack on CO 2 . 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email