Солоненко А.П., Голованова О.А., Бельская Л.В.

ТЕРМОДИНАМИЧЕСКИЙ РАСЧЕТ ПРОЦЕССА КРИСТАЛЛИЗАЦИИ ГИДРОКСИЛАПАТИТА В ПРИСУТСТВИЯ МИКРОЭЛЕМЕНТОВ


Кальций, фосфор, кислород и водород являются одними из самых распространенных элементов на поверхности нашей планеты. Данные элементы сочетаются друг с другом в различных комбинациях, это обуславливает большое разнообразие фосфатов кальция. Среди них орто-(PO43–), мета-(PO3–), пиро- (P2O74–) и полифосфаты ((PO3)nn–). Известно более десятка солей, которые могут быть образованы этими четырьмя элементами. Но наиболее чаще из всех известных фосфатов кальция как в живой, так и неживой природе встречается минерал – гидроксилапатит [Ca10(PO4)6(OH)2]. Как полезное ископаемое апатит широко используется в промышленности для получения фосфорсодержащих удобрений, фосфорной кислоты, пищевых добавок и проч. В биологических системах гидроксилапатит является основным неорганическим компонентом физиогенных (кости, зубы) и патогенных (зубные, слюнные и мочевые камни, атеросклеротические повреждения) твердых тканей человека. Стоит отметить, что гидроксилапатит практически не встречается без примесей других элементов. Например, ионы кальция могут быть частично замещены Sr, Ba, Mg, Mn, K, Na, Fe, Cu; CO32–, AsO43– и VO43– могут замещать ортофосфат-ионы; гидроксид-ионы могут быть до некоторой степени быть замещены фтором [Dorozhkin, 2007]. Замещения в структуре вещества обусловлены тем, что в среде, в которой происходит формирование основного фосфата кальция, всегда присутствуют различные примеси. Так и в организме человека кристаллизация гидроксилапатита происходит в сложной по составу среде, содержащей наряду с макрокомпонентами (ионы – Ca, Р, Na, K и некоторые другие) также ионы, концентрации которых находятся на уровне 10-6 моль/л (Cu, Fe, Al, Mn, Zn и др.).
Действие основных ионов биологических жидкостей (например, слюны) установлены, а роль микроэлементов в процессах образования биоминералов как физио-, так и патогенного характера до настоящего времени остается недостаточно изученной.

Существуют предположения о том, что микроэлементы ротовой жидкости могут либо препятствовать кристаллизации солей кальция [Wandt, Underhill, 1988] за счет адсорбции на активных центрах твердой фазы, либо выступать в качестве зародышей кристаллизации, связывая анионы раствора в малорастворимые соединения.
Важную информацию об образовании малорастворимых фаз в многокомпонентных биологических системах можно получить на основании термодинамических расчетов, в ходе которых возможно определение образования соединений в присутствии различных примесей. Ранее нами проведен термодинамический расчет возможности осаждения малорастворимых солей кальция и магния из раствора, приближенному по неорганическому составу к слюне человека [Бельская, Голованова, 2008]. Однако при этом не учитывалось влияние микроэлементов на фазообразование в растворе.
Цель работы – определение условий осаждения фосфатов кальция различного стехиометрического состава, а также возможности образования малорастворимых соединений основными микроэлементами (Zn, Fe, Cu, Al) смешанной слюны человека на основании термодинамического расчета кристаллизации.

В качестве исследуемой системы был выбран раствор, минеральный состав (неорганические макро- и микрокомпоненты), температура, ионная сила и рН (4.5–8.0) которого соответствуют составу ротовой жидкости здорового взрослого среднестатистического человека (табл.). Значения термодинамических произведений растворимости (рKºs) малорастворимых соединений, которые могут образоваться в исследуемом растворе [CaHРО4•2Н2О, Ca8H(РО4)6•5Н2О, Са10(РО4)6(ОН)2, Са10(РО4)6F2, CaCO3, Ca(OH)2, MgHPO4•3Н2О, МgNH4PO4•6Н2О, Mg3(PO4)2, MgCO3•3Н2О, Mg(OH)2, Zn3(PO4)2, ZnCO3, FeCO3, Zn(OH)2, Fe(OH)2, Al(OH)3], соответствуют данным базы констант нестойкости комплексов SC-database (SCQUERY Version 1.37 (1993)). Для определения возможности и условий образования малорастворимых фаз рассчитывали значения индексов пересыщения (SI) и энергии Гиббса кристаллизации (∆G): осаждение фазы из раствора возможно, если ∆G < 0 и SI > 0.

Таблица
Концентрации компонентов в модельной кристаллизационной среде
№ п/п Компонент Концентрация, моль/л № п/п Компонент Концентрация,
моль/л
1 Кальций 0.0015 8 Магний 0.00065
2 Натрий 0.0139 9 Фториды 1.1E-05
3 Калий 0.0262 10 Алюминий 2.60E-05
4 Аммиак 0.0036 11 Железо 4.98E-06
5 Хлориды 0.0127 12 Медь 5.38E-06
6 Карбонаты 0.0074 13 Цинк 7.28E-06
7 Фосфаты 0.0056 
  
По результатам проведенных расчетов установлено, что в исследуемых системах термодинамически вероятно образование следующих малорастворимых соединений: CaHРО4•2Н2О; Ca8H(РО4)6•5Н2О; Са10(РО4)6(ОН)2; Са10(РО4)6F2; CaCO3 (рис. 1). При рН < 7.00 по убыванию значений SI исследуемые фосфаты кальция можно ранжировать следующим образом: Са10(РО4)6F2 > Са10(РО4)6(ОН)2 > CaHРО4•2Н2О > Ca8H2(РО4)6•5Н2О; тогда как при рН > 7.0 последовательность имеет вид: Са10(РО4)6F2 > Са10(РО4)6(ОН)2 > Ca8H2(РО4)6•5Н2О > CaHРО4•2Н2О. Карбонаты кальция характеризуются отрицательными величинами энергии Гиббса кристаллизации рН ≥ 7.0, однако их значения ниже, чем для фосфатов кальция, поэтому в рамках данной модели их образование термодинамически маловероятно. Осаждение фосфатов магния в исследуемых условиях термодинамически невозможно, так как для них ∆G > 0 и SI < 0.

Помимо солей кальция также термодинамически вероятно образование Zn3(PO4)2, ZnCO3, FeCO3, Cu(OH)2, Fe(OH)3, Al(OH)3 (рис. 2). В широком диапазоне рН, включая область 4.5–8.0, возможно осаждение гидроксидов Fe (III) и Al (III). Гидроксиды кальция, магния, цинка и железа (II) образуются при рН > 8.0 (такие значения рН не характерны для ротовой жидкости). Образование малорастворимых Zn3(PO4)2, ZnCO3 и FeCO3 возможно, начиная с рН > 7.5, при этом процесс кристаллизации исследуемых солей характеризуется малыми величинами энергии Гиббса, много меньше ∆G для фосфатов кальция различного стехиометрического состава, поэтому их образование в рамках выбранной модели термодинамически маловероятно.

В целом по результатам проведенного расчета можно сказать, что с увеличением рН (в исследуемом диапазоне кислотности среды) наблюдается рост пересыщения в системе по всем малорастворимым соединениям. Однако необходимое значение данного параметра для большинства солей и гидроксидов микроэлементов, а также магния не достигается. В изучаемых условиях термодинамически наиболее вероятно осаждение фтор- и гидроксилапатита. Процесс кристаллизации брушита и октакальция фосфата характеризуется меньшими величинами ∆G и SI, следовательно, данные фазы являются метастабильными в условиях эксперимента. Образование малорастворимых соединений Zn, Сu, Fe и Al с анионами исследуемого раствора термодинамически маловероятно при малых  концентрациях  катионов  металлов в данном диапазоне рН и в присутствии
  
Рис. 1. Графики зависимости значений индексов пересыщения от рН раствора для солей кальция. 

Рис. 2. Графики зависимости значений индексов пересыщения от рН раствора для соединений Fe, Cu, Zn, Al.

катионов кальция с концентрацией, во много раз превышающей концентрации микроэлементов. Кроме того, микроэлементы Zn, Сu, Fe и Al в отличие от кальция и магния, склонны к образованию прочных хелатных комплексов с органическими компонентами ротовой жидкости, что еще больше снижает вероятность осаждения соответствующих малорастворимых соединений.

В целом по итогам термодинамического моделирования равновесий в гипотетическом растворе в присутствии ионов тяжелых металлов можно сделать вывод, что результаты проведенных расчетов хорошо согласуются с данными по изучению фазового состава коллекции зубных камней жителей Омского региона [Бельская, Голованова, 2008], согласно которым основным минеральным компонентом всех образцов является гидроксилапатит, в редких случаях встречается брушит (5–10 % от содержания апатита), а собственные фазы микроэлементов ротовой жидкости не идентифицированы. Кроме того, полученные результаты свидетельствуют о том, что ингибирующее влияние микроэлементов на кристаллизацию фосфатов кальция не является следствием связывания анионов раствора с образованием их малорастворимых соединений. Вероятно, главными причинами замедления кристаллизации солей кальция в присутствии ионов тяжелых металлов является адсорбция многозарядных ионов на поверхности растущих кристаллов, а также возможность изоморфного замещения кальция на Меn+, что приводит к дестабилизации структуры гидроксилапатита и возрастанию растворимости фазы [Вересов и др., 2004].

Литература

Бельская Л. В., Голованова О. А. Термодинамическое моделирование процесса образо¬ва¬ния зубного камня // Вестник Омского университета. 2008. № 4. С. 49–53.
Вересов А. Г., Путляев В. И., Третьяков Ю. Д. Химия неорганических биоматериалов на основе фосфата кальция // Российский химический журнал. 2004. Т. 48. № 4. С. 52–64.
Dorozhkin S. V. Calcium orthophosphates // J. Mater. Sci. 2007. V. 42. P. 1061–1095.
Wandt M. A. Е., Underhill L. G. Covariance biplot analysis of trace element concentrations in urinary stones // British Journal of Urology. 1988. V. 61. № 6. Р. 474–481.