Создание новых лекарственных веществ

Введение



Несмотря на достижения современной анестезии, продолжаются поиски менее опасных средств для наркоза, разработка различных вариантов многокомпонентного избирательного наркоза, позволяющего значительно уменьшить их токсичность и побочные отрицательные влияния.

Создание новых лекарственных веществ включает 6 стадий:

		Создание лекарственного вещества с помощью компьютерного моделирования.

		Лабораторный синтез.

		Биоскрининг и доклинические испытания.

		Клинические испытания.

		Промышленное производство.

		Продажа.

В последнее время компьютерное моделирование все более уверенно входит в практику технологии создания новых синтетических лекарственных веществ [1]. Предварительно проведенный компьютерный скрининг экономит время, материалы и силы при аналоговом поиске лекарственных препаратов. В качестве объекта исследования выбран местноанестезирующий препарат дикаин, который имеет более высокий уровень токсичности в ряду своих аналогов, но при этом не заменим в глазной и оториноларингологической практике. Для снижения и сохранения или усиления местноанестезирующего эффекта разрабатываются композиционные составы, дополнительно содержащие противогистаминные средства, содержащих аминоблокаторы, адреналин.

Дикаин относится к классу сложных эфиров п-аминобензойной кислоты (β-диметиламиноэтиловый эфир п-бутиламинобензойной кислоты гидрохлорид) [2]. Расстояние C-N в 2-аминоэтанольной группе определяет двухточечный контакт молекулы дикаина с рецептором через диполь-дипольное и ионное взаимодействие.

В основу модифицирования молекулы дикаина для создания новых анестетиков нами положен принцип введения химических группировок и фрагментов в существующий анестезиофор, которые усиливают взаимодействие вещества с биорецептором, снижают токсичность и дают метаболиты с положительным фармакодействием.

Исходя из этого нами предложены следующие варианты новых молекулярных структур:

			В бензельное кольцо введена “облагораживающая” карбоксильная группа, диметиламиногруппа замещена на более фармакоактивную диэтиламиногруппу.

			Алифатический н-бутильный радикал замещен на адреналиновый фрагмент.

			Ароматическая основа п-аминобензойной кислоты замещена на никотиновую кислоту.

			Бензольное кольцо замещено на пиперидиновое, характерное для эффективного анестетика промедол.

В работе выполнено компьютерное моделирование всех указанных структур с применением программы HyperChem. На последующих этапах компьютерного конструирования исследована биологическая активность новых анестетиков с применением программы PASS.

	
1.	Обзор литературы

	

	1.1	Лекарственные средства



Несмотря на огромный арсенал имеющихся лекарств, проблема изыскания новых высокоэффективных лекарственных средств остается актуальной. Это обусловлено отсутствием или недостаточной эффективностью лекарств для лечения некоторых заболеваний; наличие побочного действия некоторых лекарственных препаратов; ограничениями срока годности лекарственных препаратов; огромными сроками годности лекарственных препаратов или их лекарственных форм.

Создание каждого нового оригинального лекарственного вещества является результатом развития фундаментальных знаний и достижений медицинских, биологических, химических и других наук, проведения напряженных экспериментальных исследований, вложения крупных материальных затрат. Успехи современной фармакотерапии явились следствием глубоких теоретических исследований первичных механизмов гомеостаза, молекулярных основ патологических процессов, открытия и изучения физиологически активных соединений (гормоны, медиаторы, простагландины и др.) [3]. Получению новых химиотерапевтических средств способствовали достижения в изучении первичных механизмов инфекционных процессов и биохимии микроорганизмов.

Лекарственное средство – однокомпонентный или комплексный состав, обладающий профилактической и лечебной эффективностью. Лекарственное вещество – индивидуальное химическое соединение, используемое в качестве лекарственного средства [4].

Лекарственная форма – физическое состояние лекарственного средства, удобное для применения [4].

Лекарственный препарат – дозированное лекарственное средство в адекватной для индивидуального применения лекарственной форме и оптимальным оформлением с приложением аннотации о его свойствах и использовании [4].

В настоящее время каждое потенциальное лекарственное вещество проходит 3 стадии изучения: фармацевтическую, фармакокинетическую и фармакодинамическую.

На фармацевтической стадии устанавливают наличие полезного действия лекарственного вещества, после чего оно подвергается доклиническому изучению других показателей. Прежде всего определяется острая токсичность, т.е. смертельная доза для 50% опытных животных. Затем выясняется субхроническая токсичность в условиях длительного (несколько месяцев) введения лекарственного вещества в терапевтических дозах. При этом наблюдают возможные побочные эффекты и патологические изменения всех систем организма: тератогенность, влияние на репродуктивность и иммунную систему, эмбриотоксичность, мутагенность, канцерогенность, аллергенность и другие вредные побочные действия. После этого этапа лекарственное средство может быть допущено к клиническим испытаниям.

На второй стадии - фармакокинетической - изучают судьбу лекарственного вещества в организме: пути его введения и всасывания, распределение в биожидкостях, проникновение через защитные барьеры, доступ к органу-мишени, пути и скорость биотрансформации пути выведения из организма (с мочой, калом, потом и дыханием).

На третьей - фармакодинамической - стадии изучаются проблемы распознавания лекарственного вещества (или его метаболитов) мишенями и их последующего взаимодействия. Мишенями могут служить органы, ткани, клетки, клеточные мембраны, ферменты, нуклеиновые кислоты, регуляторные молекулы (гормоны, витамины, нейромедиаторы и т.д.), а также биорецепторы. Рассматриваются вопросы структурной и стереоспецифичной комплементарности взаимодействующих структур, функционального и химического соответствия лекарственного вещества или метаболита его рецептору. Взаимодействие между лекарственным веществом и рецептором или акцептором, приводящее к активации (стимулированию) или дезактивации (ингибированию) биомишени и сопровождающееся ответом организма в целом, в основном обеспечивается за счет слабых связей – водородных, электростатических, ван-дер-ваальсовых, гидрофобных [5].



	1.2	Создание и исследование новых лекарственных средств. Основное направление поиска



Создание новых лекарственных веществ оказалось возможным на основе достижений в области органической и фармацевтической химии, использования физико-химических методов, проведения технологических, биотехнологических и других исследований синтетических и природных соединений.

Общепринятым фундаментом создания теории целенаправленных поисков тех или иных групп лекарственных препаратов является установление связей между фармакологическим действием и физическими особенностями [3].

В настоящее время поиск новых лекарственных средств ведется по следующим основным направлениям.

	1.	Эмпирическое изучение того или иного вида фармакологической активности различных веществ, полученных химическим путем. В основе этого изучения лежит метод «проб и ошибок», при котором фармакологи берут существующие вещества и определяют с помощью набора фармакологических методик их принадлежность к той или иной фармакологической группе. Затем среди них отбирают наиболее активные вещества и устанавливают степень их фармакологической активности и токсичности по сравнению с существующими лекарственными средствами, которые используются в качестве стандарта.

	2.	Второе направление состоит в отборе соединений с одним определенным видом фармакологической активности. Это направление получило название направленного изыскания лекарственных средств.

	Преимущество этой системы состоит в более быстром отборе фармакологически активных веществ, а недостатком является отсутствие выявления других, может быть весьма ценных видов фармакологической активности.

	3.	Следующее направление поиска – модификация структур существующих лекарственных средств. Этот путь поиска новых лекарственных средств является теперь весьма распространенным. Химики-синтетики заменяют в существующем соединении один радикал другим, вводят в состав исходной молекулы другие химические элементы или производят иные модификации. Этот путь позволяет увеличить активность лекарственного препарата, сделать его действие более избирательным, а также уменьшить нежелательные стороны действия и его токсичность [4].

Целенаправленный синтез лекарственных веществ означает поиск веществ с заранее заданными фармакологическими свойствами. Синтез новых структур с предполагаемой активностью чаще всего проводится в том классе химических соединений, где уже найдены вещества, обладающие определенной направленностью действия на данный орган или ткань.

Для основного скелета искомого вещества могут быть выбраны также те классы химических соединений, к которым относятся естественные вещества, участвующие в осуществлении функций организма. Целенаправленный синтез фармакологических веществ труднее вести в новых химических классах соединений ввиду отсутствия необходимых первоначальных сведений о связи фармакологической активности со структурой вещества. В этом случае необходимы данные о пользе вещества или элемента.

Далее к избранному основному скелету вещества добавляют различные радикалы, которые будут способствовать растворению вещества в липидах и воде. Синтезируемую структуру целесообразно сделать растворимой одновременно и в воде, и в жирах с той целью, чтобы она могла всосаться в кровь, перейти из нее через гематотканевые барьеры в ткани и клетки и затем вступить в связь с клеточными мембранами или проникнуть через них внутрь клетки и соединиться с молекулами ядра и цитозоля [6].

Целенаправленный синтез лекарственных веществ становится удачным, когда удается найти такую структуру, которая по размеру, форме, пространственному положению, электронно-протонным свойствам и ряду других физико-химических показателей будет соответствовать живой структуре, подлежащей регулированию.

Целенаправленный синтез веществ преследует не только практическую цель - получение новых лекарственных веществ с нужными фармакологическими и биологическими свойствами, но и является одним из методов познания общих и частных закономерностей жизненных процессов. Для построения теоретических обобщений необходимо дальнейшее изучение всех физико-химических характеристик молекулы и выяснение решающих изменений в ее структуре, обусловливающих переход одного вида активности в другой.

Составление комбинированных препаратов является одним из наиболее эффективных путей поиска новых лекарственных средств. Принципы, на основе которых восставляются многокомпонентные лекарственные препараты могут быть различными и изменяются вместе с методологией фармакологии [7]. Разработаны основные принципы и правила составления комбинированных средств.

Чаще всего в комбинированные средства включаются лекарственные вещества, которые оказывают действие на этиологию заболевания и основные звенья патогенеза болезни. В комбинированное средство обычно включаются лекарственные вещества в малых или средних дозах, если между ними существуют явления взаимного усиления действия (потенцирование или суммирование).

Комбинированные средства, составленные с учетом указанных рациональных принципов, отличаются тем, что они вызывают значительный лечебный эффект при отсутствии или минимуме отрицательных явлений. Последнее их свойство обусловлено введением малых доз отдельных ингредиентов. Существенное преимущество малых доз состоит и в том что они не нарушают естественных защитных или компенсаторных механизмов организма.

Комбинированные препараты составляются также и по принципу включения в них таких дополнительных ингредиентов, которые устраняют отрицательное действие основного вещества.

Комбинированные препараты составляются с включением различных корригирующих средств, устраняющих нежелательные свойства основных лекарственных веществ (запах, вкус, раздражение) или регулирующих скорость освобождения лекарственного вещества из лекарственной формы или скорость всасывания его в кровь.

Рациональное составление комбинированных средств позволяет целенаправленно увеличить фармакотерапевтический эффект и устранить или уменьшить возможные отрицательные стороны действия лекарственных средств на организм.

При комбинировании лекарственных средств отдельные компоненты должны быть совместимы между собой в физико-химическом, фармакодинамическом и фармакокинетическом отношениях [8].



	1.3	Роль компьютера при создании новых лекарственных средств



Ежегодно химики синтезируют, выделяют и характеризуют от 100 до 200 тысяч новых веществ. Многие из этих веществ проходят первичные испытания на выявление той или иной биологической активности. Этот этап поиска лекарственного вещества называют скринингом. Скрининг проводят в биологических лабораториях на живых клетках, микроорганизмах или кусочках живых тканей, на здоровых или специально задержанных животных: на мышах, крысах, морских свинках, собаках.

При этом из сотен веществ отбираются несколько наиболее активных препаратов, которые затем передаются на углубление испытания. Если высокая активность вещества подтверждается, то его всесторонне изучают для определения токсичности и побочных эффектов, при отсутствии или незначительности которых проводят кинетические испытания на людях.

Считается необходимым, чтобы все новые синтезируемые вещества были подвергнуты первичным испытаниям. Очевидно, что возможность испытать все новые соединения на все новые соединения на все нужные виды активности пока остается малореальной. В настоящее время существует возможность определения потенциала их биоактивности путем компьютерного анализа [1]. Достаточно лишь ввести в компьютер сведения о строении вещества. По окончании компьютерного анализа оператор получает рекомендации о целесообразности или нецелесообразности испытаний данного вещества на тот или иной вид активности. Скрининг экономит время, материалы и силы при аналоговом поиске лекарственных веществ.

В настоящее время также пользуются методом химического модифицирования структуры известных синтетических и природных лекарственных веществ. Этот метод является интуитивным, умозрительным. С его помощью исходя из аналогии двух структур биоактивность известного вещества как бы переносят на новое соединение.

Метод молекулярного моделирования в сочетании с рентгеноструктурным анализом позволяет установить стехереохимические особенности молекулы лекарственного вещества и биорецептора, конфигурацию их хиральных центров, измерить расстояние между отдельными атомами, группами атомов или между зарядами в случае цвиттер-ионных структур лекарства и биорецепторного участка его захвата. Получаемые таким образом данные позволяют более целенаправленно проводить синтезы биоактивных молекул с заданными на молекулярном уровне параметрами. Этот метод был успешно использован в синтезе высокоэффективных анальгетиков – аналогов морфина, а также для получения ряда лекарственных веществ, действующих на центральную нервную систему подобно природному нейромедиатору γ – аминомасляной кислоты. Широкое развитее получил метод комбинаторной химии.

Метод комбинаторной химии возник и стал быстро развиваться в 1990-х годах, как часть общей стратегии открытия новых лекарственных веществ.

Стратегия комбинаторной химии основана на недавней разработке нескольких революционных химических и биологических методов параллельного синтеза и испытания большого числа соединений. Была создана техника, позволяющая синтезировать в растворе или на твердых подложках от сотен до нескольких тысяч новых соединений в день и быстро их тестировать в виде смесей или после выделения индивидуальных веществ. В совокупности с автоматизацией синтез целых семейств вещества требует значительно меньше затрат реагентов при огромном росте производительности [9].



	1.4	Молекулярное моделирование с помощью программы HyperChem



Молекулярное моделирование – сложная сеть различных наук, находящее применение в нанотехнологии, в молекулярной биологии, квантовой химии и биотехнологии.

Молекулярное моделирование молодая, востребованная и бурно развивающаяся наука.

На сегодняшний день методы квантовой химии и молекулярной динамики получили широкое распространение в численном моделировании электронной и атомной структур сложных систем молекулярных, кристаллических и переходных размеров. Это связано с технологическим развитием соответствующего математического обеспечения. Сейчас в мире функционирует достаточно много современных вычислительных комплексов, реализующих методы квантовой химии и молекулярной динамики. Использование многих из этих методов обеспечивается программой Hyper Chem для молекулярного моделирования.

HyperChem - комплексный программный продукт, предназначенный для задач молекулярного моделирования. Он включает в себя программы, реализующие методы молекулярной механики, квантовой химии и молекулярной динамики. Силовые поля, которые могут использоваться в HyperChem - это ММ+ (на базе ММ2), Amber, OPLS и BIO+ (на базе CHARMM). Реализованы полуэмпирические методы: расширенный метод Хюккеля, CNDO, INDO, MINDO/3, MNDO, AM1, PM3, ZINDO/1, ZINDO/S, а также возможности проведения неэмпирических расчетов и по теории возмущений Меллера-Плессета второго порядка. 

HyperChem обладает развитыми средствами визуализации, которые могут использоваться как при подготовке входной информации (структуры молекулы), так и при анализе результатов, например, рассчитанных характеристик ИК- и УФ- спектров.

Расчётные методы оказывают неоценимую помощь в создании лекарственных средств. Молекулярное моделирование входит во все области знаний и находит себе применение, порой играя одну из главных ролей. Некоторые области химии немыслимы без молекулярного моделирования. В развитых странах моделирование является современным методом изучения микроструктур.

В настоящее время для изучения реакционной способности молекул используются приближения CNDO/2, MNDO, AM1, PM3.

Метод CNDO основан на приближении нулевого дифференциального перекрывания и поэтому является одним из простейших полуэмпирических методов. Из этого факта следуют ограничения применимости метода, который из-за обедненной расчетной схемы недостаточно корректно воспроизводит многие эффекты. С появлением более совершенных версий полуэмпирических методов МО приближение CNDO все реже применяется на практике. Так, в версии 7 программного продукта МОРАС данный метод не представлен. Тем не менее, во многих случаях для быстрой оценки электронных параметров полезно использовать схему CNDO, так как вследствие резкого уменьшения количества рассчитываемых интегралов с помощью этого метода можно исследовать более сложные объекты. В целом CNDO/2 дает надежные результаты при расчете электронных распределений и свойств, зависящих от них.

Основным калибровочным параметром в CNDO является резонансный интеграл. Он подбирается так, чтобы относительный порядок энергетических уровней занятых МО и коэффициенты разложения МО в ЛKAO наилучшим образом совпали с расчетами ab initio соединений обучающей выборки.

Общим достоинством всех перечисленных версий является прежде всего сравнительно малое время расчетов и меньшие размеры занимаемой оперативной памяти по сравнению с более точными приближениями. Это дает возможность как для быстрой оценки исследуемых объектов, так и для изучения более сложных молекул, требующих длительного времени расчета и больших объемов оперативной памяти. В целом приближение CNDO хорошо описывает электростатические эффекты и полярность связи. CNDO/2 может применяться для расчета дипольных моментов и зарядов по схеме Малликена и оценки равновесной геометрии.

Недостатки приближения CNDO являются следствием усечения расчетной схемы, которая не учитывает взаимодействия между перекрываниями орбитальных зарядов. В результате многие эффекты не воспроизводятся.

Метод MNDO был разработан на основе более строгого и сложного приближения NDDO. Это позволило существенно улучшить результаты расчетов при решении многих задач. Длительное время метод рассматривался в качестве основного полуэмпирического метода квантовой химии. Его возможности позволили с достаточной степенью надежности рассчитывать физико-химические свойства, электронные структуры и реакционную способность множества молекулярных систем.

Преимущество заключается в быстродействии (по сравнению с неэмпирическими методами) программ, в которых реализована схема MNDO. Это позволяет применять ее для исследования все более сложных объектов. Недостатки связаны с тем, что точность метода не может превышать точность тех экспериментальных данных, по которым проводилась параметризация. 

В схеме MNDO используются 3 вида параметров. 

Во-первых, варьируемые параметры, значения которых определяют с помощью оптимизационной процедуры. 

Во-вторых, одноцентровые двухэлектронные интегралы, оценивающиеся из спектроскопических данных.

 В-третьих, ряд зависимых параметров, необходимых для расчета двухцентровых двухэлектронных интегралов, оценивающиеся с помощью эмпирических схем.

В настоящее время область применения метода MNDO достаточно изучена. Зная особенности расчетной схемы MNDO, его преимущества и недостатки можно с успехом применять метод для решения многих задач.

Хотя в целом метод MNDO имеет существенные преимущества перед СNDO, в некоторых случаях метод дает серьезные сбои. Это в первую очередь касается расчетов молекул с водородными связями, барьеров внутреннего вращения в -сопряженных системах и расчетов четырехчленных циклов. Поэтому в рамках метода MNDO были разработаны модифицированные варианты. 

Для расчета характеристик систем с водородными связями были разработаны методы MNDO/Н и MNDO/М, которые лучше воспроизводят экспериментальные значения энергии водородных связей в комплексах. 

Удовлетворительное описание водородных связей позволило широко использовать модифицированные варианты для исследования биологических объектов.

Таким образом, полуэмпирические квантовохимические методы можно использовать с большим практическим выходом для изучения реакционной способности различных химических соединений.



Достоинства метода MNDO

Недостатки метода MNDO

1.Быстродействующий метод, позволяющий изучать строение и свойства сложных молекул.

2.Учитывается ориентация р-орбиталей и правильно описывается отталкивание неподеленных электронных пар.

3.Значительное расширение круга доступных для расчета соединений.

4.Более корректное описание последовательности верхних молекулярных уровней.

5.Может использоваться для интерпретации фотоэлектронных спектров.

6.Преимущества проявляются в большей степени в расчетах более полярных молекул.

7.Удовлетворительно описывается строение радикалов, дает разумные результаты для катионов органических соединений.

8.Более точно рассчитываются валентные углы.

9.Более надежный расчет частот колебаний, протонных эффектов и электронного сродства.

10.При изучении химических реакций правильно описывается ППЭ и ПС.

1.Точность метода не может превышать точности экспериментальных данных.

2.Электронная корреляция учитывается дважды. Правильнее было бы внесение корреляционных поправок.

3.Недооцениваются эффекты взаимодействия через пространство.

4.Ограничения возможности расчета соединений третьего и более высоких периодов (в частности с SO и SO2-группами) вследствии пренебрежения d-орбиталями.

5.Некорректное описание водородных связей.

6.Неверен расчет барьеров внутреннего вращения в сопряженных молекулах.

7.Недооценивается энергия трехцентровых связей.

8.Неудовлетворительная оценка спиновых плотностей и констант сверхтонкого расщепления электронного парамагнитного резонанса.

9.Завышается стабильность радикалов.

10.Потенциалы ионизации атомов III периода с сильно выраженным S-характером завышены вследствие применения приближения замороженного остова.

11.Плохо воспроизводится спиновая плотность в радикалах.

12.Энтальпия одноатомных ионов может значительно отличаться от экспериментальных данных.

13.Ограниченная воспроизводимость свойств неорганических молекул. Ошибка расчета составляет для энтальпии - 10 ккал/моль, потенциала ионизации - 1 эВ, длин связей - 0,07

	
1.5	Определение спектра биологической активности с помощью программы PASS C&T (Prediction of Activity Spectra for Substances: Complex & Training)



Современная версия компьютерной системы предсказания спектра биологической активности PASS C&T (Prediction of Activity Spectra for Substances: Complex & Training) реализована в 1998 году. Она включает в себя обучающую выборку, содержащую более 45000 биологически активных веществ с известной биологической активностью, и охватывает более 400 фармакологических эффектов, механизмов действия, а также мутагенность, канцерогенность, тератогенность и эмбриотоксичность.

Работа PASS основана на анализе зависимостей «структура-активность» для веществ из обучающей выборки, содержащей более 45000 разнообразных биологически активных веществ (субстанции известных лекарственных препаратов и фармакологически активные соединения). Обучающая выборка постоянно пополняется новой информацией о биологически активных веществах, отбираемой как из публикаций в научно-технической литературе, так и из многочисленных баз данных. Химическая структура представлена в PASS в виде оригинальных MNA дескрипторов (Mulilevel Neighbourhoods of Atoms). MNA дескрипторы имеют универсальный характер и с достаточно хорошей точностью описывают разнообразные зависимости «структура-свойство». Используемый в PASS математический алгоритм был отобран путем целенаправленного анализа и сравнения эффективности для решения подобных задач большого числа различных методов. Показано, что данный алгоритм обеспечивает получение устойчивых в статистическом смысле зависимостей “структура-активность” и, соответственно, результатов прогноза. Это очень важно, поскольку включенные в обучающую выборку данные всегда обладают определенной неполнотой как в отношении охвата всех химических классов веществ, имеющих конкретный вид активности, так и в отношении изученности каждого отдельного вещества на все возможные виды активности. Средняя точность прогноза при скользящем контроле составляет свыше 85%. Скользящий контроль проводится следующим образом: из обучающей выборки поочередно удаляется одно вещество и для него делается прогноз на основе анализа оставшейся части обучающей выборки, результат сравнивается с известными экспериментальными данными. Процедура повторяется итеракивно для каждого из веществ и рассчитывается средняя точность прогноза. Точность прогноза в 85% достаточна для практического применения системы PASS с целью прогноза спектра биологической активности новых веществ, поскольку ожидаемая вероятность случайного угадывания одного из 780 видов активности составляет около 0.1% [10].

Результаты прогноза выдаются пользователю в виде списка названий вероятных видов активности с расчетными оценками вероятностей наличия (Pa) и отсутствия каждого вида активности (Pi), которые имеют значения от 0 до 1. Эти вероятности рассчитываются независимо по подвыборкам активных и неактивных соединений, и поэтому их сумма не равна единице. Pa и Pi интерпретируются как оценки меры принадлежности вещества к классам активных и неактивных соединений соответственно, либо как оценки ошибок первого и второго рода. Чем больше для конкретной активности величина Pa и чем меньше величина Pi, тем больше шанс обнаружить данную активность в эксперименте. Если при анализе прогнозируемого списка активностей для исследования выбираются те виды активности, для которых Pa>90%, то мы рискуем пропустить около 90% действительно активных соединений, но вероятность ложноположительных прогнозов при этом ничтожно мала; для Pa>80% - пропустим уже только 80% активных соединений, но и вероятность ложноположительных прогнозов будет выше, наконец, для Pa>Pi вероятности ошибок первого и второго рода равны [1, 9, 10].

На практике, однако, при отборе для исследования наиболее перспективных веществ руководствуются и другими критериями, например, критерием новизны. При этом исходят из того, что чем ближе значение Pa к единице, тем более вероятно, что вещество является близким аналогом известного препарата. Поэтому, если целью исследователя является выявление соединений с достаточно высоким уровнем новизны (New Chemical Entity, NCE), то надо выбирать вещества, для которых величина прогнозируемой вероятности Pa для требуемого вида активности несколько ниже, например, 0,5<PA< NCE.

Базируясь на данных компьютерного прогноза, исследователь может: 

		определить, какие тесты наиболее адекватны для изучения биологической активности конкретного химического соединения. 

			обнаружить новые эффекты и механизмы действия для ранее изученных веществ; 

			отобрать наиболее вероятные базовые структуры новых лекарств с требуемым биологическим действием среди доступных для скрининга химических соединений.

Система PASS позволяет получить прогноз спектра биологической активности 1000 веществ на обычном персональном компьютере менее чем за одну минуту. Поскольку прогноз выполняется по структурной формуле вещества, он может быть выполнен уже на стадии планирования синтеза [11].

	Применимость PASS для решения практических задач продемонстрирована в многочисленных экспериментах. Прогнозируемые виды активности подтверждены для веществ различных химических классов, проявляющих разнообразные эффекты: противобактериальный, антиаритмический, противоопухолевый, гепатопротекторный, антиамнестический, противоспалительный, антиоксидантный и др. 

	С применением PASS при поддержке гранта CRDF (RC1-2064) был выполнен прогноз спектра биологической активности для 250000 химических соединений.

Кроме того, если, наряду с основным действием, известен перечень нежелательных побочных эффектов, то при отборе перспективных для исследований соединений можно руководствоваться комбинированным критерием: 

		наличие в прогнозируемом спектре требуемых эффектов/механизмов; 

		отсутствие нежелательных эффектов/механизмов.

Естественно, что при рассмотрении всего списка, включающего свыше 400 прогнозируемых PASS C&T видов активности, можно составить большое количество комбинаций из требуемых и нежелательных эффектов. 

Для их анализа сотрудник Лаборатории структурно-функционального конструирования лекарств НИИ Биомедхимии РАМН А. А. Лагунин разработал специальную компьютерную систему интерпретации спектров биологической активности веществ IBIAC, основанную на знаниях об известных взаимосвязях между фармакологическими эффектами и механизмом действия биологически активных веществ (более 2000 терминов, описывающих биологическую активность). 

С использованием системы IBIAC генерация перечня эффектов, соответствующих определенному механизму действия и, наоборот, списка вероятных механизмов, ответственных за проявление определенного эффекта, осуществляется автоматически [12].

Поскольку прогноз спектра биологической активности осуществляется на основе структурной формулы химического соединения, он может быть выполнен уже на этапе планирования синтеза. 

В итоге будут синтезированы лишь некоторые из теоретически возможных производных, в наибольшей степени удовлетворяющие критериям задачи.

Необходимо отметить, что прогноз спектра биологической активности возможен для низкомолекулярных органических (drug-like) соединений, структура которых не отличается принципиально от веществ обучающей выборки. Не имеет смысла прогноз для синтетических и биополимеров, для неорганических веществ и т.п [13, 14, 15, 16].

	
1.6	Местноанестезирующие средства



Несмотря на достижения современной анестезии, продолжаются поиски менее опасных средств для наркоза, разработка различных вариантов многокомпонентного избирательного наркоза, позволяющего значительно снизить дозы используемых средств, уменьшить их токсичность и побочные отрицательные влияния. Известно, что местноанестезирующие средства понижают чувствительность окончаний афферентных нервных волокон, и угнетают проведения возбуждения по нервным волокнам. Вызывают местную потерю чувствительности. И в первую очередь они устраняют чувство боли, в связи с чем их используют главным образом для местного обезболивания (анестезии) [17].

К анестезирующим средствам предъявляют определенные требования. Прежде всего они должны иметь высокую избирательность действия, не оказывая отрицательного влияния ни на нервные элементы, ни на окружающие ткани. Короткий латентный период, высокая эффективность при разных видах местной анестезии, определенная продолжительность действия. Желательно, чтобы они сужали кровеносные сосуды. Это существенный момент, так как сужение сосудов усиливает анестезию, понижает кровотечение из тканей, а так же уменьшает возможность токсических эффектов задерживая всасывание анестетика. Так же к числу важных характеристик относятся низкая токсичность и минимальные побочные эффекты.

В настоящее время в медицинской практике используется множество местноанестезирующих веществ с различной степенью активности и разной продолжительностью действия.

По применения в клинической практике местные анестетики подразделяются на:

	1)	средства, применяемые только для поверхностной анестезии: кокаин, дикаин (тетракаин), бензокаин (анестезин), бупивакаин (пиромекаин);

	2)	средства, применяемые для инфильтрационной и проводниковой анестезии: новокаин (прокаин), тримекаин, бутивокаин, изокаин, ультракаин.

	3)	средства, применяемые для всех видов анестезии: лидокаин (кскаин).

	По химическому строению местноанестезирующие вещества можно разделить на две группы:

	1)	сложные эфиры: кокаин, дикаин, бензокаин, прокаин;

	2)	замещенные амиды кислот: лидокаин, тримекаин, бупивакаин, булекаин.

Дикаин применяется для поверхностной анестезии. Под влиянием местных анестетиков в окончании нерва и в самом нерве прекращается электрохимический процесс, осуществляющий передвижение ионов через мембрану и распространение нервных импульсов. Область окончания чувствительных нервов находится под регулирующим влиянием системы медиаторных рецепторов, синергентное взаимодействие, которых обеспечивает более эффективное развитее торможения.

Так же известно, что в присутствии противогистаминных средств (димедрол), м – холиноблокаторов (атропин) и адренанила усиливается местноанестезирующий эффект [18].

Дикаин представляет собой  - диметиламиноэтиловй эфир п – бутиламинобензойной кислоты гидрохлорид).







Работая над компьютерным дизайном местноанестезирующего препарата дикаин, был проведен информационный анализ данного препарата по сравнению с другими средствами, такими как бензокаин (анестезин), бумекаин. В ходе анализа было установлено, что дикаин является эффективным местноанестезирующим средством, значительно превосходящим по активности кокаин (примерно в 10 раз), он также превосходит кокаин по токсичности (в 2-5 раз).

Применяют дикаин в глазной практике в виде 0,1 % раствора при измерении внутриглазного давления, в виде 0,25-1% или 2% раствора при удалении инородных тел. Дикаин в отличие от кокаина не влияет на внутриглазное давление, не расширяет зрачки. При необходимости длительной анестезии используют глазные пленки с дикаином (0,75 см) изготовленные на основе биорастворимого полимера.

Анестезиофорной группой является диалкиламиноацетанилид. Расстояние между  и атомами определяет двухточечный контакт молекулы дикаина с рецептором через диполь-дипольное и ионное взаимодействие. В структуру молекулы дикаина входит фрагмент биогенного вещества коламина, производные которого оказывают противогистаминный эффект. Дикаин активнее анестезина и новокаина, но токсичнее их в несколько раз. В проведении дальнейшего анализа были предложены варианты новых структур для компьютерного анализа.


2. Методическая часть



2.1 Характеристика объекта исследования



Температура плавления 147-150. Соединение представляет белый кристаллический порошок без запаха. Легко растворим в воде и спирте. Дикаин сильное местноанестезирующее средство, обладающее высокой токсичностью. Применяют в глазной и оториноларингологической практике при некоторых оперативных вмешательствах, а также для анестезии.



2.1.1 Характеристика приборов

Методы исследования проводили с помощью программы Hyper Chem.



2.1.2 Методика выполнения

Построить молекулы при помощи программы HyperChem.

Последовательность действия:

		Запуск программы Hyper Chem.

		Создание молекулы:

- мышь устанавливают на пункт меню “Build”, щелчком по левой кнопке разворачивают меню и выбирают “Default atoms”;

- устанавливают указатель курсора в режиме построения молекулярных моделей; для этого в верхней строчке выбирают необходимый вид курсора и щелкают левой кнопкой мышки; 

- в развёрнутой на экране периодической таблице выбирают интересующий атом;

- связь между атомами обозначают, нажимая левую кнопку мыши в положение одного из атомов, и, удерживая ее, передвигают курсор к другому атому; затем кнопку отпускают, типы связей задают. Полуторные связи ароматического кольца обозначаются пунктирной линией. Их можно нарисовать при помощи двойного L-щелчка вблизи одной из внутренних сторон кольца;

- построенную модель автоматически дополняют атомами водорода; для этого мышь устанавливают на пункт меню “Build”, щелчком по левой кнопке разворачивают меню и выбирают “Add Hydrogens”.

- для устранения неточностей выполненного рисунка мышь устанавливают на пункт меню “Build”, щелчком по левой кнопке разворачивают меню и выбирают “Model build”; данная команда корректирует межатомные расстояния и углы.

С помощью программы Нурег Chem:

- курсор мыши устанавливают на пункт меню "Setup", выбирают метод молекулярной механики), устанавливают "ММ";

- запускают процесс оптимизации геометрии путем выбора пункта меню "Compute", далее выбирают "Geometry Optimize";

Определение геометрических характеристик (длин связей и валентных углов):

- выбирают курсор в виде двух концентрических окружностей, ставят этот курсор на один из интересующих атомов, нажимают левую кнопку мыши и, не отпуская ее, подводят курсор к следующему атому (для измерения длины связи) или к атому, находящемуся через один от исходного (для измерения величины валентного угла); затем кнопку отпускают. 

В нижней строке экрана появится значение длины связи (Е) или валентного угла (град.).

- в пункте меню "File, выбирают "Start Log" (создание файла отчета); файлу дают название и устанавливают "Quantum print level" = 9;

- в пункте меню "Setup", выбирают "Semiempirical methods"; в раскрывшемся окошечке устанавливают "MNDO";

- устанавливают соответствующий заряд и мультиплетность в соответствующих полях;

- запускают процесс расчета с оптимизацией геометрии путем выбора пункта меню "Compute", выбирают "Geometry Optimize";

- расчет заканчивается, когда в нижней строке окна появляется надпись Conv=YES"- закрывают файл отчета (.log file) путем выбора пункта меню "File", выбирают "Stop Log".

Графическое изображение ВЗМО и НВМО.

- для получения графического изображения молекулярной орбитали выбирают пункт меню "Compute", выбирают "Orbitals";

- выбирают номер нужной молекулярной орбитали, и устанавливают "3D";

- полученную картинку можно скопировать, используя пункт меню "Edit" и далее "Copy image".

- выбирают пункт меню "Compute", затем "Plot molecular properties";

- выбирают "electrostatic potencial" и устанавливают "3D". Положительный знак электростатического потенциала отображается зеленым цветом. В области неподеленных пар на атомах азота, кислорода и др. электростатический потенциал отрицательный, что отображается красным цветом.

- выделить в нарисованной молекуле связь;

- в меню Comput и выбрать пункт Potencial;

- задать начальную(Initial Bond Angle), конечную(Final Bond Angl) длины связи и шаг (Step) с которым будут автоматически производиться расчеты.

-после расчетов машина построит на экране график изменения потенциальной энергии выбранной связи от величины растяжения.


2.2 Определение биологической активности молекул с помощью программы PASS



2.2.1 Методика выполнения

Построить молекулы при помощи программы Isis Draw.

Последовательность действия:

Запуск программы Isis Draw.

Создание молекулы:

- мышь устанавливают на клавишу “Single Bound” для того, что бы изобразить фрагмент связи;

- устанавливают указатель курсора на клавишу “Atom” для того, что бы выбирать интересующий атом;

- построенную модель сохраняют, устанавливая курсор на пункт меню “File”, затем “Export”, устанавливают на “Molfile” и сохраняют в директорию Isis Draw.

Последовательность действия:

Запуск программы PASS.

	2.	Исследование биологической активности:

- мышь устанавливают на пункт меню “File”, затем “Open Base” выбираем “passdemo.SAR” и нажимаем “Открыть”;

- устанавливают указатель курсора на клавишу “Predict and Save Prediction as Text files”, тип файлов выбираем “mol”, нажимаем “Открыть”, затем сохраняем как “txt”, открываем сохраненный файл и изучаем биологическую активность.


3. Результаты эксперимента и их обсуждение



3.1 Разработка новых молекулярных структур на основе дикаина



Несмотря на достижения современной анестезии, продолжаются поиски менее опасных средств для наркоза, разработка различных вариантов многокомпонентного избирательного наркоза, позволяющего значительно снизить дозы используемых средств, уменьшить их токсичность и побочные отрицательные влияния [19]. В последнее время методы компьютерного моделирования все более входят в практику технологии создания новых синтетических лекарственных веществ [9]. Этот подход позволяет установить стехиометрические особенности молекулы лекарственного соединения, измерить расстояние между отдельными атомами, определить потенциал биоактивности, комплементарность взаимодействия с биорецепторным участком. Получаемые таким образом данные позволяют более целенаправленно проводить синтезы биоактивных молекул с заданными на молекулярном уровне параметрами, что значительно экономит время, материалы и силы при аналоговом поиске лекарственных веществ.

Дикаин применяется для поверхностной анестезии. Под влиянием местных анестетиков в окончании нерва и в самом нерве прекращается электрохимический процесс, осуществляющий передвижение ионов через мембрану и распространение нервных импульсов. Область окончания чувствительных нервов находится под регулирующим влиянием системы медиаторных рецепторов, синергентное взаимодействие, которых обеспечивает более эффективное развитие торможения. Известно, что в присутствии противогистаминных средств (димедрол), м-холиноблокаторов (атропин) и адреналина усиливается местноанестезирующий эффект [20].

Дикаин относится к классу сложных эфиров п-аминобензойной кислоты (β-диметиламиноэтиловый эфир п-бутиламинобеизойной кислоты гидрохлорид) [8].






Анестезиофорной группой является диалкиламиноацетанилид. Расстояние между  и  атомами определяет двухточечный контакт молекулы дикаина с рецептором через диполь-дипольное и ионное взаимодействие. В структуру молекулы дикаина входит фрагмент биогенного вещества коламина (2-аминоэтанола), производные которого оказывают противогистаминный эффект. Дикаин активнее своих аналогов (анестезина и новокаина), но и токсичнее их в несколько раз. Его используют главным образом в глазной и оториноларингологической практике.

Нами предложены варианты новых структур для компьютерного дизайна молекулы дикаина с целью снижения его токсичности с сохранением или даже усилением анестезирующих свойств.

Введение в бензольное кольцо «облагораживающей» карбоксильной группы и замена диметиламиногруппы на более фармакоактивную диэтиланиногруппу позволит снизить токсичность соединения, облегчить гидролиз сложноэфирной связи с высвобождением антигистаминного фрагмента - диэтиламиноэтанола.







Алифатический радикал н-бутил в структуре дикаина усиливает фармакологический эффект. При замене его на адреналиновый фрагмент ожидается получить более яркое анестезирующее действие.






К настоящему времени известно, что биологические системы не делают различия между плоскими кольцами, поэтому при замене ароматической основы н-аминобензойной кислоты на никотиновую (или изоникотиновую) кислоту изменяется полярность молекулы, облегчается задача введения различных заместителей в ароматическое кольцо. К тому же, аминопроизводные никотиновой кислоты (кордиамин) являются стимуляторами центральной нервной системы.







Один из наиболее эффективных анестетиков, промедол, содержит в структуре вместо ароматического пиридинового кольца пиперидиновое, что является предпосылкой для модификации молекулы дикаина.






3.2 Результаты исследований с помощью программы HyperChem



Таблица 1

Длина связи или валентный угол

Данные ММ расчёта

Данные MNDO расчёта

Справочные величины

С1-С2

С2-С3

С3-С4

С4-N5

N5-C6

C6-C7

C7-C9

C9-C11

C10-C11

C8-C10

C6-C8

C11-C12

C12-O13

C12-O14

O14-C15

C15-C16

C16-N17

N17-C19

N17-C18

C1-C2-C3

С2-С3-С4

C3-C4-N5

N5-C6-C8

N5-C6-C7

C6-C7-C9

C7-C9-C11

C9-C11-C10

C11-C10-C8

C10-C8-C6

C8-C6-C7

C9-C11-C12

C10-C11-C12

C12-O14-C15

O14-C15-C16

C15-C16-N17

C16-N17-C18

C16-N17-C19

C18-N17-C19

1,53442

1,53705

1,53626

1,44712

1,35045

1,34373

1,34335

1,34525

1,3446

1,34261

1,34354

1,36294

1,21202

1,35065

1,40995

1.53635

1,45672

1,45174

1,45291

111,792

113,518

109,67

118,219

123,605

120,397

120,124

119,658

119,913

120,666

119,242

119,948

122,115

118,556

107,759

113,916

111,47

113,282

111,148

1,53158

1,54128

1,54862

1,46777

1,41076

1,41935

1,40351

1,41482

1,41571

1,40009

1,42528

1,49637

1,22988

1,36302

1,41167

1,56222

1,47112

1,46359

1,46386

114,547

113,518

111,073

118,219

123,605

120,589

121,461

117,898

121,164

120,879

117,974

119,948

122,115

124,819

107,759

109,436

116,806

117,015

116,056

1,534

1,536

1,539

1,456

1,411

1,393

1,355

1,381

1,392

1,356

1,352

1,383

1,218

1,386

1,382

1,559

1,493

1,432

1,451

113,125

113,529

110,053

118,953

123,983

120,159

120,956

118,241

120,971

120,264

118,325

119,948

122,233

121,179

107,563

110,563

115,963

115,912

116,023



Сравнивая полученные результаты, оба метода имеют небольшие отклонения по сравнению со стандартными данными.

Расположение молекулы немного изменяется в зависимости от применяемого метода.

Положительный знак электростатического потенциала отображается зелёным цветом. В области неподелённых пар на атомах азота, кислорода электрический потенциал отрицательный, что отображается красным цветом.

Изменение потенциальной энергии связи C16-N17.

Исследуемая связь между атомами C16-N17. Задаём начальные величины начальная длина связи 0,971; конечная длина связи 2,971; шаг 0,1. Проведя расчёты, изобразим график изменения потенциальной энергии связи C16-N17 от величины растяжения.

Изменение потенциальной энергии углового напряжения C3-C4-C5.

Исследуем угол между атомами C3-C4-C5. Задаём начальные величины начальный угол 80; конечный 140; шаг 10. Проведя расчёты, изобразим график изменения энергии углового напряжения при изменении величины угла.



Таблица 2

Длина связи или валентный угол

Данные ММ расчёта

Данные MNDO расчёта

Справочные величины

C1-C2

C2-C3

C3-C4

C4-N5

N5-C6

C6-C7

C7-C9

C9-C10

C10-C11

C11-C8

C6-C8

C7-C12

C12-O13

C12-O14

C10-C15

C15-O16

C15-O17

O17-C18

C18-C19

C19-N20

N20-C23

C23-C24

N20-C21

C21-C22

C1-C2-C3

C2-C3-C4

C3-C4-N5

C4-N5-C6

N5-C6-C7

N5-C6-C8

C6-C7-C12

C7-C12-O13

C7-C12-O14

C6-C7-C9

C7-C9-C10

C9-C10-C11

C10-C11-C8

C6-C8-C11

C8-C6-C7

C11-C10-C15

C9-C10-C15

C10-C15-O16

C10-C15-O17

O16-C15-O17

C15-O17-C18

O17-C18-C19

C18-C19-N20

C19-N20-C21

C19-N20-C23

N20-C23-C24

N20-C21-C22

1,53446

1,53687

1,53769

1,46773

1,35455

1,34989

1,34842

1,3461

1,34288

1,34136

1,34515

1,36573

1,20944

1,33807

1,36377

1,21204

1,35143

1,41104

1,53791

1,46268

1,45983

1,53447

1,46076

1,53704

111,785

112,138

110,841

116,92

124,692

115,847

123,329

123,901

121,905

119,209

120,947

119,665

119,499

121,203

119,456

119,282

121,051

119,494

121,218

119,22

118,218

108,252

114,762

111,29

111,699

116,769

113,058

1,53154

1,54091

1,5495

1,47515

1,43364

1,423

1,41348

1,41233

1,41258

1,40248

1,41671

1,50204

1,22891

1,35818

1,50134

1,22774

1,35955

1,41312

1,5631

1,47233

1,47481

1,53447

1,47478

1,53801

111,785

113,622

110,841

119,144

118,391

122,522

121,784

126,174

144,416

119,516

121,22

118,829

120,542

120,848

118,99

121,044

120,117

125,793

113,418

120,418

125,334

107,378

110,057

116,349

116,243

112,593

112,593

1,532

1,539

1,544

1,483

1,423

1,312

1,342

1,401

1,362

1,341

1,395

1,502

1,211

1,348

1,551

1,221

1,355

1,413

1,531

1,452

1,453

1,534

1,460

1,538

111,785

112,953

110,841

117,526

117,318

117,286

122,463

123,901

144,411

119,291

121,231

118,973

119,513

121,203

119,361

120,964

121,051

120,169

119,432

120,111

119,213

107,314

112,953

113,561

114,246

111,548

112,987


Сравнивая полученные результаты, оба метода имеют небольшие отклонения.

Расположение молекулы немного изменяется в зависимости от применяемого метода.

Положительный знак электростатического потенциала отображается зелёным цветом. В области неподелённых пар на атомах азота, кислорода электрический потенциал отрицательный, что отображается красным цветом.

Изменение потенциальной энергии связи C6-С11.

Исследуемая связь между атомами C6-С11. Задаём начальные величины начальная длина связи 1,037; конечная длина связи 3,037; шаг 0,1. Проведя расчёты, изобразим график изменения потенциальной энергии связи C6-C11 от величины растяжения.

Изменение потенциальной энергии углового напряжения C6-C11-O14.

Исследуем угол между атомами C6-C11-O14. Задаём начальные величины начальный угол 90; конечный 150; шаг 10. Проведя расчёты, изобразим график изменения энергии углового напряжения при изменении величины угла.

Объектом является химическое соединение.



Таблица 3

Длина связи или валентный угол

Данные ММ расчёта

Данные MNDO расчёта

Справочные величины

O1-C2

C2-C3

C3-C4

C4-C5

C5-C7

C6-C7

C2-C6

C3-O8

C5-C9

C9-C10

C10-N11

N11-C12

N11-C13

C13-C14

C14-C15

C15-C16

C16-C18

C17-C18

C13-C17

C16-C19

C19-C20

C19-O21

O21-C22

C22-C23

C23-N24

N24-C25

N24-C26

O1-C2-C6

O2-C2-C3

C2-C3-C4

C3-C4-C5

C4-C5-C7

C5-C7-C6

C2-C6-C7

C3-C2-C6

O8-C3-C2

O8-C3-C4

C4-C5-C9

C7-C5-C9

C5-C9-C10

C9-C10-N11

C10-N11-C12

C13-N11-C12

C14-C13-N11

C17-C13-N11

C13-C14-C15

C14-C15-C19

C15-C16-C18

C16-C18-C17

C13-C17-C18

C17-C13-C14

C19-C16-C18

C19-C16-C15

O20-C19-C16

O21-C19-C16

O20-C19-O21

C19-O21-C22

O21-C22-C23

C22-C23-N24

C23-N24-C25

C23-N24-C26

C25-N24-C26

1,36051

1,34612

1,34441

1,3441

1,34399

1,34099

1,34297

1,36079

1,51602

1,54091

1,46417

1,45765

1,35257

1,3425

1,34235

1,34536

1,34524

1,34278

1,34206

1,36228

1,21202

1,35004

1,4099

1,53786

1,45675

1,45323

1,45163

120,802

120,483

119,222

122,579

117,432

120,669

121,368

118,714

120,184

120,594

121,181

121,385

115,402

114,762

123,704

116,621

120,198

119,72

120,071

119,877

120,004

119,927

120,013

120,082

119,383

120,612

119,581

120,826

119,469

119,174

110,679

113,688

111,433

113,311

111,143

1,35748

1,43712

1,4148

1,41582

1,41064

1,4048

1,41425

1,35924

1,51423

1,55179

1,47428

1,47004

1,43737

1,41646

1,40524

1,41186

1,41199

1,4048

1,41809

1,501

1,22707

1,36061

1,41166

1,5624

1,47048

1,46449

1,46417

123,264

118,062

119,488

121,583

118,032

121,572

120,637

118,674

117,63

122,882

121

120,916

117,262

113,931

116,04

114,642

123,731

117,39

120,376

120,696

118,949

120,676

120,373

118,878

120,928

120,123

126,438

112,643

120,919

125,616

107,249

109,754

116,591

116,835

115,996

1,359

1,381

1,368

1,388

1,347

1,358

1,413

1,362

1,514

1,551

1,469

1,462

1,395

1,561

1,395

1,354

1,456

1,400

1,346

1,447

1,210

1,369

1,469

1,537

1,475

1,462

1,495

122,645

119,472

119,265

121,583

117,463

120,873

121,145

118,714

118,139

121,246

121,012

121,213

114,786

113,746

122,612

115,621

120,198

118,365

120,376

120,519

119,789

120,581

120,373

120,566

120,491

120,573

122,469

117,721

120,651

123,125

108,452

110,956

117,854

114,651

113,786



Сравнивая полученные результаты, оба метода имеют небольшие отклонения.

Расположение молекулы немного изменяется в зависимости от применяемого метода.

Положительный знак электростатического потенциала отображается зелёным цветом. В области неподелённых пар на атомах азота, кислорода электрический потенциал отрицательный, что отображается красным цветом.

Изменение потенциальной энергии связи N11-С13.

Исследуемая связь между атомами N11-С13. Задаём начальные величины начальная длина связи 2,36; конечная длина связи 5,666; шаг 0,1. Проведя расчёты, изобразим график изменения потенциальной энергии связи N11-С13 от величины растяжения.

Изменение потенциальной энергии углового напряжения C19-О21-С22.

Исследуем угол между атомами C19-О21-С22. Задаём начальные величины начальный угол 70; конечный 130; шаг 10. Проведя расчёты, изобразим график изменения энергии углового напряжения при изменении величины угла.

Оптимизация геометрии и расчёт параметров молекулы методом молекулярной механики (ММ+ и MNDO метод).



Таблица 4

Длина связи или валентный угол

Данные ММ расчёта

Данные MNDO расчёта

Справочные величины

С1-С2

С2-С3

С3-С4

С4-N5

N5-C6

C6-C7

C7-C8

C41-C8

C9-C41

C5-C9

C7-C10

C10-O11

C10-O12

O12-C13

C13-C14

C14-N15

N15-C16

N15-C17

C1-C1-C3

C2-C3-C4

C3-C4-N5

C4-N5-C6

C4-N5-C9

N5-C6-C7

C6-C7-C8

C7-C8-C41

C8-C41-C9

C41-C9-N5

C9-N5-C6

C6-C7-C10

C8-C7-C10

C7-C10-C11

C7-C10-O12

C10-O12-O13

O12-C13-C14

C13-C14-N15

C14-N15-C16

C14-N15-C17

C16-N15-C17

1,53455

1,53762

1,5369

1,468

1,45685

1,5104

1,34294

1,34042

1,33896

1,35217

1,36168

1,21147

1,34976

1,41025

1,53651

1,45706

1,45297

1,4518

111,892

111,685

111,781

122,06

119,167

114,299

120,314

120,526

122,338

123,23

118,447

121,465

118,218

117,423

122,29

118,494

107,986

113,797

111,518

113,334

111,185

1,53153

1,54118

1,55042

1,47213

1,4767

1,51705

1,36463

1,44632

1,36782

1,39475

1,49687

1,22888

1,36285

1,41128

1,56218

1,47138

1,46379

1,46434

114,534

113,462

113,494

118,349

119,969

113,836

121,488

120,452

120,512

121,692

121,63

116,046

122,464

126,415

113,517

125,248

107,536

109,6

116,738

116,78

115,933

1,532

1,538

1,542

1,468

1,463

1,517

1,352

1,395

1,338

1,394

1,419

1,226

1,339

1,411

1,541

1,468

1,462

1,468

113,589

113,452

113,642

119,486

119,165

113,863

121,488

120,526

121,514

122,945

119,449

120,064

118,218

126,435

122,651

125,984

107,892

110,674

115,465

115,639

113,746



Сравнивая полученные результаты, оба метода имеют небольшие отклонения.

Расположение молекулы немного изменяется в зависимости от применяемого метода.

Положительный знак электростатического потенциала отображается зелёным цветом. В области неподелённых пар на атомах азота, кислорода электрический потенциал отрицательный, что отображается красным цветом.

Изменение потенциальной энергии связи С4–N5.

Исследуемая связь между атомами С4–N5. Задаём начальные величины начальная длина связи 0,972; конечная длина связи 2,972; шаг 0,1. Проведя расчёты, изобразим график изменения потенциальной энергии связи С4–N5 от величины растяжения.

Изменение потенциальной энергии углового напряжения C2-C3-С4.

Исследуем угол между атомами C2-C3-С4. Задаём начальные величины начальный угол 50; конечный 140; шаг 10. Проведя расчёты, изобразим график изменения энергии углового напряжения при изменении величины угла.

Оптимизация геометрии и расчёт параметров молекулы методом молекулярной механики (ММ+ и MNDO метод).

- Проведение оптимизации молекулы.

ММ+ метод.



Таблица 5

Длина связи или валентный угол

Данные ММ расчёта

Данные MNDO расчёта

Справочные величины

С1-С2

С2-С3

С3-С4

С4-N5

N5-C6

C6-C8

C8-C10

C9-C10

C7-C9

N5-C7

C10-C11

C11-O12

C11-O13

O13-C14

C14-C15

C15-N16

N16-C17

N16-C18

C1-C2-C3

C2-C3-C4

C3-C4-N5

C4-N5-C7

C4-N5-C6

N5-C7-C9

C7-C9-C10

C9-C10-C8

C6-C8-C10

C8-C6-N5

C6-N5-C7

C10-C11-O12

C9-C10-C11

C8-C10-C11

C10-C11-O13

O12-C11-O13

C11-O13-C14

C13-C14-C15

C14-C15-N16

C15-N16-C17

C15-N16-C18

1,53461

1,53838

1,53856

1,45536

1,45232

1,53586

1,53611

1,53558

1,536

1,4525

1,52029

1,20871

1,34376

1,40919

1,53584

1,45663

1,45283

1,45193

111,827

111,685

116,445

113,912

113,888

111,468

110,753

108,473

112,918

111,794

116,304

127,3

110,877

112,884

112,782

119,917

125,512

107,339

109,614

116,591

116,809

1,53147

1,54129

1,55152

1,47019

1,46727

1,54676

1,54834

1,54919

1,5454

1,46687

1,54228

1,22757

1,36181

1,41133

1,56262

1,47096

1,46422

1,4643

114,597

113,519

116,865

117,555

117,81

111,826

113,034

111,009

112,918

111,794

116,304

127,3

110,877

112,884

112,782

119,917

125,512

107,339

109,614

116,591

116,809

1,533

1,539

1,542

1,467

1,467

1,539

1,539

1,542

1,542

1,468

1,534

1,213

1,352

1,412

1,556

1,472

1,468

1,464

113,654

113,512

116,865

116,526

116,956

111,429

111,485

111,006

112,918

111,783

116,304

127,3

110,563

112,853

112,782

119,456

125,654

107,339

110,369

115,654

116,809



Сравнивая полученные результаты, оба метода имеют небольшие отклонения.

Расположение молекулы немного изменяется в зависимости от применяемого метода.

Положительный знак электростатического потенциала отображается зелёным цветом. В области неподелённых пар на атомах азота, кислорода электрический потенциал отрицательный, что отображается красным цветом.

Изменение потенциальной энергии связи С4–N5.

Исследуемая связь между атомами С4–N5. Задаём начальные величины начальная длина связи 0,97; конечная длина связи 2,97; шаг 0,1. Проведя расчёты, изобразим график изменения потенциальной энергии связи С4–N5 от величины растяжения.

Изменение потенциальной энергии углового напряжения C2-C3-С4.

Исследуем угол между атомами C2-C3-С4. Задаём начальные величины начальный угол 50; конечный 140; шаг 10. Проведя расчёты, изобразим график изменения энергии углового напряжения при изменении величины угла.



3.3 Исследование биологической активности с помощью программы PASS



В работе выполнено исследование биологической активности всех молекулярных структур с помощью программы PASS согласно методике п.2.2.

Дикаин







Если величина Pa близка к единице, то вещество может оказаться близким аналогом известных препаратов.

Базовые структуры лекарств, обладающие существенной новизной, целесообразно отбирать из массива доступных веществ соединения с Pa<0.7. Риск получения отрицательного результата в эксперименте тем больше, чем меньше величина Pa, однако и новизна такой структуры (при подтверждении прогноза в эксперименте) будет более высокой [12]. Pa Pi:

0.603 0.023 спазмолитик,

0.511 0.048 сосудорасширяющее средство,

0.405 0.015 антагонист кальциевых каналов,

0.350 0.107 антигипертензивный,

0.323 0.166 токсичный,

0.114 0.098 агонист β – адренорецепторов,

0.219 0.214 тератоген,

0.092 0.091 антагонист β – адренорецепторов.

	1.	Структура 1 (карбоксиструктура).

	

 Pa Pi:



0.591 0.025 спазмолитик,

0.367 0.095 сосудорасширяющее средство,

0.264 0.051 антагонист кальциевых каналов,

0.331 0.160 токсичный,

0.301 0.142 антигипертензивный,

0.211 0.144 диуретик,

0.233 0.195 тератоген,

0.113 0.101 агонист β – адренорецепторов,

0.092 0.090 антагонист β – адренорецепторов.

	2.	Структура 2 (адреноструктура).

	

 Pa Pi:


0.620 0.021 спазмолитик,

0.472 0.059 сосудорасширяющее средство,

0.362 0.020 антагонист кальциевых каналов,

0.218 0.041 агонист дофамина,

0.128 0.020 агонист Д2 дофамина,

0.291 0.188 токсичный,

0.144 0.041 агонист β1 – адренорецепторов,

0.139 0.043 агонист β – адренорецепторов,

0.243 0.182 тератоген,

0.237 0.211 антигипертензивный,

0.133 0.119 агонист α – адренорецепторов.

	3.	Структура 3 (никотиноструктура).

	

 



Pa Pi:

0.683 0.017 сосудорасширяющее средство,

0.548 0.031 спазмолитик,

0.326 0.026 антагонист кальциевых каналов,

0.364 0.098 антигипертензивный,

0.171 0.066 агонист дофамина.

	4.	Структура 4 (пиперидиноструктура).

 


Pa Pi:

0.680 0.015 спазмолитик,

0.537 0.042 сосудорасширяющее средство,

0.411 0.014 антагонист кальциевых каналов,

0.402 0.078 антигипертензивный,

0.233 0.051 антагонист α1 – адренорецепторов,

0.253 0.075 агонист α – адренорецепторов 

0.216 0.081 антагонист адреналина.

	После исследования биологической активности ряда веществ можно сделать вывод, что все структуры обладают новизной. Есть большая вероятность, что они будут обладать спазмолитической, сосудорасширяющей активностью и являются антагонистами кальциевых каналов.

	Суммарно характеристики биологической активности всех молекулярных структур предложены в таблице 6.

	

Таблица 6

Характерис-тика фармако-логической активности

Основ-ная струк-тура

Модифицированные структуры



Дикаин  

Структура 1

Карбокси-структура  

Структура 2

Адрено-структура  

Структура 3

Никотино-структура  

Структура 4

Пиперидиноструктура  

1. Спазмолитик

0,603 0,023

0,591 0,025

0,620 0,021

0,683 0,017

0,680 0,015

2. Сосудорасши-ряющее средство

0,511 0,048

0,367 0,095

0,472 0,059

0,548 0,031

0,537 0,042

3. Антагонист кальциевых каналов

0,405 0,015

0,264 0,051

0,362 0,020

0,326 0,026

0,411 0,014

4. Антигипер-тензивный

0,350 0,107

0,301 0,142

0,237 0,211

0,364 0,098

0,402 0,078

5. Агонист β-адренорецепторов

0,114 0,098

0,113 0,101

0,139 0,043





6. Токсичный

0,323 0,166

0,331 0,160

0,291 0,188





7. Тератоген

0,219 0,214

0,233 0,195

0,243 0,182





8. Антагонист β-адрено-рецепторов

0,092 0,091

0,092 0,09







9. Диуретик



0,211 0,144







10. агонист -адрено-рецепторов





0,144 0,041





11. Агонист α-адрено-рецепторов





0,133 0,119



0,253 0,075

12. Антагонист -адрено-рецепторов









0,233 0,051



Краткое описание позиций проявленной фармакологической активности.

1.	Спазмолитик.

Лекарственное средство, понижающее тонус и двигательную активность гладких мышц; применяют для предупреждения или устранения спазмов гладкомышечных органов.

По механизму действия спазмолитические средства делят на миотропные и нейротропные. Миотропные спазмолитические средства снижают тонус гладкомышечных органов путем прямого влияния на биохимические процессы в гладкомышечных клетках. Нейротропные спазмолитические средства оказывают спазмолитический эффект путем нарушения передачи нервных импульсов в вегетативных ганглиях или в области окончаний вегетативных нервов, стимулирующих гладкие мышцы [19].

2.	Сосудорасширяющее средство (α- и β-адреноблокаторы).

Лекарственное средство, вызывающее расширение кровеносных сосудов.

По принципу действия различают нейротропные, миотропные сосудорасширяющие средства, антагонисты кальция и сосудорасширяющие средства, влияющие на гуморальную регуляцию сосудистого тонуса.

К нейротропным сосудорасширяющим средствам относят препараты, влияющие на эффективную иннервацию сосудов [18].

3.	Антагонист кальциевых каналов.

Механизм сосудорасширяющего действия препаратов группы антагонистов кальция связывают с блокадой кальциевых каналов, что приводит к затруднению проникновения ионов кальция внутрь клетки и расслаблению гладкой мускулатуры. Из числа антагонистов кальция в медицинской практике широко используется верапамил и нифедипин, которые применяют в основном как антиангинальные средства [18, 20].

4.	Антигипертензивный.

Антигипертензивный – свойство вещества, препятствующего повышению гидростатического давления в полости организма, полых органах и сосудах.

Антигипертензивные вещества препятствуют развитию гипертензивного синдрома – симптомокомплекса, обусловленного стабильным или прогрессирующим поведением внутричерепного давления [8].

5.	Токсичный.

Токсичность – свойство вещества синтетического и природного происхождения при поступлении в организм в количестве, превышающем меру их фармакологической активности, что выражается в возникновении токсических эффектов разной направленности, интенсивности и продолжительности вплоть до развития отравления [20].

6.	Агонист β-адренорецепторов.

Агонист β-адренорецепторов – лекарственное вещество, которое прикрепляясь к β-адренорецептору, индуцирует эффективное конформационное изменение [3].

7.	Тератоген.

Тератоген – фактор, вызывающий развитие врожденных пороков [5].

8.	Антагонист β–адренорецепторов.

Антагонист β-адренорецепторов – лекарственное вещество, которое прикрепляется к β-адренорецептору, не индуцирует эффективного конфигурационного изменения.

β-адренолитики блокируют β-адренорецепторы, осуществляющие симпатическую иннервацию сердца (возбуждение) и торможение гладких мышц бронхов, желудка, некоторых сосудов, ресничной мышцы, поперечнополосатых мышц, а также регуляцию гликогенолиза и липолиза [7].

9.	Диуретик.

Диуретики (мочегонные средства) – лекарственные средства, увеличивающие выделение почками ионов натрия и воды и вызывающие в связи с этим уменьшение содержания жидкости в тканях и серозных полостях организма.

Основным и практически важным эффектом мочегонного средства является увеличение выделения ионов натрия.

Одновременно с выделением натрия мочегонные средства способствуют выделению других ионов [19].

10.	Агонист -адренорецепторов.

Агонист -адренорецепторов – вещество, которое посредством прикрепления к рецептору индуцирует эффективное конформационное изменение.

-адренорецепторы опосредуют влияние катехоламинов на сердце, гладкие мышцы желудочно-кишечного тракта и, возможно, липолитический эффект КА [7].

11.	Агонист α-адренорецепторов.

Агонист α-адренорецепторов – вещество, которое посредством прикрепления к α-адренорецептору индуцируют эффективное конформационное изменение.

α-адренорецепторы осуществляют возбуждение гладких мышц сосудов, гладких образований кожи, слизистых оболочек, органов брюшной полости, селезенки, сфинктеров желудочно-кишечного тракта и мочевого пузыря, мышцы, расширяющей зрачок и др. Сильное α-адренолитическое действие оказывают производные β-галоидоалкиламина, которые вызывают необратимую блокаду адренореактивных систем [7].

12.	Агонист α-адренорецепторов.

Агонист α-адренорецепторов – вещество, которое посредством прикрепления к α-адренорецептору индуцируют эффективное конформационное изменение.

α-адренорецепторы осуществляют возбуждение гладких мышц сосудов, гладких образований кожи, слизистых оболочек, органов брюшной полости, селезенки, сфинктеров желудочно-кишечного тракта и мочевого пузыря, мышцы, расширяющей зрачок и др. Сильное α-адренолитическое действие оказывают производные β-галоидоалкиламина, которые вызывают необратимую блокаду адренореактивных систем [7].

Из таблицы 6 видно, что порог ингибирования  практически для всех видов биологической активности незначителен, поэтому в дальнейшем сравнительный анализ фармакоактивности будем проводить по порогу активности . Одновременно приведем значения  программы PASS в условные проценты относительно базовой структуры – дикаина, принемая его характеристики за 100 %.


Таблица 7



Дикаин

Карбокси-структура

Адрено-структура

Никотино-структура

Пиперидно-структура

1. Спазмолитик

100 (0,603)

98,00 (0,591)

102,82 (0,620)

113,27 (0,693)

112,77 (0,680)

2. Сосудорасширя-ющее

100 (0,511)

71,82 (0,367)

92,37 (0,472)

107,24 (0,548)

105,09 (0,537)

3. Антагонист Ca каналов

100 (0,405)

65,19 (0,264)

89,38 (0,362)

80,49 (0,326)

101,48 (0,411)

4. Антигипертен-зивный

100 (0,350)

86,00 (0,301)

67,71 (0,237)

104,00 (0,364)

114,8 (0,402)

5. Агонист β-адренорецеп-торов

100 (0,114)

99,12 (0,113)

119,30 (0,139)





6. Токсичность

100 (0,323)

102,48 (0,331)

90,09 (0,291)





7. Тератоген

100 (0,219)

106,39 (0,233)

110,96 (0,243)





8. Антагонист β-адренорецеп-торов

100 (0,092)

100 (0,092)







9. Диуретик



(0,211)







10. Агонист -адренорецеп-торов





(0,144)





11 Агонист α-адренорецеп-торов





(0,133)



(0,253)

Антагонист -адренорецеп-торов









(0,233)



Сравнивая характеристики фармакологических структур и их соотношение, можно сделать следующие выводы.

		Чем больше показатель спазмолитических свойств, тем больше анестезирующий эффект.

		Чем меньше показатель сосудорасширяющего свойства, тем больше анестезирующий эффект.

		Чем больше показатель антагонист кальциевых каналов, тем больше анестезирующий эффект.

		Чем больше антигипертензивный показатель, тем меньше токсичность.

		Чем больше показатель сосудорасширяющего средства, тем меньше токсичность.

		Появление диуретических свойств снижает токсичность.

		Появление α, β-антагонистов адренорецепторов уменьшает токсичность.

В нашей работе для комплексной оценки анестезирующих и токсических свойств предлагается использовать интегральные показатели.

Расчет интегральных показателей проводили по формуле 1.



;



где -интегральный коэффициент анестезирующей активности. -порог активности каждого i – вида фармакологического действия, влияющего на анестезирующий эффект. -порог активности дикаина по соответствующему виду. n-число видов фармакологического действия, влияющего на анестезирующий эффект.







Для суммарной оценки токсических свойств предложен интегральный показатель токсичности.



;


где -интегральный коэффициент токсичности.  -порог активности каждого j-вида фармакологического действия, влияющего на токсичность.  -порог активности дикаина по соответствующему фармакологического действия.  n-число видов фармакологического действия, влияющего на токсичность.







Таблица 8

Интегральные коэфф-ы

Дикаин

Структура1

Карбокси-структура

Структура 2

Адрено-структура

Структура 3

Никотино-структура

Структура 4



Пиперидино-структура

Коэффициент анестези-рующей активности

0,5

0,443

0,419

0,543

0,560

Коэффициент токсичности

0,5

0,486

0,564

-

-



Из таблицы видно, что по анестезирующему эффекту исследуемые структуры можно ранжировать в следующий ряд:

Пиперидиноструктура>никотиноструктура>дикаин>

>карбоксиструктура>адреноструктура

Наглядно эффект анестезирующей активности представлен на диаграмме 1.

По токсичности исследуемые структуры располагаются в следующий ряд:

Таким образом, результаты компьютерного дизайна молекулы дикаина с целью снижения токсичности и усиления местноанестезирующего эффекта позволяют исследуемые структуры расположить в следующий ряд:

Никотиноструктура>пиперидиноструктура>адреноструктура>>дикаин>карбоксиструктура

	
4.	Экономическая часть

	

	4.1 Цель и база сравнения

	

	Несмотря на достижения современной анестезии, продолжаются поиски менее опасных средств для наркоза, разработка различных вариантов многокомпонентного избирательного наркоза, позволяющего значительно снизить дозы используемых средств, уменьшить их токсичность и побочные отрицательные влияния.

	В последнее время методы компьютерного моделирования все более входят в практику технологии создания новых синтетических лекарственных веществ. Полученные таким образом данные позволяют более целенаправленно проводить синтезы биоактивных молекул с заданными на молекулярном уровне параметрами, что значительно экономит время, материалы и силы при аналоговом поиске лекарственных веществ.

	

	4.2 Проведение работы связано с определенными видами затрат

	

	Затраты на проведение работы включают в себя:

	Расчет заработной платы работнику, выполняющему данную работу с окладом 3500 р/мес.

	Затраты на электроэнергию с ценой 1 кВт-1,6 р.

	Затраты на покупку компьютера и приобретение программы HyperChem.

	

	4.3 Заработная плата рассчитывается на 1 человека

	

	 Оклад составляет 3500 р/мес. Работа проводилась 4 месяца. Заработная плата за 4 месяца составляет 3500. 4=14000 р.

	Отчисления на социальные нужды:

	Пенсионный фонд .

	Фонд социального страхования .

	Фонд медицинского страхования .

	Фонд страхования от несчастных случаев  р.

	Итого: 5222 р.

	Основные производственные фонды

	Стоимость компьютера 20000 р.

	Стоимость компьютера Hyper Chem 30000 р.

	Итого: 50000 р.

	Амортизация

	

	

	

	4.4 Затраты на электроэнергию

	

	Цена за 1 кВт – 1,6 р.

	Затраты на энергоресурсы составили 0,1 кВт/ ч.

	Работа на компьютере составили 528 ч.

	

	

	

	Смета затрат

Статьи затрат

Стоимость, руб

Информационная программа HyperChem

30000

Заработная плата

19222

Амортизация

1920

Затраты на электроэнергию

844,8

ИТОГО

51986

	
Список литературы



				Поройков В.В. Компьютерное предсказание биологической активности веществ: пределы возможного. Химия в России, 1999, № 2, 8-12. 

				Кнунянц И. Л. Химическая энциклопедия. Издательство “Советская энциклопедия” Москва, 1988.

				Кукес В. Г., Стародубцева А. К. Фармакология и фармакотерапия. - М.: ГЭОТАР – МЕД, 2004.

					Беликов В. Г. Фармацевтическая химия. – М.: Высшая школа, 1985

				Харкевич Д. А. Фармакология, четвертое издание, Москва, 1993.

				Солдотенков А. Т., Колядина Н. М., Шендрик И. В. Основы органической химии лекарственных веществ. – М.: МИН, 2003.

				Аляутдин Р. Н. Фармакология. – учебник для вузов, Москва, 2004.

					Ланса Л., Лейси Ч., Голдман. М. Фармакологический справочник, Москва, 2000 г.

					Поройков В.В., Филимонов Д.А. Компьютерный прогноз биологической активности химических соединений как основа для поиска и оптимизации базовых структур новых лекарств. В сб.: Азотистые гетероциклы и алкалоиды. Москва: Иридиум-пресс, 2001, т.1, с.123-129. 

					Poroikov V.V., Filimonov D.A., Borodina Yu.V., Lagunin A.A., Kos A. Robustness of biological activity spectra predicting by computer program PASS for non-congeneric sets of chemical compounds. J. Chem. Inform. Comput. Sci., 2000, 40 (6), 1349-1355. 

					Anzali S., Barnickel G., Cezanne B., Krug M., Filimonov D., Poroikov V. Discriminating between drugs and nondrugs by Prediction of Activity Spectra for Substances (PASS). J. Med. Chem., 2001, 4 (15), 2432-2437. 

					Лагунин А.А., Филимонов Д.А., Поройков В.В. Компьютерный поиск потенциальных антигипертензивных соединений комбинированного действия. Хим.-фарм. журн., 2001, 35 (7), 28-34. 

					Filimonov D., Poroikov V., Borodina Yu., Gloriozova T. Chemical similarity assessment through multilevel neighborhoods of atoms: definition and comparison with the other descriptors. J. Chem. Inf. Comput. Sci., 1999, 39 (4), 666-670. 

					Lagunin A., Stepanchikova A., Filimonov D., Poroikov V. PASS: prediction of activity spectra for biologically active substances. Bioinformatics, 2000, 16 (8), 747-748. 

					Poroikov V., Akimov D., Shabelnikova E., Filimonov D. Top 200 medicines: can new actions be discovered through computer-aided prediction? SAR and QSAR in Environmental Research, 2001, 12 (4), 327-344. 

					Poroikov V., Filimonov D. Computer-aided prediction of biological activity spectra. Application for finding and optimization of new leads. Rational Approaches to Drug Design, Eds. H.-D. Holtje, W.Sippl, Prous Science, Barcelona, 2001, p.403-407.

					Кудрин А. Н. Фармакология, Москва “Медицина”, 1991.

					Лоуренс Д. Р., Беннетт П. Н. Браун М. Дж. Фармакология. Издание второе. Москва, 2002.

					Кудрин А. Н. Фармакология. – М.: Медицина, 2001.

					Лоуренс Д. Р., Беннетт П. Н. Фармакология Том 1. Москва, 1993.