Функции молекул в клетке, особенно пептидных и белковых, зависят от их пространственных структур. Если теоретический конформационный анализ научится надежно предсказывать такие структуры, значительный (и необходимый!) шаг на пути к направленному драг-дизайну пептидных биорегуляторов будет сделан. Членам интернационального клуба конформационных карт это было очевидно с самого начала — но их останавливала тогдашняя неполнота физических моделей молекул, начиная с неточности потенциальных функций межатомных взаимодействий и еще более приблизительного учета взаимодействий с растворителем. Отважиться на проведение конформационного расчета полной молекулы реального пептида с тем, чтобы найти все его стабильные конформации, могли лишь авантюристы. И такие появились — в нашем лице.
Объектом исследования был выбран
Рассчитанные заранее конформационные карты показали, что полный конформационный расчет брадикинина должен проанализировать как минимум пятьсот семьдесят шесть пространственных структур пептидного остова; для девятичленного пептида иной последовательности их могло бы быть порядка двадцати тысяч. Пять-шесть сотен — это тоже очень много: ведь каждую из таких конформаций следовало выбрать как начальную точку расчета межатомной энергии. Между тем при наших вычислительных ресурсах этот процесс занимал несколько часов лишь для одной конформации брадикинина. Позволить себе десятки, а то и сотни суток непрерывного счета мы не могли.
Вместо этого пришлось придумать новое упрощение — предположение, что взаимодействия аминокислот, отдаленных в пептидной последовательности друг от друга, менее влияют на конформационную стабильность молекулы в целом, чем взаимодействия отдельных фрагментов, входящих в эту последовательность. Иными словами, в первом приближении можно попытаться составить возможные конформации молекулы брадикинина, комбинируя наборы стабильных конформаций его фрагментов, — а рассчитать энергии этих наборов гораздо легче. Таким способом можно также удлинять цепочку: начать, например, с тетрапептида 1-4, затем скомбинировать его низкоэнергетические конформации со стабильными структурами тетрапептида 2-5 и провести расчет конформационных энергий пентапептида 1-5, потом добавить конформации фрагмента 3-6 и так далее. Для молекул белков, где пространственная структура зачастую определяется взаимодействиями остатков, далеко отстоящих друг от друга в аминокислотной последовательности, этот подход не годится; но для сравнительно небольших пептидов он может быть правомерен.
Полный конформационный расчет брадикинина показал возможность существования всего нескольких низкоэнергетических конформаций молекулы. Все они характеризовались резким изломом в центре пептидного остова; у некоторых к тому же боковая цепь первого в цепи аминокислотного остатка (заряженная положительно) и группа, замыкающая остов молекулы (заряженная отрицательно), близко подходили друг к другу. Оказалось, таким образом, что линейный пептид брадикинин может образовывать квазициклические конформации. Расчетная стабильность этих структур сохранялась и при нейтрализации заряженных групп, а также при введении в рассмотрение растворителя. Это была маленькая сенсация: мало того что никто в мире до тех пор не осуществил конформационный анализ пептида такой величины, было предсказано еще и новое свойство линейных пептидов — способность к самопроизвольной квазициклизации.
Главная же сенсация, однако, заключалась не в наших предсказаниях — мало ли что болтают теоретики, — а в том, что они были подтверждены независимым экспериментом. Одновременно с нашей работой московские ученые (Вадим Тихонович Иванов с сотрудниками) провели целый ряд экспериментальных измерений конформационных характеристик молекулы брадикинина. Их общий вывод был следующим: в конформационном равновесии брадикинина в растворе действительно присутствуют структуры изогнутой или псевдоциклической формы со сближенными концами пептидной цепи.