Глобулярные белки можно грубо считать как разновидность НР-сополимеров. В самом деле, мономерные звенья этих белков различаются тем, что одни аминокислотные остатки являются гидрофильными или заряженными, в то время как другие гидрофобными. Мы можем очень грубо приписать первым из них индекс Р, и индекс Н остальным. Если проанализировать первичную структуру белково-подобного сополимера, полученного таким образом, и сравнить с простой первичной структурой описанной выше, то можно сделать вывод, что белковая НР последовательность более информативна и специфична. Обычно считают, что в глобулярных белках гидрофильные Р звенья покрывают поверхность глобулы, делая её устойчивой к межмолекулярной агрегации, а гидрофобные звенья Н в основном формируют ядро глобулы. Можно считать, что такое требование является достаточно ограничивающим, то есть справедливо только для малой доли всех возможных первичных структур. Следовательно, Р/Н корреляции, определённые в этом случае, зависят от конформации глобулы в целом и их следует характеризовать как дальнодействующие.

Прежде всего, определим модель и алгоритм, используемый для моделирования эволюционного процесса. Мы рассмотрим модель сополимерной цепи, которая включает только два типа звеньев, Н (гидрофобные) и Р (гидрофильные или полярные). Для простоты зафиксируют НР состав последовательности как 1:1.

Рассмотрим континуальную модель, противоположную широко используемой решёточной модели, так как динамика последней дискретна и релаксация плотной глобулы происходит достаточно медленно. Время эволюции системы определялись уравнениями Ньютона, которые решались при помощи метода молекулярной динамики (МД). Мономерные звенья, связанные упругими связями, образуют линейный НР сополимер длины N.

Гамильтониан системы включает в себя только парные взаимодействия и имеет вид:

N-1 N-2 N N

Н = åHb(rij) + å å[Hev(rij) + Ha(rij)] + 1/2åpi

2

(3.1)

i,j = i+1 i=1 j=i+2 i=1

где Hb – потенциал деформации связи, Hev принимает во внимание исключённый объём, Ha – характеризует энергию притяжения между мономерными звеньями, и i, j изменяются в диапазоне от 1 до N. Расстояние между звеньями определяется как rij = çri - rjç, где ri обозначает расположение радиус-вектора звена в трёхмерном пространстве. Последнее значение в уравнении (3.1) – классическая кинетическая энергия цепи.

Исключённый объём между всеми несвязанными звеньями моделировался при помощи отталкивающего Леннард-Джонсовского потенциала.

4e s 12 _ s 6 1 , rij £ r0

rij rij 4

Hev(rij) = (3.2)

0, rij > r0

где s = e = 1 как для Н так и для Р мономерных звеньев и r0 = 21/6 s - расстояние «обрезки».

Параметр e, входящий в уравнение контролирует шкалу энергии, тогда как s определяет расстояние на котором взаимодействуют мономерные звенья.

Квазигармонический потенциал связанных звеньев цепи имеет вид:

Cb(1) b0 12 _ 2 b 0 6 1 , rij £ b0

rij rij

Hb (rij) = (3.3)

Cb(2 ) exp rij 2 _ 1 _ rij 2 , rij > b0

b0 b0

где b0 и rij – равновесная и мгновенная длины связей соответственно. Это квазигармоническое уравнение с константами упругости Сb(1) и Cb(2) в гамильтониане (3.1) описывает взаимодействие связанных между собой звеньев. Мы использовали комбинированный потенциал вместо обычного квадратичного, так как использование простого гармонического потенциала увеличивает время достижения равновесия. Равновесная длина связи, как и другие величины длин, измерены в еденицах d,типичное значение для реальных полимеров составляет s » 5 А°. Присваиваем b0 = 1 и Cb(1) = Cb(2) = 1. Потенциал, описывающий притяжение между несвязанными мономерными звеньями имеет вид:

_ eabs 1 _ rij 2 2 , r0 < rij £ rc

rij rc

Ha (rij) = (3.4)

0, rij > rc

Параметр eab (=eНН, eРР, eНР) устанавливают глубину минимума функции притяжения невалентно связанных звеньев от расстояния и rc = 2.8 – расстояние на котором «обрезают» потенциал, описывающий притяжение. В уравнении (3.4) параметр описывающий перекрёстные взаимодействия выбирается как: eНР = (eНН´eРР)1/2 (энергия измерена в единицах kBT). В нашем эксперименте eРР – переменный энергетический параметр. При eНН = eРР = 2 фактически мы имеем гомополимерную глобулу, при этом образуется плотноупакованная конформация. При eНН = eРР = 0 не существует притяжения между звеньями. Для эффективного вычисления использовалась достаточно короткая цепь с N = 128. Предварительные результаты с N = 512 показали, что качественно большие системы сильно не отличаются.