МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ, ЧИСЛЕННЫЕ МЕТОДЫ И КОМПЛЕКСЫ ПРОГРАММ
УДК [004.45:519.816]:[005-027.21:330.131.7]
АЛГОРИТМ ОПРЕДЕЛЕНИЯ АКТИВНЫХ ЦЕНТРОВ МЕЖМОЛЕКУЛЯРНОГО ВЗАИМОДЕЙСТВИЯ
Статья поступила в редакцию 07.04.2018, в окончательном варианте — 26.04.2018.
Жарких Леся Ивановна, Астраханский государственный университет, 414056, Российская Федерация, г. Астрахань, ул. Татищева, 20а,
кандидат технических наук, доцент, е-mail: [email protected]
Ажмухамедов Искандар Маратович, Астраханский государственный университет, 414056, Российская Федерация, г. Астрахань, ул. Татищева, 20а,
доктор технических наук, доцент, е-mail: [email protected]
В работе предложен алгоритм моделирования процесса межмолекулярного взаимодействия двух молекул для выявления активных центров, участвующих в образовании молекулярных комплексов. Использование предложенного в работе алгоритма позволяет находить активные центры межмолекулярного взаимодействия. Полученные при этом результаты могут накапливаться в специальных базах данных, что в дальнейшем позволит облегчить решение задачи поиска веществ, с помощью которых при необходимости могут быть блокированы активные центры соответствующих молекул (например, молекул биологической клеточной мембраны при воздействии на них токсикантов). Применение алгоритма проиллюстрировано на примере взаимодействия молекул пентапептида и метионина. Произведены квантово-химические расчеты с помощью программного комплекса GAMESS с использованием полуэмпирического метода PM3. Построены адсорбционные комплексы, образующиеся в результате взаимодействия метионина и пентапептида, среди которых выбраны наиболее устойчивые. Рассчитаны основные энергетические и геометрические характеристики в адсорбционных комплексах. Выделены активные центры данного взаимодействия.
Ключевые слова: математическое моделирование, молекулярное взаимодействие, квантово-химические расчеты, схема взаимодействия, активные центры, квантово-химический программный комплекс
ALGORITHM FOR DETERMINING THE ACTIVE CENTERS OF INTERMOLECULAR INTERACTION
The article was received 07.04.2018, in the final version —.26.04.2018.
Zharkikh Lesya I., Astrakhan State University, 20a Tatishchev St., Astrakhan, 414056, Russian Federation,
Cand. Sci. (Engineering), Associate Professor, е-mail: [email protected]
Azhmukhamedov Iskandar M., Astrakhan State University, 20a Tatishchev St., Astrakhan, 414056, Russian Federation,
Doct. Sci. (Engineering), Associate Professor, е-mail: [email protected]
The article gives a detailed information on a proposed simulation algorithm of intermolecular interaction of two molecules for identifying active centers involved in the formation of molecular complexes. The use of the proposed algorithm allows toidentify the active centers of intermolecular interactions. The results can be accumulated in special databases, which will facilitate problem solution of searching substances, with the help of which, if it is necessary, the active centers of the relevant molecules can be blocked (e.g., molecules of biological cell membrane under the influence of toxicants). The algorithm application is illustrated by interaction between pentapeptide and methionine molecules. The quantum-chemical calculations by means of GAMESS software package and semi-empirical PM3 method are carried out. The adsorption complexes formed by the interaction of methionine and pentapeptide are built. The most stable ones are selected among them. The article also gibes an information on major energy and geometry characteristics in adsorption complexes. The author defines the active centers of this interaction.
Keywords: mathematical modeling, molecular interaction, quantum-chemical calculations, interaction, active centers, quantum chemical program suite
Graphical annotation (Графическая аннотация)
Введение. Прогнозирование свойств химико-биологических систем является весьма актуальной проблемой. При этом внедрение компьютерной техники и математического моделирования позволяет, во-первых, облегчить проведение исследований, во-вторых, при минимальном использовании экспериментальных данных получить достаточно надежные результаты.
К настоящему времени известны отдельные математические модели и методы, описывающие взаимодействие в основном низкомолекулярных соединений. Существующий математический аппарат позволяет: смоделировать поведение атомов в системе; теоретически исследовать геометрию, электронную структуру и энергетические характеристики систем. Однако по мере усложнения структуры молекул, за счет увеличения количества атомов в системе, указанные задачи становятся практически неразрешимыми даже при использовании компьютеров. Поэтому, при решении обширного круга задач, связанных с моделированием процессов взаимодействия в макросистемах, приходится искать новые подходы, отталкиваясь от существующих моделей с привлечением численных методов и комплексов квантово-химических программ [3, 5].
При этом одной из наиболее актуальных задач является поиск активных центров взаимодействия молекул. В отличие от понятия активного центра как химической группы молекул, определяющих специфичность их действия [4] под термином активный центр взаимодействия молекул (АЦВ) будем понимать атомы, благодаря которым реализуются межмолекулярные водородные связи.
Исходя из вышесказанного, целью данной работы стала разработка методики моделирования процесса межмолекулярного взаимодействия двух молекул для выявления активных центров, участвующих в образовании молекулярных комплексов. Применение предложенного алгоритма проиллюстрировано на примере взаимодействия молекул пентапептида и метионина.
Теоретические предпосылки и программные комплексы для моделирования межмолекулярного взаимодействия. Для достижения поставленной цели необходимо было изучить существующий арсенал программных комплексов для данной предметной области и разработать алгоритм моделирования межмолекулярного взаимодействия, при необходимости дополнив существующее программное обеспечение собственными разработками.
Одной из наиболее известных программ для моделирования и анализа структуры молекул является MOPAC [7]. Обширные справочные сведения о длинах связей, геометрическом расположении атомов, электронной плотности различных молекул и т.п., встроенные в эту программу, позволяют визуально представить необходимую молекулу для дальнейшего анализа ее структуры.
Подобный анализ позволяет выбрать наиболее устойчивую структуру молекулы исходя из минимума потенциальной энергии молекулярной системы. При этом необходимым условием наличия экстремума является малое значение среднеквадратичной нормы градиента (RMS Gradient) функции потенциальной энергии молекулярной системы E от всех её независимых геометрических координат [9]:
^ (ясЛ2 (ясЛ2 ( яеЛ2^1
RMS Gradient (E) = (3и)4
Z
SE v^ j
+
dE к5У1 j
+
se
vdzi j
где Ху, у у, Zi - декартовы координаты г-го атома; суммирование производится по всем п атомам рассматриваемым в модели.
В процессе оптимизации геометрии молекулы в программе MOPAC по умолчанию значение критерия окончания поиска (отклонение величины градиента от нулевого значения) устанавливается равным 0,1 ккал/моль/ангстрем. Несмотря на ряд преимуществ, которые программный комплекс MOPAC имеет при расчете отдельных молекул, он не позволяет проводить корректные расчеты межмолекулярного взаимодействия. Кроме того, данный комплекс не обеспечивает расчет полного перечня энергетических характеристик, которые необходимы для выявления наиболее устойчивых адсорбционных комплексов (АК). А именно атомы, участвующие в образовании водородной связи в наиболее устойчивых АК и являются искомыми активными центрами взаимодействия.
Необходимые вычисления энергетических характеристик молекул удобно проводить, используя полуэмпирический метод PM3 [10], в квантово-химическом программном комплексе GAMESS (КХПК GAMESS) [8].
При этом основной энергетической характеристикой оценки устойчивости молекулы является теплота образования АН^ (Heat of Formation), которую обычно сравнивают со справочными или экспериментальными данными. Данная величина представляет собой теплоту образования соединения из составляющих его элементов в состоянии идеального газа при температуре 298 К. Теплота образования вычисляется как разность между суммой экспериментальных значений энергий образования, составляющих молекулу изолированных атомов, и энергией атомизации Eatom :
E t = E , + E , + E ,,
atom el nucl isol >
где Eel - потенциальная энергия электронов в молекуле, вычисляемая методом Хартри-Фока [3]; Enucl -энергия электростатического взаимодействия ядер; Eisol - энергии изолированных атомов, рассчитанные
полуэмпирическим методом.
Поверхность потенциальной энергии многоатомной молекулы обычно имеет большое количество локальных минимумов. Поэтому при необходимости отыскания глобального минимума (наиболее устойчивого состояния молекулярной системы) следует проводить поиск минимума потенциальной энергии из различных начальных точек. Далее производится поиск устойчивой геометрии конфигурации молекулы, отвечающей условию достижения энергетического минимума.
Исходя из вышеизложенного, предлагается следующий алгоритм определения активных центров межмолекулярного взаимодействия.
Алгоритм поиска активных центров.
1. Выбор молекул для моделирования процесса их взаимодействия.
Выход: структурные формулы молекул.
2. Описание структуры выбранных молекул с использованием программного комплекса MOPAC.
Выход: z-матрицы для каждой из молекул.
Z-матрица (или матрица внутренних координат) не предназначена для описания химических связей в молекуле. Каждая строка z-матрицы определяет лишь тип и положение атома в молекуле. При этом первый атом задается только своим символическим кодом. Второй атом считается расположенным на оси координат, проходящей через первый и второй атомы. Поэтому для второго атома достаточно определить только расстояние до первого атома.
Положение третьего атома определяется заданием: расстояния от него до второго атома и угла, образованного третьим, вторым и первым атомом.
Для однозначного определения координат всех последующих атомов необходимо задание трех параметров:
• расстояния до ближайшего атома;
• валентного угла, под которым понимается угол, образованный направлениями химических (ковалентных) связей, исходящими из предыдущего атома;
• двугранного угла, который однозначно определяет положение рассматриваемого атома относительно плоскости, проходящей через три атома, участвующие в его описании.
3. Конвертация полученных на этапе два z-матриц с использованием дополнительно разработанного программного обеспечения [2] для дальнейшего расчета в квантово-химическом программном комплексе GAMESS наиболее устойчивой (в смысле достижения минимума потенциальной энергии) структуры молекулы.
Выход: преобразованные z-матрицы молекул для расчетов с использованием КХПК GAMESS.
4. Проведение расчетов в КХПК GAMESS наиболее устойчивых структур молекул с использованием полученных на предыдущем этапе z-матриц.
Выход: z-матрицы наиболее устойчивых структур молекул.
5. Описание молекулярных взаимодействий на основе полученной на предыдущем этапе матрицы, содержащей потенциально возможные водородные связи для адсорбционных комплексов. Данная процедура достаточно сложна, и для ее реализации целесообразно использовать какую-либо программу визуализации молекул на основе выходных файлов GAMESS. Например, можно воспользоваться программой ^етСгай [6], в которой наглядно отображаются номера и названия атомов молекул; геометрические характеристики, необходимые для построения комплексов.
Также на данном этапе задаются параметры водородных связей и проводится перенумерация атомов одной из молекул. Необходимость перенумерации обусловлена тем, что при моделировании межмолекулярного взаимодействия требуется записать общую г-матрицу адсорбционных комплексов, и рассмотреть совокупность атомов молекул, участвующих во взаимодействии, как единую систему. Для этого нумерация атомов одной из молекул осуществляется как продолжение нумерации атомов другой молекулы.
Выход: множество г-матриц адсорбционных комплексов.
6. Проведение расчетов в квантово-химическом программном комплексе GAMESS с использованием полученных в предыдущем этапе г-матриц адсорбционных комплексов.
Выход: г-матрицы адсорбционных комплексов и их геометрическое и энергетическое описание.
7. Расчет дополнительных энергетических и геометрических параметров возможных АК, с использованием специально разработанного программного продукта [2].
Выход: таблица основных параметров взаимодействий, содержащая по каждому комплексу информацию об основных величинах, характеризующих образование структур:
• R - расстояние между атомами, А;
• ЛЕадс - энергия адсорбции, кДж/моль;
• Лq - разность зарядов атомов, ё - элементарный электрический заряд (и 1,602 176 6208-10-19 Кл).
8. Исключение из множества потенциально возможных АК, таких комплексов, вероятность возникновения которых, исходя из значений энергетических и геометрических параметров, полученных на предыдущем этапе, низка.
Выход: упорядоченная по возрастанию энергии адсорбции таблица со множеством АК, вероятность возникновения которых высока. Атомы, участвующие в образовании водородной связи в наиболее устойчивых АК, являются искомыми активными центрами взаимодействия.
9. Нанесение активных центров на структурную схему молекул, образующих устойчивые АК.
Выход: структурная схема молекул с указанием активных центров.
Пример применения предложенного алгоритма.
1. Выбор молекул для моделирования. Были выбраны молекулы: пентапептида (белок) и метио-нина (витамин и). Для них была принята следующая схема нумерации атомов (составлены их структурные формулы) - рисунки 1 и 2.
Рисунок 1 - Схема нумерации атомов молекулы пентапептида
74 75 93 92 91 90
н—о н н н н
\76 771 781 791 80 81 I 89
82^С-С-С —с —Б —С-н
О 83 / 861 871 881
84^ н н н
Рисунок 2 - Схема нумерации атомов молекулы метионина
2. Описание структуры выбранных молекул с использованием программного комплекса MOPAC. В результате описания структуры выбранных молекул для каждой из них были получены г-матрицы. В качестве примера на рисунке 3 приведена г-матрица молекулы метионина, полученная в программе МОРАС.
АТОМ CHEMICAL BOND LENGTH EOND ANGLE TWIST ANGLE
NUMBER (ANGSTROMS) (DEGPEES) [DEGREES)
fD NB.J ifCjlJB:KAj I HA NB NC
1 С 0.000000 0 .000000 0. 000000
2 С * 0 . 00 00 00 0. 000000 1
3 С 111 . 2 4 3 S 9 0 0. 000000 2 1
4 с * 110 .443279 •k 175. 510823 k 3 2 £
5 О 1;2151S:1 ik 130 .014317 •k 23. 65 6866 k 4 3 2
6 О ik 114 .102370 ■k -155. 962247 k 4 3 2
7 * 114 .478067 ■k 179. 975213 * 1 2 3
8 N -ti 111 .045844 ~k -66. 002866 k 3 2 ^
9 с ■k 103 . .440518 •k -70. 279975 W 7 1 2
10 н ik 109 .62 87 81 ik 55. 00 4541 1if -i 2 3
11 к 1Л.1Ш2В..... ik 110 .556149 •k -61. 722593 * 1 2 3
12 н * 110 .309891 ~k 57. 259528 ■k 2 1 7
13 н - 110 .117129 ■k -58. 745477 2 1 7
14 н ik 108 .289899 ■k 57. 157483 ik 3 2 1
15 н k 110 .04 9216 •k -178. 136506 ■if 6 4 3
16 н k 109 .747075 •k 162. 768908 k О 3 2
17 н ■k 108 .414704 •k 42. 64 652 0 k p( 3 2
IS н IJlMiijL- * 113 .118408 ■t 62. 981976 k 9 7 1
Щ н ik 112 . .901504 -k -59. 394980 k 9 7 1
20 н 1.096904 * 107 .4452 63 ~k -178. 222998 k 9 7 1
Рисунок 3 - 2-матрина молекулы метионина
3. Конвертация г-матриц с использованием дополнительного программного обеспечения для дальнейшего расчета наиболее устойчивых структур молекул в КХПК GAMESS. В качестве примера на рисунке 4 приведен пример конвертированной г-матрицы молекулы метионина. Аналогичный вид имеет г-матрица молекулы пентапептида.
CI
C2 i l. . 5185326
C3 2 i: 1 in .2300716
C4 3 Jj253668__ 2 110 .4421651 1 175. 4546054 0
05 4 _ ..3 130 .04 40 715 2 23. 6302133 0
06 4 ^555785^ 114 .1064489 2 -135. 022 62 6 3 0
S7 i 2 114 .4686758 3 175. 5837348 0
N8 3 JZ 111 .0459 630 1 -66. 028 7456 0
С 9 1 103 . 4498 888 2 -7 0. 258203 6 0
H10 i JL1£565^ 2 105 .6183775 3 55. 0202741 0
hi: i 110 .5698287 3 -61. "055944 0
HI 2 2 110 7315 110 .3365364 " 57. 2472918 0
H13 2 110 .1056161 7 -5 8. 7806020 0
HI 4 3 JL23390^ 108 .3161211 1 57. 1386551 0
HI 5 6 J.521.M-S.,, Л 110 .0547841 3 -178. 203 858 9 0
HI 6 g 0; _3 105 ."621084 2 162. 7827468 0
H17 g ■ 0024087... ..3 108 .40 69188 2 42. 650387 7 0
HI 8 £ J 113 .1157931 1 62. 5762815 0
HIS J^IiSCL, л 112 .88 40 858 1 -55. ,388 458z 0
H2 0 9 l. .0969210 i 107 .445^089 1 -178. 2157205 0
Рисунок 4 - Конвертированная Z-матрица молекулы метионина
4. Проведение расчетов в КХПК GAMESS с использованием полученных на предыдущем этапе г-матриц. 2-матрица наиболее устойчивой (в смысле достижения минимума потенциальной энергии) структуры молекулы метионина после проведенных расчетов в квантово-химическом программном комплексе GAMESS с использованием полуэмпирического метода РМ3 приведена на рисунке 5. Аналогичный вид имеет г-матрица наиболее устойчивой молекулы пентапептида.
с:
С 2 i 1. . 5135326
СЗ 2 ш . 2300716
С4 3 _2 110 .4421£51 1 175. 4946054
05 4 3 130 .04 40 715 2 23 . 6302133
06 4 7S _3 114 .1064485 2 -156. 022 6263
S7 i 1, _2 114 .4686753 3 175. 9337343
N8 3 ¿1902131„ _2 ш .0459630 1 -66. 023 7456
С9 7 103 .44 98 888 2 -7 0. 258203 6
НЮ i JU£5651„ _2 109 . 6183775 3 55 . 0202741
НИ i шжи.. _2 110 .5698287 3 -61. "05 594 4
Н12 2 110 . 3365364 7 57 . 2472513
Н13 2 110 . 1056161 " -5 8 . 730 602 0
Hl 4 3 1, _2 108 .3151211 1 57. 1336551
н: 5 6 iL ,Л 110 . 0547841 3 -17 8 . 203 358 9
Н16 О о i 09 ."621084 2 162. 7327468
н17 Q „3 103 .40 69188 2 42. 650337 7
H1S 5 L: 1 из .1157531 1 62. 5762815
Н19 9 к-: JL 112 .8340858 1 -55. 1 о S 8 4 5 S z,
Н2 0 5 1 . .0969210 - 107 . 4457035 1 -17 в . 2157205
О О О
о о о о о о о о о о о о о о
Рисунок 5 - Z-матрица наиболее устойчивой молекулы метионина
5. Описание молекулярных взаимодействий на основе полученной на предыдущем этапе матрицы, содержащей потенциально возможные водородные связи для адсорбционных комплексов. Было составлено множество z-матриц, описывающих адсорбционные комплексы, соответствующие различным формам взаимодействия исходных молекул. Для этого к z-матрице молекулы пентапептида была добавлена z-матрица молекулы метионина с соответствующей перенумерацией атомов. Выбор такой последовательности был обусловлен тем, что в молекуле пентапептида больше атомов, поэтому легче добавлять матрицу молекулы метионина. Результаты были сохранены в файлах с расширением res. Выход: множество z-матриц адсорбционных комплексов.
6. Проведение расчетов в КХПК GAMESS с использованием полученных в предыдущем этапе z-матриц АК. В результате получаются z-матрицы наиболее устойчивых комплексов и их геометрическое и энергетическое описание. В качестве примера такого описания на рисунке 6 показан фрагмент z-матрицы одного из адсорбционных комплексов.
С2 1 1. 53"384"
03 2 „JL 124 . 1234
N4 2 „JL 116 . 0524 3 -175 . 2437
С 5 1 1, J533 7MX _2 112 . 6455 3 6. 5530
N6 1 i. JMSSM. _2 107 . 5093 3 -118 . 3660
С 7 4 _ dMiSJJL _2 120. 4821 1 -168. 8142
S8 5 1, 114. 5328 2 88. 0317
С 5 1 1, ЛШ, .....л 103. 5707 2 -135. 4811
И D . Д .
Н72 38 0. 5453175 34_ .10.8, 2022 33 -173 . 8556
Н73 40 0. 9457750 37.. 110. ,3959 30 -3 . 5866
Н74 26 2. 4055677 25 „1Ж 3361 24 145 . 805 4
07 5 74 0. 5640340 26... 167. J436 25 -164. 4806
С7 6 75 1. 3435394 2SL „1JLL .6371 26 135. 8503
С7 7 76 1. 5283063 .....1.15.-: J43S "4 -1"8. 8417
С 7 S 77 1. 5358504 2L ...1.0 В, .7377 75 76. 068 7
С7 9 7 8 1. 5175132 XL „1JUL .6362 76 -168 . 0349
SB0 79 1. 8134645 UL .1.1.3, 77 174 . 1651
С 81 80 1. 8010416 Ж 103. .2.45.6 7 8 72 . 250 0
08 2 76 1. 2203619 TL, .127, J11.5 73 -104. 8360
N83 77 1. 48164 65 78.. „1JLL .0113 79 67. 4390
Н8 4 83 0. 538-307 XL, „АЖ ■JMSl 73 88. 51" 1
HS 5 83 0. 33""262 77.. ...10.9, .6253 73 179 . 4534
Н86 7 8 1. 1057165 Ж, .10.9, .5457 80 53 . 0552
HS 7 79 1. 107 63 02 ж. 81 -166. "351
HB 8 81 1. 0953351 Ж. ...1.1.2, 310 7 75 56. 7705
нз s с -1 — 1. 03620 34 Ж 112, ^МЗЗ 79 -63. 4416
Н90 81 1. 0967173 Ж „АЖ 79 1" 5 . "333
Н51 75 1. 1063883 Ж 110.: 77 -60. 5365
Н92 73 1. 1102663 21 ...1.1.0.: JLSM 80 -62 . 770 0
Н53 77 ' 1155571 "8 105 . 2583 75 -45 . 6747
Рисунок 6 - Фрагмент z-матрицы молекулярного комплекса
7. Расчет дополнительных энергетических и геометрических параметров возможных АК. Визуализация геометрического расположения атомов в молекулах и выбор оптимальных по длине водородных связей, образующихся между молекулами комплекса, проводилась с помощью программы ^етСгай. Обработка результатов расчетов по энергетическим параметрам осуществлялась с использованием специально разработанный программного продукта [2].
8. Расчет энергий адсорбции и анализ геометрии адсорбционных комплексов. В результате исключаются из множества потенциально возможных АК комплексы, вероятность возникновения которых, исходя из значений энергетических и геометрических параметров, полученных на предыдущем этапе, низка (см. табл.).
Таблица - Свойства адсорбционных комплексов
АК По связи R, А Aq, e ДЕадс, кДж/моль АК По связи R, А Aq, e ДЕадс, кДж/моль
1 O3 - H74 1.80 -0,291 -43,14 14 H72 - O82 1.81 0,016 -27,06
2 H67 - O82 1.83 -0,002 -38,75 15 S26 - H74 2.41 -0,276 -26,03
3 O21 - H84 1.87 -0,107 -38,57 16 H61 - O82 1.85 0,0004 -22,08
4 H66 - O82 1.85 0,019 -36,34 17 O28 - H84 1.87 -0,083 -20,29
5 O39 - H74 1.79 -0,296 -34,97 18 O16 - H84 1.86 -0,1003 -18,43
6 O16 - H74 1.79 -0,293 -34,86 19 H54 - O82 1.86 0,0056 -18,18
7 H46 - O82 1.85 -0,015 -32,95 20 H62 - O82 1.86 0,0097 -17,08
8 O3 - H84 1.87 -0,100 -31,38 21 H59 - O82 1.85 0,0115 -16,96
9 O21 - H74 1.79 -0,300 -30,94 22 H65 - O82 1.87 0,007 -16,51
10 H60 - O82 1.85 0,004 -30,26 23 O39 - H84 1.86 -0,101 -16,36
11 H42 - O82 1.84 0,041 -29,08 24 H48 - O82 1.87 0,005 -15,95
12 H73 - O82 1.81 0,008 -29,06 25 H45 - O82 1.86 0,012 -14,21
13 O28 - H74 1.80 -0,291 -28,63 26 S26 - H84 2.71 -0,076 -13,56
9. Нанесение активных центров на структурную схему молекул, образующих устойчивые АК. Используя данные таблицы, наносим на схему молекулы пентапептида номера соответствующих молекулярных комплексов, приведенных в порядке возрастания значения энергии адсорбции (чем меньше номер, тем прочнее соединение). Затем указываем атом метионина, с которым потенциально возможно образование водородной связи. Построенная таким образом схема активных центров представлена на рисунке 7.
Рисунок 7 - Схема активных центров при взаимодействии метионина с пентапептидом
Выводы. Использование предложенного в работе алгоритма позволяет находить активные центры межмолекулярного взаимодействия. Полученные при этом результаты могут накапливаться в специальных базах данных [6]. Это в дальнейшем позволяет облегчить решение задачи поиска веществ, с помощью которых при необходимости могут быть блокированы активные центры соответствующих молекул (например, молекул биологической клеточной мембраны при воздействии на них токсикантов).
Список литературы
1. Жарких Л. И. Генератор z-матриц молекул. Свидетельство о регистрации программы для ЭВМ № 2015616182 от 2 июня 2015 г. / Л. И. Жарких, А. С. Дегтярев.
2. Жарких Л. И. Автоматизация расчетов основных энергетических и зарядовых характеристик при моделировании межмолекулярных взаимодействий. Свидетельство о регистрации программы для ЭВМ № 2011611798 от 28.02.2011 г. / Л. И. Жарких, Ю. А. Очередко, Н. М. Алыков, А. А. Малев.
3. Степанов Н. Ф. Квантовая механика и квантовая химия / Н. Ф. Степанов. - Москва : Мир, 2001. -519 с.
4. Фирсов Н. Н. Микробиология: словарь терминов / Н. Н. Фирсов. - Москва : Дрофа, 2006. - 256 с.
5. Элькин М. Д. Молекулярное моделирование: методические аспекты / М. Д. Элькин, Г. П. Стефанова, И. А. Крутова, В. И. Коломин // Прикаспийский журнал: управление и высокие технологии. - 2012. - № 4. - С. 103112 (http://hi-tech.asu.edu.ru/files/4(20)/103-112.pdf).
6. ChemCraft. - Режим доступа: http://www.chemcraftprog.com, свободный. - Заглавие с экрана. - Яз. рус. (дата обращения: 14.04.2018).
7. MOPAC. - Режим доступа: http://old.psu.ru/science/soft/winmopac, свободный. - Заглавие с экрана. - Яз. рус. (дата обращения: 14.04.2018).
8. Schmidt M. W. The General Atomic and Molecular Electronic Structure System / M. W. Schmidt, K. K. Baldridge, J. A. Boatz, S. T. Elbert et al. // J. Comput. Chem. - 1993. - Vol. 14. - P. 1347-1363.
9. Shewchuk Jonathan Richard. Second order gradients methods / Shewchuk Jonathan Richard // School of Computer Science. - Carnegie Mellon University Pittsburg, 1994. - Vol. 7. - P. 155-163.
10. Stewart J. J. P. Optimization of Parameters for Semiempirical Methods / J. J. P. Stewart // J. Comput. Chem. -1989. - Vol. 10, № 2. - P. 209-220.
References
1. Zharkix L. I., Degtyarev A. S. Generator z-matricz molekul [Generator of z-matrices of molecules]. Svi-detel'stvo o registracii programmy' dlya E'VM №2015616182 ot 2 iyunya 2015 g. [Certificate of registration of the computer program no. 2015616182 dated June 2, 2015].
2. Zharkix L. I., Ocheredko Yu. A., Aly'kov N. M., Malev A. A. Avtomatizaciya raschetov osnovnykh e'nergeticheskikh i zaryadovykh kharakteristikpri modelirovanii mezhmolekulyarnykh vzaimodejstvij [Automation of calculations of the main energy and charge characteristics in the modeling of intermolecular interactions] Svidetel' stvo o regis-tratsii programmy' dlya E'VM №2011611798 ot 28.02.2011g. [Certificate of registration of the computer program no. 2011611798 from 28.02.2011]
3. Stepanov N. F. Kvantovaya mekhanika i kvantovaya khimiya /Quantum mechanics and quantum chemistry]. Moscow, Mir Publ., 2001. 519 p.
4. Firsov N. N. Mikrobiologiya: slovar' terminov [Microbiology: a glossary of terms]. Moscow, Drofa Publ., 2006. 256 p.
5. E'l'kin M. D., Stefanova G. P., Krutova I. A., Kolomin V. I. Molekulyarnoe modelirovanie: metodicheskie aspekty' [Molecule modelling: methodical aspects]. Prikaspijskij zhurnal: upravlenie i vy'sokie tekhnologii [Caspian Journal: Control and High Technologies], 2012, no. 4, pp. 103-112 (http://hi-tech.asu.edu.ru/files/4(20)/103-112.pdf ).
6. ChemCraft. Available at: http://www.chemcraftprog.com (accessed: 14.04.2018).
7. MOPAC. Available at: http://old.psu.ru/science/soft/winmopac (accessed: 14.04.2018).
8. Schmidt M. W., Baldridge K. K., Boatz J. A., Elbert S. T. et al. The General Atomic and Molecular Electronic Structure System . J. Comput. Chem., 1993, vol. 14, pp. 1347-1363.
9. Shewchuk Jonathan Richard. Second order gradients methods. School of Computer Science. Carnegie Mellon University Pittsburg, 1994, vol. 7, pp. 155-163.
10. Stewart J. J. P. Optimization of Parameters for Semiempirical Methods. J. Comput. Chem., 1989, vol. 10, no. 2, pp. 209-220.
УДК 519.72
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ СИСТЕМ ЗАЩИТЫ ДАННЫХ НА ОСНОВЕ ДИОФАНТОВЫХ УРАВНЕНИЙ
Статья поступила в редакцию 27.03.2018, в окончательном варианте — 12.06.2018.
Осипян Валерий Осипович, Кубанский государственный университет, 350040, Российская Федерация, г. Краснодар, ул. Ставропольская, 149,
доктор физико-математических наук, доцент, ORCID 0000-0001-6558-7998, e-mail: v.osippyan@gmail. com
Показана объективная необходимость совершенствования систем защиты информации в условиях развития информационно-телекоммуникационных технологий. Представлены математические модели систем защиты информации, разработанные на основе линейного неоднородного и квадратного однородного диофантова уравнений. Исходным сообщением служит некоторое решение заданного диофантова уравнений, а шифртекстом - его правая часть. Изучается задача нахождения диофантова представления заданного множества с целыми числовыми компонентами. Приведены теоремы, которые позволяют описать свойства параметрических решений диофантова уравнений, необходимых для