ВЫСОКОМОЛЕКУЛЯРНЫЕ СОЕДИНЕНИЯ, Серия А, 2008, том 50, № 2, с. 324-331
РАСТВОРЫ
УДК 541.64:532.77
ПОЛИАМФОЛИТНЫЕ РАСТВОРЫ В ПРИСУТСТВИИ ИОНОВ СОЛИ
w ______-t
РАЗНОЙ ВАЛЕНТНОСТИ: КОМПЬЮТЕРНОЕ МОДЕЛИРОВАНИЕ1
© 2008 г. О. А. Балашова, А. С. Павлов, П. Г. Халатур
Тверской государственный университет 170002 Тверь, Садовый пер., 35
Поступила в редакцию 30.01.2007 г.
Принята в печать 29.06.2007 г.
Методом стохастической динамики исследованы свойства разбавленных растворов регулярных по-лиамфолитов. Среднечисленная концентрация полимерных звеньев в растворе не превышала 0.038, что ниже концентрации перекрывания клубков при N = 1024. Рассмотрено влияние заряда (валентности) ионов фонового электролита на размер макромолекулы и распределение соли внутри полимерного клубка при разной температуре. Обнаружено, что для исследованного интервала температур (1 < T < 2e/kB) конденсация ионов соли происходит в случае многозарядных систем типа (+4 -4) и (+4 -2). Максимальное проникновение ионов соли в объем макромолекулы наблюдалось при добавлении фонового электролита типа (+2 -2).
ВВЕДЕНИЕ
Заряженные полимеры привлекают значительное внимание благодаря своим уникальным свойствам. К таким системам относятся полиэлектролиты и полиамфолиты. В последнее время большое внимание уделяется изучению растворов этих полимеров в присутствии многозарядных контрионов и ионов фонового электролита.
Из всего многообразия работ, посвященных исследованию поведения полиэлектролитов при добавлении ионов соли разной валентности, можно выделить работы [1-3]. Радиус инерции отдельной цепи полиэлектролита уменьшается с ростом либо концентрации полимера, либо концентрации соли [4]. Двухвалентные ионы способствуют более сильному сжатию полиэлектролитной цепи по сравнению с одновалентными [4]. Приведенные результаты получены для системы из нескольких полиэлектролитных цепей методами компьютерного моделирования. В работе [5] экспериментальными методами обнаружено, что трехвалентные ионы соли обычно конденсиру-
1 Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (коды проектов 07-03-00385 и 05-03-32952а).
E-mail: [email protected] (Балашова Ольга Андреевна).
ются на ДНК, в то время как только двухвалентные - на ^-актине и вирусах ^/-семейств, а одновалентные не конденсируются ни на одном из этих объектов. Также обнаружено [6], что муль-тивалентные ионы, конденсируясь на ДНК, приводят к формированию торов и стержней. В работе [7] показано, что при малой концентрации трехвалентные ионы вызывают компактизацию макромолекулы ДНК. Авторы отмечают, что взаимодействие изученных ими ионов носит специфический характер.
Работ по изучению полиамфолитных систем в присутствии ионов соли немного [8-12], но имеется ряд исследований, посвященных белкам [13, 14]. Системы, содержащие большое количество частиц, требуют значительной затраты компьютерного времени; чтобы этого избежать, используют экранированный потенциал Кулона [15]. Данное приближение хорошо работает для слабо заряженных макромолекул в растворе одновалентной соли, но в системах с мультивалентными противоионами и(или) с низкой диэлектрической проницаемостью становится менее применимым [13].
Экспериментально обнаружено [16], что увеличение концентрации соли ведет к значительным колебаниям стабильности белков, вероятно,
из-за изменения кажущегося гидрофобного эффекта и экранирования заряда. В теоретической работе [17] показано, что с добавлением противо-ионов увеличивается стабильность белка. Добавление ионов может сильно менять термодинамические и физические свойства белков в растворе, что подтверждается экспериментально при исследовании процессов денатурации, всаливания и высаливания белков [14].
Известно [8], что при небольшом количестве соли в полиэлектролитном режиме экранирование заряда уменьшает размер цепи. Дальнейшее добавление соли вызывает набухание макромолекулы, подтверждаемое работами [11, 12].
В экспериментальных работах [9, 10] показано, что добавление соли влияет на вязкость и растворимость полиамфолитов. Нейтральные поли-амфолиты не растворимы в чистой воде вследствие коллапса цепи, однако они становятся растворимыми при большой концентрации соли вследствие полного экранирования зарядов на цепи [9]. Известно [9], что природа электролита влияет на свойства раствора. В частности, чем больше поляризуемость аниона, тем лучше солю-билизация нейтрального полиамфолита. Следует отметить, что двухвалентные катионы более эффективны, чем одновалентные. При низкой концентрации соли уменьшается вязкость, что является типичным полиэлектролитным эффектом [10]. При дальнейшем добавлении соли вязкость увеличивается.
Имеются работы по изучению полиамфолит-ных растворов в присутствии соли методом молекулярной динамики [11, 12]. Авторы работы [12] исследовали раствор, состоящий из шести поли-амфолитных цепей длиной 32 мономерных звена со случайным распределением зарядов в присутствии положительно и отрицательно заряженных противоионов, причем количество ионов соли было таким, чтобы соблюдался принцип электронейтральности. Показано, что радиус инерции системы случайных гибкоцепных полиамфолитов уменьшается с добавлением соли в низкотемпературном режиме. В случае отдельных цепей наблюдалось незначительное увеличение размеров цепи.
Цель нашего исследования - изучение влияния мультивалентных ионов соли на поведение регулярных полиамфолитов и условий конденсации
противоионов фонового электролита на макромолекуле.
МОДЕЛЬ СИСТЕМЫ И МЕТОД МОДЕЛИРОВАНИЯ
Модель полиамфолита представляла собой молекулу линейного полимера, которая имела положительные и отрицательные заряды вдоль цепи. Заряды чередовались и распределялись регулярно (блоками по п звеньев). Противоионы отсутствовали. Расчеты проводили для макромолекулы, состоящей из 1024 мономерных звеньев при длине одноименно заряженного блока, равной 32.
Применяли модель свободносочлененной цепи, в которой мономерные звенья соединены связями фиксированной длины (1а), ван-дер-вааль-совы радиусы мономерных звеньев а считали одинаковыми. Величина заряда на мономерных звеньях полиамфолита одинакова и равна ±1 элементарного заряда. Ионы соли моделировали как заряженные сферические частицы радиуса 1а, отличие между катионом и анионом сводилось к величине и знаку заряда. Рассматривали пять случаев распределения зарядов соли: 1 - катионы имели заряд q1, равный + 4, анионы - q2 = -1, общее число частиц фонового электролита п+ = 256, п_ = 1024; 2 - 2 = ±4, п+ = 256, п_ = 256; 3 - ^ 2 = ±2, п+ = 512, п- = 512; 4 - q1 = +4, q2 = -2, п+ = 256, п- = = 512; 5 - qlJ 2 = ±1, п+ = 1024, п- = 1024. Количество ионов соли менялось в зависимости от их валентности таким образом, что суммарный заряд ионов каждого знака оставался неизменным. При увеличении заряда число частиц пропорционально уменьшалось. В целом система оставалась электронейтральной.
Короткодействующие силы притяжения между несвязанными частицами описывали потенциалом Юкавы
иа (ти) =
£ава
1-1^ Г
2-,2 су -I
Го < г1} < г с
(1)
0, г1] > г с
Параметр еар задает качество растворителя. В наших расчетах еар = 0 (для хорошего растворителя) или 0.3 (для плохого) и радиус обрезки гс = = 2.8а.
На малых расстояниях г < 21/6а кроме потенциала притяжения (1) включали потенциал Леннар-да-Джонса
иеу (Гц) = <
а)12 - (а) V"
4 £
0, Гц > Го
Гц < Го,
(2)
где а = £ = 1 и г0 = 21/6а - радиус обрезки. Параметр £, входящий в уравнение (2), контролирует шкалу энергии, в то время как а отражает равновесное расстояние между частицами.
Электростатические взаимодействия между заряженными частицами описывали кулонов-ским потенциалом
и с =
1 qiqJ
4П£о£г Гц
(3)
Здесь Гц = |гг - Гц | - расстояние между г и Ц заряженными частицами, qi - заряд данной частицы, £0 -диэлектрическая проницаемость вакуума, £Г - относительная диэлектрическая проницаемость среды.
Расчеты проводили в системе с периодическими граничными условиями в кубической ячейке с ребром 30-45а. Кулоновские взаимодействия учитывали по методу Эвальда [18, 19]. Значения параметров потенциалов (1) и (2) выбраны таким образом, что полимерная система без электростатических взаимодействий оказывается вблизи 6-точки при приведенной температуре Т = 1.0. Таким образом, учитывается изменение качества растворителя с изменением температуры.
В расчетах варьировали температуру (Т = 0.13.0) и валентность ионов соли ^ = 1-4). Общая продолжительность моделирования на стационарном участке траектории для каждого случая составляла (2-5) х 106 временных шагов. Следует отметить, что в работе используется стандартная для компьютерного моделирования система приведенных единиц, в которой длина измеряется в единицах а, энергия - в единицах £, температура Т -в единицах £/кВ. В работах по полиэлектролитной тематике единицей длины часто служит длина Бьеррума 1В = в2/£ГкВТ, представляющая собой
расстояние между двумя элементарными зарядами е в среде с диэлектрической проницаемостью £Г, взаимодействующими с энергией, равной кВТ. Рост 1В эквивалентен понижению температуры при фиксированных параметрах системы.
РЕЗУЛЬТАТЫ И ИХ ОБСУЖДЕНИЕ
Для исследования структуры системы рассматривали функции радиального распределения плотности ^(г) для мономерных звеньев и пар мономерное звено-ионы соли, методика расчета которых изложена в работе [20]. Будем обсуждать зависимости при двух температурах Т = 1.0 и 2.0, для которых наблюдается максимальное различие в поведении приведенных систем.
На рис. 1 представлены внутримолекулярные функции ^(г) звеньев типа АВ для исследуемых систем типа (+4 -1), (+4 -4), (+4 -2), (+2 -2) и (+1 -1), где числа в скобках обозначают заряды катиона и аниона соли, при Т = 1.0. Видно, что для всех систем характерно наличие двух максимумов. Первый более высокий пик, отвечающий расстоянию между частицами 1.15 а, совпадает с минимумом энергии ван-дер-ваальсовых взаимодействий и обусловлен электростатическим вкладом. Второй более низкий максимум на кривых радиального распределения плотности обусловлен связанностью звеньев в цепь. Таким образом, наблюдается только первая координационная сфера, что свидетельствует о нахождении системы в жидкоподобном состоянии, при этом звенья имеют достаточно высокую подвижность относительно друг друга. Следует отметить, что g(г) для звеньев в системе типа (+1 -1) лежит выше других кривых. Это связано с тем, что экранирование внутрицепных электростатических взаимодействий ионами низкой валентности является более слабым, несмотря на большее число ионов.
Для изучения процесса конденсации противо-ионов рассчитаны функции g(г) для пар мономер типа А-катионы и мономер А-анионы соли, показанные на рис. 2. Звенья типа А полиамфолита заряжены положительно. Представленные зависимости отражают резкий рост плотности катионов соли до г ~ 2а, а затем ее выход на плато (рис. 2а, кривая 2) или дальнейшее медленное увеличение (кривые 1, 3, 4). Вследствие того, что звенья типа А и катионы соли С(+) одноименно
Рис. 1. Функции радиального распределения плотности g(r) мономерных звеньев типа АВ при Т = 1.0 для различных систем: 1 - (+4 -1), 2 - (+4 -4), 3 - (+4 -2), 4 - (+2 -2) и 5 - (+1 -1).
заряжены, на малых расстояниях превалирует отталкивание между зарядами. На расстоянии г > 2а взаимодействие зарядов оказывается экранированным противоположно заряженными звеньями типа В, что и подтверждается ходом кривых 1-4. На рис. 2а видно, что поведение системы 5 отличается от поведения других систем. Это, вероятно, связано с тем, что среднечисленная плотность одно-однозарядной соли выше, а ее проникновение внутрь глобулы незначительно.
На рис. 26 представлена функция g(r) анионов соли относительно положительно заряженных звеньев полиамфолитной цепи. Как и следовало ожидать, противоположные заряды притягиваются, и для исследуемых систем наблюдается наличие двух максимумов при г = 1.15а и 2а. Первый максимум соответствует увеличению концентрации в первой координационной сфере, где концентрация анионов в 2-3 раза выше, чем в среднем по системе. Второй максимум обусловлен связанностью звеньев в цепь. Если расстояние между звеном А и анионом равно г, то ближайшие по цепи мономерные звенья будут на расстоянии г + 1, что и следует из рис. 26. Заметим также, что ход кривой для системы типа (+1 -1) сильно отличается от поведения других зависимостей. В отличие от систем 1-4, для кривой 5 в пер-
вой координационной сфере наблюдается пониженное содержание частиц фонового электролита. Однако с ростом г количество ионов соли возрастает до среднечисленной плотности фонового электролита в ячейке. Такое поведение g(r) анионов можно объяснить только низкой концентрацией фонового электролита внутри глобулы полиамфолита.
Для более детального исследования распределения ионов фонового электролита мы рассматривали нормированную функцию радиального распределения плотности р(г1п) ионов соли относительно центра инерции цепи при Т = 1.0 (рис. 3).
Зависимости 1 и 5 на рис. 3 ведут себя одинаково, т.е. вблизи центра инерции цепи до г < 6а плотность ионов низкая и незначительно меняется. С увеличением расстояния резко растет р до значения среднечисленной плотности ионов при их равномерном распределении в ячейке. Это свидетельствует о попадании ионов низкомолекулярного электролита только во внешние части макромолекулы. В наибольшей степени подобное поведение характерно для одно-однозарядного электролита, что объясняет отсутствие ионного окружения у мономерных звеньев полиамфо-лита (рис. 26, кривая 5). Для зависимости 3 характерно подобное поведение, однако вблизи
#АС(+)
3
6
9 г
Рис. 2. Функции g(г) для пар мономерное звено-ионы соли типа АС (+) (а) и АС (-) (б) при Т = 1.0 для различных систем: 1 - (+4 -1), 2 - (+4 -4), 3 -(+4 -2), 4 - (+2 -2) и 5 - (+1 -1).
поверхности полиамфолитной глобулы происходит незначительная конденсация ионов фонового электролита.
На кривой 2 (рис. 3) имеется максимум при г - 9а, который свидетельствует об увеличении конденсации противоионов на поверхности глобулы. Фоновый электролит типа (+2 -2) равномерно распределяется по объему ячейки, что подтверждается ходом кривой 4.
Для исследования структуры полиамфолитной макромолекулы были изучены функции р(г1п) для мономерных звеньев типа А относительно центра инерции цепи при Т = 2.0 (рис. 4). Для всех представленных систем на расстояниях г < 7 а равновесное распределение частиц выходит на пла-
то, положение которого практически не зависит от типа распределения зарядов в ионах соли. Для Т = 2.0 такой вид функции радиального распределения плотности звеньев относительно центра инерции цепи соответствует рыхлой глобуле.
Количественная информация об изменении конформации цепи может быть получена из анализа среднеквадратичного радиуса инерции Я^.
2
Рисунок 5 показывает как функцию 1В при различных типах распределения зарядов соли. Размер цепи уменьшается с ростом длины Бьер-рума (или с понижением температуры).
При 1В = 0 взаимодействие между зарядами отсутствует, и все определяется исключенным объемом, при этом цепь ведет себя подобно набухшему клубку. При низких значениях длины Бьерру-ма наблюдаются большие значения размеров цепи, т.е. молекула находится в набухшем клуб-кообразном состоянии. С увеличением 1В начинает больше влиять электростатика, способствующая сжатию макромолекулы в глобулу, при этом все зависимости совпадают в области 1В > 4, и размеры молекулы полиамфолита практически не зависят от наличия соли. Подобному поведению макромолекулы способствует и тот факт, что при 1В > 1 система оказывается в условиях плохого растворителя даже без учета электростатических взаимодействий. Из анализа этих данных можно заключить, что размер цепи зависит от заряда ионов соли. При этом следует заметить, что добавление ионов соли с валентностью больше двух значительнее меняет размер цепи по сравнению со случаем одновалентной соли. Отметим также, что наибольшее значение размеров цепи наблюдалось при добавлении фонового электролита типа (+4 -2). Вероятно, это можно объяснить значительным проникновением ионов соли в объем макромолекулы и ростом электростатического вклада.
На рис. 6 представлена "мгновенная фотография" полиамфолита в присутствии ионов соли типа (+4 -4) при Т = 1.0, соответствующей 1В = 1.0. Из анализа фотографии можно заключить, что ионы с валентностью, равной четырем, образуют ионные пары, а также более сложные агрегаты, и при этом не стремятся проникнуть внутрь полиамфолитной молекулы.
5
10
15 г
1 ' 1Г
Рис. 3. Функции радиального распределения плотности ионов соли относительно центра инерции цепи при Т = 1.0 для различных систем: 1 - (+4 -1), 2 - (+4 -4), 3 - (+4 -2), 4 - (+2 -2) и 5 - (+1 -1).
0.15
0.10
0.05
10
15
Рис. 4. Функции радиального распределения плотности мономерных звеньев типа А относительно центра инерции цепи при Т = 2.0 для различных систем: 1 - (+4 -1), 2 - (+4 -4), 3 - (+4 -2) и 4 - (+2 -2).
5
г
Из сказанного выше следует, что наибольшее различие в поведении исследуемых систем проявлялось при температурах 1.0 и 2.0. Заметим, что максимальное значение размеров цепи наблюдалось при добавлении фонового электролита типа
(+4 -2), что, вероятно, можно объяснить значительным проникновением ионов соли внутрь макромолекулы и ростом электростатического вклада. Из функций радиального распределения частиц относительно центра инерции цепи видно,
2
Рис. 5. Зависимость Rg от длины Бьеррума для Рис. 6. "Мгновенная фотография" нолиамфоли-
различных систем: 1 - (+4 -1), 2 - (+4 -4), 3 - та в нрисутствии ионов соли тина (+4 -4) нри
(+4 -2), 4 - (+2 -2) и 5 - (+1 -1). T = 1.0.
что конденсация соли происходит в случае фонового электролита типов (+4 -4) и (+4 -2); во всех остальных случаях частицы низкомолекулярного электролита попадают преимущественно во внешние части полиамфолитной глобулы.
Полученные результаты согласуются с работами других авторов [9, 12], показывающими, что при добавлении соли происходит набухание макромолекулы полиамфолита. Следует отметить, что ионы соли с валентностью больше единицы приводят к большему увеличению размеров цепи по сравнению с однозарядными ионами.
ЗАКЛЮЧЕНИЕ
Расчеты, выполненные методом стохастической динамики, для разбавленного раствора регулярного полиамфолита в присутствии фонового электролита показывают существенные изменения структурных свойств полиамфолита в зависимости от температуры и валентности ионов соли. Соответствующие заключения делаются из анализа функций радиального распределения плотностей мономеров и ионов соли относительно центра инерции макромолекулы. Найдено, что максимальное проникновение ионов соли в объем макромолекулы происходит при добавлении фонового электролита типа (+2 -2). Обнаружено, что наибольшее значение размера цепи наблюда-
ется в системе в нрисутствии ионов соли тина (+4 -2). Для исследованного интервала темнера-тур (1 < T < 2e/kB) конденсация ионов соли нроис-ходит в случае соли тинов (+4 -4) и (+4 -2).
СПИСОК ЛИТЕРАТУРЫ
1. Barrat J.-L, Joanny J.-F. // Adv. Chem. Phys. 1996. V. 94. P. 1.
2. Oosawa F. Polyelectrolytes. New York: Marcel Dekker, 1971.
3. Kuhn P.S., Barbosa M.C. // Physica A. 2005. V. 357. № 1. P. 142.
4. Liu S, Ghosh K, Muthukumar M. // J. Chem. Phys. 2003. V. 119. № 3. P. 1813.
5. Butler J.C., Angelini T, Tang J.X., Wong G.C.L. // Phys. Rev. Lett. 2003. V. 92. № 2. P. 028301.
6. Bloomfield V.A. // Curr. Opin. Struct. Biol. 1996. V. 6. № 3. P. 334.
7. Касъяненко НА., Дьяконова HE., Фрисман Э.В. // Молек. биология. 1989. Т. 23. № 4. С. 975.
8. Dobrynin A.V., Colby R.H., Rubinstein M. // J. Polym. Sci., Polym. Phys. 2004. V. 42. № 19. P. 3513.
9. Corpact J.-M., Candau F. // Macromolecules. 1993. V. 26. № 6. P. 1333.
10. Ohlemacher A., Candau F., Munch J.P., Candau S.J. // J. Polym. Sci., Polym. Phys. 1996. V. 34. № 16. P. 2747.
11. Tanaka M., Grosberg A.Yu., Tanaka T. // Langmuir. 1999. V. 15. № 12. P. 4052.
12. Tanaka M, Grosberg A.Yu., Tanaka T. // J. Chem. Phys. 1999. V. 110. № 16. P. 8176.
13. Lund M, Jönsson B. // Biophys. J. 2003. V. 85. № 5. P. 2940.
14. Ibragimova G.T, Wade R.C. // Biophys. J. 1998. V. 74. № 6. P. 2906.
15. Carlsson F., Linse P., Malmsten M. // J. Phys. Chem. B. 2001. V. 105. № 38. P. 9040.
16. Kohn W.D., Kay C.M., HodgesR.S. // J. Mol. Biol. 1997. V. 267. № 4. P. 1039.
17. York D.M., Darden T.A., Pedersen L.G., Anderson M.W. // Biochemistry. 1993. V. 32. № 6. P. 1443.
18. Allen MP, Tildesley D.J. Computer Simulation of Liquids. Oxford: Clarendon Press, 1987.
19. Pollock E.L., Glosli J. // Comput. Phys. Commun. 1996. V. 95. № 2-3. P. 93.
20. Балашова O.A., Павлов A.C., Халатур П.Г. // Вы-сокомолек. соед. А. 2007. Т. 49. № 3. С. 481.
Polyampholyte Solutions in the Presence of Salt Ions Having Different Valences: Computer Simulation
O. A. Balashova, A. S. Pavlov, and P. G. Khalatur
Tver State University, Sadovyi per. 35, Tver, 170002 Russia e-mail: [email protected]
Abstract—The properties of dilute solutions of regular polyampholytes in the presence of a low-molecular-mass salt were studied by means of stochastic dynamics simulation. The number-average concentration of polymer units in solution did not exceed 0.038, a value that is below the coil overlap concentration at N = 1024. The effect of charge (valence) of ions of an added salt on the macromolecule size and distribution of the salt inside the polymer coil at different temperatures was considered. It was found that the condensation of salt ions over the temperature range examined (1 < T < 2e/kB) takes place in the case of multiply charged systems of the (+4 -4) or (+4 -2) type. The maximal penetration of salt ions into the volume of the macromolecule was observed in the case of a (+2 -2) salt.