|
|
|
|
Unità di Farmacologia, Dipartimento Farmacobiologico (A.D.L., S.T., M.D.B., J.-F.D., D.C.C.) and Dipartimento Farmacochimico (G.L., F.C., A.S., C.F., V.T.), Facoltà di Farmacia, Università di Bari, Bari, Italy
Received December 24, 2002; accepted July 7, 2003.
| Abstract |
|---|
|
|
|---|
[N-(2,6-dimethylphenyl)pyrrolidine-2-carboxamide (To5)] or
[N-(2,6-dimethylphenyl)pyrrolidine-3-carboxamide (To9)] in a proline-like cycle and/or linked to a lipophilic benzyl moiety as in N-benzyl-tocainide (Benzyl-Toc), 1-benzyl-To5 (Benzyl-To5), and 1-benzyl-To9 (Benzyl-To9). INa were elicited with pulses to -20 mV from different holding potentials (-140, -100, and -70 mV) and stimulation frequencies (2 and 10 Hz). All compounds were voltage-dependent and use-dependent channel blockers. The presence of a proline-like cycle increased the potency; i.e., To5 was 3- and 10-fold more effective than Toc in blocking INa at the holding potential of -140 and -70 mV, respectively. The benzyl group on the amine further enhanced drug effectiveness with the following scale: Benzyl-To9
Benzyl-Toc > Benzyl-To5. At a holding potential of -100 mV and 10-Hz stimulation, Benzyl-To9 blocked INa with a half-maximal concentration of 0.5 µM, being 60 and 400 times more potent than To9 and Toc, respectively. The similar effectiveness of Benzyl-Toc and Benzyl-To9 was paralleled by a similar spatial arrangement by equilibrium geometry modeling. In addition, the latter had a higher pKa value that probably contributed to a slow kinetic during its high use-dependent behavior. Benzyl-To5 had its lowest energy level at a more folded conformation that justifies the less favorable profile among the N-benzylated analogs. The new compounds are the most potent tocainide-like sodium channel blockers so far described and have high therapeutic potentials.
subunit, as critical for LA binding and activity on Na+ channels of various excitable tissues (Ragsdale et al., 1994
and
-
interactions, respectively, with the two hydrophobic domains cited above (Dougherty, 1996
-
interaction stronger than the
-cation interaction made by the amino group (Wright et al., 1998
proline-like analog of tocainide, conventionally named To5, showed an increased stereoselectivity. In parallel, the eutomer (R)-To5 was 5- and 20-fold more potent than tocainide in producing a tonic and a 10-Hz use-dependent block of sodium currents, respectively (Talon et al., 2001
-proline analog To9) to maintain the spatial structure of the molecule and improve basicity; 2) by increasing the lipophilicity on the nitrogen atom by introducing a benzyl moiety on both tocainide (Benzyl-Toc) and its proline-like analogs (Benzyl-To5 and Benzyl-To9) (Fig. 1). The rationale for testing the effect of the original N-benzylation versus the more classic N-alkylation was aimed to combine the effect of a second aromatic moiety, able to establish high-affinity interaction with one of the two aromatic residues, with the change in basicity of the pharmacophore. This latter can influence drug potency; in fact, charged amino moiety can influence drug diffusion through the membrane and/or establish hydrogen bonds with the third site in D1S6 (Nau et al., 1999
|
We found that the N-benzyl analogs showed a remarkable gain of potency; Benzyl-To9 was up to 400 times more potent than tocainide as an inactivated use-dependent channel blocker. A conformational analysis proved that the increase in affinity parallels the spatial disposition of the two aromatic rings, corroborating their direct role in establishing specific interactions with the binding site. An optimal basicity of the amino group, influenced by its substituents, confers a high use-dependent behavior. The results shed further light on the molecular interaction of LAs with their receptor. Also, the new potent agents deserve further investigation for their wide potential therapeutic application, because previous studies indicate a good correlation between the scale of potency determined by this approach and the in vivo and in vitro antimyotonic activity of the drugs in animal models of the human disease (De Luca et al., 1997b
; Talon et al., 2001
).
| Materials and Methods |
|---|
|
|
|---|
Sodium currents were recorded using an amplifier connected via an analog-to-digital and digital-to-analog Digidata 1200 Interface (Axon Instruments, Union City, CA) to a 486 DX2/66 personal computer and stored on the hard disk. The stimulation protocols and data acquisition was driven by the Clampex program (pClamp6; Axon Instruments). The currents flowing in response to depolarizing command voltages were low-pass filtered at 10 kHz (Frequency Devices, Haverhill, MA), visualized on an oscilloscope, and sampled at 20 kHz. The acquired traces were analyzed later using Clampfit program (pClamp6 and 8.1; Axon Instruments).
Drugs and Solutions. The following solutions were used: normal physiological solution, 115 mM NaCl, 2.5 mM KCl, 1.8 mM CaCl2, 2.15 mM Na2HPO4, and 0.85 mM NaH2PO4; external solution, 77 mM NaCl, 38 mM choline-Cl, 2.5 mM KCl, 1.8 mM CaCl2, 2.15 mM Na2HPO4, and 0.85 mM NaH2PO4; internal solution, 105 mM CsF, 5 mM MOPS, 2 mM MgSO4, 5 mM EGTA, and 0.55 mM Na2ATP. The pH was adjusted to 7.2 with a standard NaOH concentrated (5 N) solution.
The compounds tested and shown in Fig. 1 were tocainide (Toc), N-benzyl-tocainide (Benzyl-Toc); the
-proline derivatives N-(2,6-dimethylphenyl)pyrrolidine-2-carboxamide (To5) and its benzyl derivative 1-benzyl-N-(2,6-dimethylphenyl)pyrrolidine-2-carboxamide (Benzyl-To5); and the
-proline derivatives N-(2,6-dimethylphenyl)pyrrolidine-3-carboxamide (To9) and its benzyl derivative 1-benzyl-N-(2,6-dimethylphenyl)pyrrolidine-3-carboxamide (Benzyl-To9).
All the compounds were prepared in our laboratories as hydrochloride or hydroiodide salts according to procedures described in details elsewhere (Franchini et al., 2000
). All compounds were racemates, with the exception of Toc and To5, synthesized as pure R- and S-enantiomers (Franchini et al., 2000
). Racemates of these latter to be tested on INa were obtained by exactly mixing equimolar solutions of the single enantiomers.
Stock solutions of Toc, To5, and To9 were prepared by dissolving the compounds in external solution, whereas stock solutions in external solution containing dimethyl sulfoxide (<1%) were used for the benzyl-derivatives. Dimethyl sulfoxide at the highest concentration used for dilution (0.2%) was without effect on the parameters recorded. All other chemicals used were of analytical grade and obtained from Sigma Chemical Co. (St. Louis, MO).
Pulse Protocols and Statistical Analysis. In agreement with previous studies (Nau et al., 1999
; De Luca et al., 2000
; Yarov-Yarovoy et al., 2002
), the voltage-dependent block exerted by the compounds was evaluated on nearly maximal INa-elicited single depolarizing 10-ms test pulses to -20 mV from two different HP: -140 and -70 mV. The evaluation of concentration-dependent drug effect, expressed as half-maximal blocking concentration (IC50), at -140 mV allowed calculation of the affinity constant for the resting state (Kr); the amount of closed-state inactivation at this potential was negligible. The IC50 value calculated from the HP of -70 mV was instead influenced both by the higher proportion of channels entering a closed-state inactivation at this potential and by the ability of the drug to modify this proportion in favor of more channels being inactivated, based on the ability of the drug to act as inactivated channel blocker. Because the inactivated state was nonconductive, the voltage-dependent block exerted by the drugs was used to estimate the affinity constant for the inactivated channel (Ki) using the above IC50 values, the relative distribution of channels in the resting and/or inactivated states at the holding potentials used (from the steady-state inactivation curve) and the equations described below.
The use-dependent block exerted by each compound was evaluated as the cumulative block of INa upon increase of the stimulation frequency. For this protocol, the membrane potential was held at -100 mV, a HP closer to the physiological values at which almost all the channels are in the resting closed state (De Luca et al., 2000
). Trains of depolarizing 10-ms pulses to -20 mV from the HP were applied at the frequency of 2 and 10 Hz for 30 s. At both frequencies, the drug-induced reduction of peak INa at the first pulse was considered tonic block (block of the channel in the resting state). The use-dependent behavior was estimated by the further reduction of the current observed in the presence, but not the absence, of use-dependent compounds that progressively cumulated over the tonic block in a frequency-dependent manner until a new equilibrium was reached. The value of the current at the equilibrium normalized with respect to the current in the absence of drug was used to calculate the absolute potency of the drug for blocking the channels under conditions of excessive stimulation (e.g., high-frequency firing). The normalization of each current trace composing the 10-Hz train in the presence of the drug with respect to the current at the first pulse after drug, allowed evaluation, for each compound, of the concentration-dependent time course to reach the steady state for use-dependent blockade. For any drug and at each drug concentration, this process was well described with a single exponential decay (
2 < 0.0005), which led to estimate the concentration-dependent time constants as well as the fractional use-dependent block (i.e., the block of INa by the use-dependent compound eliminating the tonic block component).
Steady-state inactivation (h
) curves were determined by cyclic protocol of pulse sequences. Each sequence consisted of a conditioning pulse to -140 mV for 500 ms (to have most of the sodium channels in the "activatable" state), a prepulse of variable potential of 1000-ms duration, and the 10-ms test pulse to -20 mV; after a pause of 1 s, the sequence was cyclically repeated 18 to 20 times, with the prepulse potential value increased each time by 5-mV steps (De Luca et al., 1997a
, 2000
).
Modeling Studies. All the analogs tested are weak basic amines, and it is likely that in the physiological environment, they are largely present in the protonated charged form. This plays an important role for their pharmacological profile, in that they are able to influence use-dependent behavior and gain access to the receptor through the hydrophilic pathway of the open channel by sensing the electric gradient across the membrane (Bean et al., 1983
; Liu et al., 2003
). Thus, the modeling studies were conducted on the protonated cationic form. Models were constructed by fragments, assuming the starting geometry to be similar to that found for lidocaine hydrochloride monohydrate crystal structure (Hanson and Röhrl, 1972
). Equilibrium geometries were calculated starting with a systematic conformer distribution analysis. Conformers were grouped into families based on relevant torsion angle values. The best (i.e., most stable) representative of each family was submitted to an AM1 semiempirical geometry optimization and a single-point energy calculation by a Hartree-Fock calculation at the 3-21G(*) level. On the best conformer of each analog, a 3-21G(*) Hartree-Fock geometry optimization was performed. Energies were corrected for aqueous solvation using the Cramer/Truhlar SM5.4 model (Chambers et al., 1996
). The pKa values of the newly synthesized compounds were calculated by using the approach proposed by D'Souza et al. (2000
). This method is based on the calculation of the proton affinity of a nitrogen base. Proton affinity corresponds to the energy of a reaction in which a nonbonded electron pair is replaced by a bonded pair. A good estimation of the energy for this type of reaction may be derived calculating the energy difference (
E) between the lowest energy conformations of the protonated and unprotonated forms. This energy has been shown to be related to the pKa value of the base for several nitrogen bases (D'Souza et al., 2000
). Running the same calculations on a series of nitrogen bases with different experimental pKa values, a straight-line plot is obtained when plotting
E versus the experimental pKa values. Hence, the unknown experimental pKa value for a nitrogen base may be derived by calculating the relative energy of protonation (i.e.,
E) and interpolating the latter in the curve. Thus, the 3-21G(*) estimated free energy difference (
E) between the protonated and unprotonated forms of eight pharmacologically relevant nitrogenous bases (experimental pKa values in the range 6.9-10.1) was calculated, and a straight-line plot correlating experimental pKa values and
E was derived (
E = 3.7632 pKa + 268.22, r = 0.97,
pKa = 0.3). By interpolating the
E values calculated for the newly synthesized compounds, the corresponding estimated pKa values were obtained. All calculations and graphical representations were performed by using the SPARTAN PRO software package (Wave-function, Inc., Irvine, CA).
Data Analysis and Statistics. The data were expressed as mean ± S.E.M. The estimates of S.E.M. of normalized INa values have been obtained as described previously (De Luca et al., 1997a
, 2000
). Molar concentrations of the drugs tested that produce a 50% block of INa (IC50) in the various experimental conditions were determined by using a nonlinear least-squares fit of the concentration-response curves to the following logistic equation: Effect =-100/[1 + (K/[drug])n], where Effect is the percentage change of INa, -100 is the maximal percentage block of INa, K is IC50, n is the logistic slope factor, and [drug] is the molar concentration of the compound. The h
curves have been fitted with a single Boltzmann distribution, and the potential at which 50% of the sodium channels were inactivated (Vh1/2) was calculated at the inflection point of the curves (De Luca et al., 1991
). The left shift of the h
curve produced by the drug was a function of the concentration used and of the relative affinity for the channel in the inactivated state. The difference between Vh1/2 at each drug concentration and the control Vh1/2 (
Vh1/2) measured in the same fiber, as a function of the drug concentration, were well fitted by the nonlinear equation
Vh1/2 = k ln{(1 + [drug]/Kr) (1 + [drug]/Ki)-1}, where k is the mean slope factor of the h
curves at each drug concentration [drug] (taken constant to 5.5), Kr is the affinity constant for resting state obtained from the fit of the concentration-response curve at -140 mV, and Ki is the apparent affinity of drug for the inactivated state (Bean et al., 1983
). A more detailed evaluation of absolute Ki was obtained from the voltage-dependent distribution of the channels in the resting (h) and inactivated state (1 - h) according to the h
curve, using the equation 1/K-70 = h/Kr + (1 - h)/Ki, where K-70 and Kr are the IC50 values obtained from dose-response curves at -140 and -70 mV, whereas h and (1 - h) represent the fraction of channel present at resting and inactivated state at -70 mV, respectively (Bean et al., 1983
).
Statistical significance of differences between couples mean values has been estimated by unpaired Student's t test and considered significant with p < 0.05. The statistical significance between IC50 values ± S.E. obtained from the fit was also evaluated by a Student's t distribution using a number of degrees of freedom equal to the total number of preparations determining each point of the curve minus the number of means determining the curve minus two for the free parameters (De Luca et al., 2000
).
CLogP is the calculated log value of the octanol/water partition coefficient obtained by using C Log Software v. 3.0 (Biobyte Corp., Clermont, CA), and logD is the calculated log value of distribution coefficient at physiological pH, considering the nature of weak electrolytes and the different lipophilicity of unprotonated and protonated forms. LogD is then calculated at pH 7.4 from the cLogP and pKa values of each compound according to the equation logD = LogP - log (1 + 10pKa-pH).
Correlation analysis was evaluated by fitting the experimental data points to linear regression analysis. Nonlinear equation fitting and processing for data graphics were done by Fig P Software (Bio-soft, Cambridge, UK).
| Results |
|---|
|
|
|---|
-proline derivative To5 at 100 µM produced a block of INa that was comparable with that produced by 500 µM Toc. Interestingly, the
-proline-analog To9 showed almost the same potency as To5, suggesting that the increased distance of the amino group by one methylene does not greatly modify the interaction of this constrained molecule with the receptor. In contrast, all benzylated derivatives were markedly more potent that the related unbenzylated parent compounds, with higher potency shown by Benzyl-To9, which produced a substantial block at concentrations as low as 1 µM (Fig. 2).
|
For each compound, we constructed concentration-response curves for the block produced at both HP (Fig. 2). From the fit of the experimental points, we obtained the concentrations for half-maximal block of INa (IC50) at both -140 and -70 mV. Although the IC50 value at the former HP can be considered as the affinity constant for closed-resting (Kr) channels, the IC50 at -70 mV (K-70) is strongly influenced by the relative proportion of resting channels in equilibrium with closed-inactivated ones and by the ability of the drug to influence such an equilibrium based on its affinity for the inactivated state (Table 1) (Nau et al., 1999
). In fact, as expected from inactivated channel blockers, a great increase of potency was observed when the membrane potential was held at -70 mV: the concentration-response curves of all compounds were clearly shifted to the left with respect to those obtained at -140 mV. The gain of potency for each compound at the more depolarized potential ranged between 4- and 15-fold; the smallest gain was observed with tocainide and the greatest with To5. When looking at the absolute values of IC50, we confirmed that To5 and To9 were almost equipotent, both being about 3- and 10-fold more potent than Toc at -140 and -70 mV, respectively (Table 1). As anticipated, a marked change in the potency of the compounds was observed when a benzyl group was introduced on the amine one (Table 1). Figure 3 shows the potency of each compound in comparison with that of tocainide at both holding potentials. As can be seen, the scale of potency was Benzyl-To9 > Benzyl-Toc > Benzyl-To5 > To5 = To9. Benzyl-To9 was up to 200- and 400-fold more potent than tocainide for block at -140 and -70 mV, respectively (Table 1). Focusing on the specific effect on potency caused by the simple introduction of the benzyl moiety, we found that the gain of potency observed was different for each benzylated compound although more potent than the relative nonbenzylated ones. The most dramatic increase was observed for Benzyl-Toc, which was 80- and 150-fold more potent than tocainide in determining the block at -140 and -70 mV, respectively. Benzyl-To9 showed at both potentials a similar 40-fold increase in potency compared with its parent compound (Table 1). Taking into account that To9 is already significantly more potent than Toc, its benzylated analog turned out to be the most potent tocainide-like compound described so far, with half-maximal concentration at -70 mV as low as 0.5 µM. The introduction of the benzyl group on the
-proline derivative To5 also markedly increase the potency; however, such an increase was significantly less with respect to that observed with Toc and To9 (Table 1).
|
|
The activity of the compounds as inactivated channel blockers was confirmed by evaluating their effects on the voltage dependence of steady-state channel availability by constructing h
curves in the absence and in the presence of the test compounds. All the compounds produced a shift of the h
curve and consequently of the potential for inactivating 50% of the channels (Vh1/2) toward more negative potentials, corroborating the view that the test compounds preferentially bind and stabilize the channel in the inactive state (Fig. 3). The shift of the h
curve produced by each drug was clearly concentration-dependent and related to the relative potency in blocking INa. In fact, the curve describing the concentration-dependent shift of the h
curve (taken as
Vh1/2) by Benzyl-To9 was markedly shifted to the left with respect to that of To9 (Fig. 4). The fit of the data points to the equation described under Materials and Methods, and previously used to calculate drug affinity for the inactivated state (Bean et al., 1983
), led to Ki values of 27 ± 2.2 and 0.38 ± 0.04 µM for To9 and Benzyl-To9, respectively. As can be seen, these values were in the same range but did not exactly overlap the IC50 values obtained at -70 mV (Table 1); the Ki value found for Benzyl-To9 was about 1.6-fold lower, and this trend was also observed with some other compounds (data not shown). This is related to the fact that the value obtained at -70 mV is the result of the mixed binding of the drug to both resting and inactivated channel, although the high-affinity binding to the inactivated channel tends to overwhelm (Bean et al., 1983
). To validate this view, the values of K-70 were used to calculate the absolute Ki using the equation 1/K-70 = h/Kr + 1 - h/Ki, which takes into consideration the relative distribution of channels between the two states according to the membrane potential (Bean et al., 1983
). The resulting Ki values, listed in Table 1, were almost halved with respect to K-70, as expected given the contribution of the resting-state binding. However, the relative potency between the compounds was largely maintained. By this approach, the Ki of Benzyl-To9 was as low as 0.25 µM. Although the absolute values of Ki could be better estimated with more sophisticated biophysical approaches, our results indicate that the high potency shown by the compounds at depolarized membrane potentials is related to their high affinity for the inactivated state and to their ability to markedly shift the equilibrium of closed-state inactivation.
|
Use-Dependent Block of Na+ Channels by Proline-Like and Benzylated Tocainide Derivatives. Use-dependent behavior of tested compounds has been evaluated at HP of -100 mV, closer to the physiological values, with 10-ms depolarizing pulse to -20 mV applied at the frequency of 2 and 10 Hz. At both frequencies, the first trace recorded in the presence of the compound was reduced versus the control because of the tonic block exerted by the drug. When the fiber was repetitively stimulated at high frequency, the block progressively cumulated over the initial tonic block until a new equilibrium was reached. Figure 5 shows, for each compound, a sample of the progressive use-dependent block of sodium currents at 10 Hz; the concentration-response curves comparing the block before and after steady-state, use-dependent block at 10 Hz are shown just below the sample traces. The reduction of the current at the first pulse of the train was influenced only by the membrane potential and by the relative distribution of the channels between resting and inactivated states. In fact, the curves describing the tonic block and the relative calculated IC50 values were intermediate between those obtained at -140 and -70 mV (Table 2). The use-dependent block adds cumulatively to the tonic block; thus, to evaluate the absolute drug potency in real conditions of high-frequency depolarizing trains, we constructed the dose-response curves by considering the residual current after the high-frequency stimulation with respect to that in the absence of drug. The IC50 values obtained from experimental data fitting at both 2 and 10 Hz, along with the gain of potency of each compound caused by the use-dependent behavior, are listed in Table 2. As previously shown, To5 was markedly more use-dependent with respect to Toc, with gain of potency up to 8-fold at 10 Hz; surprisingly, To9 showed use-dependent behavior similar to that of Toc, producing a similar block at 2 and 10 Hz (Table 2). Each benzylated compound showed a specific pattern of use-dependent behavior. In fact Benzyl-Toc maintained the use-dependent behavior of Toc, although it was more potent in absolute terms (Fig. 6). On the other hand, the presence of the benzyl moiety on the
-proline derivative To5 dramatically reduced use-dependent behavior, the gain of potency ranging around 2-fold at 10 Hz. In contrast, Benzyl-To9, was also remarkably use-dependent, showing gain of potency values approaching those observed with To5 (Fig. 6 and Table 2). Thus Benzyl-To9, which proved to be the most potent analog, was also the one showing almost the highest use-dependence.
|
|
|
A comparison between the IC50 values listed in Tables 1 and 2 allows us to appreciate that the high-stimulation frequency lowers the effective concentration of the compounds in a manner similar to, but not overlapping, that produced by the depolarized membrane potential. In fact, use-dependent behavior is a complex dynamic process involving the kinetics of drug binding and unbinding to the channel in relation to both state-dependent drug affinity and physicochemical properties; these properties influence the ability of the drug to access the receptor site. To gain insight into this dynamic process, the time course of reaching steady state after a 10-Hz stimulation was determined for each drug at each concentration. The temporal pattern was obtained by normalizing all subsequent current traces versus the first pulse, so as to exclude tonic block, until the steady-state fractional block was reached. Figure 6, A and B, shows the pulse-dependent current reduction at 10 Hz for two extreme compounds in terms of potency of INa block: tocainide and Benzyl-To9. For both compounds, the steady state was reached within 10 s of stimulation over the 30 s of total train duration, and this time was chosen in Fig. 6 to better show graphically the different time course between the compounds. For each concentration of these two compounds, as well as for all the others, the current reduction was well fitted with a single exponential function that allowed calculation of both the time course to steady state and the fractional use-dependent block from the residual current amplitude. Toc turned out to be the fastest compound, whereas Benzyl-To9, able to produce a notable fractional use-dependent block at much lower concentrations, reaches the steady state with a remarkably slower time course. Figure 6C shows the IC50 values of all the compounds for producing the fractional steady-state, use-dependent block. Toc, Benzyl-Toc, and Benzyl-To5 produced a fractional block at 10 Hz with IC50 values of about 400, 40, and 25 µM, respectively. These values are similar to those observed for tonic block in Table 2, whereas the gain of potency for absolute use-dependent block for these compounds at 10 Hz, calculated without excluding the contribution of the tonic block (TB), was around 2. For To5, To9, and Benzyl-To9, a significantly lowering of IC50 values versus those for tonic block was instead observed, with ratios TB/fractional use-dependent block (fUDB) of 1.8 (To9), 4.6 (To5), and 5.2 (Benzyl-To9). This phenomenon (i.e., the degree of apparent change in drug potency during the high frequency stimulation) seems influenced by the ability of the compounds to equilibrate more or less rapidly with the channel in relation to the interpulse time. This point is addressed in Fig. 6D, which compares the time constants required to reach the steady-state block for the various drugs. Because of the great differences in potency between the compounds, the time constants have been compared at equieffective points (i.e., at concentrations able to produce the same steady-state block). The time constants were inversely proportional with the amount of block and consequently with the drug concentration; thus, the higher the block, the shorter the time constant. Beside this general view, marked differences were observed between the various compounds. Generally, the compounds characterized by a high use-dependent behavior showed long time constants, with values ranging also in the tenths of seconds. However, this phenomenon was not correlated with the absolute potency for either the resting or the inactivated state, because some potent compounds, such as Benzyl-Toc and Benzyl-To5, showed a faster time constant compared with less potent compounds such as To9 and To5. This result suggests the involvement of other determinants in the kinetic of drug block and in use-dependent behavior.
Interpretation of Potency and Use-Dependence in Terms of Physicochemical Properties and Conformational Features. The physicochemical properties of the compounds are shown in Table 3. The ClogP values follow the expected increase in relation to 1) the presence of the proline-like cycle and 2) the presence of the benzyl moiety, whereas the logD undergoes more complex variations, and is highly influenced by pKa (see below). No linear correlation could be found between these parameters and drug potency, although the strong potency of the benzylated analogs was mirrored by their high lipophilicity. According to our previous hypotheses, this could be related to the presence of the second aromatic ring in a crucial pharmacophore position rather than to lipophilicity by itself (De Luca et al., 2002). In fact, Benzyl-To5, which was the compound with the best "ideal" properties in terms of logP and logD for diffusion through the membrane, was the least potent of the three benzylated analogs. To gain insight in this feature, we performed a molecular modeling study aimed at evaluating the lowest conformation geometry for these compounds. As can be seen in Fig. 8, Benzyl-Toc and Benzyl-To9 assume a similar spatial conformation, with the two benzyl groups in outstretched position. On the contrary, Benzyl-To5 attains stability in a more folded conformation, with the aromatic groups almost facing each other. The calculated distance between the centroid of the aromatic rings turned out to be higher than 8 Å for both Benzyl-Toc and Benzyl To-9, whereas it was >1 Å less for Benzyl-To5. These data strongly corroborate that the higher and similar potency of the two former compounds is related to the similar distance between the two aromatic rings, which allows a similar interaction with the hydrophobic binding pockets.
|
|
A large change in estimated pKa values was also observed. In fact, according to previous data (De Luca et al., 2000
; Talon et al., 2001
), To5 showed a higher pKa than Toc, and the homologating of the amino group by one methylene further increases the basicity. As generally observed for tertiary amines compared with secondary amines in an aqueous environment (Aue et al., 1972
), Benzyl-To5 and Benzyl-To9 were less basic than the respective nonbenzylated analogs. No linear correlation was found between pKa and drug affinities for either resting or inactivated channels, so we focused on the relationship between pKa and use-dependent behavior. In fact, as mentioned in the previous paragraph, this latter is a dynamic process resulting from state-dependent affinity of the protonated molecule and physicochemical properties that influence drug access to and egress from the binding site. We evaluated the relationship between the use-dependent properties of each compound, expressed as the ratio between TB and fUDB and pKa. As can be seen in Fig. 8A, no linear correlation could be observed when considering all the compounds (r2 = 0.219). However, the deviation from linearity is mainly caused by the high pKa value of To9; in fact, when excluding this compound from the fit, the correlation rose to r2 = 0.969. This high correlation needs to be better verified, because the points are relatively few and mostly clustered around two pKa values. However, the data might suggest that use-dependence is indeed correlated in a gaussian-like manner to the pKa, with a preferential range of pKa values leading to the highest use-dependent block. Other than the absolute use-dependent potency, pKa may influence the kinetic of drug-receptor interaction. To qualitatively test this hypothesis, we evaluated the relation between the time constants to reach the steady-state 20% fractional block of INa after a train at 10 Hz and the pKa of the various compounds (Fig. 8A). We found that the cluster of compounds with pKa of 7.4 to 7.7 were those with shorter time constants (i.e., Toc, Benzyl-Toc, and Benzyl-To5). However, remarkable differences are present between the compounds of this group; the compounds being faster with the following order Toc > Benzyl-Toc > Benzyl-To5. Interestingly, the fastest Toc was also the compound with the lowest affinity for inactivated state, which is believed to contribute to the accumulation of use-dependent blockade. Benzyl-To5, slightly less potent than Benzyl Toc, has a slightly higher pKa value (7.7 versus 7.4) which can account for its longer time constant. The increase of pKa value of one unit was paralleled by a further prolongation of the time constants with respect to the previous group, suggesting that the higher quote of ionized compound required a longer time or a higher number of openings to reach the receptor site (Liu et al., 2003
). In this case, it is worth notice that Benzyl-To9 showed a time constant value that was more than 2-fold longer that of To5, probably in relation to the 40-fold higher affinity for the inactivated state of the former. This finding supports the role of a high-affinity drug-receptor interaction in contributing to the kinetic of drug binding. Interestingly, To9, which in the previous analysis was the compound deviating from linear correlation between TB/fUDB ratio and pKa, had the same Ki values as To5 but a higher pKa value, which is probably responsible for the remarkable slowing of the time constant. These results again corroborate the idea that lower or higher pKa values can be detrimental for use-dependence, preventing proper receptor interaction or proper access to the binding site.
| Discussion |
|---|
|
|
|---|
- with a
-proline cycle did not substantially change the drug potency and produced a decrease in the use-dependent behavior, probably in relation to less favorable physicochemical properties (see below). The most interesting observation was that the introduction of a benzyl group on the pharmacophore amino terminal, a structural change that has been exploited to a limited degree in
-
interactions (Dougherty, 1996
helix at an ideal distance of 11 Å (Ragsdale et al., 1994
-cation interaction with a positively charged residue. A possible candidate is the Lys1273 in the selectivity filter, which has been found to affect LA binding (Sunami et al., 1997
-
interaction (McGaughey et al., 1998
|
Based on their high-affinity, state-dependent block, the newly synthesized compounds may have a wide and rational therapeutic use. Specific residues and gating differences account for the higher affinity of cardiac versus skeletal muscle and nerve sodium channel (Wang et al., 1996
; Li et al., 2002
). Also, others' and our results have shown that this class of compounds can selectively suppress the abnormal opening of myotonia-causing sodium channel mutants, based on a combination of their state-dependent action and mutation-related gating differences (Weckbecker et al., 2000
; Desaphy et al., 2001
; Takahashi and Cannon, 2001
). Consequently, the newly synthesized compounds may offer an improvement of the therapeutic profile and selectivity with respect to tocainide.
In particular, the N-benzylated analogs have a combined increase of potency and use-dependence and can selectively address pathological signs at very low dosages. In parallel, the increase in lipophilicity may also offer advantages in pharmacokinetic properties and tissue distribution, whereas the more constrained conformation of the proline-like analogs may help the validation of the presumed map of residues involved in drug binding.
| Footnotes |
|---|
ABBREVIATIONS: LA, local anesthetic; Toc, tocainide; Benzyl-Toc, N-benzyl-tocainide; To5, N-(2,6-dimethylphenyl)pyrrolidine-2-carboxamide; Benzyl-To5, 1-benzyl-N-(2,6-dimethylphenyl)pyrrolidine-2-carboxamide; To9, N-(2,6-dimethylphenyl)pyrrolidine-3-carboxamide; Benzyl-To9, 1-benzyl-N-(2,6-dimethylphenyl)-pyrrolidine-3-carboxamide; MOPS, 3-(N-morpholino)propanesulfonic acid; HP, holding potential; fUDB, fractional use-dependent block (percentage reduction of INa at the steady-state versus the first pulse of 10-Hz train); TB, tonic block.
1 Present address: Laboratoire de Physiologie Generale, Unite Mixte de Recherche 6018, Centre National de la Recherche Scientifique, Faculte des Sciences et des Techniques, Universite de Nantes, F-44322 Nantes, France. ![]()
Address correspondence to: Prof. Diana Conte Camerino, Unità di Farmacologia, Dipartimento Farmacobiologico, Facoltà di Farmacia, University di Bari, Via Orabona 4, Campus 70121, Bari, Italy. E-mail: conte{at}farmbiol.uniba.it
| References |
|---|
|
|
|---|
Aue DH, Webb HM, and Bowers MT (1972) A thermodynamic analysis of solvation effects on the basicities of alkylamines. An electrostatic analysis of substituent effects. J Am Chem Soc 98: 318-329.[CrossRef]
Bean PB, Cohen CJ, and Tsien RW (1983) Lidocaine block of cardiac sodium channels. J Gen Physiol 81: 613-642.
Brown ML, Zha CC, Van Dyke CC, Brown GB, and Brouillette WJ (1999) Comparative molecular field analysis of hydantoin binding to the neuronal voltage-dependent sodium channel. J Med Chem 42: 1537-1545.[CrossRef][Medline]
Cannon SC (1996) Ion-channel defects and aberrant excitability in myotonia and periodic paralysis. Trends Neurosci 19: 3-10.[CrossRef][Medline]
Catterall WA (2002) Molecular mechanisms of gating and drug block of sodium channels. Novartis Found Symp 241: 206-218.[Medline]
Chambers CC, Hawkins GD, Cramer CJ, and Truhlar DG (1996) Model for aqueous solvation based on class IV atomic charges and first solvation shell effects. J Phys Chem 100: 16385-16398.[CrossRef]
Clare JJ, Tate SN, Nobbs M, and Romanos MA (2000) Voltage-gated sodium channels as therapeutic targets. Drug Discov Today 5: 506-520.[CrossRef][Medline]
Conte Camerino D, Pierno S, De Luca A, and Bryant SH (2000) Antimyotonic effects of tocainide enantiomers on skeletal muscle fibers of congenitally myotonic goats. Neuromuscul Disord 10: 160-164.[CrossRef][Medline]
De Luca A, Natuzzi F, Desaphy JF, Loni G, Lentini G, Franchini C, Tortorella V, and Conte Camerino D (2000) Molecular determinants of mexiletine structure for potent and use-dependent block of skeletal muscle sodium channels. Mol Pharmacol 57: 268-277.
De Luca A, Natuzzi F, Falcone G, Duranti A, Lentini G, Franchini C, Tortorella V, and Conte Camerino D (1997a) Inhibition of frog skeletal muscle sodium channels by newly synthesized chiral derivatives of mexiletine and tocainide. Naunyn-Schmiedeberg's Arch Pharmacol 356: 777-787.[CrossRef][Medline]
De Luca A, Pierno S, Natuzzi F, Franchini C, Duranti A, Lentini G, Tortorella V, Jockusch H, and Conte Camerino D (1997b) Evaluation of the antimyotonic activity of mexiletine and some new analogs on sodium currents of single muscle fibers and on the abnormal excitability of the myotonic ADR mouse. J Pharmacol Exp Ther 282: 93-100.
De Luca A, Pröbstle T, Brinkmeier H, and Rüdel R (1991) The different use dependences of tocainide and benzocaine are correlated with different effects on sodium channel inactivation. Naunyn-Schmiedeberg's Arch Pharmacol 344: 596-601.[Medline]
De Luca A, Talon S, De Bellis M, Desaphy JF, Franchini C, Lentini G, Catalano A, Corbo F, Tortorella V, and Conte Camerino D (2003) Inhibition of skeletal muscle sodium currents by mexiletine analogues: specific hydrophobic interactions rather than lipophilia per se account for drug therapeutic profile. Naunyn-Schmiedeberg's Arch Pharmacol 367: 318-327.[CrossRef][Medline]
Desaphy JF, De Luca A, Tortorella P, De Vito D, George AL Jr, and Conte Camerino D (2001) Gating of myotonic Na channel mutants defines the response to mexiletine and a potent derivative. Neurology 57: 1849-1857.
Dougherty DA (1996) Cation-
interactions in chemistry and biology: a new view of benzene, Phe, Tyr, and Trp. Science (Wash DC) 271: 163-168.[Abstract]
D'Souza F, Zandler ME, Deviprasad GR, and Kutner W (2000) Acid-base properties of fulleropyrrolidines: experimental and theoretical investigations. J Phys Chem A 104: 6887-6893.[CrossRef]
Franchini C, Corbo F, Lentini G, Bruno G, Scilimati A, Tortorella V, Conte Camerino D, and De Luca A (2000) Synthesis of new 2,6-prolylxylidide analogues of tocainide as stereoselective blockers of voltage-gated Na+ channels with increased potency and improved use-dependent activity. J Med Chem 43: 3792-3798.[CrossRef][Medline]
Hanson AW and Röhrl M (1972) The crystal structure of lidocaine hydrochloride monohydrate. Acta Crystallogr Sec B 28: 3567-3571.[CrossRef]
Hille B and Campbell DT (1976) An improved Vaseline gap voltage clamp for skeletal muscle fibers. J Gen Physiol 67: 265-293.
Jurkat-Rott K and Lehmann-Horn F (2001) Human muscle voltage-gated ion channels and hereditary disease. Curr Opin Pharmacol 1: 280-287.[CrossRef][Medline]
Kuo CC, Huang RC, and Lou BS (2000) Inhibition of Na+ current by diphenhydramine and other diphenyl compounds: molecular determinants of selective binding to the inactivated channels. Mol Pharmacol 57: 135-143.
Li HL, Galue A, Meadows L, and Ragsdale DS (1999) A molecular basis for the different local anesthetic affinities of resting versus open and inactivated states of the sodium channel. Mol Pharmacol 55: 134-141.
Li RA, Ennis IL, Tomaselli GF, and Marban E (2002) Structural basis of differences in isoform-specific gating and lidocaine block between cardiac and skeletal muscle sodium channels. Mol Pharmacol 61: 136-141.
Liu H, Atkins J, and Kass RS (2003) Common molecular determinants of flecainide and lidocaine block of heart Na+ channels: evidence from experiments with neutral and quaternary flecainide analogues. J Gen Physiol 121: 199-214.
McGaughey GB, Gagnés M, and Rappé AK (1998)
-stacking interactions. J Biol Chem 273: 15458-15463.
Nau C, Wang SY, Strichartz GR, and Wang GK (1999) Point mutations at N434 in D1-S6 of Mu1 Na+ channels modulate binding affinity and stereoselectivity of local anesthetic enantiomers. Mol Pharmacol 56: 404-413.
Ragsdale DS, McPhee JC, Scheuer T, and Catterall WA (1994) Molecular determinants of state-dependent block of Na+ channels by local anesthetics. Science (Wash DC) 265: 1724-1728.
Ragsdale DS, McPhee JC, Scheuer T, and Catterall WA (1996) Common molecular determinants of local anesthetic, antiarrhythmic and anticonvulsant block of voltage-gated Na+ channels. Proc Natl Acad Sci USA 93: 9270-9275.
Rüdel R, Lehmann-Horn F, and Ricker K (1994) The non-dystrophic myotonias, in Myology (Engel AG and Franzini-Armstrong C eds) pp 1291-1302, McGraw-Hill Companies, New York.
Sunami A, Dudley SC Jr, and Fozzard HA (1997) Sodium channel selectivity filter regulates antiarrhythmic drug binding. Proc Natl Acad Sci USA 94: 14126-14131.
Takahashi MP and Cannon SC (2001) Mexiletine block of disease-associated mutations in S6 segments of the human skeletal muscle Na+ channel. J Physiol (Lond) 537: 701-714.
Talon S, De Luca A, De Bellis M, Desaphy JF, Lentini G, Scilimati A, Corbo F, Franchini C, Tortorella P, Jockusch H, et al. (2001) Increased rigidity of the chiral centre of tocainide favours stereoselectivity and use-dependent block of skeletal muscle Na+ channels enhancing the antimyotonic activity in vivo. Br J Pharmacol 134: 1523-1531.[CrossRef][Medline]
Wang DW, Nie L, George AL Jr, and Bennett PB (1996) Distinct local anesthetic affinities in Na+ channel subtypes. Biophys J 70: 1700-1708.[Medline]
Wang SY, Nau C, and Wang GK (2000) Residues in Na+ channel D3-S6 segment modulate both batrachotoxin and local anesthetic affinities. Biophys J 79: 1379-1387.[Medline]
Weckbecker K, Wurz A, Mohammadi B, Mansuroglu T, George AL Jr, Lerche H, Dengler R, Lehmann-Horn F, and Mitrovic N (2000) Different effects of mexiletine on two mutant sodium channels causing paramyotonia congenita and hyperkalemic periodic paralysis. Neuromuscul Disord 10: 31-39.[CrossRef][Medline]
Wright SN, Wang SY, and Wang GK (1998) Lysine point mutations in Na+ channel D4-S6 reduce inactivated channel block by local anesthetics. Mol Pharmacol 54: 733-739.
Yang Y-C and Kuo C-C (2002) Inhibition of Na+ current by imipramine and related compounds: different binding kinetics as an inactivation stabilizer and as open channel blocker. Mol Pharmacol 62: 1228-1237.
Yarov-Yarovoy V, Brown J, Sharp EM, Clare JJ, Scheuer T, and Catterall WA (2001) Molecular determinants of voltage-dependent gating and binding of pore-blocking drugs in transmembrane segment IIIS6 of the Na+ channel
subunit. J Biol Chem 276: 20-27.
Yarov-Yarovoy V, McPhee JC, Idsvoog D, Pate C, Scheuer T, and Catterall WA (2002) Role of amino acid residues in transmembrane segments IS6 and IIS6 of the Na+ channel
subunit in voltage-dependent gating and drug block. J Biol Chem 277: 35393-35401.
This article has been cited by other articles:
![]() |
P. J. Jones, Y. Wang, M. D. Smith, N. J. Hargus, H. S. Eidam, H. S. White, J. Kapur, M. L. Brown, and M. K. Patel Hydroxyamide Analogs of Propofol Exhibit State-Dependent Block of Sodium Channels in Hippocampal Neurons: Implications for Anticonvulsant Activity J. Pharmacol. Exp. Ther., February 1, 2007; 320(2): 828 - 836. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. M. Lipkind and H. A. Fozzard Molecular Modeling of Local Anesthetic Drug Binding by Voltage-Gated Sodium Channels Mol. Pharmacol., December 1, 2005; 68(6): 1611 - 1622. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||