Journal of Liquid Chromatography & Related Technologiesw, 31: 79–96, 2008
Copyright # Taylor & Francis Group, LLC
ISSN 1082-6076 print/1520-572X online
DOI: 10.1080/10826070701665626
Current State of the Art in HPLC
Methodology for Lipophilicity Assessment
of Basic Drugs. A Review
Costas Giaginis and Anna Tsantili-Kakoulidou
1
Department of Pharmaceutical Chemistry, School of Pharmacy,
University of Athens, Athens, Greece
Abstract: HPLC provides a user’s friendly, rapid, and compound sparing methodology, which is successfully applied to determine drug lipophilicity. Under suitable
chromatographic conditions isocratic and extrapolated retention factors correlate
well with octanol-water partition or distribution coefficients. The present review
provides an overview of the stationary and mobile phases, which are preferably used
for lipophilicity assessment mainly in the case of basic compounds. Difficulties
raised by the interference of silanophilic interactions in the partition mechanism, and
the ways proposed to face this problem are discussed. Attention has been given to
the extrapolation procedure and the standardization of conditions to obtain 1:1 correlation between extrapolated retention factors and logP or logD. Other chromatographic indices encoding information on the lipophilic behavior are briefly
presented. A separate section refers to recent advances in IAM Chromatography, its
similarities/dissimilarities with reversed phase HPLC and the octanol-water system,
as well as its potential to mimic specific interactions with phospholipids.
Keywords: Lipophilicity indices, Basic drugs, n-Octanol-water system, Reversedphase HPLC, IAM chromatography, Extrapolated retention factors
INTRODUCTION
Lipophilicity, expressed by the logarithm of octanol-water partition coefficient
logP or distribution coefficient logD, if ionized molecular species are present,
Correspondence: Anna Tsantili-Kakoulidou, Department of Pharmaceutical
Chemistry, School of Pharmacy, University of Athens, Panepistimiopolis, Zografou
Athens 157 71, Greece. E-mail: tsantili@pharm.uoa.gr
79
80
C. Giaginis and A. Tsantili-Kakoulidou
constitutes a physicochemical property of paramount importance for the
medicinal chemist. It plays an important role in ADME (Absorption, Distribution, Metabolism, and Elimination) characteristics of drugs, while
affecting their pharmacodynamic and toxicological profile, as well.[1 – 3]
Although lipophilicity is essential for penetration across biological
membranes and hydrophobic interactions with receptors, high logP/logD
values are associated with undesired drug features, like extensive and unpredictable metabolism, high plasma protein binding, or accumulation to
tissues.[4]
Basic compounds represent the major fraction in drug related databases
rendering their lipophilicity assessment as an urgent requirement in drug
design. Computed lipophilicity values are often inaccurate, especially if
they refer to ionized or partially ionized molecules and, although valuable
for screening virtual libraries, they should be replaced by measured data
early in the drug discovery process.[5,6] Different experimental protocols for
logP or logD determination have been reported in literature. The classical
shaking flask method for direct partitioning experiments is tedious and time
consuming, not suitable for degradable compounds, less amenable to automation, while it presents limitations concerning the logP/logD range which
can be reliably measured.[3,7,8] The dual phase potentiometric titration
suitable for ionisable drugs, on the other hand, requires special equipment,
not always available in an analytical laboratory.[9,10]
Reversed-phase high performance liquid chromatography (RP-HPLC)
has proven to simulate octanol-water partitioning and is considered as a
popular alternative for lipophilicity assessment. It offers several practical
advantages, including speed, reproducibility, insensitivity to impurities or
degradation products, broader dynamic range, on line detection, and reduced
sample handling and sample sizes.[11 – 13] These advantages have attracted considerable interest and the literature is rich in research articles, which investigate
the relationship of chromatographic retention with octanol-water partitioning
and the common factors underlying the two processes.[11 – 16]
The criticism towards octanol as an isotropic medium with only a superficial similarity to biomembranes and the difficulties associated with the use
of liposomes as more representative models,[17] have triggered the development of immobilized artificial membrane (IAM) stationary phases for use in
HPLC. IAM chromatography has unfolded new perspectives in the application of HPLC as a tool to mimic specific interactions with
phospholipids.[18,19]
The major part of the present review considers the chromatographic conditions that are more suitable for lipophilicity assessment focusing on basic
drugs, evaluates associated difficulties, and provides an overview on the
relation between RP-HPLC and the reference octanol-water partitioning
system. In a separate section, analogous aspects in respect to IAM chromatography and its potential as a tool for rapid evaluation of drug permeation and/
or interactions with biological membranes are discussed.
Lipophilicity Assessment of Basic Drugs
81
LIPOPHILICITY INDICES IN REVERSED-PHASE HPLC
The lipophilicity index measured by HPLC is derived by the retention time tr
that is converted to the logarithm of the retention factor log k according to
Equation (1):
tr to
log k ¼ log
ð1Þ
to
where to being the retention time of an unretained solute.
Isocratic retention factors represent a relative scale of lipophilicity. They are
preferred by some authors since they require fewer experiments.[20] However,
extrapolated retention factors logkw, corresponding to pure water as mobile
phase, are considered as more representative lipophilicity indices, their values
being of the same order of magnitude as octanol-water logP/logD.[11–16] Extrapolated logkw values are derived using the linear part of the logk/f relationships,
where f is the concentration of the organic modifier in the mobile phase. This
issue is discussed extensively in the Extrapolation Procedure section.
Both isocratic log k and extrapolated logkw values are directly correlated
to octanol-water logP/logD via Collander-type equations (Equation (2):
logP=logD ¼ a logkðwÞ þ b
ð2Þ
where a, b, constants derived by linear regression analysis.
Equations of type 2 are constructed using compounds with known logP/
logD values and can serve as calibration equations for further logP/logD calculations.[21] In many cases, when logkw values are used, a and b in Equation
(2) tend to approach 1 and 0, respectively. In such cases, retention and octanol
water partitioning are considered as homoenergetic processes. The quality of
type 2 equations, however, depends both on the chromatographic conditions
and the nature of the solutes.[20,21] Solvatochromic analysis has revealed
differences in the balance of factors involved in octanol-water partitioning
and reversed phase retention,[22,23] while conformational effects have also
been manifested.[24] In this aspect, care should be taken in the selection of
the training set of solutes for logP/logD estimation by HPLC. Nevertheless,
considerable research efforts are directed towards the standardization of chromatographic conditions, which attenuate dissimilarities between retention and
octanol-water partitioning and guarantee 1:1 correlation between logkw and
logP/logD values for structurally diverse compounds.[25 – 27]
CHROMATOGRAPHIC CONDITIONS FOR LIPOPHILICITY
ASSESSMENT
Stationary Phases
C18 silanized silica gel is the most preferred packing material for reversed
phase columns in the chromatographic analysis of basic drugs. The same
82
C. Giaginis and A. Tsantili-Kakoulidou
material is appropriate for drug lipophilicity assessment, as well. However, the
interference of silanophilic interactions in the partitioning mechanism of RPHPLC has been recognized as a serious drawback, especially in the case of
basic drugs.[28,29] Silanophilic interactions are attributed to the remaining
free silanol sites and include hydrogen bonding as well as electrostatic
forces, especially in the case of positively charged basic compounds,
producing considerable increases in retention.[28 – 32] They also depend on
the degree of ionization of the silanol groups, being less pronounced at low
pH.[33,34]
The problem of silanophilic interactions is partially faced by the development of columns with reduced free or accessible silanol sites. End-capping of
the silanol residues by trimethylchlorosilane (TMCS) or hexamethyldisilazane
(HMDS) is usually performed during the manufacturing process, leading to a
higher degree of silanization.[31,32] Hence, base deactivated silica represents a
packing material more suitable (e.g., BDS C18) for basic solutes. In addition,
recent technology has led to the development of polar embedded and polar
endcapped stationary phases, which are considered to be further deprived
from silanophilic effects.[26,35 – 37] With respect to embedded columns, a
polar functional group, such as amide, carbamate, ether, or sulfonamide, is
incorporated at the bottom of the alkyl bonded chains. This functional
group provides electrostatical shielding to the surface silanol sites. The LCABZþ and the Discovery-RP-Amide-C16 stationary phases belong to these
types of columns. However, these packing materials, which have been used
for lipophilicity determination, may exhibit other polar interactions with
analytes, such as the strong interaction between the polar embedded groups
and the phenolic analytes.[38] On the other hand, the need of a masking
agent in the case of basic compounds (see below) with these types of
columns suggests that silanophilic interactions still persist.[26] With respect
to polar endcapped columns, a second reaction is used to bond a short
carbon chain (usually C3-C4) with a polar end to the surface silanol sites. A
favourable advantage for both types is the fact that a higher degree of orientation for the alkyl chains is achieved and, thus, they can be used with mobile
phases containing high amounts of water or even pure water without the
problem of hydrophobic collapse.[39,40]
The pH limitation of the above mentioned columns lies in the range 2.5 to
7.5. Thus, in the case of basic compounds, they do not allow determination of
retention factors corresponding to the neutral form. In this case, logkw of the
neutral form can be estimated from the apparent logkapp
w by using Equation (3)
adapted from the analogous logP/logD relationship.
pKapH
logkw ¼ logkapp
Þ
w þ logð1 þ 10
ð3Þ
Equation (3) is assumed to be valid for isocratic logk values, as well.
Nevertheless, whether the effect of ionization in octanol-water and HPLC
system is similar remains to be clarified, in as much as secondary interactions
Lipophilicity Assessment of Basic Drugs
83
and the influence of organic modifier may lead to deteriorations in the logk/
pH profile (see further details in sections below).
Recently bidentate stationary phases (e.g., Zorbax-extend C218) that
include a propylene bridge, as well as surface modified silica columns (e.g.,
XTerra C218), where organic functional groups have become a constituent
of the silica backbone, allow the use of mobile phases with pH up to 12.[41]
However, the applicability of such columns in the lipophilicity assessment
of basic pharmaceuticals has not been systematically investigated yet. In a
recent publication, 1:1 correlation has been reported between logkw and
logP for 40 basic compounds measured at pH 10.5 on a Zorbax-extend
C218 column without addition of any masking agent.[42]
As an alternative choice, the polymer based octadecyl-poly(vinyl alcohol)
(ODP) stationary phase, which is completely devoid of reactive silanol groups,
has also been used for lipophilicity measurements.[15,43] The ODP column
presents stability to acidic and strongly basic conditions (at pH between 2
and 13).[44] However, it has been reported that the retention mechanism on
ODP stationary phase compared to octanol-water partitioning is controlled
by a different balance of forces. Thus, the derived data may not be so
suitable to reproduce the classical log P or logD values.[45]
Mobile Phases
The most extensively used mobile phases in RP-HPLC are mixtures of water
or buffer with an organic modifier, usually methanol, acetonitrile, or THF.
However, acetonitrile was found to produce the most asymmetrical peaks in
the analysis of organic bases.[46] This fact was attributed to the inability of
acetonitrile to form hydrogen bonds with residuals silanols, in contrast to
methanol and THF. In terms of lipophilicity assessment, methanol seems to
be the most suitable organic modifier for RP-HPLC, since it does not
disturb the hydrogen bonding network of water. Moreover, during equilibration, methanol molecules associate with the stationary phase forming a
monolayer, which provides a hydrogen bonding capability in better
agreement with n-octanol.[47]
It should be taken into account that organic modifiers are capable of
affecting the pKa of ionized solutes, as well as the acidity of the surface
silanol groups and the pH of the mobile phase. In general, the pKa of bases
decreases as the organic modifier concentration increases. However, substantial structure-dependent differences in pKa shifts for bases at a given organic
solvent composition as well as pH variations at different organic solvent compositions have been reported.[48 – 50] These effects are minimized in the extrapolation procedure, providing a further argument for the use of logkw values
instead of isocratic logk in the case of basic compounds.
The buffer composition of the aqueous component in the mobile phase
also plays an active role in the retention of protonated basic compounds,
84
C. Giaginis and A. Tsantili-Kakoulidou
which may form ion pairs with the counter ions. Morphilinepropanesulfonic
acid (MOPS), is considered as the buffer of choice for lipophilicity assessment by HPLC.[26,51,52] It exhibits a large buffering capacity coupled to
poor ion pair formation ability due to its zwitterionic nature and, thus, it
does not interfere either with solutes or with stationary phase. On the
other hand, the partitioning experiments for logD determination are
usually performed in phosphate buffer or in phosphate buffered saline
(PBS), containing NaCl and KCl at a total concentration of approximately
0.16 M. to mimic the isotonic physiological conditions.[53] Hence, this
choice is often used in HPLC as well. However, phosphate and, especially,
the chloride anions are capable of forming ion pairs with protonated
molecules with extraction constants that may differ from those in
octanol-water.
In the case of basic drugs, the addition of small amounts (0.15 –0.20% v/
v) of amines to the mobile phase is a critical prerequisite in order to suppress
silanophilic interactions, even if polar embedded or polar end capped stationary phases are used. Hydrophobic amines, such as n-decylamine and N,Ndimethyloctylamine, are considered to be the most suitable masking agents
combined with methanol as organic modifier.[15,26] The effect of hydrophobic
amines on retention is less evident with acetonitrile as organic modifier.
Acetonitrile, as a weak hydrogen bonding solvent, is not capable of
solvating the stationary phase with sufficient water, thus presumably preventing the positively charged amine to be dragged on the column and to exert its
role as a masking agent.[15]
Recently, room temperature ionic liquids of the imidazolium tetrafluoroborate family, such as 1-butyl-3-methylimidazolium (BMIM BF4), have been
reported to be suitable for suppressing silanophilic interactions for a set of
b-blockers.[54 – 56] In particular, this type of mobile phase additive combines the
silanol masking effect of the imidazolium cation with the chaotropic character
of the BF2
4 anion, providing a promising masking agent which may release
new perspectives for the optimization of lipophilicity determinations.[54 – 56]
EXTRAPOLATION PROCEDURE
If the entire organic modifier range is considered, the relationship between
retention factors and the fraction of the organic modifier f follows the
Schoenmaker’s solubility parameter model according to Equation (4).[57]
pffiffiffiffi
logk ¼ Af þ Bf2 þ E f þ logkw
ð4Þ
A, B, and E are fitting coefficients and logkw is the intercept corresponding to 100% aqueous phase. The Bf2 term accounts for the curvature
(concave) at higher organic p
modifier
concentrations partly attributed to silanoffiffiffiffi
philic interactions, while E f accounts for a curvature (concave or convex)
Lipophilicity Assessment of Basic Drugs
85
observed at water rich mobile phases (f , 0.2) due to stationary phase
solvation problems. The error in the logkw values produced as a result of
the curvature at lower fractions of organic modifier has recently been investigated by Tate et al.[58] The correct estimation of extrapolated chromatographic
indices also depends on the stationary phase. Indeed, a hydrid based polar
embedded column produced a slightly smaller error in extrapolation
procedure than a polar endcapped and a conventional non-polar endcapped
column, as a result of both lower surface area and less surface silanols.[58]
Nevertheless, N-N lone pair interactions between amide embedded groups
and the solutes containing nitrogen atoms, like basic compounds, seem to
affect the retention characteristics and thus the extrapolation accuracy.[58]
Quadratic extrapolation using the higher organic modifier concentration
range may also lead to erroneous values in respect to lipophilicity.[15,27]
Nevertheless, when methanol is used as organic modifier at fractions
.0.2 coupled to a masking agent, the linear part of Equation (4) is sufficiently
wide. Hence, it can be used to derive extrapolated logkw values according to
Snyde’s linear solvent strength model[59] via Equation (5):
logk ¼ Sf þ logkw
ð5Þ
Linear extrapolation is generally preferred to obtain logkw values representative of lipophilicity. It is assumed that linearity holds better for
modifier concentrations that produce 0 , logk , 1.[60,61]
STANDARDIZATION OF THE CHROMATOGRAPHIC
CONDITIONS FOR THE LIPOPHILICITY ASSESSMENT
OF BASIC DRUGS
Attempts to optimize the chromatographic conditions in the aim to simulate
better octanol-water partitioning have already been reported 30 years ago
by Unger et al.[62,63] These authors suggested a reversed-phase C18 packing
material as stationary phase previously coated with n-octanol and the use of
pure n-octanol saturated buffer as mobile phase. A good correlation
between logP values and logkw was obtained; however, the basic
compounds were strongly retained in the column interacting with the
surface silanol sites and disrupting the correlation.[62] A hydrophobic amine,
N,N-dimethyloctylamine, was further added to suppress the silanophilic interactions, leading to a very good correlation.[63] However, a rather limited set of
basic drugs, including phenothiazines and tricyclic antidepressants, was used.
Moreover, retention factors were less reproducible due to column instability
and bleeding. Recently, the use of n-octanol was revisited not as a principal
solvent component of the stationary phase, but as a mobile phase additive.
In fact, the addition of 0.25% n-octanol in the methanol fraction of mobile
phase coupled to n-octanol saturated MOPS buffer produced a very good
86
C. Giaginis and A. Tsantili-Kakoulidou
correlation for a C8 column.[52] Triethylamine or n-decylamine was used as
masking agent, with the latter to be advantageous concerning logP/logkw
relationships. However, the data set was limited, including nonfunctional
solutes only, while the dynamic range of logD values did not exceed the
three log units.
Based on this evidence, Lombardo et al. proposed a LC-ABZþ column as
stationary phase and mobile phase conditions similar to Minick (methanol as
organic modifier þ0.25% n-octanol, MOPS as n-octanol saturated buffer and
0.15% n-decylamine in respect to the total volume).[26] For a set of 163 structurally diverse basic and neutral drugs, a calibration curve has been obtained
reflected in Equation (6), which covers a dynamic range broader than seven
log units. To achieve the potential for automation, three isocratic logk
values were used for the extrapolation to logkw according to three lipophilicity
ranges.
logD7:4 ¼ 1:08 ð+0:02Þlogkw þ 0:20 ð+0:04Þ
n ¼ 163;
r2 ¼ 0:949;
s ¼ 0:369;
F ¼ 3000
ð6Þ
Equation (6) represents a general practically 1:1 correlation with a slope
close to unity and an intercept close to zero and is used to calculate logD
values at pH 7.4. The method, introduced as ElogD7.4, has been validated
for a large number of neutral and basic drugs.[64]
To this point, a conventional base deactivated silica column (BDS C218)
has successfully been applied for the lipophilicity assessment of 64 structurally diverse basic and neutral drugs using analogous mobile phase conditions.[27] Correlation of 1:1 with high statistics was obtained according to
Equation (7), which covers a dynamic range of six log units.
logD7:4 ¼ 1:03 ð+0:03Þ logkw þ 0:14 ð+0:07Þ
n ¼ 64;
r2 ¼ 0:937;
s ¼ 0:288;
F ¼ 908
ð7Þ
When basic drugs were analyzed separately from neutral ones, an
analogous equation was obtained (Equation 8).
logD7:4 ¼ 1:07 ð+0:04Þ logkw þ 0:00 ð+0:09Þ
n ¼ 40;
r2 ¼ 0:943;
s ¼ 0:278;
F ¼ 632
ð8Þ
If only n-decylamine is added in the mobile phase Equation (9) even
better statistics was obtained accompanied, however, by a high intercept,
which could be attributed to the presence of silanophilic interactions despite
the use of the masking agent.
logD7:4 ¼ 1:08 ð+0:03Þ logkw 0:64 ð+0:08Þ
n ¼ 40;
r2 ¼ 0:970;
s ¼ 0:201;
F ¼ 1250
ð9Þ
Lipophilicity Assessment of Basic Drugs
87
The difference in logD7.4 versus logkw relationships in the absence and
presence of n-octanol in the mobile phase is illustrated in Figure 1.
The comparison between Equations (8) and (9) further supports the
addition of n-octanol in the mobile phase as a crucial factor, favorable for
the establishment of similar energetics between retention and bulk octanolwater partitioning. In fact, as the stationary phase becomes solvated by the
mobile phase during equilibration, n-octanol, as a lipophilic component,
associates with the stationary phase providing additional masking of the
free silanols and octanol-like character in respect to hydrogen bonding capability. Nevertheless, in the presence of n-octanol, careful consideration should
be taken concerning the range of organic modifier concentration for the extrapolation procedure, since, in water rich mobile phases, it seems to affect more
markedly the linearity of the logk/f relationship leading to downward
curves.[27] Analogous findings were described earlier by Minick et al.[52]
The downward curvature has been observed for both strongly ionized basic
drugs at pH 7.4 and neutral compounds at volume percentages greater than
60% water, depending, however, on the solute as well. Therefore, to avoid
an underestimation of logkw indices, careful selection of the methanol
fraction range was proposed and the use of at least five isocratic logk values
for the extrapolation procedure.[27]
OTHER CHROMATOGRAPHIC DATA AS LIPOPHILICITY
RELEVANT EXPRESSIONS
The slope S of the linear Equation (5) is considered to encode significant information on the lipophilic behaviour of the solute. By some authors, the slope S
Figure 1. Relationships of logD7.4 values versus logkw values, in presence of decylamine (B) and in presence of decylamine þ n-octanol in the mobile phase (D) (Data
taken from reference[27]).
88
C. Giaginis and A. Tsantili-Kakoulidou
is considered to reflect the solute/solvent interactions during the retention
process and is related to the specific hydrophobic surface area.[65] The
strong influence of volume in the slope S was demonstrated for a series of substituted coumarins using PLS analysis.[24] If, within a series of compounds
there are no considerable differences in the forces involved in solute/stationary phase interactions (mainly concerning hydrogen bonding or the extent of
silanophilic interactions), a good relationship between the slope S and the
intercept logkw is anticipated (Equation (10):
S ¼ a logkw þ b
ð10Þ
The organic modifier concentration fo, which produces an equal molar
distribution between the stationary and mobile phase leading to logk ¼ 0,
has been proposed as a measure to rank lipophilicity.[66] The fo indices correspond to the quotient:
fo ¼ logkw =S
ð11Þ
Based on the fo concept, a fast gradient method has been proposed by
Valko et al. to determine the chromatographic hydrophobicity index (CHI)
as a high throughput alternative to the other lipophilicity measures.[67,68]
For this purpose, gradient retention times (tg) are measured and converted
to CHI values by means of a calibration equation, derived by a set of
standards with well determined CHI (fo) values:
CHI ¼ slope tg þ intercept
ð12Þ
The absolute magnitude of the CHI parameter depends on the values
assigned to the set of standards. The method has the advantage that, once
the calibration equation has been established, the retention parameter is
obtained from a single fast gradient run, thus saving time and solvents. The
CHI parameter has been reported to correlate satisfactorily with log P. The
reported chromatographic conditions involve acetonitrile as organic
modifier and the use of ammonium acetate as buffer, without addition of
any masking agent. It should be noted that only few basic compounds were
included in the data set and they were measured at elevated pH mainly as
uncharged species.[69] The effect of organic solvent composition on mobile
phase starting pH and on solutes pKa in gradient chromatography and its consequences in CHI indices of ionisable compounds, has been further investigated.[70] The CHI/pH profile for a number of basic drugs was established
using 2,2,2 trifluoro-ethanol as organic modifier and either ammonium
acetate or butylamine buffer as the aqueous component of the mobile
phase.[71] The focus of that study was to solve the problem of a substantial
drop in pH during gradient elution, especially at high starting pH, which
implies that the neutral form in the case of strong bases cannot be achieved.
The authors suggest the use of 50 mM butylamine as the aqueous
component of the mobile phase to overcome this drawback. In the presence
Lipophilicity Assessment of Basic Drugs
89
of butylamine, minimization of pH variation during gradient elution is
achieved, permitting the determination of the CHI index of the neutral species.
IMMOBILIZED ARTIFICIAL MEMBRANE
CHROMATOGRAPHY
IAM chromatography has been introduced as a promising alternative to
simulate liposome/water partitioning and cell membrane permeation.[17,18,72]
It is prepared by phospholipids covalently bonded to a propylamino silica
support at monolayer densities. Remaining propylamine residues are treated
in a second step to suppress an undesired basic function on the silica
backbone. Moreover, free silanol groups, although not easily accessible,
may interfere in secondary interactions. The most frequently used IAM
column is IAMPC, which contains phosphatidylcholine. In fact, three
different types of IAMPC have been introduced in the market, the single
chain IAMPC-DD, the double chain IAMPC-MG, and IAMPC-DD-2, which
differ on the way the remaining propylamine residues are treated. It is
reported that double chain IAM surfaces better simulate natural phospholipids
and the resulting chromatographic indices correlate better with permeability
data.[73,74]
IAM columns permit the use of aqueous mobile phases without addition
of organic modifier, leading to directly measured logkw values and reducing
considerably the time of analysis. The buffer of choice is phosphate
buffered saline in order to mimic physiological conditions. The pH limitations
of the column restrict measurement in the pH range 2.5 to 7.4. Many authors
prefer the use of pH 7.0, which is close to physiological pH and safer for the
column.[71,73] In the case of compounds with strong affinity for the IAM
surface, acetonitrile up to 30% is preferably added and logkw values are
obtained by linear extrapolation. The use of methanol as organic modifier is
avoided, since it affects the stability of the column, causing methanolysis of
the phospholipids. Nevertheless, the ageing of the column should be
checked from time to time, using standard compounds.[75 – 77]
According to Ong and Pidgeon,[78] partitioning seems to be the principal
retention mechanism in IAM retention, implying that besides hydrophobic
interactions, polar interactions with the solvated layer(s) of the stationary
phases and the head groups of the immobilized phospholipids should be considered. The latter constitute specific electrostatic interactions with ionized
species.[79] Such interactions are very important in the case of protonated
basic compounds, which are more strongly retained as a result of their interaction with the phosphate anions of the stationary phase. It is reported that due
to the involvement of electrostatic forces, the IAM retention of protonated
b-blockers was stronger compared to isolipophilic neutral compounds.[80] In
another study concerning structurally diverse basic and neutral compounds,
the degree of protonation had to be considered in order to obtain a good
90
C. Giaginis and A. Tsantili-Kakoulidou
correlation between logkwIAM and logD values at pH 7.4.[81] Otherwise, a
better correlation was obtained with logP values, implying that the decrease
in the retention due to ionization was compensated by the electrostatic interactions.[80] In the same study, IAM retention was compared to reversed phase
chromatographic retention. Very characteristically, the strong base
metformin, fully protonated at pH 7.4, eluted with the dead time in reversed
phase HPLC, while it was retained in IAM chromatography due to the electrostatic interactions of its positively charged center with the phosphate
anions.[81]
IAM chromatographic indices have been successfully correlated with
liposomes partitioning data; however, the balance between electrostatic and
hydrophobic interactions is considered to differ between the two systems.
Nevertheless, such studies include a rather limited number of
compounds.[82] In the case of basic drugs, silanophilic interactions have
been reported to affect the logkwIAM/pH profile as compared to the corresponding pH/partition diagram in liposomes. Thus, logkwIAM values of propranolol, measured on a double chain IAM.PC.DD2 column, was increased
between pH 6– 7, while in liposome partitioning a plateau was reached at
pH below 8.[76]
The potential of IAM chromatographic indices to predict passive
transport through various biological barriers, as well as to estimate pharmacokinetic properties and certain pharmacological activities, has recently been
reviewed by Barbato.[83] Nevertheless, in a parallel study on the similarity
between IAM columns, conventional HPLC columns, octanol-water partitioning, and biopartitioning systems by means of solvatochromic analysis,
published by Lazaro et al., the belief that IAM chromatography should be considered to be always the best choice for modelling biological processes, is
disputed.[84]
CONCLUSIONS
HPLC provides a user’s friendly, rapid, and compound sparing methodology,
which is successfully applied to determine drug lipophilicity. Although, in the
case of basic drugs, silanophilic interactions may interfere in the partition
mechanism, leading to overestimated or erroneous lipophilicity, there are
ways to reduce such secondary interactions and to obtain extrapolated
retention factors, logkw, which reproduce octanol-water logD values in a satisfactory manner. Nevertheless, protonated bases are considered to develop
specific interactions with biological membranes, which are not encoded in
octanol-water partitioning or reversed-phase chromatographic retention. The
development of IAM stationary phase has opened new perspectives in the
use of HPLC to investigate such interactions in a fast and reproducible way.
The greatest potential of IAM Chromatography is the estimation of passive
Lipophilicity Assessment of Basic Drugs
91
transport and in this aspect it may offer a high throughput screening method
for drug candidates in drug discovery and the development process.
REFERENCES
1. Lipinski, C.A.; Lombardo, F.; Dominy, B.W.; Feeney, P.J. Experimental and computational approaches to estimate solubility and permeability in drug discovery
and development settings. Adv. Drug. Deliv. Rev. 1997, 23, 3 – 25.
2. Testa, B.; Crivori, P.; Reist, M.; Carrupt, P.A. The influence of lipophilicity on the
pharmacokinetic behavior of drugs: Concepts and examples. Perspect. Drug
Discov. 2000, 17, 179–211.
3. Hansch, C. The physicochemical approach to drug design and discovery (QSAR).
Drug Dev. Res. 1981, 1, 267– 309.
4. Smith, D.A.; Jones, B.; Walker, D.K. Design of drugs involving the concepts and
theories on drug metabolism and pharmakokinetics. Med. Res. Rev. 1996, 16,
243– 266.
5. Mannhold, A.; van Waterbeemd, H. Substructure and whole molecule approaches
for calculating log P. J. Comp.-Aid. Mol. Des. 2001, 15, 337– 354.
6. Vrakas, D.; Tsantili-Kakoulidou, A.; Hadjipavlou-Litina, D. Exploring the consistency of logP estimation for substituted coumarins. Quant. Struct.-Act. Relat. 2003,
22, 622– 629.
7. Dearden, J.C.; Bresnen, G.M. The measurement of partition coefficients. Quant.
Struct.-Act. Relat. 1988, 7, 133–144.
8. Hersey, A.; Hill, A.P.; Hyde, R.M.; Livingstone, D.J. Principles of method
selection in partition studies. Quant. Struct.-Act. Relat. 1989, 8, 288– 296.
9. Avdeef, A. pH-Metric logP. 2. Refinement of partition coefficients and ionization
constants of multiprotic substances. J. Pharm. Sci. 1993, 82, 183– 190.
10. Avdeef, A.; Box, K.J. Sirius Technical Application Notes (STAN); Sirius Analytical Instruments Ltd.: Forest Row, UK, 1995; Vol. 2.
11. Braumann, T. Determination of hydrophobic parameters by reversed-phase liquid
chromatography: theory, experimental techniques, and application in studies on
quantitative structure-activity relationships. J. Chromatogr. 1986, 373, 191– 225.
12. Dorsey, J.G.; Khaledi, M.G. Hydrophobicity estimations by reversed-phase liquid
chromatography. Implications for biological partitioning processes.
J. Chromatogr. A 1993, 656, 485–499.
13. Van de Waterbeemd, H.; Kansy, M.; Wagner, B.; Fischer, H. Lipophilicity
measurement by high performance liquid chromatography (RP-HPLC). In Lipophilicity in Drug Action and Toxicology; Pilska, V., Testa, B., and van de
Waterbeemd, H. (Eds.); VCH: Weinheim, 1996; 73.
14. Kaliszan, R. High-performance liquid-chromatographic methods and procedures
of hydrophobicity determination. Quant. Struct.-Act. Relat. 1990, 9, 83 –87.
15. Bechalany, A.; Tsantili-Kakoulidou, A.; El Tayar, N.; Testa, B. Measurement of
lipophilicity indices by reversed-phase high-performance liquid chromatography:
comparison of two stationary phases and various eluents. J. Chromatogr. 1991,
541, 221– 229.
16. Slater, B.; McCormack, A.; Avdeef, A.; Comer, J.E.A. pH Metric logP. 4. Comparison of partition coefficients determined by shake-flask, HPLC and potentiometric methods. J. Pharm. Sci. 1994, 83, 1280– 1283.
92
C. Giaginis and A. Tsantili-Kakoulidou
17. Kramer, S. Liposomes/water partitioning: Theory techniques applications. In
Pharmacokinetic Optimization in Drug Research; Testa, B., van de
Waterbeemd, H., Folkers, G., and Guy, R., Eds.; Verlag Helvetica Chimica
Acta: Zürich and Wiley-VCH: Weinheim, 2001.
18. Pidgeon, C.; Ong, S.; Liu, H.; Qui, X.; Pidgeon, M.; Dantzig, A.H.; Munroe, J.;
Hornback, W.J.; Kasher, J.S.; Glunz, L.; Szczerba, T. IAM chromatography: An
in vitro screen for predicting drug membrane permeability. J. Med. Chem. 1995,
38, 590–594.
19. Ong, S.; Liu, H.; Pidgeon, C. Immobilized artificial membrane chromatography:
measurements of membrane partition coefficient and predicting drug membrane
permeability. J. Chromatogr. A 1996, 728, 113–128.
20. Yamagami, C.; Yokota, M.; Takao, N.; Alan, R. Hydrophobicity parameters determined by reversed-phase liquid chromatography. IX: Relationship between
capacity factor and water-octanol partition coefficient of monosubstituted pyrimidines. Chem. Pharm. Bull. 1994, 42, 907– 912.
21. Tsantili-Kakoulidou, A.; Antoniadou-Vyza, A. Determination of the partition coefficients of adamantyl derivatives by reversed phase TLC and HPLC. In QSAR:
Quantitative Structure-Activity Relationships in Drug Design; Alan R. Liss, Inc.:
1989; 71 –74.
22. Abraham, M.A.; Chadha, H.S.; Leo, A.J. Hydrogen bonding XXXV. Relationship
between HPLC capacity factors and water-octanol partition coefficients.
J. Chromatogr. A 1994, 685, 203– 211.
23. Abraham, M.A.; Chadha, H.S.; Leitao, R.A.E.; Mitchell, R.C.; Lambert, W.J.;
Kaliszan, R.; Nasal, A.; Haber, P. Determination of solute lipophilicity, as
logP(octanol) and logP(alkane) using poly(styrene-divinylbenzene) and immobilized artificial membrane stationary phases in reversed-phase high-performance
liquid chromatography. J. Chromatogr. A 1997, 685, 203– 211; 766, 35 – 47.
24. Vrakas, D.; Panderi, I.; Hadjipavlou-Litina, D.; Tsantili-Kakoulidou, A. Investigation of the relationships between logP and various chromatographic indices
for a series of substituted coumarins. Evaluation of their similarity/dissimilarity
using multivariate statistics. Quant. Struct.-Act. Relat. 2005, 24, 254– 269.
25. Kaliszan, R. High performance liquid chromatographic methods and procedures of
hydrophobicity determination. Quant. Struct.-Act. Relat. 1990, 9, 83 –87.
26. Lombardo, F.; Shalaeva, M.Y.; Tupper, K.A.; Gao, F. ElogDoct: A tool for lipophilicity determination in drug discovery. 2. Basic and neutral compounds. J. Med.
Chem. 2001, 44, 2490– 2497.
27. Giaginis, C.; Theocharis, S.; Tsantili-Kakoulidou, A. Contribution to the standardization of the chromatographic conditions for the lipophilicity assessment of
neutral and basic drugs. Anal. Chim. Acta 2006, 573– 574, 311– 318.
28. El Tayar, N.; van de Waterbeemd, H.; Testa, B. Measurements of protonated basic
compounds by reversed-phase high-performance liquid chromatography: II.
Procedure for the determination of a lipophilic index by reversed-phase highperformance liquid chromatography. J. Chromatogr. 1985, 320, 305– 312.
29. Kaliszan, R. High-performance liquid-chromatographic methods and procedures
of hydrophobicity determination. Quant. Struct.-Act. Relat. 1990, 9, 83 – 87.
30. El Tayar, N.; Tsantili-Kakoulidou, A.; Roethlisberger, T.; Testa, B.; Gal, J.
Different partitioning behaviour of sulphonyl-containing compounds in
reversed-phase high-performance liquid chromatography and octanol-water
systems. J. Chromatogr. 1988, 439, 237– 244.
31. Neue, U.D.; Tran, K.; Iraneta, P.C.; Alden, B.A. Characterization of HPLC
packings. J. Sep. Sci. 2003, 26, 174– 186.
Lipophilicity Assessment of Basic Drugs
93
32. Vervoor, R.J.; Debets, A.J.; Claessens, H.A.; Cramers, C.A.; de Jong, G.J. Optimisation and characterisation of silica-based reversed-phase liquid chromatographic
systems for the analysis of basic pharmaceuticals. J. Chromatogr. A 2000, 897,
1 – 22.
33. Neue, U.D. HPLC columns; Wiley-VCH: Weinheim, 1997.
34. Mendez, A.; Bosch, E.; Roses, M.; Neue, U.D. Comparison of the acidity of residual
silanol groups in several liquid chromatography columns. J. Chromatogr. A 2003,
986, 33–44.
35. Pagliara, A.; Khamis, E.; Trinh, A.; Carrupt, P.A.; Tsai, R.S.; Testa, B. Structural
properties governing retention mechanisms on RP-HPLC stationary phases used
for lipophilicity measurements. J. Liq. Chromatogr. & Rel. Technol. 1995, 18,
1721– 1745.
36. Ascah, T.L.; Kallury, K.M.R.; Szafranski, C.A.; Corman, S.D.; Liu, F. Characterization and high performance liquid chromatographic evaluation of a new AmideFunctionalized Reversed Phase Column. J. Liq. Chromatogr. & Rel. Technol.
1996, 19, 3049– 3073.
37. Stella, C.; Seuret, P.; Rudaz, S.; Carrupt, P.-A.; Gauvrit, P.-A.; Lanteri, P.;
Veuthey, J.-L. Characterization of chromatographic supports for the analysis of
basic compounds. J. Sep. Sci. 2002, 25, 1351 –1363.
38. Neue, U.D.; Cheng, Y.F.; Lu, Z.; Alden, B.A.; Iraneta, P.C.; Phoebe, C.H.;
Tran, K. Properties of reversed phase packings with an embedded polar group.
Chromatographia 2001, 54, 169– 177.
39. Layne, J. Characterization and comparison of the chromatographic performance of
conventional, polar-embedded, and polar-endcapped reversed-phase liquid chromatography stationary phases. J. Chromatogr. A 2002, 957, 149– 164.
40. O’Gara, J.E.; Alden, B.A.; Walters, T.; Petersen, J.; Niederlander, C.; Neue, U.D.
Simple preparation of a C8 HPLC stationary phase with an internal polar functional
group. Anal. Chem. 1995, 67, 3809– 3813.
41. Claessens, H.A.; van Straten, M.A. Review on the chemical and thermal stability
of stationary phases for reversed-phase liquid chromatography. J. Chromatogr. A
2004, 1060, 23 – 41.
42. Stella, C.; Galland, A.; Liu, X.; Testa, B.; Rudaz, S.; Veuthey, J.-L.; Carrupt, P.-A.
Novel RPLC stationary phases for lipophilicity measurement: Solvatochromic
analysis of retention mechanisms for neutral and basic compounds. J. Sep. Sci.
2005, 28, 2350– 2362.
43. Donavan, S.F.; Pescatore, M.C. Method for measuring the logarithm of the
octanol – water partition coefficient by using short octadecyl – poly(vinyl alcohol)
high-performance liquid chromatography columns. J. Chromatogr. A 2002, 952,
47 – 61.
44. Supelco 2001 catalog, Bellefonte, PA, USA 16823 http://www.sigmaaldrich.
com).
45. Liu, X.; Tanaka, H.; Yamauchi, A.; Testa, B.; Chuman, H. Lipophilicity measurement by reversed-phase high-performance liquid chromatography (RP-HPLC): A
comparison of two stationary phases based on retention mechanisms. Helv. Chim.
Act. 2004, 87, 2866 –2876.
46. Cimpan, G.; Irimie, F.; Gocan, S.; Claessens, H.A. Role of stationary phase and
eluent composition on the determination of log P values of N-hydroxyethylamide
of aryloxyalkylen and pyridine carboxylic acids by reversed-phase highperformance liquid chromatography. J. Chromatogr. B 1998, 714, 247– 261.
47. Karger, B.L.; Gant, J.; Hartkopf, A.; Weiner, P.H. Hydrophobic effects in
reversed-phase liquid chromatography. J. Chromatogr. 1976, 128, 65 – 78.
94
C. Giaginis and A. Tsantili-Kakoulidou
48. Buckenmaier, S.M.C.; McCalley, D.V.; Euerby, M.R. Determination of ionisation
constants of organic bases in aqueous methanol solutions using capillary electrophoresis. J. Chromatogr. A 2004, 1026, 251– 259.
49. Espinosa, A.; Bonsch, E.; Roses, M. Retention of ionizable compounds in highperformance liquid chromatography. 14. Acid-base pK values in acetonitrilewater mobile phases. J. Chromatogr. A 2002, 964, 55 – 66.
50. Espinosa, A.; Bonsch, E.; Roses, M. Retention of ionizable compounds on HPLC.
IX. Modeling in reversed phase liquid chromatography as a function of pH and
solvent composition with acetonitrile-water mobile phase. J. Chromatogr. A
2002, 947, 47 – 58.
51. Tsantili-Kakoulidou, A.; El Tayar, N.; Van De Waterbeemd, H.; Testa, B. Structural effects in the lipophilicity of di- and polysubstituted benzenes as measured by
reversed-phase high-performance liquid chromatography. J. Chromatogr. 1987,
389, 33 –45.
52. Minick, D.J.; Frenz, J.H.; Patrick, M.A.; Brent, D.A. A comprehensive method for
determining hydrophobicity constants by reversed-phase high-performance liquid
chromatography. J. Med. Chem. 1988, 31, 1923– 1933.
53. Hansch, C.; Leo, A. Exploring QSAR; ACA: Washington DC, 1995; Vol. 2.
54. Ruiz-Angel, M.J.; Card-Broch, S.; Berthod, A. Ionic liquids versus triethylamine
as mobile phase additives in the analysis of b-blockers. J. Chromatogr. A 2006,
1119, 202– 208.
55. Marszall, M.P.; Baczek, T.; Kaliszan, R. Reduction of silanophilic interactions in
liquid chromatography with the use of ionic liquids. Anal. Chim. Acta 2005, 547,
172– 178.
56. Kaliszan, R.; Marszałł, M.P.; Jan Markuszewski, M.; Baczek, T.; Pernak, J. Suppression of deleterious effects of free silanols in liquid chromatography by imidazolium tetrafluoroborate ionic liquids. J. Chromatogr. A 2004, 1030, 263–271.
57. Schoenmakers, P.J.; Billiet, H.A.H.; de Galan, L. Description of solute retention
over the full range of mobile phase compositions in reversed-phase liquid
chromatography. J. Chromatogr. 1983, 282, 107– 121.
58. Tate, P.A.; Dorsey, J.G. Column selection for liquid chromatographic estimation
of the k0 w hydrophobicity parameter. J. Chromatogr. A 2004, 1042, 37 – 48.
59. Valko, K.; Snyder, L.R.; Glajch, G.L. Retention in reversed-phase liquid
chromatography as a function of mobile-phase composition. J. Chromatogr. A
1993, 656, 501– 520.
60. Shoenmakers, P.J.; Billiet, H.A.H.; de Galan, L. Influence of organic modifiers on
the rentention behaviour in reversed-phase liquid chromatography and its consequences for gradient elution. J. Chromatogr. 1979, 185, 179– 195.
61. Hsieh, M.M.; Dorsey, J.G. Accurate determination of logkw in reversed-phased
liquid chromatography. J. Chromatogr. A 1993, 631, 63 – 78.
62. Unger, S.H.; Cook, J.R.; Hollenberg, J.S. Simple procedure for determining
octanol-aqueous partition, distribution, ionization coefficients by reversed-phase
high performance liquid chromatography. J. Pharm. Sci. 1978, 67, 1364– 1367.
63. Unger, S.H.; Chiang, G.H. Octanol-physiological buffer distribution coefficients of
lipophilic amines by reversed-phase high-performance liquid chromatography and
their correlation with biological activity. J. Med. Chem. 1981, 24, 262– 270.
64. Lombardo, F.; Shalaeva, M.Y.; Bissett, B.D.; Christokhodova, N. Physicochemical and biological profiling in drug research. ElogD7.4 20 000 compounds
later: Refinements, observations and applications. In Pharmacokinetic Profiling
in Drug Research; Testa, B., Kramer, S.D., Wunderli-Allenspach, H. and
Folkers, G. (Eds.); Wiley-VCH: Zurich, 2006; 187 –201.
Lipophilicity Assessment of Basic Drugs
95
65. Horvath, C.; Melander, W.; Molnar, I. Solvophobic interactions in liquid chromatography with nonpolar stationary phases. J. Chromatogr. 1976, 125, 129– 156.
66. Valko, K.; Slegel, P. New chromatographic hydrophobicity index (f0) based on
the slope and the intercept of the log k0 versus organic phase concentration plot.
J. Chromatogr. 1993, 631, 69 – 61.
67. Valko, K.; Bevan, C.; Reynolds, D. Chromatographic hydrophobicity index by
fast-gradient RP-HPLC: A high-throughput alternative to log P/log D. Anal.
Chem. 1997, 69, 2022–2029.
68. Du, C.M.; Valko, K.; Bevan, C.; Reynolds, D.; Abraham, M.H. Rapid gradient RPHPLC method for lipophilicity determination: A solvation equation based comparison with isocratic. Methods. Anal. Chem. 1998, 70, 4228– 4234.
69. Valko, K. Application of high-performance liquid chromatography based measurements of lipophilicity to model biological distribution. J. Chromatogr. A. 2004,
1037, 299– 310.
70. Canals, I.; Valkó, K.; Bosch, E.; Hill, A.P.; Rosés, M. Retention of ionizable
compounds on HPLC. 8. Influence of mobile-phase pH change on the
chromatographic retention of acids and bases during gradient elution. Anal.
Chem. 2001, 73, 4937– 4945.
71. Espinosa, S.; Bosch, E.; Roses, M.; Valko, K. Change of mobile phase pH during
gradient reversed-phase chromatography with 2,2,2-trifluoroethanol-water as
mobile phase and its effect on the chromatographic hydrophobicity index determination. J. Chromatogr. A. 2002, 954, 77 – 87.
72. Taillardat-Bertschinger, A.; Carrupt, P.A.; Barbato, F.; Testa, B. Immobilized artificial membrane HPLC in drug research. J. Med. Chem. 2003, 46, 655– 665.
73. Markovich, R.J.; Stevens, J.M.; Pidgeon, C. Fourier transform infrared assessed of
membrane lipids immobilized to silica: leaching and stability of immobilized artificial membrane – bonded phase. Anal. Biochem. 1989, 182, 237–244.
74. Barbato, F.; di Martino, G.; Grunetto, G.; La Rotonda, M.I. Prediction of drugmembrane interactions by IAM-HPLC: Effects of different phospholipid stationary
phases on the partition of bases. Eur. J. Pharm. Sci. 2004, 22, 261– 269.
75. Morse, K.; Pidgeon, C. Importance of mobile phase in immobilized artificial
membrane chromatography. In Pharmacokinetic Optimization in Drug
Research; Testa, B., van de Waterbeemd, H., Folkers, G., and Guy, R. (Eds.);
Verlag Helvetica Chimica Acta, Zürich and Wiley-VCH: Weinheim, 2001;
351– 381.
76. Rhee, D.; Markovich, R.; Chae, W.G.; Qiu, X.; Pidgeon, C. Chromatographic
surfaces prepared from lyse phosphatidylcholine ligands. Anal. Chim. Acta
1994, 297, 377– 386.
77. Taillardat-Bertschinger, A.; Galland, A.; Carrupt, P.A.; Testa, B. Immobilized artificial membrane (IAM)-HPLC: proposed guidelines for technical optimization of
retention measurements. J. Chromatogr. A 2002, 953, 39– 53.
78. Ong, S.; Pidgeon, C. Thermodynamics of solute partitioning into immobilized artificial membranes. Anal. Chem. 1995, 67, 2119– 2128.
79. Taillardat-Bertschinger, A.; Martinet, C.A.M.; Carrupt, P.-A.; Reist, M.;
Caron, G.; Fruttero, R.; Testa, B. Molecular factors influencing retention on
immobilized artifical membranes (IAM) compared to partitioning in liposomes
and n-octanol. Pharm. Res. 2002, 19, 729– 737.
80. Barbato, F.; di Martino, G.; Grumetto, L.; La Rotonda, M.I. Can protonated betablockers interact with biomembranes stronger than neutral isolipophilic
compounds? A chromatographic study on three different phospholipid stationary
phases (IAM-HPLC). Eur. J. Pharm.Sci. 2005, 25, 379– 386.
96
C. Giaginis and A. Tsantili-Kakoulidou
81. Vrakas, D.; Giaginis, C.; Tsantili-Kakoulidou, A. Different retention behaviour of
structurally diverse basic and neutral drugs in immobilized artificial membrane
(IAM) and reversed-phase HPLC. Comparison with octanol-water partitioning.
J. Chromatogr. A 2006, 1116, 158– 164.
82. Ottiger, C.; Wunderli-Allenspach, H. Immobilized artificial membrane (IAM)HPLC for partition studies of neutral and ionized acids and bases in comparison
with the liposomal partition system. Pharm. Res. 1999, 16, 643 –650.
83. Barbato, F. The use of immobilised artificial membrane (IAM) chromatography for
determination of lipophilicity. Curr. Comp. – Aid. Drug Design 2006, 2, 341– 352.
84. Lazaro, E.; Rafols, C.; Abraham, M.H.; Roses, M. Chromatographic estimation of
drug disposition properties by neans of IAM and C18 columns. J. Med. Chem.
2006, 48, 4861– 4870.
Received February 28, 2007
Accepted April 28, 2007
Manuscript 6172I