1
|
Růžička K, Štejfa V, Červinka C, Fulem M, Šturala J. Thermodynamic Study of N-Methylformamide and N, N-Dimethyl-Formamide. Molecules 2024; 29:1110. [PMID: 38474622 DOI: 10.3390/molecules29051110] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/10/2024] [Revised: 02/27/2024] [Accepted: 02/27/2024] [Indexed: 03/14/2024] Open
Abstract
An extensive thermodynamic study of N-methylformamide (CAS RN: 123-39-7) and N,N-dimethylformamide (CAS RN: 68-12-2), is presented in this work. The liquid heat capacities of N-methylformamide were measured by Tian-Calvet calorimetry in the temperature interval (250-300) K. The vapor pressures for N-methylformamide and N,N-dimethylformamide were measured using static method in the temperature range 238 K to 308 K. The ideal-gas thermodynamic properties were calculated using a combination of the density functional theory (DFT) and statistical thermodynamics. A consistent thermodynamic description was developed using the method of simultaneous correlation, where the experimental and selected literature data for vapor pressures, vaporization enthalpies, and liquid phase heat capacities and the calculated ideal-gas heat capacities were treated together to ensure overall thermodynamic consistency of the results. The resulting vapor pressure equation is valid from the triple point to the normal boiling point temperature.
Collapse
Affiliation(s)
- Květoslav Růžička
- Department of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, CZ-166 28 Prague, Czech Republic
| | - Vojtěch Štejfa
- Department of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, CZ-166 28 Prague, Czech Republic
| | - Ctirad Červinka
- Department of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, CZ-166 28 Prague, Czech Republic
| | - Michal Fulem
- Department of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, CZ-166 28 Prague, Czech Republic
| | - Jiří Šturala
- Department of Inorganic Chemistry, University of Chemistry and Technology, Prague, Technická 5, CZ-166 28 Prague, Czech Republic
| |
Collapse
|
2
|
Maste S, Sharma B, Pongratz T, Grabe B, Hiller W, Erlach MB, Kremer W, Kalbitzer HR, Marx D, Kast SM. The accuracy limit of chemical shift predictions for species in aqueous solution. Phys Chem Chem Phys 2024; 26:6386-6395. [PMID: 38315169 DOI: 10.1039/d3cp05471c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2024]
Abstract
Interpreting NMR experiments benefits from first-principles predictions of chemical shifts. Reaching the accuracy limit of theory is relevant for unambiguous structural analysis and dissecting theoretical approximations. Since accurate chemical shift measurements are based on using internal reference compounds such as trimethylsilylpropanesulfonate (DSS), a detailed comparison of experimental with theoretical data requires simultaneous consideration of both target and reference species ensembles in the same solvent environment. Here we show that ab initio molecular dynamics simulations to generate liquid-state ensembles of target and reference compounds, including explicitly their short-range solvation environments and combined with quantum-mechanical solvation models, allows for predicting highly accurate 1H (∼0.1-0.5 ppm) and aliphatic 13C (∼1.5 ppm) chemical shifts for aqueous solutions of the model compounds trimethylamine N-oxide (TMAO) and N-methylacetamide (NMA), referenced to DSS without any system-specific adjustments. This encompasses the two peptide bond conformations of NMA identified by NMR. The results are used to derive a general-purpose guideline set for predictive NMR chemical shift calculations of NMA in the liquid state and to identify artifacts of force field models. Accurate predictions are only obtained if a sufficient number of explicit water molecules is included in the quantum-mechanical calculations, disproving a purely electrostatic model of the solvent effect on chemical shifts.
Collapse
Affiliation(s)
- Stefan Maste
- Fakultät für Chemie und Chemische Biologie, Technische Universität Dortmund, Otto-Hahn-Straße 4a, 44227 Dortmund, Germany.
| | - Bikramjit Sharma
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany.
| | - Tim Pongratz
- Fakultät für Chemie und Chemische Biologie, Technische Universität Dortmund, Otto-Hahn-Straße 4a, 44227 Dortmund, Germany.
| | - Bastian Grabe
- Fakultät für Chemie und Chemische Biologie, Technische Universität Dortmund, Otto-Hahn-Straße 4a, 44227 Dortmund, Germany.
| | - Wolf Hiller
- Fakultät für Chemie und Chemische Biologie, Technische Universität Dortmund, Otto-Hahn-Straße 4a, 44227 Dortmund, Germany.
| | - Markus Beck Erlach
- Fakultät für Biologie und Vorklinische Medizin, Universität Regensburg, 93040 Regensburg, Germany
| | - Werner Kremer
- Fakultät für Biologie und Vorklinische Medizin, Universität Regensburg, 93040 Regensburg, Germany
| | - Hans Robert Kalbitzer
- Fakultät für Biologie und Vorklinische Medizin, Universität Regensburg, 93040 Regensburg, Germany
| | - Dominik Marx
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany.
| | - Stefan M Kast
- Fakultät für Chemie und Chemische Biologie, Technische Universität Dortmund, Otto-Hahn-Straße 4a, 44227 Dortmund, Germany.
| |
Collapse
|
3
|
Imamura K, Yokogawa D, Sato H. Recent developments and applications of reference interaction site model self-consistent field with constrained spatial electron density (RISM-SCF-cSED): A hybrid model of quantum chemistry and integral equation theory of molecular liquids. J Chem Phys 2024; 160:050901. [PMID: 38341702 DOI: 10.1063/5.0190116] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2023] [Accepted: 01/04/2024] [Indexed: 02/13/2024] Open
Abstract
The significance of solvent effects in electronic structure calculations has long been noted, and various methods have been developed to consider this effect. The reference interaction site model self-consistent field with constrained spatial electron density (RISM-SCF-cSED) is a hybrid model that combines the integral equation theory of molecular liquids with quantum chemistry. This method can consider the statistically convergent solvent distribution at a significantly lower cost than molecular dynamics simulations. Because the RISM theory explicitly considers the solvent structure, it performs well for systems where hydrogen bonds are formed between the solute and solvent molecules, which is a challenge for continuum solvent models. Taking advantage of being founded on the variational principle, theoretical developments have been made in calculating various properties and incorporating electron correlation effects. In this review, we organize the theoretical aspects of RISM-SCF-cSED and its distinctions from other hybrid methods involving integral equation theories. Furthermore, we carefully present its progress in terms of theoretical developments and recent applications.
Collapse
Affiliation(s)
- Kosuke Imamura
- Department of Molecular Engineering, Graduate School of Engineering, Kyoto University, Kyoto 615-8510, Japan
| | - Daisuke Yokogawa
- Graduate School of Arts and Science, The University of Tokyo, Komaba, Meguro-ku, Tokyo 153-8902, Japan
| | - Hirofumi Sato
- Department of Molecular Engineering, Graduate School of Engineering, Kyoto University, Kyoto 615-8510, Japan
- Fukui Institute for Fundamental Chemistry, Kyoto University, Kyoto 606-8103, Japan
| |
Collapse
|
4
|
Tran VA, Teucher M, Galazzo L, Sharma B, Pongratz T, Kast SM, Marx D, Bordignon E, Schnegg A, Neese F. Dissecting the Molecular Origin of g-Tensor Heterogeneity and Strain in Nitroxide Radicals in Water: Electron Paramagnetic Resonance Experiment versus Theory. J Phys Chem A 2023; 127:6447-6466. [PMID: 37524058 PMCID: PMC10424240 DOI: 10.1021/acs.jpca.3c02879] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2023] [Revised: 07/01/2023] [Indexed: 08/02/2023]
Abstract
Nitroxides are common EPR sensors of microenvironmental properties such as polarity, numbers of H-bonds, pH, and so forth. Their solvation in an aqueous environment is facilitated by their high propensity to form H-bonds with the surrounding water molecules. Their g- and A-tensor elements are key parameters to extracting the properties of their microenvironment. In particular, the gxx value of nitroxides is rich in information. It is known to be characterized by discrete values representing nitroxide populations previously assigned to have different H-bonds with the surrounding waters. Additionally, there is a large g-strain, that is, a broadening of g-values associated with it, which is generally correlated with environmental and structural micro-heterogeneities. The g-strain is responsible for the frequency dependence of the apparent line width of the EPR spectra, which becomes evident at high field/frequency. Here, we address the molecular origin of the gxx heterogeneity and of the g-strain of a nitroxide moiety (HMI: 2,2,3,4,5,5-hexamethylimidazolidin-1-oxyl, C9H19N2O) in water. To treat the solvation effect on the g-strain, we combined a multi-frequency experimental approach with ab initio molecular dynamics simulations for structural sampling and quantum chemical EPR property calculations at the highest realistically affordable level, including an explicitly micro-solvated HMI ensemble and the embedded cluster reference interaction site model. We could clearly identify the distinct populations of the H-bonded nitroxides responsible for the gxx heterogeneity experimentally observed, and we dissected the role of the solvation shell, H-bond formation, and structural deformation of the nitroxide in the creation of the g-strain associated with each nitroxide subensemble. Two contributions to the g-strain were identified in this study. The first contribution depends on the number of hydrogen bonds formed between the nitroxide and the solvent because this has a large and well-understood effect on the gxx-shift. This contribution can only be resolved at high resonance frequencies, where it leads to distinct peaks in the gxx region. The second contribution arises from configurational fluctuations of the nitroxide that necessarily lead to g-shift heterogeneity. These contributions cannot be resolved experimentally as distinct resonances but add to the line broadening. They can be quantitatively analyzed by studying the apparent line width as a function of microwave frequency. Interestingly, both theory and experiment confirm that this contribution is independent of the number of H-bonds. Perhaps even more surprisingly, the theoretical analysis suggests that the configurational fluctuation broadening is not induced by the solvent but is inherently present even in the gas phase. Moreover, the calculations predict that this broadening decreases upon solvation of the nitroxide.
Collapse
Affiliation(s)
- Van Anh Tran
- Max-Planck-Institut
für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Markus Teucher
- Max-Planck-Institut
für Chemische Energiekonversion, Stiftstraße 34-36, 45470 Mülheim an der Ruhr, Germany
| | - Laura Galazzo
- Department
of Physical Chemistry, University of Geneva, Quai Ernest Ansermet 30, 1211 Geneva, Switzerland
- Faculty
of Chemistry and Biochemistry, Ruhr-Universität
Bochum, 44780 Bochum, Germany
| | - Bikramjit Sharma
- Lehrstuhl
für Theoretische Chemie, Ruhr-Universität
Bochum, 44780 Bochum, Germany
| | - Tim Pongratz
- Fakultät
für Chemie und Chemische Biologie, Technische Universität Dortmund, Otto-Hahn-Str. 4a, 44227 Dortmund, Germany
| | - Stefan M. Kast
- Fakultät
für Chemie und Chemische Biologie, Technische Universität Dortmund, Otto-Hahn-Str. 4a, 44227 Dortmund, Germany
| | - Dominik Marx
- Lehrstuhl
für Theoretische Chemie, Ruhr-Universität
Bochum, 44780 Bochum, Germany
| | - Enrica Bordignon
- Department
of Physical Chemistry, University of Geneva, Quai Ernest Ansermet 30, 1211 Geneva, Switzerland
- Faculty
of Chemistry and Biochemistry, Ruhr-Universität
Bochum, 44780 Bochum, Germany
| | - Alexander Schnegg
- Max-Planck-Institut
für Chemische Energiekonversion, Stiftstraße 34-36, 45470 Mülheim an der Ruhr, Germany
| | - Frank Neese
- Max-Planck-Institut
für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
5
|
Reimann M, Kaupp M. Spin-State Splittings in 3d Transition-Metal Complexes Revisited: Toward a Reliable Theory Benchmark. J Chem Theory Comput 2023; 19:97-108. [PMID: 36576816 DOI: 10.1021/acs.jctc.2c00925] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
A new composite method for the calculation of spin-crossover energies in 3d transition-metal complexes based on multireference methods is presented. The method reduces to MRCISD+Q at the complete-basis-set (CBS) level for atomic ions, for which it gives excitation energies with a mean absolute error of only ca. 0.01 eV. For molecular complexes, the CASPT2+δMRCI composite approach corresponds to a CASPT2/CBS calculation augmented by a high-level MRCISD+Q-CASPT2 correction with a smaller ligand basis set. For a set of [Fe(He)6]n+ test complexes, the approach reproduces full MRCISD+Q/CBS results to within better than 0.04 eV, without depending on any arbitrary IPEA shifts. The high-quality CASPT2+δMRCI method has then been applied to a series of 3d transition-metal hexaqua complexes in aqueous solution, augmented by an elaborate 3D-RISM-SCF solvent treatment of the underlying structures. It provides unprecedented agreement with experiment for the lowest-lying vertical spin-flip excitation energies, except for the Fe3+ system. Closer examination of the latter case provides strong evidence that the observed lowest-energy excitation at 1.56 eV, which has been used frequently for evaluating quantum-chemical methods, does not arise from the iron(III) hexaqua complex in solution, but from its singly deprotonated counterpart, [Fe(H2O)5OH]2+.
Collapse
Affiliation(s)
- Marc Reimann
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Technische Universität Berlin, Sekr. C7, Straße des 17. Juni 135, D-10623 Berlin, Germany
| | - Martin Kaupp
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Technische Universität Berlin, Sekr. C7, Straße des 17. Juni 135, D-10623 Berlin, Germany
| |
Collapse
|
6
|
Imamura K, Yokogawa D, Higashi M, Sato H. Reference interaction site model self-consistent field with constrained spatial electron density approach for nuclear magnetic shielding in solution. J Chem Phys 2022; 157:204105. [DOI: 10.1063/5.0122326] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022] Open
Abstract
We propose a new hybrid approach combining quantum chemistry and statistical mechanics of liquids for calculating the nuclear magnetic resonance (NMR) chemical shifts of solvated molecules. Based on the reference interaction site model self-consistent field with constrained spatial electron density distribution (RISM–SCF–cSED) method, the electronic structure of molecules in solution is obtained, and the expression for the nuclear magnetic shielding tensor is derived as the second-order derivative of the Helmholtz energy of the solution system. We implemented a method for calculating chemical shifts and applied it to an adenine molecule in water, where hydrogen bonding plays a crucial role in electronic and solvation structures. We also performed the calculations of 17O chemical shifts, which showed remarkable solvent dependence. While converged results could not be sometimes obtained using the conventional method, in the present framework with RISM–SCF–cSED, an adequate representation of electron density is guaranteed, making it possible to obtain an NMR shielding constant stably. This introduction of cSED is key to extending the method’s applicability to obtain the chemical shift of various chemical species. The present demonstration illustrates our approach’s superiority in terms of numerical robustness and accuracy.
Collapse
Affiliation(s)
- Kosuke Imamura
- Department of Molecular Engineering, Graduate School of Engineering, Kyoto University, Kyoto 615-8510, Japan
| | - Daisuke Yokogawa
- Graduate School of Arts and Science, The University of Tokyo, Komaba, Meguro-ku, Tokyo 153-8902, Japan
| | - Masahiro Higashi
- Department of Molecular Engineering, Graduate School of Engineering, Kyoto University, Kyoto 615-8510, Japan
| | - Hirofumi Sato
- Department of Molecular Engineering, Graduate School of Engineering, Kyoto University, Kyoto 615-8510, Japan
- Fukui Institute for Fundamental Chemistry, Kyoto University, Kyoto 606-8103, Japan
| |
Collapse
|
7
|
Kaupp M, Schattenberg CJ, Müller R, Reimann M. Unusually Large Effects of Charge-assisted C-H⋅⋅⋅F Hydrogen Bonds to Anionic Fluorine in Organic Solvents: Computational Study of 19 F NMR Shifts versus Thermochemistry. Chemistry 2022; 11:e202200146. [PMID: 35984672 PMCID: PMC9716039 DOI: 10.1002/open.202200146] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2022] [Revised: 07/26/2022] [Indexed: 01/31/2023]
Abstract
A comparison of computed 19 F NMR chemical shifts and experiment provides evidence for large specific solvent effects for fluoride-type anions interacting with the σ*(C-H) orbitals in organic solvents like MeCN or CH2 Cl2 . We show this for systems ranging from the fluoride ion and the bifluoride ion [FHF]- to polyhalogen anions [ClFx ]- . Discrepancies between computed and experimental shifts when using continuum solvent models like COSMO or force-field-based descriptions like the 3D-RISM-SCF model show specific orbital interactions that require a quantum-mechanical treatment of the solvent molecules. This is confirmed by orbital analyses of the shielding constants, while less negatively charged fluorine atoms (e. g., in [EF4 ]- ) do not require such quantum-mechanical treatments to achieve reasonable accuracy. The larger 19 F solvent shift of fluoride in MeCN compared to water is due to the larger coordination number in the former. These observations are due to unusually strong charge-assisted C-H⋅⋅⋅F- hydrogen bonds, which manifest beyond some threshold negative natural charge on fluorine of ca. < -0.6 e. The interactions are accompanied by sizable free energies of solvation, in the order F- ≫[FHF]- >[ClF2 ]- >[ClF4 ]- . COSMO-RS solvation free energies tend to moderately underestimate those from the micro-solvated cluster treatment. Red-shifted and intense vibrational C-H stretching bands, potentially accessible in bulk solution, are further spectroscopic finger prints.
Collapse
Affiliation(s)
- Martin Kaupp
- Technische Universität BerlinInstitut für Chemie, Theoretische Chemie/QuantenchemieSekr. C7, Strasse des 17. Juni 13510623BerlinGermany
| | - Caspar J. Schattenberg
- Technische Universität BerlinInstitut für Chemie, Theoretische Chemie/QuantenchemieSekr. C7, Strasse des 17. Juni 13510623BerlinGermany
| | - Robert Müller
- Technische Universität BerlinInstitut für Chemie, Theoretische Chemie/QuantenchemieSekr. C7, Strasse des 17. Juni 13510623BerlinGermany
| | - Marc Reimann
- Technische Universität BerlinInstitut für Chemie, Theoretische Chemie/QuantenchemieSekr. C7, Strasse des 17. Juni 13510623BerlinGermany
| |
Collapse
|
8
|
Weinreich J, Lemm D, von Rudorff GF, von Lilienfeld OA. Ab initio machine learning of phase space averages. J Chem Phys 2022; 157:024303. [PMID: 35840379 DOI: 10.1063/5.0095674] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Equilibrium structures determine material properties and biochemical functions. We here propose to machine learn phase space averages, conventionally obtained by ab initio or force-field-based molecular dynamics (MD) or Monte Carlo (MC) simulations. In analogy to ab initio MD, our ab initio machine learning (AIML) model does not require bond topologies and, therefore, enables a general machine learning pathway to obtain ensemble properties throughout the chemical compound space. We demonstrate AIML for predicting Boltzmann averaged structures after training on hundreds of MD trajectories. The AIML output is subsequently used to train machine learning models of free energies of solvation using experimental data and to reach competitive prediction errors (mean absolute error ∼ 0.8 kcal/mol) for out-of-sample molecules-within milliseconds. As such, AIML effectively bypasses the need for MD or MC-based phase space sampling, enabling exploration campaigns of Boltzmann averages throughout the chemical compound space at a much accelerated pace. We contextualize our findings by comparison to state-of-the-art methods resulting in a Pareto plot for the free energy of solvation predictions in terms of accuracy and time.
Collapse
Affiliation(s)
- Jan Weinreich
- Faculty of Physics, University of Vienna, Kolingasse 14-16, AT-1090 Wien, Austria
| | - Dominik Lemm
- Faculty of Physics, University of Vienna, Kolingasse 14-16, AT-1090 Wien, Austria
| | | | | |
Collapse
|
9
|
Reimann M, Kaupp M. Reaction Entropies in Solution from Analytical Three-Dimensional Reference Interaction Site Model Derivatives with Application to Redox and Spin-Crossover Processes. J Phys Chem A 2022; 126:3708-3716. [PMID: 35652546 DOI: 10.1021/acs.jpca.2c02317] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
An analytical approach to compute the excess entropy of solvation at constant pressure in three-dimensional reference interaction site model (3D-RISM) calculations is presented. It includes the changes in the macroscopic dielectric constant of the solvent upon variation of temperature and density. The approach is exact within the framework of force-field descriptions of the solute and gives reasonable results for self-consistently determined electrostatics as used in the 3D-RISM-self-consistent field approach, particularly for entropy differences. The new method is applied to simple examples of reaction entropies of iron complexes in aqueous solution, for which simple gas-phase calculations and many other approaches give unreliable estimates. For both redox half-reactions and spin-crossover processes, (semi)quantitative agreement with experimental reaction entropies can be achieved out of the box.
Collapse
Affiliation(s)
- Marc Reimann
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Technische Universität Berlin, Sekr. C7, Straße des 17.Juni 135, Berlin D-10623, Germany
| | - Martin Kaupp
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Technische Universität Berlin, Sekr. C7, Straße des 17.Juni 135, Berlin D-10623, Germany
| |
Collapse
|
10
|
Fuchs O, Trunschke S, Hanebrink H, Reimann M, Seitz O. Enabling Cysteine-Free Native Chemical Ligation at Challenging Junctions with a Ligation Auxiliary Capable of Base Catalysis. Angew Chem Int Ed Engl 2021; 60:19483-19490. [PMID: 34165893 PMCID: PMC8457107 DOI: 10.1002/anie.202107158] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2021] [Indexed: 12/14/2022]
Abstract
Ligation auxiliaries are used in chemical protein synthesis to extend the scope of native chemical ligation (NCL) beyond cysteine. However, auxiliary-mediated ligations at sterically demanding junctions have been difficult. Often the thioester intermediate formed in the thiol exchange step of NCL accumulates because the subsequent S→N acyl transfer is extremely slow. Here we introduce the 2-mercapto-2-(pyridin-2-yl)ethyl (MPyE) group as the first auxiliary designed to aid the ligation reaction by catalysis. Notably, the MPyE auxiliary provides useful rates even for junctions containing proline or a β-branched amino acid. Quantum chemical calculations suggest that the pyridine nitrogen acts as an intramolecular base in a rate-determining proton transfer step. The auxiliary is prepared in two steps and conveniently introduced by reductive alkylation. Auxiliary cleavage is induced upon treatment with TCEP/morpholine in presence of a MnII complex as radical starter. The synthesis of a de novo designed 99mer peptide and an 80 aa long MUC1 peptide demonstrates the usefulness of the MPyE auxiliary.
Collapse
Affiliation(s)
- Olaf Fuchs
- Department of Chemistry, Humboldt-Universität zu Berlin, Brook-Taylor-Str. 2, 12489, Berlin, Germany
| | - Sebastian Trunschke
- Department of Chemistry, Humboldt-Universität zu Berlin, Brook-Taylor-Str. 2, 12489, Berlin, Germany
| | - Hendrik Hanebrink
- Department of Chemistry, Humboldt-Universität zu Berlin, Brook-Taylor-Str. 2, 12489, Berlin, Germany
| | - Marc Reimann
- Department of Chemistry, Technische Universität Berlin, Straße des 17. Juni 135, 10623, Berlin, Germany
| | - Oliver Seitz
- Department of Chemistry, Humboldt-Universität zu Berlin, Brook-Taylor-Str. 2, 12489, Berlin, Germany
| |
Collapse
|
11
|
Ein zur Basenkatalyse befähigtes Ligationsauxiliar ermöglicht die cysteinfreie native chemische Ligation an anspruchsvollen Verknüpfungsstellen. Angew Chem Int Ed Engl 2021. [DOI: 10.1002/ange.202107158] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
|
12
|
Biomolecular Simulations with the Three-Dimensional Reference Interaction Site Model with the Kovalenko-Hirata Closure Molecular Solvation Theory. Int J Mol Sci 2021; 22:ijms22105061. [PMID: 34064655 PMCID: PMC8151972 DOI: 10.3390/ijms22105061] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2021] [Revised: 05/08/2021] [Accepted: 05/10/2021] [Indexed: 11/17/2022] Open
Abstract
The statistical mechanics-based 3-dimensional reference interaction site model with the Kovalenko-Hirata closure (3D-RISM-KH) molecular solvation theory has proven to be an essential part of a multiscale modeling framework, covering a vast region of molecular simulation techniques. The successful application ranges from the small molecule solvation energy to the bulk phase behavior of polymers, macromolecules, etc. The 3D-RISM-KH successfully predicts and explains the molecular mechanisms of self-assembly and aggregation of proteins and peptides related to neurodegeneration, protein-ligand binding, and structure-function related solvation properties. Upon coupling the 3D-RISM-KH theory with a novel multiple time-step molecular dynamic (MD) of the solute biomolecule stabilized by the optimized isokinetic Nosé-Hoover chain thermostat driven by effective solvation forces obtained from 3D-RISM-KH and extrapolated forward by generalized solvation force extrapolation (GSFE), gigantic outer time-steps up to picoseconds to accurately calculate equilibrium properties were obtained in this new quasidynamics protocol. The multiscale OIN/GSFE/3D-RISM-KH algorithm was implemented in the Amber package and well documented for fully flexible model of alanine dipeptide, miniprotein 1L2Y, and protein G in aqueous solution, with a solvent sampling rate ~150 times faster than a standard MD simulation in explicit water. Further acceleration in computation can be achieved by modifying the extent of solvation layers considered in the calculation, as well as by modifying existing closure relations. This enhanced simulation technique has proven applications in protein-ligand binding energy calculations, ligand/solvent binding site prediction, molecular solvation energy calculations, etc. Applications of the RISM-KH theory in molecular simulation are discussed in this work.
Collapse
|