Проблема укладки генома в вирусном капсиде стоит на стыке физики, химии и вирусологии. Соседние работы в архиве исследуют смежные аспекты: как бактериальные гены собираются в живые капли, как ионные транзисторы работают в нанопорах, как квантовые эффекты проявляются в биологических системах. Данная статья выделяется тем, что впервые систематически включает в молекулярную динамику одновременно явные ионы, явную воду, дискретные зарядовые сайты капсида и полную кулоновскую электростатику. Ранее такие расчёты либо использовали неявный растворитель, либо упрощённые модели полимера, либо фиксированные ионные условия. Здесь же ионное экранирование возникает естественно из динамики ионов, что позволяет проверять предсказания теорий среднего поля на молекулярном уровне.
Авторы построили крупнозернистую модель, в которой РНК представлена линейной цепью из 1122 бусинок трёх типов (фосфат, сахар, сахар) на основе силового поля MARTINI-2. Капсид смоделирован как двойная сферическая оболочка с внутренним радиусом около 12 нм, на внутренней поверхности которой закреплены 180 положительно заряженных сайтов, имитирующих N-концевые домены белков. Вода описана трёхбусинной поляризуемой моделью, ионы натрия и хлора — стандартными частицами MARTINI. Электростатические взаимодействия рассчитывались методом PPPM с диэлектрической проницаемостью 2.5. Моделирование проводилось в пакете LAMMPS при температуре 310 K с ланжевеновским термостатом. Для проверки модели измеряли персистентную длину свободной цепи при разных концентрациях соли (10–100 мМ) и сравнивали с экспериментальными данными для одноцепочечной РНК. Затем полимер помещали внутрь капсида и анализировали радиальные и угловые профили плотности мономеров после достижения равновесия.
Модель воспроизвела экспериментальную зависимость персистентной длины одноцепочечной РНК от концентрации соли: при росте соли от 10 до 100 мМ персистентная длина уменьшается с ~2.5 до ~1.0 нм, что согласуется с данными Chen et al. для поли(rU). Внутри капсида повышение концентрации соли от 10 до 100 мМ систематически сдвигает радиальное распределение генома от стенки к центру, как и предсказывают теории среднего поля. Однако время установления равновесия растёт в семь раз: от 64 нс при 10 мМ до 470 нс при 100 мМ. Независимое варьирование заряда капсида и плотности зарядовых сайтов показало, что суммарная электростатическая связь управляет глобальной организацией, а геометрические детали модулируют локальный контакт генома со стенкой. Угловые профили выявили двойной пик вокруг каждого зарядового сайта при 180 сайтах, который исчезает при увеличении числа сайтов до 1470, когда расстояние между ними становится меньше персистентной длины цепи. Представление N-концевых доменов одной бусинкой (выступающей) расщепляет одиночный максимум на два, что указывает на оборачивание генома вокруг выступа, а не прямое прилипание.
Работа стоит особняком среди соседних исследований: в отличие от работ по бактериальным генам или ионным транзисторам, она напрямую проверяет предсказания теорий среднего поля для укладки полиэлектролита в заряженной сферической оболочке с явным учётом молекулярных деталей. Соседние работы используют либо континуальные модели, либо неявный растворитель, либо фиксированные ионные условия. Здесь же впервые показано, что дискретные ионы и выступающие зарядовые сайты приводят к количественным отклонениям от среднего поля, особенно в угловом распределении генома. Это означает, что континуальные приближения верны для грубой радиальной картины, но не улавливают локальную организацию, важную для понимания специфичности упаковки.
Естественное развитие — включение вторичной структуры РНК (спаривания оснований) и разветвлённой топологии, что может изменить относительную роль суммарного заряда и распределения зарядов. Также важно расширить модель на многовалентные ионы (Mg2+ и др.), которые сильнее влияют на сворачивание РНК. Длинные гибкие N-концевые домены, выходящие за рамки одной-двух бусинок, могут добавить новые конформационные эффекты. Наконец, переход от равновесной организации внутри готового капсида к кинетике сборки позволит напрямую сравнивать с экспериментами по самосборке вирусов.
Результаты важны для вирусологии (понимание упаковки генома), физики полиэлектролитов (поведение заряженных цепей в ограниченном объёме) и нанотехнологий (дизайн капсидов для доставки терапевтических нуклеиновых кислот).
Ближайшие шаги — включение вторичной структуры РНК, моделирование сборки капсида вокруг генома и тестирование эффектов многовалентных ионов.
В физике полиэлектролитов остаётся открытым вопрос о точной роли дискретных ионов и корреляций в ограниченной геометрии. Данная работа показывает, что эти эффекты заметны уже при моновалентной соли и умеренных зарядах, что ставит задачу разработки более точных теорий, выходящих за рамки среднего поля. Для вирусологии ключевая нерешённая проблема — связь между последовательностью РНК, её вторичной структурой и эффективностью упаковки в капсид; данная модель создаёт платформу для её изучения.
Этот раздел адресован автору работы и его коллегам. Соседние работы нашёл смысловой поиск по нашему архиву — не по совпадению слов, а по смыслу. Мы не рецензируем: мы показываем поле вокруг.
Что сделано в работе
Авторы построили крупнозернистую модель одноцепочечной РНК внутри вирусного капсида с явным учетом ионов, воды и электростатики. Модель воспроизвела экспериментальную зависимость персистентной длины РНК от концентрации соли и показала, что повышение соли сдвигает геном внутрь капсида, при этом время установления равновесия растет в семь раз.
В чём сила
Сила работы — в систематическом включении молекулярных деталей (дискретные заряды, явные ионы, поляризуемая вода) в электростатическую модель упаковки генома, что позволяет проверять предсказания теорий среднего поля на молекулярном уровне.
Куда можно двигаться
Что можно попробовать
Почему это может быть важно
Понимание того, как ионные условия и дискретность зарядов управляют упаковкой генома, может помочь в разработке стратегий, направленных на нарушение сборки вирусов или стабилизацию капсидов для биотехнологических применений.
Работы рядом в нашем архиве
Вокруг этой работы в архиве пусто: ближайший разбор всего на 0.51, при обычных 0.64. Похоже, рядом почти не бурили.
Время установления равновесия в модели растёт в семь раз при увеличении концентрации соли от 10 до 100 мМ — это значит, что в более солёной среде геном «дольше ищет» свою окончательную укладку. Для сравнения, при 10 мМ равновесие достигается за 64 наносекунды, а при 100 мМ — почти за полмикросекунды.