Лекции Рубина (1123233), страница 21
Текст из файла (страница 21)
В реальных условиях амплитуды изгибных флуктуаций могут достигатьнесколькихангстрем,авременарелаксациилежатвмикросекундном диапазоне.Таким образом, в белке имеется набор связанных элементов, существенно различающихся по временам релаксации. Наиболее быстрые имелкомасштабные флуктуации присущи боковым группам. Эти группыобразуют жидкоподобную опушку вокруг спиральных участков полипептидного каркаса и играют роль демпфирующей среды.
Иерархия во временах релаксаций позволяет представить динамику белковой глобулыкак флуктуации в жидкоподобной капле, армированной упругим полипептидным каркасом. Зависимость скорости реакции в белках от вязкостирастворителяможнопонятьнаосновеэтихпредставлений.Диффузия лигандов внутри глобулы происходит лишь при образованиифлуктуационных полостей или "дырок". Появление "дырки" внутри глобулы может быть инициировано образованием ее вначале в растворителена поверхности глобулы. Вероятность этого процесса обратно пропор-113циона вязкости растворителя.
За счет конформационных движений поверхностная группа белка заполняет "дырку" в растворителе. Тем самым"дырка" теперь оказывается уже в наружном слое белка. Далее, за счетдвижений групп второго и третьих слоев "дырка" диффундирует внутрьглобулы, обеспечивая появление дополнительных флуктуационных полостей. Форма этих полостей в белке, как мы видели, не произвольна, аимеет вид флуктуирующих щелей, параметры которых определяютсягеометрией жестких элементов белкового каркаса. Решение диффузионных уравнений позволяет вычислить скорость диффузии частиц черезсистему таких флуктуирующих щелей. Так, при диффузии в миоглобинелиганд СО должен пройти несколько "ворот", которые открываются засчет конформационных движений.
Скорость диффузии зависит от диаметра лиганда, амплитуды флуктуации и времени релаксации щели,котораявсвоюочередьопределяетсяжесткостьюимикровязкостьюстенок. Для миоглобина расчетные и экспериментальные данные совпали. В частности, оказалось, что общее время прохождения СО в миоглобинесоставляет710-сисоответствуетсложениювременконформационных релаксаций нескольких ворот в глобуле. Отметим,что в жесткой молекуле белка, где отсутствуют внутримолекулярныедвижения и структурные флуктуации, диффузия лиганда должна бытьсопряжена с преодолением больших активационных барьеров (до100ккал/моль). Эти барьеры настолько замедлят движение лиганда, что онопрактически станет бесконечно медленным в масштабах биологическоговремени.
Иными словами, в реальных биополимерах с плотной упаковкой именно структурные флуктуации делают возможным перенос лигандов внутри молекулы, что важно для ее функциональной активности.Более сложные виды кооперативных релаксационных процессов вбелке мы рассмотрим позднее в связи с механизмами ферментативногокатализа.Численное моделирование динамики белка-сравнительноновое направление в молекулярной биофизике. В предыдущих лекцияхбьшо показано, что конформационная энергия белка определяется атоматомными взаимодействиями и описывается специальными потенциальными функциями. В результате можно получить энергетические карты,на которыхвидныусредненныекоординатыатомов,соответствующиеусловиям минимума общего конформационного потенциала.
Подобнуюкартину дает и метод рентгеноструктурного анализа, с помощью которого определяют среднестатистические положенияатомов вбелковойструктуре. Однако такими способами невозможно проследить за движениями и флуктуациями положения отдельных атомов, которые лежат воснове конформационных флуктуаций и переходов в белках. В методечисленного моделирования динамики белка для отдельных атомов непосредственнорешаютсяклассическиеуравнениядвижения,вкоторыхдвижущие силы определены из известных потенциальных функций атоматомных взаимодействий. Исходные координаты тяжелых атомов (неводородных) задаются по рентгеноструктурным данным, причем в на-114чальный момент времени скорости их движения считаются равными вразличных, случайно выбранных направлениях.
Начальные ускорениядвижения вычисляются как первые производные сил, действующих наатомы в исходнойструктуре.Затем,на электронно-вычислительныхмашинах производится решение системы уравнений движения для всехатомов. Численное решение этих уравнений представляет собой траектории движенияотдельных атомови дает совокупнуюкартину развитияструктурных флуктуаций во времени. Пока еще удается проследить задвижениями атомов на коротком отрезке времени отt = от О до ~ 200 пс12(1 пс= 10- с), что обусловлено в основном большим объемом вычислений.В качестве примера приведем результаты моделирования внутренней динамики белка-ингибитора трипсина (ИТ) панкреатической железы, молекула которого содержит454тяжелых атома.
Оказалось, чтореальные флуктуации положений атомов в белке по отношению к усредненной во времени структуре составляют для а-углеродных 0,6 А и 0,75А для всех остальных атомов. Наблюдаются также флуктуации в значениях двугранных углов ер и \jf вращения в пептидной цепи в пределах200и для углаroв пределахзатухают в течение1-27-9°.10-Эти флуктуации положений быстропс. Однако имеются и долгоживущие, до20пс,флуктуации в положении а-углеродных атомов, которые, по-видимому,отражают конформационные переходы в белке.
Регулярность флуктуационных движений нарушается тем значительнее, чем чаще атомнаягруппа испытывает столкновения с другими атомами своего микроокружения. В пределах общего широкого конформационного минимума вбелке совершаются спонтанные переходы из одного микросостояния вдругое за счет тепловой энергии, например, вращение ароматическогокольца тирозина в молекуле ИТ. Моделирование на ЭВМ этого процессапоказало, что сам переход через потенциальный барьер происходит самопроизвольно, а не за счет сильных активационных соударений с атомами микроокружения кольца.
Кольцо тирозина пересекает барьер завремя ~ 1 пс по определенной траектории, а толчки микроокружениятолько стремятся "отвести" кольцо от барьера и "сбить" его с естественной траектории спонтанного перехода. Флуктуации положений отдельных атомов в белке коррелируют друг с другом, что может привести кбольшим по масштабу структурным сдвигам и конформационным перестройкам. Флуктуационные "дрожания" атомов создают условия и предпосьшкидляфункциональнонаправленныхконформационныхпереходов в белках. Мы еще пока далеки от построения детальной картины динамики белка. Однако уже сейчас можно сделать некоторыеобщие выводы, основанные на сопоставлении теоретических и экспериментальных результатов по внутримолекулярной подвижности и ее ролив функциональной активности белка. Атомные группы в белке испытывают на себе действие различных сил (кулоновские, ван-дер-ваальсовывзаимодействия), а также случайных "тепловых" толчков со сторонысоседних групп.
Кроме того, они могут участвовать и в нормальных ко-115лебаниях, когда небольшие смещения атомов от положения равновесияносят гармонический характер. Наблюдаются также сравнительно редкие "прыжки" с преодолением активационных барьеров. Вообще говоря,каждая атомная группа может участвовать в различных движениях, такчто общее смещение атома зависит сложным образом от вкладов различных сил. В большинстве белков амплитуда флуктуаций растет при удалении от центра к периферии молекулы. Для атомовсреднеквадратичные смещения составляютных боковых цепях - околоосновной цепи20,4 А, а для атомов на длин1,5 А.
Элементы вторичной структуры (а иимеющие водородные связи, характеризуются меньшими флуктуа13),циями, чем другие, "неорганизованные", части белка. Широкие вариацииамплитуд флуктуаций указывают на гетерогенность белка и наличиеучастков, различающихся по подвижности и гибкости. В отличие отостатков, находящихся внутри молекулы белка, движение боковых цепей,особенно содержащих заряженные группы, подвержено сильному влиянию окружающей среды (растворителя). Само по себе движение атомовносит анизотропный характер, когда случайные флуктуации сильно отличаются(иногда впо2амплитудамвзависимостиотнаправлениясмещенияраза). Направление наибольшего смещения атома определяется особенностями крупномасштабных коллективных взаимодействийэтого атома и соседних групп.
Такие кооперативные движения происходят за долгое время(~0,2::::>:10пс по сравнению с локальными смещениямипс). Корреляция мелких флуктуационных смещений приводит кпоявлению кооперативных крупномасштабных движений, включающихбольшое число атомов. Мелкие смещения играют еще роль "смазки",облегчающей крупномасштабные смещения, когда происходит вовлечение белковых доменов уже за более длинные ("физиологические") времена.Именноколлективныедвиженияважнывпроцессефункционирования белка, при котором происходят направленные смещения его структурных элементов, необходимые для достижения активной конформации.116ЛекцияМИГРАЦИЯ ЭНЕРГИИ И11.ПЕРЕНОС ЭЛЕКТРОНА В БИОСТРУКТУРАХДо сих пор мы рассматривали внутримолекулярную подвижностьбелков, которая носит самопроизвольный характер и обусловлена тепловым движением.
Однако участие белков в метаболических процессах имеханизмы функциональной активности всегда связаны с изменением ихэлектронного состояния. Именно этот фактор индуцирует в белке определенные конформациюнные переходы, направленные к достижениюнового равновесия с минимумом конформационной энергии, котороесоответствуетизмененномуэлектронномусостояниюмакромолекулы.Сами электронные переходы и изменение электронного состояния биополимеров происходят намного быстрее, чем вызванные ими конформационные перестройки. Поэтому в первом приближении электронные иконформационные переходы можно рассматривать отдельно. Электронные переходы в биополимерах имеют самостоятельное значение в целомряде важнейших биологических процессов.