Изобретение относится к способам оптимизации условий кристаллизации белков для последующего выращивания белкового кристалла в оптимизированных условиях и расшифровки структуры белка с высоким разрешением с помощью рентгеноструктурного анализа и может быть применено в области белковой кристаллографии. Техническим результатом является удешевление и ускорение этапа оптимизационного скрининга условий кристаллизации белков. Для его достижения предложен способ оптимизации условий кристаллизации белков с применением метода молекулярной динамики, заключающийся в том, что выделяют олигомеры, содержащиеся в предкристаллизационном растворе, к полученным олигомерам применяют вычисление молекулярной динамики и изменения свободной энергии при формировании олигомера в различных условиях рН, типа осадителя, концентрации осадителя, температуры согласно следующему протоколу, заключающемуся в том, что для каждых моделируемых условий кристаллизации сначала задают рН системы путем определения состояний протонирования аминокислотных остатков в составе олигомеров с помощью программы PROPKA, после чего все расчеты и подготовку структур к ним выполняют в программном пакете GROMACS, подготовленные структуры помещают в центр ячейки моделирования, при этом минимальное расстояние между ее краем и молекулой белка составляет 1 нм, при этом остальное пространство ячейки заполняют водой, в систему добавляют компоненты осадителя в необходимой концентрации, осуществляют минимизацию энергии и NVT- и NPT-уравновешивание систем, проводят вычисление продуктивной МД в NPT-ансамбле, при этом интегрирование производят по стандартному алгоритму leap-frog с временным шагом интегрирования 2 фс и длительностью рассчитанных траекторий 100 нс, при этом перед анализом полученных траекторий устраняют артефакты, возникающие из-за квазибесконечных периодических граничных условий, по данным полученных траекторий строят графики среднеквадратичных флуктуаций атомов С? (RMSF) и вычисляют изменения энергий при образовании олигомеров белка, при этом для расчетов изменения свободной энергии при образовании олигомеров используют метод MM/GBSA, проводят вычисление свободной энергии для каждой траектории по результатам МД-моделирования, с помощью программы GROMACS, на 1000 фреймах с применением модуля gmx_MMPBSA версии 1.4.0 в сочетании со скриптом MMPBSA.py и пакетом AmberTools20, устанавливают оптимальные для кристаллизации исследуемого белка параметры, при которых наблюдаются наименьшие среднеквадратичные флуктуации атомов олигомеров и наибольшее уменьшение свободной энергии. Способ оптимизации условий кристаллизации белков с применением метода молекулярной динамики, заключающийся в том, что выделяют олигомеры, содержащиеся в предкристаллизационном растворе, к полученным олигомерам применяют вычисление молекулярной динамики и изменения свободной энергии при формировании олигомера в различных условиях рН, типа осадителя, концентрации осадителя, температуры согласно следующему протоколу, заключающемуся в том, что для каждых моделируемых условий кристаллизации сначала задают рН системы путем определения состояний протонирования аминокислотных остатков в составе олигомеров с помощью программы PROPKA, после чего все расчеты и подготовку структур к ним выполняют в программном пакете GROMACS, подготовленные структуры помещают в центр ячейки моделирования, при этом минимальное расстояние между ее краем и молекулой белка составляет 1 нм, при этом остальное пространство ячейки заполняют водой, в систему добавляют компоненты осадителя в необходимой концентрации, осуществляют минимизацию энергии и NVT- и NPT-уравновешивание систем, проводят вычисление продуктивной МД в NPT-ансамбле, при этом интегрирование производят по стандартному алгоритму leap-frog с временным шагом интегрирования 2 фс и длительностью рассчитанных траекторий 100 нс, при этом перед анализом полученных траекторий устраняют артефакты, возникающие из-за квазибесконечных периодических граничных условий, по данным полученных траекторий строят графики среднеквадратичных флуктуаций атомов С? (RMSF) и вычисляют изменения энергий при образовании олигомеров белка, при этом для расчетов изменения свободной энергии при образовании олигомеров используют метод MM/GBSA, проводят вычисление свободной энергии для каждой траектории по результатам МД-моделирования, с помощью программы GROMACS, на 1000 фреймах с применением модуля gmx_MMPBSA версии 1.4.0 в сочетании со скриптом MMPBSA.py и пакетом AmberTools20, устанавливают оптимальные для кристаллизации исследуемого белка параметры, при которых наблюдаются наименьшие среднеквадратичные флуктуации атомов олигомеров и наибольшее уменьшение свободной энергии.