1
|
Levine HA, McGee MP, Serna S. Diffusion and reaction in the cell glycocalyx and the extracellular matrix. J Math Biol 2009; 60:1-26. [PMID: 19274464 DOI: 10.1007/s00285-009-0254-y] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/03/2008] [Revised: 01/27/2009] [Indexed: 11/28/2022]
Abstract
Many biologically important macromolecular reactions are assembled and catalyzed at the cell lipid-surface and thus, the extracellular matrix and the glycocalyx layer mediate transfer and exchange of reactants and products between the flowing blood and the catalytic lipid-surface. This paper presents a mathematical model of reaction-diffusion equations that simply describes the transfer process and explores its influence on surface reactivity for a prototypical pathway, the tissue factor (Tf) pathway of blood coagulation. The progressively increasing friction offered by the matrix and glycocalyx to reactants and to the product (coagulation factors X, VIIa and Xa) approaching the reactive surface is simulated and tested by solving the equations numerically with both, monotonically decreasing and constant diffusion profiles. Numerical results show that compared to isotropic transfer media, the anisotropic structure of the matrix and glycocalyx sharply decreases overall reaction rates and significantly increases the mean transit time of reactants; this implies that the anisotropy modifies the distribution of reactants. Results also show that the diffusional transfer, whether isotropic or anisotropic, influences reaction rates according to the order at which the reactants arrive at the boundary. Faster rates are observed when at least one of the reactants is homogeneously distributed before the other arrives at the boundary than when both reactants transfer simultaneously from the boundary.
Collapse
Affiliation(s)
- Howard A Levine
- Department of Mathematics, Iowa State University, Ames, IA 50011, USA.
| | | | | |
Collapse
|
2
|
Affiliation(s)
- Soohyung Park
- Institute of Chemistry and the Fritz Haber Research Center, The Hebrew University, Jerusalem 91904, Israel
| | - Noam Agmon
- Institute of Chemistry and the Fritz Haber Research Center, The Hebrew University, Jerusalem 91904, Israel
| |
Collapse
|
3
|
Park S, Agmon N. Theory and Simulation of Diffusion-Controlled Michaelis−Menten Kinetics for a Static Enzyme in Solution. J Phys Chem B 2008; 112:5977-87. [DOI: 10.1021/jp075941d] [Citation(s) in RCA: 38] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Soohyung Park
- Institute of Chemistry and the Fritz Haber Research Center, The Hebrew University, Jerusalem 91904, Israel
| | - Noam Agmon
- Institute of Chemistry and the Fritz Haber Research Center, The Hebrew University, Jerusalem 91904, Israel
| |
Collapse
|
4
|
Chiam KH, Tan CM, Bhargava V, Rajagopal G. Hybrid simulations of stochastic reaction-diffusion processes for modeling intracellular signaling pathways. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2006; 74:051910. [PMID: 17279942 DOI: 10.1103/physreve.74.051910] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/16/2006] [Indexed: 05/13/2023]
Abstract
In the intracellular environment, signaling takes place in a nonideal environment that is spatially heterogeneous and that is noisy, with the noise arising from the low copy numbers of the signaling molecules involved. In this paper, we model intracellular signaling pathways as stochastic reaction-diffusion processes and adapt techniques commonly used by physicists to solve for the spatiotemporal evolution of the signaling pathways. We then apply it to study two problems of relevance to the modeling of intracellular signaling pathways. First, we show that, in the limit of small protein diffusion which is typically the case for proteins in the cytosol crowded by other macromolecules, the extent of diffusion control, in the transient regime, on reactions is greater than previous predictions. Second, we show that the presence of scaffold proteins can modify the phosphorylation activity of a mitogen-activated protein kinase cascade, and explain how this activity is modulated by the scaffold protein concentration.
Collapse
Affiliation(s)
- K-H Chiam
- Bioinformatics Institute, 30 Biopolis Street, Singapore 138671, Singapore.
| | | | | | | |
Collapse
|
5
|
Tzafriri AR, Edelman ER. On the validity of the quasi-steady state approximation of bimolecular reactions in solution. J Theor Biol 2005; 233:343-50. [PMID: 15652144 DOI: 10.1016/j.jtbi.2004.10.013] [Citation(s) in RCA: 18] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2004] [Revised: 10/07/2004] [Accepted: 10/08/2004] [Indexed: 10/26/2022]
Abstract
Two-step binding kinetics are extensively used to study the relative importance of diffusion in biochemical reactions. Classical analysis of this problem assumes ad hoc that the encounter complex is at quasi-steady state (QSS). Using scaling arguments we derive a criterion for the validity of this assumption in the limit of irreversible product formation. We find that the QSS approximation (QSSA) of two-step binding is only valid if the total ligand and receptor concentrations are much smaller than (k2+k-1)/k1, where k1 and k-1 are, respectively, the forward and reverse diffusion encounter rate constants and k2 is the chemical association rate constant. This criterion can be shown to imply that the average time between encounters is much longer than the half-life of the encounter complex and also guarantees that the concentration of the encounter complex is negligible compared to the reactant and product concentrations. Numerical examples of irreversible and reversible cases corroborate our analysis and illustrate that the QSS may be invalid even if k-2<<k2. Our analysis of the irreversible case is shown to carry through to the more rigorous framework of the Smoluchowski theory of diffusion-controlled reactions. This work underscores the need for exercising greater caution in invoking the QSSA.
Collapse
Affiliation(s)
- A Rami Tzafriri
- Harvard-MIT Division of Health Sciences and Technology, Massachusetts Institute of Technology, Room 16-343, Cambridge, MA 02139, USA.
| | | |
Collapse
|
6
|
Høyrup P, Callisen TH, Jensen MØ, Halperin A, Mouritsen OG. Lipid protrusions, membrane softness, and enzymatic activity. Phys Chem Chem Phys 2004. [DOI: 10.1039/b314146b] [Citation(s) in RCA: 31] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
7
|
Bednarczyk D, Ekins S, Wikel JH, Wright SH. Influence of molecular structure on substrate binding to the human organic cation transporter, hOCT1. Mol Pharmacol 2003; 63:489-98. [PMID: 12606755 DOI: 10.1124/mol.63.3.489] [Citation(s) in RCA: 84] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022] Open
Abstract
Organic cation transporters play a critical role in the elimination of therapeutic compounds in the liver and the kidney. We used computational quantitative structure activity approaches to predict molecular features that influence interaction with the human ortholog of the organic cation transporter (hOCT1). [(3)H]tetraethylammonium uptake in HeLa cells stably expressing hOCT1 was inhibited to varying extents by a diverse set of 30 molecules. A subset of 22 of these was used to produce, using Catalyst, a pharmacophore that consisted of three hydrophobic features and a positive ionizable feature. The correlation coefficient of observed versus predicted IC(50) was 0.86 for this training set, which was superior to calculated logP alone (r = 0.73) as a predictor of hOCT1 inhibition. A descriptor-based quantitative structure-activity relationship study using Cerius(2) resulted in an equation relating five molecular descriptors to log IC(50) with a correlation coefficient of 0.95. Furthermore, a group of phenylpyridinium and quinolinium compounds were used to investigate the spatial limitations of the hOCT1 binding site. The affinity for hOCT was higher for 4-phenylpyridiniums > 3-phenylpyridiniums > quinolinium, indicating that substrate affinity was influenced by the distribution of hydrophobic mass. In addition, supraplanar hydrophobic mass was found to increase the affinity for binding hOCT1. These results indicate how a combination of computational and in vitro approaches may yield insight into the binding affinity of transporters and may be applicable to predicting these properties for new therapeutics.
Collapse
Affiliation(s)
- Dallas Bednarczyk
- Department of Physiology, University of Arizona, Tucson, Arizona, USA
| | | | | | | |
Collapse
|
8
|
Vissing T, Ipsen JH, Callisen TH. Vesicles of one monoglyceride and two phospholipids: phase behavior and susceptibility to hydrolysis by glyceride lipase and phospholipase A2. Colloids Surf B Biointerfaces 2002. [DOI: 10.1016/s0927-7765(01)00313-7] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
9
|
Nelsestuen GL, Stone M, Martinez MB, Harvey SB, Foster D, Kisiel W. Elevated function of blood clotting factor VIIa mutants that have enhanced affinity for membranes. Behavior in a diffusion-limited reaction. J Biol Chem 2001; 276:39825-31. [PMID: 11517221 DOI: 10.1074/jbc.m104896200] [Citation(s) in RCA: 21] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022] Open
Abstract
Blood clotting factor VIIa is involved in the first step of the blood coagulation cascade, as a membrane-associated enzyme in complex with tissue factor (TF). Factor VIIa is also an important therapeutic agent for hemophilia where its function may include TF-independent as well as TF-dependent mechanisms. This study compared the activity of wild type factor VIIa (WT-VIIa) with that of a mutant with elevated affinity for membrane (P10Q/Q32E, QE-VIIa). Phospholipid and cell-based assays showed the mutant to have up to 40-fold higher function than WT-VIIa in both TF-dependent and TF-independent reactions. Tissue factor-dependent reactions displayed the maximum enhancement when binding had reached equilibrium in competition with another TF-binding protein. In liposome-based assays, the association rate of WT-VIIa with TF occurred at a physical maximum and could not be improved by site-directed mutagenesis. A practical consequence was identical function of WT-VIIa and QE-VIIa in assays that depended entirely on assembly kinetics. Thus, factor VIIa mutants provided unique reagents for probing the mechanism of factor VIIa action. They may also offer superior agents for therapy.
Collapse
Affiliation(s)
- G L Nelsestuen
- Department of Biochemistry, Molecular Biology and Biophysics, University of Minnesota, 1479 Gortner Ave., St. Paul, MN 55108, USA.
| | | | | | | | | | | |
Collapse
|
10
|
Boskovic DS, Krishnaswamy S. Exosite binding tethers the macromolecular substrate to the prothrombinase complex and directs cleavage at two spatially distinct sites. J Biol Chem 2000; 275:38561-70. [PMID: 10984491 DOI: 10.1074/jbc.m006637200] [Citation(s) in RCA: 39] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022] Open
Abstract
The prothrombinase complex, composed of the proteinase, factor Xa, bound to factor Va on membranes, catalyzes thrombin formation by the specific and ordered proteolysis of prothrombin at Arg(323)-Ile(324), followed by cleavage at Arg(274)-Thr(275). We have used a fluorescent derivative of meizothrombin des fragment 1 (mIIaDeltaF1) as a substrate analog to assess the mechanism of substrate recognition in the second half-reaction of bovine prothrombin activation. Cleavage of mIIaDeltaF1 exhibits pseudo-first order kinetics regardless of the substrate concentration relative to K(m). This phenomenon arises from competitive product inhibition by thrombin, which binds to prothrombinase with exactly the same affinity as mIIaDeltaF1. As thrombin is known to bind to an exosite on prothrombinase, initial interactions at an exosite likely play a role in the enzyme-substrate interaction. Occupation of the active site of prothrombinase by a reversible inhibitor does not exclude the binding of mIIaDeltaF1 to the enzyme. Specific recognition of mIIaDeltaF1 is achieved through an initial bimolecular reaction with an enzymic exosite, followed by an active site docking step in an intramolecular reaction prior to bond cleavage. By alternate substrate studies, we have resolved the contributions of the individual binding steps to substrate affinity and catalysis. This pathway for substrate binding is identical to that previously determined with a substrate analog for the first half-reaction of prothrombin activation. We show that differences in the observed kinetic constants for the two cleavage reactions arise entirely from differences in the inferred equilibrium constant for the intramolecular binding step that permits elements surrounding the scissile bond to dock at the active site of prothrombinase. Therefore, substrate specificity is achieved by binding interactions with an enzymic exosite that tethers the protein substrate to prothrombinase and directs cleavage at two spatially distinct scissile bonds.
Collapse
Affiliation(s)
- D S Boskovic
- Joseph Stokes Research Institute, Children's Hospital of Philadelphia Philadelphia, Pennsylvania 19104, USA
| | | |
Collapse
|
11
|
Nielsen LK, Risbo J, Callisen TH, Bjørnholm T. Lag-burst kinetics in phospholipase A(2) hydrolysis of DPPC bilayers visualized by atomic force microscopy. BIOCHIMICA ET BIOPHYSICA ACTA 1999; 1420:266-71. [PMID: 10446309 DOI: 10.1016/s0005-2736(99)00103-0] [Citation(s) in RCA: 93] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/16/2022]
Abstract
The lag-burst phenomenon in the phospholipase A(2) mediated hydrolysis of phospholipid bilayers is for the first time demonstrated in an atomic force microscopy (AFM) study. Simultaneous AFM measurements of the degree of bilayer degradation and the physical-chemical state of the membrane reveals growing nanoscale indentations in the membrane during the lag phase. It is argued that these indentations are domains of hydrolysis products (lysoPC/PC) which eventually trigger the burst. The rate of the rapid hydrolysis following the burst is found to be proportional to the length of the edge between membrane adsorbed and desorbed to the mica base. The observed maximal rate of membrane degradation is approx. 0.2 mmol lipid/min/mol lipase in solution.
Collapse
Affiliation(s)
- L K Nielsen
- Department of Chemistry, Technical University of Denmark, DK-2800, Lyngby, Denmark
| | | | | | | |
Collapse
|
12
|
Affiliation(s)
- S Ekins
- Department of Drug Disposition, Lilly Research Laboratories, Eli Lilly and Co., Lilly Corporate Center, Indianapolis, Indiana 46285, USA
| | | |
Collapse
|
13
|
Nelsestuen GL. Significance of reduced dimensionality in reaction kinetics: impact of multi-site particles. Chem Phys Lipids 1999; 101:37-44. [PMID: 10810923 DOI: 10.1016/s0009-3084(99)00053-5] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
Abstract
This review examines novel kinetic properties of enzymes on membrane surfaces or states of restricted diffusion. A leading feature is the presence of multiple enzymes and/or substrates per particle. In these states, enzymes can be influenced by parameters such as the number of substrates or enzymes per particle, particle size, the rates of exchange of substrate or enzyme from the particle, or substrate diffusion to the particle. These steps are independent of the enzyme site parameters which are described by classical enzymology. The results make it clear that non-classical behaviors are important to biological systems, are the basis for some enzyme expression levels and are determinants of cellular design. To identify more unique functions of these states, descriptions of catalysis in the non-solution state should become a part of kinetic education in biology.
Collapse
Affiliation(s)
- G L Nelsestuen
- Department of Biochemistry, University of Minnesota, St. Paul 55108, USA.
| |
Collapse
|
14
|
Martinez MB, Flickinger MC, Nelsestuen GL. Steady-state enzyme kinetics in the Escherichia coli periplasm: a model of a whole cell biocatalyst. J Biotechnol 1999; 71:59-66. [PMID: 10483101 DOI: 10.1016/s0168-1656(99)00015-2] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
Abstract
This study provided analysis of in vivo enzyme kinetics in a model system which consisted of alkaline phosphatase in the periplasm of Escherichia coli. Modeling of complete substrate titration curves was achieved for a wide range of intraperiplasmic enzyme levels and outer membrane permeabilities. The results helped to identify the features most important to optimize in vivo reaction velocity. For many situations, a surprising finding was that maximum enzyme expression was not a major concern. For example, for moderate enzyme expression levels and moderate substrate levels (ca 0-5 mM), the limiting step for the enzyme in the periplasm was substrate (para-nitrophenylphosphate) diffusion through the outer membrane. In vivo reaction velocity was directly proportional to substrate concentration, outer membrane permeability, and the cell concentration. Velocity was also quite insensitive to a potent inhibitor of the enzyme. Even though diffusion-limited, periplasmic reaction velocity was quite sensitive to temperature, suggesting that the conformation of porin proteins in the E. coli outer membrane governed the average size of the pore. This model system therefore defined important features of bacterial whole cell biocatalyst design, which may also apply to other reactors using intact cells as catalysts.
Collapse
Affiliation(s)
- M B Martinez
- Department of Biochemistry, University of Minnesota, St. Paul 55108, USA
| | | | | |
Collapse
|
15
|
Sinniger V, Merton RE, Fabregas P, Felez J, Longstaff C. Regulation of tissue plasminogen activator activity by cells. Domains responsible for binding and mechanism of stimulation. J Biol Chem 1999; 274:12414-22. [PMID: 10212215 DOI: 10.1074/jbc.274.18.12414] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022] Open
Abstract
A number of cell types have previously been shown to bind tissue plasminogen activator (tPA), which in some cases can remain active on the cell surface resulting in enhanced plasminogen activation kinetics. We have investigated several cultured cell lines, U937, THP1, K562, Molt4, and Nalm6 and shown that they bind both tPA and plasminogen and are able to act as promoters of plasminogen activation in kinetic assays. To understand what structural features of tPA are involved in cell surface interactions, we performed kinetic assays with a range of tPA domain deletion mutants consisting of full-length glycosylated and nonglycosylated tPA (F-G-K1-K2-P), DeltaFtPA (G-K1-K2-P), K2-P tPA (BM 06.022 or Reteplase), and protease domain (P). Deletion variants were made in Escherichia coli and were nonglycosylated. Plasminogen activation rates were compared with and without cells, over a range of cell densities at physiological tPA concentrations, and produced maximum levels of stimulation up to 80-fold with full-length, glycosylated tPA. Stimulation for nonglycosylated full-length tPA dropped to 45-60% of this value. Loss of N-terminal domains as in DeltaFtPA and K2P resulted in a further loss of stimulation to 15-30% of the full-length glycosylated value. The protease domain alone was stimulated at very low levels of up to 2-fold. Thus, a number of different sites are involved in cell interactions especially within finger and kringle domains, which is similar to the regulation of tPA activity by fibrin. A model was developed to explain the mechanism of stimulation and compared with actual data collected with varying cell, plasminogen, or tPA concentrations and different tPA variants. Experimental data and model predictions were generally in good agreement and suggest that stimulation is well explained by the concentration of reactants by cells.
Collapse
Affiliation(s)
- V Sinniger
- National Institute for Biological Standards and Control, Blanche Lane, South Mimms, Hertfordshire EN6 3QG, United Kingdom
| | | | | | | | | |
Collapse
|
16
|
Zhu G, Allende ML, Jaskiewicz E, Qian R, Darling DS, Worth CA, Colley KJ, Young WW. Two soluble glycosyltransferases glycosylate less efficiently in vivo than their membrane bound counterparts. Glycobiology 1998; 8:831-40. [PMID: 9639544 DOI: 10.1093/glycob/8.8.831] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
Many Golgi glycosyltransferases are type II membrane proteins which are cleaved to produce soluble forms that are released from cells. Cho and Cummings recently reported that a soluble form of alpha1, 3-galactosyltransferase was comparable to its membrane bound counterpart in its ability to galactosylate newly synthesized glycoproteins (Cho,S.K. and Cummings,R.D. (1997) J. Biol. Chem., 272, 13622-13628). To test the generality of their findings, we compared the activities of the full length and soluble forms of two such glycosyltransferases, ss1,4 N-Acetylgalactosaminyltransferase (GM2/GD2/ GA2 synthase; GalNAcT) and beta galactoside alpha2,6 sialyltransferase (alpha2,6-ST; ST6Gal I), for production of their glycoconjugate products in vivo . Unlike the full length form of GalNAcT which produced ganglioside GM2 in transfected cells, soluble GalNAcT did not produce detectable GM2 in vivo even though it possessed in vitro GalNAcT activity comparable to that of full length GalNAcT. When compared with cells expressing full length alpha2,6-ST, cells expressing a soluble form of alpha2,6-ST contained 3-fold higher alpha2,6-ST mRNA levels and secreted 7-fold greater alpha2,6-ST activity as measured in vitro , but in striking contrast contained 2- to 4-fold less of the alpha2,6-linked sialic acid moiety in cellular glycoproteins in vivo . In summary these results suggest that unlike alpha1,3-galactosyltransferase the soluble forms of these two glycosyltransferases are less efficient at glycosylation of membrane proteins and lipids in vivo than their membrane bound counterparts.
Collapse
Affiliation(s)
- G Zhu
- Departments of Biological and Biophysical Sciences and Biochemistry and Molecular Biology, Schools of Dentistry and Medicine and James G. Brown Cancer Center, University of Louisville, Louisville, KY 40292, USA
| | | | | | | | | | | | | | | |
Collapse
|