| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Department of Physiology, Michigan State University, East Lansing, Michigan 48824
Correspondence: Address reprint requests to Patrick F. Dillon, Tel.: 517-355-6475, ext. 1102; Fax: 517-355-5125; E-mail: dillon{at}msu.edu.
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
In testing these questions we ascertained the effect of electric fields on the small molecule-small molecule complex NE-morphine sulfate (MS), the small molecule-protein complex epinephrine-bovine serum albumin (Epi-BSA), and the protein-protein complex insulin-glucagon (IG). The known binding of NE and MS (2
) may play a role in both the therapeutic (3
) and abusive (4
) effects of morphine. The known binding of Epi to albumin (5
) protects Epi from oxidation and enzymatic degradation as it passes through the circulatory system. Yet, there still must be dissociation of the complex for Epi to have its hormonal effects. The data show that for the measured dissociation constant of NE-Asc binding, a significant fraction of NE will circulate bound to Asc (1
). In contrast, it has been shown that insulin and glucagon can bind to one another (6
) but with a dissociation constant determined here that is well above both their physiological concentrations. Thus, although this is a useful model system for measuring protein-protein binding, it can also preclude the binding of insulin and glucagon under physiological conditions.
The positive results of the experiments using these complexes allowed us to make several unique conclusions. These include the quantitative expression of the molecular shielding constant, the degree to which a molecular complex shields its binding sites from an electric field based on the association energy of the complex and the molecular radius of the complex components. Also, having ascertained the behavior of different sized molecules in the electric field, we were able to conclude that complexes of very large proteins have such a large amount of shielding and limited approach to the membrane that they are not substantially affected by the membrane electric field. The consequences of our results are discussed below.
| METHODS |
|---|
|
|
|---|
Solutions
All solutions were prepared in the buffer which was used in the particular capillary electrophoresis (CE) experiment, 25 mM sodium borate 10-hydrate (Na2B4O7·10 H2O) at pH 9.4. NE, Asc, MS, Epi, BSA, insulin, and glucagon were purchased from Sigma-Aldrich Chem. Co. (St Louis, MO). The BSA is Sigma A7409, Initial Fractionation by Heat Shock, 35% in 0.85% sodium chloride. All chemicals were made as stock solutions and diluted to the required concentration in the CE injection vials. The NE-Asc dissociation constants were determined using 1 mM NE and altering the concentration of Asc between 0 and 100 mM. The NE-MS dissociation constants were determined using 1 mM NE at 204 and 255 V/cm and 0.1 mM at 153 V/cm and altering the concentration of MS between 0 and 40 mM. The IG dissociation constants were determined using a constant concentration of 20 µM glucagon and altering the concentration of insulin between 0 and 1.6 mM. The reciprocity of the IG binding was demonstrated using 20 µM insulin and 0, 19, and 190 µM glucagon. The Epi-BSA dissociation constants were determined using a constant concentration of 40 µM Epi and altering the concentration of BSA between 0 and 2.32 mM.
CE procedures
Samples were vacuum injected into a 100-µm diameter, 98-cm length capillary tube (volume 7.7 µl) for 2 s (injection volume 8.6 nl/s). The capillary tube had a detection window 66 cm from the injection site in an ISCO (Lincoln, NE) model 3850 electropherograph with an absorbance detector at 195 nm, absorbance maximum from 0.020.2 absorbance units, and a rise time of 3.2 s. The driving voltage was varied between 5 and 25 kV, corresponding to electric fields of 51 and 255 V/cm. The peaks were recorded on a chart recorder (The Recorder Company, Houston, TX) at 1 cm/min and 1 V full scale.
Molecular radius
Estimates of molecular radius are based on the model of the NE-Asc complex (7
) and the CPK molecular model system (8
) of 1.25 cm/Å. For the 345 g/mol NE-Asc complex, the complex volume is estimated to be 0.63 nm3. The relationship between molecular weight (MW) and volume is therefore 0.00183 (MW) = nm3. For a sphere of volume 4
r3/3, this corresponds to a radius of 0.53 nm for the NE-Asc complex. The radii for the other molecules and complexes are calculated in the same way.
Data analysis
The constant molecule (NE, Epi, or glucagon) was always prepared and run as a standard. For each complex, the initial electric field dependent dissociation constant was calculated by plotting the log of the varied molecule concentration against the log of the constant-varied complex/free constant ratio and applying a least-squares best fit for the linearization and 95% confidence interval using the Axum data processing system. The initial dissociation constant at a particular electric field, the KE, was estimated from the linearization at the 0.5 complex-0.5 free point. The initial KE = [free constant][free varied]/[constant-varied complex]. The [free constant] equals the [constant-varied complex] at the KE. Thus, the true KE is the [free varied], which is equal to the [total varied] 0.5[total constant], and this was the value used for the KE. The estimate of the dissociation constant at zero electric field, the KD, was made by extrapolating the log of the KE values at a given electric field to zero electric field. The binding energy of the complexes was calculated by multiplying the log(KD) by 2.303RT for the determination of J/mol at 298 K.
Membrane dissociation
The membrane potential of a membrane decreases exponentially with a Debye-Hückel length constant of 1 nm and a voltage change of 30 mV (9
11
). The electric field is the length derivative of the voltage decay (10
,11
). Having ascertained the relationship of slope (
) (converting V/cm into mV/nm) and KD between electric field and dissociation constants for a complex pair, the dissociation constant at a given distance (x) from a membrane with a change in voltage (V0) and a length constant (
) (1
) is
![]() |
This equation can be rearranged to calculate the distance at which different changes in KD occur. This is done in Table 1 for 10-, 100-, and 1000-fold changes in the dissociation constant as the four different complexes approach the membrane.
|
| RESULTS |
|---|
|
|
|---|
|
|
|
|
-
bond (7
5 kJ/mol (12
|
As the change in dissociation produced by the membrane requires that the complex approach the membrane, very large molecular complexes will not be able to approach the membrane sufficiently closely to increase dissociation significantly. Fig. 6 shows the log of the complex radius plotted against the log of the distance necessary for a 10-fold decrease in the dissociation constant. At the x intercept, the radius, calculated from the log(radius), will be the largest size a complex can be and still have a 10-fold change in dissociation. This size corresponds approximately to an 800,000 MW complex.
|
|
|
| DISCUSSION |
|---|
|
|
|---|
The multiple examples of molecular coupling show that the CE method of measuring dissociation constants is general. When holding one molecule concentration constant, changes in the concentration of its complementary molecule result in peak changes to the first molecule. CE causes molecules to separate in an applied electric field. If two molecules are bound at zero electric field, and separated at high electric field, there must be an electric field at which they are half-bound and half-free. Further, as the electric field increases, increasing the force of separation, then the concentration of complementary molecule needed to maintain half of the first molecule in the complexed form must also increase. Fig. 5 shows that this is the case for all these combinations of molecules tested here. Extrapolating to zero electric field defines the dissociation constant, KD, in free solution.
The insulin-glucagon and NE-MS electropherograms in Figs. 2 and 3 show that peak changes can be reciprocal. Just as insulin causes the glucagon peak to decrease, so too does glucagon cause the insulin peak to increase. In the case of NE-MS, however, just as the NE peak decreases with an increase in MS, the rightmost peak of the MS complex decreases in the presence of NE. Thus, one cannot make a priori judgments on the direction a peak change will take when its complement is present. Molecular coupling would normally be expected to produce changes in the peak characteristics, although there could be circumstances in which this does not occur.
We have shown that NE and Asc are connected by four hydrogen bonds and one
-
interaction (7
). Hydrogen bonds and
-
bonds have similar energy (13
). In aqueous solutions, hydrogen bonds must compete with water molecules, and hydrogen bonds between nonwater molecules have a strength of 5 kJ/mol (12
). The binding energy we calculate for the NE-Asc binding of 25.2 kJ/mol is consistent with the five interactions previously shown (7
). The IG and Epi-BSA bindings may have a different number or type of bonds but produce a similar binding energy. The NE-MS association energy is consistent with having stronger bonds than the NE-Asc binding. NE-MS binding involves a stronger
-
bond than NE-Asc, hydrogen bonds and an ionic bond (2
), resulting in a stronger association energy.
Even a cursory analysis of Fig. 5 shows that those complexes containing a protein have much shallower electric field dependence than those with only small molecules. Since increasing electric field strength increases dissociation, proteins must reduce the amount of dissociation as the electric field increases. This cannot be due to association energy differences, as the protein association energies are in the same range as the NE-Asc small molecule pair association energy. Therefore, qualitatively, increasing the size of molecules must protect their binding sites from the effects of the electric field. The molecular shielding constant provides a starting point for estimates of the effect of an electric field on a complex. Since the radii for most molecules can be estimated fairly well, the slope/association energy can be predicted. The units of the molecular shielding constant, cm2/V, are also intuitively satisfying. An applied voltage is spread over an area, and the greater the area the more the voltage is distributed, lessening its effects. Larger molecules will have greater areas, lessening the dissociative effects of the electric field.
Given the assumptions made regarding the spherical shapes of the different molecules and their subunits (insulin is a dimer (14
) and glucagon a trimer (15
)), it is remarkable that this model of the data, using only the radii and association energies, should give such a good correlation. The radii used in these calculations are based solely on the MW of the molecules, assuming a spherical shape. No estimate of the nonspherical shape of particular molecules or of molecular complexes was included. That using the simplest assumptions produced such robust results leads us to conclude that although nonspherical shapes may alter the quantitation of the molecular shielding constant, this alteration will be small, as will changes in the estimates of the dissociation distance as a complex approaches the membrane. The spherical assumption may work, in part, because molecules are constantly rotating in the electric field. Calculations of the effect of the electric field on the rotational energy of a molecule show that this effect is minimal, with thermal forces causing random orientations far exceeding any orientation of the complex produced by the field (16
). The change in dipole orientation equilibrium will be <1% at 104 V/cm. Thus, any quantitative effects produced by the nonspherical nature of a molecule or complex will be minimized by the randomization of its orientation within the electric field.
Because the membrane potential, in millivolts, exponentially decays over a distance of nanometers in the surrounding fluid, the electric field near the membrane, in mV/nm, or 104 V/cm, will be much greater than the field needed to produce dissociation, which is in the 102 V/cm range (see Fig. 5). The exponential decay of the field means this field will only exist within 10 nm or less from the membrane. But within this range, as shown in Table 1, the membrane electric field is more than sufficient to produce dissociation of molecular complexes. Very large molecular complexes of proteins may be unable to approach the membrane sufficiently close to cause dissociation, as shown in Fig. 6. Large protein complexes will have both molecular shielding and limited membrane approach, minimizing the effect of the electric field on these complexes. This will not be the case for small molecule complexes.
The distance over which the membrane electric field could affect agonist binding covers the range of membrane protein projection into the interstitial fluid. The electric field (E) is the length derivative of the voltage. As previously cited, the dielectric constant of water produces a Debye-Hückel length constant (
) of 1 nm (9
11
) at physiological ionic strength. The equation governing the determination of the electric field is
![]() |
This phenomenon will not be limited to the interstitial space. The electric field will also be significant for several nanometers on the intracellular side of the membrane, where the electric field will mirror the electric field on the outside of the membrane (20
). So too, molecular complexes circulating in the plasma, such as Epi-BSA, will dissociate as they near the capillary membrane. Other factors not addressed here experimentally should also be considered. The glycocalyx extends the membrane potential by several nanometers on the outside (but not the inside) of the cell membrane (21
). The membrane potential has its greatest change (and largest electric field) just beyond the glycocalyx. Thus, in cells with a large glycocalyx, the effect of the electric field on agonist binding will have this additional factor. This will be a further consideration when considering the degree to which the water shell around an ion such as Na+ is stripped by the electric field as it approaches an ion channel. If the extent of the pore is within the range of the glycocalyx, the stripping effect of the electric field will be reduced.
It may be possible to apply electric dissociation in a specific manner for practical purposes such as drug delivery. Bonds with energies in this range absorb energy in the infrared spectrum (13
). Capillary electrophoresis uses DC voltage to effect dissociation. Within the DC field are all the AC frequencies with wavelengths shorter than the length of the DC field, in this case 98 cm. The power input in our system for the 204 V/cm field is 1.58 J/s. Application of one or several frequencies specific to the complementary pair will split the pair with a lower total energy input. Differences in the number of H-bonds,
-
bonds, and frequency differences produced by bond stretching (13
) may produce unique frequency patterns that allow selective dissociation of a specific complex, if these frequencies can be determined and applied. Thus, it should be possible to separate molecular complexes with much lower energy input than using conventional CE. If these frequencies can penetrate skin sufficiently, they can be used to dissociate complexes in situ.
In summary, capillary electrophoresis can be used to measure the dissociation constants of many complementary pairs; the large size of proteins will reduce the effect of the electric field on the dissociation of bound molecules; membrane electric fields will cause dissociation of all complexes except the largest protein-protein complexes within 10 nm of the membrane; and agonists can be delivered to membrane receptor sites bound to carrier molecules and dissociate in the vicinity of the receptor.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Submitted on August 3, 2005; accepted for publication October 27, 2005.
| REFERENCES |
|---|
|
|
|---|
2. Root-Bernstein, R. S. 1987. Catecholamine bind to enkephalins, morphiceptin, and morphine. Brain Res. Bull. 18:509532.[CrossRef][Medline]
3. Motsch, J., and M. Kamler. 1997. Alpha 2-adrenergic agonists. Use in chronic paina meta analysis. Schmerz. 11:339344.[CrossRef][Medline]
4. Georges, F., and G. Aston-Jones. 2003. Prolonged activation of mesolimbic dopaminergic neurons by morphine withdrawal following clonidine: participation of imidazoline and norepinephrine receptors. Neuropsychopharmacology. 28:11401149.[Medline]
5. Patil, N., P. Fraundorfer, and P. K. Dutta. 1996. Stereo selective modification of circular dichroism spectra of rat lung ß-adrenoceptor protein preparation by enantiomers of epinephrine. Chirality. 8:43634365.
6. Root-Bernstein, R. S., and C. Dobbelstein. 2001. Insulin binds to glucagon forming a complex relationship. Autoimmunity. 33:153169.[Medline]
7. Root-Bernstein, R. S., and P. F. Dillon. 2002. Fostering venture research: a case study of the discovery that ascorbate enhances adrenergic drug activity. Drug Dev. Res. 57:5874.[CrossRef]
8. Harte, R. A. 1969. Molecules in Three Dimensions: A Guide to the Construction of Models of Biochemically Interesting Compounds with CPK Models. American Society of Biological Chemists, Bethesda, MD.
9. Delahay, P. 1966. Double Layer and Electrode Kinetics. John Wiley and Sons, New York.
10. Moore, W. J. 1972. Physical Chemistry, 4th ed. Prentice-Hall, Englewood Cliffs, NJ.
11. Tien, H. T. 1974. Bilayer Lipid Membranes (BLM): Theory and Practice. Marcel Dekker, New York.
12. Tinoco, I., K. Sauer, and J. C. Wang. 1995. Physical Chemistry: Principles and Applications in Biological Sciences. Prentice Hall, Upper Saddle River, NJ.
13. Chang, R. 1977. Physical Chemistry with Application to Biological Systems. MacMillan, New York.
14. Falconi, M., M. T. Cambria, A. Cambria, and A. Desideri. 2001. Structure and stability of the insulin dimer investigated by molecular dynamics simulation. J. Biomol. Struct. Dyn. 18:761772.[Medline]
15. Oomen, R. P., and H. Kaplan. 1990. Binding of glucagon to lipid bilayers. Biochem. Cell Biol. 68:284291.[Medline]
16. Dörr, F. 1983. Structure determination of biomolecules by physical methods. In Biophysics. W. Hoppe, W. Lohmann, H. Markl, and H. Ziegler, editors. Springer Verlag, Berlin. 168169.
17. Faller, M., M. Niederweis, and G. E. Schulz. 2004. The structure of mycobacterial outer membrane channel. Science. 303:11891192.
18. Fuchs, H., U. Lüchen, R. Tauben, A. Engel, and R. Gessner. 1998. Structure model of phospholipids-reconstitutes human transferring derived by electron microscopy. Structure. 6:12351243.[Medline]
19. Woldin, C. N., F. S. Hing, L. Jongsoon, P. F. Pilch, and G. G. Shipley. 1999. Structural studies of the detergent-solubilized and vesicle-reconstituted insulin receptor. J. Biol. Chem. 274:3498134992.
20. McLaughlin, S., and H. Harahy. 1974. Phospholipid flip-flop and the distribution of surface charges in excitable membranes. Biophys. J. 14:200208.
21. Glaser, R. 2001. Biophysics. Springer-Verlag, Berlin.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |