help button home button Biophys. J.
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Gabdoulline, R. R.
Right arrow Articles by Wade, R. C.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Gabdoulline, R. R.
Right arrow Articles by Wade, R. C.
Biophysical Journal 85:1421-1428 (2003)
© 2003 The Biophysical Society

Concerted Simulations Reveal How Peroxidase Compound III Formation Results in Cellular Oscillations

Razif R. Gabdoulline *, Ursula Kummer {dagger}, Lars F. Olsen {ddagger} and Rebecca C. Wade *

* Molecular and Cellular Modeling Group, European Media Laboratory, Heidelberg, Germany; {dagger} Bioinformatics and Computational Biochemistry Group, European Media Laboratory, Heidelberg, Germany; and {ddagger} CelCom, Institute of Biochemistry and Molecular Biology, University of Southern Denmark, Odense, Denmark

Correspondence: Address reprint requests to Rebecca C. Wade, Molecular and Cellular Modeling Group, European Media Laboratory, Schloss-Wolfsbrunnenweg 33, D-69118 Heidelberg, Germany. Tel.: 49-622-153-3247; Fax: 49-622-152-2298; E-mail: rebecca.wade{at}eml.villa-bosch.de.


   ABSTRACT
TOP
 ABSTRACT
INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 APPENDIX A
 APPENDIX B
 ACKNOWLEDGEMENTS
 REFERENCES
 
A major problem in mathematical modeling of the dynamics of complex biological systems is the frequent lack of knowledge of kinetic parameters. Here, we apply Brownian dynamics simulations, based on protein three-dimensional structures, to estimate a previously undetermined kinetic parameter, which is then used in biochemical network simulations. The peroxidase-oxidase reaction involves many elementary steps and displays oscillatory dynamics important for immune response. Brownian dynamics simulations were performed for three different peroxidases to estimate the rate constant for one of the elementary steps crucial for oscillations in the peroxidase-oxidase reaction, the association of superoxide with peroxidase. Computed second-order rate constants agree well with available experimental data and permit prediction of rate constants at physiological conditions. The simulations show that electrostatic interactions depress the rate of superoxide association with myeloperoxidase, bringing it into the range necessary for oscillatory behavior in activated neutrophils. Such negative electrostatic steering of enzyme-substrate association presents a novel control mechanism and lies in sharp contrast to the electrostatically-steered fast association of superoxide and Cu/Zn superoxide dismutase, which is also simulated here. The results demonstrate the potential of an integrated and concerted application of structure-based simulations and biochemical network simulations in cellular systems biology.


   INTRODUCTION
TOP
 ABSTRACT
 INTRODUCTION
MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 APPENDIX A
 APPENDIX B
 ACKNOWLEDGEMENTS
 REFERENCES
 
There are a burgeoning number of programs aimed at creation of a "virtual cell" by simulation of cellular biochemical and signaling dynamics. Macroscopic kinetic simulations are used to model the complex multistep dynamics. However, experimental data on kinetic parameters are often missing, and this hinders construction and validation of models. With advances in protein structure determination, in particular arising from structural genomics programs, it is increasingly likely that macromolecular structural data will be available for proteins in the simulated cellular pathways. Here we demonstrate how this protein structure data can be exploited, by means of atomic-detail molecular simulations, to compute kinetic parameters that are fed into macroscopic biochemical network simulations. We show the value of this concerted simulation approach by applying it to the PO reaction.

The PO reaction (Scheeline et al., 1997Go), which displays oscillatory dynamics important for immune response (Olsen et al., 2003Go), involves many elementary steps and five different oxidation states of the peroxidase enzyme (see Fig. 1). It thus provides a good model system for complex biochemical behavior. In addition to the classical peroxidase reaction, peroxidase catalyzes the PO reaction. This involves the oxidation of an organic electron donor (typically NADH) by molecular oxygen. Referring to Scheeline et al. (1997)Go, the overall reaction can be written as

(1)
The electron transfer can be mediated by aromatic compounds such as melatonin that react with the enzyme to form radicals that in turn oxidize NADH (Olsen et al., 2001Go). When NADH and O2 are supplied continuously to a stirred aqueous solution at pH 5.0–6.5 containing peroxidase, a suitable aromatic compound, and methylene blue, the concentrations of the reactants and five enzyme intermediates oscillate (Nakamura et al., 1969Go) (see Fig. 2). The reaction intermediates include hydrogen peroxide and superoxide (Scheeline et al.). These reactive oxygen species were previously considered undesirable in cellular metabolism, but have recently been suggested to function as secondary messengers in certain cell signaling processes (Amit et al., 1999Go). The oscillatory dynamics appear to protect peroxidase from the harmful effects of the superoxide radical (Hauser et al., 2001Go) and, in the case of myeloperoxidase (MPO), to be involved in the activation of neutrophilic leukocytes and, therefore, of immune response (Olsen et al., 2003Go). The importance of the oscillatory dynamics due to MPO is highlighted by the observation that changes in the oscillations in activated neutrophils are associated with medical conditions such as septic shock and myocardial infarction (Kindzelskii and Petty, 2002Go). Therefore, an understanding of the determinants of cellular oscillations that reaches from the macroscopic level down to the detailed level of protein structure would aid rational, targeted drug design.



View larger version (12K):
[in this window]
[in a new window]
  		FIGURE 1  Schematic representation of the elementary steps involving enzymatic intermediates in the peroxidase-oxidase reaction. Several additional nonenzymatic steps also take place. AH represents NADH or certain aromatic compounds such as melatonin. The enzyme can adopt one of five states: Per2+ and Per3+, and compounds (co) I to III with a bound ligand. The step shown by the thick arrow (binding of superoxide to form compound III) is simulated by Brownian dynamics to compute the second-order rate constant. The complete set of reactions is given, together with rate constants, in Table A1 of Appendix A.

 


View larger version (24K):
[in this window]
[in a new window]
  		FIGURE 2  Time series from biochemical network simulations showing oscillation of the concentrations of NAD(P)H and superoxide using the parameters given in Table A1. Initial concentrations of myeloperoxidase and melatonin were 200 and 300 µM, respectively. All other initial concentrations were zero. The bottom figure shows the envelope of the oscillations of NAD(P)H plotted against the second-order rate constant for the association of ferric myeloperoxidase and superoxide.

 

View this table:
[in this window]
[in a new window]
  		TABLE A1  Elementary reaction steps and rate constants of the model of the peroxidase-oxidase reaction in activated neutrophils

 
The rate constants of most of the elementary steps of the PO reaction have been measured experimentally and, by fitting the few unknown constants, it is possible to perform simulations of the biochemical network (Scheeline et al., 1997Go). Such simulations have been performed primarily for horseradish peroxidase (HRP), but recently also for MPO, the focus of the present study, in the context of neutrophil activation (Olsen et al., 2003Go). Thus, the PO reaction is suitable for investigating the incorporation of kinetic parameters estimated from protein structure-based simulations. Experimental measurements of kinetic parameters can be costly and time-consuming, especially when they are fast and involve highly reactive compounds like the radical intermediates in this reaction. Therefore, we focus here on using structure-based BD simulations to estimate the association rate constant for superoxide binding to peroxidase. The rate constant for this elementary step has been determined experimentally for HRP (Scheeline et al., 1997Go), and for MPO at pH 7.8 (Kettle et al., 1988Go) but not at the physiological pH of 5. So far, a fitted value has been used for this rate constant when simulating the overall network (Olsen et al., 2003Go).

The rate of diffusional association of two molecules places an upper limit on the rate of a reaction between them. The rate constant for diffusional association can be computed by combining an analytical solution to the diffusion equation for two spherical particles (Smoluchowski, 1917Go) with data on the probability of reaction of the two molecules from BD simulations (Northrup et al., 1983Go). The BD simulations provide a means to account for the nonuniform shape, charge distribution, and reactivity of biomolecules when computing association rate constants. In the present BD simulations, the diffusion of superoxide to the enzyme active site was simulated.

We first carried out BD simulations for superoxide dismutases from two species, bovine (SOD-B) and Photobacterium leiognahi (SOD-P). These SODs are diffusion-controlled and well-characterized experimentally. Consequently, BD simulations for these SODs permitted calibration of our simulation method and provided a point of reference for the simulations with peroxidases. We then performed BD simulations for three peroxidases: MPO, the homologous lactoperoxidase (LPO), and the nonhomologous HRP. The rate constants computed for the peroxidases were entered into macroscopic biochemical network simulations. These showed oscillatory dynamics consistent with experimental measurements. Comparison of the determinants of the rate constants in the five enzymes revealed a novel enzymatic control mechanism: negative steering by the electrostatic potential of MPO depressed the rate constant for its association with superoxide, bringing the rate constant for this reaction into the range necessary for oscillatory behavior.


   MATERIALS AND METHODS
TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
RESULTS AND DISCUSSION
 APPENDIX A
 APPENDIX B
 ACKNOWLEDGEMENTS
 REFERENCES
 
Biochemical network model
The model is shown in Table A1 in Appendix A and was derived as described in Olsen et al. (2003)Go. Simulations were performed with the MADONNA software (www.berkeleymadonna.com) using the Rosenbrock routine for the numerical solution of stiff differential equations.

Molecular modeling
Protein crystal structures were taken from the PDB for bovine (2sod) and photobacterium leiognahi (1yai) Cu/Zn SOD and for HRP (7atj) and human MPO (1d2v). The structure of human LPO has not yet been determined experimentally. It was modeled by homology using the SwissModel automated modeling server (www.expasy.ch/swissmod). The LPO sequence has 67% and 56% identity to chains A and C, respectively, of MPO in the structure files 1d2v (human), 1mhl (human), and 1myp (dog), which were used as templates. The SODs were modeled as homodimers, each monomer having one Cu2+ and one Zn2+ ion. MPO was modeled as a homodimer with a heme and one bound Ca2+ ion (Shin et al., 2001Go). HRP was modeled as a monomer with a heme and two bound Ca2+ atoms (Shiro et al., 1986Go). LPO was modeled as a homodimer with a heme but without bound ions.

Polar hydrogen atom coordinates were added using the WHATIF software (Hooft et al., 1996Go) after assigning titration states as follows. For SOD, His41 was doubly protonated and His61, which ligates both Cu and Zn ions was unprotonated (-1 e charge). For the peroxidases, protonation states were assigned on the basis of pKa values computed with the UHBD program (Madura et al., 1994Go) using a Poisson-Boltzmann electrostatic model (Demchuk and Wade, 1996Go; Raquet et al., 1997Go); see Table B1 of Appendix B.


View this table:
[in this window]
[in a new window]
  		TABLE B1  Computed pKa values and assigned net charges for titratable residues in the peroxidases

 
Brownian dynamics simulations
Simulations were carried out with the SDA software (Gabdoulline and Wade, 1998Go) using forces derived from protein electrostatic potential grids computed with the UHBD program (Madura et al., 1994Go) by numerical solution of the finite difference linearized Poisson-Boltzmann equation.

OPLS (Jorgensen and Tirado-Rives, 1988Go) atomic radii and partial charges were assigned to protein atoms. The charges on the metal ions and their ligating residues in SOD were assigned following Shen et al. (1990)Go. Parameters for ferric heme were taken from Kuczera et al. (1990)Go.

The superoxide ion was modeled as two spheres, each having a charge of -0.5 e and radius of 1.5 Å, at 1.5 Å separation. The motion of the ion in the presence of a fixed enzyme was simulated. Superoxide was assigned a translational diffusion constant of 0.1283 Å2/ps, which corresponds to a hydrodynamic radius of 1.62 Å, as used in Sergi et al. (1994)Go. The rotational diffusion constant assigned of 3.9 10-2 rad2/ps was derived with this hydrodynamic radius using the Stokes-Einstein equation.

Electrostatic forces were computed from the protein electrostatic potential grids and the charges on superoxide atoms, taking into account the desolvation penalties (Elcock et al., 1999Go) that occur when a charged ligand approaches a protein's low dielectric cavity. The parameters used were the same as in a recent study of protein-protein association (Gabdoulline and Wade, 2001Go).

The enzyme was fixed and the diffusion of superoxide, in a continuum ionic solution exerting frictional and stochastic forces on the solute, was simulated by Brownian dynamics (Ermak and McCammon, 1978Go). Association was monitored by measuring the distance between the oxygen atoms of superoxide and the Cu or Fe ions.


   RESULTS AND DISCUSSION
TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
APPENDIX A
 APPENDIX B
 ACKNOWLEDGEMENTS
 REFERENCES
 
Validation of the Brownian dynamics protocol for the calculation of diffusional association rate constants
The rate of the reaction of superoxide with the diffusion-controlled enzyme, SOD, has been studied extensively by BD simulation (Folcarelli et al., 1999Go; Sergi et al., 1994Go; Sines et al., 1990Go) and experiment (Getzoff et al., 1992Go; Tainer et al., 1983Go). Consequently, we first performed simulations of superoxide association with SODs from two species, bovine (SOD-B) and Photobacterium leiognahi (SOD-P), that show different arrangements of the SOD monomer in the active homodimer (Bourne et al., 1996Go). This permitted us to define a reaction criterion for rate constant determination from our BD simulations that reproduced the experimental rates for SOD well (see Table 1). This "reaction criterion" is consistent with those that we have used previously (Gabdoulline and Wade, 2001Go). The rate constants are calculated by requiring approach of a superoxide atom to within 8 Å of the metal ion for a BD trajectory to be considered reactive. The computed rate constant is then divided by a scaling factor of 1.3. A distance of 8 Å to the Cu of SOD-B, for example, corresponds to a distance of <6 Å to one of the side-chain hydrogen bond donors of Arg141, Lys134, Thr135, and some of the main-chain atoms. The "reaction distance" of 8 Å to the Cu differs significantly from the distance of 4 Å used in some other SOD simulations (Sergi et al., 1994Go). This difference is due to several factors. First, we use desolvation penalties (Elcock et al., 1999Go) while computing the interaction of superoxide with a protein. Second, a 2-atom rather than 1-atom representation of superoxide is used; this effectively makes superoxide larger and hinders its approach to the copper ion. Third, we model charge redistribution between the metal ions and their ligands. This weakens the electrostatic steering of the superoxide specifically to the Cu, thus increasing the derived "reaction distance" value. Test simulations with these factors eliminated show that experimental rates are reproduced with a reaction distance of 4–4.5 Å from the Cu.


View this table:
[in this window]
[in a new window]
  		TABLE 1  Computed and experimental rate constants for the binding of superoxide to different enzymes

 
Application of the Brownian dynamics protocol to calculate diffusional association rate constants for peroxidases
We then applied an identical procedure and reaction criterion to simulate the association of superoxide with the three peroxidases: HRP, MPO, and LPO. The results are reported in Table 1. The rate constant for HRP at 100 mM ionic strength, which is ~2 orders-of-magnitude lower than that for SOD at 20 mM ionic strength, was correctly computed as 1.2 x 107 M-1s-1. The rate constant for MPO at pH 8 was computed as 8 x 105 M-1s-1, in agreement with the experimental determination (Kettle et al., 1988Go). The estimated rate constant for MPO at pH 5 was almost the same as that for HRP (1.1 x 107 M-1s-1), whereas the rate constant for LPO was higher (2.6 x 108 M-1s-1).

Dependence of the association rate constant on protein electrostatics
Comparison of the rate constants computed at low ionic strength and at infinite ionic strength (modeled by neglecting electrostatic interactions) shows that the rate of association of superoxide and SOD is enhanced by electrostatic steering, reaching 8.5 x 109 M-1s-1 at 20 mM ionic strength. This electrostatic steering has been seen in previous BD simulations for SOD (Folcarelli et al., 1999Go; Sergi et al., 1994Go; Sines et al., 1990Go) and other diffusion-influenced enzymes such as triose phosphate isomerase (Wade et al., 1998Go). HRP and LPO, on the other hand, show only modest electrostatic enhancement of their association rate with superoxide. This is compatible with their much lower rate constants. MPO has approximately the same rate as HRP at 100 mM but, surprisingly, shows electrostatic depression of its rate, which is an order-of-magnitude lower at 100 mM ionic strength than in the absence of electrostatic interactions.

Dependence of oscillatory behavior on the rate constant for peroxidase compound III formation
Simulation of the time-dependence of all the concentrations in the biochemical network model for MPO displayed in Table A1 of Appendix A reveals that the depression of the rate constant for superoxide binding to MPO has a crucial impact on the overall behavior of the system. In this model, MPO is located in a different compartment from one of its substrates (NAD(P)H). The rate constant computed at pH 5 is within the range of (0.4–8) x 107 M-1s-1 for oscillatory dynamics to occur, whereas the rate at pH 8 is below the lower border of this range and the rate in the absence of electrostatic interactions is well above the upper border (see Fig. 2).

Analysis of the sensitivity of oscillatory behavior to all the individual kinetic parameters reveals that the rate constant for association of superoxide to peroxidase (k4) is one of the most crucial parameters for oscillatory behavior (see Table 2). Only two out of the other 10 kinetic parameters show similar importance (the reaction of compound II with melatonin, k3, and the reaction of the NADP-radical with oxygen, k8). This comparison also reveals the robustness of the oscillatory behavior in this system to the other parameters, most of which can be changed by many orders of magnitude without destroying the system's ability to display oscillatory behavior. Therefore, a few crucial parameters have to be more strictly controlled (as is shown here for the case of superoxide association to peroxidase), whereas other parameters can be adjusted flexibly according to the needs of other parts of the system.


View this table:
[in this window]
[in a new window]
  		TABLE 2  Sensitivity of the oscillatory behavior toward the individual kinetic constants

 
Structural comparison of the active sites of SOD and the peroxidases
Despite both being peroxidases, HRP and MPO have negligible sequence identity. Thus, it appears puzzling that the two proteins can achieve approximately the same rate constant for the superoxide binding step. However, even though the overall structures are very different, the active sites show key similarities in the region around the superoxide binding site. This suggests that postdiffusional-encounter events may be similar (see Fig. 3). In particular, MPO, HRP, and SOD-B have a catalytic arginine with a guanadinium nitrogen atom ~6 Å from the metal ion. This arginine is thought to be crucial for the dismutation reaction in SOD (Fisher et al., 1994Go; Folcarelli et al., 1999Go; Wade et al., 1998Go).



View larger version (41K):
[in this window]
[in a new window]
  		FIGURE 3  Superposition of the active sites of Cu/Zn SOD-B (red) and the MPO (yellow) and HRP (blue) heme proteins. A conserved water site is shown by a sphere in the SOD and MPO structures. A cyanide molecule is shown in the same location in the HRP structure where a recent crystal structure shows superoxide binds (Berglund et al., 2002Go). Despite lack of sequence similarity, there are similar features in the active sites, such as the location of the guanadinium group of an arginine. The distal histidines of MPO and HRP superimpose well, even though the entrance channels to the active sites of these two peroxidases have completely different shapes and orientations. The channel in MPO extends into the page behind and to the right of the distal histidine, whereas that in HRP extends out of the page in front of the distal histidine.

 
Sequence and structural comparison of the peroxidases to reveal features important for determining rate constants for compound III formation
A structural superposition also yields insights into the different factors governing diffusional association rates. The active site entrance channels are oriented in different directions in the two proteins. MPO appears to have a larger entrance channel than HRP, and therefore the rates of diffusional association of uncharged substrates to MPO are larger. However, due to the presence of negatively charged residues near the active site entrance of MPO, the association rate for the negatively charged superoxide is depressed by a factor of 10 (see Fig. 4). These negatively charged residues are also present in eosinophil peroxidase, which has a 70% sequence identity to MPO and is likely to function in eosinophils in a similar way to MPO in neutrophils. On the other hand, these negatively charged residues are not conserved in other structurally similar peroxidases, such as LPO with a 56% sequence identity. Indeed, as a result of the more positive charge in this region (+5 e more), BD simulations show that the rate constant for association of superoxide with LPO should be approximately an order-of-magnitude greater than for MPO (see Table 1). This is consistent with the observation that LPO catalyzes an oscillating reaction in vitro at very low enzyme concentrations (Kummer et al., 1996Go).




View larger version (197K):
[in this window]
[in a new window]
  		FIGURE 4  Electrostatic potential of the MPO homodimer at pH 5 (a) and pH 8 (b) contoured at + 0.5 (blue) and - 0.5 (red) kT/e. The black arrows show the locations of channels into the active sites. The yellow contours represent the reaction criterion for superoxide binding and are 8 Å from the Fe atom. The difference in the region of negative contours near the active sites at the two pHs is due to a histidine residue that changes protonation state between these pHs.

 
The discovery of the use of negative electrostatic steering to depress a reaction rate is unlikely to be limited to MPO alone. Investigation of how widespread such an electrostatic control mechanism is should be facilitated by the increasing availability of kinetic parameters from high-throughput experiments and the use of combined structural and kinetic simulation techniques.


   APPENDIX A
TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 APPENDIX A
APPENDIX B
 ACKNOWLEDGEMENTS
 REFERENCES
 


   APPENDIX B
TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 APPENDIX A
 APPENDIX B
ACKNOWLEDGEMENTS
 REFERENCES
 


   ACKNOWLEDGEMENTS
TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 APPENDIX A
 APPENDIX B
 ACKNOWLEDGEMENTS
REFERENCES
 
We are grateful to the Klaus Tschira Foundation and the Bundesministerium für Bildung und Forschung (BMBF) for financial support.


   FOOTNOTES
 
Abbreviations used: MPO, myeloperoxidase; HRP, horseradish peroxidase; LPO, lactoperoxidase; PO reaction, peroxidase-oxidase reaction; BD, Brownian dynamics; SOD, superoxide dismutase; SOD-B, bovine SOD; SOD-P, SOD from Photobacterium leiognahi.

Submitted on February 6, 2003; accepted for publication April 21, 2003.


   REFERENCES
TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 APPENDIX A
 APPENDIX B
 ACKNOWLEDGEMENTS
 REFERENCES
 
Amit, A., A. L. Kindzelskii, J. Zanoni, J. N. Jarvis, and H. R. Petty. 1999. Complement deposition on immune complexes reduces the frequencies of metabolic, proteolytic and superoxide oscillations of migrating neutrophils. Cell. Immunol. 194:47–53.[Medline]

Berglund, G. I., G. H. Carlsson, A. T. Smith, H. Szöke, A. Henriksen, and J. Hajdu. 2002. The catalytic pathway of horseradish peroxidase at high resolution. Nature. 417:463–467.[Medline]

Bourne, Y., S. M. Redford, H. M. Steinman, J. R. Lepock, J. A. Tainer, and E. D. Getzoff. 1996. Novel dimeric interface and electrostatic recognition in bacterial Cu,Zn superoxide dismutase. Proc. Natl. Acad. Sci. USA. 93:12774–12779.[Abstract/Free Full Text]

Demchuk, E., and R. C. Wade. 1996. Improving the continuum dielectric approach to calculating pKas of ionizable groups in proteins. J. Phys. Chem. 100:17373–17387.

Elcock, A. H., R. R. Gabdoulline, R. C. Wade, and J. A. McCammon. 1999. Computer simulation of protein-protein association kinetics: acetylcholinesterase-fasciculin. J. Mol. Biol. 291:149–162.[Medline]

Ermak, D. L., and J. A. McCammon. 1978. Brownian dynamics with hydrodynamic interactions. J. Chem. Phys. 69:1352–1360.

Fisher, C. L., D. E. Cabelli, J. A. Tainer, R. A. Hallewell, and E. D. Getzoff. 1994. The role of arginine 143 in the electrostatics and mechanism of Cu,Zn superoxide dismutase: computational and experimental evaluation by mutational analysis. Proteins. 19:24–34.[Medline]

Folcarelli, S., F. Venerini, A. Battistoni, P. O'Neill, G. Rotilio, and S. Desideri. 1999. Toward the engineering of a super efficient enzyme. Biochem. Biophys. Res. Commun. 256:425–428.[Medline]

Gabdoulline, R. R., and R. C. Wade. 1998. Brownian dynamics simulation of protein-protein diffusional encounter. Methods. 14:329–341.[Medline]

Gabdoulline, R. R., and R. C. Wade. 2001. Protein-protein association: investigation of factors influencing association rates by Brownian dynamics simulations. J. Mol. Biol. 306:1139–1155.[Medline]

Getzoff, E. D., D. E. Cabelli, C. L. Fisher, H. E. Parge, M. S. Viezzoli, L. Banci, and R. A. Hallewell. 1992. Faster superoxide dismutase mutants designed by enhancing electrostatic guidance. Nature. 358:347–351.[Medline]

Hauser, M. J. B., U. Kummer, A. Z. Larsen, and L. F. Olsen. 2001. Oscillatory dynamics protect enzymes and possibly cells against toxic substances. Faraday Discuss. 120:215–227.[Medline]

Hooft, R. W. W., C. Sander, and G. Vriend. 1996. Positioning hydrogen atoms by optimizing hydrogen bond networks in protein structures. Proteins. 26:363–376.[Medline]

Jorgensen, W. L., and J. Tirado-Rives. 1988. The OPLS force field for proteins. Energy minimizations for crystals of cyclic peptides and crambin. J. Am. Chem. Soc. 110:1657–1666.

Kettle, A. J., D. F. Sangster, J. M. Gebicki, and C. C. Winterbourn. 1988. A pulse radiolysis investigation of the reactions of myeloperoxidase with superoxide and hydrogen peroxide. Biochim. Biophys. Acta. 956:58–62.[Medline]

Kindzelskii, A. L., and H. R. Petty. 2002. Apparent role of traveling metabolic waves in oxidant release by living neutrophils. Proc. Natl. Acad. Sci. USA. 99:9207–9212.[Abstract/Free Full Text]

Kuczera, K., J. Kuriyan, and M. Karplus. 1990. Temperature dependence of the structure and dynamics of myoglobin. A simulation approach. J. Mol. Biol. 213:351–373.[Medline]

Kummer, U., K. R. Valeur, G. Baier, K. Wegmann, and L. F. Olsen. 1996. Oscillations in the peroxidase-oxidase reaction: a comparison of different peroxidases. Biochim. Biophys. Acta. 1289:397–403.[Medline]

Madura, J. D., M. E. Davis, M. K. Gilson, R. C. Wade, B. A. Luty, and J. A. McCammon. 1994. Biological applications of electrostatic calculations and Brownian dynamics simulations. Comp. Chem. Rev. 5:229–267.

Nakamura, S., K. Yokota, and I. Yamazaki. 1969. Sustained oscillations in a lactoperoxidase, NADPH and O2-system. Nature. 222:794.[Medline]

Northrup, S. H., S. A. Allison, and J. A. McCammon. 1983. Brownian dynamics simulations of diffusion influenced bimolecular reactions. J. Chem. Phys. 80:1517–1524.

Olsen, L. F., U. Kummer, A. L. Kindzelskii, and H. R. Petty. 2003. A model of the oscillatory metabolism of activated neutrophils. Biophys. J. 84:69–81.[Abstract/Free Full Text]

Olsen, L. F., A. Lunding, F. R. Lauritsen, and M. Allegra. 2001. Melatonin activates the peroxidase-oxidase reaction and promotes oscillations. Biochem. Biophys. Res. Commun. 284:1071–1076.[Medline]

Raquet, X., V. Lounnas, J. Lamotte-Brasseur, J. M. Frere, and R. C. Wade. 1997. pKa calculations for Class A ß-lactamases: methodological and mechanistic implications. Biophys. J. 73:2416–2426.[Abstract/Free Full Text]

Scheeline, A., D. L. Olson, E. P. Williksen, G. A. Horras, M. L. Klein, and R. Larter. 1997. The peroxidase-oxidase oscillator and its constituent chemistries. Chem. Rev. 97:739–756.[Medline]

Sergi, A., M. Ferrario, F. Poltichelli, P. O'Neill, and A. Desideri. 1994. Simulation of superoxide-superoxide dismutase association rate for six natural variants. Comparison with the experimental catalytic rate. J. Phys. Chem. 98:10554–10557.

Shen, J., C. F. Wong, S. Subramaniam, T. A. Albright, and J. A. McCammon. 1990. Partial electrostatic charges for the active site center of Cu,Zn superoxide dismutase. J. Comp. Chem. 11:346–350.

Shin, K., H. Hayasawa, and B. Lonnerdal. 2001. Mutations affecting the calcium-binding site of myeloperoxidase and lactoperoxidase. Biochem. Biophys. Res. Comm. 281:1024–1029.[Medline]

Shiro, Y., M. Kurono, and I. Morishima. 1986. Presence of endogenous calcium ion and its functional and structural regulation in horseradish peroxidase. J. Biol. Chem. 261:9382–9390.[Abstract/Free Full Text]

Sines, J. J., S. A. Allison, and J. A. McCammon. 1990. Point charge distributions and electrostatic steering in enzyme/substrate encounter: Brownian dynamics of modified copper/zinc superoxide dismutases. Biochemistry. 1990:9403–9412.

Smoluchowski, M. V. 1917. Versuch einer mathematischen theorie der koagulationskinetik kolloider loeschungen. Z. Phys. Chem. 92:129–168.

Tainer, J. A., E. D. Getzoff, J. S. Richardson, and D. C. Richardson. 1983. Structure and mechanism of copper, zinc superoxide dismutase. Nature. 306:284–287.[Medline]

Wade, R. C., R. R. Gabdoulline, S. K. Luedemann, and V. Lounnas. 1998. Electrostatic steering and ionic tethering in enzyme-ligand binding: insights from simulations. Proc. Natl. Acad. Sci. USA. 95:5942–5949.[Abstract/Free Full Text]





This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Gabdoulline, R. R.
Right arrow Articles by Wade, R. C.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Gabdoulline, R. R.
Right arrow Articles by Wade, R. C.

HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Copyright © 2003 by the Biophysical Society.