Application Notes
Does it possible to explain MAC (minimum alveolar concentration) of inhaled anesthetics from molecular structure? The minimum alveolar concentration is one of the index to count inhaled anesthetics strength of anesthesia.
Do it yourself (DIY): It was found that the solubility parameter has played a surprisingly large role to estimate MAC. This treatment is not only the position in three-dimensional space of Hansen parameters, but it is necessary to use the results of YMB calculations. So, I will explain in the DIY corner. Try it by yourself!
2012.5.4
HSPiP Team Senior Developer, Dr. Hiroshi Yamamoto
see article on Wiki Anesthetics
Inhalation agents (inhaled anesthetics) are absorbed from the lungs and excretion is very early, so they are widely used.
Once, ether, laughing gas (N2O), Halothane were used.
To evaluate of the strength of these inhalation anesthetic, MAC (minimum alveolar concentration) is used. MAC is the lung inhalation anesthetic concentration that 50% of patients do not show any response to skin stimulation. The effects of inhalation anesthetics is stronger MAC is small. Partition coefficient blood/gases is the ratio of the concentration of inhalation anesthetic in the blood on the concentration of inhaled anesthetics has reached a state of equilibrium. This is the indicator of recovery from anesthesia with the introduction of inhaled anesthetics.
Ether is a flammable, take a long time to deployment and awakening, and now is not used.
Laughing gas is a colorless and odorless gas, soluble in water. Analgesic action is strong, so weak anesthetic action by itself. Laughing gas is not used now.
Halothane is a volatile anesthetic and not flammable. Anesthetic action is quite strong. Analgesic, muscle relaxant effect is so weak, that combination of muscle relaxants and nitrous. Because there is a vasodilatory and myocardial depressant effect, easy to drop in blood pressure. Because it enhances the adrenaline sensitivity of myocardial conduction system, arrhythmia-prone. In some cases, may cause liver dysfunction or malignant hyperthermia.
There is little isoflurane myocardial depression. Liver failure is low because they are not metabolized in the body. introduced from halothane, arousal is faster. Has the effect of an increase in cerebral blood flow. There is a dose-dependent respiratory depressant effect.
Sevoflurane has been used most often in our country. Fast to deploy and arousal. There is a dose-dependent respiratory depressant effect. There is a strong analgesic effect. Is degraded by soda lime adsorbent and Bararaimu CO2, occurs in compound A nephrotoxicity.
I found one interesting literature that summarize MAC of inhaled anesthetics.
Polyhalogenated Methyl Ethyl Ethers: Solubilities and Anesthetic Properties
International Anesthesia Research Society Anesth Analg 1999;88:1161–7
For more information please read the above literature.
I summarized the molecules below table.
Hcode | Compound | name | CAS | log(KMAC) |
22329 | CCIF2OCCIFCF3 | |||
21997 | CCIF2OCF2CCIF2 | (Chlorodifluoromethyl)(2-chloro-1,1,2,2-tetrafluoroethyl) ether | ||
21998 | CCIF2OCCI2CF3 | 1,1-Dichloro-1-[chloro(difluoro)methoxy]-2,2,2-trifluoroethane | 32778-09-9 | 2.120573931 |
21999 | CCIF2OCF2CCI2 F | 1,1-Dichloro-2-[chloro(difluoro)methoxy]-1,2,2-trifluoroethane | 37136-24-6 | 2.269512944 |
22330 | CCI2FOCF2CCIF2 | 1-chloro-2-[dichloro(fluoro)methoxy]-1,1,2,2-tetrafluoro-ethane | 2.260071388 | |
22331 | CCIF2OCFHCF3 | 2-[Chloro(difluoro)methoxy]-1,1,1,2-tetrafluoroethane | 2.463892989 | |
22332 | CF2HOCCIFCF3 | 1-chloro-1-(difluoromethoxy)-1,2,2,2-tetrafluoro-ethane | 2.507855872 | |
22333 | CF2HOCF2CCIF2 | (1,1,2,2-Tetrafluoro-2-chloroethyl)(difluoromethyl) ether、(Difluoromethyl)(2-chlorotetrafluoroethyl) ether | 2.781036939 | |
21984 | CCIF2OCH2CF3 | 2-[Chloro(difluoro)methoxy]-1,1,1-trifluoroethane | 33018-78-9 | 2.457881897 |
3571 | CF2HOCCIHCF3 | isofluran | 26675-46-7 | 1.161368002 |
3572 | CF2HOCF2CCIFH | Enflurane | 13838-16-9 | 1.342422681 |
22334 | CCIF2OCCIHCF3 | 2-Chloro-2-[chloro(difluoro)methoxy]-1,1,1-trifluoroethane | 1.686636269 | |
22335 | CF2HOCCI2CF3 | 1,1-Dichloro-1-(difluoromethoxy)-2,2,2-trifluoroethane | 32778-07-7 | 1.992111488 |
22336 | CCIF2OCF2CCIFH | 1.477121255 | ||
22337 | CF2HOCF2CFCI2 | 1,1-Dichloro-2-(difluoromethoxy)-1,2,2-trifluoroethane | 1.954724791 | |
21979 | CH3OCF2CHCI2 | Methoxyflurane | 76-38-0 | 0.431363764 |
21994 | CF2HOCF2CF3 | 1-(Difluoromethoxy)-1,1,2,2,2-pentafluoroethane | 3.751279104 | |
3268 | CF3OCFHCF3 | 1,2,2,2-Tetrafluoroethyl trifluoromethyl ether | 2356-62-9 | 3.292256071 |
3106 | CF2HOCFHCF3 | Desflurane | 57041-67-5 | 1.892094603 |
3272 | CF2HOCH2CF3 | 2,2,2-trifluorethyl difluoromethyl ether | 1885-48-9 | 2.042181595 |
21980 | CFH2OCFHCF3 | 1,1,1,2-Tetrafluoro-2-(fluoromethoxy)ethane | 1.666517981 | |
3273 | CFH2OCF2CF2 H | fluoromethyl 1,1,2,2-tetrafluoroethyl ether | 37031-31-5 | 1.624282096 |
3574 | CH3OCF2CCIFH | 2-Chloro-1,1,2-trifluoroethyl methyl ether | 425-87-6 | 1.201397124 |
22338 | CF2HOCBrCICF3 | 1.178976947 | ||
22339 | CF2HOCF2CBrCIF | 1.176091259 | ||
21983 | CF2HOCBrHCF3 | 1-Bromo-2,2,2-trifluoroethyl difluoromethyl ether | 32778-10-2 | 0.716003344 |
21978 | CH3OCF2CBrFH | ROFLURANE | 679-90-3 | 0.838849091 |
255 | EtOEt | EtOEt | 60-29-7 | 1.278753601 |
3546 | CF3CHBrCl | Halothane | 151-67-7 | 0.875061263 |
3281 | CH2FOCH(CF3)2 | Sevoflurane | 28523-86-6 | 1.77815125 |
649 | Cl2C=CClH | TriChloroethylene | 79-01-6 | 0.230448921 |
YMB (Properties Estimation function) has been implemented into HSPiP ( Hansen Solubility Parameters in Practice) and with using this function, I calculated all these molecules. It is possible to estimate the various properties of molecules from Smiles structure notation. Then I performed QSAR (Quantitative Structure Activity Relationship) analysis. (I used log(MAC*1000)=log(KMAC) as the target property)
At first, I checked the rough correlation log(KMAC) with physical properties of each compound. For example, I made the graph log(KMAC) to boiling point of compounds.
There is some tendency between log(KMAC) and boiling point. If I assignd log(KMAC)<2.0 as strong inhaled anesthetics, then boiling point range of 292K-375K compounds are the answer. We can easily know this result is strange. Because almost all common solvents are in this range.
Then what other additional factors determin anesthesia strength?
One factor is Hansen solubility parameters, dD (dispersion term) as show below graph.
If I assignd log(KMAC)<2.0 as strong inhaled anesthetics, then the value of dD is between 13.3-18.
In addition, molecular volume is very important. If I assignd log(KMAC)<2.0 as strong inhaled anesthetics, then the adequate molecular volume is 90-140.
Then I made QSAR model that estimate log(KMAC) from boiling point, dD and molecular volume. I show the result below.
When I made this scheme, I use only halogenated ether data (Blue diamond), and using this scheme, I predict diethyl ether, Halothane, Sevoflurane, Trichloroethylene (Red square). The result is fairly good.
Then I search HSPiP database with the condition 292K<boiling point<375K, 13.3< dD < 18, 90< molecular volume <140. HSPiP data base have over 10,000 compounds but there is only 373 compounds that meat this ranges.
Why there are so few molecules that meat this range?
Larger boling point and smaller molecular volume is completely confricted. Halogenated compounds that include Cl atom or Br, are selected because Cl or Br atom size are not so large but they increase dD and Boiling point. F atom adjust volume, dD and boiling point.
YMB have the function that estimate Antoine parameters to calculate vapor pressure. I checked these Antoine parameters applicability for thse haloganated compounds. In this paper, they listed vapor pressure of these compounds. The results is fairly good as shown in below chart even YMB calculate vapor pressure from only chemical structure,
If we have the experimental boiling point for these compounds, the YMB calculation will increase accuracy.
These anesthetics are absorbed from the lungs, absorbed into blood then circulate over the body. I am just the Polymer Chemist and do not know how these chemicals work in body, but if delivered anesthetics dissolve into the nerve to act, then solubility to nerve become very important. HSPiP have he function to determine HSP of target that absorb anesthetics. The basic concept of this thought is “likes dissolve likes”. When we say “likes” we means “the HSP vector similarity”.
We calculate vectors similarity with next scheme.
HSP distance = (4*(dDsolv - dDtarget)^2 + (dPsolv -dPtarget)^2 + (dHsolv- dHtarget)^2)^0.5
If we apply Sphere Calc. Data that implemented in Sphere GA search, program will determine Target HSP as [12.9, 3.1, 2.3]. The CCIF2OCFHCF3 is shortest HSP distance. What is this region?
By chance, I found interesting patent that Pfizer took. JP 2000-26367
They made adamantane that have 3 fluorine atoms in R2-R4 position. They calmed this fluorinated adamantane deliver drugs into central nervous system (CNS) of patients who suffer Alzheimer's or Parkinson's disease. So I try to calculate HSP of fluorinated adamantane with YMB function.
The fluorinated adamantane’s HSP become [13.9, 0, 3.2] . This HSP is very similar to CCIF2OCFHCF3 [12.9, 3.1, 2.3] or CF3OCFHCF3 [13.2, 3.3, 2.7].
Is there anyone who want to make the drug attached these structure instead of fluorinated adamantane?
Let's see these materials in the Hansen 3D space.
If you are using HTML5 enable browser such as Chrome or Safari,FireFox4 (IE9 have problems), you can see the Canvas. Drag=Rotate, Drag+Shift key=Mag/Red, Drag+alt key(option key for Mac)=Translate. If you pick solvent, solvent name will appear.
Green Large sphere show CNS, Cyan large Sphere show fluorinated adamantane. Red small spheres show haloganate ethers and Blue small spheres show traditional anesthetics.
The other day, I read newspaper that said Alzheimer's disease start some protein dissolve nerve and spread whole nerves. If this protein’s HSP is near to [12.9, 3.1, 2.3], what we can do? In electrone beam etching field, they made adamanatane attached (Meth)acrylate as regist polymer, If we made fluorine modefied adamanatane poly-(meth)acrylate, this polymer catch this protein instead of nerve like oil absorbent paper?
I searched HSPiP database with the condition of 12.0 <dD <14, dP <4, dH is <4 Boiling Point> 350K, only 69 compounds exist. These chemicals also become good lead compound I think.