В данной работе мы описываем протокол оптимизации и параметризации аминокислотных остатков, модифицированных реакционноспособными карбонильными формами, адаптируемыми к белковым системам. Этапы протокола включают в себя проектирование и оптимизацию структуры, назначение зарядов, построение параметров и подготовку белковых систем.
Карбонилирование белков реакционноспособными альдегидами, полученными в результате перекисного окисления липидов, приводит к сшивке, олигомеризации и агрегации белков, вызывая внутриклеточное повреждение, нарушение функций клеток и, в конечном итоге, их гибель. Он был описан при старении и нескольких возрастных хронических состояниях. Тем не менее, основа структурных изменений, связанных с потерей функции у белковых мишеней, до сих пор не совсем понятна. Таким образом, описан путь к построению in silico новых параметров для аминокислот, карбонизированных реакционноспособными карбонильными формами, полученными в результате окисления жирных кислот. Были построены аддукты Майкла для Cys, His и Lys с 4-гидрокси-2-ноненалом (HNE), 4-гидрокси-2-гексеналом (HHE) и формой фуранового кольца для 4-Oxo-2-ноненаля (ONE), в то время как малоновый диальдегид (MDA) был непосредственно присоединен к каждому остатку. В протоколе описываются детали построения, оптимизации геометрии, назначения зарядов, недостающих связей, углов, двугранных углов, а также их валидация для каждой модифицированной структуры остатков. В результате, структурные эффекты, индуцированные карбонилированием с этими липидными производными, были измерены с помощью моделирования молекулярной динамики различных белковых систем, таких как фермент тиоредоксин, бычий сывороточный альбумин и мембранный домен Zu-5-анкирина с использованием среднеквадратичного отклонения (RMSD), среднеквадратичной флуктуации корня (RMSF), структурного вторичного прогнозирования (DSSP) и анализа площади поверхности, доступного растворителю (SASA). в том числе.
В постоянном стремлении понять молекулярное поведение белков с окислительными модификациями, вычислительная химия стала фундаментальной опорой в широкой области научных исследований. Это опирается на использование теоретических моделей, способных интерпретировать физические явления в электронных системах, используя математические уравнения для описания атомного поведения молекул. В этом ландшафте вычислительное моделирование белков выделяется как важнейший инструмент для анализа атомного поведения молекулярных систем. Основываясь на оценке структурного поведения, энергетических расчетах и конформационныхсостояниях, эти методы становятся стратегическими союзниками в прогнозировании поведения биомолекулярных систем.
Эти симуляции специализируются на изучении структурных изменений и оценке потери или усиления биологических функций в белковых системах. Тем не менее, вычислительные подходы показали значительные ограничения при применении к белковым системам, содержащим модифицированные остатки, образованные ковалентными посттрансляционными модификациями в последовательности. Это связано с тем, что многим доступным методам не хватает ресурсов с параметрами, адаптируемыми к силовым полям, совместимыми с наиболее распространенными пакетами программ для молекулярно-динамического моделирования белков 2,3,4,5,6. Таким образом, стандартизация адаптивных параметров силового поля, совместимых с вычислительным программным обеспечением, имеет важное значение для облегчения точной связи топологий и атомных координат с уравнением, управляющим потенциальнойэнергией системы.
В ответ на эти вызовы с использованием методов ab initio был разработан протокол, адаптируемый к новым модифицированным аминокислотным остаткам с альдегидами, полученными в результате перекисного окисления липидов. В этом смысле оптимизация структурной геометрии новых остатков позволяет присвоить адаптивные заряды новым параметрам связи, угла и двугранника, которые могут быть запущены в общих силовых полях, таких как ЯНТАРЬ. Последующая валидация этих параметров позволяет определить согласованность и надежность метода, применимого к моделированию молекулярной динамики.
Одним из заметных преимуществ этого метода является его способность адаптироваться к различным посттрансляционным модификациям, от карбонилирования до фосфорилирования, ацетилирования и метилирования, среди прочего. Эта универсальность не ограничивается только белковыми системами, но распространяется и на макромолекулярные структуры, обеспечивая связь с атомными топологиями и координатами. В отличие от этого, предыдущие исследования показывают, что стандартная параметризация посттрансляционных модификаций подходит только для конкретного типа модификации и может быть получена только из опубликованных репозиториев, не имеющих возможности создаватьновые структуры.
В настоящее время проблемы прогнозирования и проектирования структуры белка становятся все более очевидными при моделировании структур с посттрансляционными модификациями. Нехватка параметров, описывающих изменения в конкретных аминокислотных сайтах, подчеркивает настоятельную необходимость разработки и применения вычислительных методов, которые могут быть адаптированы к стандартным параметризациям. Целью данного протокола является обеспечение маршрута для построения in silico новых параметров для аминокислот, ковалентно модифицированных реактивными карбонильными формами, полученными в результате окисления жирных кислот. Эти модифицированные аминокислоты распознаются общим янтарным силовым полем (GAFF) и, следовательно, могут быть использованы для оценки in silico структурных и функциональных эффектов, которые этот вид карбонилирования оказывает на их целевые белки.
Одним из важнейших шагов в разработке протокола параметризации AMBER была квантовая оптимизация новых аминокислотных остатков, модифицированных производными перекисного окисления липидов, из-за энергетической изменчивости, связанной с минимизацией и способом назначения зарядов RESP в ?…
The authors have nothing to disclose.
Эта работа была поддержана исследовательским грантом с кодом 1107-844-67943 от Ministerio de Ciencia, Tecnología e Innovación (Minciencias) и Университета Картахены (Колумбия) на грант для поддержки исследовательских групп 2021 года и Acta 017-2022.
AmberTools16 or Upper | The Amber Project | Amber is a suite of biomolecular simulation programs | |
Gaussian 09 or Upper | Gaussian Inc | Draw and optimize structures | |
Linux Ubuntu | GNU/Linux | Platform for AmberTools | |
NVIDIA GPUs GTX 1080 or Upper | Nvidia | Compatible with PMEMD |