Построение 3D-моделей циклических молекул в естественных переменныхСтатьи / Построение 3D-моделей циклических молекул в естественных переменныхСтраница 1
1. Введение
Интерес к геометрическому строению циклических молекул, интенсивно изучаемых как экспериментальными, так и расчетными методами, определяется не только их важнейшей ролью в органической химии и биохимии, но также сложностью и практически неисчерпаемым количеством соответствующих конформационных вариантов, особенно в случае гетероциклических соединений. Для построения модели (т.е. вычисления 3N декартовых координат) N-атомной молекулы в общем случае достаточно задать 3N-6 значений структурных параметров - межъядерных расстояний, валентных углов и углов внутреннего вращения, называющихся также внутренними или естественными переменными и легко оценивающихся по имеющимся эмпирическим закономерностям [1]. Оставшиеся 6 степеней свободы связаны с выбором положения и ориентацией молекулы в пространстве. Тем не менее, число структурных параметров, описывающих строение N-атомных моноциклических молекул, равно 3N (N межъядерных расстояний, N валентных углов и N углов внутреннего вращения). Из этих параметров лишь 3N-6 являются независимыми, и их значения можно выбирать произвольно (в пределах условия замыкания цикла). Оставшиеся 6 параметров называются зависимыми и определяются значениями независимых параметров.
Отметим, что пространственное строение нециклических молекул полностью описывается заданием значений N-1 межъядерных расстояний, N-2 валентных углов и N-3 углов внутреннего вращения. Замыкание цепи атомов в цикл увеличивает на единицу количество независимых межъядерных расстояний. При этом количество независимых угловых переменных уменьшается и становится недостаточным для непосредственного использования ранее рассмотренных алгоритмов построения нециклических молекул [2].
Алгоритмы построения циклических молекул по естественным переменным можно разделить на две группы.
Для итерационных методов (методы "стягивающего потенциала" и Шераги) характерна слабая чувствительность к качеству стартового приближения значений структурных параметров. Однако низкое быстродействие делает их малоэффективными при решении задач, требующих многократного построения модели молекулы (решение обратной задачи при поиске структурных параметров в дифракционных методах исследования, уточнение геометрии в методах молекулярной механики и квантовой химии, конформационный поиск и т.д.).
Алгоритмы построения геометрической модели молекулы неитерационными методами (метод Нордландера) опираются на вспомогательные геометрические построения, отличаются способом выбора 3N-6 назависимых параметров из общего их количества, работают значительно быстрее методов первой группы, однако требуют аккуратного выбора значений независимых геометрических параметров, не противоречащих условию замыкания цикла.
При сравнении алгоритмов полезно иметь в виду, что точность задания структурных параметров на основании эмпирических закономерностей существенно падает в ряду "межъядерные расстояния", "валентные углы", "торсионные углы", и включение в число зависимых параметров максимального количества торсионных углов является предпочтительным.
2. Метод "стягивающего потенциала" [3]
Очевидно, что линейная цепь атомов может быть неотличима от циклической, если подобрать соответствующие значения геометрических параметров. Подбор осуществляется итерационно, так, чтобы расстояние между концами линейной цепочки атомов постепенно приближалось к длине соответствующей химической связи. Для этого к обычному минимизируемому функционалу прибавляется так называемый "стягивающий потенциал", исчезающи по мере приближения расстояния между концами цепи к эталонному значению.
3. Метод Шераги
Авторам [4] удалось включить в набор независимых структурных параметров все N межъядерных расстояний и N валентных углов. Теперь лишь N-6 углов внутреннего вращения требуется задавать во входных данных. Оставшиеся шесть зависимых торсионных углов должны удовлетворять системе из шести уравнений, формулирующих условия замыкания цикла, сводящейся к уравнению с весьма громоздкими коэффициентами, решаемому итерационно.
4. Метод Нордландера [5]
Строится линейная цепь из N-1 атома. Если расстояние между концами этой цепи не превышают суммы длин двух оставшихся связей, то замыкание легко обеспечивается достраиванием последнего атома между этими концами. Метод формально неитерационный, но обеспечить отмеченное требование к стартовому набору структурных параметров практически невозможно без итерационного подбора.