Биоизостеризм

ОГЛАВЛЕНИЕ

ВВЕДЕНИЕ.....................................................................................................3

  1. Развитие концепции биоизостеризма.....................................................4
    1. Работы Ирвинга Ленгмюр................................................................4
    2. Работы Г. Гримма...............................................................................5
    3. Работы Ганса Эрленмейера...............................................................6
    4. Работы Корвина Ганча......................................................................8
  2. Классификация биоизостеров................................................................10
  3. Биоизостеризм пиретрина .....................................................................13
  4. Биоизостеризм анальгетиков на основе морфина..............................16
  5. Биоизостеризм цефалоспаринов............................................................17

ЗАКЛЮЧЕНИЕ............................................................................................19

СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ......................................20

 

ВВЕДЕНИЕ

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

 

  1. Развитие биоизостеризма
    1. Работы Ирвинга Ленгмюра

История возникновения концепции биоизостеризма началась в 1919 г., когда американский химик Ирвинг Ленгмюр выдвинул идею об изостеризме [2]. Развивая предложенную им ранее октетную теорию валентности, он заключил, что в соединениях, имеющих одинаковое количество атомов и одно и то же общеечисло электронов, расположение последних также одинаково. При этом Ленгмюр обнаружил, что такие молекулы имеют удивительно похожие физико-химические свойства. В работе “Изоморфизм, изостеризм и ковалентность” Ленгмюр сравнил молекулы закиси азота и диоксида углерода и показал, что более десяти их физико-химических констант (включая вязкость, теплопроводность, растворимость в воде и т. п.) оказались чрезвычайно близкими. Такие соединения Ленгмюр назвал “изостерическими соединениями или изостерами”, подчеркнув, что это понятие применимо не только к молекулам, но и к радикалам и к группам атомов, имеющим обобществленные пары электронов. Таким образом, по Ленгмюру, изостеры – это молекулы или ионы, содержащие одинаковое число атомов, а также имеющие одинаковое количество и расположение электронов. Группы таких изостеров, предложенные в таблице 1.1

Таблица 1.1 − Группы изостеров по Ленгмюру [2]

№ группы

Изостер

№ группы

Изостер

1

H-, He, Li+

12

NO2-, O3

O2-, F-, Ne, Na+, Mg2+, Al3+

13

HF, OH-

3

S2-, Cl-, Ar, K+, Ca2+

14

ClO4-, SO42-, PO43-

4

Cu+, Zn2+

15

ClO3-, SO32-, PO33-

5

Br-, Kr, Rb+, Sr2+

16

SO3, PO3-

6

Ag+, Cd2+

17

S2O62-, P2O64-

7

I-, Xe, Cs+, Ba2+

18

S2O72-, P2O74-

8

N2, CO, CN-

19

SiH4, PH4+

9

CH4, NH4+

20

MnO4-, CrO42-

10

CO2, N2O, N3-, CNO-

21

SeO42-, AsO43-

11

NO3-, CO32-

   



Очевидно, что только изостеры, имеющие одинаковый общий заряд, обладают сходными физико-химическими свойствами. В таблице 1.1 представлены только три пары таких соединений: N2, CO (8 группа); CO2, N2O (10 группа) и N3–, CNO– (10 группа). Сходство в свойствах для соединений первых двух групп было показано Ленгмюром в [2]. Сравнение свойств цианатов и азидов в то время было затруднено из-за отсутствия данных об их физико-химических свойствах, но по предположению Ленгмюра последние должны были быть очень близки. Кроме того, на основании предложенной октетной теории и концепции изостеризма Ленгмюр предположил, что схожими свойствами должны обладать и неизвестные в то время молекулы диазометана (H2C=N=N) и кетена (H2C=C=O), что подтвердилось впоследствии, когда эти соединения были синтезированы.

Распределение изостеров на группы (таблица 1.1) позволило находить определенные соотношения между ними. Ленгмюр заключил, что если аргон является изостером иона калия (3 группа), а метан – изостером иона аммония (9 группа), то свойства ионов калия и аммония должны быть похожи, так как аргон и метан имеют очень близкие физико-химические свойства. Такие соотношения, по мнению Ленгмюра, делали теорию изостеризма чрезвычайно плодотворной.

Важно подчеркнуть, что говоря о близости характеристик изостерических соединений, Ленгмюр имел в виду в основном их физико-химические свойства (в том числе, кристаллические формы), но не принимал во внимание их физиологическую активность*(Через двадцать лет после работы Ленгмюра была установлена идентичность свойств углекислого газа и закиси азота (10 группа изостеров) как обратимых анестетиков для некоторых организмов [3].).

    1. Работы Г. Гримма

Дальнейшее развитие концепция изостеризма получила в работах Г. Гримма [4–7]. Гримм вывел правило (получившее название закона гидридного замещения) для описания сходства между группами соединений, которые имеют одинаковое количество электронов, но различное количество атомов. Закон гидридного замещения заключается в следующем: “Атомы в любом месте периодической системы в пределах четырех клеток до инертного газа при присоединении к ним от одного до четырех атомов водорода изменяют свои свойства таким образом, что образующиеся комбинации ведут себя как псевдоатомы, аналогичные элементам в группах, находящихся справа от них от одной до четырех клеток соответственно” (таблица 1.2) [4, 7]. Согласно этому закону, каждый вертикальный столбец в таблице 1.2 объединяет изостерические комбинации, т.е. свойства однотипных соединений, содержащих, например фрагменты –О–, –NH–, –CH2– (VI группа), должны быть похожими.

Таблица 1.2 − Некоторые группы изостеров по Гримму [4, 7]

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

    1. Работы Ганса Эрленмейера

Дальнейшее развитие концепция изостерической замены получила в многочисленных работах Ганса Эрленмейера [8–14], который расширил классификацию Ленгмюра и Гримма и определил изостеры как “атомы, ионы или молекулы, в которых наружные электронные оболочки могут считаться идентичными” [8] (таблица 1.3). Основой теории Эрленмейера был тот факт, что сходства в свойствах элементов, находящихся в вертикальных группах периодической таблицы, могут быть отнесены к образуемым из этих элементов молекулам или молекулярным фрагментам. Некоторым подтверждением справедливости вышесказанного было то, что пары, такие, как бензол и тиофен, тиофен и фуран, и даже бензол и пиридин обладали сходством по многим физическим и химическим характеристикам [9].

Важнейшей заслугой Г. Эрленмейера, является то, что при сопоставлении свойств изостеров он сравнивал не только их физико-химические свойства и реакционную способность [9, 11], но и физиологическую (фармакологическую) активность [8, 10, 13, 14]. Так, исследуя свойства ряда изостеров (–O–, –NH–, –CH2–), предложенных Гриммом (таблица 1.2, VI группа), Эрленмейер сравнил активность соединений 1–3 в иммунологических тестах и показал, что она практически одинакова [8, 14].

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

 Таблица 1.3 − Группы изостеров, полученные на основании идентичности количества электронов на внешних оболочках

Количество электронов во внешнем слое

4

5

6

7

8

N+

P

S

Cl

ClH

P+

As

Se

Br

BrH

S+

Sb

Te

I

IH

As+

 

PH

SH

SH2

Sb+

   

PH2

PH3



 

 

 

 

 

 

 

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

Термин “биоизостеризм” был предложен в 1951 г. Г. Фридманом, называвшим биоизостерами все соединения (атомы и молекулы), которые “удовлетворяют самому широкому определению изостеров и имеют тот же тип биологической активности” [18]. При этом Фридман допускал, что активность может быть и антагонистической, а самое важное заключается в том, что биоизостеры воздействуют на одну и ту же биологическую мишень в организме, и, следовательно, имеют сопоставляемые физиологические свойства. Впоследствии некоторые авторы [19, 20] попытались ввести понятие о так называемых “частичных биоизостерических группах”, т.е. тех, которые превращают агонистов в антагонистов. Однако эта идея оказалась не слишком удачной, так как замена, являющаяся биоизостерической в одной серии соединений, не обязательно будет таковой в другой серии. Таким образом, “антагонистическая” группа одной молекулы изменяет на обратную активность “агонистической” группы другой молекулы, только если агонистические группы в обеих сериях выполняют одну и ту же функцию. Отметим, что в настоящее время термин “частичные биоизостерические группы” не употребляется.

Для неклассических биоизостеров, как правило, не проводили деления на подтипы, однако в некоторых случаях в отдельную группу выделяли циклические соединения, заменяющие нециклические структуры [15]. Как видно из таблиц 2.1 и 2.2, классическая изостерическая замена предполагает замещение подобного на подобное, принимая во внимание количество атомов, валентность, степень ненасыщенности и ароматичности, и становится биоизостерической, если при этом сохраняется биологическая активность. Почему же сохраняют активность неклассические биоизостеры?

    1. Работы Корвина Ганча

Поиск ответа на этот вопрос обусловил дальнейшее развитие концепции биоизостеризма в работах Корвина Ганча – одного из основателей методологии QSAR (количественные соотношения структура – активность) [22]. Ганч определил биоизостеры как соединения, “вызывающие идентичный биохимический или фармакологический ответ в стандартной тест-системе”. Этой системой может быть мембрана клетки, фермент или рецептор, а также подопытное животное, человек и т.п. Ганч использовал также термин “частичные биоизостеры” (практически не употребляющийся в настоящее время), т. е. соединения, вызывающие одинаковый качественный (но не количественный) биохимический или фармакологический ответ.

В своей работе [23] Ганч попытался объяснить причину одинаковой активности неклассических биоизостеров, не “выискивая сходство в их химических формулах”, а используя предложенный им ранее математический подход, описывающий соотношения структура–активность в терминах трех основных различий в свойствах физиологически активных соединений: стерических, электронных и гидрофобных. Корреляционное уравнение, предложенное Ганчем для определенного вида активности в данном ряду химических соединений, имеет вид:

lg1/c = ao + a1logP – a2(logP)2 + a3σ + a4Es, (1)

где с – любая экспериментальная величина, характеризующая биологическую активность, Р – коэффициент распределения вещества между липидной и водной фазами (обычно в системе октанол–вода), параметры σ (константа Гаммета) и Es (константа Тафта) отражают электронные и стерические влияния заместителей, ai – параметры, полученные при обработке экспериментальных данных методом наименьших квадратов.

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

lg1/c = 3,3 (±0,1) + 0,84 (±0,07)logP, (2)

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

(r2 = 0,96) показывает, что различия данных по активности на 96% объясняются разной липофильностью тестируемых веществ. Таким образом, для этого конкретного случая биоизостерические соединения будут “изолипофильными”, как, например, структуры, представленные в таблице 1.4.

Ганч отмечал, что “традиционный путь поиска общих структурных черт” в парах соединений типа представленных в таблице 1.4, оказывается бесполезен. “Это не означает, – писал он, – что сходства вообще нет. Для ротовой ганглии указанные ионы кажутся одинаковыми, и только наши искусственные формулы и номенклатура, не учитывающая их липофильный характер, делают их совершенно различными” [23]. Традиционный подход, по мнению Ганча, ведет лишь к увеличению списка “неклассических биоизостеров”, в то время как предложенное им более общее определение биоизостеризма, учитывающее гидрофобные, электронные и стерические эффекты атомов или групп атомов, оказывается существенно более плодотворным (В своей работе Ганч упомянул также о необходимости учета еще одного важного фактора, а именно метаболизма, так как большая подверженность метаболической атаке какого-либо аналога из серии веществ с близкими основными характеристиками может “свести на нет его биоизостеричность”) [23]. Подход, предложенный Ганчем, широко используется в настоящее время для дизайна физиологически активных веществ.

Таблица 1.4 − Анионы, вызывающие увеличение на 20 мВ мембранного потенциала ротовой ганглии моллюска (из [23])

  1. Классификация биоизостеров

Традиционно биоизостеры разделяют на две группы: классические изостеры, имеющие приблизительно одинаковый размер, форму и конфигурацию внешнего электронного слоя ([20, 21], таблица 2.1) и неклассические биоизостеры (таблица 2.2), т. е. атомы или группы, не обладающие характеристиками классических изостеров, однако имеющие сходную биологическую активность.

Фактически классические (био)изостеры соответствовали определению Эрленмейера для изостеров. Традиционно их разделили на несколько категорий:

 

Таблица 2.1 − Примеры классических изостеров

Одновалентные атомы или группы

–CH3; –NH2; –OH; –F; –Cl

–PH2; –SH

–Br; -изопропил

–I; -третбутил

Бивалентные атомы или группы

–CH2–; –NH–; –O–;–S–; –Se–

–COCH2–; –CONH–; –COO–; –COS–

Заместители сложноэфирной группы

–S(OOH)N(R) –; –C(R)=NO–;

1,2-оксазол-3,5-диил,

1,3-оксазол-2,5-диил,

1,3,4-оксадиазол-2,5-диил;

Трехвалентные атомы или группы

=C(H) –; =N–;

=P–; As–;

Четырехвалентные атомы

=C=; =Si=;

=N+=; =P+=

Заместители ГАМК

β-циклогесилиден ГАМК;

пирацетам;

1,2,3,6-тетрагидропиридин-4-карбоновая  кислота;

Заместители пептидных групп

–NH-CH2C(O)NHCH2C(O)–;

–NH-CH2C(S)NHCH2C(O)–;

–NH-CH2C(O)NHC(R)2C(O)–;

–NH-CH2C(O)NHN(R)C(O)–;

–NH-CH2C(O)OCH2C(O)–;

–NH-CH2CH=CHCH2C(O)–;

Заместители спейсера –(CH2)3

1,4-фенилиден

Эквиваленты кольцевых радикалов

 

Заместители циклопентильного радикала

–СН2–;

Заместители тетрагидрофурильного радикала

–O–; –-СН2–;

Заместители пирролидинильного радикала

–NH–; –СН2–;

Заместители тиофенильного радикала

–S–;

Заместители фенильного радикала

–CH=CH–; –CH=; –NH=; пиридил

Заместители пиридильного радикала

–CH=CH–; –CH=; –NH=; фенил

Эквивалентные кольца

 
 
 

Продолжение таблицы 2.1

Заместители бензольного ядра

пиридильное, пиримидиновое, тиофеновое, пиррольное, фурановое, изотиазольное, изоксазольное, пиразольное

Заместители пиридинового ядра

нитрофенильное, тиофеновое, пиррольное, фурановое, пиридиниевое, R3N+C6H4;

Заместители фенола

фторбензол, индол, имидазол

Заместители1, 2-дигидроксифенола (катехола)

бензимидазол, бензопиразол,

3-гидрокси-1-R-4-пиридон,

N-гидрокси-2-пиридон;

Заместители индола

фенол, 4,5-бензоизоксазол, 4,5-бензопиразол, 1-азаиндолизин, хинолин;

Заместители карбазола

9-замещённый флуорен, дифенилметан


 

Таблица 2.2 − Примеры некоторых “неклассических” биоизостеров

  1. Биоизостеризм пиретрина

Группа пиретроидных инсектицидов сформировалась на основе упрощения структуры-лидера и введения разнообразных, как оказалось биоизостерных, групп. Систематический подход в постепенном изменении его структуры позволил получить несколько поколений синтетических пиретроидных инсектицидов. Первое поколение появилось после замены в спиртовой компоненте природного пиретрина I пентадиенильного фрагмента сначала на химически более устойчивый пропенильный (аллетрин, 2). Позднее была осуществлена более принципиальная «хирургическая» операция: вместо вырезанной циклопентеноновой части молекулы инсектицида провели химическую имплантацию 2-бензилфурилметильного блока (ресметрин, 3), а затем и его заменили на фенилоксбензильную группировку (фенотрин, 4). В результате появилось второе поколение синтетических пиретроидов с повышенной фотоустойчивостью. Следующим успешным этапом на пути химической модификации оказалось варьирование заместителей в циклопропильном кольце, наряду с введением цианогруппы в спиртовую компоненту. Это привело к появлению препаратов третьего (соединения 5) и четвертого (соединения 6 и 7) поколения. Появление следующего поколения инсектицидов, имеющих большую активность, более широкий спектр действия и повышенную устойчивость, также обязано аналоговому синтезу. Новизна их строения связана, во-первых, с раскрытием циклопропанового кольца (флуцитринат, 8), во-вторых, с переходом от сложных эфиров к простым (препараты 9). Затем упразднили диалкилоксидный спейсер, связывающий две арильные группировки, заменив его на алкен (вещество 10), и, наконец, ввели атом кремния в диарилалкильную цепочку (силафлуофен 11). Последние новые группы высокоэффективных пестицидов (8-11) имеют лишь отдаленное структурное сходство (соединения 8,9) или не имеют его вовсе (соединения 10,11) со своим природным родоначальником (1). Тем не менее, все инсектициды (1-11) обладают близким механизмом физиологического контактного действия на насекомых, которое проявляется в блокировании передачи нервных импульсов[26].

Рисунок 3.1 – Модификация молекулы пиретрина[26]

 

  1. Биоизостеризм анальгетиков на основе морфина

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

Рисунок 4.1 – Структурные формулы оксилорфана и буторфанола

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

Рисунок 4.2 – Изменение строения морфина

В результате удаления циклогексанового кольца были синтезированы вещества группы бензоморфанов (пентазоцин и др.), которые обладают меньшими анальгетическими и наркотическими эффектами по сравнению с морфином. Затем в результате следующего этапа умозрительного и экспериментального дизайна − «вырезания» бензильной СН2-группы,− был осуществлён переход к большому семейству мощных анальгетиков 4-арилпиперидинового ряда, обладающих выраженным наркотическим действием[26].

  1. Биоизостеризм цефалоспоринов

Цефалоспорины I покаления

Рисунок 5.1 – Структурная формула цефалоридина

Цефалоспорины II покаления

Рисунок 5.2 – Структурная формула цефуроксима

Цефалоспорины III покаления

Рисунок 5.3 – Структурная формула цефтриаксона

Цефалоспорины IV покаления

 Рисунок 5.4 – Структурная формула цефпирама

 

ЗАКЛЮЧЕНИЕ

  1. Прослежена эволюция понятия “биоизостеризм”, изменившего за пять десятилетий свое содержание от предложенного Ленгмюром термина “изостеризма” до максимально обобщенного определения Ганча, принятого и в настоящее время.
  2. Рассмотрена классификация биоизостеров на основе физических и биологических параметров, согласно которой биоизостеры разделили на классические и неклассические
  3. Рассмотрена разработка лекарственных средств на основе цефалоспоринового ряда антибиотиков, анальгетиков на основе морфина, инсектицида пиретрина. 
    СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ

1. Зефирова О.Н., Зефиров Н.С. // Вестн. Моск. ун-та. Сер.2. Хи-

мия. 2000. 41. С. 103.

2. Langmuir I. // JACS. 1919. 41. P. 1543.

3. Seifriz W. // Science. 1948. 107. P. 15.

4. Grimm H.G. // Z. Electrochem. 1925. 31. S. 474.

5. Grimm H.G. // Z. Electrochem. 1928. 34. S. 430.

6. Grimm H.G. // Z. Electrochem. 1934. 47. S. 594.

7. Grimm H.G. // Naturwissenschaften. 1929. 17. S. 557.

8. Erlenmeyer H., Leo M. // Helv. Chim. Acta. 1932. 15. S. 1171.

9. Erlenmeyer H., Leo M. // Helv. Chim. Acta. 1933. 16. S. 1381.

10. Erlenmeyer H., Berger E., Leo M. // Helv. Chim. Acta. 1933. 16.

S. 733.

11. Erlenmeyer H., Leo M. // Helv. Chim. Acta. 1933. 16. S. 57.

12. Erlenmeyer H., Willi E. // Helv. Chim. Acta. 1935. 18. S. 740.

13. Erlenmeyer H., Berger E. // Biochemische Z. 1933. 262. S. 196.

14. Erlenmeyer H., Berger E. // Biochemische Z. 1932. 252. S. 22.

15. Patani C. A., LaVoie E. J. // Chem. Rev. 1996. 96. P. 3147.

16. Schatz V.B. // Medicinal Chemistry, 2nd Edn., ed. Burger A., Wiley-

Interscience. N.Y., 1960.

17. Thornber C. W. // Chem. Soc. Rev. 1979. 8. P. 563.

18. Friedman H.L. // National Academy of Sciences–National

Research Council Publication N. 206. Washington, 1951. P. 295.

19. Ariens E. J. // Drug Design. Ed. E. J. Ariens. N.Y., 1971.

20. Korolkovas A. // Essentials of Molecular Pharmacology:

Background for Drug Design. N.Y., 1970.

21. Burger A. // Medicinal Chemistry, 3rd Edn., ed. Burger A., Wiley-

Interscience. N.Y., 1970.

22. Hansch C., Fujita T. // J. Am. Chem. Soc. 1964. 86. Р. 1616.

23. Hansch C. // Intra-Science Chem. Rep. 1974. 8. P. 17.

24. Lipinski C.A. // Annu. Rep. Med. Chem. 1986. 21. P. 283.

25. Burger A. // Prog. Drug. Res. 1991. 37. P. 287.

26. Основы дизайна и синтеза  лекарственных соединений / Солдатенков А.Т. – Москва, РУДН, 2012. – 203 с.

27. Фармакология/ Д. А. Харкевич. – Москва, Издательская группа «Геотар-Медиа», 2008. – 440 с.

 


Биоизостеризм