Квантово-механічний конформаційний аналіз молекули 2'-дезоксицитидилової кислоти — структурної ланки ДНК
На рівні теорії MP2/6-311G++(d,p)//DFT B3LYP/6-31G(d,p) вперше отримано повне сімейство конформерів молекули 2'-дезоксицитидилової кислоти. Проаналізовано розподіли значень її торсійних кутів і порівняно з такими в молекулі 2'-дезоксицитидину. Виявлено кореляцію між амплітудою псевдооберта...
Saved in:
| Published in: | Доповіді НАН України |
|---|---|
| Date: | 2010 |
| Main Authors: | , |
| Format: | Article |
| Language: | Ukrainian |
| Published: |
Видавничий дім "Академперіодика" НАН України
2010
|
| Subjects: | |
| Online Access: | https://nasplib.isofts.kiev.ua/handle/123456789/30409 |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| Journal Title: | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
| Cite this: | Квантово-механічний конформаційний аналіз молекули 2'-дезоксицитидилової кислоти — структурної ланки ДНК / Т.Ю. Ніколаєнко, Д.М. Говорун // Доп. НАН України. — 2010. — № 9. — С. 173-184. — Бібліогр.: 15 назв. — укр. |
Institution
Digital Library of Periodicals of National Academy of Sciences of Ukraine| _version_ | 1860235661647282176 |
|---|---|
| author | Ніколаєнко, Т.Ю. Говорун, Д.М. |
| author_facet | Ніколаєнко, Т.Ю. Говорун, Д.М. |
| citation_txt | Квантово-механічний конформаційний аналіз молекули 2'-дезоксицитидилової кислоти — структурної ланки ДНК / Т.Ю. Ніколаєнко, Д.М. Говорун // Доп. НАН України. — 2010. — № 9. — С. 173-184. — Бібліогр.: 15 назв. — укр. |
| collection | DSpace DC |
| container_title | Доповіді НАН України |
| description | На рівні теорії MP2/6-311G++(d,p)//DFT B3LYP/6-31G(d,p) вперше отримано повне сімейство конформерів молекули 2'-дезоксицитидилової кислоти. Проаналізовано розподіли значень її торсійних кутів і порівняно з такими в молекулі 2'-дезоксицитидину. Виявлено кореляцію між амплітудою псевдообертання фуранозного кільця і syn-орієнтацією нуклеотидної основи. Охарактеризовано ''тонкі'' конформаційні параметри, що описують непланарність кільця нуклеотидної основи і деяких його бічних груп.
The full set of 2'-deoxycytidilic acid molecule conformers has been obtained by ab initio quantum mechanical calculations at the MP2/6-311G++(d,p)//DFT B3LYP/6-31G(d,p) theory level. The distributions of torsion angles are analyzed and compared with those of a 2'-deoxycytidine molecule. The strong correlation between the sugar puckering amplitude and the syn-orientation of the nucleotide base unit is revealed. 'Fine' conformational parameters defining a nonplanarity of the base ring and its side groups are characterized.
|
| first_indexed | 2025-12-07T18:23:32Z |
| format | Article |
| fulltext |
УДК 539.194
© 2010
Т.Ю. Нiколаєнко, член-кореспондент НАН України Д.М. Говорун
Квантово-механiчний конформацiйний аналiз молекули
2
′-дезоксицитидилової кислоти — структурної ланки
ДНК
На рiвнi теорiї MP2/6-311G++(d,p)//DFT B3LYP/6-31G(d,p) вперше отримано повне
сiмейство конформерiв молекули 2′-дезоксицитидилової кислоти. Проаналiзовано розпо-
дiли значень її торсiйних кутiв i порiвняно з такими в молекулi 2′-дезоксицитидину.
Виявлено кореляцiю мiж амплiтудою псевдообертання фуранозного кiльця i syn-орiєн-
тацiєю нуклеотидної основи. Охарактеризовано “тонкi” конформацiйнi параметри, що
описують непланарнiсть кiльця нуклеотидної основи i деяких його бiчних груп.
Молекула 2′-дезоксицитидилової кислоти (рис. 1) є нуклеотидом з електронейтральною
фосфатною групою в 5′-положеннi i цитозиновою основою. Бiологiчна важливiсть вивчення
її конформацiйного рiзноманiття зумовлена тим, що функцiонування ДНК, складовою час-
тиною якої є ця молекула, прямо пов’язане зi змiнами конформацiй її структурних ланок.
Крiм цього, порiвняння конформацiйних властивостей нуклеотиду i нуклеозиду [1] з однiєю
i тiєю ж нуклеотидною основою дозволить виявити внутрiшнньомолекулярнi механiзми, що
визначають основнi закономiрностi геометричної структури складових частин ДНК.
Iнтерес до вивчення 2′-дезоксицитидилової кислоти зумовлений ще й тим, що саме
GC-збагаченi дiлянки ДНК мають “схильнiсть” до B–Z переходу [2].
Вiдомi експериментальнi методи [3, 4] дослiдження окремих нуклеотидiв не в змозi на
сьогоднi дати вичерпної iнформацiї про конформацiйнi властивостi цих об’єктiв, оскiльки
в експериментах спостерiгаються лише кiлька енергетично “найвигiднiших” конформацiй.
У той же час провести аналiз усiх можливих конформерiв (включно з тими, що реалiзуються
в складi ДНК [5]), незалежно вiд їхньої вiдносної енергiї, дозволяють методи обчислюваль-
ної квантової механiки.
Рис. 1. Структура молекули 2′-дезоксицитидилової кислоти та позначення її атомiв
ISSN 1025-6415 Доповiдi Нацiональної академiї наук України, 2010, №9 173
Об’єкти i методи дослiдження. Для опису конформацiї молекули 2′-дезоксицити-
дилової кислоти (див. рис. 1) використано вiсiм конформацiйних параметрiв: загально-
прийнятi [6] торсiйнi змiннi β, γ, δ1, ε i χ = O4′C1′N1C2 та торсiйнi кути α′ = C5′O5′POP,
ζ1 = O5′POP1HP1 i ζ2 = O5′POP2HP2. Для трьох атомiв кисню фосфатної групи використано
такi позначення: не з’єднаний з атомом водню — OP, два iнших — вiдповiдно OP1 i OP2, якщо
дивитися вздовж зв’язку O5′−P i рухатися вiд OP за годинниковою стрiлкою. Зазначимо,
що введений таким чином кут α′ вiдрiзняється вiд номенклатурного кута α [6] на постiйну
величину. Кути ζ1 i ζ2 описують орiєнтацiї гiдроксильних груп OP1HP1 i OP2HP2 вiдповiдно.
Таким чином, кожен конформер описувався значеннями восьми торсiйних кутiв.
Пошук повного сiмейства конформерiв молекули 2′-дезоксицитидилової кислоти про-
водили таким чином. До кожного з конформерiв молекули 2′-дезоксицитидину [1] замiсть
атома H5′ приєднували фосфатну групу −POPOP1HP1OP2HP2 у кожному з трьох положень:
α′ = +60◦, −60◦ та 180◦ (кути ζ1 i ζ2 мали фiксованi значення −106◦ i 177◦ вiдповiдно). Гео-
метрiї одержаних у такий спосiб 273 конформерiв молекули 2′-дезоксицитидилової кислоти
було оптимiзовано на рiвнi теорiї MP2/6-311G++(d,p)//DFT B3LYP/6-31G(d,p) з викорис-
танням програмного пакета Gaussian 03w для платформи win32 [7] iз налаштуваннями то-
чностi обчислень за замовчуванням. Для кожного конформера розраховували в гармонiй-
ному наближеннi частоти нормальних коливань. Цi данi використовували для вилучення
нестiйких конформацiй (у коливальних спектрах яких присутня хоча б одна уявна часто-
та) та для розрахунку енергiї Гiббса конфоромерiв (при значеннях температури 298,15 К
i тиску 1 атм).
У результатi було одержано 266 стiйких конформерiв. Проведений аналiз показав, що
розподiл значень кута χ є бiмодальним (усi його значення лежать у секторах 70◦±13◦ (syn),
i −142◦±35◦ (anti)), а кутiв γ i ε — тримодальним (усi значення γ попадають у сектори 45◦±
±25◦ (gouche+), −72◦±46◦ (gouche−) i 173◦±18◦ (trans), а ε — у сектори 53◦±23◦ (gouche+),
−79◦±21◦(gouche−) i −175◦±15◦ (trans)). Усi значення кута псевдообертання P знаходяться
в областях 167◦±48◦ (S-сiмейство) i 0◦±55◦ (N -сiмейство). Розподiл кута β виявився майже
неперервним, крiм областi −45◦. . .+50◦, в яку не потрапило жодного значення. У розподiлах
значень кожного з кутiв α′, ζ1 i ζ2 наявнi один сектор i одна майже суцiльна область значень.
Таким чином, уже на першому етапi у комфомерах нуклеотиду було виявлено змiни
статистичного розподiлу значень кута β порiвняно з конформерами вiдповiдного нуклео-
зиду [1].
Пiсля цього будь-якi два конформери одержаного сiмейства, у яких значення кутiв χ, γ
i ε попадали в один i той самий сектор, кута α′ — вiдрiзнялися менше нiж на 5◦, а кутiв β,
ζ1 i ζ2 — менше нiж на 60◦, вважалися однаковими, що звузило кiлькiсть “унiкальних”
конформерiв до 215. Це сiмейство доповнювали в такий спосiб, щоб у ньому при однакових
iнших параметрах iснували конформери, у яких кут χ набуває як syn-, так i anti -значень,
а кут γ — кожного з трьох значень: gouche+, gouche− i trans, для чого було згенеровано 587
нових конформерiв. Їх геометрiї було оптимiзовано на тому ж рiвнi теорiї DFT B3LYP/6-
31G(d,p), що i на першому етапi.
Доповнене таким чином сiмейство мiстило 778 оптимiзованих стiйких конформерiв. Вра-
ховуючи неминучi обчислювальнi похибки, що виникають у процесi оптимiзацiї геометрiї,
будь-якi два конформери, у яких вiдповiднi значення кутiв χ, α′, β, γ, ε, ζ1 i ζ2 вiдрiз-
нялися не бiльше нiж на 10◦, а кута P — не бiльше нiж на 5◦, вважалися однаковими
1Кут δ однозначно пов’язаний з кутом P i амплiтудою νmax псевдообертання фуранозного кiльця [6].
174 ISSN 1025-6415 Reports of the National Academy of Sciences of Ukraine, 2010, №9
(вiдмiннiсть мiж електронними енергiями таких “однакових” конформерiв не перевищила
0,14 ккал/моль, що помiтно менше, нiж kT ∼ 0,6 ккал/моль при T = 298 К). Пiсля вилу-
чення “однакових” структур залишилося 326 конформерiв.
Для того щоб усунути неповноту, обумовлену тим, що при побудовi вихiдного сiмейства
орiєнтацiї гiдроксильних груп фосфатної групи (кути ζ1 i ζ2) були фiксованi i для них не
проводився “перебiр”, для кожного з наявних стiйких конформерiв створювався новий, що
вiдрiзнявся замiною ζ1 ↔ ζ2
2. Геометрiї новостворених конформерiв оптимiзувалися i до-
лучалися до наявного сiмейства, пiсля чого виконувалася така ж процедура вiдсiювання
однакових структур, як i на попередньому етапi. У результатi двократного проведення та-
кого усунення неповноти кiлькiсть унiкальних конформерiв склала 481.
Останнiм кроком було доповнення цього сiмейства в такий спосiб, щоб при однакових
iнших параметрах у ньому iснували конформери, у яких кут ε набуває кожного з трьох
значень: gouche+, gouche− i trans. Пiсля проведення оптимiзацiї геометрiї новостворених
конформерiв i вилучення однакових було одержано множину з 613 конформерiв, яку вва-
жаємо повною.
Результати дослiдження та їхнiй аналiз. Структурнi i енергетичнi параметри ви-
браних3 конформерiв молекули 2′-дезоксицитидилової кислоти наведено в табл. 1. Значення
їхнiх енергiй Гiббса (при T = 298,15 К i P = 1 атм) лежать у дiапазонi 0–14,1 ккал/моль;
глобальний мiнiмум i максимум енергiї Гiббса — вiдповiдно конформери 1 i 40 (див. табл. 1).
Syn-орiєнтацiя (χ ∈ 75◦ ± 19◦) нуклеотидної основи спостерiгається у 291 конформера,
anti - (χ ∈ −148◦ ± 42◦) — у 318 конформерiв, “промiжна” (χ = 142◦, 145◦, 141◦, 104◦) —
у конформерiв 10, 11, 14 i 364. “Пiвнiчну” (P ∈ 0◦ ± 89,9◦) конформацiю фуранозного
кiльця мають 292 конформери (з яких 200 — класичну C3′-endo (P ∈ 0◦. . .36◦)), “пiвденну”
(P ∈ 180◦ ± 90◦) — 321 корформер (з яких 141 — класичну C2′-endo (P ∈ 144◦. . .180◦)).
Кореляцiї торсiйних кутiв. Доволi велика кiлькiсть одержаних конформерiв дозво-
лила провести статистичний аналiз — обчислити за методикою [8] коефiцiєнти лiнiйної
кореляцiї мiж усiма конформацiйними параметрами молекули 2′-дезоксицитидилової кис-
лоти. Серед 45 цих коефiцiєнтiв лише 4 мають значення, бiльшi вiд 0,2: r(δ, P ) = 0,96,
r(P, νmax) = 0,59, r(χ, νmax) = 0,34 i r(ζ1, ζ2) = 0,30.
Близькiсть r(δ, P ) до одиницi пов’язана з тим, що кут δ однозначно визначається ку-
том P i амплiтудою νmax псевдообертання фуранозного кiльця, причому всi значення νmax
лежать у досить вузькому дiапазонi (8◦. . .43◦: найбiльший νmax виявлено в конформерi 15
(див. табл. 1)).
Менш очiкуваним є те, що глiкозидний кут χ найбiльше корелює саме з νmax. Це зумов-
лено тим (рис. 2, а), що “пласкi” цукри (0◦ < νmax < 30◦; найменшi значення νmax виявлено
у конформерiв 26, 27 i 33; див. табл. 1) завжди супроводжуються syn-орiєнтацiєю нуклео-
тидної основи. Аналогiчна закономiрнiсть має мiсце i для нуклеозиду (див. рис. 2, б ), але
не була ранiше вiдмiчена [1].
Вiрогiдною фiзичною причиною появи “пласких” конформацiй фуранозного кiльця в ну-
клеотидах i нуклеозидах є наявнiсть внутрiшньомолекулярних водневих зв’язкiв з енергiя-
ми, достатнiми для компенсацiї енергетичної “невигiдностi” затiнених конформацiй C–C
2Такий метод було обрано з огляду на те, що в розподiлах кутiв ζ1, ζ2 вiдсутнi будь-якi чiтко окресленi
сектори, а отже, не було можливостi вибрати скiльки-небудь обгрунтований крок перебору їх значень.
3Файли структур повного сiмейства конформерiв можна отримати у авторiв за запитом.
4“Промiжнi” значення глiкозидного кута обумовленi сильними водневими зв’язками (OP1−H · · ·N3 або
OP2−H · · ·N3) мiж гiдроксильними групами фосфатної групи i атомами нуклеотидної основи.
ISSN 1025-6415 Доповiдi Нацiональної академiї наук України, 2010, №9 175
зв’язкiв кiльця (яка i є причиною непланарностi фуранозного кiльця [6]). Так, у конформе-
ра 33 шляхом аналiзу просторового розподiлу густини електронного заряду було виявлено
водневi зв’язки O3′H · · ·OP i C5′H · · ·O2 (рис. 3).
Розподiли значень торсiйних кутiв. На рис. 4, а наведено розподiли значень конформа-
цiйних параметрiв одержаних конформерiв молекули 2′-дезоксицитидилової кислоти. Для
кутiв χ, γ i ε вони виявляють чiтку “секторну” структуру: бiмодальну для χ (усi значення
Таблиця 1. Конформацiйнi (град) i енергетичнi (ккал/моль) характеристики вибраних конформерiв моле-
кули 2′-дезоксицитидилової кислоти
Конформер β γ δ ε χ α′ ζ1 ζ2 P νmax
∆G,
ккал/моль
1 −90,1 55,5 81,0 −92,3 66,7 9,3 −44,7 −74,7 29,0 33,0 0,00
2 −90,3 55,2 81,1 −92,6 67,3 9,0 −43,2 121,8 28,4 32,8 0,32
3 148,1 46,2 100,0 −158,0 88,5 −48,0 −139,3 57,1 338,4 28,3 2,08
4 110,7 41,9 87,7 −153,6 65,3 −7,6 76,1 51,7 45,4 29,1 2,08
5 106,2 47,9 89,0 −102,6 66,8 85,8 53,8 −24,9 44,0 31,1 2,13
6 106,2 49,3 87,8 −148,2 67,4 85,0 52,6 −24,0 44,9 30,1 2,17
7 112,1 41,6 88,0 −150,4 65,3 −7,8 −122,9 49,9 44,9 29,0 2,24
8 111,0 39,5 89,3 −93,2 64,6 −5,4 80,1 52,7 44,4 30,2 2,40
9 112,1 39,3 89,3 −96,7 64,7 −6,1 −122,7 51,5 43,5 30,0 2,43
10 162,7 43,3 150,9 174,8 141,7 −43,8 −150,3 72,9 224,6 35,2 2,72
11 165,9 40,7 150,7 55,3 145,1 −42,3 −149,8 71,9 229,2 39,1 2,93
12 −84,7 179,9 140,0 172,6 −149,4 −34,5 −55,8 105,5 153,8 38,1 2,96
13 132,7 49,8 133,9 −176,9 69,1 −15,7 −131,0 54,2 147,5 36,1 2,97
14 162,8 43,2 155,2 −71,9 140,6 −43,4 −150,6 72,5 224,6 35,5 3,52
15 75,4 51,9 118,8 −175,8 −111,2 −116,6 −64,1 133,0 125,8 42,7 3,85
16 163,4 58,5 146,7 170,1 −172,7 68,6 53,0 145,7 189,8 30,4 5,51
17 161,8 51,5 141,2 168,6 −161,3 51,4 −104,5 156,4 161,7 34,6 5,95
18 −172,4 56,2 137,2 171,9 −144,9 −172,3 −106,1 176,9 153,6 35,4 6,42
19 −175,0 55,2 138,2 170,3 −156,1 −174,7 −106,9 175,4 155,4 35,5 6,49
20 −173,0 56,6 82,3 −171,2 −165,8 56,3 −110,1 157,1 13,9 34,7 6,60
21 177,2 58,6 81,4 −171,1 −163,8 −170,7 −104,2 176,2 14,4 35,5 6,62
22 −177,2 56,5 81,3 −175,1 −165,0 −179,5 −175,8 103,2 14,0 35,6 6,74
23 153,4 45,3 146,9 170,0 −176,5 168,2 −78,4 143,5 187,5 31,4 6,78
24 157,8 50,6 91,6 −168,1 −179,0 62,0 42,8 147,7 351,8 31,6 6,87
25 163,4 58,5 150,8 −68,3 −172,5 68,3 53,3 146,0 192,3 30,5 6,88
26 −56,0 −55,7 118,1 −58,2 70,6 72,5 −109,1 156,1 344,4 8,9 6,89
27 −56,9 −54,5 118,7 −58,6 70,3 68,5 32,6 149,1 342,5 8,3 6,92
28 158,9 50,9 140,5 168,5 −159,2 174,1 179,5 106,2 160,0 34,8 7,06
29 164,9 51,0 144,8 −68,3 −158,1 49,5 −102,2 157,4 163,8 33,5 7,32
30 168,8 50,0 152,0 −69,1 −177,2 36,5 −91,6 −46,4 199,9 30,9 7,43
31 136,8 48,0 143,3 −67,5 −160,3 −56,8 −153,0 −45,0 159,5 34,2 7,49
32 −170,1 59,7 81,4 −169,0 −163,9 58,8 −134,3 −59,1 15,8 35,3 7,83
33 −56,7 −55,6 118,6 −56,3 70,1 70,2 −132,4 −59,6 345,7 8,1 8,45
34 −123,5 55,8 142,1 −67,9 −137,4 −19,9 −174,4 97,6 157,6 34,3 8,57
35 −158,3 −177,5 85,7 −176,4 −172,1 −56,4 −158,4 103,8 10,0 33,5 9,07
36 173,4 37,0 147,6 54,9 104,1 148,5 122,0 −59,2 217,7 34,6 11,29
37 60,2 −121,9 84,2 −168,2 71,5 −97,1 −82,0 137,2 40,6 35,3 11,86
38 55,5 −120,5 86,2 −170,9 69,9 −91,3 −85,9 −38,3 38,6 33,5 12,22
39 56,0 −114,3 81,9 172,8 −170,1 −88,9 −89,9 −38,3 22,4 38,6 13,10
40 −177,9 44,3 130,6 60,6 71,4 51,4 −101,5 159,9 150,8 35,5 14,13
176 ISSN 1025-6415 Reports of the National Academy of Sciences of Ukraine, 2010, №9
Рис. 2. Розподiл конформацiйних параметрiв у координатах (νmax, χ) для молекул 2′-дезоксицитидилової
кислоти (а) та 2′-дезоксицитидину [1] (б )
попадають в областi syn i anti) i тримодальну для γ та ε (усi значення попадають в областi
gouche+, gouche− i trans). У конформерiв 37, 38, 39 кут γ набуває “нетипових” значень,
близьких до −120◦. В iншому розподiли значень кутiв γ, ε i χ якiсно не вiдрiзняються вiд
вiдомих у нуклеозидах [1].
Секторна структура вiдсутня в розподiлi значень кута β, за винятком областi [−45. . .
+45◦] (“прогалини”), у яку не потрапляє жодного значення5. Її походження обумовлене,
на нашу думку, тим, що при β ≈ 0◦ атоми фосфатної групи опиняються занадто близько
до атомiв фуранозного кiльця, призводячи до значних стеричних ускладнень (йдеться про
вiдстанi, що не перевищують 1,0 Å, вiд атома OP (або OP1 чи OP2) до C2′ фуранозного
кiльця та/або до атомiв O2 i C6 нуклеотидної основи). Такi стеричнi обмеження, якi, на-
певне, мають мiсце i в полiнуклеотидах, принципово не могли бути виявленi при аналiзi
нуклеозиду [1], адже у нього вiдсутня фосфатна група, а кут β задає орiєнтацiю лише ато-
ма H5′ , вiдстань вiд якого до атомiв фуранозного кiльця i нуклеотидної основи є достатньо
великою при будь-яких значеннях β.
Розподiл значень кута α′ також не виявляє чiткої секторної структури, але характери-
зується “згущеннями” поблизу значень +45, −45 i 180◦ (останнє особливо виражене). Цей
розподiл є симетричним вiдносно α′ = 180◦, що обумовлено симетричним розташуванням
гiдроксильних груп OP1HP1 i OP2HP2 фосфатної групи вiдносно площини O5′POP. У скла-
дi полiнуклеотиду така симетрiя, очевидно, втрачається i можна очiкувати iстотної змiни
розподiлу значень кута α′ у них.
Проявом згаданої симетрiї фосфатної групи є також перехiд розподiлу ζ1 у розподiл ζ2
при вiднiманнi 120◦ вiд усiх значень першого. Це стає особливо помiтним, якщо орiєнтацiї
атомiв HP1, HP2 i O5′ фосфатної групи описувати, як це зроблено в [10], торсiйними кутами,
вiдрахованими вiд атома OP: τ1 = OPPOP1HP1, τ2 = OPPOP2HP2 i α′ вiдповiдно. У такому
випадку розподiли всiх цих трьох кутiв є якiсно однаковими (див. рис. 4, а).
Водночас, такий зв’язок розподiлiв є суто статистичною їх властивiстю, адже в кожно-
му окремо взятому конформерi значення кутiв ζ1 i ζ2 є майже незалежними (коефiцiєнт
лiнiйної кореляцiї r(ζ1, ζ2) = 0,30).
5Така ж особливiсть розподiлу кута β була виявлена у дослiдженнi, виконаному з використанням потен-
цiалу твердих сфер, який враховував лише сили мiжатомного вiдштовхування [9].
ISSN 1025-6415 Доповiдi Нацiональної академiї наук України, 2010, №9 177
Рис. 3. Структура вибраних конформерiв молекули 2′-дезоксицитидилової кислоти (цифри вiдповiдають
номеру конформера в табл. 1). Пунктиром позначено водневi зв’язки
Кути α′, ζ1, ζ2 є ступенями вiльностi, характерними для нуклеотиду, у нуклеозидiв вони
вiдсутнi.
У розподiлi значень кута псевдообертання P спостерiгаються згущення поблизу 18◦ (що
вiдповiдає конформацiї C3′-endo) i 180◦ (конформацiї C2′-endo та C3′-exo). Лише нечисленнi
значення потрапляють в East- (P = 90◦ ± 45◦) та West- (P = 270◦ ± 45◦) областi. “Прога-
лини” обумовленi тим, що значення P = 70◦ i P = 270◦ вiдповiдають бар’єрам у циклi
псевдообертання фуранозного кiльця [6 (рис. 47), 11].
Деталiзацiї потребує i поняття syn- i anti -конформацiї. Вiдповiдними дiапазонами до-
цiльно називати не формальнi геометричнi областi, а одержанi сектори в розподiлi кута χ,
а саме — syn = {45◦ < χ < 105◦}, anti = {−180◦ < χ < −100◦}.
Конформацiйнi вiдмiнностi нуклеотидiв i нуклеозидiв. У нуклеозидi [1] розподiли всiх
торсiйних кутiв мають чiтко виражену “секторну” структуру. Приєднання фосфатної групи
iстотно змiнило цю властивiсть: розподiли кутiв β, ζ1, ζ2 i α′ стали майже неперервними
(за винятком “прогалин”, походження яких описано нижче). Ми пов’язуємо цi змiни з iсну-
ванням “iєрархiї жорсткостей” рiзних торсiйних ступенiв вiльностi. Так, ширина кожного
178 ISSN 1025-6415 Reports of the National Academy of Sciences of Ukraine, 2010, №9
Рис. 4. Розподiл конформерiв молекули 2′-дезоксицитидилової кислоти за конформацiйними параметра-
ми (а) та енергiями Гiббса (б — при T = 0 К; в — при T = 298,15 К)
ISSN 1025-6415 Доповiдi Нацiональної академiї наук України, 2010, №9 179
iз секторiв є результатом балансу мiж крутильною жорсткiстю щодо вiдносного повороту
атомних груп “1” i “2” навколо одинарного хiмiчного зв’язку 1–2 (бар’єром обертання) та
“моментом сил”, прикладеним до цих груп. Величина останнього визначається як силовими
характеристиками — енергiями водневих зв’язкiв, у яких беруть участь атоми груп 1 i 2,
так i геометричними — розмiрами атомних груп 1 i 2.
Другий чинник особливо важливий у випадку кута β, адже в молекулi 2′-дезоксици-
тидилової кислоти крутильний момент щодо осi C5′−O5′ обумовлений силами притягання,
прикладеними до атомiв водню гiдроксильних груп фосфатної групи i зумовленими участю
цих атомiв у порiвняно сильних водневих зв’язках. Плече цих сил становить не менше 1,60 Å
(довжина зв’язку O5′−P), тодi як у 2′-дезоксицитидину [1] — не бiльше 0,97 Å (довжина
зв’язку O5′−H5′ : крутильний момент створювали сили, прикладенi до атома H5′).
Цiлком аналогiчно, наявнiсть “прогалин” у розподiлах кутiв ζ1, ζ2 свiдчить про те, що
орiєнтацiї фосфатної групи та гiдроксильних груп у її складi, що вiдповiдають ζ1 ∈ gouche+
i ζ2 ∈ gouche−, вiддiленi вiд iнших порiвняно високими потенцiальними бар’єрами обер-
тання.
Сектори в розподiлi значень кута γ стали ширшими (а прогалини — вужчими), нiж
у молекулi 2′-дезоксицитидину [1], але не зникли взагалi. Це свiдчить про велику (порiвняно
з кутом β) жорсткiсть цього ступеня вiльностi.
Енергетичнi властивостi. Розподiл одержаних конформерiв за енергiями Гiббса про-
iлюстровано на рис. 4, б (T = 0 К) та рис. 4, в (T = 298,15 К, p = 1 атм).
Глобальний мiнiмум енергiї Гiббса досягається у конформерiв 1 i 2 (див. табл. 1), струк-
тури яких (див. рис. 3) вiдрiзняються тiльки орiєнтацiєю гiдроксильної групи OP2HP2 (ку-
том ζ2). Цей мiнiмум є “квазiвиродженим”, оскiльки рiзниця енергiй G2−G1 = 0,3 ккал/моль
є малою, порiвняно з величиною енергетичної “щiлини” (1,8 ккал/моль) мiж конформерами
1 i 2 та енергетично найвигiднiшим з решти (3).
Враховуючи, що iмовiрнiсть pG реалiзацiї конформацiї молекули з енергiєю Гiббса G
пропорцiйна exp(−G/(kBT )), ми розраховували, що сумарна iмовiрнiсть експерименталь-
ного спостереження в газовiй фазi конформерiв 1 i 2 при T = 298,15 К становить 86%.
Сумарна iмовiрнiсть спостереження конформерiв 1–13 (див. табл. 1) становить 97%; кон-
формерiв з syn-орiєнтацiєю нуклеотидної основи — 98%; конформерiв з P ∈ [270 . . . 90◦] —
96% (з яких P ∈ C3′-endo — 86 %); з P ∈ C2′-endo — 2,9%.
Найменшу енергiю Гiббса за умови χ ∈ anti в одержаному сiмействi має конформер 12
(∆G = 3,0 ккал/моль (див. табл. 1)), за умови χ ∈ syn — конформер 1 (∆G=0,0 ккал/моль).
Конформер 1 одночасно є енергетично найвигiднiшим з P ∈ C3′-endo; найменшу енергiю
Гiббса за конформацiї цукру C2′-endo має конформер 12 (∆G = 3,0 ккал/моль).
Конформацiйна мiнливiсть. Крiм торсiйних ступенiв вiльностi, якi є найменш жорстки-
ми, було проаналiзовано так званi тонкi конформацiйнi змiннi, що описують непланарнiсть
кiльця нуклеотидної основи i деяких його бiчних груп. У табл. 2 наведено середнi значення
Таблиця 2. Ендоциклiчнi торсiйнi кути кiльця нуклеотидної основи (град)
Параметр
Кут
N3C2N1C6 C2N1C6C5 N1C6C5C4 C6C5C4N3 C5C4N3C2 C4N3C2N1
〈x〉 −0,5 0,3 0,0 −0,1 −0,1 0,4
δx 4,2 3,1 1,5 1,5 0,7 3,0
max(x− 〈x〉) 24,3 10,1 8,1 11,5 3,8 21,2
180 ISSN 1025-6415 Reports of the National Academy of Sciences of Ukraine, 2010, №9
〈x〉 (з урахуванням перiодичностi кутових змiнних [8]), дисперсiї δx та найбiльшi вiдхилення
значень вiд середнiх для ендоциклiчних торсiйних кутiв кiльця нуклеотидної основи, якi
є мiрою його непланарностi.
Згiдно з одержаними даними, кiльце нуклеотидної основи є плоским лише “в середньо-
му”: хоча середнi значення всiх шести ендоциклiчних кутiв близькi до нуля, але амплiтуда
змiни кута N3C2N1C6 досягає 24◦. Кут C5C4N3C2 є “найжорсткiшим”: чотири атоми, що
його утворюють, фактично завжди лежать в однiй площинi (виходять з неї не бiльше нiж
на C5C4N3C2 6 3,8◦). Саме вiд атомiв C5, C4, N3, C2 вiдраховувалися торсiйнi кути, що
описують вихiд бiчних радикалiв нуклеотидної основи з площини її кiльця (табл. 3).
Атоми N4, H5, H6 фактично завжди лежать у площинi кiльця основи (дисперсiї вiд-
повiдних торсiйних кутiв не перевищують 3◦, амплiтуди змiн — 10◦). Участь атома O2
у внутрiшньомолекулярних водневих зв’язках обумовлює його максимальне вiдхилення вiд
площини основи (кут N4N3C2O2) у 17◦.
Значно лабiльнiшими є атоми H41 i H42 амiногрупи: дисперсiї кутiв N3C4N4H41
i C5C4N4H42 перевищують модулi їх середнiх значень. Це свiдчить про те, що амiдний
фрагмент C4N4H41H42 не тiльки не є плоским [5] навiть “у середньому”, але його структура
в рiзних конформерах зазнає вiддзеркалення вiдносно площини кiльця основи.
Атом C1′ (глiкозидний зв’язок N1C1′) вiдхиляється вiд площини кiльця нуклеотидної
основи найбiльше. Так, максимальне вiдхилення кута N3C2N1C1′ вiд середнього значення
(−178,6◦) становило 34◦ (!) i було виявлено у конформера 36 (див. табл. 1, рис. 3). Ми пов’я-
зуємо це з наявнiстю внутрiшньомолекулярних водневих зв’язкiв OP1H · · ·N3 i C2′H · · ·O2.
У той же час довжина глiкозидного зв’язку змiнюється мало (див. δl/〈l〉 у табл. 4), порiв-
няно з довжинами хiмiчних зв’язкiв гiдроксильних груп OP2HP2 i OP1HP1 фосфатної групи.
Ця властивiсть також вiдрiзняє 2′-дезоксицитидилову кислоту вiд 2′-дезоксицитидину.
ДНК-подiбнi конформери. Конформацiя нуклеотиду в складi ДНК визначається мiнiму-
мом вiльної енергiї всiєї макромолекули. Ця енергiя є сумою енергiй окремих нуклеотидiв
та енергiї їхньої взаємодiї. Для того щоб розрiзнити внески цих складових, з одержано-
го сiмейства конформерiв молекули 2′-дезоксицитидилової кислоти нами було вибрано тi,
у яких значення торсiйних кутiв є подiбними до наявних у ДНК [12–15]. Так, значення
торсiйних кутiв у конформерiв 20, 21, 22, 24 i 32 (див. табл. 1) подiбнi до тих, що реа-
лiзуються в ДНК у A-I формi, 35 — у A-II формi, 16, 17, 18, 19, 23 i 28 — у B-I формi,
34 — у Z-I формi, 25, 29, 30 i 31 — у Z-II формi. При цьому структура конформера 25
подiбна до виявленої в кристалi моногiдрату дезоксицитидин-5′-монофосфату [3]. Енерге-
тично найвигiднiший (∆G = 5,5 ккал/моль) з ДНК-подiбних конформерiв (16) вiдповiдає
B-I формi ДНК.
Таким чином, методами обчислювальної квантової механiки вперше одержано повне сi-
мейство конформерiв молекули 2′-дезоксицитидилової кислоти. Показано, що залежно вiд
конформацiї вiдносна енергiя Гiббса молекули змiнюється в дiапазонi 0–14,1 ккал/моль,
причому iснує два енергетично “найвигiднiшi” конформери, якi вiддiленi вiд решти енерге-
тичною “щiлиною” в 1,8 ккал/моль.
Присутнiсть фосфатної групи iстотно вiдрiзняє розподiл кута β вiд наявного у 2′-дезо-
ксицитидинi, хоча не вносить якiсних змiн у розподiли кутiв γ i ε. При цьому спостерiгається
розширення меж секторiв, у яких знаходяться значення всiх торсiйних кутiв.
Визначальним є вплив на конформацiйнi властивостi молекули 2′-дезоксицитидило-
вої кислоти внутрiшньомолекулярних водневих зв’язкiв, зокрема тих, де донором протона
є OH-групи фосфатної групи. Якщо акцептором протона при цьому є атом N3 нуклеотидної
ISSN 1025-6415 Доповiдi Нацiональної академiї наук України, 2010, №9 181
Таблиця 3. Параметри, що описують вiдхилення бiчних радикалiв (град) вiд площини кiльця нуклеотидної основи
Параметр
Радикал
C1′ N4 H41 H42 O2 H5 H6
C5C6N1C1′ N3C2N1C1′ N3C5C4N4 N3C4N4H41 C5C4N4H42 N4N3C2O2 N3N4C5H5 N4C5C6H6
〈x〉 178,6 −178,6 179,2 −5,3 8,7 −179,6 −179,3 −179,6
δx 3,5 4,6 1,0 7,3 11,8 2,5 2,1 1,7
max(x − 〈x〉) 20,4 34,1 2,3 18,0 28,8 17,4 8,5 10,4
182
IS
S
N
1
0
2
5
-6
4
1
5
R
epo
rts
o
f
th
e
N
a
tio
n
a
l
A
ca
d
em
y
o
f
S
cien
ces
o
f
U
kra
in
e,
2
0
1
0
,
№
9
Таблиця 4. Середнi значення i дисперсiї (Å) довжин деяких хiмiчних зв’язкiв
Параметр
Зв’язок
OP2HP2 OP1HP1 O5′P POP1 C1′N1 POP2 O4′C1′
〈l〉 0,972 0,972 1,612 1,611 1,475 1,610 1,421
δl 0,009 0,008 0,012 0,011 0,010 0,010 0,009
max(l − 〈l〉) 0,033 0,027 0,031 0,035 0,035 0,030 0,037
δl/〈l〉, % 0,88 0,85 0,77 0,69 0,67 0,63 0,63
Пр и м i т ка . δl =
√
〈(l − 〈l〉)2〉.
основи, то такi зв’язки призводять до “нетипових” значень глiкозидного кута χ i помiтного
вiдхилення глiкозидного зв’язку вiд площини кiльця нуклеотидної основи.
Виявлено кореляцiю мiж планарнiстю фуранозного кiльця i syn-орiєнтацiєю нуклеотид-
ної основи.
Автори висловлюють щиру вдячнiсть Р.О. Журакiвському (Iнститут молекулярної
бiологiї i генетики НАН України) за наданi поради щодо пошуку повного сiмейства кон-
формерiв; а також корпорацiї “Gaussian Inc.” (США) за люб’язно наданий Д.М. Говоруну
грант — пакет “GAUSSIAN03” для платформи win32.
1. Журакiвський Р.О., Говорун Д.М. Вичерпний конформацiйний аналiз канонiчного нуклеозиду 2′-де-
зоксицитидину квантово-механiчним методом функцiоналу густини // Фiзика живого. – 2006. – 14,
№ 3. – С. 35–48.
2. Fuertes M., Cepeda V., Alonso C., Perez J. Molecular Mechanisms for the B-Z Transition in the Example
of Poly[d(G-C)d(G-C)] Polymers. A Critical Review // Chem. Rev. – 2006. – 106. – P. 2045. – 2064.
3. Viswamitra M.A., Reddy B. S., Lin G.H.Y., Sundaraligam M. Stereochemistry of Nucleic Acids and Their
Constituents. XVII. la Crystal and Molecular Structure of Deoxycytidine 5′-Phosphate Monohydrate. lb A
Possible Puckering for the Furanoside Ring in B-Deoxyribonucleic Acidic // J. Amer. Chem. Soc. – 1971. –
93, issue 18. – P. 4565–4573.
4. Pandit J., Seshardi T.P., Viswamitra M. A. Structure of Disodium Deoxycytidine 5′-Phosphate Hepta-
hydrate, C9H12N3O7P
2− · 2Na+ · 7H2O // Acta Cryst. – 1983. – C39. – P. 342–345.
5. Shishkin O.V., Gorb L., Zhikol O.A., Leszczynski J. Conformational Analysis of Canonical 2-Deoxyri-
bonucleotides. 1. Pyrimidine Nucleotides // J. Biomol. Structure & Dynamics. – 2004. – 21, Issue 4. –
P. 537–553.
6. Зенгер В. Принципы структурной организации нуклеиновых кислот. – Москва: Мир, 1987. – 584 с.
7. Frisch M. J., Trucks G.W., Schlegel H.B. et al. Gaussian 03, Revision С. 02 / Gaussian, Inc., Wallingford
CT, 2004.
8. Kitamura K., Wakahara A., Mizuno H. et al. Conformationally “Concerted” Changes in Nucleotide Struc-
tures. A New Description Using Circular Correlation and Regression Analyses // J. Amer. Chem. Soc. –
1981. – 103. – P. 3899–3904.
9. Murthy V. L., Srinivasan R., Draper D.E., Rose G.D. A Complete Conformational Map for RNA //
J. Mol. Biol. – 1999. – 291. – P. 313–327.
10. Булавiн Л.А., Нiколаєнко Т.Ю., Говорун Д.М. Структурна нежорсткiсть молекули ортофосфорної
кислоти: неемпiричне квантово-механiчне дослiдження // Доп. НАН України. – 2009. – № 10. – С. 80–
87.
11. Foloppe N., Nilsson L., MacKerell A. Ab initio Conformational Analysis of Nucleic Acid Components:
Intrinsic Energetic Contributions to Nucleic Acid Structure and Dynamics // Biopolymers. – 2002. – 61,
Issue 1. – P. 61–76.
12. Berman H. Crystal studies of B-DNA: The answers and the questions // Ibid. – 1997. – 44, Issue 1. –
P. 23–44.
13. Schneide B., Neidle St., Berman H.M. Conformations of the Sugar-Phosphate Backbone in Helical DNA
Crystal Structures // Ibid. – 1997. – 42. – P. 113–124.
14. Svozil D., Kalina J., Omelka M., Schneider B. DNA conformations and their sequence preferences // Nucl.
Acids Res. – 2008. – 36, No 11. – P. 3690–3706.
ISSN 1025-6415 Доповiдi Нацiональної академiї наук України, 2010, №9 183
15. Gessner R.V., Frederick C.A., Quigley G. J. et al. The Molecular Structure of the Left-handed Z-DNA
Double Helix at 1.0-Å Atomic Resolution // J. Biol. Chem. – 1989. – 264, No 14. – P. 7921–7935.
Надiйшло до редакцiї 19.11.2009Київський нацiональний унiверситет
iм. Тараса Шевченка
Iнститут молекулярної бiологiї i генетики
НАН України, Київ
T.Yu. Nikolaienko, Corresponding Member of the NAS of Ukraine D.M. Hovorun
Quantum-mechanical conformational analysis of 2′-deoxycytidilic acid
molecule — the DNA structural unit
The full set of 2′-deoxycytidilic acid molecule conformers has been obtained by ab initio quantum
mechanical calculations at the MP2/6–311G++(d,p)//DFT B3LYP/6–31G(d,p) theory level. The
distributions of torsion angles are analyzed and compared with those of a 2′-deoxycytidine molecule.
The strong correlation between the sugar puckering amplitude and the syn-orientation of the nucleo-
tide base unit is revealed. ‘Fine’ conformational parameters defining a nonplanarity of the base
ring and its side groups are characterized.
184 ISSN 1025-6415 Reports of the National Academy of Sciences of Ukraine, 2010, №9
|
| id | nasplib_isofts_kiev_ua-123456789-30409 |
| institution | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
| issn | 1025-6415 |
| language | Ukrainian |
| last_indexed | 2025-12-07T18:23:32Z |
| publishDate | 2010 |
| publisher | Видавничий дім "Академперіодика" НАН України |
| record_format | dspace |
| spelling | Ніколаєнко, Т.Ю. Говорун, Д.М. 2012-02-02T19:45:26Z 2012-02-02T19:45:26Z 2010 Квантово-механічний конформаційний аналіз молекули 2'-дезоксицитидилової кислоти — структурної ланки ДНК / Т.Ю. Ніколаєнко, Д.М. Говорун // Доп. НАН України. — 2010. — № 9. — С. 173-184. — Бібліогр.: 15 назв. — укр. 1025-6415 https://nasplib.isofts.kiev.ua/handle/123456789/30409 539.194 На рівні теорії MP2/6-311G++(d,p)//DFT B3LYP/6-31G(d,p) вперше отримано повне сімейство конформерів молекули 2'-дезоксицитидилової кислоти. Проаналізовано розподіли значень її торсійних кутів і порівняно з такими в молекулі 2'-дезоксицитидину. Виявлено кореляцію між амплітудою псевдообертання фуранозного кільця і syn-орієнтацією нуклеотидної основи. Охарактеризовано ''тонкі'' конформаційні параметри, що описують непланарність кільця нуклеотидної основи і деяких його бічних груп. The full set of 2'-deoxycytidilic acid molecule conformers has been obtained by ab initio quantum mechanical calculations at the MP2/6-311G++(d,p)//DFT B3LYP/6-31G(d,p) theory level. The distributions of torsion angles are analyzed and compared with those of a 2'-deoxycytidine molecule. The strong correlation between the sugar puckering amplitude and the syn-orientation of the nucleotide base unit is revealed. 'Fine' conformational parameters defining a nonplanarity of the base ring and its side groups are characterized. Автори висловлюють щиру вдячнiсть Р.О. Журакiвському (Iнститут молекулярної бiологiї i генетики НАН України) за наданi поради щодо пошуку повного сiмейства конформерiв; а також корпорацiї “Gaussian Inc.” (США) за люб’язно наданий Д.М. Говоруну грант — пакет “GAUSSIAN03” для платформи win32. uk Видавничий дім "Академперіодика" НАН України Доповіді НАН України Біофізика Квантово-механічний конформаційний аналіз молекули 2'-дезоксицитидилової кислоти — структурної ланки ДНК Quantum-mechanical conformational analysis of 2'-deoxycytidilic acid molecule — the DNA structural unit Article published earlier |
| spellingShingle | Квантово-механічний конформаційний аналіз молекули 2'-дезоксицитидилової кислоти — структурної ланки ДНК Ніколаєнко, Т.Ю. Говорун, Д.М. Біофізика |
| title | Квантово-механічний конформаційний аналіз молекули 2'-дезоксицитидилової кислоти — структурної ланки ДНК |
| title_alt | Quantum-mechanical conformational analysis of 2'-deoxycytidilic acid molecule — the DNA structural unit |
| title_full | Квантово-механічний конформаційний аналіз молекули 2'-дезоксицитидилової кислоти — структурної ланки ДНК |
| title_fullStr | Квантово-механічний конформаційний аналіз молекули 2'-дезоксицитидилової кислоти — структурної ланки ДНК |
| title_full_unstemmed | Квантово-механічний конформаційний аналіз молекули 2'-дезоксицитидилової кислоти — структурної ланки ДНК |
| title_short | Квантово-механічний конформаційний аналіз молекули 2'-дезоксицитидилової кислоти — структурної ланки ДНК |
| title_sort | квантово-механічний конформаційний аналіз молекули 2'-дезоксицитидилової кислоти — структурної ланки днк |
| topic | Біофізика |
| topic_facet | Біофізика |
| url | https://nasplib.isofts.kiev.ua/handle/123456789/30409 |
| work_keys_str_mv | AT níkolaênkotû kvantovomehaníčniikonformacíiniianalízmolekuli2dezoksicitidilovoíkislotistrukturnoílankidnk AT govorundm kvantovomehaníčniikonformacíiniianalízmolekuli2dezoksicitidilovoíkislotistrukturnoílankidnk AT níkolaênkotû quantummechanicalconformationalanalysisof2deoxycytidilicacidmoleculethednastructuralunit AT govorundm quantummechanicalconformationalanalysisof2deoxycytidilicacidmoleculethednastructuralunit |