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