УДК: 539.27+544.18
СТРОЕНИЕ МОЛЕКУЛЫ ТРИПТОФАНА: ЭЛЕКТРОНОГРАФИЧЕСКОЕ И КВАНТОВО-ХИМИЧЕСКОЕ ИССЛЕДОВАНИЕ
В.В. Дунаева, Г.В. Гиричев, Н.И. Гиричева
Валерия Валерьевна Дунаева*, Георгий Васильевич Гиричев
Кафедра физики, Ивановский государственный химико-технологический университет, Шереметевский
просп., 7, Иваново, Россия, 153000
E-mail: tyunina.vv@gmail.com*, g.v.girichev@mail.ru
Нина Ивановна Гиричева
Кафедра органической и физической химии, Ивановский государственный университет, ул. Ермака, 39, Иваново, Россия, 153025 E-mail: n.i.giricheva@mail.ru
В представленной работе теоретическими (квантово-химическиерасчеты) и экспериментальными (газовая электронография) методами проведено исследование структуры молекулы и конформационных свойств триптофана. Расчеты выполнены в программе Gaussian 03 методом B3LYP с базисным набором cc-pVTZ. Они показали наличие шести кон-формеров, которые существуют в паре при температуре эксперимента (T = 495 K). Теоретически оценена чувствительность данных электронографии к конформационным изменениям в структуре молекулы триптофана. Конформеры имеют различную ориентацию как карбоксильной и аминогруппы, так и основной цепи и индольного фрагмента относительно друг друга. На основании этого их можно разделить на две группы: различимые (с разным торсионным углом C(OOH)-C(HNH2)-C(H2)-C(ind)) и слаборазличимые (с разными торсионными углами H-N-C-C и H-O-C-C) в рамках электронографического эксперимента. Определены молекулярные параметры конформеров триптофана. Установлено, что в конформерах этой аминокислоты присутствует внутримолекулярная водородная связь H2N"HO. Анализ электронографических данных проведен на основе предположения, что насыщенный пар над триптофаном при температуре эксперимента T = 495 K представляет собой смесь как минимум двух конформеров с наименьшими энергиями. Показано, что наилучшее соотношение этих конформеров между собой I : IIявляется 50: 50, соответственно. Проведено сравнение геометрических параметров молекул аминокислот (глицин, аланин, триптофан), полученных из данных газовой электронографии в этой работе и ранее. Установлено влияние внутримолекулярной водородной связи на структурные параметры основной цепи молекул указанных аминокислот.
Ключевые слова: триптофан, аминокислота, газовая электронография, квантово-химические расчеты, конформер, структура молекулы
MOLECULAR STRUCTURE OF TRYPTOPHAN: GAS-PHASE ELECTRON DIFFRACTION
AND QUANTUM-CHEMICAL STUDIES
V.V. Dunaeva, G.V. Girichev, N.I. Giricheva
Valeriya V. Dunaeva*, Georgy V. Girichev
Department of Physics, Ivanovo State University of Chemistry and Technology, Sheremetevskiy ave., 7, Ivanovo, 153000, Russia
E-mail: tyunina.vv@gmail.com*, g.v.girichev@mail.ru Nina I. Giricheva
Department of Organic and Physical Chemistry, Ivanovo State University, Ermak st., 39, Ivanovo, 153025, Russia E-mail: n.i.giricheva@mail.ru
The molecular structure and conformational properties of tryptophan have been investigated by gas-phase electron diffraction and theoretical methods. Quantum chemical calculations realized by program Gaussian 03 (B3LYP/cc-pVTZ) have been predicted the existence of six con-formers at the temperature of experiment (T = 495 K). The ability of gas-phase electron diffraction
method to distinguish the structure of conformers was estimated theoretically. Conformers have different orientations of carboxylic and amine group, backbone and indole fragment to each other. These conformers can be divided on two groups: distinguishable parameters (with different torsion angle C(OOH)-C(HNH2)-C(H2)-C(ind)) and weekly distinguishable ones (with different torsion angles H-N-C-C and H-O-C-C) by gas-phase electron diffraction. The molecular parameters of the conformers were determined. The conformers have intramolecular hydrogen bonding of the HN"HO. The analysis of the gas-phase electron diffraction data have been carried out assuming the saturated vapor of tryptophan at T = 495 K consists of mixture at least of two conformers with lowest energy. It was shown that optimal ratio between conformers I : II was 50 : 50, respectively. The geometrical parameters of amino acids molecules (glycine, alanine, tryptophan) obtained by gas-phase electron diffraction were compared. The influence of the intramolecular hydrogen bond was established onto the structural parameters of the backbone of the above amino acids molecules.
Key words: tryptophan, amino acid, gas-phase electron diffraction, quantum chemical calculations, conformer, molecular structure
Для цитирования:
Дунаева В.В., Гиричев Г.В., Гиричева Н.И. Строение молекулы триптофана: электронографическое и квантово-хи-мическое исследование. Изв. вузов. Химия и хим. технология. 2020. Т. 63. Вып. 3. С. 37-45 For citation:
Dunaeva V.V., Girichev G.V., Giricheva N.I. Molecular structure of tryptophan: gas-phase electron diffraction and quantum-chemical studies. Izv. Vyssh. Uchebn. Zaved. Khim. Khim. Tekhnol. [Russ. J. Chem. & Chem. Tech.]. 2020. V. 63. N 3. P. 37-45
INTRODUCTION
Amino acids are the building blocks of peptides and proteins having fundamental chemical, biochemical and biological significance. They take part in many biochemical processes (neuroregulation, enzyme catalysis, synthesis of bioactive substances) [1]. In the last years, the synthetic peptides with amino acid residues have gained high interest as new drugs to overcome bacterial resistance to classical antibiotics [2]. Peptides have been used in bionanotechnology (biosensor, molecular electronics) [3]. Amino acid molecules have a unique structure. Zwitterionic forms of amino acids are stabilized in the crystalline state and in solution whereas the neutral structures are found in the gas phase [4]. Therefore, the structure of amino acids are still the objects of careful study in both of these media.
In the gaseous phase the aromatic amino acids, tryptophan (Trp), phenylalanine (Phe), tyrosine (Tyr), have been studied by UV spectroscopy due to presence of UV chromophores in their side chains [5], matrixisolation spectroscopy with quantum chemical calculations [6,7], UV and IR ion-dip spectroscopy combined with quantum chemical calculations [8], microwave spectroscopy [9] and other methods.
Amino acid molecules have low symmetry and many rotational degrees of freedom. The conformational behavior of amino acids in the gas phase is considerably complicated, particularly for the aromatic amino acids with high rotational freedom or plausible intramolecular hydrogen bonds (H-bonds). For this reason a large number of different low-energy conformers are expected on the surface of potential energy. Therefore,
it is difficult to carry out a structural study of the amino acid molecule in the gas phase by experimental methods. Nevertheless, the molecular structures of some amino acids and bioactive compounds have been experimentally investigated by the gas-phase electron diffraction (GED) [4, 10-14]. The structures and preferred conformations of aromatic amino acids in gas phase are still indefinite. The molecular structures of the most stable conformers of neutral glycine (Gly), alanine (Ala), and proline (Pro) have been determined by the joint analysis of electron diffraction data and rotational constants, in conjunction with the results of an ab initio calculation [4, 13, 14]. Study of the Trp molecule was carried out by ultrafast electron diffraction in work [15]. The structure of Trp conformer without intramolecular hydrogen bond was obtained by experiment and quantum chemical calculations. However, the conformer obtained was no minimum on the potential energy surface in the quantum chemical studies of gaseous molecule of Trp.
The aim of our study was to determine the geometrical parameters and conformational composition of gaseous Trp using electron diffraction and quantum chemistry methods. The theoretical radial distribution functions are compared for different types of conformers in electron diffraction. It should be noted that the presented study of the Trp led to the identification of highly accurate molecular structures, which demonstrated the reliability of the GED method associated with high-level quantum-chemical calculations for determining the structure of the amino acid. The main goal of this work is experimental and theoretical investigations of the structure
and properties of the aromatic amino acid, Trp, and their conformational dependency by synchronous GED and mass spectrometry (GED/MS) along with high level quantum chemical calculations. We propose conformational classification for the neutral form of Trp with a density functional theory using the Gaussian 03 program system [16]. We investigate the differences in structure and charge distribution between the subgroups of Trp.
EXPERIMENTAL PART
Chemicals
Commercially available L-tryptophan [CASRN 73-22-3] with purity of 99%, was obtained from Al-drich Chemical Co., and was used without further purification. Before use, the sample was carefully dried under vacuum during 48 h at the temperature 323 K and then stored over P2O5 in desiccators, no further purification being followed.
Computational details
In this work, quantum chemical calculations of geometry and electronic properties of Trp's conform-ers were carried out by Gaussian 03 [16] at the B3LYP/ cc-pVTZ level of theory. The charge distribution was analyzed using the natural bond orbital (NBO) analysis
[17] using the Gaussian 03 package [16].
GED/MS experiment
Electron diffraction patterns were recorded simultaneously with the mass spectra using special techniques. The electron diffraction experiment was carried out using the combined GED/MS apparatus
[18], consisting of the EMR-100 gas-phase electron diffraction and the APDM-1 monopole mass spectro-metric units. The conditions of the GED/MS experiment are shown in Table 1.
The temperature of stainless steel effusion cell was measured by a W/Re-5/20 thermocouple, which was standardized using the melting points of Sn and Al. Polycrystalline ZnO was used to establish the wavelength of the fast electrons recorded before and after the GED experiments on Trp.
Table 1
Experimental conditions for GED/MS study of Trp Таблица 1. Экспериментальные условия синхронного электронографического и масс-спектрометри-
The Trp sample was evaporated at T = 495±5 K from a stainless steel effusion cell. Microphotometric measurements were carried out by means of automatic techniques [19]. The scattered electrons were collected on the Kodak Electron Image films SO-163 of 9x 12 cm at two camera distances, 598 and 338 mm. The electron diffraction patterns were scanned by a modified [19] computer-controlled microdensitometer Carl Zeiss Jena MD-100. For each film, a rectangular area of 10 x 130 mm was scanned at a diagonal direction (33 equidistant lines with a step of 0.1 mm along each line). The molecular intensities were obtained in the ranges 2.628.5 Á-1 (short camera) and 1.2-15.0 Á-1 (long camera).
The mass spectra of the effusing molecular beam were recorded simultaneously with recording diffraction patterns. These data indicate that the monomeric molecules of Trp are presented in the vapor. No volatile impurities were detected at a level of 1%, at least. Mass-spectrum is the same as a described in [20].
RESULTS AND DISCUSSION
Quantum chemical calculations
Several theoretical studies have been carried out on the conformational analyses of neutral trypto-phan [21, 22]. In case of structural analyses of electron diffraction experiment, the data of geometry of molecule, force field, vibrational amplitudes, and vibra-tional corrections of intermolecular distances were necessary. We perform the conformational analysis of the Trp molecule at the B3LYP/cc-pVTZ level of electronic structure theory using the Gaussian 03 program [16]. The calculations predict 40 stable conformations of neutral tryptophan, energies of which lie within 10 kJ/mol. No imaginary frequencies were observed for any of the conformers examined, confirming that the optimized structures were minima on the potential energy surface. All conformers can be divided into 6 groups in relation to different orientations of car-boxylic and amine group, backbone and indole fragment to each other. Relative energies and structures of conformer groups are given in Supplementary (Table S1 and Fig. S1 http://journals.isuct.ru/ctj/arti-cle/view/1985/1128). Structures of six most stable conformers with the lowest-energy are shown in Fig. 1. Each of them has an intramolecular hydrogen bond.
Because intramolecular hydrogen bonds play an important role in the molecular structures and con-formational preferences of amino acids, for the six low-energetic conformers of neutral Trp, NBO-analy-sis at the B3LYP/cc-pVTZ level of theory was performed. The most stable Trp conformers exhibit an intramolecular hydrogen bond between the hydrogen atom of the hydroxyl group and the lone pair of electrons on the nitrogen atom of the amino group (N- • H-O) (Fig. 1).
ческого исследования Trp
Nozzle-to-plate distance, mm 338 598
Electron beam current, ^A 1.50 0.70
Ionizing voltage, V 50 50
Accelerating voltage, kV 94.0 90.8
Temperature of effusion cell, °C 222 222
Residual gas pressure diffraction chamber/MS-block, torr 2 10-6/ 6 10-7 3 • 10-6/ 5 10-7
Scattering angles 5mm - .w, Â-1 3.7-31.0 1.2-17.2
I
Erei = 0 kJ/mol
II
Erei = 1.80 kJ/mol
III
/•;,,i = 4.02 kJ/mol
IV
Erei = 4.81 kJ/mol
V
Erei = 4.90 kJ/mol Fig. 1. Low-energetic conformers of Trp (Erei - relative energy) Рис. 1. Конформеры Trp с наименьшей энергией (Erei - относительная энергия)
VI
Erei = 6.19 kJ/mol
As in some other amino acids [23], the hydrogen bond between the carboxyl hydrogen and the nitrogen in the amino group is the strongest. The average energies of the interaction, £'(LP(N) ^ o*(O-H)), between atomic orbital of N with a lone electron pair (LP(N)) as a donor and an "antibonding" orbital of the O-H bond (o*(O-H)) as acceptor, and transfer charge values (qr) are shown in table 2. The amount of transfer charge (qtr) is higher than 0.01е [24].
Table 2
Energies of orbital interactions LP(N)-g*(0-H) and amounts of charge transfer (qtr) for main conformers of Trp based on NBO theory (at B3LYP/cc-PVTZ level of electronic structure theory) Таблица 2. Энергии взаимодействия орбиталей LP(N) с с*(О-Н) и величины перенесенного заряда
Conformer £(LP(N)^-ct*(0-H)), kJ/mol qtr
I 63.89 0.033
II 57.20 0.029
III 61.92 0.031
IV 56.57 0.029
V 60.25 0.031
VI 66.36 0.035
As it is known, the relative energy values of these conformers reflect their concentrations in the gas
phase at Texp. The six of the 40 Trp conformers obtained have a significant concentration in a saturated vapor (72% at T = 495 K). The concentrations of the other conformers are negligible. These minor species could not be detected in the GED study. Thus, the main contributions into the electron scattering intensities are made by the six conformers during the electron diffraction experiment.
Sensitivity of GED method to distinguish the structure of conformers
Among 6 selected conformers of Trp the conformers II and III are indistinguishable for GED because of them differs only by hydrogen orientation in amino group. The same reason explains the closeness of GED data of conformers IV and VI. Conformers I and V have the same structure of backbone and differ by indole ring orientation. Radial distribution functions f(r) for the Trp conformers with different orientation of indole ring and backbone were compared (Fig. 2). Electron diffraction data were found to be more sensitive to the rotation of the backbone than the indole ring.
Tryptophan conformers may be distinguishable, since the discrepancy factor Rf between theoretical functions f(r) for the conformers exceeds the level determined by the experimental noise of photo plate. Thus, the pair of conformers I and V, II and IV, III and
tr
VI may be integrated. Taking into account the indistin-guishing by GED conformers II and III, IV and VI, the six conformers can be reduced to I and II.
f(r)
f)
0 ' 4 ' 8 12
a
f(r)
fr)
I-1-1-1-1-1-1
0 4 8 A 12
b
Fig. 2. Radial distribution functions of Trp. f(r) - theoretical radial distributions, Af(r) - difference curve a) for I and V conform-ers: different orientation of indole ring Rf = 5.3%; b) for I and II
conformers: rotation of backbone Rf = 7.7% Рис. 2. Функции радиального распределения Trp. f(r) - теоретические радиальные распределения, Af(r) - разностная кривая a) для I и V конформеров: различная ориентация индоль-ного кольца Rf = 5,3%; b) для I и II конформеров: вращение основной цепи Rf = 7,7%
Structure analysis
Structural analysis was performed under the assumption that there are 2 conformers in the Trp vapour — I and II with different backbone locations relative to Ca-Cp. A conventional least squares analysis of sM(s) was carried out using a modified version of KCED program [25]. The starting values for bond distances and angles were taken from the B3LYP/cc-pVTZ calculations for I and II. Vibrational corrections, Ar = rhi - ra, and the starting values for vibrational amplitudes were derived from the B3LYP/cc-pVTZ force field, using the curvilinear approach of Sipachev (program SHRINK) [25]. According to Fig. 2, the indole ring has a weak influence on radial distribution functions. The two conformers, I and II, were only used for structural analysis.
The molecular structures of each conformer were constructed with 75 independent structural para-
sM(s)
s, A-1
Fig. 3. Experimental (dots) and theoretical (solid line) molecular intensities sM(s) curves at long (1) and short (2) nozzle-to-plate distances and the difference curves AsM(s) at long (3) and short
(4) nozzle-to-plate distances for Trp Рис. 3. Экспериментальная (точки) и теоретическая (сплошная линия) функции молекулярной составляющей интенсивности рассеяния sM(s) для длинного (1) и короткого (2) расстояния сопло-фотопластинка и разностные кривые AsM(s) для длинного (3) и короткого (4) расстояния сопло-фотопластинка для молекулы Trp
meters: 28 bond distances, 25 bond angles and 22 torsion angles. In our investigation an attempt to reduce the complexity of model was done by means of decreasing the number of parameters. The following assumptions were used: (i) the 4 bond distances (N1-H16, C2-C3, C10-C3, C12-O13) were chosen independently variable (the difference between the none-quivalent distances of the same type was hold as obtained in quantum chemical calculations); (ii) the 9 bond angles (N3-C2-C1, C4-N1-C2, C6-C5-C4, C10-C3-C2, C11-C10-C3, O13-C12-C11, O14-C12-C11, N15-C11-C10, N15-C11-C12) independently modified (the differences between all angles fixed were constrained to calculated quantum chemical values); (iii) the 3 torsion angles (C1-C10-C3-C2, C12-C11-C10-C3, O13-C12-C11-C10) were chosen independently variable. With these restrictions, the number of parameters reduces to 16. The independent rh parameters were used for the geometrical constructions. Vibra-tional amplitudes were refined in groups with fixed differences. The experimental molecular intensities and the radial distributing curves are shown in Fig. 3 and Fig. 4, respectively. The molecular intensities sM(s) were obtained in the ranges 3.7-31.0 A-1 and 1.2-17.2 A-1 for the short and long camera distance, respectively (Fig. 3). The molecular intensity was calculated by the formula
sM(s) = s[/(s) - G(s)]/G(s).
The results of the least squares analysis of experimental data showed that at least, two conformers of tryptophan are present as the most stable conformers
Table 3
Experimental (conformer I) and optimized geometric parameters (conformers I and II) of Trp Таблица 3. Экспериментальные (конформер I) и теоретические структурные параметры (конформеры I
0 4 8 r, А
Fig. 4. Experimental (dots) and theoretical (line) radial distribution curves f(r) for Trp and the difference curves Af(r) for various
ratio of Trp conformers Рис. 4. Экспериментальная (точки) и теоретическая (линия) кривые радиального распределения f(r) для Trp и разностная кривая hf(r) для различного соотношения конформеров Trp
in the gas phase in ratio as I: II = 50:50. Under the experimental condition of T = 495 K, about 50% of the gas exists in the I and V conformers form and 50% of tryptophan in the gas is probably a mixture of II, III, IV, and VI conformers. It should be noted that the molecular constants of these conformers have no significant differences (Table 3). Final results of the least squares analysis of GED data are given in Table 3. The best fit of the experimental scattering intensities (Rf = 4.02%) was obtained for the conformer ratio 50:50. The experimental structural parameters are found to be in good agreement with the results of theoretical calculations.
Vibrational amplitudes and correction terms, Ar = rhi-ra, were derived from the theoretical force field taking into account nonlinear kinematic effects at the level of the first order perturbation theory with the use of the program SHRINK [25]. Some of the most important amplitudes and correction terms are listed in Table 4.
Conformer I I II
Parameters: r, Â; Z, ° GED, rh1 B3LYP/ cc-pVTZ, re B3LYP/ cc-pVTZ, re
N1-H16 1.014(4) " pi 1.003 1.003
C2-H17 1.087(4)(p) 1.077 1.078
C10-H22 1.099(4)(p) 1.089 1.093
O14-H25 0.994(4)p) 0.984 0.982
C2-C3 1.374(3) p2 1.370 1.370
N1-C2 1.380(3)p) 1.376 1.378
C10-C3 1.498(5) p3 1.500 1.495
C12-O13 1.195(6) p4 1.203 1.202
N3-C2-C1 110.1 (10) b p5 110.0 110.1
C4-N1-C2 110.6(10) pe 109.4 109.2
C6-C5-C4 118.7(2) p7 118.6 118.8
C10-C3-C2 125.9(20) ps 126.4 126.6
C11-C10-C3 112.5(13) pg 113.9 113.6
O13-C12-C11 123.9(25) pw 123.3 122.9
O14-C12-C11 110.0 (25) pn 113.7 113.9
N15-C11-C10 118.2(50) p!2 115.9 110.7
N15-C11-C12 111.3 (18) p13 109.3 109.9
C1-C10-C3 -C2 82(6) b pi4 84.5 108.2
C12-C11-C10-C3 -71 (5) pi5 -69.9 176.4
O13-C12-C11-C10 -35 (12) pi6 -37.7 -73.1
mol.% 50 50 50
Rf, % 4.017 — —
Notes: a -uncertainty in rhi- distances were calculated as a=(aSc2+(2.5aLs)2)1/2 (asc =0.002r, ols -standard deviation in LS-analysis); b uncertainty for angles is 3ols Примечания: a -погрешность в rhi- расстояниях вычислялась как a=(asc2+(2,5aLs)2)1/2 (asc =0,002r, ols -стандартное отклонение в МНК-анализе); b погрешность в величинах углов принята равной 3ols
Table 4
Bond distances, experimental and calculated vibration amplitudes and vibrational corrections (without non-bonded
distances involving hydrogen atoms) a Таблица 4. Межъядерные расстояния, экспериментальные и теоретические амплитуды колебаний и колеба-
Interaction c r b ' a l ь lexp lcalc Ar = rh1 Г a Refinement group
1 2 3 4 5 6
O14- H25 0.9935 0.074(2) 0.072 0.0007 1
N 1- H16 1.0116 0.071 0.069 0.0020 1
C 2- H17 1.0854 0.076 0.074 0.0015 1
C 8- H20 1.0903 0.077 0.075 0.0014 1
C 6- H18 1.0915 0.077 0.075 0.0014 1
C10-H22 1.0984 0.078 0.076 0.0011 1
O14- C12 1.3282 0.046(1) 0.046 0.0005 2
C 3- C 2 1.3730 0.044 0.044 0.0012 2
C 4- N 1 1.3819 0.046 0.046 -0.0037 2
C 2- N 1 1.3843 0.046 0.046 -0.0039 2
C 6- C 7 1.3864 0.045 0.045 0.0009 2
C 9- C 8 1.3876 0.045 0.045 0.0006 2
a
1 2 3 4 5 6
C 4- C 9 1.3974 0.046 0.046 -0.0000 2
C 5- C 4 1.4212 0.047 0.047 0.0019 2
C 3- C 5 1.4390 0.049 0.049 0.0051 2
C11- N15 1.4743 0.052 0.052 -0.0001 2
C10- C 3 1.4973 0.051 0.051 0.0002 2
C11-C10 1.5394 0.055 0.056 0.0005 2
C11-C12 1.5411 0.056 0.056 -0.0004 2
O14...O13 2.2496 0.057(3) 0.055 0.0008 3
C 2...C 4 2.2620 0.054 0.053 0.0054 3
C 5...N 1 2.2283 0.055 0.053 -0.0038 3
C11...H25 2.1611 0.145 0.144 0.0052 3
C 3...C 4 2.3307 0.056 0.055 0.0067 3
C 5...C 7 2.4074 0.058 0.057 0.0032 3
O13...C11 2.4144 0.069 0.068 0.0052 3
C 4...C 6 2.4278 0.060 0.059 0.0062 3
C12...N15 2.4852 0.072 0.070 0.0044 3
O14...N15 2.5531 0.108 0.107 -0.0071 3
C 4...C 7 2.7602 0.065(4) 0.066 0.0072 4
C 6...C 9 2.8329 0.065 0.066 0.0067 4
O13...C10 2.8552 0.126 0.127 -0.0009 4
C 3...N15 3.0789 0.156 0.157 0.0125 4
C10...C 6 3.3218 0.114(7) 0.119 0.0136 5
O13...C 2 3.2611 0.352 0.357 0.0077 5
C 5...N15 3.5843 0.256(6) 0.263 0.0189 6
C 2...C 9 3.6329 0.056 0.063 0.0117 6
C 3...C 9 3.6569 0.055 0.062 0.0121 6
C10...N 1 3.7005 0.060 0.067 0.0051 6
C 8...N 1 3.7228 0.059 0.066 0.0076 6
C 6...N15 3.8193 0.415 0.423 0.0430 6
C10...C 4 3.7852 0.062 0.069 0.0153 6
O14...C 3 3.8775 0.265 0.272 0.0456 6
C11...C 6 4.0970 0.248(5) 0.249 0.0357 7
C 7...N 1 4.1090 0.068 0.069 0.0079 7
C12...N 1 4.4686 0.279(4) 0.275 0.0203 8
C 4...N15 4.6176 0.281 0.278 0.0285 8
C 2...C 8 4.5780 0.070 0.066 0.0166 8
O13...N 1 4.5343 0.388 0.385 0.0274 8
N15...N 1 4.7524 0.300 0.296 0.0368 8
O13...C 5 4.7228 0.273 0.270 0.0575 8
C12...C 4 5.0018 0.252(10) 0.259 0.0445 9
O14...N 1 4.8924 0.467 0.474 0.0704 9
C12...C 6 5.1567 0.292(14) 0.296 0.0719 10
O13...C 4 5.3020 0.338 0.342 0.0573 10
C 9...N15 5.6223 0.330 0.334 0.0332 10
O13...C 9 6.6539 0.371(32) 0.362 0.0775 11
O14...C 7 6.7211 0.551 0.542 0.1267 11
Notes: a -all values in A; b - total error limit in interatomic distances was estimated by formula g=(gsc2+(2.5gls)2)1/2 (Osc =0.002r, ols -standard deviation in LS-analysis); uncertainty for amplitudes are 3ols, c Numbering of atoms see at Fig. S2 http://journals.isuct.ru/ctj/article/view/1985/1128.
Примечания: a - все величины в A; b - погрешность в межатомных расстояниях оценивалась по формуле c=(csc2+ +(2,5cls)2)12 (Osc =0,002r, ols -стандартное отклонение в МНК-анализе); погрешность для амплитуд составила 3ols, с Нумерация атомов дана на рис. S2 http://journals.isuct.ru/ctj/article/view/1985/1128
Despite the great complexity of the structure of the amino acid molecule due to the need to take into account many of its conformers during the study, the
experimental data obtained are in good agreement with the theoretical calculations.
since amino acid molecules include the backbone and the various side chains, a certain structural features of the aliphatic and aromatic molecules might be expected. For three amino acids molecules, such as Gly and Ala with aliphatic side fragment and Trp with indole ring in side chain, the backbone deviates from the plane differently depending on the side chain: Gly is a flat structure, Ala is a strong distortion, Trp is the average between them.
The most stable forms of Trp molecule have the O=C-O-H dihedral angle in the trans orientation. As it was shown before [6, 8, 26-29], the stabilization of conformers of this type is characteristic of aromatic amino acids, contrary to the aliphatic family, where the conformers with the carboxylic group in the cis orientation and the bifurcated N-H.. ,O=C bond are of lower energy [4, 14, 30]. In other words, the side chains influence the intramolecular H-bond formation. For the aliphatic (Gly [4], Ala [14]) and the hydroxyl-contain-ing (ser, Thr) [31] amino acids the conformers with hydrogen-bond between amino and carboxylic groups do not have the lowest energies, whereas the presence of aromatic fragment (Tyr, Phe [6] and Trp) in molecule structure is promoted to stabilize the conformers with one H-bond.
The intramolecular hydrogen bonds affect onto the structural parameters of the backbone of the above amino acids molecules. As can be seen from the data Table 5, in case of Trp molecule, Z(Cm-C-N) increase and Z(C-C-N) decrease, but r(C-N) longer and r(C-O) shorter on the contrary to analogical parameters of Gly and Ala molecules.
Table 5
Some geometrical parameters of amino acids by GED data Таблица 5. Некоторые геометрические параметры
CONCLUSIONS
The first study on structure determination of tryptophan molecule was performed employing the gas electron diffraction method. Geometrical parameters and conformational composition of neutral tryptophan were determined based on joint analysis of GED data and quantum chemical force constants. The DFT
method with the B3LYP/cc-pVTZ level was used to optimize geometries and to obtain corrections to the total energy and Gibbs free energies to verify that all structures had real minima. A total of 40 minima have been located on the conformational space of Trp, within the 0-10 kJ/mol energy range. The conforma-tional distribution of tryptophan at the temperature of sublimation of the compound (495 K) was calculated, and six conformers were predicted to have abundances higher than 6%. The composition of saturated vapor of tryptophan mainly contains conformers I and II. The assignment of the experimental data was based on the calculated molecular structures of these 2 conformers.
The GED refined structure of tryptophan was found to be close to the predicted one from quantum chemical calculation. The well agreement of observed and calculated molecular intensity curves was shown.
An NBO analysis was used to get quantitative evaluations of interaction energies of atomic pairs involved in intramolecular hydrogen bond. The results obtained support the presence of intramolecular hydrogen bond in all six low-energetic conformers of neutral Trp. The observed conformers for Trp are stabilized by an O-H--N interaction between the carboxylic and amino groups. Comparison of GED data of the aliphatic amino acids (Gly, Ala) with the aromatic amino acid (Trp) allowed to indicate the influence of intramolecular hydrogen bond on the structural parameters of the backbone of the molecules.
references
ЛИТЕРАТУРА
1. Boldyreva E.V. Crystalline amino acids - a link between chemistry, materials sciences and biology. In: Models, Mysteries, and Magic of Molecules. Ed. by J.C.A. Boeyens, J.F. Ogilvie. springer Verlag: Berlin. 2007. P. 169-194.
2. Wesolowski M., Erecinska J. Relation between chemical structure of amino acids and their thermal decomposition J. Therm. Anal. Cal-orim. 2005. V. 82. N 2. P. 307-313. DOI: 10.1007/s10973-005-0895-z.
3. Gasit E. self-assembled peptide nanostructures: the design of molecular building blocks and their technological utilization. Chem. Soc. Rev. 2007. V. 36. N 8. P. 1263-1269. DOI: 10.1039/B605536M.
4. Iijima K., Tanaka K., Onuma S. Main conformer of gaseous glycine: molecular structure and rotational barrier from electron diffraction data and rotational constants. J. Mol. Struct. 1991. V. 246. N 3-4. P. 257-266. DOI: 10.1016/0022-2860(91)80132-N.
5. Lee Y.H., Jung J.W., Kim B., Butz P., Snoek L.C., Kroemer R.T., Simons J.P. Alanyl side chain folding in phenylalanine: conformational assignments through ultraviolet rotational band contour analysis. J. Phys. Chem. A. 2004. V. 108. N 1. P. 69-73. DOI: 10.1021/jp0368280.
6. Kaczor A., Reva ID., Proniewicz L.M., Fausto R. Importance of entropy in conformational equilibrium of phenylalanine: A matrixisolation infrared spectroscopy and density functional theory study. J. Phys. Chem. A. 2006. V. 110. N 7. P. 2360-2370. DOI: 10.1021/jp0550715.
7. Kaczor A., Reva I.D., Proniewicz L.M., Fausto R. Matrix-isolated monomelic tryptophan: Electrostatic interactions as nontrivial factors stabilizing conformers. J. Phys. Chem. A. 2007. V. 111. N 15. P. 2957-2965. DOI: 10.1021/jp070097c.
аминокислот по данным газовой электронографии
r, Â; Z, ° Gly [4] Ala [14] Trp
C-Cm — 1.536(11) 1.540(5)
C-N 1.467(5) 1.453(2) 1,474(3)
C-C 1.526(3) 1.527(11) 1.541(5)
C=O 1.205(1) 1.197(1) 1.195(6)
C-O 1.355(2) 1.341(2) 1.329(6)
C-C-Cm — 111.9(2) 109.6(13)
C-C-N 112.1(5) 112.9(3) 111.3 (18)
Cm-C-N — 110.0(2) 118.2(50)
C-C=O 125.1(5) 125.7(3) 123.9(25)
C-C-O 111.6(5) 110.3(2) 110.0(25)
8. Snoek L.C., Robertson E.G., Kroemer R.T., Simons J.P. Conformational landscapes in amino acids: infrared and ultraviolet ion-dip spectroscopy of phenylalanine in the gas phase. Chem. Phys. Lett. 2000. V. 321. N 1-2. P. 49-56. DOI: 10.1016/S0009-2614(00)00320-1.
9. Lesarri A., Cocinero E.J., Lopez J.C., Alonso J.L. Shape of4(S)-and 4(R)-hydroxyproline in gas phase. J. Am. Chem. Soc. 2005. V. 127. N 8. P. 2572-2579. DOI: 10.1021/ja045955m.
10. Petrov V.M., Girichev G.V., Oberhammer H., Giricheva N.J., Bardina A.V., Petrova V.N., Ivanov S.N. Molecular structure and conformations of 2-nitrobenzenesulfonyl fluoride: Gas-phase electron diffraction and quantum chemical calculations study. J. Mol. Struct. 2010. V. 978. N 1-3. P. 97-103. DOI: 10.1016/j.mol-struc.2009.12.005.
11. Vogt N., Khaikin L.S., Grikina O.E., Rykov A.N. A benchmark study of molecular structure by experimental and theoretical methods: Equilibrium structure of uracil from gas-phase electron diffraction data and coupled-cluster calculations. J. Mol. Struct. 2013. V. 1050. P. 114-121. DOI: 10.1016/j.molstruc.2013.07.014.
12. Marochkin Ы, Altova E.P., Rykov A.N., Shishkov IF. Molecular structure of tryptamine in gas phase according to gas electron diffraction method and quantum chemistry calculations. J. Mol. Struct. 2017. V. 1148. P. 179-184. DOI: 10.1016/j.molstruc.2017.07.016.
13. Belyakov A.V., Gureev M.A., Garabadzhiu A.V., Losev V.A., Rykov A.N. Determination of the molecular structure of gaseous proline by electron diffraction, supported by microwave and quantum chemical data. Struct. Chem. 2015. V. 26. N 5-6. P. 1489-1500. DOI: 10.1007/s11224-015-0589-5.
14. Iijima K., Nakano M. Reinvestigation of molecular structure and conformation of gaseous l-alanine by joint analysis using electron diffraction data and rotational constants. J. Mol. Struc. 1999. V. 485-486. P. 255-260. DOI: 10.1016/S0022-2860(99)00047-2.
15. Gahlmann A. Ultrafast electron diffraction: pulsed laser desorption enables time-resolved structural determination of thermally labile chromophores. Pasadena, California: California institute of technology. 2011. 274 p.
16. Frisch M.J., Trucks G.W., Schlegel H.B., Scuseria G.E., Robb M.A., Cheeseman J.R., Montgomery J.A., Vreven Jr.T., Kudin K.N., Burant J.C., Millam J.M., Iyengar S.S., Tomasi J., Barone V., Mennucci B., Cossi M., Scalmani G., Rega N., Petersson G.A., Nakatsuji H., Hada M., Ehara M., Toyota K., Fukuda R., Haseg-awa J., Ishida M., Nakajima T., Honda Y., Kitao O., Nakai H., Klene M., Knox X., Li J.E., Hratchian H.P., Cross J.B., Adamo C., Jaramillo J., Gomperts R., Stratmann R.E., Yazyev O., Austin A.J., Cammi R., Pomelli C., Ochterski J.W., Ayala P.Y., Morokuma K., Voth G.A., Salvador P., Dannenberg J.J., Zakrzewski V.G., Dapprich S., Daniels A.D., Strain M.C., Farkas O., Malick D.K., Rabuck A.D., Raghavachari K., Foresman J.B., Ortiz J.V., Cui Q., Baboul AG., Clifford S., Cioslowski J., Stefa-nov B.B., Liu G., Liashenko A., Piskorz P., Komaromi I., Martin R.L., Fox D.J., Keith T., Al-Laham M.A., Peng C.Y., Nanayak-kara A, Challacombe M., Gill P.M.W., Johnson B., Chen W., Wong M.W., Gonzalez C., Pople J.A Gaussian, Inc. Pittsburgh PA. 2003. 15. SDD. Gaussian 03. Revision B.03.
17. Weinhold F., Landis C.R. Discovering chemistry with natural bond orbitals. Wiley. 2012. 319 p.
18. Girichev G.V., Utkin A.N., Revichev Yu.F. Modernization of the EMR-100 electronograph for gas research. Pribory Tekhn. Ek-sperim. 1984. V. 27. N 2. P. 187-190 (in Russian). Гиричев Г.В., Уткин АН., Ревичев Ю.Ф. Модернизация электронографа ЭМР-100 для исследования газов. Приборы и техника эксперимента. 1984. Т. 27. № 2. С. 187-190.
19. Zakharov A.V., Zhurko G.A., Girichev E.G., Zhabanov Yu.A., Girichev G.V. Methods for processing electron diffraction data and their software implementations. Part I. Initial data processing. Izv. Vyssh. Uchebn. Zaved. Khim. Khim. Tekhnol. 2013. V. 56. N 1. P. 65-70 (in Russian).
Захаров А.В., Журко Г.А., Гиричев Е.Г., Жабанов Ю.А., Гиричев Г.В. Методики обработки электронографических данных и их программные реализации. Часть I. Первичная обработка электрнографических данных. Изв. вузов. Химия и хим. технология. 2013. Т. 56. Вып. 1. P. 65-70.
20. Tyunina V.V., Krasnov A.V., Tyunina E.Yu., Badelin V.G., Girichev G.V. Enthalpy of sublimation of natural aromatic amino acids determined by Knudsen's effusion mass spectrometric method. J. Chem. Thermodynamics. 2014. 74. P. 221-226. DOI: 10.1016/j.jct.2014.02.003.
21. Kaczor A., Reva ID., Proniewicz L.M., Fausto R. Matrix-Isolated Monomelic Tryptophan: Electrostatic Interactions as Nontrivial Factors Stabilizing Conformers. J. Phys. Chem. A. 2007. V. 111. N 15. P. 2957-2965. DOI: 10.1021/jp070097c.
22. Purushotham U., Narahari Sastry G. A comprehensive conformational analysis of tryptophan, its ionic and dimeric forms. J. Computational. Chem. 2014. V. 35. N 8. P. 595-610. DOI: 10.1002/jcc.23482.
23. Huang Z., Yu W., Lin Z. Exploration of the full conformational landscapes of gaseous aromatic amino acid phenylalanine: An ab initio study. J. Mol. Struct. Theochem. 2006. V. 758. N 2-3. P. 195-202. DOI: 10.1016/j.theochem.2005.10.043.
24. Weinhold F. Nature of H-bonding in cluster, liquids and enzymes: an ab initio, natural bond orbital perspective. J. Mol. Struct. Theochem. 1997. V. 398-399. P. 181-197. DOI: 10.1016/S0166-1280(96)04936-6.
25. Zakharov A.V., Zhurko G.A., Girichev G.V. Methods for processing electron diffraction data and their software implementations. Part II. Calculation of root-mean-square amplitudes of vibrations and vibrational corrections to internuclear distances. Least-squares analysis. Izv. Vyssh. Uchebn. Zaved. Khim. Khim. Tekhnol. 2013. V. 56. N 1. P. 71-74 (in Russian).
Захаров А.В., Журко Г.А., Гиричев Г.В. Методики обработки электронографических данных и их программные реализации. Часть II. Расчет среднеквадратичных амплитуд колебаний и колебательных поправок к межъядерным расстояниям. МНК-анализ. Изв. вузов. Химия и хим. технология. 2013. Т. 56. № 1. С. 71-74.
26. Schemer S., Kar T., Pattanayak J. Comparison of various types of hydrogen bonds involving aromatic amino acids. J. Am. Chem. Soc. 2002. V. 124. N 44. P. 13257-13264. DOI: 10.1021/ja027200q.
27. Sanz M.E., Cabezas C., Mata S., Alonsod J.L. Rotational spectrum of tryptophan. J. Chem. Phys. 2014. V. 140. N 20. P. 2043081 - 204308-6. DOI: 10.1063/1.4876001.
28. Snoek L.C., Kroemer R.T., Hockridge M.R., Simons J.P. Con-formational landscapes of aromatic amino acids in the gas phase: Infrared and ultraviolet ion dip spectroscopy of tryptophan. Phys. Chem. Chem. Phys. 2001. V. 3. N 10. P. 1819-1826. DOI: 10.1039/b101296g.
29. Ramaekers R., Pajak J., Rospenk M., Maes G. Matrix-isolation FT-IR spectroscopic study and theoretical DFT(B3LYP)/6-31++ G** calculations of the vibrational and conformational properties of tyrosine. Spectrochim. Acta. Part A. Molec. Boio-molec. Spectroscopy. 2005. V. 61. N 7. P. 1347-1356. DOI: 10.1016/j.saa.2004.10.003.
30. Lambie B., Ramaekers R., Maes G. On the contribution of intramolecular H-bonding entropy to the conformational stability of alanine conformations. Spectrochim. Acta. Part A. Molec. Biomolec. Spectroscopy. 2003. V. 59. N 6. P. 1387-1397. DOI: 10.1016/S1386-1425(02)00353-0.
31. Tyunina V.V., Krasnov A.V., Badelin V. G., Girichev G.V. Enthalpy of sublimation of hydroxyl-containing amino acids: Knud-sen's effusion mass spectrometric study. J. Chem. Thermodynamics. 2016. 98. P. 62-70. DOI: 10.1016/j.jct.2016.02.021.
Поступила в редакцию (Received) 11.03.2019 Принята к опубликованию (Accepted) 26.09.2019