Тем не менее первые параграфы практического руководства по расчету пространственной структуры белковой глобулы в рамках любовно изложенного нами попарно-аддитивного приближения можно уже написать. Как и положено настоящему руководству, оно будет начинаться с оценки сил и средств: прикинем для начала, хотя бы примерно, объем предстоящих расчетов. Если количество атомов в белке, как говорилось, исчисляется тысячами, то счет внутримолекулярным атом-атомным взаимодействиям, как нетрудно убедиться, пойдет на миллионы (ничего не поделаешь — каждый с каждым). Следовательно, чтобы вычислить энергию внутримолекулярных взаимодействий в какой-либо конформации, необходимо рассчитать все атом-атомные расстояния (их миллионы!), найти по каждому из них энергию невалентных взаимодействий, затем, если нужно, еще и электростатических, и все просуммировать. (Вообще-то говоря, кроме этих двух типов взаимодействий, во внутримолекулярную энергию вносят вклад и другие слагаемые: например, водородные связи — особые взаимодействия, о которых, мы упоминали, говоря о структуре белковых α-спиралей, и т. д. Однако наличие этих дополнительных эффектов лишь несколько увеличивает объем вычислений, не внося ничего принципиально нового в излагаемую нами картину.)
Итак, формула, определяющая величину энергии внутримолекулярных взаимодействий в белке, представляет собой сумму нескольких миллионов слагаемых, каждое из которых зависит от расстояния между некоторой парой атомов; в свою очередь, эти расстояния зависят от углов поворота вокруг одинарных связей. (Мы уже говорили о том, что различные конформации молекулы определяются именно значениями этих (как их еще называют) углов внутреннего вращения: выражения «молекула находится в данной конформации» и «молекула обладает данным набором углов внутреннего вращения» эквивалентны.) Написать такую формулу попросту невозможно — и не только по причине существующего дефицита бумаги. Точнее, формулу-то написать можно, но она окажется такой громоздкой и запутанной, что работать с ней будет совершенно немыслимо.
Намного легче, конечно, составить программу численного расчета энергии для данного набора значений углов внутреннего вращения и обратиться за помощью к ЭВМ. Применение ЭВМ для расчета внутримолекулярной энергии — принципиальная особенность теоретического конформационного анализа белковых молекул (наш учебник вполне мог бы называться и так), и для специалиста в этой области ЭВМ такое же орудие производства, как рояль для музыканта или угольный комбайн для шахтера.
И все же на этот раз проблема состоит не в объеме вычислений, необходимых для расчета энергии, соответствующей какой-то одной конформации, — современные вычислительные машины с такими задачами справляются довольно просто. Но ведь наша конечная цель не ограничивается расчетом энергии, мы хотим найти ту трехмерную структуру, в которую белковая молекула сворачивается под влиянием различных внутренних взаимодействий. Как же это сделать, умея вычислять энергию таких взаимодействий в любой конформации молекулы?
К искреннему сожалению авторов, путь к ответу на этот вопрос вновь пролегает через дремучие дебри физики — на сей раз уже не квантовой механики, а термодинамики и статистической физики. Обе упомянутые дисциплины в трогательном согласии утверждают, что всякая система молекулярных размеров, будучи предоставлена самой себе, стремится перейти в состояние, которому соответствует наименьшая энергия. Не будем придираться к нашему пересказу формулировки, который и вправду не грешит излишней точностью: что означают, например, слова «стремится принять»?
Рассмотрим для начала случай, когда система уже находится в состоянии с наинизшей энергией. Очевидно, для того, чтобы его покинуть, ей нужен какой-то приток энергии извне. Откуда же взять эту энергию?
Как сказано, система предоставлена «самой себе», так что получить энергию как будто неоткуда, и выходит, что оставаться ей, системе, в этом самом состоянии с наинизшей энергией до конца своих дней. Однако, с другой стороны, полностью «самой себе» система никогда не бывает предоставлена: она находится в некой среде, и, если температура среды отлична от абсолютного нуля, при котором тепловое движение молекул прекращается, система постоянно испытывает толчки со стороны молекул окружения, которые, вполне возможно, смогут вывести ее из состояния с наинизшей энергией.
Здесь важно, что мы заботливо оговорили нашей системе именно молекулярные размеры. Будь ее детали существенно побольше, никакого влияния на состояние системы тепловые толчки не оказали бы. Может быть, следовало строгости ради сказать «почти не оказали бы» или «скорей всего не оказали бы» — ведь тепловое движение молекул совершенно беспорядочно, и в принципе нельзя исключить возможность того, что в некоторый момент на одну сторону какого-нибудь, скажем, поршня, вставленного в цилиндр, придется в десять раз больше толчков, чем на другую (вот, кстати, реальный путь к заветной мечте человечества — вечному двигателю).
Все тот же А. Сент-Дьёрдьи обратился однажды к своему знакомому, нобелевскому лауреату по физике, с вопросом, поверил ли бы он его рассказу очевидца о том, как письменный стол неожиданно сам по себе поднялся в воздух. Ведь молекулы стола также находятся в тепловом движении, и теоретически вполне возможно, что все они одновременно (или хотя бы большая их часть) двинулись в одном направлении — вверх. Физик ответил, что подобному рассказу он тем не менее не поверил бы: если вычислить вероятность такого (в принципе возможного) события, то она наверняка окажется намного меньше вероятности того, что А. Сент-Дьёрдьи ошибся… (Примерно ту же мысль, но более лапидарно выражает поговорка «Врет, как очевидец».)
Действительно, вероятность самопроизвольного взлета стола по упомянутой причине просто-таки невероятно мала — право же, нет никакой возможности сравнить эту величину с чем бы то ни было. Однако от нуля она все же отлична. Это обстоятельство, кстати, привело к некоторому ренессансу современной теологии, стимулом для которого послужило утверждение «Чудо — это крайне маловероятное событие». Нет оснований сомневаться в чудесах, творимых в свое время Иисусом и святыми, говорят теологи, ведь вот и наука утверждает, что, хоть малая, а вероятность таких событий есть…
Вернемся, однако, в мир, где «чудеса» возможны: к системам молекулярных размеров. Здесь энергия отдельных тепловых толчков оказывается вполне сравнимой с энергией, необходимой для перехода системы из одного состояния в другое. Если эта энергия, например, в два раза больше средней энергии теплового толчка, значит, такой переход может совершиться при одновременном получении двух толчков одинакового направления, если в три раза — трех и т. д. Ясно, что вероятность такого стечения обстоятельств будет убывать в геометрической прогрессии по мере роста количества требуемых толчков, или, что то же самое, с увеличением разности в энергиях исходного и конечного состояний.
Таким образом, для всякой молекулярной системы вероятность оказаться (под воздействием тепловых возмущений) в некотором состоянии резко падает с увеличением энергии этого состояния. Иначе это можно сформулировать следующим образом: если мы будем длительное время наблюдать за нашей системой, то окажется, что она тем реже пребывает в каждом из возможных состояний, чем большая энергия ему соответствует.
И если какое-нибудь состояние существенно выделяется среди прочих низкой энергией, система будет пребывать в нем почти все время, а если изредка и покинет его, то очень ненадолго.
Надо ли говорить, что рассказ обо всех этих подробностях бытия систем молекулярных размеров нам понадобился с вполне определенной целью. Совершенно очевидно, что под системой подразумевается белковая молекула, а различные ее состояния — это различные конформации, которым соответствуют те или иные значения внутримолекулярной энергии…