Силовое поле CHARMM

Силовое поле CHARMM — это набор параметров и функциональных форм, используемых в молекулярной механике для описания потенциальной энергии молекулярных систем. CHARMM (Chemistry at HARvard Macromolecular Mechanics) является одним из наиболее широко применяемых силовых полей для моделирования биомолекул, таких как белки, нуклеиновые кислоты, липиды и углеводы, а также для изучения их взаимодействий с лигандами и растворителями. Разработка силового поля ведется с 1980-х годов в Гарвардском университете под руководством Мартина Карплуса и его коллег, и на протяжении десятилетий оно претерпело несколько поколений уточнений и расширений.

История

Разработка силового поля CHARMM началась в 1983 году, когда была опубликована первая версия программы CHARMM, включавшая базовый набор параметров для аминокислот и нуклеотидов. Изначально силовое поле было ориентировано на моделирование белков и пептидов в вакууме. В 1990-х годах, с ростом вычислительных мощностей, началось активное включение явного растворителя (воды) и липидных мембран. Ключевым этапом стало появление параметризации CHARMM22 (1998 год) для белков, которая стала стандартом де-факто для многих исследований. В 2000-х годах были выпущены версии CHARMM27 (для нуклеиновых кислот и липидов) и CHARMM36 (2012 год), которая существенно улучшила описание углеводов, липидных бислоев и ионных взаимодействий. На сегодняшний день наиболее актуальной является версия CHARMM36m (2017 год), оптимизированная для моделирования мембранных белков и гибких петель.

Функциональная форма

Потенциальная энергия в силовом поле CHARMM представляется как сумма нескольких вкладов, описывающих внутримолекулярные и межмолекулярные взаимодействия. Основное уравнение имеет вид:

\[ E_{\text{total}} = E_{\text{bond}} + E_{\text{angle}} + E_{\text{dihedral}} + E_{\text{improper}} + E_{\text{vdW}} + E_{\text{elec}} \]

Внутримолекулярные взаимодействия

• Валентные связи (\(E_{\text{bond}}\)): описываются гармоническим потенциалом \(E_{\text{bond}} = \sum_{\text{bonds}} k_b (r - r_0)^2\), где \(k_b\) — силовая константа, \(r\) — текущая длина связи, \(r_0\) — равновесная длина.
• Валентные углы (\(E_{\text{angle}}\)): также гармонический потенциал \(E_{\text{angle}} = \sum_{\text{angles}} k_\theta (\theta - \theta_0)^2\).
• Двугранные (торсионные) углы (\(E_{\text{dihedral}}\)): описываются рядом косинусов \(E_{\text{dihedral}} = \sum_{\text{dihedrals}} \sum_n k_n [1 + \cos(n\phi - \delta)]\), где \(n\) — кратность, \(\phi\) — угол, \(\delta\) — фаза.
• Неправильные торсионные углы (\(E_{\text{improper}}\)): используются для поддержания планарности (например, в ароматических кольцах) или хиральности, часто задаются гармоническим потенциалом.

Межмолекулярные взаимодействия

• Ван-дер-ваальсовы взаимодействия (\(E_{\text{vdW}}\)): описываются потенциалом Леннарда-Джонса 12-6: \(E_{\text{vdW}} = \sum_{i,j} 4\varepsilon_{ij} \left[ \left(\frac{\sigma_{ij}}{r_{ij}}\right)^{12} - \left(\frac{\sigma_{ij}}{r_{ij}}\right)^{6} \right]\), где \(\varepsilon\) — глубина потенциальной ямы, \(\sigma\) — расстояние нулевой энергии. Для смешанных пар используются правила комбинирования (обычно Лоренца-Бертло).
• Электростатические взаимодействия (\(E_{\text{elec}}\)): описываются законом Кулона \(E_{\text{elec}} = \sum_{i,j} \frac{q_i q_j}{4\pi\varepsilon_0 r_{ij}}\), где \(q_i\) — парциальные заряды атомов. Для дальнодействующих взаимодействий применяется метод суммирования Эвальда (PME) или методы усечения с переключающими функциями.

Параметризация и типы атомов

Силовое поле CHARMM использует детальную классификацию атомов по типам (atom types), которые определяют их химическое окружение, гибридизацию и функциональные группы. Например, для углерода существуют типы: CT (sp3-гибридизованный), CA (ароматический), C (карбонильный), C (альфа-углерод в пептидах) и другие. Параметры (силовые константы, равновесные длины, заряды) получают из квантово-химических расчетов высокого уровня (например, MP2/6-31G) и экспериментальных данных (рентгеноструктурный анализ, спектроскопия ЯМР, термодинамические измерения). Особенностью CHARMM является использование так называемого «химически разумного» подхода: параметры переносимы между близкими по структуре молекулами, что позволяет моделировать широкий спектр соединений без повторной параметризации.

Модели воды

Для моделирования водного окружения в CHARMM используются несколько моделей воды, совместимых с силовым полем. Наиболее распространены:

• TIP3P (Transferable Intermolecular Potential 3P) — трехточечная модель, стандартная для CHARMM. Хорошо воспроизводит плотность и диэлектрическую проницаемость, но завышает коэффициент самодиффузии.
• TIP4P-Ew — четырехточечная модель, улучшенная для использования с методом Эвальда. Дает более точные термодинамические свойства.
• CHARMM-modified TIP3P — версия с модифицированными параметрами Леннарда-Джонса для лучшей совместимости с биомолекулами.

Применение

Силовое поле CHARMM широко используется в молекулярной динамике (МД) и методах Монте-Карло для изучения:

• Структуры и динамики белков — предсказание конформационных изменений, сворачивания, стабильности.
• Взаимодействия лиганд-рецептор — докинг, свободная энергия связывания (методы FEP, TI).
• Липидных мембран — моделирование бислоев, липидных рафтов, взаимодействия с мембранными белками.
• Нуклеиновых кислот — изучение ДНК и РНК, их гибкости, взаимодействия с белками.
• Углеводов — конформационный анализ, глико-протеиновые взаимодействия.
• Материаловедения — полимеры, наночастицы, углеродные нанотрубки (с расширениями CHARMM).

Программная реализация

Исходное название CHARMM относится к программному пакету, разработанному в Гарварде. Однако силовое поле CHARMM поддерживается и другими популярными пакетами молекулярного моделирования, такими как:

• NAMD (NAnoscale Molecular Dynamics) — высокопроизводительный пакет для МД, разработанный в Университете Иллинойса.
• GROMACS — свободный пакет, в котором реализованы версии CHARMM36 и CHARMM36m.
• AMBER — поддерживает импорт параметров CHARMM через интерфейс.
• OpenMM — библиотека для GPU-ускоренного моделирования.
• Desmond (Schrödinger) — коммерческий пакет.

Сравнение с другими силовыми полями

CHARMM отличается от других популярных силовых полей (AMBER, OPLS, GROMOS) по ряду параметров:

• Парциальные заряды — в CHARMM они получаются из квантово-химических расчетов с использованием метода RESP (Restrained Electrostatic Potential) или его модификаций, в то время как в AMBER часто применяется подход HF/6-31G* с масштабированием.
• Типы атомов — CHARMM имеет более детальную классификацию (более 200 типов для биомолекул), что повышает точность, но увеличивает сложность параметризации.
• Потенциал двугранных углов — CHARMM использует ряды косинусов, что позволяет точнее описывать торсионные профили.
• Ван-дер-ваальсовы параметры — в CHARMM они оптимизированы для конденсированной фазы, в то время как AMBER изначально ориентирован на газовую фазу.

Ограничения и критика

Несмотря на широкое распространение, силовое поле CHARMM имеет ряд ограничений:

• Параметризация — требует больших вычислительных ресурсов и экспертных знаний для новых типов молекул.
• Поляризация — классическое CHARMM не учитывает поляризацию электронного облака, что может приводить к ошибкам в описании сильных электростатических взаимодействий (например, ион-ион). Существует расширение CHARMM Drude (с использованием осцилляторов Друде), но оно менее распространено.
• Точность для гибких молекул — стандартные версии могут недооценивать энтропийные эффекты для больших конформационных переходов.
• Воспроизводимость — различия в реализации в разных программных пакетах могут приводить к небольшим расхождениям в результатах.

Интересные факты

• Название CHARMM изначально было аббревиатурой, но впоследствии стало самостоятельным брендом; программа CHARMM распространяется под лицензией для академических и коммерческих организаций.
• Силовое поле CHARMM лежит в основе многих успешных предсказаний структуры белков в рамках конкурса CASP (Critical Assessment of Structure Prediction).
• В 2013 году Мартин Карплус получил Нобелевскую премию по химии за разработку многомасштабных моделей сложных химических систем, в том числе за вклад в создание CHARMM.

Источники

• Brooks, B. R., et al. (1983). "CHARMM: A program for macromolecular energy, minimization, and dynamics calculations." Journal of Computational Chemistry, 4(2), 187-217.
• MacKerell, A. D., et al. (1998). "All-atom empirical potential for molecular modeling and dynamics studies of proteins." Journal of Physical Chemistry B, 102(18), 3586-3616.
• Huang, J., & MacKerell, A. D. (2013). "CHARMM36 all-atom additive protein force field: Validation based on comparison to NMR data." Journal of Computational Chemistry, 34(25), 2135-2145.
• Best, R. B., et al. (2012). "Optimization of the additive CHARMM all-atom protein force field targeting improved sampling of the backbone φ, ψ and side-chain χ1 and χ2 dihedral angles." Journal of Chemical Theory and Computation, 8(9), 3257-3273.
• Vanommeslaeghe, K., et al. (2010). "CHARMM general force field: A force field for drug-like molecules compatible with the CHARMM all-atom additive biological force fields." Journal of Computational Chemistry, 31(4), 671-690.