Äîêóìåíò âçÿò èç êýøà ïîèñêîâîé ìàøèíû. Àäðåñ îðèãèíàëüíîãî äîêóìåíòà : http://library.biophys.msu.ru/PDF/3347.pdf Äàòà èçìåíåíèÿ: Wed Jun 7 18:46:30 2006 Äàòà èíäåêñèðîâàíèÿ: Mon Oct 1 19:42:30 2012 Êîäèðîâêà: Ïîèñêîâûå ñëîâà: molecules

I

NSTITUTE OF

PHYSICS PUBLISHING

PHYSICAL B

IOLOGY

Phys. Biol. 3 (2006) 121­129

doi:10.1088/1478-3975/3/2/004

Direct simulation of plastocyanin and cytochrome f interactions in solution
I B Kovalenko1,AMAbaturova1, P A Gromov2,DMUstinin1, E A Grachev2, G Yu Riznichenko1 and ABRubin1
1 2

Biological Faculty, Lomonosov Moscow State University, Leninskie Gory, Moscow 119992, Russia Physical Faculty, Lomonosov Moscow State University, Leninskie Gory, Moscow 119992, Russia

E-mail: kovalenko78@mail.ru

Received 25 December 2005 Accepted for publication 11 May 2006 Published 5 June 2006 Online at stacks.iop.org/PhysBio/3/121 Abstract Most biological functions, including photosynthetic activity, are mediated by protein interactions. The proteins plastocyanin and cytochrome f are reaction partners in a photosynthetic electron transport chain. We designed a 3D computer simulation model of diffusion and interaction of spinach plastocyanin and turnip cytochrome f in solution. It is the first step in simulating the electron transfer from cytochrome f to photosystem 1 in the lumen of thylakoid. The model is multiparticle and it can describe the interaction of several hundreds of proteins. In our model the interacting proteins are represented as rigid bodies with spatial fixed charges. Translational and rotational motion of proteins is the result of the effect of stochastic Brownian force and electrostatic force. The Poisson­Boltzmann formalism is used to determine the electrostatic potential field generated around the proteins. Using this model we studied the kinetic characteristics of plastocyanin­cytochrome f complex formation for plastocyanin mutants at pH 7 and a variety of ionic strength values.

1. Introduction
In green plants and algae light absorption and charge separation followed by photosynthetic electron and proton transfer take place in the photosynthetic membrane of chloroplasts. Protein complexes photosystem 1, photosystem 2, cytochrome bf complex and ATP-synthase are embedded in a photosynthetic membrane (figure 1). A mobile carrier plastocyanin, a small copper protein, transfers electrons from cytochrome f, a subunit of the cytochrome bf complex, to P700, the reaction center of photosystem 1 (Hope 2000, Gross 1996). Plastocyanin molecules are assumed to diffuse for long distances in the luminal space of chloroplasts, shuttling electrons between granal and stromal areas of the thylakoid membrane (Hope 2000, Albertsson 2001). A shift in the plastocyanin concentration between granal and stromal areas of the membrane during light­dark transition was interpreted as favoring the long distance transport between these compartments (Haehnel et al 1989). In native ° chloroplasts the luminal space is narrow (40­100 A) and contains proteins protruding through the thylakoid membrane. ° Since the size of plastocyanin (40 â 28 â 30 A) is about
1478-3975/06/020121+09$30.00

the thickness of the lumen, the diffusion of plastocyanin in the lumen is probably hindered by the membrane protein complexes. There is a conceptual difficulty in explaining the mechanism of plastocyanin diffusion in the lumen. It is not clear whether plastocyanin can rapidly diffuse over a distance of hundreds of nanometers, between grana stacks and stroma lamellae. The 30­50 µs lag phase of cytochrome f oxidation, after a saturating pulse of light, shows that it takes some time for the oxidized plastocyanin to diffuse from P700 to the nearest cytochrome bf complex (Haehnel et al 1980). This lag is followed by multicomponent relaxation kinetics that may correspond to different plastocyanin diffusion distances. It is widely recognized that electrostatic interactions play a crucial role in the binding of plastocyanin to its reaction partners, photosystem 1 and cytochrome f (Durell et al 1990, Gross 1996). Long-range electrostatic interactions govern the diffusion of plastocyanin to its docking sites at protein complexes and thus influence the rate of association, whereas short-range electrostatic interactions make plastocyanin molecules orient in the electrostatic field of the protein complex (Gross 1996). 121

© 2006 IOP Publishing Ltd Printed in the UK


I B Kovalenko et al

Figure 1. Schematic view of photosynthetic electron transport. Two thylakoid membranes and the luminal space between them are shown. Multiprotein complexes photosystem 1, photosystem 2, and cytochrome bf are embedded in the membrane. Mobile electron carrier protein plastocyanin (Pc) diffuses in the lumen. The arrows denote electron transport.

The limited diffusion of plastocyanin and its interaction with photosystem 1 and cytochrome bf complexes in the narrow thylakoid lumen are difficult to study by experimental techniques. A computer simulation model can give some information on plastocyanin diffusion and interaction in lumen. On one hand, this simulation should take into account the geometry of the luminal space packed with many protein molecules. On the other hand, it is necessary to consider electrostatic interactions of plastocyanin with its reaction partners in a thylakoid membrane. Brownian dynamics (BD) models are usually used to simulate the interaction of two individual protein molecules to investigate the process of complex formation, to predict the structure of the complex and to estimate the rate constant of protein association. In BD simulations the interaction of two proteins is investigated in solution, the proteins are treated as rigid or semi-rigid bodies, their geometric shape is considered at atomic resolution and electrostatic interactions of the proteins are presented in great detail. Simulation of plastocyanin and cytochrome f association by the Brownian dynamics method was done by Pearson and Gross (1998), Gross and Pearson (2003), Rienzo et al (2001) and Ullmann et al (1997). In Ullmann et al (1997), BD simulations were followed by molecular dynamics and quantum chemistry studies. These BD simulations cannot provide the concentration time courses for different biochemical species which are registered experimentally and which are usually used as a source of information about rate constants of biochemical reactions. In our studies we develop models of the limited diffusion of plastocyanin and its interaction with photosystem 1 and cytochrome bf complex in a narrow thylakoid lumen. In this paper we present the results on the simulation of spinach plastocyanin and turnip cytochrome f diffusion and interaction in solution, but not in the lumen. This is the first step in the development of the complete luminal model, where 122

we consider many plastocyanin molecules that compete to form complexes with PS1 and cytochrome bf. This study develops the method of the direct 3D simulation of electron transport processes in a thylakoid, which was described previously (Kovalenko et al 2003). In contrast to models described in Pearson and Gross (1998), Gross and Pearson (2003), Rienzo et al (2001) and Ullmann et al (1997), our model is multiparticle and includes many plastocyanin and cytochrome f molecules. We do not describe complex formation in detail (such as hydrophobic interactions and conformational changes of proteins which are important for the final complex formation). Our model uses a continuum electrostatic approach that describes molecules at the atomic level using a macroscopic description. From crystallographic studies, we know the structure of plastocyanin­cytochrome f encounter complex (Ubbink et al 1998, PDB ID is 2PCF). We use this structure in our simulation.

2. Description of the model
2.1. The model scene The model scene is a cubic 3D reaction volume (figure 2). Initially, the protein molecules are distributed randomly in the reaction volume. In our simulation, protein molecules are considered as Brownian particles. These particles move stochastically in a viscous medium guided by random Brownian force and electrostatic force. As shown in Doi and Edwards (1986), one can use the Langevin equation for the mathematical description of such a process. The Langevin equation for translational motion has the following form: x dx tr = fx (t ) + Fx . dt x where x is the coordinate, tr is the viscous friction coefficient at x, fx (t ) and Fx are the projections of the random Brownian force and electrostatic force onto the x axis, respectively.


Direct simulation of plastocyanin and cytochrome f interactions in solution

Figure 2. Visualization of the model scene of the computer simulation of protein plastocyanin and cytochrome f in solution. Protein molecules are visualized by ellipsoids of revolution.

Figure 3. Description of plastocyanin and cytochrome f by ellipsoids of revolution for calculation of the viscous friction coefficients of the molecules.

2.2. Interactions of proteins: simulation of the geometric shape of the proteins and protein collisions The interaction of two protein­electron carriers can be divided into four stages (Gross 1996): (1) convergence due to diffusion and electrostatic attraction; (2) interaction of separate parts of the molecules determined by their charge configuration, mutual orientation of the molecules, with the possible formation of a predocking complex; (3) final binding at the site where electron transfer can occur; and (4) electron transfer. In our simulation we consider the diffusion of molecules, electrostatic interactions and formation of a predocking complex. In our model proteins are considered as rigid bodies. They diffuse in the simulated reaction volume and can come into close proximity to each other but without overlapping. The 3D geometric surfaces of plastocyanin and cytochrome f are known at atomic resolution based on crystal structures of these proteins. We approximated the surfaces of the molecules by a set of spheres. With this approximation, two molecules overlap if the respective sets of spheres overlap. Geometric interpretation of the surface of the molecule by a small number of spheres (10­50 spheres), on one hand, ensures a sufficiently realistic representation of the molecular surface. On the other hand, it is less time consuming than calculation of each atom overlap, because the number of spheres is by two orders of magnitude less than the number of atoms. To verify the overlapping of two spheres is very easy--it is enough to know the radii and coordinates of their centers. 2.3. Simulation of the electrostatic interactions of cytochrome f and plastocyanin Electrostatic interactions play a key role in the binding of plastocyanin to its reaction partners, photosystem 1 and cytochrome f (Durell et al 1990, Hope 2000, Rienzo et al 2001, Gross 1996). Charged amino acid residues and partial charges of the proteins produce a heterogeneous electrostatic field. This field is screened by the polar water molecules and thus the electrostatic potential rapidly decreases with distance from the protein surface. If a protein is far away from other proteins, it moves under the effect of Brownian force but not electrostatic force, and hence its motion is determined only 123

The Langevin equation for rotational motion is
x rot

d = mx (t ) + Mx , dt

x where rot is the viscous friction coefficient for rotational motion around axis x, mx (t ) and Mx are the moments of the random Brownian force and electrostatic force relative to the x axis, respectively. The random Brownian force has a normal x 2kT tr . Here distribution with zero mean value and dispersion t k is the Boltzmann constant and T is the temperature. To simplify the calculation of viscous friction coefficients, the molecules of plastocyanin and cytochrome f were represented as ellipsoids of revolution with axes a, b and c (a is the axis of revolution, b = c). The axes and sizes of ellipsoids were chosen in such a way that the moment of inertia relative to the axis of revolution (axis a) was minimal and equal to the moment of inertia of the protein molecule. The obtained ellipsoids of revolution are apt representations of the shape of plastocyanin and cytochrome f (figure 3). Translational and rotational viscous friction coefficients of an ellipsoid are calculated as described in Bloomfield (2000). For the description of the position of each object we chose the Cartesian coordinate system with mirror boundary conditions at the edges (the scene coordinate system). Every object is assigned its local Cartesian coordinate system, which coincides with the axes of the respective ellipsoid. Such a choice of local coordinate system is very convenient because the viscous friction coefficients for the ellipsoid in this system are well known (Bloomfield 2000). The Langevin equation is solved numerically in the local coordinate system: on each time step, the displacements of the ellipsoid along its axes and rotations around its axes are calculated. The displacements and rotations found in the local coordinate system for each object are then recalculated in the scene coordinate system. The position of every object in the model is specified by the displacement and rotation of the local coordinate system of the object with respect to the scene coordinate system.


I B Kovalenko et al

Figure 4. The grid for calculation of the electrostatic potential field around a protein. Dark gray represents the cells occupied with the protein, light gray represents the cells of separation, white cells represent water.

by Brownian diffusion. As the protein comes close to other proteins, its behavior is directed by the electrostatic field of the latter so that a favorable position may be reached for the formation of a predocking complex. In our simulation a protein is presented as a low dielectric area ( = 2) with spatially fixed partial charges. The solution is presented as a high dielectric area ( = 80) with mobile charges (ions). When the proteins are separated from each other by more ° than 100 A, the electrostatic interactions are very weak because of the dielectric medium (water with ions) screening. In our simulation electrostatic interactions are considered only when proteins converge to some distance which we call the radius of the electrostatic interaction. Since a single protein is influenced by all other proteins inside the radius of electrostatic interaction, it is necessary to calculate the electrostatic force value and its moment. For this reason we chose a 3D rectangular region around each protein, set up a rectangular grid in the region and mapped the charges, dielectric constants and ionic strength to the grid. Then we calculate the electrostatic potential value for a given ionic strength and pH in every cell of the grid. We calculate this electrostatic potential map around all types of the objects in our model scene (both for plastocyanin and cytochrome f in the oxidized and reduced forms). The step of the grid is ° 2 A. The formal charges on the proteins were calculated as described in Durell et al (1990). The dielectric constant of the cells inside the protein was assigned 2; inside the solution-- 80; at the boundary protein solution--40; and the value of the ionic strength inside the protein was assigned 0 (figure 4). A Poisson­Boltzmann calculation (Ullmann and Knapp 1999) is used to determine the electrostatic potential grid around each type of the object. In the finite difference representation, the electrostatic potential of a given cell on the grid in the nth iteration is calculated by an iteration formula (1) depending on the potential values in the neighboring cells 124

Figure 5. Equipotential surfaces ­6.5 mV (light gray) and + 6.5 mV (dark gray) for reduced cytochrome f and oxidized plastocyanin (wild-type and mutant E59K/E60Q), calculated by the Poisson­Boltzmann equation; ionic strength is 100 mM, pH = 7, sol = 80, prot = 2. The lines connect the amino acid residues on plastocyanin and cytochrome f that were used in the simulation for calculation of the distance between the proteins. For mutant plastocyanin the black points represent atoms with changed charges relative to the wild type.

in the previous time step and the entire charge of the cell. In the nth iteration the electrostatic potential value in the cell is =
6 i =1 hi i 6 i =1 hi

+4q + h3 1 2
2 K

(1)

2 =

8NA e2 I , kB T

I=

cibulk Zi2 .
i =1

Here is the electrostatic potential, is the dielectric constant, is the charge density, cibulk is the concentration of the ith ion in the solution, Zi is the charge of the ith ion, e is the electron charge, T is the temperature (K), NA is the Avogadro number, I is the solution ionic strength and h is the step of the grid. Figure 5 demonstrates equipotential surfaces of the proteins cytochrome f and plastocyanin. After the calculation of the potential grid of each type of object in the model scene, we ran a Brownian dynamics simulation. In calculating the electrostatic force and its moment applied to a charge in a protein, it is necessary to calculate the potential gradient created by other proteins at the point where the charge is situated. To calculate the electrostatic force and its moment applied to the entire protein molecule, we geometrically sum the forces and their moments applied to every charge of the molecule. 2.4. Simulation of the plastocyanin­cytochrome f complex formation In the model, plastocyanin and cytochrome f molecules move randomly, guided by Brownian and electrostatic forces. As two molecules approach each other they are directed by the electrostatic field and may occupy a favorable position for


Direct simulation of plastocyanin and cytochrome f interactions in solution

docking. In the model it means that the distances between certain amino acid residues of the molecules become less than certain parameters of the model, which we call the docking distances. The corresponding amino acid residues are in a close contact in the crystal structure of the plastocyanin­ cytochrome f complex. From crystallographic studies we know the structure of the plastocyanin­cytochrome f encounter complex (we use the structure of the complex described in Ubbink et al (1998), PDB ID is 2PCF). In the model the distances between these amino acid residues determine the position and orientation of one molecule with respect to another. The program monitors the distances between chosen amino acid residues at every time step in the process of diffusion. If at some time step the distances are less than the docking distances, then the molecules can form a complex with probability p, which is a model parameter. We do not explicitly consider the process of complex formation. Instead, we simulate this event implicitly by introducing the probability, p, of the accidental complex formation. The accidental nature of a complex formation in every particular case is represented by a random choice of a number from the segment [0; 1]. At each time step, the algorithm generates a random number, uniformly distributed from 0 to 1. If the generated number is less than the given docking probability p, it means that the probability of the choice of this random number is the same as the probability of the complex formation. In other words, it means that the complex formation occurs. If the generated number is greater than the given docking probability p, the program flows to the next step. Provided the proteins still satisfy the docking conditions, the described procedure is repeated. As two molecules form a complex, it is considered that electron transfer occurs immediately and the molecules change their redox states. In our model we consider only electrostatic interactions and do not simulate explicitly the final complex formation and electron transfer. We do not consider hydrophobic interactions, possible complex rearrangements and any conformational changes of proteins which may be important for the final complex formation. We think that probability p is a convenient parameter that can take into account all these processes implicitly. Introducing probability p, we are able to simulate not only diffusion-controlled reactions but also the reactions in which the reaction rate is influenced by diffusion and some other processes. A recent paper concerning a cyanobacterial plastocyanin (Schlarb-Ridley et al 2005) suggests that the reaction between plastocyanin and cytochrome f is activation controlled, because the increase in the solution viscosity by addition of high molecular mass viscogens (Ficoll 70 or Dextran 70) had little or no effect on the second-order rate constant of the protein interaction. On the other hand, when the viscosity of the solution was raised by addition of low molecular mass viscogens, a marked inhibition of second-order rate constant was observed, as would be expected for a reaction in which diffusion plays an important part (a diffusion-controlled reaction). We suppose that the reaction is diffusion controlled because the true behavior of plastocyanin in a solution of high

Figure 6. Kinetic curve of cytochrome f oxidation, obtained from the simulation, and its approximation by the mass action law (see the formula in the figure). N0 is the total number of cytochrome f molecules, k is a second-order rate constant.

molecular mass viscogens such as Ficoll or Dextran is not known in detail.

3. Simulation algorithm parameters
Calculating the number of oxidized (PcII) and reduced (PcI) plastocyanin and cytochrome f (cyt f ) molecules at each time step, we obtain the kinetic curves of oxidation of cytochrome f and reduction of plastocyanin (figure 6). We can use the kinetic curves to compare the simulation results with experimental data and identify the model parameters. The second-order reaction rate constant can be determined by fitting kinetic curves by the mass action law. For this purpose we consider a simple kinetic scheme of a second-order reaction: PcII +cyt f - PcI +cyt f - , where k is a second-order rate constant. In all computer experiments the number of plastocyanin and cytochrome molecules was N0 = 30, the reaction volume was cubic with 70 nm on a side and the time step was 100 ps. Every simulation was repeated ten times and the obtained kinetic curves were averaged. 3.1. Rate constant dependence on the model parameters The protein association rate constant depends on numerous factors: the type of the proteins, diffusion coefficients, temperature, pH, ionic strength, etc. In our model, the second-order rate constant k is a function of the parameters of the model such that the probability of final complex formation, docking distances, time step, ionic strength, pH and temperature is k = F(p, ri,dt, I,pH, T). In this section we study the rate constant dependence on the docking distance and docking probability. In section 4 we provide the dependence of k on I. We studied the dependence of the rate of the plastocyanin and cytochrome f reaction on the docking distance. The docking distance is the model parameter that determines 125
k


I B Kovalenko et al

Table 1. Distances between contacting atoms in the complex plastocyanin­cytochrome f, obtained by the NMR method in Ubbink et al (1998) and distances chosen in our simulation as the condition of predocking complex formation. Amino acid residues of plastocyanin and cytochrome f, which are in close contact in plastocyanin­cytochrome f complex Asp42­Arg209 (D42-R209) Glu43­Lys187 (E43-K187) Asp44­Lys185 (D44-K185) Glu60­Lys58 (E60-K58) Glu59­Lys65 (E59-K65) His87­Tyr1 (H87-Y1) His87­Phe4 (H87-F4) ­ Numbers and names of the atoms to calculate the distances between Plastocyanin 591, OD2 607, HB 618, OD2 842, HA 825, OE2 1270, HE2 1270, HE2 1435, Cu Cytochrome f 3278, HH2 2930, HE 2895, HZ 912, HZ 1014, HE 21, HE2 65, CE2 3881, Fe R,A distance in a complex (experiment) 4.1 1.34 3.56 4.35 3.4 2.45 3.67 10.7­11.3 R,A distance in a complex (simulation) 18 18 18 25 ­ ­ ­ 40

Figure 7. Dependence of the rate constant of the reaction plastocyanin + cytochrome f on the docking distance (the distance between metal centers) at I = 100 mM for different docking probabilities p.

Figure 8. Dependence of the reaction rate on the docking probability p for the wild-type and mutant plastocyanins at I = 100 mM, pH = 7.0. The changes of charge on amino acid residues for different mutations and docking distances are presented. The values of distances between contacting amino acid residues are given in table 1.

the distance between the proteins when they are ready to form a complex capable for electron transfer. As the docking distance we chose the distance between metal centers (Cu in plastocyanin and Fe in cytochrome f). Figure 7 demonstrates the obtained dependences for three values of docking probabilities p (1, 0.01 and 0.001) and 100 mM ionic strength. As seen in figure 7, the rate of the plastocyanin and cytochrome f reaction increases with an increase in the docking distance value. The dependence is nonlinear: at large docking distances the rate of the reaction increases faster than at small docking distances. In fact, due to the electrostatic interactions and steric factors, diffusing protein molecules find themselves much more frequently at large distances from each other than at short distances. We also studied the dependence of the reaction rate on docking probability. In the simulation we changed the docking probability, while docking distances were fixed (see table 1). Figure 8 presents the dependence of the reaction rate on the docking probability for a wild type and some mutants of plastocyanin. The mutations modify the charge distribution on the surface of plastocyanin and thus significantly influence the process of protein docking, which is mainly determined 126

by electrostatic interactions. Some mutations lead to an increase in the rate constant, while some mutations actually decrease it. As seen in figure 8, there is a saturation dependence of the reaction rate on the docking probability. At high docking probabilities the rate of the reaction is almost constant. Proteins stay close to each other in the position ready to form a complex for a time period that depends on the diffusion coefficients and the time step value. At low docking probabilities they leave this position before the complex is formed. However, at high docking probabilities they have enough time to form a complex, which is observed as shown in figure 8, at p > 0.1. Figure 8 shows that for the plastocyanin mutant Q88E with high negative potential, the rate constant is higher than for the wild-type plastocyanin. This demonstrates that high electrostatic attraction between the negative patch on plastocyanin and the positive patch on cytochrome f leads to the increase in the rate constant. This is in accordance with experimental results (Kaant et al 1996) and other simulations (Rienzo et al 2001). All other plastocyanin mutations that we simulated have a lower negative potential at the binding


Direct simulation of plastocyanin and cytochrome f interactions in solution

(a )

Table 2. The experimental (Kaant et al 1996) and calculated rates of reaction plastocyanin and cytochrome f at I = 100 mM, pH = 6, for different mutants of plastocyanin. Model parameters are given in table 1; docking probability is 0.01. Charge change of the respective k,106 (M s) amino acid residue experiment ­ 0 -1 -1 -1 -2 -2 -3 -1 0 0 +1 0 +1 +1
-1

Mutation Wt Q88E D42N E43N E43K E43Q/D44N E59K/E60Q E59K/E60Q/ E43N

k,106 (M c) simulation 198 (±4) 345 (±7) 94 (±3) 80 (±3) 44 (±2) 40 (±3) 27 (±2) 6.4 (±0.9)

-1

185 (±20) 220 (±30) 76.5 (±1.5) 56.1 (±1) 29.3 (±0.7) 26.9 (±0.6) 13.8 (±0.2) 5.56 (±0.01)

(b )

Figure 9. The logarithm of the second-order rate constant k dependence on the root square of the ionic strength I for the wild-type and mutant plastocyanin and cytochrome f: (a) experimental data reprinted from Kaant et al 1996 Theroleof acidic residues of plastocyanin in its interaction with cytochrome f Biochim. Biophys. Acta 1277 115­26, Copyright (1996), with permission from Elsevier and from Crowley et al 2003 The parsley plastocyanin-turnip cytochrome f: a structurally distorted but kinetically functional acidic patch Biochem. J. 378(1) 45­57, reproduced with permission from the Biochemical Society; (b) simulation results (the parameter values are given in table 1, docking probability is 0.01). k is in (M s)-1, the ionic strength I is in M.

of estimated docking distances for the reaction in solution are given in table 1. On the other hand, we estimated the values of docking distances for which the calculated second-order rate constants for the wild-type and mutant plastocyanins at the ionic strength of 100 mM were close to the experimental values in Kaant et al (1996). As seen in table 1, the estimated docking distances are several times greater than the experimentally determined ones between plastocyanin and cytochrome f arranged in a unique complex. This is because in the model the representation of the surface of the protein molecule is quite rough (due to the approximation by a set of spheres) and it is impossible for the proteins in the model to come close at a distance of several angstroms. In the model the values of docking distances correspond to the position of separate molecules, cytochrome f and plastocyanin, in solution as they approach each other to form a complex. The value of docking probability p was estimated by comparison of calculated and experimental reaction rates. The good agreement with experimental rates for a wide range of mutants was achieved for p = 0.01. The comparison of experimental and model data at the ionic strength of 100 mM is represented in figure 9 and in table 1.

4. Results and discussion
On the basis of model parameters (docking distances and probability, figure 8 and table 1) we calculated rate constant values which fit the ex