Medical, Pharma, Engineering, Science, Technology and Business

^{1}Department of Physics and Astronomy, University of Southern California, Los Angeles, CA, USA

^{2}Department of Biomedical Engineering, University of Southern California, Los Angeles, CA, USA

^{3}Graduate College of Biomedical Sciences, Western University of Health Sciences, Pomona, CA, USA

- *Corresponding Author:
- Albash T

Department of Physics and Astronomy

University of Southern California

Los Angeles, CA, USA

**Tel:**+1 213 740 7492

**E-mail:**[email protected]

**Received date:** October 14, 2013; **Accepted date:** November 04, 2013; **Published date:** November 15, 2013

**Citation:** Albash T, Bouteiller JMC, Berger TW, Baudry M, Haas S (2013) Back Action on Neurotransmitters by Receptor Binding Reveals an Optimal Receptor Density Profile. J Comput Sci Syst Biol 6:327-336. doi:10.4172/0974-7230.1000129

**Copyright:** © 2013 Albash T, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Visit for more related articles at** Journal of Computer Science & Systems Biology

We discuss how integration of back action into coupled rate equations describing dynamical biophysical processes can lead the identification of optimized structural features. This approach is applied to analyze neural receptor binding and function. In functional receptor studies, the influence of ligand binding to the receptor on free ligand concentration in the synaptic cleft is rarely considered, especially when the number of ligand molecules vastly exceeds the number of receptors. Here we evaluate the role of ligand binding/unbinding to the receptor on ligand concentration and the resulting change in receptor dynamics using the example of glutamate interaction with the AMPA receptor subtype of glutamate receptors. We find a significant difference for AMPA receptor-mediated current between the free diffusion case, where binding/unbinding is neglected, and the case when glutamate binding to AMPA receptors is taken into account for evaluating free ligand concentration. Furthermore, taking into account receptor binding/unbinding reveals new properties of the receptor/neurotransmitter system, and in particular, indicates the existence of an optimum receptor density profile with an optimal radius where the total charge and peak current are maximal, a property that cannot be captured by the free diffusion case. This may provide an explanation for the disposition of AMPA receptors and the synaptic geometry based on the optimization of the receptor-mediated current.

Neurotransmitting receptors; Diffusion

Coupled rate equations are commonly used to describe biophysical time evolution processes on a high course-grained level. The rate constants are typically treated as fitting parameters, which contain information on the microscopic details underlying the dynamical processes. Typically, the rate equations form a differential linear system of equations given by

(1)

Where external sources are treated as being unperturbed by the evolution of the variables and are absorbed into the matrix A(t). In this work we consider a generalization of this approach, which includes the effect of the back action of the variables on the source terms. The resulting set of coupled differential equations is non-linear, and it encodes temporal feedback of the system on its environment.

In this study, we illustrate the effect of such back action on structural optimization by focusing on a particular example of postsynaptic receptors. Transmission of information through chemical synapses consists of three parts: release of chemical messengers (neurotransmitters) from one neuron into the extracellular space (synaptic cleft), neurotransmitter diffusion in the synaptic cleft, and neurotransmitter binding to receptors located on the receiving (postsynaptic) neuron. The process of neurotransmitter binding and unbinding (to receptors) is usually assumed not to modify the concentration of free neurotransmitter in the synaptic cleft, although this assumption is only valid when the number of neurotransmitter molecules is very high relative to the number of receptors and the diffusion rate is fast relative to the binding kinetics, which is not necessarily the case. How the dynamics of diffusing ligands are altered in the presence of interactions (such as binding/unbinding to ligand binding sites or consumption via oxidation) has been extensively studied both experimentally and theoretically. Very early examples involve the study of oxygen diffusion through tissue and its consumption via oxidation [1,2]. Early experiments [3-5] involving neurotransmitter diffusion indicated that the binding of neurotransmitters to receptors can effectively slow the diffusion of neurotransmitters, a phenomenon often referred to as buffered diffusion, especially when binding kinetics are much faster than diffusion rates, and additional experiments [6] have confirmed other pre- dictions of buffered diffusion models. There have also been numerous theoretical works, involving deterministic studies using partial differential equations [7–11] and stochastic studies using Monte Carlo simulations [12,13], which have revealed the importance of incorporating receptor-neurotransmitter interactions in evaluating free neurotransmitter concentration. However, most of these studies have focused on the changes in the temporal characteristics of ligand-receptor interactions due to buffered diffusion rather than on its resulting functional implications which we investigate in this study.

Incorporating receptor-neurotransmitter interactions is expected to have a significant effect when the number of receptors can be much larger than the number of ligand molecules and the kinetics of binding/ unbinding can be very fast relative to the free diffusion rate [5]. On the other hand, when the number of neurotransmitter molecules far exceeds the number of receptors and/or when binding is slow relative to diffusion rates, the effects of binding/unbinding are expected to be minimal. In these cases, it is generally assumed that the concentration of free neurotransmitter is not significantly modified by the binding of the neurotransmitter to the receptors. For instance, if the number of postsynaptic glutamate receptors is small relative to the number of glutamate molecules being released by a presynaptic event, it is logical to expect that the binding/unbinding of glutamate to its receptors will minimally affect free glutamate concentration [14,15] where binding/ unbinding is neglected in the modelling). This assumption is further supported by experimental results indicating that a single vesicle releases approximately 3000 glutamate molecules, and that the number of postsynaptic glutamate receptors is on the order of 100 receptors [16]. It would therefore seem that the number of receptors is sufficiently small relative to the number of glutamate molecules to neglect the effect of glutamate binding to its receptor on free glutamate concentration. However, since the association rate constant of glutamate to AMPA receptors is of the same order of magnitude as the diffusion rate, glutamate binding to receptors could significantly decrease the number of glutamate molecules as the distance from the release site increases, suggesting that it is not obvious that this assumption is valid.

The goal of this study is to apply a rate equation approach to determine whether neurotransmitter/receptor interactions would be different with a free diffusion model compared to a receptor-ligand diffusion model. We will investigate the temporal characteristics of ligand-receptor interactions due to buffered diffusion as well as its resulting quantitative functional implications on global synaptic AMPA receptor function. We consider a classical receptor kinetic model for the AMPA receptors, which has been described in numerous publications [17-19] using a glutamate diffusion profile unaffected by binding (free dif- fusion case). We then determine the dynamics of the receptor-limited diffusion model, this time considering the impact of glutamate binding/unbinding on free glutamate concentration profile in the cleft. We discover that there is a clear difference between the free diffusion and the receptor-limited diffusion model, thereby underscoring that the free diffusion case, which is most commonly used, may lead to inaccuracies in the calculation of transmitter/ receptor interactions. Moreover, when we take into consideration actual dimensions of the postsynaptic densities, our results indicate the existence of an optimal receptor density for a given radial distance from the release site, a feature that cannot be accounted for with the free diffusion model. Interestingly, this optimal density is on the order of the experimentally determined value. These results suggest that the dimensions of postsynaptic densities and the receptor density are set by some basic physical properties of ligand-receptor interactions.

**Glutamate-AMPA receptor model**

An illustration of the glutamate/receptor system we are considering is shown in **Figure 1**. Glutamate is released from the presynaptic terminal in the synaptic cleft in the form of a disk of radius r_{Glu} and height h, which corresponds respectively to the size of vesicle, and the distance between the pre- and post- synaptic elements (the transverse length of the cleft) [20,21]. Diffusion of glutamate takes place over a disk of receptors of radius r_{PSD}. Glutamate molecules are considered to be small enough to be described in the continuum limit, and we assume they are uniformly distributed throughout the cylinder. We assume that the height h is sufficiently small for the dynamics along the height direction to be negligible, such that any changes in glutamate concentration occur instantaneously along this direction. This has been shown to be a valid approximation for typical cleft heights [22]. In turn, the only relevant diffusion occurs in the plane transverse to the height.

The dynamics of AMPA receptors are governed by the kinetic schema for the 16-state AMPA receptor model [17], depicted in **Figure 2**. The model describes the opening of ion channels due to the binding of Glutamate to the receptor. The states in the schema are labeled by R,D,E,O, where D refers to being singly desensitized, E to doubly desensitized, and O to the opening of the receptor-associated ion channel. Furthermore, states are followed by a number 1,2,3,4 referring to how many glutamate molecules are bound. The states take values from 0 to 1 corresponding to the fraction of the total number of receptors in that state. Moving to the right along the diagram requires the binding of glutamate (measured in mM). Although a simpler model could be used to study the phenomena of interest, we selected this model because results generated by this model provide a good fit with numerous experimental data under a variety of experimental conditions. Furthermore, as will be discussed below, the existence of multiple binding sites for each molecule of receptors plays a significant role in the difference between the 2 cases. The parameters for the model given in **Table 1** were adapted from the original values used in Robert and Howe [17] to additionally fit the experimental results from Kessler et al. [23]. This optimization was presented in Bouteiller et al. [19]. We write the coupled rate equations associated with the kinetic schemes as follows:

(2)

Parameter | Value |
---|---|

k_{1} |
10 mM^{-1} ms^{-1} |

δ_{0} |
3.3×10^{-6} ms^{-1} |

δ_{1} |
0.42 ms^{-1} |

β |
0.55 ms^{-1} |

δ_{2} |
0.2 ms^{-1} |

k_{1} |
10 mM^{-1} ms^{-1} |

k_{-1} |
7 ms^{-1} |

γ_{0} |
0.001 ms^{-1} |

γ_{1} |
0.017 ms^{-1} |

α |
0.3 ms^{-1} |

γ_{2} |
0.035 ms^{-1} |

k_{-2} |
0.00041 ms^{-1} |

**Table 1:** Parameters used in the AMPA receptor model [19].

where x is the spatial coordinate, t is time, and for conciseness the states in the schema are labelled by yi, with i=0,..., 15 and y_{16} corresponds to glutamate concentration. We explain the form of the rate equation for glutamate and the function F (y(x, t)) below.

**Rate equation for glutamate concentration**

The number of glutamate molecules in the system is given by:

(3)

Where is the fraction of receptors in the state i with J glutamate molecules bound to receptors N_{cleft} is the number of molecules in the cleft (which are not bound to a receptor) and N is the number of receptors. Since the dynamics we are considering do not create or destroy glutamate molecules, taking a time derivative of equation (3) gives zero on the left-hand side, which then gives us an equation for the time evolution of the molecules in the cleft:

(4)

The derivative acting on the y_{i}‘s can be replaced with the coupled rate equations in equation (2). We can now convert this equation to describe the concentration of glutamate at a given location x as:

(5)

where ρ_{rec} is the density of receptors (in mM) associated with the spatial distribution of receptors. For notational conciseness, we define:

(6)

For simplicity we assume a “smearing” of the receptors such that they can be described by a continuous distribution. For example, a uniform distribution of receptors in a disk of radius r_{PSD}, ρ_{rec} would be given by:

(7)

Where h_{D} is the Debye length (μm) the appropriate conversion constant from μm^{-3} to mM:

(8)

Here, the Debye length corresponds to an interaction distance, representing the effective length along the height of the cylinder of Figure 1 from which the receptors can bind/unbind glutamate. To include diffusion, we simply include the diffusion term:

(9)

Where *D*_{Glu} is the (free) diffusion constant for glutamate in the cleft.

**Numerical methods**

The system of partial differential equations (n+1 in total, where n=16) we have to solve is given by:

(10)

Since the system and initial conditions are rotationally symmetric about the origin (this means we have Neumann boundary conditions for glutamate concentration at the origin), the angular dependence in the diffusion term can be ignored:

(11)

Where r denotes the radial coordinate. We discretize r with a grid of m elements with a maximum radius of 1 *μm*, and then we consider the time evolution of the *m×(n+1)* dependent variables. Since we are discretizing in the radial direction, our grid is formed from annular slices of width Δr =1μm / m, where we take *m*=250 for our simulations. At the maximum radius of 1 *μm*, we impose Dirichlet boundary conditions (an absorbing boundary) on the glutamate concentration to ensure that any glutamate molecule that has diffused to this point is irrecoverable. To solve this problem, we use the numerical method TR-BDF2 with an adaptive step size [24,25].

We focus on the case where receptors have a uniform distribution within a disk of fixed radial dimension, r_{PSD}. We set the values of the following parameters to:

(12)

Furthermore, we assume that the 3000 glutamate molecules are released as a disk of uniform concentration with radius *r*=0.026 *μm*.

Several physical parameters related to receptor function can be computed, including the current density *J _{i}(r,t)* generated by each open state:

(13)

Where is the density of receptors

α_{i} is the conductance of the *i*-th open state,

*V _{m}* is the voltage across the membrane.

Here we assume that V_{m} is a constant. We use the following values for the conductance [17]:

(14)

The total current density *J(t)* generated by the receptors is simply the sum of the individual *J _{i}(t)*. Since it is difficult to experimentally resolve the radial dependency of the current density, we also will be considering the current

(15)

The total current *I(t)* generated by the receptor is simply the sum of the individual *I _{i}(t)*. We give examples of both

**Figure 3:** Current density, current, and glutamate time evolution for a particular configuration of AMPA receptors. (A) Changes in current density as a function of time
and radius of the PSD. (B) Changes in current as a function of time for rPSD=102 nm and number of AMPA receptors=100. (C) Changes in the glutamate concentration
as a function of time and radius of the PSD.

Since the current *I(t)* is generically maximal for the single vesicle release event, a useful quantity to define is the maximum current *I _{peak}* attained at a time

I_{peak}= I(T_{peak})=max(I(T)), ∀ (16)

The total charge generated by the receptor distribution is given by:

(17)

with the total charge *Q* being the sum of the individual charges.

While there is no constraint for the receptor density in the equations we use above, the dimensions of the receptors and of the postsynaptic density impose some boundary conditions for this parameter. Thus, from electron microscopic studies, AMPA receptors are generally assumed to be cylindrical proteins with a diameter of approximately 15 nm [26], and the maximum receptor density that can be physically achieved in a postsynaptic density is about 5.7×10^{−3} *nm ^{−2}*:

(18)

Therefore, this value is highlighted in all our results and graphs (in the form of a dashed line). Furthermore, although this physical value is fixed, it is important to keep in mind that the dynamics of the system depend only on the ratio of *N/h _{D}*. Therefore, if we double our choice of the Debye length, we would see the same dynamics as our original Debye length with twice the receptor number. Note that in this study, we ignored the phenomenon of glutamate uptake (by astrocytic or neuronal glutamate transporters [27-29]), which significantly affects glutamate concentration in the synaptic cleft. However, incorporating this phenomenon would have a similar impact for both the free diffusion and the receptor-limited diffusion case (in fact, it would further emphasize the differences we observe, as shown below).

**Free diffusion**

In order to understand the differences introduced by the receptor binding/unbinding, we first consider the case in which the receptors do not modify the free glutamate concentration. This amounts to setting *F(y*) to zero in equation (9). We refer to this as the free diffusion case since glutamate diffuses as if no receptors are present. The free diffusion case can be understood as the limit of *N→0* of the rate equations (10). There cannot be any current generated in the free diffusion case since the number of receptors is technically zero. However, one way to proceed is to simply consider equations (10) with *N=0* but then use the definitions of *J _{i}(r,t)* and

Under these conditions, the behavior of *I(t)* and *Q* are shown in **Figure 4**. The linear relationship between current and charge and receptor density is simply due to the fact that the current is proportional to the number of receptors activated. For a fixed density of receptors, as radial dimension increases, the number of receptors that contribute to the current increases, thereby increasing the rate of change of current (and charge), which accounts for the increasing slope we observe. However, for sufficiently large radial dimension, the receptors located far from the release site experience a small enough glutamate concentration that the open states are barely populated. Therefore, the slope value should eventually asymptote to some maximum slope value.

**Figure 4:** Free diffusion case: Effects of receptor density on peak current and total charge. Peak current (A) and total charge (B) were calculated as a function of
receptor density. The graphs illustrate the relationship for the following radial dimensions: rPSD=42 nm (purple), rPSD=62 nm (blue), rPSD=82 nm (cyan), rPSD=102 nm
(green), rPSD=122 nm (yellow), rPSD=142 nm (orange), and rPSD=162 nm (red). The dashed line corresponds to the maximal receptor density for the AMPA receptor.

It is also interesting to consider the individual contribution of the various open states to the total current and charge (**Figure 5**). The key point is that, depending on the radial dimension, the open state that contributes most to the total peak current varies. For example, for a 42 *nm* radial dimension (a relatively small radial dimension), the O3 state contributes most, and then the O4 state, and finally the O2 state, independently of receptor density over the range we study. However, at 162 *nm*, the O2 state contributes the most, then the O3 state, and finally the O4 state. This is of course not surprising since for sufficiently large radial dimension, most of the receptors lie further away from the center, where glutamate concentration is lower, and the majority of the receptors are using the O2 state to generate current, a state in which binding of only two molecules of glutamate is sufficient to induce channel opening. A similar behavior occurs for the individual open state contribution to the charge, where for small radial dimension, the O3 state makes the largest contribution, while at much larger radial dimensions, and the O2 state has the largest contribution, independently of receptor density. Furthermore, we observe that the behavior of *I _{O2}* and

**Figure 5:** Free diffusion case: Contribution of the receptor’s open states on peak current and total charge. The individual contributions of the three open states,
O2, O3 and O4, to the peak current and total charge were calculated as a function of receptor density. The graphs illustrate the relationships for the following radial
dimensions: rPSD=42 nm (purple), rPSD=62 nm (blue), rPSD=82 nm (cyan), rPSD=102 nm (green), rPSD=122 nm (yellow), rPSD=142 nm (orange), and rPSD=162 nm (red). The
dashed line corresponds to the maximal receptor density for the AMPA receptor. (A) O2 Current. (B) O2 Charge. (C) O3 Current. (D) O3 Charge. (E) O4 Current.
(F) O4 Charge.

**Receptor-limited diffusion**

We then incorporate the receptor-limited diffusion to perform the same calculations as for the free diffusion case. An example of a direct comparison between the free diffusion and receptor-limited diffusion for peak current and total charge is shown in **Figure 6**. Interestingly, even for a relatively small number of receptors (<50), there is already a significant difference between the two results. In this example we use a postsynaptic density with a radius of 102 *nm*, and for *N=20* there is already about a 30% difference (corresponding to a difference of 1.3 *pA* for *V _{m}*=−70

**Figure 6:** Comparison of peak current and total charge between free diffusion and receptor-limited diffusion models. Peak current (A) and total charge (B) were
calculated with the free diffusion (blue) or receptor-limited diffusion (red) equations for a radial dimension of rPSD=102 nm. The dashed line corresponds to the
maximal receptor density for the AMPA receptor.

**Figure 7** shows the peak total current and the total charge for different radial dimensions. For very small receptor density, we observe the same linear relationship as in the free diffusion case, which is expected since the limit of small density corresponds to the free diffusion case. However, with increasing receptor densities, we observe very strong deviations from the behavior obtained with the free diffusion case. For all radial dimensions studied, both peak current and charge reach a maximum; although at somewhat different values of receptor density (the charge reaches a maximum at a higher receptor density). As receptor density increases, *I _{peak}* and

**Figure 7:** Receptor-limited diffusion case: Effects of receptor density on peak current and total charge. Peak current (A) and total charge (B) were calculated as
a function of receptor density. The graphs illustrate the relationship for the following radial dimensions: rPSD=42 nm (purple), rPSD=62 nm (blue), rPSD=82 nm (cyan),
rPSD=102 nm (green), rPSD=122 nm (yellow), rPSD=142 nm (orange), and rPSD=162 nm (red). The dashed line corresponds to the maximal receptor density for the
AMPA receptor.

Considering the critical physical receptor density, it appears that the only receptor distributions that can generate a maximum peak current are those with radial dimensions approximately greater than 90 nm, while those that can generate a maximum charge are those with radial dimensions approximately greater than 100 *nm*. These values correspond to 145 receptors for the former (90 *nm*) and 180 receptors for the latter (100 *nm*). (Recall the discussion at the end of section 4 about the close relationship of this value with the Debye length.) Furthermore, for large radial dimensions, the optimal receptor density decreases. This means that the optimal configuration of AMPA receptors for large radial dimensions will deviate significantly from a tight packing configuration. This is not unexpected since for larger radial configurations, glutamate concentration far away from the release site is small, and a decrease in receptor density will reduce competition between receptors for the limited number of glutamate molecules.

Again it is interesting to consider the individual contribution of the various open states to peak current and total charge (**Figure 8**). The results are strikingly different from those shown in Figure 5. One very important difference is that the O2 state can contribute more than the O3 state. This is an important result, since the O2 state is the only state that shows an increasing contribution with increasing receptor density (up to a point). This is perhaps not entirely surprising since any increase in receptor density would mean an immediate competition between receptors, which would result in a decline in O3 and O4 states’ contributions, and allowing the O2 state to play a bigger role. It is exactly this interplay that results in extreme values for peak current and charge.

**Figure 8:** Receptor-limited diffusion case: Contribution of the receptor’s open states on peak current and total charge. The individual contributions of the three open
states, O2, O3 and O4, to the peak current and total charge were calculated as a function of receptor density. The graphs illustrate the relationships for the following
radial dimensions: rPSD=42 nm (purple), rPSD=62 nm (blue), r=82PSD nm (cyan), rPSD=102 nm (green), rPSD=122 nm (yellow), rPSD=142 nm (orange), and rPSD=162 nm
(red). The dashed line corresponds to the maximal receptor density for the AMPA receptor. (A) O2 Current. (B) O2 Charge. (C) O3 Current. (D) O3 Charge. (E) O4
Current. (F) O4 Charge.

It is also interesting to analyze the unique patterns observed for *I _{O2}* and

**Comparison to experimental results**

We compare the results obtained with our two models to published experimental data involving one or two release events (Figure 1A) [30]. Experimental results show that for a single release event, the peak current (at a potential of −75 mV) is 5 pA, and the peak current doubles to 10 pA for two simultaneous release events. Unfortunately, available experimental data do not provide the number of glutamate molecules released, the density of AMPA receptors, and the radial extent of AMPA receptors. Thus, we choose an AMPA receptor configuration that allows us to fit the experimental data with both the free diffusion and receptor-limited model. We use *r _{PSD}*=102 nm for both models, and the number of AMPA receptors

**Figure 9:** Comparison of current response between the free diffusion model and the receptor-limited diffusion model. The free diffusion model (N=28, rPSD=102 nm,
blue curve) and the receptor-limited diffusion model (N=40, rPSD=102 nm, red curve) are compared with experimental data (black curve, extracted from reference
[30]) after one and two release events. (A) One release event, NGlu=3000. (B) Two release events. NGlu=6000. See text for details about the temporal discrepancy.

Although these results suggest that the receptor-limited model provides a better match with experimental data, this finding may not be true for all AMPA receptor configurations. For example, if for a particular configuration, both models result in low saturation of receptors after a single release event, it would be expected that both models would provide similar responses for two release events. Therefore, an important way to differentiate the two models is by controlling the number of receptors, i.e., changing the number or density of receptors while holding the number of glutamate molecules released fixed. As shown in Figure 6, our simulation results predict a substantial difference in the observed peak current as the density of receptors is scaled.

Our results indicate that taking into account the feedback controlled receptor-limited diffusion to calculate free glutamate concentration in the cleft significantly changes the patterns of important functional parameters of receptor function, such as total current and total charge through the receptor channel. The physical dimensions at which these effects becomes relevant are within the values experimentally observed (for a reasonable choice of Debye length), indicating that such consideration needs to be incorporated when modeling neurotransmitter/receptor interactions. Furthermore, our results indicate that considering the effect of binding/unbinding shed new lights on the understanding of receptor density and receptor distribution profile within the postsynaptic density. In particular, we discover that, for a uniform distribution of AMPA receptors on a disk, and for a given radial dimension, there is an optimal radial density profile of receptors where peak current or total charge generated by the receptors is maximal, which is a property that cannot be observed using a free diffusion model. This is a direct result of the competition between receptors for the limited number of glutamate molecules, and as expected, when the dimension of the postsynaptic density increases, the optimal receptor density decreases.

Using a typical size for the AMPA receptor, we determine a maximum receptor density limit, and our results indicate that for a given receptor density below this limit, there is a corresponding radial dimension for which peak current or peak charge is maximal. These values are about 90 *nm* for the former and 100 *nm* for the latter. While these values are reasonably close to the value of a typical postsynaptic density (between 50 and 500 *nm*), they are clearly on the low side compared to experiments. However, it is important to note that we have only considered the presence of AMPA receptors in the postsynaptic density and inclusion of NMDA receptors (which have higher affinity for glutamate than AMPA receptors) as well as glutamate transporters will modify this optimum value. The density of receptors at which the maximum occurs corresponds to approximately 145 or 180 receptors (for a Debye length of 1 *nm*) depending on whether we consider the maximum of peak current or charge.

We made a number of simplifying assumptions throughout this study. First, we assumed that the dynamics of glutamate diffusion perpendicular to the plane (i.e., across the cleft) plays a negligible role. This means that we assumed an instantaneous travel for glutamate upon release across the cleft. In previous studies, it was shown that this effect is negligible [22], and this assumption also greatly simplifies the numerical calculation. However, including these effects would make the role of the Debye length far more realistic, since the receptor would now experience a concentration gradient transverse to the plane. Second, the smearing of the receptor was a necessary assumption in order to solve the problem in the continuum limit. We tried to take into account the physical size of the receptor by focusing our analysis on realistic receptor densities, but it would be interesting to construct models that include finite-size corrections, while retaining the simplicity of the continuum limit. Third, the glutamate-receptor interaction is strictly local, meaning that a receptor at position *x* (in the plane) can only interact with glutamate molecules (therefore concentration) at position *x*. However, the finite Debye length implies that the receptor can interact with glutamate concentration at *x ± hD* as well. It would be interesting to study the effect of such non-local interactions in these models.

An important step forward will be to study a more realistic configuration for the receptors. It is generally assumed that AMPA receptors are distributed in a mosaic pattern at the majority of excitatory synapses [31]. The approach we have used in this study is general enough to accommodate such a distribution, and therefore studying the system integrating NMDA receptors, AMPA receptors, as well as metabotropic receptors and glutamate transporters is an important future direction for our work. Since AMPA receptors in this more realistic situation will be exposed to even less glutamate molecules, the effect of receptor binding/unbinding should be quite relevant. Furthermore, our results indicate that the optimum receptor distribution is likely not uniform. From our results, it is more likely that a tightly packed configuration with AMPA receptors close to the release site and decreasing density away from the release site might correspond to the optimum distribution. This result appears to agree with experimental observations for receptor distributions [32]. It would be very interesting to demonstrate that the integration of the receptorlimited diffusion provides a physical solution for the actual optimal receptor distribution in the postsynaptic density and to evaluate how general this finding could be in different neurotransmitter/receptor systems at different synapses. Furthermore, it would be interesting to adapt our techniques to better understand the clustering of receptors observed experimentally [32]. Finally, these results could also account for the generally admitted view that synaptic plasticity is the result of movements of AMPA receptors closer or further from the center of the postsynaptic density [33].

This work was supported in part by the National Institute of Biomedical Imaging and Bio Engineering (NIBIB, grant number 5 P41 EB001978-25) and the Office of Naval Research (ONR, grant number N00014-10-1-0685) for Theodore W. Berger, by the National Institute of Neurological Disorders and Stroke (NINDS, grant number 1R01N057128-01A2) for Michel Baudry, Theodore W. Berger, and Stephan Haas.

- Hill AV (1928) The diffusion of oxygen and lactic acid through tissues. Proceedings of the Royal Society of London Series B, Containing Papers of a Biological Character 104: 39-96.
- Gerard RW (1931) Oxygen diffusion into cells. Biological Bulletin 60: 245-268.
- Colquhoun D, Henderson R, Ritchie JM (1972) The binding of labelled tetrodotoxin to non-myelinated nerve fibres. J Physiol 227: 95-126.
- Katz B, Miledi R (1973) The binding of acetylcholine to receptors and its removal from the synaptic cleft. J Physiol 231: 549-574.
- Hartzell HC, Kuffler SW, Yoshikami D (1975) Post-synaptic potentiation: interaction between quanta of acetylcholine at the skeletal neuromuscular synapse. J Physiol 251: 427-463.
- Glavinovic MI, Law Min JC, Kapural L, Donati F, Bevan DR (1993) Speed of action of various muscle relaxants at the neuromuscular junction binding vs. buffering hypothesis. J Pharmacol Exp Ther 265: 1181-1186.
- Junge W, McLaughlin S (1987) The role of fixed and mobile buffers in the kinetics of proton movement. Biochim Biophys Acta 890: 1-5.
- Sala F, Hernández-Cruz A (1990) Calcium diffusion modeling in a spherical neuron. Relevance of buffering properties. Biophys J 57: 313-324.
- Aquino G, Endres RG (2010) Increased accuracy of ligand sensing by receptor diffusion on cell surface. Phys Rev E Stat Nonlin Soft Matter Phys 82: 041902.
- Holmes WR (1995) Modeling the effect of glutamate diffusion and uptake on NMDA and non-NMDA receptor saturation. Biophys J 69: 1734-1747.
- Boucher J, Kröger H, Sík A (2010) Realistic modelling of receptor activation in hippocampal excitatory synapses: analysis of multivesicular release, release location, temperature and synaptic cross-talk. Brain Struct Funct 215: 49-65.
- Stiles JR, Van Helden D, Bartol TM Jr, Salpeter EE, Salpeter MM (1996) Miniature endplate current rise times less than 100 microseconds from improved dual recordings can be modeled with passive acetylcholine diffusion from a synaptic vesicle. Proc Natl Acad Sci U S A 93: 5747-5752.
- Dilger JP (2010) Monte Carlo simulation of buffered diffusion into and out of a model synapse. Biophys J 98: 959-967.
- Rusakov DA, Kullmann DM (1998) Extrasynaptic glutamate diffusion in the hippocampus: ultrastructural constraints, uptake, and receptor activation. J Neurosci 18: 3158-3170.
- Xu-Friedman MA, Regehr WG (2004) Structural contributions to short-term synaptic plasticity. Physiol Rev 84: 69-85.
- Matsuzaki M, Ellis-Davies GC, Nemoto T, Miyashita Y, Iino M, et al. (2001) Dendritic spine geometry is critical for AMPA receptor expression in hippocampal CA1 pyramidal neurons. Nat Neurosci 4: 1086-1092.
- Robert A, Howe JR (2003) How AMPA receptor desensitization depends on receptor occupancy. J Neurosci 23: 847-858.
- Massicotte G, Baudry M (2004) Brain plasticity and remodeling of AMPA receptor properties by calcium-dependent enzymes. Genet Eng (N Y) 26: 239-254.
- Bouteiller JM, Allam SL, Hu EY, Greget R, Ambert N, et al. (2011) Integrated multiscale modeling of the nervous system: predicting changes in hippocampal network activity by a positive AMPA receptor modulator. IEEE Trans Biomed Eng 58: 3008-3011.
- Xie X, Liaw JS, Baudry M, Berger TW (1997) Novel expression mechanism for synaptic potentiation: Alignment of presynaptic release site and postsynaptic receptor. Proceedings of the National Academy of Sciences 94: 6983-6988.
- Ghaffari-Farazi T, Liaw JS, Berger TW (1999) Consequence of morphological al- terations on synaptic function. Neurocomputing 26-27: 17 - 27.
- Savtchenko LP, Rusakov DA (2007) The optimal height of the synaptic cleft. Proc Natl Acad Sci U S A 104: 1823-1828.
- Kessler M, Suzuki E, Montgomery K, Arai AC (2008) Physiological significance of high- and low-affinity agonist binding to neuronal and recombinant AMPA receptors. Neurochemistry International 52: 1383-1393.
- Coughran J WM, Grosse E, Rose D (1983) Cazm: A circuit analyzer with macro- modeling. Electron Devices, IEEE Transactions on 30: 1207-1213
- Bank R, Coughran JWM, Fichtner W, Grosse E, Rose D, et al. (1985) Transient simulation of silicon devices and circuits. Computer-Aided Design of Integrated Circuits and Systems, IEEE Transactions on 4: 436-451.
- Chen X, Winters C, Azzam R, Li X, Galbraith JA, et al. (2008) Organization of the core structure of the postsynaptic density. Proc Natl Acad Sci U S A 105: 4453-4458.
- Tong G, Jahr CE (1994) Block of glutamate transporters potentiates postsynaptic excitation. Neuron 13: 1195-1203.
- Bergles DE, Jahr CE (1998) Glial contribution to glutamate uptake at Schaffer collateral-commissural synapses in the hippocampus. J Neurosci 18: 7709-7716.
- Anderson CM, Swanson RA (2000) Astrocyte glutamate transport: review of properties, regulation, and physiological functions. Glia 32: 1-14.
- Pankratov YV, Krishtal OA (2003) Distinct quantal features of AMPA and NMDA synaptic currents in hippocampal neurons: implication of glutamate spillover and receptor saturation. Biophys J 85: 3375-3387.
- Fukazawa Y, Shigemoto R (2012) Intra-synapse-type and inter-synapse-type relationships between synaptic size and AMPAR expression. Curr Opin Neurobiol 22: 446-452.
- Tarusawa E, Matsui K, Budisantoso T, Molnár E, Watanabe M, et al. (2009) Input-specific intrasynaptic arrangements of ionotropic glutamate receptors and their impact on postsynaptic responses. J Neurosci 29: 12896-12908.
- Choquet D (2010) Fast AMPAR trafficking for a high-frequency synaptic transmission. Eur J Neurosci 32: 250-260.

Select your language of interest to view the total content in your interested language

- Advanced DNA Sequencing
- Algorithm
- Animal and Tissue Engineering
- Applications of Bioinformatics
- Artificial Intelligence Studies
- Artificial intelligence
- Artificial neural networks
- Big data
- Bioinformatics Algorithms
- Bioinformatics Databases
- Bioinformatics Modeling
- Bioinformatics Tools
- Biology Engineering
- Biostatistics: Current Trends
- Cancer Proteomics
- Chemistry of Biology
- Clinical Proteomics
- Cloud Computation
- Cluster analysis
- Comparative genomics
- Comparative proteomics
- Computational Chemistry
- Computational Sciences
- Computational drug design
- Computer Science
- Computer-aided design (CAD)
- Current Proteomics
- Data Mining Current Research
- Data algorithms
- Data mining applications in genomics
- Data mining applications in proteomics
- Data mining in drug discovery
- Data mining tools
- Data modelling and intellegence
- Data warehousing
- Ethics in Synthetic Biology
- Evolution of social network
- Evolutionary Optimisation
- Evolutionary algorithm
- Evolutionary algorithm in datamining
- Evolutionary computation
- Evolutionary science
- Experimental Physics
- Findings on Machine Learning
- Gene Synthesis
- Genome annotation
- Genomic data mining
- Genomic data warehousing
- Handover
- Human Proteome Project Applications
- Hybrid soft computing
- Industrial Biotechnology
- Knowledge modelling
- Machine Learninng
- Mapping of genomes
- Mass Spectrometry in Proteomics
- Mathematical Modeling
- Mathematics for Computer Science
- Meta genomics
- Microarray Proteomics
- Models of Science
- Molecular and Cellular Proteomics
- Multi Objective Programming
- Neural Network
- Ontology Engineering
- P4 medicine
- Physics Models
- Protein Sequence Analysis
- Proteogenomics
- Proteome Profiling
- Proteomic Analysis
- Proteomic Biomarkers
- Proteomics Clinical Applications
- Proteomics Research
- Proteomics Science
- Proteomics data warehousing
- Python for Bioinformatics
- Quantitative Proteomics
- Robotics Research
- Scientific Computing
- Simulation Computer Science
- Soft Computing
- Statistical data mining
- Studies on Computational Biology
- Swarm Robotics
- Swarm intelligence
- Synthetic Biology
- Synthetic Biology medicine
- Synthetic Biotechnology
- Synthetic Genomics
- Synthetic biology drugs
- Systems Biology
- Technologies in Computer Science
- Theoretical Chemistry
- Theoretical Computer Science
- Theoretical Issues in Ergonomics Science
- Theoretical Methods
- Theoretical and Applied Science

- International Conference on
**Metabolomics**and**Systems Biology**

June 11-12, 2018 London, UK - International Conference on
**Computational Biology**and**Bioinformatics**

Sep 05-06 2018 Tokyo, Japan - International Conference on Advancements in
**Bioinformatics**and**Drug Discovery**

November 26-27, 2018 Dublin, Ireland

- Total views:
**11909** - [From(publication date):

December-2013 - May 23, 2018] - Breakdown by view type
- HTML page views :
**8126** - PDF downloads :
**3783**

Peer Reviewed Journals

International Conferences 2018-19