Бесплатный автореферат и диссертация по биологии на тему
Внутримолекулярное плавление ДНК
ВАК РФ 03.00.03, Молекулярная биология

Автореферат диссертации по теме "Внутримолекулярное плавление ДНК"

АИ 00

МИНИСТЕРСТВО МЕДИЦИНСКОЙ И МИКРОБИОЛОГИЧЕСКОЙ ПРОМЫШЛЕННОСТИ СССР ВСЕСОЮЗНЫЙ НАУЧНО-ИССЛЕДОВАТЕЛЬСКИЙ ИНСТИТУТ ГЕНЕТИКИ И СЕЛЕКЦИИ ПРОМЫШЛЕННЫХ МИКРООРГАНИЗМОВ

На правах рукописи

УДК 547.963.3 577.323.4

ЛЮБЧЕНКО ЮРИИ ЛЬВОВИЧ

ВНУТРИМОЛЕКУЛЯРНОЕ ПЛАВЛЕНИЕ ДНК

(Специальность 03.00.03 - Молёкулярная биология)

Автореферат диссертации на соискание ученой степени доктора биологических наук

Москва - 1988

Работа выполнена в Институте молекулярной генетики Академии Наук СССР

Официальные оппоненты:

доктор биологических наук, профессор И.Н.Тодоров доктор физико-математических наук, профессор Ю.П.Влагой доктор биологических наук Ю.И.Козлов

Ведущая организация:

ЫежФакультетская проблемная научно-исследовательская лаборатория молекулярной биологии и биоорганической химии им. А.Н.Белозерского при МГУ им. М.В.Ломоносова

Защита состоится " " 1988 г. в часов

на заседании Специализированного ученого совета Д.025.01.01 при Всесоюзном Научно-исследовательском институте генетики и селекции промышленных микроорганизмов по адресу: 113545, Москва 1-й, Дорожный пр., I.

С диссертацией мокно ознакомиться в библиотеке Инетитута1 ВНИИ Генетика.

Автореферат разослан " н 1988 г.

Ученый секретарь Специализированного совета кандидат биологических наук

В.И.Щербакова

ВВЕДЕНИЕ

Актх§льность_проблемы

Плавление - одно из наиболее важных явлений, наблюдаемых в молекулах ДНК, и интерес к нему двоякий. С одной стороны, как Физическое явление, внутримолекулярное плавление ДНК занимает значительное место в молекулярной физике вообще, и в Физике полимеров в частности. С другой стороны, детальное понимание процесса плавления важно для изучения функцишаль'ных свойств ДНК, поскольку в клетке эта молекула не инертна. В процессе акционирования в ДНК происходят конФормационные превращения, возникают новые структуры, и, как правило, все такие'перестройки идут с разрушением исходной спиральной структуры молекулы.

Понимание всех особенностей плавления ДНК важно и при использовании явления внутримолекулярного плавления ДНК как метода, позволяющего анализировать особенности в строении ДНК.

Изучение плавления ДНК началось вскоре после того, как Уотсо-ном и Криком была установлена ее-спиральная структура. Экспериментальные исследования проводились с применением различных Физических методов. Было показано, что при повышении температуры спиральная структура разрушается, и происходит такой процесс в интервале температур порядка нескольких градусов.

Вслед за первыми экспериментальными данными появились теоретические работы, а спустя 10 лет была сформулировайа теоретическая модель, удовлетворяющая накопленным экспериментальным данным.

На основании теоретической модели плавления можно предсказывать юведение ДНК при температурах, лежащих ниже интервала плавления в условиях, близких к физиологическим, выяснять влияние как внешних факторов, так и свойств самой молекулы на структурные перестройки в ней. Понятно, что такое широкое использование теории возможно лишь после ее всесторонней проверки.

Здесь стоит отметить, что в основном адекватность теоретической модели проверялась лишь по одному параметру, характеризующему переход спираль-клубок в.ДНК,' а именно, по ширине интервала перехода. Конечно, теоретические'расчеты позволяли получать и более детальные характеристики плавления ДНК, но отсутствие соответствующего экспериментального материала не позволяло провести необходимое сопоставление.

Ситуация коренным образом изменилась с середины 70-х годов, когда выяснилось, что кривая плавления ДНК не гладкая, как считалось ранее, а содержит ряд ступенек и перегибов. При дифференцировании таких крйвых по температуре выявляются узкие пики шириной в несколько десятых долей градуса. По аналогии со спектроскопией этот эффект получил название тонкой структуры кривых плавления ДНК и, конечно, он сразу же привлек внимание исследователей к проблеме плавления ДНК.

В исследовании ДНК наступил качественно новый этап, названный известным японским ученым А.Вадой "эрой тонкой структуры кривых плавления". За небольшой промежуток времени было опубликовано более 100 работ, в которых исследовалось плавлений с высоким разрешением ДНК из различных источников, от простейших прокариот до иукариот. В таком виде экспериментальные данные содержат значительно больше информации.

Возникающие при этом вопросы можно условно разбить на три группы:

1. Получение экспериментально обоснованного объяснения эффекта тонкой структуры кривых плавления ДНК.

V.. Всестороннее изучение явления внутримолекулярного плавления ДНК I: целью получения его полной Физической картины, включающей помимо термодинамики еще и основные кинетические характеристики процесса плавления ДНК, полная проверка адекватности теоретичес- ' ких моделей плавления ДНК.

'Л. Выяснение новых возможностей плавления ДНК как метода исле-

дяьании ее строения.

11<_Ч]ь. Уё1'ТоаШй0_Вйботы состояла в том, чтобы с помощью метода плавления ДНК с высоким разрешением решить сформулированные задачи.

На эчщиту выносятся следующие основные результаты работы:

[. Проведено экспериментальное исследование плавления ДНК из разных источников, на основании которого получено эксперимен-

- и -

тальное обоснование об'яснения природы эффекта тонкой структуры кривых плавления ДНК.

2. Разработан принципиально новый подход к изучению плавления ДНК, позволивший впервые получить прямые- данные о характере плавления ДНК при любой температуре в интервале плавления.

3. Получены полные данные о плавлении ДНК плазмиды Со1Е1, связавшие карты плавления с пиками профиля плавления. На основании этих данных получено экспериментальное обоснование модели возникновения необратимости плавления ДНК. Предложен метод прямого определения Фактора кооперативное™ плавления ДНК и определено с высокой точностью его значение. '

4. Проведена прямая проверка теории перехода спираль-клубок в ДНК* ¡путем сопоставления с результатами теории профилей и карт плавления ДНК с известной нуклео.тидной последовательностью.

5. Исследованы неравновесные стадии плавления ДНК и выяснено их влияние на профиль плавления ДНК. Проведен детальный анализ особенностей протекания процессов денатурации и ренатурации отдельных доменов ДНК.

6. Исследована возможность использования плавления ДНК с высоким разрешением как для выяснения ряда свойств ДНК, так и для анализа особенностей в ее первичной структуре. Выявлены и локализованы АТ-богатые участки в ДНК Фага Т7. Показана высокая чувствительность профиля плавления ДНК к локальным изменениям в ее нук-леотидной последовательности.

Перечисленные вопросы охватывают практически полностью всю проблему внутримолекулярного плавления ДНК, и это дает основание говорить о построении завершенной картины этого Физического явления.

Решение поставленных задач потребовало не только кардинального совершенствования и автоматизации самого метода плавления ДНК, основанного на широком использовании вычислительной техники, но

я

и создания новых методов и подходов, позволивших получить полную информацию о плавлении ДНК.

Науцная_новиэна_и_ыуактическая_значимости

Все приведенные в работе данные получены автором впервые.

Прежде всего это относится к экспериментальному обоснованию причини возникновения самого эффекта тонкой структуры кривых плавления, полученному на основании сопоставления данных по плавлению ДНК из различных источников с расчетами.

Разработан принципиально новый экспериментальный подход к изучению плавления ДН|$- получение карт плавления, основанный на химической модификации ДНК глиоксалем в интервале плавления ДНК и последущим электронномикроскопическим анализом частично денатурированных молекул. ,Благодаря возможности проведения такой процедуры для любой стадии в интервале плавления ДНК удалось впервые проследить за ходом плавления ДНК от начала до конца и сопоставить данные по локализации на молекуле ДНК расплавленных участков (топография плавления) с характером их плавления, получаемым из профиля плавления ДНК.

Наличие такой полной информации о плавлении ДНК плазмиды СоХЕ1 позволило разработать метод прямого определения одного из важных параметров теории перехода спираль-клубок в ДНК -фактора коопе-ративности плавления б. Прямое сопоставление экспериментальных данных с расчетными дало возможность впервые измерить значение параметра о для ДНК с высокой точностью.

Сопоставление карт плавления с данными по денатурации и рена-турации ДНК Со1Ё1 позволило идентифицировать стадии неравновесного плавления ДНК и обосновать ранее предложенную модель возникновения необратимости плавления ДНК.

Наличие полного экспериментального материала-карт и профилей плавления, полученных для ДНК с известной нуклеотидной последовательностью, позволило впервые провести детальную проверку теории. На основании этих данных выяснены возможности основной теоретической модели перехода спираль-клубок в ДНК и установлена ограниченность применимости модели двух соотояний.

Детально иссследован один из наиболее важных вопросов плавления ДНК-неравновесность плавления. Впервые экспериментально доказана з&в< симость проФиля плавления ДНК от скорости нагревания и пока-

зано, что эффект обусловлен неравновесно плавящимися участками ДНК. Сопоставление данных по плавлению и ренатурации ДНК с известной топографией плавления позволило получить количественные оценни времен релаксации отдельно плавящихся участков, выяснить принципиальные различия в температурной зависимости двух путей ренатурации участков-ренатурации с концов и нуклеации внутри.

С применением эффекта тонкой структуры впервые удалось выявить, а затем, используя специальный подход, локализовать в ДНК Т7 АТ-богатые участки.

Исследование плавления ДНК с локальным изменением нуклеотидной последовательности выявило необычайно высокую чувствительность профиля плавления ДНК к нуклеотидным заменам, позволяющую регистрировать делеции даже одной'нуклеотидной ^ары.

Полученные в работе данные, а также разработанные методы и подходы могут оказаться полезными в ряде приложений, например, для изучения строения рекомбинантных ДНК, выявления и локализации на ДНК схожих и, наоборот, отличающихся участков с их детальной характеристикой, исследования взаимодействия ДНК с белками и низкомолекулярными лигандами.

Проведенные исследования охватывают практически все вопросы, относящиеся к равновесной термодинамике плавления ДНК, и,в основном, завершают исследования по этой проблеме, значительно улучшают наше понимание процесса протекания отдельных стадий плавления ДНК и открывают новые перспективы использования плавления ДНК как метода исследования.

Структ^Ба_работы

Диссертация состоит из введения и 6 частей. Во введении обоснована актуальность темы исследования, сформулирована цель работы, ее новизна и практическая значимость. Первая часть представляет собой обзор литературных данных, относящихся к исследованиям плавления ДНК до обнаружения эффекта тонкой Структуры.В ней содержатся основные экспериментальные данные о плавлении ДНК, а также теоретические модели. Вторая часть - методическая, и в ней, наряду с подробным описанием применявшихся методов и приемов, даны характеристики препаратов ДНК и ферментов. В автореферате изложено содержание последующих частей диссертации, в которых отражены результаты исследований. Каждой части предшествует

короткий обзор данных ПО данной проблеме. В конце диссертации после заключения приводятся выводы.

Апро^ауия^аооты/

Результаты работ донладын.инюь на Втором всесоюзном биохимическом с'езде (г.Ташкент, ¡Уб9г\), на П Всесоюзной конференции по конформационным изменениям опопшшморов в растворах (г.Тбилиси, 1975г.), на X Всесоюзной конференции по электронной микроскопии (Таллин, 1979г.), на У Всесоюзном совещании "Плазмида" (Г.Тарту, 1979г.), на Симпозиуме <№Б0 по ДНК (г.Прага, 1979г.) на У Всесоюзном совещании но конфирмационным изменениям биополимеров в растворе (г.Тбилиси. 1980г.), на конференции СЭВ "Биофизика нуклеиновых кислот и белков" (г.Райнхайнбрун, 1980г.),' на Симпозиуме СЭВ но бчтЬи:тш) нуклеиновых кислот и нуклеопротеидов (г.Таллин, 1981г.), на Всесоюзном биофизическом съезде (г.Москва, 1982г.), на ХП Всесоюзной конференции по электронной микроскопии (г.Сумы, 1982г.), на Симпозиуме СЭВ "Физико-химические свойства биополимеров в растворе и клетках" (г.Пущино, 1985г.), нн Всесоюзном совещании по конформационным изменениям биополимеров н растворах ¡Тбилиси, 1985г.), на Международном симпозиуме "Физико-иимия ДНК и молекулярные механизмы функционирования генома (Тбиписи, 1987г.).

В этой части рассмотрены данные по экспериментальному изучению плавлиния ДИН, полученные разними методами до обнаружения эффекта гонкой структуры кривых плавления. Коротко рассмотрены теоретические модели плавления полинуклеотидов.

Приводится подробное описание использованных в работе методов. включая методику измерения температуры; методы получения дифференциальных кривых плавления как путем обработки исходных нриаих плаьленяя (графические дифференцирование и использование коипгЛ)теров), так и путем прямой регистрации дифференциальных крин: > плавления. Изложена подробно процедура модификации ДНК ь им:е^ьале плавления глиоксалем и последующий электронномикро-скоги''лч:ниЙ анмиз полученных препаратов. Отдельная глава посвя-

щнча харлктериляции препаратов ДНК и Фсрментои.

ЭфФект тонной структуры кривых плавления ДНК демонстрирует рис. I, на.котором приведена дифференциальная кривая плавления ДНК Фат. Т", полученная дифференцированием исходной кривой плавления по температуре. Видно, что кривая не гладкая, а содержит отчетливые узкие пики и ступеньки, и именно эти особенности плавлении ДНК получили название эффекта тонкой структуры.

Т °С

Рис.1 Профиль плавления ДНК фага Т7 в 0,1х55С. Символы • и о отвечают двум независимым экспериментам.

Впервые тонкая структура кривых плавления была обнаружена в работе Гомез и Ланга (1972), в которой была приведена дифференциальная кривая плавления ДНК Т7, но удовлетворительного об'яс-нения этому эффекту в работе не приводилось. Аналогичные особенности наблюдались также на кривых плавления ДНК из митохондрий дрожжей (Мишель и др.,1974), эти данные Также указывали на наличие интересного эффекта, который и был детально исследован в настоящей работе.

Вернемся вновь к рис.1, на котором приведена дифференциальная кривая плавления (в дальнейшем такие кривые будем называть профилями плавления) ДНК фага Т7 в 0.1х$£С На этой кривой вое особенности тонкой структуры перенумерованы. Обращает на себя вни-

мание наличие очень узкил пиков, так, например, пик 6 имеет ширину О,4°С.

Независимо от первых экспериментальных исследований аналогичный эффект бйп получен при чеоретических расчетах кривых плавле-

т-ъ

Рис.2 Профили плавления, рассчитанные для двунитевых полинуклео-тидов из 30 ООО случайно чередующихся пар АТ и вС. Сплошная линия и пунктир отвечают двум различным последовательностям.

Вологодским и Франк-Каменецким были проведены расчеты в рамках модели плавления двунитевого гетерополимера, учитывающей образование петель при плавлении ДНК (основная модель, Веденов и др., 1971). На рис.2 приведены профили плавления для двух гетерополи-меров со случайным чередованием нуклеотидов из 30 . ООО звеньев. По своему виду эти расчетные кривые похожи на экспериментальную кривую для ДНК Т7 (см.рис.1), на них также видны пики шириной 0,5°С. Помимо профилей плавления, расчеты дают и другую важную характеристику-топографию плавления. Сопоставление этих данных с пиками на профилях плавления позволило выдвинуть объяснение происхождения эффекта тонкой структры кривых плавления ДНК.

Оказалось, ,что в ДНК,из-за кооперативного характера ее плавления, идет выплавление участков из нескольких сотен пар нуклеотидов^ происходит это в интервале температур в десятые доли градуса, что и приводит к появлению узких пиков на профилях плавления.

Наглядное об'яснение может быть получено, если обратиться к упрощенному представлению о плавлении ДНК, основанному на том, что каждый участок ДНК может находиться либо в полностью спиральном, либо в полностью расплавленном состоянии. Тогда его плавление представляет собой переход по принципу "все или ничего", и в таком приближении для участка из п пар нуклеотидов получается следующее выражение для ширины интервала плавления:

г2

АТ= 4 Цо. П&Н

(I)

где Т0 - температура плавления участка (температура пика на профиле плавления) идН- энтальпия плавления в расчете на одно звено (энтальпия пары нуклеотидов для плавления данного участка с учетом его среднего СС-содержания).

Из формулы (I) следует, что для значениЙ4Н=9ккал/моль и Т„= 350°К для участка из 400 п.н. величина дТ=0<3°С. Следовательно^ денатурации такого участка должен отвечать узкий пик на профиле плавления. Интеграл под этим пиком равен отношению длины участка (« ) к длине всей ДНК, т.е. для ДНК изА^-звеньев эта величина равна п/Ы . Для Л^=40 ООО п.н. пШ=1%, т.е. эффект невелик, и он тем меньше, чем длиннее ДНК. К тому же,с увеличением длины ДНК возрастает число кооперативно плавящихся участков, что при той же ширине интервала плавления всей ДНК приведет к наложению пиков, т.е. сглаживанию всей кривой.

Итак, из расчетов следует, что на профилях плавления ДНК должны появляться узкие пики (-~0,5°С)-тонкая структура кривых плавления. Эффект должен проявляться тем ярче, чем короче ДНК.

Для проверки • этих выводов были проведены исследования на ряде ДНК различной длины. На рис.3 приведены профили плавления трех ДНК длиной от 104до 17.104п.н. Кривые А и В обладают тонкой структурой, и наиболее ярко она выражена для самой короткой молекулы- ДНК фага РМ2 {10 600 п.н.). На профиле плавления

этой ДНК имеется отчетливый узкий пик в начале (Т=68°С). Он имеет ширину 0,4°С, и интеграл под ним отвечает денатурации участка из £>00 п.н,

Видно, что увеличение длины генома (ср. рис.1 и рис.За) приводит к заметному сглаживанию профиля плавления ДНК. Еще слабее тонкая структура проявляется доя ДНК фага Сд (72 ООО п.н.),и она не регистрируется совсем на профиле плавления ДНК Фага Т2 (170 ООО п.н.).

Таким образом, экспериментальные данные показывают, что профили плавления ДНК геномов размером не более 100 ОООп.н. обладают тонкой структурой. характеризующейся наличием узких пиков. Тонкая структура тем ярче выражена, чем короче геном.

Все это находится в соответствии с тем, что предсказывает теория, тем самым подтверждается предложенное об'яснение природы эффекта тонкой структуры.

т«с

jtyc.3 Профили плавления в 0,Ix SSC ДНК из различных источников: Aj) ДНК фага РМ2: В) ДНК фага Сд и С) ДНК фага Т2

Предлагались так же другие об»яснения, но они оказались ошибочными. Так, в частности,в работах Вады и др.(1976) наличие узких ' гужов на. профилях плавления ДНК связывалось с денатурацией протяженных участков с регулярной нуклеотидной последовательностью. Этот вывод основывался на результатах расчетов профилей плавления, проведенных в рамках модели, учитывающей влияние на стабильность пары нуклеотидов соседней (т.н. эффект гетерогенности стэкинг-взаимодействий). Безусловно,учет этого эффекта является .правильным шагом на пути улучшения теории, однако Вологодским и Франк-Каменецким (1977) было показано, что в работах Вады и др. именнЛ'эдесь была допущена ошибка, которая привела к неправильной трактовке природы эффекта тонкой структуры кривых плавления.

Безусловно,сопоставление'экспериментальных данных, полученных на природных ДНК,с результатами расчетов, проведенных на модельных последовательностях, чаще со случайным чередованием звеньев, нельзя считать вполне удовлетворительным. Развитие методов сек-венирования ДНК дало возможность начать исследования ДНК с известной нуклеотидной последовательностью. Первой была расшифрована нуклеотидная последовательность генома фага ¡Ш74, и на репликативной двунитевой форме этой ДНК впервые была проведена работа по сопоставлению данных теории и эксперимента.

На рис.4 приведены экспериментальный (а) и теоретический (б) профили плавления двунитевой РФП формы ДНК 0X174 в 0,1х5£С Кривая (б) быларассчиуана Вологодским и Франк-Каменецким в рамках основной теоретической модели с использованием строгого метода расчета, предложенного в работе Фиксмана и др.(1977).

На рис.4 можно выделить две области. Ниже 72° С экспериментальная и теоретические кривые похожи друг на друга по положению пиков и относительной амплитуде, хотя и отличаются по абсолютной амплитуде. Наличие теоретических карт плавления дает возможность проанализировать полученные данные. Укажем на некоторые моменты. ТвЯ, основной пик при 71,6°С отвечает плавлению участка из-800 п,«.. Несколько менее интенсивный пик при 70,9°С связан с денатурацией двух соседних участков из~300 п.н. Эти данные соот-ветптнумт изложенный выше представлениям о плавлении ДНК.

гТ

Рис.4 Профили плавления двунитевой репликативной формы (РФП) ДНК фага 0X174 в 0,1x556 а) экспериментальная кривая, в и о отвечают двум независимым экспериментам: б) теоретическая кривая, рассчитанная по известной нуклеотидной последовательности.

Выше 72,7°С наблюдается сильное отличие экспериментальной кривой от теоретической. Позже аналогичные данные были получены в группе Вады (1978), и из них следовало, что наряду со случаями хорошего совпадения расчетных профилей плавления с теоретическими, имеются и сильные расхождения. На другое сильное отличие данных теории и эксперимента указывалось в работе Визарда и Айсеви-

на С1978)( в которой сопоставлялись с результатами расчетов экспериментальные данные по плавлению ДНК ФХ174 в различных ионных условиях. . Оказалось/ что'в отличие от наблюдающихся па эксперименте сильных изменений профилей плавления в зависимости от ионной силы,теоретические профили не меняются.

Таким образом, теоретическое и экспериментальное изучение плавления ДНК с известной нуклеотидной последовательностью, подтвердив в основном теоретическую модель, поставило ряд принципиальных вопросов. Этоотносится как к теории, так и к эксперименту. Так, например, в теории не все параметры известны с хорошей точностью. Кро^е того, ее можно совершенствовать, усложняя модель и вводя новые параметры. Что же касается эксперимента, то здесь возникает вопрос о равновесности наблюдаемого процесса

I

плавления ДНК.

Ясно, что для ответа на поставленные вопросы нужны дополнительные исследования, и решающую роль в них сыграл принципиально новый подход, позволивший получить такую важную информацию о плавлении ДНК,как карты плавления.

тенссмнтикс/*с

Рис.5 Кривая плавления (а) и дифференциальная кривая плавления Ш лилейной молекулы ДНК ллаэмиды Со1Е1 в 0,1х£5С. Стрелками обозначены те точки в интервале плавления, при которых проводились фиксация.

ПВЗМО§_Ц§блюдение_плавления_ДНКр:.

Этот подход был развит в настоящей работе и состоит в фиксации глиоксалем расплавленных участков ДНК с последующим электрон-номикроскопическим анализом полученных препаратов ДНК. Была разработана процедура фиксации, позволяющая осуществить модификацию оснований в уде расплавленных участках, предотвращая их ренату-рацию при охлаждении. Причем сама методика позволяла осуществить такую фиксацию при любой температуре в интервале плавления ДНК.

Для проведения картирования бьша'выбрана ДНК СоСЕ1, обладающая ярко выраженной -тонкой структурой кривой плавления. На рис.5 приведены интегральная (а) и дифференциальная (б) кривые плавления линейной формы ДНК Со1Е1 в 0,1x55^. Видно, что профиль плавления содержит много отчетливых пиков. Обращает на себя внимание первые два совершенно изолированных пика шириной 0,5°С. Для выяснения вопроса, каким участкам ДНК Со1Е1 отвечает появление пиков на профиле плавления, была осуществлена фиксация в 9 точках в интервале плавления. Температуры фиксации обозначены на рис.5 стрелками и видно, что они, как 'правило, отвечали минимумам на профилях плавления. Сразу же после окончания процедуры фиксации готовились препараты для электронной микроскопии.

Рис. 6 Микрофотографии молекул ДНК Со1Е1, отвечающих температурам фиксации 65°С (а) и 72,2°С (б), расплавленные участки обозначены стрелками,

На рис.6 приведен электронномикроскспический снимок для двух точек в интервале плавления. На них отчетливо видны расплавленные участки как на.концах молекулы, так и внутри (на Фотографиях расплавленные участки обозначены стрелками). Такие молекулы измеряли, длины и координаты расплавленных и спиральных участков вводили в калькулятор НР9825А (Хьюлетт-Паккард) и по массиву молекул для каждой из точек строили карты плавления.

Рис.7 Трехмерная диаграмма плавления ДНК Со1Е1. Положение плавящегося участка на молекуле ДНК отложено по горизонтальной оси, степень денатурации- по вертикальной.

На рис.7 приведены в виде трехмерной Диаграммы карты плавления, отвечающие 9 точкам фиксации. На этой диаграмме по горизонтальной оси отложена координата плавящегося участка на молекуле ДНК, по вертикальной- степень денатурации каждого участка, а по наклонной йгихгвмпература. Для наглядности приведен также профиль плавления ДпК СоШ. Из карт видно, что первому узкому пику от-

аечает денатурация участка из~400 п.н. на левом конце молекулы, следующему, также узкому пику, отвечает денатурация участка из -600 п.н. на правом конце. Затем плавится участок в спиральной части молекулы (координаты 74,1-79,1). Отметим также противоположную ситуацию: плавление участка между двумя расплавленными (например, участок из 360 п.н. с координатами 61,0-66,7). Одновременно с ним плавится участок приблизительно такой же длины слева от ранее расплавленного, в результате чего возникает в правой части молекулы петля размером -1500 п.н. .Наиболее тугоплавкая область ДНК РсЛЕ1 располагается в левой части молекулы,и она удерживает комплементарные нити от расхождения.

Таким образом, полученные данные наглядно и однозначно демонстрируют одно из важнейших свойств плавления ДНК - коопера-тивность: мы видим, что в молекуле происходит постадийное выплавление участков ДНК из нескольких сотен пар нуклеотидов. Сопоставление карт плавления с профилем плавления показывает, что каждый из таких участков плавится в интервале температур в десятые доли градуса. Тем самым, полученные данные являются прямым и надежным подтверждением раннее выдвинутых теоретических

представлений о Физике внутримолекулярного плавления ДНК.

Наличие таких полных данных о ходе плавления даже одной ДНК из плазмиды Со1Е1, позволяет решить ряд общих задач внутримолекулярного плавления ДНК.

Неравновесность плавления ДНК.

До сих пор, говоря о плавлении ДНК, мы молчаливо предполагали, что оно представляет собой термодинамически равновесный процесс. Слабая зависимость профиля плавления ДНК от условий проведения эксперимента казалось бы указывала на равновесность плавления (см:Вада и др.,1976). Однако в работе ХофФа и Роса (1972) были получены экспериментальные данные, свидетельствующие о существо-* вании необратимых стадий плавления ДНК, которые проявляются в наличии гистерезиса: несовпадении кривых частичного плавления и последующей ренатурации. ХоАФом и Росом была сформулирована модель, объясняющая появление гистерезиса. Согласно модели воз-

никновение стадий неравновесного плавления связано с плавлением участков ДНК; располагающихся между уже расплавленными, поскольку обратный процесс, представляющий собой ренатурацию протяженной области, в условиях низкой ионной силы может быть сильно замедли»1

Рис.8 Профили частичного плавления (сплошная линия) и последующей ренатурации (пунктир) ДНК Со1Е1 в 0,1х55С. Данные работы Вады и др.(1980).

Наглядно наличие гистерезиса регистрируется для ДНК, обладающей ярко выраженной тонкой структурой кривых плавления. На рис. 8 приведены данные по плавлению и ренатурации ДНК СоСEI в O.IxSSC (данные взяты из работы Вады и др. 1980). Сплошной линией обозначены профили плавления, а пунктиром - ренатурации. Видно, что вплоть до появления на профиле плавления пика при Т= 74°С,все стадии плавления ДНК Cof-EI обратимы: профили ренатурации и денатурации совпадают. Наличие гистерезиса на последнем графике указывает на То, что на стадии плавления ДНК CoCEI, отвечающей появлению пика при Т=74°С,плавление протекает неравновесно.

Полученные нами полные данные по плавлению ДНК CoCEI могут ответить на вопрос, в чем причина возникновения необратимости. Из рис. 7 видно, что вплоть до пика при температуре 74°С все предыдущие соответствуют денатурации участков либо на конце спирали, либо внутри спирали, и лишь при Т=74°С наблюдается плавление участка между двумя расплавленными. Сопоставление этих данных с результатами, Приведенными на рис.8,показывает, что необратимость плавления ДНК CoIEI, наблюдавшуюся в работах Вады и др. (1980), следует связать с денатураций участка ДНК, плавящегося между ранее расплавленными областями, что находится в полном соответствии с моделью Хоффа и Роса.

Дальнейшие исследования, проведенные Перельройзеным и др. (1982) по такой же схеме, когда данные по ренатурации и денатурации других ДНК сопоставлялись с картами их плавления, показали, что необратимость возникает и в тех случаях, когда плавится участок на конце молекулы рядом с расплавленной областью.

Однако эти данные не дают (Тгвета на вопрос, насколько сильно неравновесность искажает профиль плавления ДНК. К этому вопросу мы вернемся позжеf а здесь отметим, что, как показали дальнейшие исследования, неравновесность,несмотря на сильный гистерезисный эффектов малой степени искажает профиль плавления ДНК, и различия в профилях плавления ДНК, получающиеся при разной ионной силе, об' ясняются другими причинами.

0ОВ§аедение^актдЕа_кодп^

■ При анализе данных по плавлению ДНК ФХ 174 отмечалось, что не все теоретические параметры известны с достаточно высокой точностью, и,в частности, это относилось к такому важному параметру как (Ьактор кооперативности плавления ДНК. Анализ данных по плавлению ДНК СоГЕ1 позволил разработать подход по прямому измерению этого параметра. -Ъ/цт

Фактор кооперативности € , где ^ - свободная энергия, затрачиваемая на образование двух границ между спиральной и расплавленной областями, и для кооперативно плавящегося участка она может быть определена из разности температур одного и того же участка, возникающей из-за того, что в одном случае он находится внутри молекулы, а в_другом -на ее конце.

ТСИРСЙЯТиЯС/ "С __;

Рис.9 ,_Профад;Гплавления в "6,1x55(1 молекул ДНК СоШ, получающихся при расщеплении рестриктазами Есой! (А) и Г (В).

В качестве экспериментальной системы была выбрана ДНК Со1Е1, у которой, как видно из рис. 7, первый плавящийся участок из 380 п.н.' располагается на конце молекулы в том случае, когда кольцевая ДНК Со1Е1 расщепляется рестриктазой ЕсоИ1. Но при расщеплении рестриктазой 5/гм2этот же участок уже плавится внутри спирали ДНК.

На рис.9 А и В показаны профили плавления, полученные для таких препаратов ДНК. Видно, что основное отличие в этих кривых, связано с положением первого пика: его температура на~1°С ниже для ЕсоК - расщепленной ДНК по сравнению со случаем, когда линеаризация молекул проводилась рестриктазой $ггн!

В.5

В.Ч

>

• ИЗ ч

ни

I

0.1

е.в

5В 6В 62 64 66 ТЕИРЕЯЯТиЯЕ/ ШС

Рис. 10 Частичные профили плавления в 0,1х53Ссмеси линейных молекул ДНК СоШ обоих типов.

Наглядно это видно из рис.10, на котором приведены частичные профили плавления препарата ДНК, содержащего оба типа молекул.

Первые два пика соответствуют плавлению одного и того же участка ДНК, а второй очен1 интенсивный пик обусловлен денатурацией второго участка, котофг^ в обоих типах молекул ДНК плавится практически при одной и той же температуре.

Аналогичные данные были получены для различных ионных условий, и для 0.2 М Мч+(при этих ионных условиях плавление первого участка молекул ДНК, получающихся расщеплением рестриктазой , происходит равновесно) полученные экспериментальные данные сопоставлялись с расчетными. Наилучшему совпадению с экспериментом соответствуют данные, отвечающие значению (Т=(2,575)-Ю-".

Во время выполнения этой работы было обнаружено, что в линейных молекулах Со1Е1 плавление даже первого участка в том случае, когда он располагается внутри спирали Ота -расщепление) в низкой ионной силе,происходит равновесно. Видимо, обнаруженная необычная ситуация проявления Неравновесности обусловлена особенностями в нуклеотидной последовательности в этой части ДНК Со1Е1. В ней самый первый плавящийся участок йшанкирован двумя очень АТ-богатыми последовательностями, препятствующими ренатурации всей области с концов. Тем самым,здесь мы имеем ситуацию, подобную модели Хоффа и Роса.

Неравновесность плавления ДНК СоШ не позволила определить б для других значений ионной силы. Это было сделано в работе Ко-зявкина и др. (1987), в которой определение проводилось по такой же схеме* Но на другой ДНК. Для 0,2 М значение совпало с величиной ¡определенной нами, уменьшение же ионной силы приводило к падению о , так,что при десятикратном снижении,»уменьшилась в 30 раз. Это обстоятельство необходимо учитывать при сопоставлении экспериментальных данных с расчетами при различной ионной ейле. В. работе Гото (1983) были проведены расчеты с разной величиной . О, и значительным варьированием этого параметра удалось уменьшить расхождение между расчетными и экспериментальными профилями плавления ДНК ФХ17^ отвечающими различной ионной-силе.

ЧастьДУ л„ДНК_с_известкой_нунлеотидноЯ _последова'£ел ьностью^

Наличие полной информации о плавлении ДНК, получаемой из детальных профилей плавления и соответствующих им карт плавления, позволило провести всестороннюю проверку теоретических моделей плавления ДНК.

Ранее мы приводили данные пс сопоставлению экспериментального профиля плавления ДНК ФХ174 с расчетным. Расчеты проводились ь рамках строгого алгоритма Поланда-Фиксмана-Фрейера (1977), позволяющего точно и быстро проводить расчеты даже для длинных последовательностей. При этом отмечалось общее сходство расчетных профилей плавления с экспериментальны! и, однако полного совпаде-

Рис.П Экспериментальный (А) и теоретический (В) профили плавления ДНК рВ1?322 в 0,1х£$С. Стрелками обозначены точки, для которых получались карты плавления.

Аналогичные данные были нами получены и для двух плаэмидных ДНК - рВ!?322 и рАОЗ. Обе плаэмидн имеют ярко выраженную тонкую структуру кривых плавления. На рис.11 показаны результаты расчетов и эксперимента для ДНК рВВ322. Зидно, что з очличие о"1 пксл>-

риментапьной кривой, обладающей отчетливой тонкой структурой, расчетные кривые флее гладкие. Такой же результат был получен нами на ДНК плаэмиды рАОЗ, а также в работе Вады и др. (1977). Таким образом,получается, что расчеты , проведенные в рамках основной модели, не позволяют получить кривые плавления, совпадающие с экспериментальными.

Другой характеристикой плавления ДНК, которую можно получить, являются карты плавления, и результаты анализа карт плавления приведены ниже.

' МИР 1= □ 5 I -Г I □ N

Рис.12 Экспериментальные (сплошная линия) и расчетные (точки) карты плавления ДНК рВ1?322.

Экспериментальные карты получались также, как это делалось для ДНК Со1Е1, и они строились для температур за отдельным пиком на профиле плавления ДНК. Теоретические карты рассчитывались для той же степени денатурации. На рис.12 приведены карты плавления ДНК плазмиды рВ1?322. Экспериментальные карты получены для 4 точек в интервале плавления, обозначенных на рис.11 стрелками. На рис.12 они показаны сплошной линией. На этом же рис.12 приведены теоретические карты плавления, они обозначены точками. Мы видим хорошее совпадение теории с экспериментом: совпадают как размеры расплавленных и спиральных областей при всех степенях денатурации, так и их размер. Отсутствие на экспериментальных картах 1,2 и 3 расплавленного участка на левом конце молекулы объясняется трудностями регистрации с помощью электронной микроскопии расплавленных участков размером менее 100 п.н., т.е. 2% длины ДНК рВ1?322 или 2 ед. карты на рис.12. Именно такого размера участок можно ожидать, исходя из теоретических данных.

Аналогичные данные были получены для других ДНК (плазмида рАОЗ и фаг ФХ174),и во всех случаях было хорошее совпадение карт плавления экспериментальных с расчетными. Таким образом,строгие расчеты , выполняемые в рамках основной модели, дают вполне удовлетворительное совпадение с экспериментатьными картами плавления ДНК, но соответствие профилей плавления намного хуже.

Нами был также проверен другой теоретический подход, именуе мый моделью двух состояний. Согласно этому подходу предполагается, что плавление каждого участка ДНК происходит по принципу "все или ничего", т.е. участок*-ДНН может находиться лишь в двух состояниях: полностью спиральном или полностью расплавленном. Кривая плавления для каждого участка описывается простой аналитической функцией. Полная кривая плавления ДНК получается суперпозицией таких функций. Простота теоретического описания привлекла к модели двух состояний внимание значительного числа исследователей. В большинстве случаев оно диктовалось естественным желанием разложить полный профиль плавления на отдельные пики, но в отличие от первых работ Ансевина и др. (1976) и Габарро и др. (1979),данное разбиение основывалось на физической модели плавления ДНК.

Теоретические основы модели двух состояний в применении к анализу плавления ДНК были детально разработаны Азбелем (19791981). Разбиение ДНК на отдельные участки осуществлялось путем разложения исходного профиля плавления на отдельные пики и подбиралось такое разбиение, которое соответствовало данной нуклео-тидной последовательности. Азбелем (1980) было детально исследовано плавление РФП ДНК ФХ174 , для которой, на основании нашей экспериментальной кривой, была получена полная топография плавления. Сопоставление полученных результатов Азбеля с экспериментальными картами плавления представляет собой единственный путь' проверки выбранной теоретической модели.

-НИ а .

1жЩ

юн

~1Г

ЗМИИ

игапд - .

в

• ■Г Г«III I I1

С I

Рис.13 Карты плавления двунитевой ДНК /5X174: А) эксперименталь- . ные данные; В) карты< построенные по данным Азбеля (1980) и (С)-карты, полученные строгими расчетами .

На рис.13 приведены экспериментальные карты плавления РФ ДНК ФХ174 (А), которые сравниваются с результатами приближенного (В), выполненного для этой ДНК Азбелем (1980), и строгого расчета (с). Видно, что несмотря на полное совпадение профилей плавления, карты плавления, полученные Азбелем, неправильно описывают ход плавления этой ДНК. В особенности это относится к началу и концу плавления * хотя следовало ожидать полного совпадения. Наоборот, строгие расчеты, несмотря на значительное отличие в профилях плавления,дают практически правильную топографию плавления этой ДНК. Поэтому можно сказать; что модель двух состояний для детального описания плавления ДНК не годится. Последующий анализ показал, что причина кроется именно в ошибочности самого исходного предположения о плавлении участков ДНК по принципу "все или ничего": имеется очень большое число промежуточных состояний которые существенно влияют на характер плавления ДНК. Таким образом, прямое сопоставление экспериментальных и расчетных данных по плавлению ДНК с известной нуклеотидной последовательностью приводит к следующему. Расчеты, проводимые в рамках основной теоретической модели, в целом правильно описывают ход плавления. Однако эти же расчеты не дают удовлетворительного соответствия в профилях плавления. Ясно8 что необходимо совершенствование теории.

Ранее нами (см.гл.III) высказывалось предположение о необходимости учета в теории влияния на стабильность данной пары соседней, т.н. гетерогенность стэкинг-взаимодействия. Отсутствие необходимых параметров не позволяло учесть этот эффект, но если наблюдаемое различие в расчетных и экспериментальных профилях плавления связано именно с необходимостью учета гетерогенности стэкинг-взаимодействияí то соответствующие параметры можно определить путем подгойки расчетных данных под экспериментальные.

Впервые это было сделано в работах Гото и др". (1981), а затем Вологодским и др.(1983). Было показано, что можно подобрать необходимые параметры, характеризующие гетерогенность стэкинг-взаимодействия ¡'которые приближают теоретические кривые к экспериментальным, однако таким способом полного совпадемля теории с экспериментом добиться не удается.

Таким образом, этап исследования плавления ДНК с применением эффекта тонкой структуры оказался чрезвычайно плодотворным. В результате проведейных исследований были получены ответы практически на все основные вопросы, относящиеся к термодинамике перехода спираль-клубок в ДНК. Однако оказалось, что полные экспериментальные данные дают возможность изучить особенности кинетики протекания отдельных стадий плавления ДНК.

Часть_УА_Кинетические^

Ранее, в части-Ш, мы указывали, что в плавлении ДНК может возникать неравновесность, приводящая .к сильным гистерезисным эффектам. Били получены данные, подтверждающие модель возникновения необратимости, однако модель не объясняет как суть неравновесности, так и связанный с этим важный вопрос о том, насколько неравновесность искажает профиль плавления ДНК. Так, из общих соображений ясно, что в тех случаях, когда наблюдается гистерезис, профиль плавления ДНК должен меняться при изменении скорости увеличения температуры. Однако на опыте такого эффекта не наблюдалось. Более того, в работе Вады .и др. (1976) утверждалось, что зависимости от скорости нагревания нет,, и что плавление - равновесный процесс, а неравновесна обратная стадия - ре-натурация. Видно, что ситуация противоречивая, и прояснилась она лишь после работ Перельройзена и др. (1983) и Аншелевича и др. (1985), в которых были детально проанализированы особенности плавления участка со свободными концами. Оказалось, что денатурация такого участка,из-за медленности процесса денатурации при равновесной температуре плавления .происходит при более высокой температуре, и она тем выше, чем быстрее скорость изменения температуры. Однако,из-за очень сильной зависимости константы скорости денатурации от температуры, эффект сдвига наблюдаемой температуры плавлении очень мал: даже при десятикратном изменении температуры он не превышает десятых долей градуса.

Итак, сформулировано теоретическое представление, объясняющее аффекты неравновесности при плавлении ДНК, и одним из важных вы-ьодоа теории медленных релаксационных процессов является пред-

сказание слабой зависимости профиля плавления ДНК от скорости нагревания, которое и было проверено на эксперименте.

3§§ЙСЙМ25ТЬ_пвофиля_плавления_ЙНК_от

Регистрация столь слабого эффекта предъявляла строгие требования к эксперименту как в отношении спектрофотометрической точности, так и в особенности к измерению температуры. В наших опытах она составляла 0,02°С. Точность регистрации профилей плавления демонстрирует рис.14^ на котором приведены данные двух независимых экспериментов для ДНК Со1Е1. Видно * что кривые практически совпадаютj причем отметим, что наложение кривых проводилось без какого-либо относительного сдвига} их по температурной шкале и без трансформации по оси ординат. Кривые просто совмещались по оси температур. .

Рис.14 Профили плавления ДНК СоШ в 1х5$С, полеченные при скорости нагревания 0,05 град/мин. Приведены результаты двух независимых экспериментов,-

77

79

в:

аз

аэ

Н£ В 7 ВЯ 91 ТЕМРЕЯЯТиЯЕ/»С

Рис. 1Ь Профили плавления рестрикционного фрагмента С1 ДНК фага

Т7, пш!ученные при скоростях нагревания 0,05 град/мин (сплошная кринйл) и 0,5 град/мин (точки). А) 0,1х35С(на вставке приведена часть кривой ь увеличением масштабе) и В) I . -s.se

На рис.15 Л приведены профили плавления рестрикционного фрагмента 01 ДНК фага Т7 в 0,1х5$С, полученные при различных значениях скопости нагревания. Видно, что изменение скорости нагревания не меняет положение первых двух пиков, а положение следующего, третьего пика, меняется на 0,3°С при десятикратном изменении скорости нагревания. Полученные нами данные по электронномикрос-копическому картированию отдельных стадий плавления этой ДНК (на рис.16 изображены схематически такие карты) показали, что именно пик 3 отвечает денатурации неравновесно плавящегося участка, и после пика 3 наблюдается гистерезис: профиль ренатурации отличается от профиля денатурации.

(1) с

(2) С

13)1

I

15 ОЬр

250 Ьр

У////

500Ьр

I_I

500 1000

1*пд1Ь I» Ьи» р«1г»

1500

Рис.16 Схема*иллюстрирующая топографию плавления фрагмента ДНК С1. Заштрихованы расплавленные участки.

На рис.15В приведены аналогичные данные для высокой ионной силы (1х5$С). Здесь также отчетливо видны первые три пика, но теперь уже третий пик не меняет своего положения на профиле плавления в зависимости от скорости нагревания. Этот результат отвечает ожидаемому на основании полученных данных по плавлении и ренатурации.

В то *е самое время плавление последнего участка ДНК, приводящее к расхождению комплементарных1 нитей, независимо от ионных условий,из-за малой скорости реассоциации разошедшихся цепей всегда должно происходить неравновесно. Это наблюдается и на опыте. На рис.17 приведены данные по плавлению с разной скоростью в IxSSC рестрикционного фрагмента ДНК CoIEI. Видно, что отличие в положении последнего пика, отвечающего денатурации участка из 180 п.н., составляет 0,3°С.

Полученные экспериментальные данные собраны в табл.1, в которой они сопоставлены с результатами> ожидаемымй на основании теории мздленных релаксационных процессов.

Таблица I. Сдвиг в температурах плавления неравновесно плавящихся участков при 10-кратном изменении скорости нагревания.

ДНК (пик) O.IxSSC 1х S5C

&Т эксп Ьтк ът и1эксп &Тк

CI-фрагмент ДНК Т7, пик 3 0,32^.0,05 0,33±.0,03 0±0,05 0

EcoRI-фрагмент ДНКСое* I (1300 п.н.) пик4 0,24±0,05 0,42*0,04 0,39±0,05 0,42±0,04

бТ„„,~раэиица в температурах пиков при разной скорости нагрева-

оКСП

ния, 01^- ожидаемое различие, которое определяется исходя из со-отновения БТК=80^, где Л/- размер плавящегося участка, который определяется по относительному гиперхромному эффекту, отвечающему данному пику на профиле плавления. Ошибка определенияЙТ^ос-т&вляат а5Т -ЮЯ, и она связана с ошибкой в определении М .

Данные табл.1 демонстрируют вполне удовлетворительное согласие теории с экспериментом для Фрагмента С1. Несколько хуже совпадение для ДНК СоСЕГ, но и здесь значения довольно близки. Таким образом, полученные данные однозначно демонстрируют наличие зависимости профиля плавления от скорости нагревания, но эффект очень мал, и поэтому его не удавалось обнаружить раньше. Удовлетворительное количественное совпадение экпериментальных данных с теорией медленных релаксационных процессов свидетельствует о правильности теоретических представлений.

тсмрекипше:/ «с.

Рис.17 Профиль плавления рестрикционного фрагмента ДНК Со1В1 (1291 п.н.) в 1х эбс при скоростях нагревания 0,05 град/мин (.....) и 0,5 град/чин (-).

Нами, с целью получения количественных характеристик отдельных стадий денатурации и ренатурации ДНК, было подробно исследовано плавление и ренатурация фрагмента С1 при различных ионных условиях. На рис.18 приведены серии экспериментальных данных по денатурации и ренатурации ДНК С1 со скоростями 0,5 и 0^05 град/мин. Ренатурацию проводили от точки за третьим пиком, который отвечает неравновесному плавлению. Эти данные показывают влияние условий эксперимента на равновесность плавления. Так, в 0,4х взс плавление протекает равновесно лишь при очень малой скорости нагревания, но уже • уменьшение ионной силы в 2 раза настолько замедляет процесс ренатурации, что даже при скорости нагревания 0,05 град/мин плавление становится сильно неравновесным.

0.25

13.20

т

^0.15

ч

е и

\

ш0. 10 и

I

0.05

0.00

73 75 77 73 В1 ТЕМРЕКНТиЯЕ/ <С

Рис]8. Профили ренатурации и денатурации фрагмента С1 в 0,4хз5£г 1А): 0,1к ЭБС (В) и 0,05х-БЗС(С) при скоростях изменения температуры: +0,05 град/мин(-), +0,5 Град/мин(- - -), -0,5 град/мин

и -0,05 град/мин(.....К

Наличие карт плавления фрагмента С1 ДНК Т7 позволило провести детальный анализ процессов ренатурации всех трех участков при различной температуре. Анализ полученных данных показал , что два различных процесса ренатурации рассматриваемого участка ".зашивание" с концов и нуклеацию-можно разделить. И тогда получилось, что нуклеация, подобно реассоциации,характеризуется довольно слабой зависимостью от температуры. Наоборот, "зашивание" с концов, подобно процессу денатурации, представляет собой медленный процесс, сильно зависящей от температуры. Особенно наглядно это видно из рис.18С, на котором приведены данные для очень низкой ионной силы. В этих условиях ренатурация всей расплавленной области, состоящей из трех участков (см. рис.16), происходит целиком (ей отвечает узкий пик при Ь8°С) вслед за ренатурацией с правого конца участка I. Увеличение скорости охлаждения в 10 раз сдвигает пик на 0,6°С влево.

Приведенные данные позволили также получить оценки характерных времзн релаксации отдельных областей. Так, времена релаксации первого и второго участков фрагмента С1 не превосходят десятков секунд во всем исследованном диапазоне значений ионной силы. Для третьего участка ситуация иная: время его денатурации резко уменьшается с температурой. Независимо от ионных условий, время денатурации уменьшается в 10 раз при уменьшении температуры на 0,3°С. Это приводит к тому, что при равновесной температуре плавления время релаксации третьего участка уменьшается на три порядка (от секунд до десятков часов) при уменьшении концентрации ионов в 10 раз (от 1х ЭБСдо 0,1х&зс). Оказалось, что ренатурация путем нуклеации имеет совершенно иные характеристики: скорость нуклеации слабо уменьшается с понижением температуры, но резко растет с увеличением ионной силы. Таким образом, видно, что сопоставление результатов денатурации и ренатура-цим позволяет получить важные количественные закономерности в отношении плавления и ренатурации отдельных областей ДНК.

Из данных по плавлению с различной скоростью ДНК Со1Е1 следует, что даже в 1х зге у этой ДНК до полного расхождения нитей имеются стадии неравновесного плавления. Они соответствуют ситуации, когда после денатурации одного из участков образуется расплавленная область из~2000 п.н. Видимо здесь мы имеем ситуацию, похожую на расхождение нитей, когда поиск комплементарных участков идет медленно. Следовательно, необходимо для каждой из ДНК проверять равновесность плавления даже при высокой ионной силе.

Часть У1. Некоторые примеры использования эффекта тонкой струк-

В этой части диссертации рассмотрены два примера использования эффекта тонкой структуры для анализа особенностей нуклеотид-ной последовательности ДНК.

туры кривых плавления.

«л

о

1-1-I-I-1-Г-1-1-1-Г

66 68 70 72 74

тшшши ( * *0)

76 78

Рис.19 Профиль плавления ДНК фага Т7 в 0,1х БЗС

i.Siif-HiSiSSÜilS^iSl^^iSXyi-iiyäSIiiSS-S^iiiii-i'äiä-IZ •

На рис.19 приеден профиль плавления ДНК фага Т7 в 0,lxssc,

и на нем обращаем'на себя внимание самый первый пик (Т=68°С). Он отвечает плавление участков, денатурирующих заметно раньше остальной части ДНК Т7, суммарная длина которых составляет 400 п. н. (IX от всей ДНК). Из данных по плавлению ДНК Т7, полученных для длин волн 260 и 280 нм (при последней очень мал вклад в ги-перхромный эффект АТ-пар), удалось оценить среднее gc -содержание легкоплавких участков. Оно оказалось равным 28-2%, и эта величина заметно меньше среднего gc-содержания (48М.

Средний размер АТ-богатых участков можно оценить из ширины пика I. Для кооперативно плавящегося участка ширина кривой плавления определяется формулой (I). Если считать, что все АТ-бога-тые участки одинаковы и плавятся по принципу "все или ничего", то из формулы (I) следует8 что п=100-150 п.н., т.е. в молекуле ДНК Т7 должно быть 3-4 AT-богатых участка.

Для прямого определения числа таких участков в ДНК Т7, а также их локализации на ДНК, был использован следующий подход. АТ-богатые участки фиксировали в денатурированном состоянии глиок-салем, а затем по этим участкам ДНК гцфолизовали нуклеазой 51, атанующей денатурированные участки ДНК, и тем самым расщепляли ее на определенное число фрагментов.

На рис.20 приведены данные, полученные по такой схеме. Отчетливо видны три узких полосы, отвечающие размерам фрагментов 22j 12 и 5,6 т.п.н.. Причем, как видно из рис.20А и 20В, получающиеся фрагменты по своей гомогенности не отличались огрестрикцион-ных фрагментов ДНК. Это подтверждается также прямым измерением длин 51-фрагментов с помощью электронной микроскопии и аналитического центрифугирования.______ __ _________

Из данных по кинетике SI-расщепления ДНК Т7 по АТ-богатым участкам (см.рис.20С) следует, что фрагменты А, В и С расположены именно в таком порядке. Сопоставление данных по расщеплению ДНК Т7 дикого типа и делеционного мутанта DIU позволило ориентировать карту расположения АТ-богатых участков с генетической картой ДНК Т7. Из таких данных следовало, что "ранние" гены фага располагаются в наиболее коротком С-фрагменте.

\~А

1-е

-17 ОНА

Т7 ОМА

А* В. А

— В* С

В -С

Рис.20 Электрофоретическое разделение в 0,7% агарозном геле фрагментов, получающихся после расщепления ДНК Т7 по АТ-богатым местам (А), рестрикционных фрагментов ДНК фага Л (В) и фрагментов неполного расщепления ДНК фага Т? по АТ-учаоткам/С/,

Итак, предложенным способом в ДНК Т7 выявлено два АТ-богатых укастка. Из приведенных выше оценок следовало, что их должно быть 3-4, и в том^случае, когда АТ-богатые участки располагаются вблизи концов молекулы или рядом с обнаруженными, нужны более точные методы их выявления.

Для поиска дополнительных АТ-богатых участков были использованы рестрикционный анализ неразделенных ¿1- фрагментов и выделенных, гель-электрофорез рестрикционных фрагментов ДНК с фиксированными денатурированными участками и денатурационное картирование. Эти способы позволили выявить еще два АТ-богатых участка и результаты полученных данных представлены на рис.21, где также приведена карта наиболее важных в функциональном отношении областей генома Т7.

0 Ю 20 90 40 50 (0 70 во 90 ЮО

1 ' I ' \ ' I ' '1—' I ' I ■ I ' I '' I"

low-melting regions

и* г* т»* с* г* о» т>*

щ " |f п0 ■ ■ L Ц-—-IH Stand Saa 3fl

В GO E Я F С

i-u—i—i-u-f s»u 3fl

IGflTC)

flOftl promoter»

|leUltl I--сииш-,eenetlw(J

primary origin of rtpticalfon

О Юг0904050С0 70 80 90К» POSltiofl (T7units)

£ис. 21 Карта расположения АТ-богатых участков в ДНК фага Т7, которые показаны стрелками. Ниже приведены карты для рестрикта-зы ,'§аи ЗА и для нуклеазы 51 и положения наиболее важных "функциональных сайтов фага Т7,-

Видно, что АТ-богатые участки приходятся на функционально важные участки генома (промоторы АО, А1, участки начала репликации). Следует отметить, что данные по локализации АТ-богатых участков совпали с данными, полученными на основе нуклеотидной последовательности ДНК (Стадиер,1982), что указывает на адекватность предложенного подхода к выявлению легкоплавких участков ДНК на основе данных по плавлению.

Влияние локальных изменений в нуклеотидной последовательности ДНК на профили денатурации.

Наличие хорошо разрешенной тонкой структуры кривых плавления ДНК, когда отдельному пику отвечает плавление определенного участка, дает возможность обнаруживать .саке отдельные нуклеотид-ные замены в молекуле ДНК. На такую высокую чувствительность профиля плавления ДНК к нуклеотидной последовательности указывалось в нашей работе по плавлению ДНК <Ш74, а также в работах Азбеля (1979,1980). Замена АТ-пары на более стабильную сс-лару должна приводить к увеличению среднего ^с -содержания соответствующего участка ДНК, тем самым увеличивая его температуру, плавления. Замена бс-пары на АТ - наоборот уменьшает температуру плавления. Однако в обоих случаях эффект негелик и составляет не более 0,2°С.

III п , I •

А) 4 *

а_1—ь

... ААвбвТАТСбСббААТТСАТбГ^. ТС-7 а) ... ААвСбТАТСбСббААТТ.АТб... ОТ-19 ... ШМТМСЬСв......Щ... Щ-25

Рис.22 Рестрикционные карты ДНК плазмид ТС-7, ДТ-19 и ДТ-25 (А). Стрелкой указано положение сайтов, последовательности которых приведены ниже (В)

Исследования проводились на достаточно длинных ДНК (~П50 и.к.), обладающих ярко выраженной тонкой структурой кривых плавления. Все ДНК был(1 устроены таким образом, что нуклеотидные замены затрагивали СЬределенний участок ДНК, плавлению которого отвечал отдельный пик на дифференциальной кривой, положение которого и должно меняться при нуклеотидных заменах. Отметим, что такой подход, основанный на измерении сдвига одного из пиков по отношению к остальным, был ранее применен нами для определения фактора неоперативности плавления ДНК и изучения зависимости проФиля плавления ДНК от скорости нагревания. Во всех случаях измеряемый эффект составлял 0,2*0,5°С, т.е. был сравним с ожидаемым в случае нуклеотидных замен.

Для исследования были выбраны фрагменты ДНК плазмид ТС-7, ДТ-19 и ДТ-25. На рис.22 схематически изображено строение указанных плазмид. Отличия в их строении относятся линь к одному сайту на ДНК, указанному стрелкой. Ниже приведены нуклеотидные последовательности фрагмента ДНК, включающего все замены. Видно, что ДНК ДТ-19 отличается от ДНК ТС-7 делецией одной бС-пары. Другая ДНК ДТ-25 от ДНК ТС-7 - делецией более протяженного участка, содержащего четыре АТ-пары и две ос-пары.

На рис.23 приведены на одном графике профили плавления всех трех ДНК. Видно, что профили плавления отличаются липь положением второго пика, остальные пики совпадают как по положению, так и по амплитуде. В табл.2 приведены данные, показывающие отличия в положении пика 2 для всех трех ДНК.

Таблица 2. Разница в температурах плавления второго пика (вЩ для ДНК ТС-7, ДТ-19 и ДТ-25.

ДНК (ДТ-19)-(ТС-71 (ДТ-25)-(ТС-7) (ДТ-25ЫДТ-19)

8Т2(°С) -0,14+ 0,03 +0(14+ 0,03 , +0,28+ 0,03

ТЕМРЕРНТиРЕГ/ "С

Рис.23 Профили плавления в 1,6хЗБСДНК плазмид ТС-7 (------),

ДТ-25 <- - -) и ДТ-19 (......).

Из табл.2 следует, что делеция одной 6С пари приводит к сдви» гу пика 2 в сторону меньших температур на 0,12° С (см. рис.23, профили плавления ДНК ТС-7 и ДТ-198). Такий же аффект наблюдается и при делеции 6 п.н. - двух пар сс и четырех АТ. Но как это следует из рис.23, такие изменения в нуклеотидной последовательности приводят к сдвигу пика вправо. В ту же.сторону сдвигается пик 2 при делеции одной 6С пары и четырех АТ (см. на рис.23 профили плавления для ДНК ДТ-19 и ДТ-25), но здесь эффект вдвое больше. Оценки показали, что полученный эффект сдвига пиков на профилях плавления обусловлен изменением среднего сс-содержания соответствующего участка ДНК. Безусловно все исследованные нук-леотидные замены меняют и стэкинг-взаимодействие, и этот эффект должен сказываться на сдвиге пика, но его вклад, видимо, сравним с ошибкой измерений - ±0,02°С.

Таким образом,, в настоящей работе нами показано, что профили плавления ДИК очень чувствительны н нуклеотидным заменам, так что удается регистрировать делецию одной пары нуклеотидов. Отметим, что из проведенных оценок следует, что замена одной пары нуклеотидов.например, АТ на GC должна приводить к вдвое большему эффекту, чем делеция.одной пары.

Другой способ, которым можно обнаружить точечные замены, является гель-электрофорез в градиенте денатурирующих агентов. Этот метод был предложен в работах Лермана и др. (1979,1980), и он основан на эффекте резкого торможения ДНК в геле в условиях ее частичной денатурации. Суть эффекта заключается в том,, что образование денатурированной области в ДНК резко тормозит ее движение при электрофорезе в мелкопористом геле. Нами было показано» что эффект торможения сильно зависит от размера денатурированной области так* что участок из 300 п.н. уменьшает подвижность всей ДНК в~20 раз. Регистрируется эффект задерживания благодаря тому, что в геле создается градиент концентрации денатурирующих агентов, и подбираются такие условия, которые обеспечивают плавление ДНК. Вариации в нуклеотидной последовательности» приводящие к изменению температуры плавления участка, ответственного за торможение, , меняют положение фрагмента на алектрофоретограмме. В работах Лермана (1980,1985) показано, что в ряде случаев таким' способом удается регистрировать отдельные нуклеотидные замены. Поскольку представляет интерес сопоставление возможностей обоих методов, мы провели такое сопоставление на одних и тех жё объектах. -

На рис.24 приведены результаты разделения фрагментов ДНК ДТ-19 и ТС-7. По данным плавления здесь должен быть эффект 0,1<Гс и оказалось, что расстояние между профилями денатурации, действительно отвечает ожидаемому.

Рис.24 Гель-электрофорез фрагментов ДНК ТС-7 и ДТ-19 в градиенте концентрации денатурирующих агентов (формамид с мочевиной). Движение ДНК происходит сверху вниз. Концентрация денатурантов возрастает слева направо.

Полученные данные показывают, что оба метода, электрофорети-ческий и спектрофотометрический, обладают очень высокой чувствительностью к нуклеотидной последовательности ДНК. Однако, информация, получаемая каждым из методов различна: электрофорез дает в основном качественную картину, и он очень эффективен при сопоставлении сразу нескольких препаратов ДНК. Спектрофотоыетри-ческая регистрация плавления позволяет количественно охарактеризовать изменения в нуклеотидной последовательности. Такие возможности метода плавления ДНК могут оказА*ься полезными в ряде случаев, когда необходимо выявить лональные изменения в куклео-

тиднсй последовательности ДНК. Тан, я работах Майерса и др. (1985) здектро<{оретичоский метод бил применен для выявления мутаций, припоикщих : появлению нлслецстг.енных заболеваний у человека. Безусловно ось метода не могут конкурировать с прямым сек-ьенирпванием, дающим полную информацию об изменениях в последо-гатьлыюпти, но при исследовании одних и тех же объектов такое в;шное обстоятельство, что с помощью плавления удается получить информации гораздо бистрее, может оказаться полезным. В Ы В О л н

1. Получены профили плавления с высоким разрешением ДНК из различных источников. Показано, что ДНК генома, размер которого не превышает 10^ п.н., обладает тонкой структурой кривой плавления, характеризующейся наличием на дифференциальной кривой пиков шириной 0,5°С, и что эффект тем ярче выражен, чем короче ДНК. Проведено сопоставление полученных данных с расчетными, что позволило экспериментально обосновать теоретическую модель, объясняющую эффект тонкой структуры.

2. Разработана методика избирательной модификации глиоксалем в ДНК участков, расплавленных при выбранной температуре в интервале плавления ДНК, в результате которой предотвращается их рена-турация. Сочетание данной методики с последующим электронномикро-скопическик анализом ДНК позволило разработать принципиально новый экспериментальный подход к анализу плавления ДНК, с помощью которого удается получить информацию о характере плавления ДНК при любой температуре в интервале плавления.

3. Для ДНК плазмиды СЫЕ1,. обладающей ярко выраженной тонкой структурой кривой плавления, получены электронноыикроскопические Карты плавления| отвечающие 9 выбранным точкам в интервале плавления. Проведено сопоставление плавящихся участков ДНИ с пиками ча профиле плавления. Анализ полученнноА полной картины плавления ДНК Соепривел к .следующий выводам:

* Процесс плавления ДИК представляет собой поэтапное выплавление. ё, ней участков размером от-200 до-500 п.н., каждый участок пла-¿итей в^тервале температур в несколько десятых долей градуса, Зто является прямым подтверждением модели возникновения тонкой

структуры кривых плавления ДНК.

- Экспериментально обоснована модель возникновения неравновео-ности при плавлении ДНК. Показано, что она связана с образованием стадий, приводящих к уменьшению числа границ между раиплав-леннши и спиральными участками.

4. Предложен подход пряного определения одного из важных параметров теории перехода спираль-клубок - фактора кооперативности, основанный на изменении температуры плавления определенного участка ДНК при изменении состояния его границ.

Измерена разница и температурах плавления одного и того же участка ДНК СоСЕI в тех случаях, когда он расположен внутри спирали ДНК (кольцевая ДНК расш,еплялась рестриктазой Зша1 ) и на ее конце (расщепление рестриктазой ЕсоИ) лри различной концентрации N5?". Путем прямого сравнения результатов эксперимента с расчетами впервые надежно определено значение фактора кооперативности

к

б0, которое оказалось равным (2,5-5) ГО .

5. Проведена прямая проверка теории перехода спираль-клубок в ДНК. Для этой цели для ряда ДНК с известном нуклеотидной последовательностью были получены профили и соответствующие им карты плавления. Проведено сопоставление экспериментальных данных с теоретическими. Показано, что полного совпадения распетых профилей с экспериментальными получить не удается. В то сниои время получено вполне удовлетворительное совпадение расчетных карт плавления ДНК с экспериментальными, что указывает нп адекватность теории в отношении предсказания топографии плавления ДНК. Сформулированы возможные "причини расхождения теории и опыта.

Проведен анализ применимости к описании) плавления ДНК модели двух состояний. Показана ограниченность этой модели в отношении детального описания плавления ДНК.

6. Исследованы кинетические эффекты при плавлении ДНК, вызванные денатурацией неравновесно плавящихся участков. Полученные результаты можно суммировать следующим образом:

-Обнаружено влияние скорости нагревания на профили плавления ДНК Показано, что десятикратное изменение скорости нагревания приводит к смещению в область высоких темперят;»,/ лпыь тех пиков, ко-

торые соответствуют стадиям неравновесного плавления ДНК. Величина эфФекта( о'.- составлял десятые доли градуса) оказалась в удовлетворительней* согласии с теорией.

-Детально исследован процесс возникновения неравновесности плав- . ления участка, окруженного расплавленными областями. Проведен анализ процессов денатурации и ренатурации, в различной степени удаленных от равновесия. Получены количественные оценки времен денатурации и ренатурации отдельных участков при различной ионной силе.

-Показано, что скорости процессов как денатурации участков, окруженных расплавленными областями, так и ренатурации расплавленных областей с концов резко зависят от температуры. Наоборот, скорость образования спирали путем нуклеации расплавленной области, слабо зависит от температуры.

7. Исследованы возможности метода плавления ДНК с высоким разрешением для изучения особенностей в первичной структуре ДНК: - Проведен детальный анализ ранних стадий плавления ДНК фага Т7. Предсказано наличие в этой ДНК нескольких АТ-богатых участков размером-ЮО п.н. Разработан подход к прямой локализации этих участков, основанный на расщеплении молекулы ДНК нуклеэой 51 по этим участкам после их предварительной модификации глиоксалем. Электронномикроскопический и рестрикционный анализ получающихся. Фрагментов позволил выявить и локализовать на молекуле ДНК Фага .Т7 4 АТ-богатых участка. Показано, что они располагаются в важных в функциональном отношении.местах фагового генома. -Изучен экспериментально вопрос о влиянии замен в нуклеотидной последовательности ДНК на профиль плавления. Получены профили плавления ДНК плазмид, характеризующихся строго локализованными делециями в несколько пар нуклеотидов. Показано, что Даже деле-ция одной пары приводит к смещению определенного пика на профиле плавления ДНК. В этом случае эффект составил 0,12+0,03°С.

Перечень работ, в которых отражены основные научные результаты диссертации.

1. Любченко Ю.Л., Трифонов Э.Н., Лазуркин Ю.С., Франк-Каменец-кий М.Д. Обнаружение легкоплавких мест в молекуле ДНК Лага Т2.- Мол.биол., 1971, т.5, вып.5, с.772-779.

2. Баев А.С., Любченко Ю.Л., Лазуркин Ю.С., Трифонов Э.Н., Франк-Каменецкий М.Д. Изучение легкоплавких участков ДНК фага Т2 с помощью электронной микроскопии и кинетического Формальдегидного метода.-Мол.биол., 1972, т.6, вып.:, с.760-765.

3. Любченко Ю.Л., Баев А.С., Трифонов Э.Н., Лазуркин Ю.1 Франк-Каменецкий М.Д. Исследование ранних стадий термодгн.ч-турации ДНК фага Т2.-С6. "Конформж¡ионные изменения биополимеров в растворах". М.:Наука, 1973, с.15-22.

4. Lazurkin Yu.S^ Pavlov V.M,t Lyubchenko Yu.L,, Frank-Kuif-.. netskii M.D., Eerestetskaya I.V. On the parameter:-, ¡>t' DNA melting curves in the low Na ion concentration ran^.-Biopolyoers, 1975, v.I4, n.7, pp.I55I-I552.

5. Lyubchenko Yu.L., Vologodskii A.V., Lazurkin Yu.S., Frank-Kamenetskii M.D., Gause G.G. Fine structure of DNA melt curves,-Eiopolymers, 1976, v.15, no 6.FP. I0I9-I036.

6. Pavlov V.M., Lyubchenko Yu, L.f Eorovik A.S.^, Lazurkin Yu.S. Specific fragmentation of T7 phage DNA at low-melting range,-Nucleic Acids Res., 1977, v.4, n.II, pp.4053-4062.

7. Pavlov V;Mlt Lyubchenko Yu.L. A new method for recording of DNA differential melting curves.- Biopolymers, 1978,

v,17, no 7, PP.775-778.

8. Lyubchenko Yu.L., Eorovik A.S., Kalambet Yu.A., Shitov V.T., Golovanov E.I. Mapping of different stages of i.oiicln EI plasmid DNA melting. In: DNA-recombination, interaction:' and repair. -Zadrazil S,, Sponar J. Eds., Perg. Pre ;.s; Oxford N.Y 1980,.ppЛ55-162.

9. Eorovik A.S., Kalambet Yu.A., Shitov V. ., Golovanov E.j Lyubchenko Y.u.L-. ; Equilibrium melting ColEI DNA: Fie< !r n microscopic visualisation .-Nucleic Acids Res., 19Я0. v.8 n 9, pp.4165-4184.

Ю. Lyubchenko Yu.L., Vologodskii A.V Frank-Kamenetskii M.D. Direct, comparison of theoretical and experimental melting profiles for H? H ФХ174 DMA.- Nature, 1978, v.271, pp.2831 .

1!. Ti kr.honenko T.t., Dubichev A.G., Lyubchenko Yu.L., Kvitko N.P., Chaplygina N.M., Kalinina T.I., Dreizin R.S., Narodit яку F.S The distributIon of guanine-cytosine pairs in adenovirus DNAs .-J.Gen.Virol., T98I, v.54, pp.425-429.

12. Lyubchenko Yu. L., Kalambet Yu. A.( Lyamichev V.I., Eorovik A . 3 . A comparison of experimental and theoretical melting Tips for replicative form of ФХ174 DNA. -Nucleic Acids Res., 1982, v.IO, n.6, pp.1867-1876.

13. Lyamichev V.I., Panyutin I iG. , Cherny D.I., Lyubchenko Yu.L. Localization of low-melting regions in phage T7 DNA.-Nuc-leic Acids Res., 1983, v'.II, n7, pp .2165-2176.

14. Пюбченко Ю.Л., Завальский Л.Ю., Боровик А.С. Реакция взаимодействия рестриктазы EcoRI с ДНК плазмиды ColEl.-Мол.биол ., 1982, т.16, вып.4, с.871-877..

[5. Kozyavkin S.A., Lyubchenko Y .L. The nonequilibrium character of DNA melting: effects of the heating rate on the fine structure of melting curves.- Nucl. Acids Res., 1984, v.-12, no 10, 4329-4349: ,' . '

16. Боровик А.С., Любченко Ю.Л., Народйцкий B.C., Тихоненко Т.И. Электронномикроскопические карты плавления ДНК SA7.- Мол. биол., 1984, т.18, вып.6, с.1634-1638.

17. Kalambet Yu.A.( Borbvik A.S., Lyamichev V.X., Lyubchenko Yu L. Electron microscopy of the melting of sequenced DNA,-Bio polymers, 1985, v.24, p.359->377.

18. Kozyavkin S.A., Naritsin.D .E,, ' Lyubchenko Yu.L. The kinetics of DNA helix-coll subtransitionsJourn. Biomolec. Struct. Dyn., 1986, v.3, n.4, pp.689-704.

Lb. Razlutskii I.V.f 'Shlyakhtenko L.S., Lyubchenko Yu.L. The effect of nucleotide substitutions on DNA denaturation -profiles.-. Nucl. Acids Res., 1987, v.15, n.I6, pp.66fi5-6676.

20. Нариций Д.В., ЛгобченКо Ю.Л. Плавление коротких фрагментов ДНК. > Определение»зависимости константы нуклеации спирального участка ДНК от ионной силы. - Колек.биол., 1988, т.22,

вып.1, с.242-248.

21. Амосова O.A., Черный Д.И., Любченко В.Д., Дужак А.Б., Лобова H.H., Подгорный В.Ф. Сравнительный анализ строения дпу-спиральных РНК киллерных штаммов дрожжей ■ n.-uwy. » cerevisiae.- Молек.биол., 1988, т.22, вып.1, с.224-230.

22. Разлуцкий И.В., Шляхтенко J1 .С., Любченко ЮЛ. Влияние i ук ■ леотидных замен на профили денатурации ДНК,- Биополимори и клетка, 1988, т.4, вып.5, с.225-227.