Studied bile acids
1. Introduction
Beside of separation and identification of chemical compounds, chromatography can be used to obtain molecular parameters that reflect their structural characteristics – molecular descriptors. Most often it is a parameter of hydrophobicity (lipophilicity) of molecules (ions) which is obtained by using retention parameters of reversed phase liquid chromatography or thin layer chromatography of high resolutions (RPHPLC and RPTLC) [1,2]. Hydrophobicity of molecules is an important feature in medical chemistry, and arbitrarily connected with the logarithm of the solute partition (distribution) coefficient log
Bile acids are amphiphilic molecules that have peculiar structure, because molecular descriptors that are obtained on the basis of the molecular graph or fragmentation methods often do not reflect their true structural features [4]. Therefore, the bile acid chromatographic lipophilicity play an important application in obtaining QSA(P)R models that connect biological and pharmacological or other physical-chemical properties (solubility, critical micelle concentration, critical micelle temperature, etc.) for their structure [5]. In the QSA(P)R models chromatographic parameters are independent variables.
In the following section presents the introduction chromatographic parameters (TLC and HPLC) that are used to represent the hydrophobicity of compounds, and presents the main structural features of bile acids.
2. Chromatographic parameters
The basic characteristic of the position of analytes spot on TLC chromatograms is the
If the solute spends more time in the mobile phase, then its chromatographic spot position is closer to the front of solvent and its
The ratio of time spent by the solute in the stationary phase
The connection between
The logarithm of the above expression is:
and the retention parameter log
In the above equation
RM0 parameters value depends on the type of organic modifier on the basis of which is determined by the function (1.5). RM0 is usually in the good correlation of the compounds lipophilicity. In recent times to describe the molecules lipophilicity also used the chromatographic parameters φ0 which is the ratio of RM0 and the slope S (Fig. 2) [7,8].
The chromatographic parameters φ0 corresponding volume ratio of organic modifier in aqueous mobile phase in which the same amount of solute in the mobile phase and in the stationary phase. Indeed solute in the above mobile phase composition, during the chromatographic process, spends at the same time in the stationary phase and the mobile phase, therefore capacitive factor (2) is k = 1, respectively Rf value (3) is 0.5. This Rf values corresponding to RM chromatographic parameter (4) whose value is zero. Which means that equation (5) is: 0=RM0 + Sφ0 from which follows the expression (6). If the solute is more hydrophobic then a larger quantity of organic modifier is needed to equalize the amount of solute in two phases, i.e. the chromatographic parameter φ0 has higher value. With a high efficient reverse phase liquid chromatography, stationary phase is also hydrophobic environment. Solute is characterized by retention time ts, which represents the elapsed time from injection to the occurrence of the same solute in detectors, i.e. the retention time of solute in the column. Chemical compounds in RP-HPLC analysis is usually characterized by the retention coefficient (capacity factor):
where t0 is the retention time of solvent from the mobile phase (Fig. 1 (a)). If a chemical compound is more hydrophobic, then more time is spent in the hydrophobic stationary phase, i.e. it takes more time to pass the column and the retention time or retention factor even greater. The logarithm of the retention coefficient logk is used as a parameter of hydrophobicity of chemical compound. An important chromatographic parameter of lipophilicity is logkw which were obtained an extrapolation of the linear equations (8) to the zero volume fraction φ of organic modifier in aqueous mobile phase.
Similar as in the RP TLC analysis, in RP HPLC also can defined the chromatographic indices φ0 (with the same meaning) as the ratio of logkw and the slope from logk = f(φ) [9,10].
There is a possibility of application principal component analysis (PCA) on chromatographic data to the thin layer and the liquid chromatography. PCA is applied to the data matrix D of retention parameters RM, logk (or k). The columns of the matrix D corresponding retention parameters for different volume fraction φ of organic modifier (columns represent the organic modifier) while the matrix rows represent the different types of chemical compounds. PCA is usually applied directly on the covariance matrix DTD retention parameters. As a result of PCAs mathematical procedure (orthogonal diagonalization) are obtained orthonorms PC score vectors of whose number is 2 (PC1, PC2) or 3 (PC1, PC2, PC3), depending on the percentage of variance explained from the data matrix D. Accordingly objects (tested molecules) from a multidimensional space matrix of retention parameters D mapped to 2d or 3d space of PC. In the graphs of PC scores can be found congeneric group of chemical compounds by their lipophilicity. As a parameter of lipophilicity of chemical compounds in the QSA(P)R studies can be applied also to the scores of PC [3,6].
3. The hydrophobic effect
Generally, the dissolution of an amphiphilic object in water is accompanied by the disruption of the hydrogen bonds between water molecules and formation of a hydration sheath (hydration layer) around the particles of the dissolved substance. If we observe such a system (amphiphilic solution), which consists of two subsystems: hydrophobic part (a) and hydrophilic part (b), then the thermodynamic functions can be considered separately for each subsystem (Fig. 4). For both subsystems it holds that in the formation of the solvation sheaths (around the amphiphilic objects) approximately the same number of hydrogen bonds are formed as existed between the water molecules in the bulk water, i.e. before introducing the amphiphilic object. Hence the change of the enthalpy for each of the subsystems is equal to zero: ΔH(a)≈ΔH(b)≈0. Also, in the formation of the solvation sheath, the entropy (translational and rotational) of both subsystems decreases, that is: ΔS(a)≈ΔS(b)<0. However, the water molecules from the hydrophilic side (b) of the amphiphilic molecule form additional hydrogen bonds; hence for this subsystem there is an additional negative enthalpy (ΔHɛ<0), which is dissipated as heath in the environment (bulk of the solution), thus giving rise to a positive change of the entropy (of the environment). Therefore, the free enthalpy change for the hydrophilic subsystem is: ΔG(b)<0, on the basis of which the water molecules from the hydrophilic side (b) of the amphiphil can be denoted as stabilized water molecules (SWM), while the water molecules from the hydrophobic side (a) are nonstabilized water molecules (NSWM) [4,11].
The ratio of the hydrophilic-to-hydrophobic surface area of an amphiphilic molecule determines the overall change of the Gibbs energy of formation of the hydration sheath (ΔG), and, since the hydrophobic surface of the amphiphilic molecule is larger, then ΔG>0. The larger the amount of the amphiphil present in the solution, the more water molecules participate in the formation of the hydration sheath, and the more negative is the overall entropy change. This results in the changes in the system (solution) due to the passing of NSMW from the amphiphil hydration sheath to the bulk of the solution, giving rise to the system entropy. Amphiphilic molecules (ions) are oriented so that their (formally) desolvated hydrophobic sides are to the smallest extent exposed to water molecules. This is realized either in the form of self-association of the amphiphils via their hydrophobic surfaces (Fig. 4) (entropy decrease due to the amphiphil self-association is smaller than the increase in the entropy due to their release from the hydration layer) or by binding of the amphiphil to the hydrophobic surface of the system: hydrophobic stationary phase, proteins, distribution in the organic solvent, etc. [11]. Generally, the hydrophobicity of a particle (molecule, ion) can be expressed via the logarithm of the coefficient of its partition (logP) between 1-octanol and water (P stands for the ratio of the equilibrium concentration of the particle of the same electronic structure in 1-octanol and its equilibrium concentration in water). If logP > 1, the particle is hydrophobic (lipophilic), and if logP < 1, the particle is hydrophilic [12]. Apart from the partition coefficient, the hydrophobicity of a particle can be also expressed via chromatographic parameters, either of the normal-phase or reversed-phase chromatography.
4. Structure of bile acids
Bile acid enzymatic produced in the liver of human and other mammals are the primary bile acids, and their intestinal microbial transformation to obtain the secondary bile acids. Both groups of bile acids are hydroxy derivatives of 5β-cholanic acid (cholic acid and chenodeoxycholic acid) [13]. The geometry of the 5β-cholanic acid largely determines of bile acids properties. The steroid skeleton of 5β-cholanic acid can be distinguished convex β and concave α surface (Fig. 5) [5].
Hydroxyl groups of bile acids are mostly oriented towards α side of steroid skeleton. X ray diffraction is demonstrated that the carboxyl group of the side chain from C17 carbon is also located on the α side of the steroid ring system. A concave (α) of bile acids steroid skeleton is the polar surface – hydrophilic side, while the convex surface (β) is the nonpolar surface – hydrophobic side. Simultaneous presence of hydrophobic and hydrophilic regions in the molecule of bile acids is named as amphiphilic. Since the carboxyl group ionized under physiological conditions, bile acids in biochemical systems belong to the ionic amphiphile. Also, bile acids are a special group of amphiphilic compounds, the group of biplanar molecules [4.5]. This is most apparent in cholic acid molecules (Fig. 6) because the molecules oxygen atoms from the α OH groups are in the same plane, the polar plane. In chenodeoxycholic and deoxycholic acid cannot be talking about the true planar polarity, but the literature often uses the term molecule with a hydrophilic edge.
If it oxidation of C7 or C12 α axial hydroxy group in cholic acid (C) (or deoxycolic acid (DC) or chenodeocycholic acid (CDC)) molecule leads to oxo derivative whose oxigen atom is switched for 60° (Newman projection) in relation to its' starting axial orientation i.e forms an angle of 30° with steroid skeleton mean plain (SSMP) (Newman projection, Fig. 7) [11,14]. Oxigen atom from C7 or C12 oxo group has the same steric orientation as has the equatorial C6 OH group from hyodeoxycholic acid (HDA).
Convergence of oxigen core (from oxo group) toward steroid sceleton β side means that there is the rise in number of stabilised water molecules (SWM, stabilised by hydrogen bonds) in hydratation layer of bile acids' steroid sceleton while hydrophobic surface of the convex side of the steroid sceleton decreases (Fig. 7). Thus, the number of unstabilised water molecules (NSWM) declines [11].
If bile acids' hydrophobic surface deceases (B), i.e. the number of NSWM is lower, the observed molecule bonds to the hydrophobic stationary phase with lower intensity, i.e. reaction (R1) is moved to the right. This change in hydrophobic surface reflects in the values of retention coefficient k (RPHPLC) and retention parameter RM0 (RPTLC).
5. Discussion: bile acids and the hromatographic parametars (indices)
In the RPHPLC(TLC), there is an equilibrium between the bile acid molecules (ions) from the polar mobile phase BA(aq) and the bile acid BA(s) adsorbed on the hydrophobic stationary phase.
This equilibrium is characterized by the equilibrium constant Kad = [BA(s)] / [BA(aq)], which is connected with the retention factor (capacity factor) (7) via the following relation [15]:
where Vhc represents the volume of hydrocarbon in the stationary phase, whereas Vpf is the volume of the polar phase. Their ratio is denoted as φ, and it represents a characteristic of the column. The change of the standard Gibbs energy
In both RPHPLC and RPTLC, in the adsorption of bile acids on the hydrophobic stationary phase, the water molecules from the solvation sheath that are not stabilized by hydrogen bonds (NSWM) pass to the hydrophilic mobile phase, giving a positive entropy contribution. Beside that, during the return of these water molecules in the bulk solution, NSWMs surface energy is released, which is then dissipated as a heat in the environment (solvophobic theory: according to this theory change in surface energy is a key factor in the reversed phase chromatographic processes energy). The formally dehydrated hydrophobic surfaces of bile acid molecules (ions) bind to the surface of the adsorbent by hydrophobic interactions (London dispersion forces, and dipol-induced dipole interactions). This makes a negative enthalpic contribution to the adsorption process, which is also dissipated as heat into the mobile phase (environment), giving rise to a positive change of the entropy of the environment.
Hence, the overall entropy change is positive, i.e. the Gibbs energy is lowered. The decrease in the free enthalpy of adsorption (ΔadGº) of bile acids (molecules or ions) on the hydrophobic stationary phase is more pronounced if the solvation layer of the steroid skeleton contains a greater number of NSWM, since then, there exists a largest increase in the overall entropy. In the binding of bile acids (molecules or ions) to the hydrophobic stationary phase, their steroid skeletons orient in such a way that the formally desolvated surfaces are to the smallest possible extent exposed to the solvent molecules from the hydrophilic mobile phase (methanol and water). Since bile acids are biplanar compounds, their molecules bind to the hydrophobic stationary phase by the side of the steroid skeleton that is most hydrophobic. With natural bile acids and their oxo derivatives, this is usually the convex β side of the steroid skeleton (an exception is the enantiomer of cholic acid with β configurations of all three OH groups, where α side is more hydrophobic than the β side of the steroid skeleton). Thus, the lnk chromatographic parameter has the liner dependence of ΔadGº, which means that lnk has also linear depending on the hydrophobic surface of the β side of steroid nucleus. Similar is true for RM chromatographic parameter (RPTLC). Therefore, lnk could serve as a molecular descriptor for the hydrophobicity of bile acids. But the studies dealing with the hydrophobicity of bile acids as expressed via the retention coefficients (k), there is a problem of mutual comparability of the lnk values. Because of that it is usually accustomed to give an order of hydrophobicity of the molecules, while the lnk values, though determined under identical experimental conditions, may vary because of the characteristic of the column (RPHPLC), which is in equation (9) denoted with φ. Because of that, Heuman introduced for bile acids a relative retention coefficient (kr) which is defined as the ratio of the retention coefficient of the given bile acid (k(BA)) and the retention coefficient of taurocholic acid (TC) (k(TC)) (Fig. 8) [16].
Taurocholic acid (37) was chosen as a reference because its retention coefficient in the pH range of the mobile phase from 2.9 to 9.0 is practically constant (complete ionization) [16]. By combining equations (10) and (11), the free enthalpy of adsorpion of bile acid (BA) on the hydrophobic stationary phase is:
In this way are eliminated the individual characteristics of the column, so that the relative retention coefficient depends only on the composition of the mobile phase and the degree of hydrophobicity of the stationary phase.
Besides, Heuman defined also the hydrophobicity index HI as the ratio of the logarithm of the relative retention coefficient of the given bile acid (BA) and of the relative retention coefficient of taurolithocholic acid (TL) (lnkr(TL)) (Fig. 7):
By introducing equation (13) into (12), one obtains the following equation:
from which is evident that there is a linear relation between the hydrophobicity index of a bile acid and the Gibbs energy of adsorption on the hydrophobic stationary phase.
The HI value for taurocholic acid (TC) is zero since its relative retention coefficient equation (11) is equal to one, and the same value (one) holds also for taurolithocholic acid (TL). If a certain bile acid is more hydrophilic than taurocholic acid (TC), its retention constant (k) is smaller than the k value of taurocholic acid (TC), so that the relative retention coefficient of that bile acid is smaller than unity (lnk<0), which means that, according to equation (13), the HI(BA) is negative. In the opposite case, if a certain bile acid is more hydrophobic than taurocholic acid (TC), then the value of its retention coefficient is greater than the retention coefficient of taurochollic acd (37), and the relative retention coefficient for that bile acid is larger than unity (lnk>0), that is HI>0 (Fig. 9).
Besides of Heumans hydrophobic index (based on data retention) logk is used to obtain predictive equations where logk is the independent variable while bile acid properties associated with its hydrophobicity is the dependent variable. Armstrong and Carey were given a linear regression equation, where logk is associated with the solubilized molar fraction of cholesterol ϑ (solubilized by micellar solutions of bile acid salts in the equilibrium conditions) [17]:
where
Natalini et al. as an experimental molecular descriptor used chromatographic parameter φ0 (6) for modeling the negative logarithm of the critical micellar concentration (pCMC = -logCMC) [9,10]. The study was included 20 non-conjugated and conjugated bile acids. For the full set of bile acids were obtained following linear regression equation:
R = 0.8246
Regression coefficient is improved if the equations of linear regression obtained separately for conjugated (18) and for non-conjugated bile acids (19):
R = 0.9643
R = 0.9380
Equations (18) and (19) have a good predictiveness in their congeneric groups, also with φ0 is well described hydrophobicity of the bile acids steroid skeleton β side. The mutual association of bile acid salts, especially around the first critical micellar concentration, which is through to their hydrophobic β side of steroid skeleton (Fig. 11) indicated that the chromatographic parameter φ0 adequately be used to predict its critical micellar concentration. Advantages of chromatographic parameters φ0 (6) with respect to the parameters RM0 (5) and logkw (8) is reflected in the absence of appropriate extrapolation of retention coefficients. Extrapolation error can be large if the linear area of equations (5) and (8) are away from the mobile phase without an organic modifier.
Sarbu et al. in RPTLC experiments examined which chromatographic parameter the best describes the hydrophobicity of bile acids. In the studied bile acids (Table 1) between RM0 parameter and the slope S (Equation 5) there is a good correlation, suggesting the presence of congeneric groups among chemical compounds in the tested assembly [3].
Lithocholic acid | 3α | |
Deoxycholic acid | 3α, 12α | |
Chenodeoxycholic acid | 3α, 7α | |
Ursodeoxycholic acid | 3α, 7β | |
Hyodeoxycholic acid | 3α, 6α | |
Hyocholic acid | 3α, 6α, 7α | |
Cholic acid | 3α, 7α, 12α | |
Glycolithocholic acid | 3α | |
Taurolithocholic acid | 3α | |
Glycochenodeoxycholic acid | 3α, 7α | |
Taurochenodeoxycholic acid | 3α, 7α | |
Glycodeoxycholic acid | 3α, 12α | |
Taurodeoxycholic acid | 3α, 12α | |
Glycocholic acid | 3α, 7α, 12α | |
Taurocholic acid | 3α, 7α, 12α |
Application of PCA on the matrix of RM values (D) of the studied bile acids are obtained space of principal components with smaller dimensions then matrix D. Screen plot showed that the significant first and second principal components (PC1 and PC2) that explained 99.75% of variance in the starting matrix of retention parameters. Bartlett’s statistic shows that PC1, PC2 and PC3 are significant which explained of 99.92% of the total variance. In the set of studied bile acids (Table 1) there is a good correlation between RM0 chromatographic parameters and scores of PC1:
R = - 0.9203
5.33 | 0.195 | 0.122 | 0.080 | ||||
4.78 | 0.346 | 0.179 | 0.111 | 3.50 | 2.65 | ||
5.01 | 0.340 | 0.177 | 0.082 | 3.28 | 2.25 | ||
3.65 | 0.545 | 0.226 | 0.114 | 3.00 | 2.20 | ||
3.88 | 0.488 | 0.208 | 0.110 | 3.08 | 2.28 | ||
3.66 | 0.559 | 0.229 | 0.107 | 2.80 | 1.84 | ||
3.84 | 0.533 | 0.217 | 0.105 | 2.02 | 1.10 | ||
3.39 | 0.674 | 0.246 | 0.119 | ||||
2.12 | 1.000 | 0.249 | 0.077 | ||||
2.87 | 0.819 | 0.237 | 0.088 | 2.12 | 0.45 | ||
2.27 | 1.249 | 0.194 | 0.082 | ||||
2.90 | 0.850 | 0.247 | 0.087 | 2.25 | 0.80 | ||
2.30 | 1.165 | 0.191 | 0.090 | ||||
2.30 | 1.102 | 0.186 | 0.090 | 1.65 | -0.40 | ||
2.29 | 1.524 | 0.137 | 0.125 |
Sarbu for assessing whether RM0 or score of PC adequate to describe lipophilicity of bile acids using the experimental data of particion coefficients between 1-octanol and water of bile acids in molecular (HA) and ionized (A) forms. Between the logarithm partition coefficients and chromatographic parameters on the basis of Table 1 are obtained by the following equation of linear regression:
R = 0.8559
R = 0.8751
R = 0.8675
R = 0.9379
R = 0.8715
R = 0.9395
Linear regression equations (23 and 24) in which the chromatographic parameter RM0 replaced with a score of PC1 has correlation coefficients with larger values compared to the equations (21 and 23). Correlation coefficient is slightly increased in the multiple linear regression equation for logP (molecular and ionized form) when the equations involved scores of PC1 and PC2. Increasing the coefficient of correlation in the above equations can be explained with Wolds assumption according to which measurement errors are contained in the principal components with small eigenvalues, so the scores of PC1 is purified from the errors. In Fig. 12 bile acids are shown from the Table 2 in the plane of the scores of principal components. It can be seen to form three linear congeneric groups: non- conjugated bile acids (BA), glyco-conjugate (G-BA) and tauro-conjugated (T-BA). Grouping of bile acids on the basis of the scores of PC as the appropriate regression equation for logP (from 23 to 26), indicating that the PC1 scores is adequate replacement for RM0 chromatographic parameters to describe the lipophilicity (hydrophobicity) of bile acids, namely the scores of PC1 can be used as molecular descriptors for hydrophobicity of the steroid skeleton β (bile acids).
If increase the set of bile acids from Table 1 by including oxo and acetoxy derivatives, RM0 chromatographic parameter also included in the multiple linear regressions for the logPHA and the logPA with molecular descriptors such as polarizability and molar specific sum of absolute values of the charges on each atom of the molecule (electronic parameter), molecular volume, the third order of the connectivity index and Wiener index (shape parameter). Based on the loadings of PC can be concluded that RM0 carries different information about the structure of the studied bile acids than the above in silico molecular descriptors [6].
Effect of temperature on the chromatographic process in reverse phase chromatography can also be applied to obtain chromatographic parameters that have properties of molecular descriptors [15]. Namely, the reaction (R2) is exothermic (Had <0), then based on the van' t Hoff equation, the equilibrium constant of reaction (R2) decreases
This means that the retention factor k (9) also decreases. Thus, the increase in temperature has the same effect on the value of k as there are increasing the volume fraction of organic modifier in aqueous mobile phase (reversed phase chromatography). In the experiment, the influence of temperature with an the retention factor at constant mobile phase included 25 different bile acids (conjugated and non-conjugated) (Table 3) [18].
For each tested bile acids between retention factor k and temperature t there is a linear relation:
The linear model (28) explains from 96% to 99% of the variance at tested bile acids (Table 4). In the set of the studied bile acids (Table 3) between the parameters of equation (28) (a and b) there is a good correlation, indicating the formation of congeneric groups in the plane of a and b:
R=0.9996
Lithocholic acid | 3α-OH | |
Deoxycholic acid | 3α-OH, 12α-OH | |
Chenodeoxycholic acid | 3α-OH, 7α-OH | |
Cholic acid | 3α-OH, 7α-OH, 12α-OH | |
Ursodeoxycholic acid | 3α-OH, 7β-OH | |
Hyodeoxycholic acid | 3α-OH, 6α-OH | |
Hyocholic acid | 3α-OH, 6α-OH, 7α-OH | |
Glycochenodeoxycholic acid | 3α-OH, 7α-OH | |
Taurodeoxycholic acid | 3α-OH, 12α-OH | |
Glycocholic acid | 3α-OH, 7α-OH, 12α-OH | |
Glycodeoxycholic acid | 3α-OH, 12α-OH | |
Taurolithocholic acid | 3α-OH | |
Taurochenodeoxycholic acid | 3α-OH, 7α-OH | |
Glycolithocholic acid | 3α-OH | |
12-Monoketocholic acid | 3α-OH, 7α-OH, 12α-Oxo | |
7-Monoketocholic acid | 3α-OH, 7α-Oxo, 12α-OH | |
7,12-Diketocholic acid | 3α-OH, 7α-Oxo, 12α-Oxo | |
3,7-Diketocholic acid | 3α-Oxo, 7α-Oxo, 12α-OH | |
3,12-Diketocholic acid | 3α-Oxo, 7α-OH, 12α-Oxo | |
Dehydrocholic acid | 3α-Oxo, 7α-Oxo, 12α-Oxo | |
12-Monoketodeoxycholic acid | 3α-OH, 12α-Oxo | |
Diketodeoxycholic acid | 3α-Oxo, 12α-Oxo | |
7-Monoketochenodeoxycholic acid | 3α-OH, 7α-Oxo | |
Diketochenodeoxycholic acid | 3α-Oxo, 7α-Oxo | |
6-Monoketohyodeoxycholic acid | 3α-OH, 6α-Oxo |
Fig. 13 shows that the tested bile acids in the plane of parameters of equation (23) formed by three congeneric groups. The first group (I) formed by bile acids with oxo groups in the steroid nucleus and the bile acids that in addition to C3 have one more equatorial OH groups. The α equatorial OH groups have the same position relative to the steroid skeleton mean plane such have the α oxo groups. However, in this group also includes the cholic acid and the glycocholic acid which is structurally not fit into the above mentioned group. The second group (II) formed deoxycholic acid and chenodeoxycholic acid with its conjugates. The third group (III) forms lithocholic acid and its conjugates. Question may arise as whether can be improved with grouping of bile acid gained by in the plane of a and b, especially in terms of group I which provides the vendor cholic acid? The answer may be the application of methods of principal components of the matrix of the data of retention factor dependence on temperature. The first principal component (PC1) explains 99.99715% of variance in the starting data matrix, while the remaining variance explained by PC2. In Table 4 presents the scores of principal components of bile acids studied.
43.0501 | 0.7200 | 0.9879 | 5.37637 | 0.023150 | ||
21.8893 | 0.3518 | 0.9906 | 1.79015 | -0.051081 | ||
18.4260 | 0.3021 | 0.9807 | 1.05575 | 0.020646 | ||
8.4873 | 0.1319 | 0.9828 | -0.66245 | 0.016292 | ||
4.4884 | 0.0680 | 0.9843 | -1.40668 | 0.031446 | ||
5.6203 | 0.0804 | 0.9538 | -1.11512 | -0.034642 | ||
5.5477 | 0.0864 | 0.9836 | -1.23205 | 0.041543 | ||
16.6862 | 0.2651 | 0.9946 | 0.85970 | -0.024730 | ||
19.5503 | 0.3165 | 0.9802 | 1.30724 | -0.001736 | ||
7.8569 | 0.1213 | 0.9841 | -0.77443 | 0.014749 | ||
20.6286 | 0.3336 | 0.9786 | 1.51098 | -0.012582 | ||
39.1628 | 0.6544 | 0.9728 | 4.65836 | 0.029036 | ||
16.3808 | 0.2671 | 0.9796 | 0.70668 | 0.007651 | ||
41.1016 | 0.6926 | 0.9750 | 4.95053 | 0.026326 | ||
3.1340 | 0.0402 | 0.9704 | -1.58139 | -0.004996 | ||
2.8202 | 0.0353 | 0.9676 | -1.64238 | -0.002356 | ||
0.7141 | 0.0078 | 0.9769 | -2.12906 | 0.056420 | ||
0.6100 | 0.0054 | 0.9974 | -2.11715 | 0.056756 | ||
0.6237 | 0.0060 | 0.9752 | -2.12981 | 0.056122 | ||
0.5042 | 0.0056 | 0.9738 | -2.17086 | 0.066570 | ||
7.7659 | 0.1040 | 0.9652 | -0.57871 | -0.116428 | ||
4.1490 | 0.0511 | 0.9681 | -1.32315 | -0.048688 | ||
6.1746 | 0.0817 | 0.9665 | -0.91824 | -0.078586 | ||
3.4944 | 0.0423 | 0.9704 | -1.46618 | -0.031508 | ||
6.3614 | 0.0913 | 0.9667 | -0.96809 | -0.039373 |
In the plane of principal component scores are obtained by analogous groups such as Fig. 14. However in the group I the oxo derivatives and the derivatives with α equatorial OH groups forming linear congeneric group, while based on Cooks distance cholic acid and glycocholic acid may be considered as outliners. With loadings of principal component can be determined by mutual correlation of different in silico descriptors (topological and electronic) and obtained experimental descriptors a, b, and PC1 on the set of 25 bile acids (Table 4). It may be concluded that chromatographic parameters a and PC1 are more or less orthogonal to most topological descriptors (Winers index, Balabans index, shape attributes, konectivities, etc.), while PC1 and a are mutually parallel so that explain the same structural characteristics of bile acids. Namely, as flat molecules nitrazepan guest molecule is incorporated in the hydrophobic domain of micelles. Therefore, if the β side of bile acids steroid skeleton more hydrophobic then nitrazepans partition coefficient has a higher value. Using the method of partial least square (PLS) regression equation was obtained for the nitrazepans partition coefficients in which experimentally determined variables a and PC1 are the most significant with respect to in silico descriptors. Which means that PC1 and a can be used as molecular descriptors of bile acids lipophilicity, and that the temperature dependence of retention factors may obtain information regarding the structure of bile acids (lipophilicity).
6. Conclusions
Chromatographic data (parameters) can be used for prediction of following feaatures of bile acids:
Critical Micellar Concentration
Haematolityc potentitial
Binding for Albumin
Solubilization of cholesterol, lecitin, vitamin E etc.
Mixed micelles parameters (for example partition coefficient of nitrazepan in bile acids' micelles)
Application of retention parameters for prediction of bile acids' properties is convenient for finding the size of the starting set of molecules. If the set is too small (less than ten molecules), by using in silico descriptors multiple regression equations are obtained which are overfitted, while by using chromatographic parameters that posses appropriate structural informations about bile acids, simple linear equations are derived with acceptable predictive power.
Acknowledgement
This study was supported by the Secretariat for Science and Technological Development of AP Vojvodina, Republic of Serbia (Grant No. 114-451-2113/2012).
Abbreviations
QSA(P)R = Quantitative Structure Activity (Property) Relationship
RPHPLC = Reversed Phase High Pressure (Resolutions) Liquid Chromatography
RPTLC = Reversed Phase Thin Layer Chromarography
PCA = Principal Component Analysis
PC = Principal Component
SWM = Stabilised Water Molecules
NSWM = Non Stabilised Water Molecules
SSMP = steroid skeleton mean plane
References
- 1.
Camile W. 2003 The Practice of Medicinal Chemistry. Oxford: Academic Press. - 2.
King F.D 2002 Medicinal Chemistry Principles and Practice. Cambridge: R.S.C. - 3.
Sârbu C. Kuhajda K. Kevrešan S. 2001 Evaulation of the lipofilicity of bile acids and their derivatives by thin-layer chromatography and principal component analysis. J. Chromatogr. A917 361 366 - 4.
Poša M. 2011 QSPR study of the effect of steroidal hydroxy and oxo substituents on the critical micellar concentrations of bile acids. Steroids 76(1-2):85 EOF 93 EOF - 5.
Poša M. Kuhajda K. 2010 Hydrophobiciti and haemolytic potential of oxo derivatives of cholic, deoxycholic and chenodeoxycholic acids. Steroids75 6 424 431 - 6.
Sârbu C. Onişor C. Poša M. Kuhajda K. Kevrešan S. 2008 Modeling and prediction (correction) of partition coefficients of bile acids and their derivatives by multivariate regression methods. Talanta75 651 657 - 7.
Onişor C. Poša M. Sârbu C. 2011 Modeling of chromatographic lipophilicity indices of bile acids and their derivatives. Rev. Roum. Chim.56 3 217 229 - 8.
Onişor C. Poša M. Kevrešan S. Kuhajda K. Sârbu C. 2010 Modeling chromatographic lipophilicity of bile acids and their derivatives estimated by HP-RPTLC. J. Sep. Sci.33 20 3110 3118 - 9.
Natalini B. Sardella R. Camaioni E. Gioiello A. Pellicciari R. 2007 Correlation between CMC and chromatographic index: simple and effective evaluation of the hydrophobic/hydrophilic balance of bile acids. Anal. Bioanal. Chem.388 1681 1688 - 10.
Natalini B. Sardella R. Camaioni E. Macchiarulo A. Gioiello A. Carbone G. Pellicciari R. 2009 Derived chromatographic indices as effective tools to study the self-aggregation process of bile acids. J. Pharm. Biomed. Anal.50 613 621 - 11.
Poša M. Pilipović A. Lalić M. 2010 The influence of NaCl on hydrophobicity of selected, pharmacologically active bile acids expressed with chromatographic retention index and critical micellar concentration. Colloids Surf. B: Biointerfaces81 336 343 - 12.
Roda A. Minutello A. Angellotti M. A. Fini A. 1990 Bile acid structure-activity relationships: evaluation of bile acids lipophilicity using 1-octanol/water partition coefficient and reverse phase HPLC. J. Lipid Res.31 1433 1443 - 13.
Mikov M. Fawcett J. P. 2007 Bile Acids. Geneva: Medishet Publisher. - 14.
Poša M. Farkaš Z. 2010 Cholesterol solubilization by oxo derivatives of selected bile acids and their osmotic resistance. Collect. Czech. Chem. Commun.75 8 767 784 - 15.
Poša M. Pilipović A. Lalić M. Popović J. 2010 Hydrophobicity and retention coefficient of selected bile acid oxo derivatives. Acta Chim. Slov.57 4 828 835 - 16.
Heuman M.D 1989 Quantitative estimation of the hydrophilic-hydrophobic balance of mixed bile salt solutions. J. Lipid Res.30 719 730 - 17.
Armstrong J.M, Carey C.M 1982 The hydrophobic-hydrophilic balance of bile salts. Inverse correlation between reverse-phase high performance liquid chromarographic mobilities and micellar cholesterol-solubilizing capacities. J. Lipid Res.23 70 80 - 18.
Poša M. Pilipović A. Lalić M. Popović J. 2011 Determination and importance of temperature dependence of retention coefficient (RPHPLC) in QSAR model of nitrazepams’ partition coefficient in bile acid micelles. Talanta83 5 1634 1642