1
|
Unke OT, Stöhr M, Ganscha S, Unterthiner T, Maennel H, Kashubin S, Ahlin D, Gastegger M, Medrano Sandonas L, Berryman JT, Tkatchenko A, Müller KR. Biomolecular dynamics with machine-learned quantum-mechanical force fields trained on diverse chemical fragments. SCIENCE ADVANCES 2024; 10:eadn4397. [PMID: 38579003 DOI: 10.1126/sciadv.adn4397] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/10/2023] [Accepted: 02/29/2024] [Indexed: 04/07/2024]
Abstract
The GEMS method enables molecular dynamics simulations of large heterogeneous systems at ab initio quality.
Collapse
Affiliation(s)
- Oliver T Unke
- Google DeepMind, Tucholskystraße 2, 10117 Berlin, Germany and Brandschenkestrasse 110, 8002 Zürich, Switzerland
- Machine Learning Group, Technische Universität Berlin, 10587 Berlin, Germany
- DFG Cluster of Excellence "Unifying Systems in Catalysis" (UniSysCat), Technische Universität Berlin, 10623 Berlin, Germany
| | - Martin Stöhr
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| | - Stefan Ganscha
- Google DeepMind, Tucholskystraße 2, 10117 Berlin, Germany and Brandschenkestrasse 110, 8002 Zürich, Switzerland
| | - Thomas Unterthiner
- Google DeepMind, Tucholskystraße 2, 10117 Berlin, Germany and Brandschenkestrasse 110, 8002 Zürich, Switzerland
| | - Hartmut Maennel
- Google DeepMind, Tucholskystraße 2, 10117 Berlin, Germany and Brandschenkestrasse 110, 8002 Zürich, Switzerland
| | - Sergii Kashubin
- Google DeepMind, Tucholskystraße 2, 10117 Berlin, Germany and Brandschenkestrasse 110, 8002 Zürich, Switzerland
| | - Daniel Ahlin
- Google DeepMind, Tucholskystraße 2, 10117 Berlin, Germany and Brandschenkestrasse 110, 8002 Zürich, Switzerland
| | - Michael Gastegger
- Machine Learning Group, Technische Universität Berlin, 10587 Berlin, Germany
- DFG Cluster of Excellence "Unifying Systems in Catalysis" (UniSysCat), Technische Universität Berlin, 10623 Berlin, Germany
- BASLEARN - TU Berlin/BASF Joint Lab for Machine Learning, Technische Universität Berlin, 10587 Berlin, Germany
| | - Leonardo Medrano Sandonas
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| | - Joshua T Berryman
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| | - Alexandre Tkatchenko
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| | - Klaus-Robert Müller
- Google DeepMind, Tucholskystraße 2, 10117 Berlin, Germany and Brandschenkestrasse 110, 8002 Zürich, Switzerland
- Machine Learning Group, Technische Universität Berlin, 10587 Berlin, Germany
- Department of Artificial Intelligence, Korea University, Anam-dong, Seongbuk-gu, Seoul 02841, Korea
- Max Planck Institute for Informatics, Stuhlsatzenhausweg, 66123 Saarbrücken, Germany
- BIFOLD - Berlin Institute for the Foundations of Learning and Data, Berlin, Germany
| |
Collapse
|
2
|
Otaki H, Ishiuchi SI, Fujii M, Sugita Y, Yagi K. Similarity scores of vibrational spectra reveal the atomistic structure of pentapeptides in multiple basins. Phys Chem Chem Phys 2024; 26:9906-9914. [PMID: 38477212 DOI: 10.1039/d4cp00064a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/14/2024]
Abstract
Vibrational spectroscopy combined with theoretical calculations is a powerful tool for analyzing the interaction and conformation of peptides at the atomistic level. Nonetheless, identifying the structure becomes increasingly difficult as the peptide size grows large. One example is acetyl-SIVSF-N-methylamide, a capped pentapeptide, whose atomistic structure has remained unknown since its first observation [T. Sekiguchi, M. Tamura, H. Oba, P. Çarçarbal, R. R. Lozada-Garcia, A. Zehnacker-Rentien, G. Grégoire, S. Ishiuchi and M. Fujii, Angew. Chem., Int. Ed., 2018, 57, 5626-5629]. Here, we propose a novel conformational search method, which exploits the structure-spectrum correlation using a similarity score that measures the agreement of theoretical and experimental spectra. Surprisingly, the two conformers have distinctly different energy and geometry. The second conformer is 25 kJ mol-1 higher in energy than the other, lowest-energy conformer. The result implies that there are multiple pathways in the early stage of the folding process: one to the global minimum and the other to a different basin. Once such a structure is established, the second conformer is unlikely to overcome the barrier to produce the most stable structure due to a vastly different hydrogen bond network of the backbone. Our proposed method can characterize the lowest-energy conformer and kinetically trapped, high-energy conformers of complex biomolecules.
Collapse
Affiliation(s)
- Hiroki Otaki
- Center for Bioinformatics and Molecular Medicine, Graduate School of Biomedical Sciences, Nagasaki University, 1-14 Bunkyo, Nagasaki, Nagasaki 852-8521, Japan
| | - Shun-Ichi Ishiuchi
- Department of Chemistry, School of Science, Tokyo Institute of Technology, 2-12-1 Ookayama, Meguro-ku, Tokyo, 152-8550, Japan
| | - Masaaki Fujii
- School of Life Science and Technology, Tokyo Institute of Technology, 4259 Nagatsuta-cho, Midori-ku, Yokohama, Kanagawa 226-8503, Japan
| | - Yuji Sugita
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan.
- Computational Biophysics Research Team, RIKEN Center for Computational Science, 7-1-26 Minatojima-Minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
- Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, 1-6-5 Minatojima-Minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| | - Kiyoshi Yagi
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan.
| |
Collapse
|
3
|
Schultz M, Parker SL, Fernando MT, Wellalage MM, Thomas DA. Diserinol Isophthalamide: A Novel Reagent for Complexation with Biomolecular Anions in Electrospray Ionization Mass Spectrometry. JOURNAL OF THE AMERICAN SOCIETY FOR MASS SPECTROMETRY 2023; 34:745-753. [PMID: 36975839 DOI: 10.1021/jasms.3c00010] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/18/2023]
Abstract
Transferring biomolecules from solution to vacuum facilitates a detailed analysis of molecular structure and dynamics by isolating molecules of interest from a complex environment. However, inherent in the ion desolvation process is the loss of solvent hydrogen bonding partners, which are critical for the stability of a condensed-phase structure. Thus, transfer of ions to vacuum can favor structural rearrangement, especially near solvent-accessible charge sites, which tend to adopt intramolecular hydrogen bonding motifs in the absence of solvent. Complexation of monoalkylammonium moieties (e.g., lysine side chains) with crown ethers such as 18-crown-6 can disfavor structural rearrangement of protonated sites, but no equivalent ligand has been investigated for deprotonated groups. Herein we describe diserinol isophthalamide (DIP), a novel reagent for the gas-phase complexation of anionic moieties within biomolecules. Complexation is observed to the C-terminus or side chains of the small model peptides GD, GE, GG, DF-OMe, VYV, YGGFL, and EYMPME in electrospray ionization mass spectrometry (ESI-MS) studies. In addition, complexation is observed with the phosphate and carboxylate moieities of phosphoserine and phosphotyrosine. DIP performs favorably in comparison to an existing anion recognition reagent, 1,1'-(1,2-phenylene)bis(3-phenylurea), that exhibits moderate carboxylate binding in organic solvent. This improved performance in ESI-MS experiments is attributed to reduced steric constraints to complexation with carboxylate groups of larger molecules. Overall, diserinol isophthalamide is an effective complexation reagent that can be applied in future work to study retention of solution-phase structure, investigate intrinsic molecular properties, and examine solvation effects.
Collapse
Affiliation(s)
- Madeline Schultz
- Department of Chemistry, University of Rhode Island, Kingston, Rhode Island 02881, United States
| | - Sarah L Parker
- Department of Chemistry, University of Rhode Island, Kingston, Rhode Island 02881, United States
| | - Maleesha T Fernando
- Department of Chemistry, University of Rhode Island, Kingston, Rhode Island 02881, United States
| | - Miyuru M Wellalage
- Department of Chemistry, University of Rhode Island, Kingston, Rhode Island 02881, United States
| | - Daniel A Thomas
- Department of Chemistry, University of Rhode Island, Kingston, Rhode Island 02881, United States
| |
Collapse
|
4
|
Kotobi A, Schwob L, Vonbun-Feldbauer GB, Rossi M, Gasparotto P, Feiler C, Berden G, Oomens J, Oostenrijk B, Scuderi D, Bari S, Meißner RH. Reconstructing the infrared spectrum of a peptide from representative conformers of the full canonical ensemble. Commun Chem 2023; 6:46. [PMID: 36869192 PMCID: PMC9984374 DOI: 10.1038/s42004-023-00835-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2022] [Accepted: 02/08/2023] [Indexed: 03/05/2023] Open
Abstract
Leucine enkephalin (LeuEnk), a biologically active endogenous opioid pentapeptide, has been under intense investigation because it is small enough to allow efficient use of sophisticated computational methods and large enough to provide insights into low-lying minima of its conformational space. Here, we reproduce and interpret experimental infrared (IR) spectra of this model peptide in gas phase using a combination of replica-exchange molecular dynamics simulations, machine learning, and ab initio calculations. In particular, we evaluate the possibility of averaging representative structural contributions to obtain an accurate computed spectrum that accounts for the corresponding canonical ensemble of the real experimental situation. Representative conformers are identified by partitioning the conformational phase space into subensembles of similar conformers. The IR contribution of each representative conformer is calculated from ab initio and weighted according to the population of each cluster. Convergence of the averaged IR signal is rationalized by merging contributions in a hierarchical clustering and the comparison to IR multiple photon dissociation experiments. The improvements achieved by decomposing clusters containing similar conformations into even smaller subensembles is strong evidence that a thorough assessment of the conformational landscape and the associated hydrogen bonding is a prerequisite for deciphering important fingerprints in experimental spectroscopic data.
Collapse
Affiliation(s)
- Amir Kotobi
- grid.7683.a0000 0004 0492 0453Deutsches Elektronen-Synchrotron DESY, Hamburg, Germany
| | - Lucas Schwob
- Deutsches Elektronen-Synchrotron DESY, Hamburg, Germany.
| | - Gregor B. Vonbun-Feldbauer
- grid.6884.20000 0004 0549 1777Hamburg University of Technology, Institute of Advanced Ceramics, Hamburg, Germany
| | - Mariana Rossi
- grid.469852.40000 0004 1796 3508Max Planck Institute for the Structure and Dynamics of Matter, Hamburg, Germany
| | - Piero Gasparotto
- grid.5991.40000 0001 1090 7501Scientific Computing Division, Paul Scherrer Institute, Villigen, Switzerland
| | - Christian Feiler
- grid.24999.3f0000 0004 0541 3699Helmholtz-Zentrum Hereon, Institute of Surface Science, Geesthacht, Germany
| | - Giel Berden
- grid.5590.90000000122931605Radboud University, Institute for Molecules and Materials, FELIX Laboratory, Nijmegen, The Netherlands
| | - Jos Oomens
- grid.5590.90000000122931605Radboud University, Institute for Molecules and Materials, FELIX Laboratory, Nijmegen, The Netherlands
| | - Bart Oostenrijk
- grid.7683.a0000 0004 0492 0453Deutsches Elektronen-Synchrotron DESY, Hamburg, Germany ,grid.9026.d0000 0001 2287 2617The Hamburg Centre for Ultrafast Imaging, Hamburg, Germany
| | - Debora Scuderi
- grid.503243.3Institut de Chimie Physique, CNRS UMR8000, Université Paris-Saclay, Orsay, France
| | - Sadia Bari
- Deutsches Elektronen-Synchrotron DESY, Hamburg, Germany. .,The Hamburg Centre for Ultrafast Imaging, Hamburg, Germany. .,Zernike Institute for Advanced Materials, University of Groningen, Groningen, The Netherlands.
| | - Robert H. Meißner
- grid.24999.3f0000 0004 0541 3699Helmholtz-Zentrum Hereon, Institute of Surface Science, Geesthacht, Germany ,grid.6884.20000 0004 0549 1777Hamburg University of Technology, Institute of Polymers and Composites, Hamburg, Germany
| |
Collapse
|
5
|
Villard J, Kılıç M, Rothlisberger U. Surrogate Based Genetic Algorithm Method for Efficient Identification of Low-Energy Peptide Structures. J Chem Theory Comput 2023; 19:1080-1097. [PMID: 36692853 PMCID: PMC9933449 DOI: 10.1021/acs.jctc.2c01078] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
Identification of the most stable structure(s) of a system is a prerequisite for the calculation of any of its properties from first-principles. However, even for relatively small molecules, exhaustive explorations of the potential energy surface (PES) are severely hampered by the dimensionality bottleneck. In this work, we address the challenging task of efficiently sampling realistic low-lying peptide coordinates by resorting to a surrogate based genetic algorithm (GA)/density functional theory (DFT) approach (sGADFT) in which promising candidates provided by the GA are ultimately optimized with DFT. We provide a benchmark of several computational methods (GAFF, AMOEBApro13, PM6, PM7, DFTB3-D3(BJ)) as possible prescanning surrogates and apply sGADFT to two test case systems that are (i) two isomer families of the protonated Gly-Pro-Gly-Gly tetrapeptide (Masson, A.; J. Am. Soc. Mass Spectrom.2015, 26, 1444-1454) and (ii) the doubly protonated cyclic decapeptide gramicidin S (Nagornova, N. S.; J. Am. Chem. Soc.2010, 132, 4040-4041). We show that our GA procedure can correctly identify low-energy minima in as little as a few hours. Subsequent refinement of surrogate low-energy structures within a given energy threshold (≤10 kcal/mol (i), ≤5 kcal/mol (ii)) via DFT relaxation invariably led to the identification of the most stable structures as determined from high-resolution infrared (IR) spectroscopy at low temperature. The sGADFT method therefore constitutes a highly efficient route for the screening of realistic low-lying peptide structures in the gas phase as needed for instance for the interpretation and assignment of experimental IR spectra.
Collapse
|
6
|
Lin K, TomHon P, Lehmkuhl S, Laasner R, Theis T, Blum V. Density Functional Theory Study of Reaction Equilibria in Signal Amplification by Reversible Exchange. Chemphyschem 2021; 22:1947-1957. [PMID: 34549869 DOI: 10.1002/cphc.202100204] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2021] [Revised: 05/19/2021] [Indexed: 11/07/2022]
Abstract
An in-depth theoretical analysis of key chemical equilibria in Signal Amplification by Reversible Exchange (SABRE) is provided, employing density functional theory calculations to characterize the likely reaction network. For all reactions in the network, the potential energy surface is probed to identify minimum energy pathways. Energy barriers and transition states are calculated, and harmonic transition state theory is applied to calculate exchange rates that approximate experimental values. The reaction network energy surface can be modulated by chemical potentials that account for the dependence on concentration, temperature, and partial pressure of molecular constituents (hydrogen, methanol, pyridine) supplied to the experiment under equilibrium conditions. We show that, under typical experimental conditions, the Gibbs free energies of the two key states involved in pyridine-hydrogen exchange at the common Ir-IMes catalyst system in methanol are essentially the same, i. e., nearly optimal for SABRE. We also show that a methanol-containing intermediate is plausible as a transient species in the process.
Collapse
Affiliation(s)
- Kailai Lin
- Department of Chemistry, Duke University, Durham, NC 27708, USA
| | - Patrick TomHon
- Department of Chemistry, North Carolina State University, Raleigh, NC 27606, USA
| | - Sören Lehmkuhl
- Department of Chemistry, North Carolina State University, Raleigh, NC 27606, USA
| | - Raul Laasner
- Thomas Lord Department of Mechanical Engineering and Materials Science, Duke University, Durham, NC 27708, USA
| | - Thomas Theis
- Department of Chemistry, North Carolina State University, Raleigh, NC 27606, USA.,Joint Department of Biomedical Engineering, UNC, Chapel Hill, and NC State University, Raleigh, NC 27606, USA.,Department of Physics, North Carolina State University, Raleigh, NC 27606, USA
| | - Volker Blum
- Department of Chemistry, Duke University, Durham, NC 27708, USA.,Thomas Lord Department of Mechanical Engineering and Materials Science, Duke University, Durham, NC 27708, USA
| |
Collapse
|
7
|
Adasme-Carreño F, Caballero J, Ireta J. PSIQUE: Protein Secondary Structure Identification on the Basis of Quaternions and Electronic Structure Calculations. J Chem Inf Model 2021; 61:1789-1800. [PMID: 33769809 DOI: 10.1021/acs.jcim.0c01343] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The secondary structure is important in protein structure analysis, classification, and modeling. We have developed a novel method for secondary structure assignment, termed PSIQUE, based on the potential energy surface (PES) of polyalanine obtained using an infinitely long chain model and density functional theory calculations. First, uniform protein segments are determined in terms of a difference of quaternions between neighboring amino acids along the protein backbone. Then, the identification of the secondary structure motifs is carried out based on the minima found in the PES. PSIQUE shows good agreement with other secondary structure assignment methods. However, it provides better discrimination of subtle secondary structures (e.g., helix types) and termini and produces more uniform segments while also accounting for local distortions. Overall, PSIQUE provides a precise and reliable assignment of secondary structures, so it should be helpful for the detailed characterization of the protein structure.
Collapse
Affiliation(s)
- Francisco Adasme-Carreño
- Departamento de Bioinformática, Centro de Bioinformática, Simulación y Modelado (CBSM), Facultad de Ingeniería, Universidad de Talca, Campus Talca, 1 Poniente No. 1141, Casilla 721, Talca 3460000, Chile
| | - Julio Caballero
- Departamento de Bioinformática, Centro de Bioinformática, Simulación y Modelado (CBSM), Facultad de Ingeniería, Universidad de Talca, Campus Talca, 1 Poniente No. 1141, Casilla 721, Talca 3460000, Chile
| | - Joel Ireta
- Departamento de Química, División de Ciencias Básicas e Ingeniería, Universidad Autónoma Metropolitana-Iztapalapa, A.P. 55-534, Ciudad de Mexico 09340, Mexico
| |
Collapse
|
8
|
Colell JFP, Logan AWJ, Zhou Z, Lindale JR, Laasner R, Shchepin RV, Chekmenev EY, Blum V, Warren WS, Malcolmson SJ, Theis T. Rational ligand choice extends the SABRE substrate scope. Chem Commun (Camb) 2020; 56:9336-9339. [PMID: 32671356 PMCID: PMC7443256 DOI: 10.1039/d0cc01330g] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Here we report on chelating ligands for Signal Amplification By Reversible Exchange (SABRE) catalysts that permit hyperpolarisation on otherwise sterically hindered substrates. We demonstrate 1H enhancements of ∼100-fold over 8.5 T thermal for 2-substituted pyridines, and smaller, yet significant enhancements for provitamin B6 and caffeine. We also show 15N-enhancements of ∼1000-fold and 19F-enhancements of 30-fold.
Collapse
Affiliation(s)
| | | | - Zijian Zhou
- Department of Chemistry, Duke University, Durham, NC 27708, USA
| | | | - Raul Laasner
- Department of Mechanical Engineering and Materials Science, Duke University, Durham, NC 27708, USA
| | - Roman V. Shchepin
- Department of Chemistry, Biology, and Health Sciences, South Dakota School of Mines and Technology, Rapid City, South Dakota 57701, USA
| | - Eduard Y. Chekmenev
- Russian Academy of Sciences, Leninskiy Prospekt 14, 119991 Moscow, Russia
- Department of Chemistry, Integrative Biosciences (IBio), Karmanos Cancer Institute (KCI), Wayne State University, Detroit, MI 48202, USA
| | - Volker Blum
- Department of Chemistry, Duke University, Durham, NC 27708, USA
- Department of Mechanical Engineering and Materials Science, Duke University, Durham, NC 27708, USA
| | - Warren S. Warren
- Department of Chemistry, Duke University, Durham, NC 27708, USA
- Departments of Physics, Radiology and Biomedical Engineering, Duke University, Durham, NC 27707, USA
| | | | - Thomas Theis
- Department of Chemistry, Duke University, Durham, NC 27708, USA
- Department of Chemistry, North Carolina State University, Raleigh, NC 27695
- Joint Department of Biomedical Engineering University of North Carolina at Chapel Hill and North Carolina State University, Chapel Hill, NC, USA
| |
Collapse
|
9
|
Bravo MF, Palanichamy K, Shlain MA, Schiro F, Naeem Y, Marianski M, Braunschweig AB. Synthesis and Binding of Mannose‐Specific Synthetic Carbohydrate Receptors. Chemistry 2020; 26:11782-11795. [DOI: 10.1002/chem.202000481] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/28/2020] [Revised: 03/15/2020] [Indexed: 12/16/2022]
Affiliation(s)
- M. Fernando Bravo
- Advanced Science Research Center at the Graduate Center City University of New York 85 St Nicholas Terrace New York NY 10031 USA
- Department of Chemistry and Biochemistry Hunter College 695 Park Ave New York NY 10065 USA
- The PhD Program in Chemistry The Graduate Center of the, City University of New York 365 5th Ave New York NY 10016 USA
| | - Kalanidhi Palanichamy
- Advanced Science Research Center at the Graduate Center City University of New York 85 St Nicholas Terrace New York NY 10031 USA
- Department of Chemistry and Biochemistry Hunter College 695 Park Ave New York NY 10065 USA
| | - Milan A. Shlain
- Advanced Science Research Center at the Graduate Center City University of New York 85 St Nicholas Terrace New York NY 10031 USA
- Department of Chemistry and Biochemistry Hunter College 695 Park Ave New York NY 10065 USA
| | - Frank Schiro
- Advanced Science Research Center at the Graduate Center City University of New York 85 St Nicholas Terrace New York NY 10031 USA
- Department of Chemistry and Biochemistry Hunter College 695 Park Ave New York NY 10065 USA
| | - Yasir Naeem
- Advanced Science Research Center at the Graduate Center City University of New York 85 St Nicholas Terrace New York NY 10031 USA
- Department of Chemistry and Biochemistry Hunter College 695 Park Ave New York NY 10065 USA
| | - Mateusz Marianski
- Advanced Science Research Center at the Graduate Center City University of New York 85 St Nicholas Terrace New York NY 10031 USA
- Department of Chemistry and Biochemistry Hunter College 695 Park Ave New York NY 10065 USA
- The PhD Program in Chemistry The Graduate Center of the, City University of New York 365 5th Ave New York NY 10016 USA
- The PhD Program in Biochemistry The Graduate Center of the, City University of New York 365 5th Ave New York NY 10016 USA
| | - Adam B. Braunschweig
- Advanced Science Research Center at the Graduate Center City University of New York 85 St Nicholas Terrace New York NY 10031 USA
- Department of Chemistry and Biochemistry Hunter College 695 Park Ave New York NY 10065 USA
- The PhD Program in Chemistry The Graduate Center of the, City University of New York 365 5th Ave New York NY 10016 USA
- The PhD Program in Biochemistry The Graduate Center of the, City University of New York 365 5th Ave New York NY 10016 USA
| |
Collapse
|
10
|
Thomas DA, Chang R, Mucha E, Lettow M, Greis K, Gewinner S, Schöllkopf W, Meijer G, von Helden G. Probing the conformational landscape and thermochemistry of DNA dinucleotide anions via helium nanodroplet infrared action spectroscopy. Phys Chem Chem Phys 2020; 22:18400-18413. [PMID: 32797142 DOI: 10.1039/d0cp02482a] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Isolation of biomolecules in vacuum facilitates characterization of the intramolecular interactions that determine three-dimensional structure, but experimental quantification of conformer thermochemistry remains challenging. Infrared spectroscopy of molecules trapped in helium nanodroplets is a promising methodology for the measurement of thermochemical parameters. When molecules are captured in a helium nanodroplet, the rate of cooling to an equilibrium temperature of ca. 0.4 K is generally faster than the rate of isomerization, resulting in "shock-freezing" that kinetically traps molecules in local conformational minima. This unique property enables the study of temperature-dependent conformational equilibria via infrared spectroscopy at 0.4 K, thereby avoiding the deleterious effects of spectral broadening at higher temperatures. Herein, we demonstrate the first application of this approach to ionic species by coupling electrospray ionization mass spectrometry (ESI-MS) with helium nanodroplet infrared action spectroscopy to probe the structure and thermochemistry of deprotonated DNA dinucleotides. Dinucleotide anions were generated by ESI, confined in an ion trap at temperatures between 90 and 350 K, and entrained in traversing helium nanodroplets. The infrared action spectra of the entrained ions show a strong dependence on pre-pickup ion temperature, consistent with the preservation of conformer population upon cooling to 0.4 K. Non-negative matrix factorization was utilized to identify component conformer infrared spectra and determine temperature-dependent conformer populations. Relative enthalpies and entropies of conformers were subsequently obtained from a van't Hoff analysis. IR spectra and conformer thermochemistry are compared to results from ion mobility spectrometry (IMS) and electronic structure methods. The implementation of ESI-MS as a source of dopant molecules expands the diversity of molecules accessible for thermochemical measurements, enabling the study of larger, non-volatile species.
Collapse
Affiliation(s)
- Daniel A Thomas
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany.
| | | | | | | | | | | | | | | | | |
Collapse
|
11
|
Bannwarth C, Caldeweyher E, Ehlert S, Hansen A, Pracht P, Seibert J, Spicher S, Grimme S. Extended
tight‐binding
quantum chemistry methods. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2020. [DOI: 10.1002/wcms.1493] [Citation(s) in RCA: 218] [Impact Index Per Article: 54.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Christoph Bannwarth
- Department of Chemistry and The PULSE Institute Stanford University Stanford California USA
| | - Eike Caldeweyher
- Mulliken Center for Theoretical Chemistry Rheinische Friedrich‐Wilhelms‐Universität Bonn Bonn Germany
| | - Sebastian Ehlert
- Mulliken Center for Theoretical Chemistry Rheinische Friedrich‐Wilhelms‐Universität Bonn Bonn Germany
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry Rheinische Friedrich‐Wilhelms‐Universität Bonn Bonn Germany
| | - Philipp Pracht
- Mulliken Center for Theoretical Chemistry Rheinische Friedrich‐Wilhelms‐Universität Bonn Bonn Germany
| | - Jakob Seibert
- Mulliken Center for Theoretical Chemistry Rheinische Friedrich‐Wilhelms‐Universität Bonn Bonn Germany
| | - Sebastian Spicher
- Mulliken Center for Theoretical Chemistry Rheinische Friedrich‐Wilhelms‐Universität Bonn Bonn Germany
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry Rheinische Friedrich‐Wilhelms‐Universität Bonn Bonn Germany
| |
Collapse
|
12
|
Warnke S, Ben Faleh A, Pellegrinelli RP, Yalovenko N, Rizzo TR. Combining ultra-high resolution ion mobility spectrometry with cryogenic IR spectroscopy for the study of biomolecular ions. Faraday Discuss 2020; 217:114-125. [PMID: 30993271 PMCID: PMC6657637 DOI: 10.1039/c8fd00180d] [Citation(s) in RCA: 37] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
We explore the capability of SLIM-based IMS for isomer selectivity in combination with cryogenic, messenger-tagging IR spectroscopy.
Double-resonance spectroscopic schemes in combination with cryogenic ion traps are the go-to techniques when isomer-specific high-resolution spectra are required for analysis of molecular ions. Their limitation lies in the requirement for well-resolved, isomer-specific absorption bands as well as in the potentially time-consuming steps to identify each isomer. We present an alternative approach where isomeric species are readily separated using ion mobility spectrometry (IMS) and selected prior to cryogenic spectroscopic analysis. To date, most IMS approaches suffer from relatively low resolution, however, recent technological developments in the field of travelling-wave ion mobility using structures for lossless ion manipulation (SLIM) permit the use of extremely long drift paths, which greatly enhances the resolution. We demonstrate the power of combining this type of ultra-high resolution IMS with cryogenic vibrational spectroscopy by comparing mobility-resolved IR spectra of a disaccharide to those acquired using IR–IR double resonance. This new approach is especially promising for the investigation of larger molecules where spectral congestion interferes with double resonance techniques.
Collapse
Affiliation(s)
- Stephan Warnke
- Laboratoire de Chimie Physique Moléculaire, École Polytechnique Fédérale de Lausanne, EPFL SB ISIC LCPM, Station 6, CH-1015 Lausanne, Switzerland.
| | - Ahmed Ben Faleh
- Laboratoire de Chimie Physique Moléculaire, École Polytechnique Fédérale de Lausanne, EPFL SB ISIC LCPM, Station 6, CH-1015 Lausanne, Switzerland.
| | - Robert P Pellegrinelli
- Laboratoire de Chimie Physique Moléculaire, École Polytechnique Fédérale de Lausanne, EPFL SB ISIC LCPM, Station 6, CH-1015 Lausanne, Switzerland.
| | - Natalia Yalovenko
- Laboratoire de Chimie Physique Moléculaire, École Polytechnique Fédérale de Lausanne, EPFL SB ISIC LCPM, Station 6, CH-1015 Lausanne, Switzerland.
| | - Thomas R Rizzo
- Laboratoire de Chimie Physique Moléculaire, École Polytechnique Fédérale de Lausanne, EPFL SB ISIC LCPM, Station 6, CH-1015 Lausanne, Switzerland.
| |
Collapse
|
13
|
Janke SM, Qarai MB, Blum V, Spano FC. Frenkel-Holstein Hamiltonian applied to absorption spectra of quaterthiophene-based 2D hybrid organic-inorganic perovskites. J Chem Phys 2020; 152:144702. [PMID: 32295353 DOI: 10.1063/1.5139044] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
For the prototypical two-dimensional hybrid organic-inorganic perovskites (2D HOIPs) (AE4T)PbX4 (X = Cl, Br, and I), we demonstrate that the Frenkel-Holstein Hamiltonian (FHH) can be applied to describe the absorption spectrum arising from the organic component. We first model the spectra using only the four nearest neighbor couplings between translationally inequivalent molecules in the organic herringbone lattice as fitting parameters in the FHH. We next use linear-response time-dependent density functional theory (LR-TDDFT) to calculate molecular transition densities, from which extended excitonic couplings are evaluated based on the atomic positions within the 2D HOIPs. We find that both approaches reproduce the experimentally observed spectra, including changes in their shape and peak positions. The spectral changes are correlated with a decrease in excitonic coupling from X = Cl to X = I. Importantly, the LR-TDDFT-based approach with extended excitonic couplings not only gives better agreement with the experimental absorption line shape than the approach using a restricted set of fitted parameters but also allows us to relate the changes in excitonic coupling to the underlying geometry. We accordingly find that the decrease in excitonic coupling from X = Cl to Br to I is due to an increase in molecular separation, which in turn can be related to the increasing Pb-X bond length from Cl to I. Our research opens up a potential pathway to predicting optoelectronic properties of new 2D HOIPs from ab initio calculations and to gain insight into structural relations from 2D HOIP absorption spectra.
Collapse
Affiliation(s)
- Svenja M Janke
- Department of Mechanical Engineering and Materials Science, Duke University, Durham, North Carolina 27708, USA
| | - Mohammad B Qarai
- Chemistry, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Volker Blum
- Department of Mechanical Engineering and Materials Science, Duke University, Durham, North Carolina 27708, USA
| | - Frank C Spano
- Chemistry, Temple University, Philadelphia, Pennsylvania 19122, USA
| |
Collapse
|
14
|
Pracht P, Bohle F, Grimme S. Automated exploration of the low-energy chemical space with fast quantum chemical methods. Phys Chem Chem Phys 2020; 22:7169-7192. [PMID: 32073075 DOI: 10.1039/c9cp06869d] [Citation(s) in RCA: 810] [Impact Index Per Article: 202.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
We propose and discuss an efficient scheme for the in silico sampling for parts of the molecular chemical space by semiempirical tight-binding methods combined with a meta-dynamics driven search algorithm. The focus of this work is set on the generation of proper thermodynamic ensembles at a quantum chemical level for conformers, but similar procedures for protonation states, tautomerism and non-covalent complex geometries are also discussed. The conformational ensembles consisting of all significantly populated minimum energy structures normally form the basis of further, mostly DFT computational work, such as the calculation of spectra or macroscopic properties. By using basic quantum chemical methods, electronic effects or possible bond breaking/formation are accounted for and a very reasonable initial energetic ranking of the candidate structures is obtained. Due to the huge computational speedup gained by the fast low-cost quantum chemical methods, overall short computation times even for systems with hundreds of atoms (typically drug-sized molecules) are achieved. Furthermore, specialized applications, such as sampling with implicit solvation models or constrained conformational sampling for transition-states, metal-, surface-, or noncovalently bound complexes are discussed, opening many possible applications in modern computational chemistry and drug discovery. The procedures have been implemented in a freely available computer code called CREST, that makes use of the fast and reliable GFNn-xTB methods.
Collapse
Affiliation(s)
- Philipp Pracht
- Mulliken Center for Theoretical Chemistry, Universität Bonn, Beringstr. 4, 53115 Bonn, Germany.
| | | | | |
Collapse
|
15
|
Abstract
Estimating the range of three-dimensional structures (conformations) that are available to a molecule is a key component of computer-aided drug design. Quantum mechanical simulation offers improved accuracy over forcefield methods, but at a high computational cost. The question is whether this increased cost can be justified in a context in which high-throughput analysis of large numbers of molecules is often key. This chapter discusses the application of quantum mechanics to conformational searching, with a focus on three key challenges: (1) the generation of ensembles that include a good approximation to a molecule's bioactive conformation at as prominent a ranking as possible; (2) rational analysis and modification of a pre-established bioactive conformation in terms of its energetics; and (3) approximation of real solution-phase conformational ensembles in tandem with NMR data. The impact of QM on the high-throughput application (1) is debatable, meaning that for the moment its primary application is still lower-throughput applications such as (2) and (3). The optimal choice of QM method is also discussed. Rigorous benchmarking suggests that DFT methods are only acceptable when used with large basis sets, but a trickle of papers continue to obtain useful results with relatively low-cost methods, leading to a dilemma that the literature has yet to fully resolve.
Collapse
|
16
|
Stöhr M, Tkatchenko A. Quantum mechanics of proteins in explicit water: The role of plasmon-like solute-solvent interactions. SCIENCE ADVANCES 2019; 5:eaax0024. [PMID: 31853494 PMCID: PMC6910842 DOI: 10.1126/sciadv.aax0024] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/13/2019] [Accepted: 09/18/2019] [Indexed: 05/05/2023]
Abstract
Quantum-mechanical van der Waals dispersion interactions play an essential role in intraprotein and protein-water interactions-the two main factors affecting the structure and dynamics of proteins in water. Typically, these interactions are only treated phenomenologically, via pairwise potential terms in classical force fields. Here, we use an explicit quantum-mechanical approach of density-functional tight-binding combined with the many-body dispersion formalism and demonstrate the relevance of many-body van der Waals forces both to protein energetics and to protein-water interactions. In contrast to commonly used pairwise approaches, many-body effects substantially decrease the relative stability of native states in the absence of water. Upon solvation, the protein-water dispersion interaction counteracts this effect and stabilizes native conformations and transition states. These observations arise from the highly delocalized and collective character of the interactions, suggesting a remarkable persistence of electron correlation through aqueous environments and providing the basis for long-range interaction mechanisms in biomolecular systems.
Collapse
|
17
|
Stöhr M, Van Voorhis T, Tkatchenko A. Theory and practice of modeling van der Waals interactions in electronic-structure calculations. Chem Soc Rev 2019; 48:4118-4154. [PMID: 31190037 DOI: 10.1039/c9cs00060g] [Citation(s) in RCA: 68] [Impact Index Per Article: 13.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
The accurate description of long-range electron correlation, most prominently including van der Waals (vdW) dispersion interactions, represents a particularly challenging task in the modeling of molecules and materials. vdW forces arise from the interaction of quantum-mechanical fluctuations in the electronic charge density. Within (semi-)local density functional approximations or Hartree-Fock theory such interactions are neglected altogether. Non-covalent vdW interactions, however, are ubiquitous in nature and play a key role for the understanding and accurate description of the stability, dynamics, structure, and response properties in a plethora of systems. During the last decade, many promising methods have been developed for modeling vdW interactions in electronic-structure calculations. These methods include vdW-inclusive Density Functional Theory and correlated post-Hartree-Fock approaches. Here, we focus on the methods within the framework of Density Functional Theory, including non-local van der Waals density functionals, interatomic dispersion models within many-body and pairwise formulation, and random phase approximation-based approaches. This review aims to guide the reader through the theoretical foundations of these methods in a tutorial-style manner and, in particular, highlight practical aspects such as the applicability and the advantages and shortcomings of current vdW-inclusive approaches. In addition, we give an overview of complementary experimental approaches, and discuss tools for the qualitative understanding of non-covalent interactions as well as energy decomposition techniques. Besides representing a reference for the current state-of-the-art, this work is thus also designed as a concise and detailed introduction to vdW-inclusive electronic structure calculations for a general and broad audience.
Collapse
Affiliation(s)
- Martin Stöhr
- Physics and Materials Science Research Unit, University of Luxembourg, L-1511 Luxembourg, Luxembourg.
| | | | | |
Collapse
|
18
|
Scutelnic V, Perez MAS, Marianski M, Warnke S, Gregor A, Rothlisberger U, Bowers MT, Baldauf C, von Helden G, Rizzo TR, Seo J. The Structure of the Protonated Serine Octamer. J Am Chem Soc 2018; 140:7554-7560. [DOI: 10.1021/jacs.8b02118] [Citation(s) in RCA: 48] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Valeriu Scutelnic
- Laboratory of Molecular Physical Chemistry, Ecole Polytechnique Fédérale de Lausanne, Station 6, CH-1015 Lausanne, Switzerland
| | - Marta A. S. Perez
- Laboratory of Computational Chemistry and Biochemistry, Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| | - Mateusz Marianski
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| | - Stephan Warnke
- Laboratory of Molecular Physical Chemistry, Ecole Polytechnique Fédérale de Lausanne, Station 6, CH-1015 Lausanne, Switzerland
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| | - Aurelien Gregor
- Laboratory of Computational Chemistry and Biochemistry, Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| | - Ursula Rothlisberger
- Laboratory of Computational Chemistry and Biochemistry, Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| | - Michael T. Bowers
- Department of Chemistry and Biochemistry, University of California Santa Barbara, Santa Barbara, California 93106, United States
| | - Carsten Baldauf
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| | - Gert von Helden
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| | - Thomas R. Rizzo
- Laboratory of Molecular Physical Chemistry, Ecole Polytechnique Fédérale de Lausanne, Station 6, CH-1015 Lausanne, Switzerland
| | - Jongcheol Seo
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| |
Collapse
|
19
|
Hoffmann W, Marianski M, Warnke S, Seo J, Baldauf C, von Helden G, Pagel K. Assessing the stability of alanine-based helices by conformer-selective IR spectroscopy. Phys Chem Chem Phys 2018; 18:19950-4. [PMID: 27398722 DOI: 10.1039/c6cp03827a] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Polyalanine based peptides that carry a lysine at the C-terminus ([Ac-AlanLys + H](+)) are known to form α-helices in the gas phase. Three factors contribute to the stability of these helices: (i) the interaction between the helix macro dipole and the charge, (ii) the capping of dangling C[double bond, length as m-dash]O groups by lysine and (iii) the cooperative hydrogen bond network. In previous studies, the influence of the interaction between the helix dipole and the charge as well as the impact of the capping was studied intensively. Here, we complement these findings by systematically assessing the third parameter, the H-bond network. In order to selectively remove one H-bond along the backbone, we use amide-to-ester substitutions. The resulting depsi peptides were analyzed by ion-mobility and m/z-selective infrared spectroscopy as well as theoretical calculations. Our results indicate that peptides which contain only one ester bond still maintain the helical conformation. We conclude that the interaction between the charge and the helix macro-dipole is most crucial for the formation of the α-helical conformation and a single backbone H-bond has only little influence on the overall stability.
Collapse
Affiliation(s)
- Waldemar Hoffmann
- Freie Universität Berlin, Takustraße 3, 14195 Berlin, Germany. and Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| | - Mateusz Marianski
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| | - Stephan Warnke
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| | - Jongcheol Seo
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| | - Carsten Baldauf
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| | - Gert von Helden
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| | - Kevin Pagel
- Freie Universität Berlin, Takustraße 3, 14195 Berlin, Germany. and Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| |
Collapse
|
20
|
Bocklitz S, Suhm MA. Polymer Segments at the Folding Limit: Raman Scattering for the Diglyme Benchmark. Chemphyschem 2017; 18:3570-3575. [DOI: 10.1002/cphc.201701169] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/26/2017] [Indexed: 11/11/2022]
Affiliation(s)
- Sebastian Bocklitz
- Institut für Physikalische Chemie; Georg-August-Universität Göttingen; Tammannstraße 6 37077 Göttingen Germany
| | - Martin A. Suhm
- Institut für Physikalische Chemie; Georg-August-Universität Göttingen; Tammannstraße 6 37077 Göttingen Germany
| |
Collapse
|
21
|
Haeffner F, Irikura KK. N-Protonated Isomers and Coulombic Barriers to Dissociation of Doubly Protonated Ala 8Arg. JOURNAL OF THE AMERICAN SOCIETY FOR MASS SPECTROMETRY 2017; 28:2170-2180. [PMID: 28699065 DOI: 10.1007/s13361-017-1719-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/07/2016] [Revised: 05/15/2017] [Accepted: 05/15/2017] [Indexed: 06/07/2023]
Abstract
Collision-induced dissociation (or tandem mass spectrometry, MS/MS) of a protonated peptide results in a spectrum of fragment ions that is useful for inferring amino acid sequence. This is now commonplace and a foundation of proteomics. The underlying chemical and physical processes are believed to be those familiar from physical organic chemistry and chemical kinetics. However, first-principles predictions remain intractable because of the conflicting necessities for high accuracy (to achieve qualitatively correct kinetics) and computational speed (to compensate for the high cost of reliable calculations on such large molecules). To make progress, shortcuts are needed. Inspired by the popular mobile proton model, we have previously proposed a simplified theoretical model in which the gas-phase fragmentation pattern of protonated peptides reflects the relative stabilities of N-protonated isomers, thus avoiding the need for transition-state information. For singly protonated Ala n (n = 3-11), the resulting predictions were in qualitative agreement with the results from low-energy MS/MS experiments. Here, the comparison is extended to a model tryptic peptide, doubly protonated Ala8Arg. This is of interest because doubly protonated tryptic peptides are the most important in proteomics. In comparison with experimental results, our model seriously overpredicts the degree of backbone fragmentation at N9. We offer an improved model that corrects this deficiency. The principal change is to include Coulombic barriers, which hinder the separation of the product cations from each other. Coulombic barriers may be equally important in MS/MS of all multiply charged peptide ions. Graphical Abstract ᅟ.
Collapse
Affiliation(s)
- Fredrik Haeffner
- Chemical Sciences Division, National Institute of Standards and Technology, Gaithersburg, MD, 20899-8320, USA
- Department of Chemistry, Boston College, 2609 Beacon Street, Chestnut Hill, MA, 02467-3860, USA
| | - Karl K Irikura
- Chemical Sciences Division, National Institute of Standards and Technology, Gaithersburg, MD, 20899-8320, USA.
| |
Collapse
|
22
|
Schneider M, Masellis C, Rizzo T, Baldauf C. Kinetically Trapped Liquid-State Conformers of a Sodiated Model Peptide Observed in the Gas Phase. J Phys Chem A 2017; 121:6838-6844. [DOI: 10.1021/acs.jpca.7b06431] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Affiliation(s)
- Markus Schneider
- Theory
Department, Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany
| | - Chiara Masellis
- Laboratoire
de Chimie Physique Moléculaire, EPFL SB ISIC LCPM, Ecole Polytechnique Fédérale de Lausanne, Station 6, CH-1015 Lausanne, Switzerland
| | - Thomas Rizzo
- Laboratoire
de Chimie Physique Moléculaire, EPFL SB ISIC LCPM, Ecole Polytechnique Fédérale de Lausanne, Station 6, CH-1015 Lausanne, Switzerland
| | - Carsten Baldauf
- Theory
Department, Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany
| |
Collapse
|
23
|
Supady A, Hecht S, Baldauf C. About Underappreciated Yet Active Conformations of Thiourea Organocatalysts. Org Lett 2017; 19:4199-4202. [PMID: 28786673 DOI: 10.1021/acs.orglett.7b01782] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Conformational dynamics can define the function of organocatalysts. While the accepted mechanism of Schreiner's catalyst features a double hydrogen bond to the substrate that only forms with the anti-anti conformation of its central thiourea group, our electronic-structure theory study reveals that binding of the model substrate methyl vinyl ketone prefers syn-anti conformations. We find a new mechanism featuring π stacking interactions and highlight the need for extensive structure searches for flexible molecules, especially when aiming for structure-based design of catalytic activity.
Collapse
Affiliation(s)
- Adriana Supady
- Fritz-Haber-Institut der Max-Planck-Gesellschaft , 14195 Berlin, Germany
| | - Stefan Hecht
- Department of Chemistry, Humboldt-Universität zu Berlin , 12489 Berlin, Germany
| | - Carsten Baldauf
- Fritz-Haber-Institut der Max-Planck-Gesellschaft , 14195 Berlin, Germany
| |
Collapse
|
24
|
Hermann J, DiStasio RA, Tkatchenko A. First-Principles Models for van der Waals Interactions in Molecules and Materials: Concepts, Theory, and Applications. Chem Rev 2017; 117:4714-4758. [PMID: 28272886 DOI: 10.1021/acs.chemrev.6b00446] [Citation(s) in RCA: 255] [Impact Index Per Article: 36.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Abstract
Noncovalent van der Waals (vdW) or dispersion forces are ubiquitous in nature and influence the structure, stability, dynamics, and function of molecules and materials throughout chemistry, biology, physics, and materials science. These forces are quantum mechanical in origin and arise from electrostatic interactions between fluctuations in the electronic charge density. Here, we explore the conceptual and mathematical ingredients required for an exact treatment of vdW interactions, and present a systematic and unified framework for classifying the current first-principles vdW methods based on the adiabatic-connection fluctuation-dissipation (ACFD) theorem (namely the Rutgers-Chalmers vdW-DF, Vydrov-Van Voorhis (VV), exchange-hole dipole moment (XDM), Tkatchenko-Scheffler (TS), many-body dispersion (MBD), and random-phase approximation (RPA) approaches). Particular attention is paid to the intriguing nature of many-body vdW interactions, whose fundamental relevance has recently been highlighted in several landmark experiments. The performance of these models in predicting binding energetics as well as structural, electronic, and thermodynamic properties is connected with the theoretical concepts and provides a numerical summary of the state-of-the-art in the field. We conclude with a roadmap of the conceptual, methodological, practical, and numerical challenges that remain in obtaining a universally applicable and truly predictive vdW method for realistic molecular systems and materials.
Collapse
Affiliation(s)
- Jan Hermann
- Fritz-Haber-Institut der Max-Planck-Gesellschaft , Faradayweg 4-6, 14195 Berlin, Germany
| | - Robert A DiStasio
- Department of Chemistry and Chemical Biology, Cornell University , Ithaca, New York 14853, United States
| | - Alexandre Tkatchenko
- Fritz-Haber-Institut der Max-Planck-Gesellschaft , Faradayweg 4-6, 14195 Berlin, Germany.,Physics and Materials Science Research Unit, University of Luxembourg , L-1511 Luxembourg, Luxembourg
| |
Collapse
|
25
|
De S, Musil F, Ingram T, Baldauf C, Ceriotti M. Mapping and classifying molecules from a high-throughput structural database. J Cheminform 2017; 9:6. [PMID: 28203290 PMCID: PMC5289135 DOI: 10.1186/s13321-017-0192-4] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2016] [Accepted: 01/17/2017] [Indexed: 11/10/2022] Open
Abstract
High-throughput computational materials design promises to greatly accelerate the process of discovering new materials and compounds, and of optimizing their properties. The large databases of structures and properties that result from computational searches, as well as the agglomeration of data of heterogeneous provenance leads to considerable challenges when it comes to navigating the database, representing its structure at a glance, understanding structure-property relations, eliminating duplicates and identifying inconsistencies. Here we present a case study, based on a data set of conformers of amino acids and dipeptides, of how machine-learning techniques can help addressing these issues. We will exploit a recently-developed strategy to define a metric between structures, and use it as the basis of both clustering and dimensionality reduction techniques-showing how these can help reveal structure-property relations, identify outliers and inconsistent structures, and rationalise how perturbations (e.g. binding of ions to the molecule) affect the stability of different conformers.
Collapse
Affiliation(s)
- Sandip De
- National Center for Computational Design and Discovery of Novel Materials (MARVEL), Lausanne, Switzerland.,Laboratory of Computational Science and Modelling, Institute of Materials, Ecole Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
| | - Felix Musil
- National Center for Computational Design and Discovery of Novel Materials (MARVEL), Lausanne, Switzerland.,Laboratory of Computational Science and Modelling, Institute of Materials, Ecole Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
| | - Teresa Ingram
- Theory Department of the Fritz Haber Institute, Faradayweg 4-6, 14195 Berlin-Dahlem, Germany
| | - Carsten Baldauf
- Theory Department of the Fritz Haber Institute, Faradayweg 4-6, 14195 Berlin-Dahlem, Germany
| | - Michele Ceriotti
- National Center for Computational Design and Discovery of Novel Materials (MARVEL), Lausanne, Switzerland.,Laboratory of Computational Science and Modelling, Institute of Materials, Ecole Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
| |
Collapse
|
26
|
Marianski M, Supady A, Ingram T, Schneider M, Baldauf C. Assessing the Accuracy of Across-the-Scale Methods for Predicting Carbohydrate Conformational Energies for the Examples of Glucose and α-Maltose. J Chem Theory Comput 2016; 12:6157-6168. [DOI: 10.1021/acs.jctc.6b00876] [Citation(s) in RCA: 73] [Impact Index Per Article: 9.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Mateusz Marianski
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany
| | - Adriana Supady
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany
| | - Teresa Ingram
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany
| | - Markus Schneider
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany
| | - Carsten Baldauf
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany
| |
Collapse
|
27
|
Ratcliff LE, Mohr S, Huhs G, Deutsch T, Masella M, Genovese L. Challenges in large scale quantum mechanical calculations. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2016. [DOI: 10.1002/wcms.1290] [Citation(s) in RCA: 70] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Affiliation(s)
- Laura E. Ratcliff
- Argonne Leadership Computing Facility Argonne National Laboratory Lemon IL USA
| | - Stephan Mohr
- Department of Computer Applications in Science and Engineering Barcelona Supercomputing Center (BSC‐CNS) Barcelona Spain
| | - Georg Huhs
- Department of Computer Applications in Science and Engineering Barcelona Supercomputing Center (BSC‐CNS) Barcelona Spain
| | - Thierry Deutsch
- University Grenoble Alpes INAC‐MEM Grenoble France
- CEA, INAC‐MEM Grenoble France
| | - Michel Masella
- Laboratoire de Biologie Structurale et Radiologie, Service de Bioénergétique, Biologie Structurale et Mécanisme Institut de Biologie et de Technologie de Saclay, CEA Saclay Gif‐sur‐Yvette Cedex France
| | - Luigi Genovese
- University Grenoble Alpes INAC‐MEM Grenoble France
- CEA, INAC‐MEM Grenoble France
| |
Collapse
|
28
|
Ropo M, Blum V, Baldauf C. Trends for isolated amino acids and dipeptides: Conformation, divalent ion binding, and remarkable similarity of binding to calcium and lead. Sci Rep 2016; 6:35772. [PMID: 27808109 PMCID: PMC5093913 DOI: 10.1038/srep35772] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2016] [Accepted: 10/03/2016] [Indexed: 12/16/2022] Open
Abstract
We derive structural and binding energy trends for twenty amino acids, their dipeptides, and their interactions with the divalent cations Ca2+, Ba2+, Sr2+, Cd2+, Pb2+, and Hg2+. The underlying data set consists of more than 45,000 first-principles predicted conformers with relative energies up to ~4 eV (~400 kJ/mol). We show that only very few distinct backbone structures of isolated amino acids and their dipeptides emerge as lowest-energy conformers. The isolated amino acids predominantly adopt structures that involve an acidic proton shared between the carboxy and amino function. Dipeptides adopt one of two intramolecular-hydrogen bonded conformations C5 or . Upon complexation with a divalent cation, the accessible conformational space shrinks and intramolecular hydrogen bonding is prevented due to strong electrostatic interaction of backbone and side chain functional groups with cations. Clear correlations emerge from the binding energies of the six divalent ions with amino acids and dipeptides. Cd2+ and Hg2+ show the largest binding energies-a potential correlation with their known high acute toxicities. Ca2+ and Pb2+ reveal almost identical binding energies across the entire series of amino acids and dipeptides. This observation validates past indications that ion-mimicry of calcium and lead should play an important role in a toxicological context.
Collapse
Affiliation(s)
- M. Ropo
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany
- Department of Physics, Tampere University of Technology, Finland
- COMP, Department of Applied Physics, Aalto University, Finland
| | - V. Blum
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany
- Department of Mechanical Engineering and Materials Science, Duke University, Durham, NC, USA
| | - C. Baldauf
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, D-14195 Berlin, Germany
| |
Collapse
|
29
|
Schwing K, Gerhards M. Investigations on isolated peptides by combined IR/UV spectroscopy in a molecular beam – structure, aggregation, solvation and molecular recognition. INT REV PHYS CHEM 2016. [DOI: 10.1080/0144235x.2016.1229331] [Citation(s) in RCA: 39] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
|
30
|
Otaki H, Yagi K, Ishiuchi SI, Fujii M, Sugita Y. Anharmonic Vibrational Analyses of Pentapeptide Conformations Explored with Enhanced Sampling Simulations. J Phys Chem B 2016; 120:10199-10213. [DOI: 10.1021/acs.jpcb.6b06672] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
| | | | - Shun-ichi Ishiuchi
- Laboratory
for Chemistry and Life Science, Institute for Innovative Research, Tokyo Institute of Technology, 4259 Nagatsuta-cho, Midori-ku, Yokohama 226-8503, Japan
| | - Masaaki Fujii
- Laboratory
for Chemistry and Life Science, Institute for Innovative Research, Tokyo Institute of Technology, 4259 Nagatsuta-cho, Midori-ku, Yokohama 226-8503, Japan
| | - Yuji Sugita
- RIKEN Advanced Institute for Computational Science, 7-1-26 Minatojima-Minamimachi,
Chuo-ku, Kobe, Hyogo 650-0047, Japan
- RIKEN Quantitative Biology Center, 1-6-5 Minatojima-Minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| |
Collapse
|
31
|
Ireta J. Polyalanine α-helix microsolvation: assessing the energy of the peptide desolvation penalty with density functional theory. Theor Chem Acc 2016. [DOI: 10.1007/s00214-016-1981-y] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
32
|
Voronina L, Masson A, Kamrath M, Schubert F, Clemmer D, Baldauf C, Rizzo T. Conformations of Prolyl–Peptide Bonds in the Bradykinin 1–5 Fragment in Solution and in the Gas Phase. J Am Chem Soc 2016; 138:9224-33. [DOI: 10.1021/jacs.6b04550] [Citation(s) in RCA: 50] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Affiliation(s)
- Liudmila Voronina
- Laboratoire
de Chimie Physique Moléculaire, École Polytechnique Fédérale de Lausanne, EPFL SB ISIC LCPM, Station 6, CH-1015 Lausanne, Switzerland
| | - Antoine Masson
- Laboratoire
de Chimie Physique Moléculaire, École Polytechnique Fédérale de Lausanne, EPFL SB ISIC LCPM, Station 6, CH-1015 Lausanne, Switzerland
| | - Michael Kamrath
- Laboratoire
de Chimie Physique Moléculaire, École Polytechnique Fédérale de Lausanne, EPFL SB ISIC LCPM, Station 6, CH-1015 Lausanne, Switzerland
| | - Franziska Schubert
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, D-14195 Berlin, Germany
| | - David Clemmer
- Department
of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Carsten Baldauf
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, D-14195 Berlin, Germany
| | - Thomas Rizzo
- Laboratoire
de Chimie Physique Moléculaire, École Polytechnique Fédérale de Lausanne, EPFL SB ISIC LCPM, Station 6, CH-1015 Lausanne, Switzerland
| |
Collapse
|
33
|
Theis T, Ortiz GX, Logan AWJ, Claytor KE, Feng Y, Huhn WP, Blum V, Malcolmson SJ, Chekmenev EY, Wang Q, Warren WS. Direct and cost-efficient hyperpolarization of long-lived nuclear spin states on universal (15)N2-diazirine molecular tags. SCIENCE ADVANCES 2016; 2:e1501438. [PMID: 27051867 PMCID: PMC4820385 DOI: 10.1126/sciadv.1501438] [Citation(s) in RCA: 181] [Impact Index Per Article: 22.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/12/2015] [Accepted: 02/05/2016] [Indexed: 05/17/2023]
Abstract
Conventional magnetic resonance (MR) faces serious sensitivity limitations which can be overcome by hyperpolarization methods, but the most common method (dynamic nuclear polarization) is complex and expensive, and applications are limited by short spin lifetimes (typically seconds) of biologically relevant molecules. We use a recently developed method, SABRE-SHEATH, to directly hyperpolarize (15)N2 magnetization and long-lived (15)N2 singlet spin order, with signal decay time constants of 5.8 and 23 minutes, respectively. We find >10,000-fold enhancements generating detectable nuclear MR signals that last for over an hour. (15)N2-diazirines represent a class of particularly promising and versatile molecular tags, and can be incorporated into a wide range of biomolecules without significantly altering molecular function.
Collapse
Affiliation(s)
- Thomas Theis
- Department of Chemistry, Duke University, Durham, NC 27708, USA
- Corresponding author. E-mail: (W.S.W.); (Q.W.); (T.T.)
| | | | | | | | - Yesu Feng
- Department of Chemistry, Duke University, Durham, NC 27708, USA
| | - William P. Huhn
- Department of Mechanical Engineering and Materials Science, Duke University, Durham, NC 27708, USA
| | - Volker Blum
- Department of Mechanical Engineering and Materials Science, Duke University, Durham, NC 27708, USA
| | | | - Eduard Y. Chekmenev
- Departments of Radiology and Biomedical Engineering, Vanderbilt University, Institute of Imaging Science, Nashville, TN 37232, USA
| | - Qiu Wang
- Department of Chemistry, Duke University, Durham, NC 27708, USA
- Corresponding author. E-mail: (W.S.W.); (Q.W.); (T.T.)
| | - Warren S. Warren
- Department of Chemistry, Duke University, Durham, NC 27708, USA
- Department of Physics, Duke University, Durham, NC 27708, USA
- Departments of Radiology and Biomedical Engineering, Duke University, Durham, NC 27708, USA
- Corresponding author. E-mail: (W.S.W.); (Q.W.); (T.T.)
| |
Collapse
|
34
|
First-principles data set of 45,892 isolated and cation-coordinated conformers of 20 proteinogenic amino acids. Sci Data 2016; 3:160009. [PMID: 26881946 PMCID: PMC4755128 DOI: 10.1038/sdata.2016.9] [Citation(s) in RCA: 41] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2015] [Accepted: 01/15/2016] [Indexed: 11/24/2022] Open
Abstract
We present a structural data set of the 20 proteinogenic amino acids and their amino-methylated and acetylated (capped) dipeptides. Different protonation states of the backbone (uncharged and zwitterionic) were considered for the amino acids as well as varied side chain protonation states. Furthermore, we studied amino acids and dipeptides in complex with divalent cations (Ca2+, Ba2+, Sr2+, Cd2+, Pb2+, and Hg2+). The database covers the conformational hierarchies of 280 systems in a wide relative energy range of up to 4 eV (390 kJ/mol), summing up to a total of 45,892 stationary points on the respective potential-energy surfaces. All systems were calculated on equal first-principles footing, applying density-functional theory in the generalized gradient approximation corrected for long-range van der Waals interactions. We show good agreement to available experimental data for gas-phase ion affinities. Our curated data can be utilized, for example, for a wide comparison across chemical space of the building blocks of life, for the parametrization of protein force fields, and for the calculation of reference spectra for biophysical applications.
Collapse
|
35
|
Goyal B, Srivastava KR, Kumar A, Patwari GN, Durani S. Probing the role of electrostatics of polypeptide main-chain in protein folding by perturbing N-terminal residue stereochemistry: DFT study with oligoalanine models. RSC Adv 2016. [DOI: 10.1039/c6ra22870d] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
Energetics of folding (ΔHE→F, in kcal mol−1) from the extended (E) structure to the folded (F) structure for Ia and Ib critically depend on the geometrical relationship between the backbone peptide units of the polypeptide structure.
Collapse
Affiliation(s)
- Bhupesh Goyal
- Department of Chemistry
- Indian Institute of Technology Bombay
- Mumbai-400076
- India
| | | | - Anil Kumar
- Department of Chemistry
- Indian Institute of Technology Bombay
- Mumbai-400076
- India
| | - G. Naresh Patwari
- Department of Chemistry
- Indian Institute of Technology Bombay
- Mumbai-400076
- India
| | - Susheel Durani
- Department of Chemistry
- Indian Institute of Technology Bombay
- Mumbai-400076
- India
| |
Collapse
|
36
|
Baldauf C, Rossi M. Going clean: structure and dynamics of peptides in the gas phase and paths to solvation. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2015; 27:493002. [PMID: 26598600 DOI: 10.1088/0953-8984/27/49/493002] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
The gas phase is an artificial environment for biomolecules that has gained much attention both experimentally and theoretically due to its unique characteristic of providing a clean room environment for the comparison between theory and experiment. In this review we give an overview mainly on first-principles simulations of isolated peptides and the initial steps of their interactions with ions and solvent molecules: a bottom up approach to the complexity of biological environments. We focus on the accuracy of different methods to explore the conformational space, the connections between theory and experiment regarding collision cross section evaluations and (anharmonic) vibrational spectra, and the challenges faced in this field.
Collapse
Affiliation(s)
- Carsten Baldauf
- Fritz Haber Institute, Faradayweg 4-6, 14195 Berlin, Germany
| | | |
Collapse
|
37
|
Supady A, Blum V, Baldauf C. First-Principles Molecular Structure Search with a Genetic Algorithm. J Chem Inf Model 2015; 55:2338-48. [DOI: 10.1021/acs.jcim.5b00243] [Citation(s) in RCA: 68] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Adriana Supady
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, 14195 Berlin, Germany
| | - Volker Blum
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, 14195 Berlin, Germany
- Department of Mechanical Engineering & Materials Science, Duke University, Durham, North Carolina 27708, United States
| | - Carsten Baldauf
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, 14195 Berlin, Germany
| |
Collapse
|