Header

UZH-Logo

Maintenance Infos

Potential energy surface-based automatic deduction of conformational transition networks and its application on quantum mechanical landscapes of d-glucose conformers


Satoh, Hiroko; Oda, Tomohiro; Nakakoji, Kumiyo; Uno, Takeaki; Tanaka, Hiroaki; Iwata, Satoru; Ohno, Koichi (2016). Potential energy surface-based automatic deduction of conformational transition networks and its application on quantum mechanical landscapes of d-glucose conformers. Journal of Chemical Theory and Computation, 12(11):5293-5308.

Abstract

This paper describes our approach that is built upon the potential energy surface (PES)-based conformational analysis. This approach automatically deduces a conformational transition network, called a conformational reaction route map (r-map), by using the Scaled Hypersphere Search of the Anharmonic Downward Distortion Following method (SHS-ADDF). The PES-based conformational search has been achieved by using large ADDF, which makes it possible to trace only low transition state (TS) barriers while restraining bond lengths and structures with high free energy. It automatically performs sampling the minima and TS structures by simply taking into account the mathematical feature of PES without requiring any a priori specification of variable internal coordinates. An obtained r-map is composed of equilibrium (EQ) conformers connected by reaction routes via TS conformers, where all of the reaction routes are already confirmed during the process of the deduction using the intrinsic reaction coordinate (IRC) method. The postcalculation analysis of the deduced r-map is interactively carried out using the RMapViewer software we have developed. This paper presents computational details of the PES-based conformational analysis and its application to d-glucose. The calculations have been performed for an isolated glucose molecule in the gas phase at the RHF/6-31G level. The obtained conformational r-map for α-d-glucose is composed of 201 EQ and 435 TS conformers and that for β-d-glucose is composed of 202 EQ and 371 TS conformers. For the postcalculation analysis of the conformational r-maps by using the RMapViewer software program we have found multiple minimum energy paths (MEPs) between global minima of 1C4 and 4C1 chair conformations. The analysis using RMapViewer allows us to confirm the thermodynamic and kinetic predominance of 4C1 conformations; that is, the potential energy of the global minimum of 4C1 is lower than that of 1C4 (thermodynamic predominance) and that the highest energy of those of all the TS structures along a route from 4C1 to 1C4 is lower than that of 1C4 to 4C1 (kinetic predominance).

Abstract

This paper describes our approach that is built upon the potential energy surface (PES)-based conformational analysis. This approach automatically deduces a conformational transition network, called a conformational reaction route map (r-map), by using the Scaled Hypersphere Search of the Anharmonic Downward Distortion Following method (SHS-ADDF). The PES-based conformational search has been achieved by using large ADDF, which makes it possible to trace only low transition state (TS) barriers while restraining bond lengths and structures with high free energy. It automatically performs sampling the minima and TS structures by simply taking into account the mathematical feature of PES without requiring any a priori specification of variable internal coordinates. An obtained r-map is composed of equilibrium (EQ) conformers connected by reaction routes via TS conformers, where all of the reaction routes are already confirmed during the process of the deduction using the intrinsic reaction coordinate (IRC) method. The postcalculation analysis of the deduced r-map is interactively carried out using the RMapViewer software we have developed. This paper presents computational details of the PES-based conformational analysis and its application to d-glucose. The calculations have been performed for an isolated glucose molecule in the gas phase at the RHF/6-31G level. The obtained conformational r-map for α-d-glucose is composed of 201 EQ and 435 TS conformers and that for β-d-glucose is composed of 202 EQ and 371 TS conformers. For the postcalculation analysis of the conformational r-maps by using the RMapViewer software program we have found multiple minimum energy paths (MEPs) between global minima of 1C4 and 4C1 chair conformations. The analysis using RMapViewer allows us to confirm the thermodynamic and kinetic predominance of 4C1 conformations; that is, the potential energy of the global minimum of 4C1 is lower than that of 1C4 (thermodynamic predominance) and that the highest energy of those of all the TS structures along a route from 4C1 to 1C4 is lower than that of 1C4 to 4C1 (kinetic predominance).

Statistics

Citations

Dimensions.ai Metrics
7 citations in Web of Science®
7 citations in Scopus®
7 citations in Microsoft Academic
Google Scholar™

Altmetrics

Downloads

1 download since deposited on 16 Dec 2016
0 downloads since 12 months
Detailed statistics

Additional indexing

Item Type:Journal Article, refereed, original work
Communities & Collections:07 Faculty of Science > Department of Chemistry
Dewey Decimal Classification:540 Chemistry
Language:English
Date:2016
Deposited On:16 Dec 2016 10:06
Last Modified:02 Feb 2018 10:56
Publisher:American Chemical Society (ACS)
ISSN:1549-9618
OA Status:Closed
Publisher DOI:https://doi.org/10.1021/acs.jctc.6b00439
PubMed ID:27673598

Download