Оптимизация процессов кристаллизации биоматериалов тема автореферата и диссертации по физике, 01.04.07 ВАК РФ
Безбах, Илья Жанович
АВТОР
|
||||
кандидата физико-математических наук
УЧЕНАЯ СТЕПЕНЬ
|
||||
Калуга
МЕСТО ЗАЩИТЫ
|
||||
2007
ГОД ЗАЩИТЫ
|
|
01.04.07
КОД ВАК РФ
|
||
|
На правах рукописи УДК 548 5 532 5
БЕЗБАХ Илья Жанович
ОПТИМИЗАЦИЯ ПРОЦЕССОВ КРИСТАЛЛИЗАЦИИ БИОМАТЕРИАЛОВ
Специальность 01 04 07 - Физика конденсированного состояния
АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата физико-математических наук
Москва-2007
003174179
Работа выполнена в государственном образовательном учреждении высшего профессионального образования "Московский государственный технический университет им Н Э Баумана" и в Научно-исследовательском центре "Космическое материаловедение" Института кристаллографии им А В Шубникова Российской академии наук
Научный руководитель
доктор физико-математических наук Стрелов Владимир Иванович
Официальные оппоненты
доктор технических наук, профессор Андреев Владимир Викторович кандидат физико-математических наук Плотников Федор Алексеевич
Ведущая организация
Филиал ФГУП "Научно-исследовательский физико-химический институт им Л Я Карпова"
Защита состоится «31» октября 2007 г. в 16.00 часов на заседании диссертационного совета Д 212141 17 при ГОУ ВПО "Московский государственный технический университет им Н Э Баумана" по адресу 248600, г Калуга, ул Баженова, д 2, МГТУим НЭ Баумана, Калужский филиал
Ваши отзывы в двух экземплярах, заверенные печатью, просьба высылать по указанному адресу
С диссертацией можно ознакомиться в библиотеке ГОУ ВПО МГТУ им Н Э Баумана, Калужский филиал (г Калуга, ул. Баженова, д 2)
Автореферат разослан <<^у> с-ектЯиТрд 2007 г
Ученый секретарь диссертационного совета канд техн. наук, доцент
С/ / Лоскутов С А
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Актуальность работы Кристаллизация биоматериалов в настоящее время используется для удовлетворения растущей потребности структурной биологии для определения кристаллографических параметров органических молекул, что в дальнейшем позволяет проводить как синтез новых веществ с требуемыми свойствами, так и решать некоторые фундаментальные вопросы функционирования живых систем в целом Важнейшим фактором, ограничивающим эти исследования, являются процессы роста биокристаллов, оптимизация условий которых в настоящее время проводится, в основном, эмпирическими методами
В связи с этим, изучение закономерностей роста биокристаллов и разработка новых методов их получения являются актуальной научной задачей, что подтверждается и ростом числа публикаций в этой области исследований
Кристаллы белков выращиваются из растворов путем доведения раствора до определенных значений пересыщения Пересыщение сг определяется как отношение текущей концентрации белка к его растворимости В общем случае пересыщение может достигаться путем изменения концентрации оса-дителя, белка и нейтральных добавок, рН, температуры и других параметров
Как и в случае неорганических соединений, при достижении пересыщения система, содержащая белок, движется к состоянию равновесия, при котором вещество распределяется между раствором и твердой фазой
Рост кристаллов белков можно разделить на три существенных этапа зародышеобразование, рост, достижение равновесия и прекращение роста Кристаллы белков зарождаются при высоком уровне пересыщения, от 3 до 10, растут же при а> 1, а при а< 1 происходит растворение твердой фазы Кристаллы белков растут очень медленно - их скорость роста на 2-3 порядка меньше, чем у неорганических кристаллов
В настоящее время для получения биокристаллов используются такие методы, как объемная кристаллизация, микродиализ, свободная диффузия через поверхность раздела белка и осадителя, метод диффузии паров Одним из важных факторов, объединяющих эти процессы, является тепломассопе-ренос (ТМП), определяющий качество получаемых кристаллов Изучение особенностей и оптимизация процессов ТМП при кристаллизации белковых молекул позволяет определить лимитирующие условия процессов получения совершенных кристаллов, разработать методы управления процессом кристаллизации В связи с этим, целью работы являлась оптимизация процесса выращивания в земных и космических условиях высокосовершенных кристаллов белков размером не менее 100 мкм на основе математического моделирования и экспериментального исследования процессов ТМП в растворе белка, а также разработка метода и аппаратуры для температурного управления ростом кристаллов на примере модельного белка лизоцима
Объекты и методы исследования Объектом исследований были выбраны кристаллы белка (фермента) лизоцима — широко известного представителя биоматериалов, структура и свойства которого достаточно исследованы, что позволяло проводить сравнительный анализ результатов разработанного метода кристаллизации В качестве методов исследования использовались математическое моделирование и экспериментальные исследования особенностей кристаллизации биоматериалов (в т ч для определения граничных условий), разработка экспериментальной аппаратуры для получения высокосовершенных кристаллов белков, а также методы рентгеноструктурного анализа для изучения совершенства структуры полученных кристаллов
Задачи исследования
1 Определить основные факторы, влияющие на процесс кристаллизации белков
2 Установить закономерности процессов кристаллизации биоматериалов на примере кристаллизации белка лизоцима для земных и космических условий, в зависимости от интенсивности конвективных течений (в т ч конвекции Марангони) и внешних энергетических воздействий на раствор белка Определить условия преобладания диффузионного ТМП
3 Определить оптимальные условия и разработать метод выращивания высокосовершенных кристаллов белков, пригодных для рентгеноструктурного анализа (размером не менее 100 мкм)
4 Разработать методику и алгоритм (программу) управления процессом зародышеобразования и роста биокристаллов
Научная новизна В ходе выполнения диссертационной работы впервые получены следующие результаты
1 Установлены закономерности роста кристаллов биоматериалов при задании и прецизионном управлении температурой всего раствора белка и в локальной точке капилляра для обеспечения образования единичных центров кристаллизации
2 Разработана математическая модель управления процессом кристаллизации, позволяющая определить зависимость количества и размеров образующихся кристаллов от величины пересыщения в растворе, величины критического пересыщения раствора и скорости кристаллизации
3 Разработан алгоритм и оригинальная методика управления температурой как всего объема раствора белка, так и в локальной точке при кристаллизации биоматериалов
4 На основе теоретических расчетов процессов ТМП в растворах белков разработан метод и экспериментальная установка для получения высокосовершенных кристаллов биоматериалов, в которой осуществляется не только управление температурой с точностью регулирования не хуже ±0,1 °С, но и реализована возможность создания пересыщения раствора в заданном месте капилляра за счет точечного подвода или отвода тепла (т е путем создания локального градиента температуры)
Практическая значимость работы состоит в разработке метода управляемой кристаллизации биоматериалов, обеспечивающего раздельное управление процессом как на этапе зародышеобразования, так и в процессе роста образовавшихся кристаллов Этот способ температурного управления процессом кристаллизации биоматериалов является технологичным и эффективным для выращивания высокосовершенных кристаллов, при этом конвекция в растворе практически исключается, а влияние вибраций на процесс кристаллизации минимизируется и, таким образом, обеспечиваются преимущественно диффузионные условия массопереноса белка, т е условия стабильности для самоорганизации молекул белка и встраивания их в кристалл
На защиту выносятся
1 Результаты математического моделирования процесса взаимной диффузии молекул белка и осадителя для реальных условий кристаллизации модельного белка лизоцима
2 Результаты исследования влияния диффузионных и конвективных процессов на кристаллизацию белков в земных и космических условиях с учетом термогравитационной и термокапиллярной (Марангони) конвекции
3 Математическая модель процессов ТМП в растворе белка для управления процессом кристаллизации при создании в растворе локального градиента температуры
4 Разработанный метод выращивания высокосовершенных кристаллов белков в градиенте температуры в ячейках капиллярного типа
Апробация работы Основные результаты, представленные в диссертации, были доложены и обсуждены на следующих конференциях
1 Всероссийская научно-техническая конференция "Прогрессивные технологии, конструкции и системы в приборо- и машиностроении " (Калуга, 2004 г )
2 XIV Международная конференция по росту кристаллов "Fourteenth International Conference on Crystal Growth" (ICCG-14) (Гренобль, Франция, 2004 г )
3 VI Международная конференция "Рост монокристаллов и тепломассоперенос" (ICSC-2005) (Обнинск, 2005 г )
4 XI и XII Национальные конференции по росту кристаллов (Москва, 2004, 2006 гг )
5 Российский симпозиум "Космическое материаловедение" (КМ-2007) (Калуга, 2007 г )
Результаты работы использованы в НИЦ "Космическое материаловедение" ИК РАН (г Калуга) для экспериментальной кристаллизации модельного белка лизоцима с высоким совершенством структуры при прецизионном управлении градиентом температуры в ячейках капиллярного типа
Публикации По теме диссертации опубликовано 12 печатных работ
Объем работы Диссертационная работа включает шесть глав, выводы, список литературы и три приложения Она изложена на 126 страницах текста, содержит 23 рисунка, 2 таблицы, 71 библиографических названий
Благодарности Автор выражает искреннюю благодарность доктору технических наук, профессору В Г Косушкину (КФ МГТУ им НЭ Баумана, г Калуга) за неоценимую поддержку при подготовке диссертации, ценные советы и замечания, также благодарит коллектив Лаборатории белковой кристаллографии и лично И П Куранову (ИК РАН, г Москва) за предоставление кристаллизационных растворов и проведение рентгеноструктурных исследований
СОДЕРЖАНИЕ РАБОТЫ Во введении обосновывается актуальность задачи, указываются основные цели исследования и положения, выдвигаемые на защиту
Первая глава посвящена анализу литературных данных по теме работы В ней иллюстрируются особенности структуры биоматериалов, указывается проблема сложности кристаллизации, приводится обзор применяемых методов выращивания биокристаллов в наземных и космических условиях На примере изучения особенностей оборудования, применявшегося для выращивания биокристаллов в космосе, иллюстрируется низкая стабильность тепловых условий и невозможность динамического управления процессом кристаллизации в большинстве из описываемой кристаллизационной аппаратуры Анализируются существующие методы математического моделирования процессов роста биокристаллов, указываются особенности различных подходов к моделированию моделирование диффузионного массопереноса в кристаллизационных растворах и средах с поверхностью раздела газ-жидкость, моделирование диффузионно-конвективного тепломассопереноса, моделирование зародышеобразования, моделирование кинетики роста биокристаллов, моделирование температурно-зависимых процессов кристаллизации биоматериалов Обращается внимание на отсутствие общепринятой модели, описывающей рост биокристаллов На основе анализа литературных данных формулируются задачи исследования
Во второй главе диссертации проводится исследование процесса взаимной диффузии для молекул белка и осадителя в условиях реальных экспериментов по выращиванию кристаллов белков на примере модельного протеина лизоцима Моделировались изменения во времени концентраций соли, белка и пересыщения, результаты сравнивались с данными практических экспериментов, выполненных по двум методикам 1) диализа - камеры с растворами соли и белка разделены мембраной, проницаемой для молекул соли и непроницаемой для молекул белка, 2) свободной диффузии - мембрана между камерами с растворами соли и белка отсутствует Нестационарный перенос массы, количества движения и плотности в растворе описывался системой уравнений Навье-Стокса в безразмерной форме и приближении Буссинеска при постоянной температуре
Приводятся результаты одномерного расчета динамики распределения пересыщения раствора белка в кристаллизационной камере, демонстрирующие хорошее согласие полученных результатов с результатами практических экспериментов Максимальное различие в количественных значениях относительных пересыщений составляет ~4%
Отмечается, что описанная модель позволяет в первом приближении решать основные задачи моделирования процесса диффузионного массопе-реноса при росте биокристаллов в земных и космических условиях Однако исследовать особенности процессов тепло- и массопереноса для кристаллизационных растворов в земных и космических условиях требуется с учетом влияния на них конвективных и вибрационных воздействий
Третья глава диссертации посвящена исследованию влияния конвективных процессов на процесс кристаллизации для случая нестационарного режима ТМП Решается задача определения характера (структуры и скоростей) конвективных потоков в зависимости от интенсивностей термогравитационной и термокапиллярной (Марангони) конвекций при выращивании кристаллов биоматериалов в земных и космических условиях на примере модельного белка лизоцима
Учитывалось, что конвекция Марангони является отрицательным фактором при кристаллизации методом диффузии паров, поскольку в этом методе присутствует свободная поверхность (жидкость-газ) Также и в методе жидкостной диффузии возможно присутствие открытой поверхности раствора
Одним из основных факторов, влияющих на процесс кристаллизации, является стабильность поддержания температуры кристаллизационного процесса, которая на практике, как правило, не превышает ±1°С
Для описания процессов, происходящих в растворах при выращивании кристаллов биоматериалов, использовалась математическая модель на основе системы уравнений Навье-Стокса в приближении Буссинеска Расчетная область имела цилиндрическую форму с осью симметрии на левой ее границе Расчет проводился для экспериментально используемой концентрации белка лизоцима в водном растворе - 50 мг/мл Перепад температур задавался на верхней границе расчетной области
Проведенные расчеты показали, что при отсутствии конвекции Марангони в земных условиях реализуется режим термогравитационной конвекции, конкурирующей с диффузионным массопереносом (рис 1) При этом конвективные потоки развиваются до половины кристаллизационного объема по высоте Скорость конвективных потоков в верхней части составляет ~10~5 мм/с Эти расчеты подтверждают общепринятое представление о доминировании конвекционного переноса над диффузионным при единичной гравитации В космических же условиях, как показывают расчеты, термогравитационная конвекция значительно ослабляется - скорость конвективных потоков уменьшается до ~10~9 мм/с (рис 2) Диффузионный массоперенос в таких условиях является преобладающим
I. Мп = 0; К-КГ5, мм/с II. МпфО-У, мм/с
Рис. 1. Структура (а) и изолинии (б) скорости V конвективных потоков в растворе белка при наличии (II) и отсутствии (I) конвекции Марангони в земных условиях, ДГ= 1°С
Рис. 2. Структура (а) и изолинии (б) скорости V конвективных потоков в растворе белка при наличии (II) и отсутствии (I) конвекции Марангони в космических условиях, АТ = 1°С
Совсем другая картина наблюдается при наличии конвекции Марангони В этом случае даже при небольшом перепаде температур на открытой поверхности раствора существенным образом меняются структура и интенсивность конвективных течений для значения АТ= 1°С скорость конвективных потоков существенно возрастает и составляет ~1 мм/с как для земных, так и для космических условий, однако снижаясь до -КГ1 мм/с лишь при понижении перепада температуры до ДГ=0,1оС В этих случаях вблизи свободной поверхности раствора наблюдается преобладание конвективного массопереноса
Для широко применяемого метода лежачей капли проведенные расчеты показали, что при наличии перепада температур (который всегда возможен в этом методе, что вызывает конвекцию Марангони) как для земных, так и для космических условий реализуется режим доминирования конвекции Марангони Конвективные потоки развиваются по высоте всей капли, придавая кристаллизационному раствору вращательное движение и приводя тем самым к нарушению диффузионного режима массопереноса Для обычных значений ДТ= 1°С их скорость имеет порядок ~1 мм/с и не зависит от размера капли При понижении перепада температуры до ДГ = 0,1°С скорость снижается до -КГ1 мм/с, те имеет пропорциональную зависимость от градиента температур, как и для случая объемного раствора Проведенные расчеты показывают, что для метода лежачей капли конвективный массоперенос в объеме раствора всегда будет преобладать над диффузионным
Проведенные теоретические расчеты показали существенное влияние не-стабильностей температуры и возникающей конвекции Марангони на процессы ТМП при выращивании биоматериалов как в земных, так и в космических условиях При этом показано, что причиной их возникновения является наличие температурного градиента на свободной поверхности раствора Полученные результаты доказывают необходимость достижения высокой стабильности тепловых условий при выращивании кристаллов биоматериалов
Четвертая глава диссертации посвящена исследованиям влияния вибраций на процесс кристаллизации белков, так как энергетические воздействия существенно влияют на стабильность процесса кристаллизации. При этом исследования влияния вибрационных воздействий на процесс кристаллизации биоматериалов основаны на предварительных расчетно-экспериментальных исследованиях влияния вибрационных воздействий нарост неорганических кристаллов Установлено, что при выращивании кристаллов полупроводников (на примере ве, легированного ва) поступательные колебания, осуществляемые вдоль градиента температуры (по нормали к фронту кристаллизации), в широком диапазоне частот (0,5 - 200 Гц) и амплитуд микроускорений (до 2 10"' Що) не оказывают влияния на процессы тепломассопереноса в расплаве и, соответственно, на однородность распределения легирующей примеси в кристалле, что согласуется с проведенными теоретическими расчетами Однако к существенным нарушениям в распределении примеси приводят торсионные колебания (осуществляемые вокруг оси выращиваемого кристалла) при наличии свободной поверхности расплава
При открытой поверхности расплава (присутствие конвекции Маран-гони) к образованию в кристаллах микронеоднородности (полосчатости) приводят в области частот (0,5 — 200 Гц) колебания с уровнем амплитуд > 10~2§о При амплитудах < 5 10~3»о в выращенных кристаллах полосы роста отсутствуют Очевидно, что возникновение полос роста связано с тем, что при достижении критического уровня вибраций в непосредственной близи от фронта кристаллизации возникают сдвиговые смещения расплава, которые изменяют стационарный характер конвективных течений
Влияние вибраций на процесс кристаллизации биоматериалов исследован для двух случаев вибрационных воздействий, характерных для условий роста белков в космосе на Международной космической станции (МКС) (£/§о> Ю-6) вибрационные возмущения, перпендикулярные градиентам концентрации белка и соли, и параллельные к градиентам концентрации
Расчеты для земных условий показали, что если эксперимент проводить в гелевой среде, то вибрации не оказывают влияния на характер массоперено-са в растворе, и процесс идет в режиме диффузионного массопереноса
В условиях же микрогравитации действуют периодические, высокочастотные возмущения (10~6 < ^ёо < Ю-2, 0,1 -300 Гц по данным цитируемой работы), практически определенные при реализации экспериментов по кристаллизации модельного белка лизоцима на МКС
Расчеты для первого случая показали, что при кристаллизации из-за образования градиентов концентрации соли и белка в белковой камере при вибрационных возмущениях, перпендикулярных градиентам концентрации, возникают конвекционные участки Для второго случая установлено, что при вибрационных возмущениях, параллельных градиентам концентрации, возникновение таких участков строго ограничено зонами, в которых происходит зародышеобразование или осаждение образовавшихся кристаллов Область, подверженная конвекции при наличии таких вибрационных возмущений, перемещается с фронтом зародышеобразования, тогда как в оставшейся части белковой камеры конвекции нет, поэтому в ней следует ожидать улучшения качества кристаллов
В обоих случаях (при росте неорганических и органических кристаллов) установлено существенное влияние вибрационных воздействий на процессы массопереноса и, соответственно, совершенство получаемых кристаллов в определенном диапазоне частот и амплитуд микроускорений Однако при этом влияние вибрационных воздействий на процесс кристаллизации биоматериалов носит более сложный характер
В пятой главе диссертации исследована возможность и показано преимущество использования температуры как средства управления кристаллизацией белков и представлен разработанный метод управления процессом кристаллизации Показано, что в качестве регулирующих параметров могут использоваться концентрация и разные типы осадителей, внесение или изначальное присутствие нейтральных добавок (буферы, гели, масла), затравливание (в т ч материалами с особыми свойствами), наложение электрических
и магнитных полей Однако для практической реализации управляемого выращивания биокристаллов необходимо было выбрать пригодный для физического моделирования параметр Исследована возможность раздельного управления кристаллизацией белка в капилляре с помощью температуры отвод тепла в одной точке раствора (71,) с термостатированием всего объема (Г2) (рис 3)
Тг
Г
Рис 3 Схематическое изображение кристаллизационной ячейки с капилляром
В этом методе становится возможным разделить зародышеобразование и непосредственно рост путем повышения пересыщения <т = Cв/S'(7,) (Ся -концентрация белка) в конкретном месте капилляра для стимулирования заро-дышеобразования (сг> 3) и, после образования единичных центров кристаллизации, понижения пересыщения в этом месте до оптимального для роста (<т>1)
Для описания температурно-управляемого процесса кристаллизации белка была разработана математическая модель, включающая в себя уравнения неразрывности, количества движения и переноса концентрации белка, решаемые в естественных переменных Для описания встраивания белковых молекул в кристалл потребовалась специфическая постановка граничных условий - сток массы белка, вызываемый только изменением Т Отметим, что влияние Г осуществляется через зависимость кривой растворимости на левой границе, зародышеобразование на которой наиболее вероятно
Система уравнений имеет следующий вид
(ди ЗиЛ др д2и о\ — + и— =—— + М—2 \ч дх) дх дх
V Ы дх ) дх дх
ЁР + дРЦ_ = о
д( дх
где р - плотность раствора, и - скорость, / - время, х - координата, р — давление, /л — динамическая вязкость белка, р* = рб/р - концентрация белка, ря - плотность белка, О - коэффициент диффузии белка Граничные условия При х - 0 ставятся следующие условия
для скорости и =
О, если рб<р„ас\
--■^,есшрб>р11ас(Т);
р дх
если
Рб^РнаЛТ)'Т0 Рб=Р„ас(Т)>
дРб
иначе = 0;
дх
(2)
(3)
р = ре+рв; (4) где икр — скорость вывода массы белка в кристалл, ре - плотность воды;
ршс(т) — плотность белка, соответствующая насыщению при данной температуре; Т = Т^ (7) — температура на левом конце, значение может поддерживаться или меняться по линейному закону. При х = Ь и = 0 и —= 0; Т = Т2.
дх
В расчетах использовались реальные размеры капилляра: ¿ = 30 мм, й= 1,5 мм (в таких капиллярах проводился процесс кристаллизации). В связи с наличием преобладающего размера (¿/¿/ = 20) была рассмотрена одномерная модель гидродинамики.
Начальные условия. При 1 = 0 концентрация белка в растворе р6 =50 мг/мл, г/ = 0. Используемая для вычислений зависимость растворимости модельного белка лизоцима от Т показана на рис. 4.
5, мг/мл 60
Т, °С
Рис. 4. График аппроксимированной зависимости растворимости лизоцима, использованный для вычислений
Проводились расчеты процессов ТМП при кристаллизации белка лизоцима с наложением градиента температуры на ячейку капиллярного типа (рис. 3). Предварительно методом математического моделирования процессов ТМП было установлено, что при таких температурных условиях тепловая конвекция в капилляре отсутствует.
т2, °с
т;, °с
т, мг-10"
9 10 11 12 13 14 15 16 / час
Т. = Г,, °с
2 3 4 5 6 7
т, мг-Ю
9 10 11 12 13 14 15 16 I, час
Т2, °с т{, °с
т, мг-10~
1 2 3 4 5
10 11 12 13 14 15 16 час
Рис. 5. Временные зависимости выведенной в кристалл массы белка (т, мг, ■) и создаваемого при этом в растворе пересыщения (сг, *) от значений температуры в точке отвода тепла (7|, °С, ♦) и в объеме (Т2, °С, *)
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 Г, Час
Рис 6 Временные зависимости создаваемого в растворе пересыщения ег пунктирная линия (а) - неуправляемый эксперимент с разностью температур, точечная линия (б) - эксперимент без разности температур, сплошная линия (в) - эксперимент с управлением температурой
Расчетная область описывает только жидкую фазу, процесс кристаллизации характеризуется значением пересыщения на границе и количеством выведенной в кристалл массы белка Сетка из 60 точек, равномерная Шаг по времени 1 с
В результате проведенных расчетов были получены временные зависимости выведенной в кристалл массы белка и создаваемого при этом в растворе пересыщения (для характеризации процесса кристаллизации). Меняющимся параметром в расчетах являлась температура в точке отвода тепла Тх при стабилизации температуры объема раствора Т2 В первом случае она поддерживалась постоянной и достаточной для зародышеобразования и кристаллизации, вызываемой разностью температур (Тх = 19°С), во втором - равной температуре объема раствора, что недостаточно для зародышеобразования и кристаллизации (7^ = Тг = 32°С), в третьем же - менялась по линейному закону со скоростью 1,5°С/час (0,025°С/мин), начиная с расчетного момента зародышеобразования (<т>3) Отметим, что тк требуется снижение пересыщения до сг>1, необходимо повышение растворимости, и, соответственно, температуры
На рис 5 представлены результаты расчетов для трех различных типов экспериментов, а рис 6 резюмирует расчетные значения пересыщения для этих случаев Проведенные расчеты показывают, что в первом случае (рис 5а), (рис 6а) кристаллизация имеет место, но пересыщение возрастает до ег » 5, что является избыточным и будет приводить к массовому зародышеоб-разованию Во втором случае (рис 56), (рис 66) пересыщение не достигает значения ег > 1, и зародышеобразования не происходит В третьем же случае (рис 5в), (рис 6в) управление температурой по заданному закону на определенном этапе, после образования единичных зародышей, позволяет предот-
вратить дальнейшее возрастание пересыщения и плавно понизить его до значения а> 1, необходимого для роста образовавшихся зародышей При этом показано, что расчетное значение выведенной в кристалл массы белка /и» 0,18 мг соответствовало размерам (0,8 мм х 0,9 мм х 0,5 мм) образовавшегося в реальных экспериментах кристалла лизоцима При таком варианте проведения эксперимента (с введением теплового управления) следует ожидать улучшения совершенства образующегося кристалла
Полученные результаты теоретически доказали возможность управления с помощью температуры процессом кристаллизации биоматериалов в земных и, соответственно, в космических условиях при создании локального градиента температуры в растворе На их основе была поставлена и реализована задача создания кристаллизационной ячейки для выращивания в капиллярах кристаллов биоматериалов с помощью управления температурой
В шестой главе диссертации представлены результаты разработки метода и кристаллизационной аппаратуры для выращивания высокосовершенных кристаллов биоматериалов Рассматривается методика получения биокристаллов с управляемым процессом кристаллизации Приводятся результаты экспериментов по кристаллизации модельного белка лизоцима
Описывается конструкция экспериментальной аппаратуры, которая применяется для температурного управления процессом кристаллизации Она способна поддерживать и управлять температурой в пределах 4-40°С с точностью ±0,1 °С Запечатанный капилляр, содержащий кристаллизационный раствор, находится в одной части термостатированного корпуса, а в другой части располагается управляющая электроника Капилляр в определенном месте контактирует с термоуправляемым конусом, создающим перепад температуры в заданной точке Управляющими параметрами в процессе кристаллизации являются две температуры всего корпуса-термостата и конуса, которые могут быть установлены в заданных пределах с шагом 0,1 °С При практических экспериментах корпус работает как термостат, поддерживая температуру в капилляре с общим объемом раствора 10-20 мкл с точностью ±0,1 °С, а конус исполняет роль источника или стока тепла в зависимости от режима работы Процесс зародышеобразования и кристаллизации наблюдается через микроскоп, для чего объем раствора подсвечивается свето-диодом Это позволяет на определенном этапе обнаружить появление единичных зародышей вблизи точки контакта и своевременно начать изменение температуры конуса (и, соответственно, пересыщения в растворе), предотвращая дальнейшее массовое зародышеобразование
Лизоцим из белка куриных яиц (производитель Sigma) кристаллизовался из водного раствора (концентрация 50 мг/мл) с использованием в качестве осадителя соли NaCl (концентрация 8%, рН = 4,6) Для оптимизации экспериментов по кристаллизации использовалась известная зависимость растворимости белка от температуры
На рис 7 показаны кристаллы лизоцима, выращенные в капилляре без управления температурой, т е обычным методом свободной диффузии Ис-
ходные водные растворы содержали 50 мг/мл белка и 8 % соли 1ЧаС1 при рН = 4,6 соответственно. Продолжительность процесса кристаллизации -1 сутки; в результате эксперимента было получено большое число кристаллов малого размера (0,1-0,2 мм) и низкого качества, при этом температура корпуса была постоянной и равной 20°С, с точностью поддержания ±0,1 °С.
На рис. 8 показаны кристаллы лизоцима, выращенные в разработанной экспериментальной аппаратуре с применением теплового управления. Капилляр с кристаллами белка сфотографирован через микроскоп. Состав исходного раствора и капилляр те же. Продолжительность процесса до появления кристаллических зародышей - 1 сутки; роста образовавшихся кристаллов - 5 суток. Получено три кристалла большого размера (0,6-0,8 мм). Температура корпуса устанавливалась равной 32°С, с точностью поддержания ±0,1 °С, температура конуса в начале эксперимента составляла 7°С и увеличивалась до 19°С после обнаружения зародышей.
Рис. 7. Фотография кристаллов лизоцима, выращенных без управления температурой методом свободной диффузии
Рис. 8. Фотография кристаллов лизоцима, выращенных в капилляре с применением метода управления температурой
Приводятся результаты рентгеновского исследования выращенных кристаллов. Снятые 75 рентгенограмм позволили собрать набор с полнотой 97,8%
в интервале углов дифракции, позволяющих получение разрешения до 1,5 А Исследованный кристалл лизоцима принадлежал к пространственной группе Р4 с параметрами элементарной ячейки а = 79,131 А, 6= 79,131 А, с = 37,992 А, мозаичность кристалла 0,204, средний /?-фактор по эквивалентным отражениям - 5,63% в зоне разрешения до 1,6 А
Полученные характеристики набора интенсивностей указали на высокое, с современных позиций, качество выращенного кристалла, которое превосходит качество кристаллов, обычно выращиваемых в наземных изотермических экспериментах
Полученные результаты позволяют сделать вывод, что управление температурой дает возможность не только управлять процессом на этапах за-родышеобразования и дальнейшей кристаллизации, но и создать оптимальные условия для выращивания высококачественных кристаллов белков
В 3-х приложениях приводятся материалы, иллюстрирующие программную реализацию основных аспектов диссертационной работы (алгоритм расчета взаимной диффузии белка и осадителя в практических кристаллизационных экспериментах, алгоритм расчета процесса кристаллизации лизоцима с управлением температурой, реализованный в математическом пакете Matlab на встроенном языке программирования, алгоритм обмена данными и командами между управляющим персональным компьютером и разработанной кристаллизационной аппаратурой, реализованный в среде разработки Delphi)
ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ
1 Установлены закономерности процессов тепломассопереноса при росте в капиллярах кристаллов биоматериалов, в том числе при действии возмущающих воздействий (конвективных и вибрационных), определяющие условия преобладания диффузионного массопереноса
2 Разработана и программно реализована математическая модель, описывающая процесс зарождения и роста кристаллов белка из раствора Математическая модель описывает образование зародышей кристаллов и их рост в зависимости от локального значения перенасыщения
3 С помощью предложенной модели определены оптимальные условия (на примере выращивания кристаллов лизоцима из водного раствора) кристаллизации белков размером не менее 100 мкм стабильность поддержания температуры должна быть не хуже ±0,1 °С, локальный градиент температуры - не более 1 °С/см Созданная модель вполне адекватна моделируемому явлению и может быть использована для параметрических исследований и прогнозных расчетов процесса кристаллизации белков при переменных температурных полях
4 На основе проведенных теоретических расчетов исследована и показана возможность кристаллизации в капиллярах белков с высоким совершенством их структуры, недостижимом при использовании традиционно используемых в настоящее время методов кристаллизации, заключающаяся в реализации метода управляемой кристаллизации белков, обеспечиваю-
щего раздельное управление процессом роста кристаллов как на этапе их зародышеобразования, так и в процессе кристаллизации
5 Показано, что управление процессом зарождения и роста может осуществляться путем а) задания и прецизионного поддержания требуемой температуры всего раствора белка в капилляре, б) поддержания в локальной точке капилляра с раствором соответствующей температуры для обеспечения необходимого пересыщения для зарождения единичных (1-2) центров кристаллизации, в) изменения температуры в локальной точке капилляра и, соответственно, пересыщения в процессе разращивания кристалла
6 Разработан опытный лабораторный образец экспериментального блока-кристаллизатора биокристаллов В этой системе осуществляется не только управление температурой раствора с точностью регулирования ±0,1°С, но и реализована возможность создания пересыщения белка в соответствующем месте капилляра за счет точечного подвода или отвода в этом месте тепла (т е создания в этом месте локального температурного градиента), что исключает возможность образования множества зародышей по всему объему раствора, находящегося в контролируемых изотермических условиях, что исключает и многоцентровую спонтанную кристаллизацию
ОСНОВНОЕ СОДЕРЖАНИЕ ДИССЕРТАЦИИ ОТРАЖЕНО В СЛЕДУЮЩИХ РАБОТАХ
1 Математическое моделирование и оптимизация роста биокристаллов / M Ж Безбах, В1 Косушкин, В M Стрелов и др // Прогрессивные технологии, конструкции и системы в приборо- и машиностроении Материалы Всероссийской научно-технической конференции - M, 2004 - Т 1 - С 277-278
2 Numerical models for biocrystal growing ! VP Gmkin В G Zakharov, IZh Bezbakh et al // Abstracts of the 14th International Conference on Crystal Growth - Grenoble (France), 2004 - P 335
3 Управление процессами кристаллизации биоматериалов в земных и космических условиях ! И Ж Безбах, Б Г Захаров В И Стрелов ыдр //Тезисы докладов XI Национальной конференции по росту кристаллов - M, 2004 ~ С 355
4 Математическое моделирование и экспериментальные исследования влияния вибраций и конвекции Марангони на микрооднородность кристаллов полупроводников / В И Стрелов, Б Г Захаров, И Ж Безбах и др //Поверхность -2005 -№10 - С 80-87
5 Моделирование процесса управления ростом кристаллов белков в земных и космических условиях / И Ж Безбах, В И Стрелов, В Г Косушкин и др // Рост монокристаллов и тепломассоперенос Труды 6-й международной конференции - Обнинск, 2005 -T3 -С 545-551
6 Математическое моделирование и экспериментальные исследования влияния температурных градиентов на процессы кристаллизации для земных и космических условий / В И Стрелов, Б Г Захаров, И Ж Безбах и др //Кристаллография -2005 -Т50,№3 -С 536-544
I Mathematical modeling and expenmental investigation of the effect of temperature gradients on crystallization processes under terrestrial and space conditions /VI Strelov, В G Zakharov, IZh Bezbakh et al II Crystallography Reports -2005 - V 50, No 3 -P 490-498
8 Прецизионное температурное управление процессом кристаллизации биоматериалов / И Ж Безбах, В И Стрелов, H И Сосфетв и др // Тезисы докладов XII Национальной конференции по росту кристаллов - M, 2006 -С 330
9 Оптимизация роста кристаллов белков с применением метода теплового управления ! И Ж Безбах, В Г Косушкин, В И Стрелов и др // Методы исследования и проектирования сложных технических систем Сборник статей (Труды МГТУ №592) -М,2006 -С 18-26
10 Кристаллизация белка лизоцима в прецизионно-управляемом градиенте температуры / В И Стрелов, Б Г Захаров, И Ж Безбах и др //Кристаллография -2007 -Т52,№6 -С 1134-1139
II ЬезбахИЖ, Стрелов ВИ, ЗахаровБГ Управление процессом зародышеобразования и роста биокрисггал-лов // Тезисы докладов Российского симпозиума "Космическое материаловедение" - Калуга, 2007 - С 38
12 Моделирование роста биокристаллов под воздействием управляющего теплового поля / В П Гшкин И Ж Безбах, В И Стрелов и др П Тезисы докладов Российского симпозиума "Космическое материаловедение" -Калуга, 2007 - С43
Безбах Илья Жанович
Оптимизация процессов кристаллизации биоматериалов
Автореферат диссертации на соискание ученой степени кандидата физико-математических наук
Подписано в печать 25 09 2007 Бумага офсетная 80г/м2 Формат 60x84 1/16 Уел печ л 1 Уч -изд л 1 Тираж 100 экз Заказ № 176
Отпечатано в Редакционно-издательском отделе
КФМГТУим НЭ Баумана 248000, г Калуга, ул Баженова, 2, тел 57-31-87
СПИСОК ИСПОЛЬЗУЕМЫХ СОКРАЩЕНИЙ.
Введение.
Глава 1. Анализ литературы.
1.1. Общая характеристика проблемы кристаллизации биоматериалов.
1.2. Обзор методов выращивания биокристаллов, применяемых в наземных и космических условиях.
1.3. Обзор средств управления кристаллизацией.
1.4. Выводы к главе 1.
Глава 2. Исследование процесса взаимной диффузии для молекул белка и осадителя в условиях реальных экспериментов по выращиванию кристаллов белков на примере модельного белка лизоцима.
2.1. Способы описания диффузионных процессов.
2.2. Описание математической модели.
2.2.1. Уравнения тепломассопереноса.
2.2.2. Начальные и граничные условия.
2.2.3. Моделирование осаждения кристаллов.
2.3. Результаты расчета процесса взаимной диффузии молекул белка и осадителя.
2.4. Выводы к главе 2.
Глава 3. Исследование влияния конвективных воздействий на характер процесса кристаллизации для случая нестационарного тепломассопереноса в земных и космических условиях с учетом плотностной, термогравитационной и термокапиллярной (Марангони) конвекций.
3.1. Описание возможных причин возникновения конвективных процессов и методов кристаллизации, подверженных их воздействию.
3.2. Описание примененного подхода к математическому моделированию.
3.3. Результаты расчетов структуры и скоростей конвективных потоков при кристаллизации в объеме раствора.
3.4. Результаты расчетов структуры и скоростей конвективных потоков при кристаллизации методом лежачей капли.
3.5. Выводы к главе 3.
Глава 4. Исследование влияния вибрационных возмущений на процесс кристаллизации.
4.1. Описание особенностей кристаллизации белков в невесомости.
4.2. Влияние вибраций на процесс выращивания кристаллов полупроводниковых материалов.
4.3. Анализ влияния вибраций на процесс кристаллизации биоматериалов.
4.4. Выводы к главе 4.
Глава 5. Разработка метода управления кристаллизацией.
5.1. Анализ существующих методик осуществления управляемой кристаллизации.
5.2. Выбор методики управления кристаллизацией путем создания и регулирования температурного градиента в растворе.
5.3. Создание математической модели для описания вызываемой градиентом температуры кристаллизации биоматериала.
5.4. Результаты расчетов и теоретическая иллюстрация потенциального преимущества температуры для осуществления управляемой кристаллизации в земных и космических условиях.
5.5. Выводы к главе 5.
Глава 6. Разработка экспериментального оборудования для управляемого температурой выращивания высокосовершенных кристаллов биоматериалов.
6.1. Известные предпосылки для осуществления температурно-управляемой кристаллизации биоматериалов.
6.2. Разработка конструкции ростовой аппаратуры для выращивания биокристаллов в условиях температурного градиента в ячейке капиллярного типа при прецизионном управлении температурой.
6.3. Результаты экспериментов по управляемой кристаллизации биоматериалов на примере модельного белка лизоцима.
6.4. Результаты рентгенодифракционных исследований выращенных кристаллов модельного белка лизоцима.
6.5. Выводы к главе 6.
Актуальность работы. Кристаллизация биоматериалов в настоящее время используется для удовлетворения растущей потребности структурной биологии для определения кристаллографических параметров органических молекул, что в дальнейшем позволяет проводить как синтез новых веществ с требуемыми свойствами, так и решать некоторые фундаментальные вопросы функционирования живых систем в целом. Важнейшим фактором, ограничивающим эти исследования, являются процессы роста биокристаллов, оптимизация условий которых в настоящее время проводится, в основном, эмпирическими методами.
В связи с этим, изучение закономерностей роста биокристаллов и разработка новых методов их получения являются актуальной научной задачей, что подтверждается и ростом числа публикаций в этой области исследований.
Кристаллы белков выращиваются из растворов путем доведения раствора до определенных значений пересыщения. Пересыщение <т определяется как отношение текущей концентрации белка к его концентрации в насыщенном растворе. В общем случае пересыщение может достигаться путем изменения концентрации осадителя, белка и нейтральных добавок, рН, температуры и других параметров.
Как и в случае неорганических соединений, при достижении пересыщения система, содержащая белок, движется к состоянию равновесия, при котором вещество распределяется между раствором и твердой фазой.
Рост кристаллов белков можно разделить на три существенных этапа: зародышеобразование, рост, достижение равновесия и прекращение роста. Кристаллы белков зарождаются при высоком уровне пересыщения, от 3 до 10, растут же при сг > 1, а при сг < 1 происходит растворение твердой фазы.
Кристаллы белков растут очень медленно - их скорость роста на 2-3 порядка меньше, чем у неорганических кристаллов.
В настоящее время для получения биокристаллов используются такие методы, как объемная кристаллизация, микродиализ, свободная диффузия через поверхность раздела белка и осадителя, метод диффузии паров. Фактором, объединяющим эти процессы, является тепломассоперенос (ТМП), определяющий качество получаемых кристаллов. Изучение особенностей переноса тепла и массы при кристаллизации белковых молекул позволяет определить лимитирующие условия процессов получения совершенных кристаллов, разработать методы управления кристаллизацией. В связи с этим, целью работы являлась оптимизация условий выращивания совершенных кристаллов белков размером не менее 100 мкм в условиях наземных и космических экспериментов на основе экспериментального изучения и математического моделирования процессов ТМП и разработка метода теплового управления ростом кристаллов на примере модельного белка лизоцима.
Объекты и методы исследования. Объектом исследований были выбраны кристаллы белка (фермента) лизоцима - широко известного представителя биоматериалов, структура и свойства которого достаточно исследованы, что позволяло проводить сравнительный анализ результатов разработанного метода кристаллизации. В качестве методов исследования использовались математическое моделирование и экспериментальные исследования особенностей кристаллизации биоматериалов (в т.ч. для определения граничных условий), разработка экспериментальной аппаратуры для получения высокосовершенных кристаллов белков, а также методы рентгеноструктурного анализа для изучения совершенства структуры полученных кристаллов.
Задачи исследования:
1. Определить основные факторы, влияющие на процесс кристаллизации белков.
2. Установить закономерности процессов кристаллизации биоматериалов на примере кристаллизации белка лизоцима для земных и космических условий, в зависимости от интенсивности конвективных течений (в т.ч. конвекции Марангони) и внешних энергетических воздействий на раствор белка. Определить условия преобладания диффузионного ТМП.
3. Определить оптимальные условия и разработать метод выращивания высокосовершенных кристаллов белков, пригодных для рентгеноструктурного анализа (размером не менее 100 мкм).
4. Разработать методику и алгоритм (программу) управления процессом зародышеобразования и роста биокристаллов.
Научная новизна. В ходе выполнения диссертационной работы впервые получены следующие результаты:
1. Установлены закономерности роста кристаллов биоматериалов при задании и прецизионном управлении температурой всего раствора белка и в локальной точке капилляра для обеспечения образования единичных центров кристаллизации.
2. Разработана математическая модель управления процессом кристаллизации, позволяющая определить зависимость количества и размеров образующихся кристаллов от величины пересыщения в растворе, величины критического пересыщения раствора и скорости кристаллизации.
3. Разработан алгоритм и оригинальная методика управления температурой как всего объема раствора белка, так и в локальной точке при кристаллизации биоматериалов.
4. На основе теоретических расчетов процессов ТМП в растворах белков разработан метод и экспериментальная установка для получения высокосовершенных кристаллов биоматериалов, в которой осуществляется не только управление температурой с точностью регулирования не хуже ±0,1 °С, но и реализована возможность создания пересыщения раствора в заданном месте капилляра за счет точечного подвода или отвода тепла (т.е. путем создания локального градиента температуры).
Практическая значимость работы состоит в разработке метода управляемой кристаллизации биоматериалов, обеспечивающего раздельное управление процессом как на этапе зародышеобразования, так и в процессе роста образовавшихся кристаллов. Этот способ температурного управления процессом кристаллизации биоматериалов является технологичным и эффективным для выращивания высокосовершенных кристаллов, при этом конвекция в растворе практически исключается, а влияние вибраций на процесс кристаллизации минимизируется и, таким образом, обеспечиваются преимущественно диффузионные условия массопереноса белка, т.е. условия стабильности для самоорганизации молекул белка и встраивания их в кристалл.
На защиту выносятся:
1. Результаты математического моделирования процесса взаимной диффузии молекул белка и осадителя для реальных условий кристаллизации модельного белка лизоцима.
2. Результаты исследования влияния диффузионных и конвективных процессов на кристаллизацию белков в земных и космических условиях с учетом термогравитационной и термокапиллярной (Марангони) конвекции.
3. Математическая модель процессов ТМП в растворе белка для управления процессом кристаллизации при создании в растворе локального градиента температуры.
4. Разработанный метод выращивания высокосовершенных кристаллов белков в градиенте температуры в ячейках капиллярного типа.
Апробация работы. Основные результаты, представленные в диссертации, были доложены и обсуждены на следующих конференциях:
1. Всероссийская научно-техническая конференция "Прогрессивные технологии, конструкции и системы в приборо- и машиностроении " (Калуга, 2004 г.).
2. XIV Международная конференция по росту кристаллов "Fourteenth International Conference on Crystal Growth" (ICCG-14) (Гренобль, Франция, 2004 г.)
3. VI Международная конференция "Рост монокристаллов и тепломассоперенос" (ICSC-2005) (Обнинск, 2005 г.).
4. XI и XII Национальные конференции по росту кристаллов (Москва, 2004, 2006 гг.).
5. Российский симпозиум "Космическое материаловедение" (КМ-2007) (Калуга, 2007 г.).
Публикации. По теме диссертации опубликовано 12 печатных работ.
Объем работы. Диссертационная работа включает шесть глав, выводы, список литературы и три приложения. Она изложена на 126 страницах текста, содержит 23 рисунка, 2 таблицы, 71 библиографических названий.
Общие выводы и заключения
1. Установлены закономерности процессов тепломассопереноса при росте в капиллярах кристаллов биоматериалов, в том числе при действии возмущающих воздействий (конвективных и вибрационных), определяющие условия преобладания диффузионного массопереноса.
2. Разработана и программно реализована математическая модель, описывающая процесс зарождения и роста кристаллов белка из раствора. Математическая модель описывает образование зародышей кристаллов и их рост в зависимости от локального значения перенасыщения.
3. С помощью предложенной модели определены оптимальные условия (на примере выращивания кристаллов лизоцима из водного раствора) кристаллизации белков размером не менее 100 мкм: стабильность поддержания температуры должна быть не хуже ±0,1 °С; локальный градиент температуры - не более 1°С/см. Созданная модель вполне адекватна моделируемому явлению и может быть использована для параметрических исследований и прогнозных расчетов процесса кристаллизации белков при переменных температурных полях.
4. На основе проведенных теоретических расчетов исследована и показана возможность кристаллизации в капиллярах белков с высоким совершенством их структуры, недостижимом при использовании традиционно используемых в настоящее время методов кристаллизации, заключающаяся в реализации метода управляемой кристаллизации белков, обеспечивающего раздельное управление процессом роста кристаллов как на этапе их зародышеобразования, так и в процессе кристаллизации.
5. Показано, что управление процессом зарождения и роста может осуществляться путем:
- задания и прецизионного поддержания требуемой температуры всего раствора белка в капилляре;
- поддержания в локальной точке капилляра с раствором соответствующей температуры для обеспечения необходимого пересыщения для зарождения единичных (1-2) центров кристаллизации;
- изменения температуры в локальной точке капилляра и, соответственно, пересыщения в процессе разращивания кристалла.
6. Разработан опытный лабораторный образец экспериментального блока-кристаллизатора биокристаллов. В этой системе осуществляется не только управление температурой раствора с точностью регулирования ±0,1 °С, но и реализована возможность создания пересыщения белка в соответствующем месте капилляра за счет точечного подвода или отвода в этом месте тепла (т.е. создания в этом месте локального температурного градиента), что исключает возможность образования множества зародышей по всему объему раствора, находящегося в контролируемых изотермических условиях, что исключает и многоцентровую спонтанную кристаллизацию.
86
Благодарности
Автор выражает искреннюю благодарность доктору технических наук, профессору В.Г. Косушкину (КФ МГТУ им. Н.Э. Баумана, г. Калуга) за неоценимую поддержку при подготовке диссертации, ценные советы и замечания; также благодарит коллектив Лаборатории белковой кристаллографии и лично И.П. Куранову (ИК РАН, г. Москва) за предоставление кристаллизационных растворов и проведение рентгеноструктурных исследований.
1. Куранова И.П. Кристаллизация белков на земле и в невесомости // Поверхность. — 2004. — №6. — С.4-12.
2. Chayen N.E. Turning protein crystallisation from an art into a science // Current Opinion in Structural Biology. — 2004. — V.14. — P.577-583.
3. McPherson A. Protein crystallization in the structural genomics era // Journal of Structural and Functional Genomics. — 2004. — V.5. — P.3-12.
4. Chayen N.E. Tackling the bottleneck of protein crystallization in the post-genomic era // Trends in Biotechnology. — 2002. — V.20. — P.98-104.
5. Protein crystal growth in space, past and future / LJ.DeLucas, K.M. Moore, MM. Long et al. II Journal of Crystal Growth. — 2002. — V.237. — P.l646-1650.6. www.pdb.org.
6. Вашштейн Б.К., Симонов В.И. Структурная кристаллография. — М.: Наука, 1992. —293 с.
7. Thaimattam R., Jaskolski М. Synchrotron radiation in atomic-resolution studies of protein structure // Journal of Alloys and Compounds. — 2004. — V.362. — P. 12-20.
8. Phenomena and mechanisms of mixed crystal formation in solutions II. Mechanism of interface processes / A.E. Voloshin, S.I. Kovalev, E.B. Rudneva et al. II Journal of Crystal Growth. — 2004. — V.261. — P.105-117.
9. Crystallization of chicken egg white lysozyme from assorted sulfate salts / E.L. Forsythe, E.H. Snell, C.C. Malone et al. II Journal of Crystal Growth. — 1999. — V.l 96. — P.332-343.
10. Galkin O., Vekilov P.G. Nucleation of protein crystals: critical nuclei, phase behavior, and control pathways // Journal of Crystal Growth. — 2001. — V.232. — P.63-76.
11. McPherson A. Virus and protein crystal growth on Earth and in microgravity // Journal of Physics D: Applied Physics. — 1993. — V.26. — P. 104-114.
12. DeLucas L.J. Recent results and new hardware developments for protein crystal growth in microgravity // Journal of Crystal Growth. — 1994. — V.135. — P.183-194.
13. Haas C., Drenth J. Understanding protein crystallization on the basis of the phase diagram // Journal of Crystal Growth. — 1999. — V.196. — P.447-449.
14. Temperature dependence of protein solubility determination and application to crystallization in X-ray capillaries / F. Rosenberger, S.B. Howard, J. W. Sowers et al. II Journal of Crystal Growth. — 1993. — V.129. — P.l-12.
15. Galkin O., Vekilov P.G. Control of protein crystal nucleation around the me-tastable liquid-liquid phase boundary // Proceedings of the National Academy of Sciences, USA. — 2000. — V.97. — P.6277-6281.
16. Durbin S.D., Feher G. Simulation of lysozyme crystal growth by the Monte Carlo method // Journal of Crystal Growth. — 1991. — V. 110. — P.41-51.
17. Chernov A. A. Crystal growth science between the centuries // Journal of Materials Science Materials in Electronics. — 2001. — V.12. — P.437-449.
18. Otalora F., Garcia-Ruiz J.M. Crystal growth studies in microgravity with the APCF. I. Computer simulation of transport dynamics // Journal of Crystal Growth. — 1997. — V. 182. — P. 141-154.
19. Experimental observations and numerical modelling of diffusion-driven crystallisation processes / L. Carotenuto, C. Piccolo, D. Castagnolo et al. II Acta Crystallographica Section D. — 2002. — V.58. — P. 1628-1632.
20. Convective-diffusive transport in protein crystal growth / Я Lin, F. Rosenberger, J.I.D. Alexander et al. II Journal of Crystal Growth. — 1995. — V.151. — P.153-162.
21. Nonlinear dynamics of layer growth and consequences for protein crystal perfection / P.G. Vekilov, F. Rosenberger, H.Lin et al. I I Journal of Crystal Growth. — 1999. —V.196. —P.261-275.
22. Lee C.P., Chernov A.A. Solutal convection around growing protein crystals and diffusional purification in Space // Journal of Crystal Growth. — 2002.1. V.240.—P.531-544.
23. Savino R., Monti R. Convection induced by residual-g and g-jitters in diffusion experiments // International Journal of Heat and Mass Transfer. — 1999.1. V.42. — P.111-126.
24. Savino R., Monti R. Buoyancy and surface-tension-driven convection in a hanging drop protein crystallizer // Journal of Crystal Growth. — 1996. — V.165. — P.308-318.
25. Investigation of Marangoni and natural convection during protein crystal growth / M. Kawaji, O. Gamache, D.H. Hwang et al. II Journal of Crystal Growth. — 2003. — V.258. — P.420-430.
26. Математическое моделирование и экспериментальные исследования влияния температурных градиентов на процессы кристаллизации для земных и космических условий / В.И. Стрелов, Б.Г. Захаров, И.Ж. Безбах и др. II Кристаллография. — 2005. — №3. — С.536-544.
27. Nucleation and crystallization of globular proteins what we know and what is missing / F. Rosenberger, P.G. Vekilov, M. Muschol et al. II Journal of Crystal Growth. — 1996. — V. 168. — P. 1-14.
28. Garcia-Ruiz J.M., Santos A., Alfaro E.J. Oscillatory growth rates in single crystal growing under diffusional control // Journal of Crystal Growth. — 1996. —V.84. —P.555-564.
29. Проблемы и перспективы получения в условиях микрогравитации монокристаллов полупроводников с высокой микрооднородностью свойств / Б.Г. Захаров, П.К. Волков, Ю.А. Серебряков и др. II Поверхность. — 2001.—№9. —С.49-56.
30. Гершуни Г.З., Жуховицкий Е.М. Вибрационная тепловая конвекция в невесомости // Гидромеханика и процессы переноса в невесомости. — Свердловск, 1983. — Р.86-105.
31. The control of Bridgman crystal growth with the help of temperature gradient / V.I. Strelov, VS. Sidorov, B.G. Zakharov et al. II Single crystal growth and heat & mass transfer: Third international conference. — Obninsk, 1999. — P.180-191.
32. Luft J.R., RakD.M., DeTitta G.T. Microbatch macromolecular crystallization on a thermal gradient // Journal of Crystal Growth. — 1999. — V.196. — P.447-449.
33. Vekilov P.G., Rosenberger F. Protein crystal growth under forced solution flow: experimental setup and general response of lysozyme // Journal of Crystal Growth. — 1998. — V. 186. — P.251-261.
34. Saikumar M. V., Glatz C.E., Larson MA. Lysozyme crystal growth and nucleation kinetics // Journal of Crystal Growth. — 1998. — V.187. — P.277-288.
35. Wakayama N.I. Quantitative study of crystallization kinetics of hen egg-white lysozyme using magnetic orientation // Journal of Crystal Growth. — 1998. — V.191. —P.199-205.
36. Why is the osmotic second virial coefficient related to protein crystallization? / B.L. Neal, D. Asthagiri, O.D. Velev et al. II Journal of Crystal Growth. — 1999. —V.196. —P.377-387.
37. Piazza R. Interactions in protein solutions near crystallisation: a colloid physics approach // Journal of Crystal Growth. — 1999. — V.196. — P.415-423.
38. Chernov A.A. Estimates of internal stress and related mosaicity in solution grown crystals: proteins // Journal of Crystal Growth. — 1999. — V.196. — P.524-534.
39. Crystallization of the E. coli polyamine-induced protein: a novel procedure based on the concept of ionic strength reducers / Y. Papanikolau, R. Gessmann, K. Petratos et al. II Journal of Crystal Growth. — 2000. — V.210. — P.761-766.
40. The effect of ordering of internal water in thaumatin and lysozyme crystals as revealed by Raman method / A.B. Kudryavtsev, G. Christopher, S.B. Mirov et al. II Journal of Crystal Growth. — 2000. — V.219. — P.102-114.
41. Perez C., Griebenow K. Fourier-transform infrared spectroscopic investigation of the thermal denaturation of hen egg-white lysozyme dissolved in aqueous buffer and glycerol // Biotechnology Letters. — 2000. — V.22. — P.l 899-1905.
42. Biasutti M.A., Soltermann A.T., Garcia N.A. Photodynamic effect in lysozyme: a kinetic study in different micellar media // Journal of Peptide Research. — 2000. — V.55. — P.41-50.
43. Zhong C., Wakayama N.I. Effect of a high magnetic field on the viscosity of an aqueous solution of protein // Journal of Crystal Growth. — 2001. — V.226. — P.327-332.
44. Qi J., Wakayama N.I., Ataka M. Magnetic suppression of convection in protein crystal growth processes // Journal of Crystal Growth. — 2001. — V.232. — P.132-137.
45. Kulkarni A., Zukoski C. Depletion interactions and protein crystallization //Journal of Crystal Growth. — 2001. — V.232. — P.l 56-164.
46. Thomas B.R., Chernov A.A. Acetylated lysozyme as impurity in lysozyme crystals: constant distribution coefficient // Journal of Crystal Growth. — 2001. — V.232. — P.237-243.
47. Real time evolution of concentration distribution around tetragonal lysozyme crystal: case study in gel and free solution / W.B. Hou, A.B. Kudryavtsev, T.L. Bray et al. II Journal of Crystal Growth. — 2001. — V.232. — P.265-272.
48. Burke M. W., Judge R.A., Pusey M.L. The effect of solution thermal history on chicken egg white lysozyme nucleation // Journal of Crystal Growth. — 2001.1. V.232. — P.301-307.
49. Sumida J.P., Forsythe E.L., Pusey M.L. Preparation and preliminary characterization of crystallizing fluorescent derivatives of chicken egg white lysozyme // Journal of Crystal Growth. — 2001. — V.232. — P.308-316.
50. Influence of impurities on protein crystal perfection / M.C. Robert, B. Capelle, B. Lorber et al. II Journal of Crystal Growth. — 2001. — V.232.1. P.489^97.
51. Wang L., Zhong C., Wakayama N.I. Damping of natural convection in the aqueous protein solutions by the application of high magnetic fields // Journal of Crystal Growth. — 2002. — V.237. — P.312-316.
52. Effects of pressure on the crystallization and the solubility of proteins in agarose gel / A. Kadri, B. Lorber, G. Jenner et al. II Journal of Crystal Growth.2002. — V.245. — P. 109-120.
53. Pusey M.L., Chayen N.E. Well connected: protein-protein interactions and crystallization // Trends in Biotechnology. — 2002. — V.20. — P.322-323.
54. Bakk A., Hoye J.S., Hansen A. Specific heat upon aqueous unfolding of the protein interior: a theoretical approach // Physica A. — 2002. — V.304. — P.355-361.
55. Kundrot C.E. Which strategy for a protein crystallization project? // Cellular and Molecular Life Sciences. — 2004. — V.61. — P.525-536.
56. Savino R., Lappa M. Dissolution and solutal convection in partially miscibile liquid systems I I International Journal of Heat and Mass Transfer. — 2004. — V.47. — P.601-612.
57. Numerical models for biocrystal growing / V.P. Ginkin, B.G. Zakharov, I.Zh. Bezbakh et al II Abstracts of the 14th International Conference on Crystal Growth. Grenoble (France), 2004. - P.335.
58. Управление процессами кристаллизации биоматериалов в земных и космических условиях / И.Ж. Безбах, Б.Г. Захаров, В.И. Стрелов и др. II Тезисы докладов XI Национальной конференции по росту кристаллов. -М., 2004.-С.355.
59. Математическое моделирование и экспериментальные исследования влияния вибраций и конвекции Марангони на микрооднородность кристаллов полупроводников / В.И. Стрелов, Б.Г. Захаров, И.Ж. Безбах и др. И Поверхность. 2005. - №10. - С.80-87.
60. Прецизионное температурное управление процессом кристаллизации биоматериалов / И.Ж. Безбах, В.И. Стрелов, Н.И. Сосфенов и др. II Тезисы докладов XII Национальной конференции по росту кристаллов. М., 2006.-С.330.
61. Оптимизация роста кристаллов белков с применением метода теплового управления / И.Ж. Безбах, В.Г. Косушкин, В.И. Стрелов и др. II Методы исследования и проектирования сложных технических систем: Сборник статей (Труды МГТУ №592). М., 2006. - С.18-26.
62. Кристаллизация белка лизоцима в прецизионно-управляемом градиенте температуры / В.И. Стрелов, Б.Г Захаров, И.Ж. Безбах и др. Н Кристаллография. -2007. -Т.52, №6. С.1134-1139.
63. Безбах И.Ж., Стрелов В.И, Захаров Б.Г. Управление процессом зародышеобразования и роста биокристаллов // Тезисы докладов Российского симпозиума "Космическое материаловедение". Калуга, 2007. - С.38.
64. Моделирование роста биокристаллов под воздействием управляющего теплового поля / В.П. Гинкин, И.Ж. Безбах, В.И. Стрелов и др. II Тезисы докладов Российского симпозиума "Космическое материаловедение". -Калуга, 2007.-С.43.