![]() |
|
|
Vol. 54, Issue 1, 129-138, July 1998
1-Acid Glycoprotein: Generation of a Three-Dimensional
Quantitative Structure-Activity Relationship Model for Drug Binding to
the A Variant
Service Hospitalo-Universitaire de Pharmacologie de Paris XII, Centre Hospitalier Intercommunal de Créteil, F-94010 Créteil Cedex, France (F.H., J.-C.D., J.-P.T.), Institut de Chimie Thérapeutique Section de Pharmacie, Université de Lausanne, CH-1015 Lausanne-Dorigny, Switzerland (G.C., P.G., N.A.R., A.T.-K., P.-A.C., B.T.), and Laboratoire d'Informatique Médicale, Faculté de Médecine de Dijon, F-21033 Dijon Cedex, France (P.d'A.)
| |
Summary |
|---|
|
|
|---|
Human
1-acid glycoprotein (AAG) is a mixture of at least
two genetic variants: the A variant and the F1 and/or S variant or
variants, which are encoded by two different genes. In a continuation of previous studies indicating specific drug transport roles for each
AAG variant according to its separate genetic origin, this work was
designed to (1) determine the affinities of the two main gene products
of AAG (i.e., the A variant and a mixture of the F1 and S variants) for
35 chemically diverse drugs and (2) to obtain meaningful 3D-QSARs for
each binding site. Affinities were obtained by displacement
experiments, leading to qualitative indications about binding site
characteristics. In particular, drugs binding selectively to the A
variant displayed some common structural features, but this was not
seen for the F1*S variants. Three-dimensional QSAR analyses using the
CoMFA method yielded a steric model for binding to the A variant, from
which a simplified haptophoric model was derived. In contrast, no
statistically sound model was found for the F1*S variants, possibly due
(among other reasons) to an insufficient number of high affinity
ligands in the set.
| |
Introduction |
|---|
|
|
|---|
AAG
is a
protein whose heterogeneity has significant effects on drug transport;
it is the main carrier of basic drugs and other ligands in plasma
(Kremer et al., 1988
). AAG exists as a mixture of two or
three genetic variants (i.e., the A variant and the F1 and/or S
variants) in the plasma of most individuals (Eap and Baumann, 1989
;
Yuasa et al., 1993
). The AAG polymorphism is explained by
the presence of two different genes coding for the protein (Dente
et al., 1987
), of which the AAG-A gene encodes the F1 and S variants and the AAG-B/B' gene encodes the A
variant (Tomei et al., 1989
). The AAG-B/B' gene
is structurally similar to the AAG-A gene but contains 22 base substitutions. Accordingly, the A variant and the F1*S variants
should differ in their amino acid sequences by (at least) 22 residues,
among a total of 181 residues (Schmid, 1975
). In addition, the F1 and S
variants encoded by two alleles of the AAG-A gene (Yuasa
et al., 1993
) should differ by only a few amino acids (less
than five) (Dente et al., 1987
).
Recently, a fractionation method was developed for the AAG variants
(Hervé et al., 1992
, 1993a
), and it was used on a
preparative scale to purify large amounts of the A variant and of a
mixture of the F1 and S variants. The source was a commercial AAG of
human origin containing similar proportions of the three variants.
Large differences in the binding of various drugs have been
demonstrated between the A and F1*S variants, indicating specific drug
transport roles for each AAG variant, according to their separate
genetic origin (Hervé et al., 1993b
, 1996
).
AAG is an important determinant of the plasma binding of many basic
drugs (Routledge, 1986
). In this context, the existence of functional
heterogeneity between the AAG variants, together with variable
concentrations of these variants between individuals (Eap et
al., 1990
), could be responsible for interindividual differences in drug binding, such as those reported by some investigators (Tinguely
et al., 1985
; Eap et al., 1990
). Changes in the
expression of the genetic variants of AAG during inflammation (Eap
et al., 1991
) or illness might also result in altered plasma
drug binding. Finally, the drug binding differences demonstrated
between the A and the F1*S variants indicate that AAG, considered as a
single protein, would present not one but at least two separate and
fractional drug binding sites with different binding specificities. In
addition, because the A-to-F1*S variant ratio varies among individuals, so would the proportions of these sites. The fact that some drugs are
bound only by one of these sites might result in a binding lower than
expected from the total AAG concentration and in an increased risk of
drug interactions due to higher site occupancy.
Given the potential clinical and biological significance of AAG
polymorphism, it would be useful to gain insight into the nature and
topography of the binding site or sites likely to be present on each
AAG variant, with the ultimate goal of predicting affinities and
anticipating possible risk factors. Hence, the first objective of this
work was to improve our understanding of the binding specificity of the
AAG variants using a variety of drugs belonging to different
pharmacological and chemical classes. To this end, displacement
experiments were performed using [3H]imipramine
and [14C]warfarin as selective and high
affinity ligands for the A and F1*S variants, respectively. In the
second part of the work, the binding data collected here and previously
(Hervé et al., 1996
) were used in a 3D-QSAR analysis
using CoMFA, with the objective of obtaining a three-dimensional model
of the binding sites.
| |
Experimental Procedures |
|---|
|
|
|---|
Materials
The A variant and a mixture of the F1 and S variants in the
native (sialylated) form were separated by chromatography from commercial preparations of human AAG [from Cohn Fraction VI; Sigma Chemical, St. Louis, MO (batches 90H9317 and 13H9336], as described previously (Hervé et al., 1992
, 1993a
). The
proportions of the F1, S, and A variants in commercial AAG preparations
are nearly constant (~40%, ~30%, and ~30%, respectively)
(Hervé et al., 1997
). After chromatography, >95% of
each variant was recovered. The composition of the isolated variants
was checked by polyacrylamide gel isoelectric focusing (Eap and
Baumann, 1988
), after small amounts of the different samples had been
desialylated with neuraminidase (Hervé et al., 1993b
).
The F1*S mixture consisted of ~60% F1 and ~40% S, as determined
by laser densitometry. It also contained trace amounts (<4%) of the A
variant. The A variant sample was devoid of any F1*S variants. To
remove possible endogeneous inhibitors, the AAG variant samples were
treated with charcoal, pH 7.4, as described previously (Hervé
et al., 1993b
).
[3H]Imipramine (25 Ci/mmol, 925 Gbq/mmol) and [14C]warfarin (46 mCi/mmol, 1.70 Gbq/mmol) were purchased from Amersham International (Buckinghamshire, UK). The radiochemical purity of the drugs was >98% by thin layer chromatography. Binedaline was a gift from Cassenne (Paris-La Défense, France). Bornaprolol was a gift from Rhône-Poulenc Rorer (Neuilly/Seine, France). Imipramine, desipramine, and clomipramine were gifts from from Ciba-Geigy (Rueil-Malmaison, France). Isradipine was a gift from Sandoz (Basel, Switzerland). Minaprine was a gift from Clin Midy-Sanofi (Paris, France). Propafenone was a gift from Knoll France (Levallois-Perret, France). Tertatolol was a gift from Servier (Jidy, France). Cetirizine was a gift from UCB Pharma France S.A. (Nanterre, France). Desmethylclomipramine was a gift from Dr. C. B. Eap (University Psychiatric Hospital, Prilly-Lausanne, Switzerland). Amitriptyline, auramine O, capsaicin (8-methyl-N-vanillyl-6-nonenamide), chlorpheniramine, diazepam, diphenhydramine, chlorcyclizine, maprotiline, nortriptyline, prazosin, promethazine, (rac)-propranolol, (S)-(-)propranolol, (R)-(+)-propranolol, pyrilamine maleate, quinidine, thioridazine, trazodone, viloxazine, and warfarin were purchased from Sigma.
Binding Experiments
Methods.
AAG concentrations of the charcoal-extracted A and
F1*S variant samples were measured by an immunonephelometric method,
with a Beckman assay kit and apparatus (model 7571, ARRAY TM Protein System; Beckman Instruments, Fullerton, CA). To calculate the molar
concentration of AAG, a molecular mass of 40,000 kDa was assumed
(Kremer et al., 1988
).
|
Evaluation of the displacement data. Two different expressions of a model for the competitive inhibition of the marker (imipramine or warfarin) binding by displacer were used to analyze the displacement data, to determine (1) the ligand inhibitory potencies (IC50) and (2) the ligand association constants (here abbreviated K') for binding to each AAG variant. Each parameter was estimated with its standard deviation by nonlinear regression analysis using a Gauss-Newton algorithm.
Determination of the IC50 parameter. The relationship between the percentage displacement of the respective marker by an inhibitor and the concentration of the inhibitor was modeled with eq. 1, assuming a pseudo-Hill coefficient of 1:
|
(1) |
Determination of the K' parameter.
These
equations were used to determine the association constant
(K') of each inhibitor (Hervé et al.,
1996
):
|
(2) |
|
(3) |
CoMFA Study
Binding data. Input data (dependent variables) were pIC50 values (i.e., minus the logarithmic value of IC50). For the chiral ligands studied as racemates, the same affinity was attributed to both enantiomers. This assumption was necessary in the absence of information on the binding stereoselectivity to the separate AAG variants. The effects of this assumption on CoMFA results are discussed separately for each variant (see below).
Technical specification.
All molecular modeling calculations
were performed with the Sybyl software (ver. 6.3; Tripos Associates,
St. Louis, MO) running on Silicon Graphics Indy R4400 and
O2 R5000 workstations. Ligand starting geometries
were built in their neutral form by the Concord algorithm (Tripos
Associates). Energy minimization was performed with the Tripos force
field including the electrostatic energy term calculated using
Gasteiger and Marsili (1980)
atomic charges. The method of Powell was
used for minimizations, with convergence being reached when the
gradient decrease was <0.001 kcal/mol/Å. The DISCO module of Sybyl
was used in combination with quenched molecular dynamics as described
previously (Caron et al., 1997
).
= 2) to retain models with q2 > 0.4 (Gaillard et al., 1994
= 0). The other options were chosen according to
published standards (Simon, 1992Alignment of ligands. For the A variant, manual ligand alignments were performed using the optimized geometries modified by manual geometrical fitting. All molecules were aligned on the template amitriptyline. Amitriptyline was chosen because it was the most rigid among compounds with high binding affinity. The geometry of the template was optimized by energy minimization.
The alignment of ligands was based on their aromatic moieties (often two) and basic nitrogen: (1) the aromatic moieties were aligned by their centroid and the normal to their plane, and (2) the free-electron lone pairs of the basic nitrogen were pointed in the same direction. Due to their small structural similarity with other ligands, diazepam and progesterone could not be aligned satisfactorily and had to be excluded from the CoMFA study. For the F1*S variants, the DISCO approach was used to find alignments (Martin et al., 1993Presentation of Results
The graphic results present the most relevant regions of space where the variations of a chosen statistical field are the largest. In each grid point, a statistical field is defined as the product of its coefficient in the PLS equation by the standard deviation.
| |
Results |
|---|
|
|
|---|
Ligand Displacement Studies
The results of the displacement of [3H]imipramine bound to the A variant and of [14C]warfarin bound to the F1*S variant mixture showed that the displacing ligands displayed very different inhibitory potencies. Examples of displacement of imipramine from the A variant and of warfarin from the F1*S variants are shown in Figs. 1, A and B, respectively. The theoretical curves representing a competitive inhibition of the marker binding fitted the data well, with a few exceptions discussed below. The IC50 values of the ligands calculated from the displacement curves are listed in Table 1. The values were in the range of 5.6 µM (11, propafenone) to 673 µM (32, minaprine) for [3H]imipramine binding to the A variant and in the range of 5.0 µM (31, thioridazine) to 814 µM (22, desipramine) for [14C]warfarin binding to the F1*S variants.
|
Amitriptyline, a tricyclic antidepressant drug (19;
IC50 = 6.1 µM), was one of the most
potent displacers of [3H]imipramine. The
tricyclic antidepressants as a group (compounds 19 and
21-26) were potent displacers of imipramine but
had poor or no inhibitory potency toward warfarin (Table 1). Several
antihistaminic drugs (compounds 4, 6, and 7) also seemed to displace preferentially imipramine bound to the A variant. In contrast, other ligands, such as
blockers (compounds 14-16; see Table 1), had comparable inhibitory potencies for the A and F1*S variants. Plotting the ligand
inhibitory potencies to F1*S variants versus those to the A variant
revealed the absence of any correlation between these two sets of data
(result not shown).
A few drugs did not inhibit competitively the binding of the marker (Table 1). The small displacement of bound radiomarker measured after the addition of these drugs (i.e., the percentage of displacement did not exceed 16 ± 6% at the highest inhibitor concentrations) (data not shown) was assigned to nonspecific binding.
To compare further the binding specificities of the A and F1*S variants, the binding association constants (K') of the inhibitors were calculated from the displacement data using a model for competitive binding at varying numbers of identical protein sites (see Experimental Procedures). The calculated values of K' are shown in Table 1. This table also shows a selectivity factor (S), namely the ratio of affinities of a ligand for the A variant to those for the F1*S variants. The values of K' ranged from 4.18 × 106 liter/mol (11, propafenone) to 0.015 × 106 liter/mol (32, minaprine) for the A variant and from 6.00 × 106 liter/mol (31, thioridazine) to 0.018 × 106 liter/mol (22, desipramine) for the F1*S variants. It is of interest to note that the ranking of the ligands according to the values of K' in Table 1 closely parallels the ranking according to the IC50 values. Furthermore, back-calculation from the binding parameters and total reagent concentrations gave theoretical marker ligand displacements that correlated closely with the experimental values (r > 0.95), showing that the model used was able to account for the experimental results. The only exceptions concerned the few low affinity ligands indicated in Table 1.
CoMFA
A variant. A number of alignments of ligands were tested. The alignment shown in Fig. 2 was retained for the final analyses because it afforded the best statistical results and the best superposition of haptophoric elements. The many CoMFA models obtained were compared by means of a plot of q2 versus N (Fig. 3). This indispensable procedure identified the best CoMFA models to be submitted to final PLS analysis. As shown in Fig. 3, the electrostatic field was poorly correlated with affinity and decreased q2 when combined with other fields. The final model therefore was free from electrostatic contributions. The lipophilic field alone gave no model with q2 > 0.4 and even lowered the q2 of the steric field alone when combined with it. This led us to believe that the variability in affinity is explainable mainly by steric properties of the ligands. However, and according to general CoMFA rules, a final analysis was performed for all models with a value of q2 > 0.4. To avoid possible artifacts due to the mixing of several molecular fields, the lipophilic model also was submitted to a final analysis.
|
|
|
F1*S variants. Should the F1*S variants be characterized by a single binding mode different from that of the A variant, a CoMFA alignment successful in modeling binding to the A variant would fail for the F1*S variants? This was tested with the manual alignment that was successful for the A variant used for the F1*S binding data (i.e., steric model 1). No model with q2 > 0.4 was obtained.
The DISCO approach was used to find other alignments (Martin et al., 1993| |
Discussion |
|---|
|
|
|---|
Binding Studies
Our previous studies with smaller sets of ligands have
demonstrated differences and similarities in the ligand-binding
properties of the A and F1*S variants (i.e., the two main
AAG gene products) (Hervé et al., 1993b
,
1996
). The results presented here both extend the knowledge of the AAG
variants binding properties and allowed us to model the high affinity
site of the A variant.
Although the IC50 values calculated here are
independent of the binding model, the K' values and the site
modeling assume competitive binding at a single-site class with
insignificant nonspecific binding. This is justified on several
grounds. Theoretical calculations from known binding parameters for
both high and low affinity sites showed that under the conditions used
here, <6% of the binding of marker ligand would be at the low
affinity sites. The good correlation between the experimentally
determined displacements and theoretical marker ligand displacements
back-calculated from the binding parameters and total reagent
concentrations demonstrated that the model adequately fitted the
experimental results. Further evidence for the validity of the model is
shown by the parallel ranking of the inhibitor K' and
(model-independent) IC50 values. Also, previous
binding studies with some of the ligands used here as inhibitors (i.e.,
compounds 1, 9, 10, 18,
30, and 33 described in Table 1) showed that the
stochiometry of the interactions between the A and the F1*S variants
with their respective specific ligands was ~1, suggesting that these
ligands share a single common binding site on each of these variants
(Hervé et al., 1996
). Furthermore, the association constants (K') estimated for these ligands from the
displacement experiments based on the assumption of competitive
inhibition were comparable to the association constants
(Ka) determined in the direct binding
studies. For the current work, the direct binding of (radiolabeled)
propranolol also was determined (data not shown), and it gave
association constants close to the K' values in Table 1
(Ka = 0.22 and 0.21 × 106 liter/mol for the A and F1*S variants,
respectively).
The current results allow insight into the nature and topography of the binding site or sites likely to be present on each AAG variant. Using a large number of medicinal ligands, we bring evidence that the drugs with selective affinity for the A variant share marked structural similarities, most of them containing a basic amino group linked by a short carbon chain to two aromatic rings that are either bridged to form a tricyclic structure (e.g., amitriptyline, nortriptyline, imipramine, desipramine, maprotiline, and promethazine) or unbridged (e.g., disopyramide, methadone, and diphenhydramine). However, other analogues of imipramine (e.g., clomipramine and desmethylclomipramine) or diphenhydramine (e.g., chlorcyclizine, cetirizine, and chlorpheniramine) or promethazine (e.g., chlorpromazine) had little or no selectivity for the A variant.
These results suggest that substitution of an aromatic ring with a chlorine atom or the presence of a piperidine ring may be factors leading to decreased ligand binding selectivity to the A variant. These findings indeed are made explicit by the CoMFA model (see below).
In contrast, no clear analogies were found for the selective ligands of
the F1*S variants (e.g., dipyridamole, warfarin, binedaline, and
prazosin), except for a relatively large hydrophobic ring system
(Hervé et al., 1996
). This suggests that the ligand
binding site of the F1 and S variants could be a relatively large
hydrophobic pocket able to accommodate a variety of chemical
structures, whereas the A variant binding site seems to be of smaller
size and of greater ligand specificity. In addition, the A and F1*S
variants display similar affinities for a variety of other ligands,
suggesting that these share some common structural features or
physicochemical properties with other ligands selectively bound by the
A or F1*S variants. This preliminary hypothesis then was challenged by
a 3D-QSAR study using the CoMFA approach.
CoMFA
Graphic representation of the model for the A variant. The graphic results of the best CoMFA model obtained (Table 2, model 1) are represented in Fig. 4 with (R)-bornaprolol (a high affinity ligand) in position and in Fig. 5 with minaprine (a low-affinity ligand) in position. A number of sterically unfavorable regions (dark gray) can be seen located around and below the aromatic ring R1 (Fig. 5A) and below the basic nitrogen (Fig. 5B), indicating forbidden positions. In contrast, there are favorable influences (light gray) above the second ring R2 (Fig. 5C) and behind the basic nitrogen (Fig. 5D). This is reasonable given that (1) zone A was deduced from the low affinity of compounds containing a chlorine atom as substituent in the aromatic ring, (2) zone B explains the low affinity of ligands containing a piperidine ring, (3) zone C was deduced from the relative orientation of the two rings present in high affinity tricyclic drugs, and (4) zone D explains the high affinity of N-methylated amines.
|
|
|
F1*S variants.
The impossibility of finding a statistically
acceptable (q2 > 0.4) 3D-QSAR model for binding
to the F1*S variants using the same alignment as for the A variant
confirms that their binding site is different from that of the A
variant, as found by Raes et al. (1987)
. Furthermore, the
fact that no other alignment gave a 3D-QSAR model suggests a number of
different, partly incompatible explanations. The first is that the F1
and S variants have different binding sites and that the experimental
affinity data are composite values. However, the existence of both a
high degree of primary sequence homology (Dente et al.,
1987
) and similar ligand binding properties between the individual F1
and S variants (Hervé et al., 1993a
and 1993b
) does
not support this hypothesis.
Conclusion
To the best of our knowledge, this study represents the first
successful attempt to generate a homogeneous, large set of data on the
binding of drugs to variants of a transport macromolecule and to obtain
a 3D-QSAR model from some of these data. Until the three-dimensional
structure of the AAG variants is revealed through direct investigations
(e.g., NMR or X-ray crystallography), the CoMFA approach is one the
best tools available to obtain indirect information on binding sites.
The results reported here add to our fundamental knowledge of AAG
variants, and they may well be of clinical significance given the
genetic polymorphism existing for AAG in humans (Schmid, 1975
; Dente
et al., 1987
; Yuasa et al., 1993
) and the clear
importance of this protein in the plasma binding of drugs (Routledge,
1986
).
| |
Footnotes |
|---|
Received August 1, 1997; Accepted March 16, 1998
1 Current affiliation: Department of Chemistry, University of Malaya, 50603 Kuala Lumpur, Malaysia.
B.T. and P.A.C. are indebted to the Swiss National Science Foundation for support. F.H., J.C.D., and J.P.T. are indebted to the Ministère de l'Education Nationale (EA 427) and to the Réseau de Pharmacologie Clinique for support. N.A.R. acknowledges receipt of a JWT Jones Travelling Fellowship given by the Royal Society of Chemistry. A.T.K. is indebted to the Foundation Herbette (University of Lausanne) for a travel grant. F.H. and G.C. contributed equally to this study.
Send reprint requests to: Prof. Bernard Testa, School of Pharmacy, BEP, University of Lausanne, CH-1015 Lausanne, Switzerland. E-mail: bernard.testa{at}ict.unil.ch
| |
Abbreviations |
|---|
AAG, human
1-acid
glycoprotein;
CoMFA, comparative molecular field analysis;
DISCO, distance comparison;
PLS, partial least squares;
QSAR, quantitative
structure-activity relationships;
SAMPLS, sample-distance partial least
squares;
3D-QSAR, three-dimensional quantitative structure-activity
relationships.
| |
References |
|---|
|
|
|---|
1-acid glycoprotein: evidence for drug-binding differences between the variants and for the presence of two separate drug-binding sites on
1-acid glycoprotein.
Pharmacogenetics
6:
403-415[Medline].
1-acid glycoprotein preparation using chromatography on immobilized metal affinity absorbent and on hydroxyapatite.
J Chromatogr B
688:
35-46.
1-acid glycoprotein in the native form by chromatography on an immobilized copper(II) affinity absorbent: heterogeneity of the separate variants by isoelectrofocusing and by concanavalin A affinity chromatography.
J Chromatogr
615:
47-57[Medline].
3 receptor.
J Med Chem
37:
4109-4117[Medline].
1-acid glycoprotein (orosomucoid) variants.
Hum Genet
84:
89-91[Medline].This article has been cited by other articles:
![]() |
G. L. Hortin, N. Seam, and G. T. Hoehn Bound Homocysteine, Cysteine, and Cysteinylglycine Distribution between Albumin and Globulins Clin. Chem., December 1, 2006; 52(12): 2258 - 2264. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Nishi, N. Fukunaga, and M. Otagiri CONSTRUCTION OF EXPRESSION SYSTEM FOR HUMAN {alpha}1-ACID GLYCOPROTEIN IN PICHIA PASTORIS AND EVALUATION OF ITS DRUG-BINDING PROPERTIES Drug Metab. Dispos., October 1, 2004; 32(10): 1069 - 1074. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Taheri, L. P. Cogswell III, A. Gent, and G. R. Strichartz Hydrophobic and Ionic Factors in the Binding of Local Anesthetics to the Major Variant of Human alpha 1-Acid Glycoprotein J. Pharmacol. Exp. Ther., January 1, 2003; 304(1): 71 - 80. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Boffito, K. Sciole, R. Raiteri, S. Bonora, P. G. Hoggard, D. J. Back, and G. Di Perri alpha 1-Acid Glycoprotein Levels in Human Immunodeficiency Virus-Infected Subjects on Antiretroviral Regimens Drug Metab. Dispos., July 1, 2002; 30(7): 859 - 860. [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||