Введение

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

Создание новых лекарственных веществ включает 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, в некоторых случаях метод дает серьезные сбои. Это в первую очередь касается расчетов молекул с водородными связями, барьеров внутреннего вращения в p-сопряженных системах и расчетов четырехчленных циклов. Поэтому в рамках метода 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,07A

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.