Optimal plasticity for memory maintenance during ongoing synaptic change

  1. Dhruva V Raman  Is a corresponding author
  2. Timothy O'Leary  Is a corresponding author
  1. Department of Engineering, University of Cambridge, United Kingdom

Abstract

Synaptic connections in many brain circuits fluctuate, exhibiting substantial turnover and remodelling over hours to days. Surprisingly, experiments show that most of this flux in connectivity persists in the absence of learning or known plasticity signals. How can neural circuits retain learned information despite a large proportion of ongoing and potentially disruptive synaptic changes? We address this question from first principles by analysing how much compensatory plasticity would be required to optimally counteract ongoing fluctuations, regardless of whether fluctuations are random or systematic. Remarkably, we find that the answer is largely independent of plasticity mechanisms and circuit architectures: compensatory plasticity should be at most equal in magnitude to fluctuations, and often less, in direct agreement with previously unexplained experimental observations. Moreover, our analysis shows that a high proportion of learning-independent synaptic change is consistent with plasticity mechanisms that accurately compute error gradients.

Introduction

Learning depends upon systematic changes to the connectivity and strengths of synapses in neural circuits. This has been shown across experimental systems (Moczulska et al., 2013; Lai et al., 2012; Hayashi-Takagi et al., 2015) and is assumed by most theories of learning (Hebb, 1949; Bienenstock et al., 1982; Gerstner et al., 1996).

Neural circuits are required not only to learn, but also to retain previously learned information. One might therefore expect synaptic stability in the absence of an explicit learning signal. However, many recent experiments in multiple brain areas have documented substantial ongoing synaptic modification in the absence of any obvious learning or change in behaviour (Attardo et al., 2015; Pfeiffer et al., 2018; Holtmaat et al., 2005; Loewenstein et al., 2015; Yasumatsu et al., 2008; Loewenstein et al., 2011).

This ongoing synaptic flux is heterogeneous in its magnitude and form. For instance, the expected lifetime of dendritic spines in mouse CA1 hippocampus has been estimated as 1–2 weeks (Attardo et al., 2015). Elsewhere in the brain, over 70% of spines in mouse barrel cortex are found to persist for 18 months (Zuo et al., 2005), although these persistent spines exhibited large deviations in size over the imaging period (on average, a >25% deviation in spine head diameter).

The sources of these ongoing changes remain unaccounted for, but are hypothesised to fall into systematic changes associated with learning, development and homeostatic maintenance, and unsystematic changes due to random turnover (Rule et al., 2019; Mongillo et al., 2017; Ziv and Brenner, 2018). A number of experimental studies have attempted to disambiguate and quantify the contributions of different biological processes to overall synaptic changes, either by directly interfering with synaptic plasticity, or by correlating changes to circuit-wide measurements of ongoing physiological activity (Nagaoka et al., 2016; Quinn et al., 2019; Yasumatsu et al., 2008; Minerbi et al., 2009; Dvorkin and Ziv, 2016). Consistently, these studies find that the total rate of ongoing synaptic change is reduced by only 50% or less in the absence of neural activity or when plasticity pathways are blocked.

Thus, the bulk of steady-state synaptic changes seem to arise from fluctuations that are independent of activity patterns at pre/post synaptic neurons or known plasticity induction pathways. As such, it seems unlikely that their source is some external learning signal or internal reconsolidation mechanism. This is surprising, because maintenance of neural circuit properties and learned behaviour would intuitively require changes across synapses to be highly co-ordinated. To our knowledge, there is no theoretical account or model prediction that explains these observations.

One way of reconciling stable circuit function with unstable synapses is to assume that ongoing synaptic changes are localised to ‘unimportant’ synapses, which do not affect circuit function. While this may hold in particular circuits and contexts (Mongillo et al., 2017), at least some of the ongoing synaptic changes are likely associated with ongoing learning, which must somehow affect overall circuit function to be effective (Rule et al., 2020). Furthermore, this model does not account for the dominant contribution of fluctuations among those synapses that do not remain stable over time.

In this work we explore another, non-mutually exclusive hypothesis that active plasticity mechanisms continually maintain the overall function of a neural circuit by compensating changes that degrade memories and learned task performance. This fits within the broad framework of memory maintenance via internal replay and reconsolidation, a widely hypothesised class of mechanisms for which there is widespread evidence (Carr et al., 2011; Foster, 2017; Nader and Einarsson, 2010; Tronson and Taylor, 2007).

Compensatory plasticity can be induced by external reinforcement signals (Kappel et al., 2018), interactions between different brain areas and circuits (Acker et al., 2018), or spontaneous, network-level reactivation events (Fauth and van Rossum, 2019). Either way, we can conceptually divide plasticity processes into two types: those that degrade previously learned information, and those that protect against such degradation. We will typically refer to memory-degrading processes as ‘fluctuations’. While these may be stochastic in origin, for example due to intrinsic molecular noise in synapses, we do not demand that this is the case. Fluctuations will therefore account for any synaptic change, random or systematic, that disrupts stored information.

The central question we address in this work is how compensatory plasticity should act in order to optimally maintain stored information at the circuit level, in the presence of ongoing synaptic fluctuations. To do this, we develop a general modelling framework and conduct a first-principles mathematical analysis that is independent of specific plasticity mechanism and circuit architectures. We find that the rate of compensatory plasticity should not exceed that of the synaptic fluctuations, in direct agreement with experimental measurements. Moreover, fluctuations should dominate as the precision of compensatory plasticity mechanisms increases, where ‘precision’ is defined as the quality of approximation of an error gradient. This provides a potential means of accounting for differences in relative magnitudes of fluctuations in different neural circuits. We validate our theoretical predictions through simulation. Together, our results explain a number of consistent but puzzling experimental findings by developing the hypothesis that synaptic plasticity is optimised for dynamic maintenance of learned information.

Results

Review of key experimental findings

To motivate the main analysis in this paper we begin with a brief survey of quantitative, experimental measurements of ongoing synaptic dynamics. These studies, summarised in Table 1, provide quantifications of the rates of systematic/activity-dependent plasticity relative to ongoing synaptic fluctuations.

Table 1
Synaptic plasticity rates across experimental models, and the effect of activity suppression.
ReferenceExperimental systemTotal baseline synaptic change% synaptic change that is activity / learning-independent
Pfeiffer et al., 2018Adult mouse hippocampus40% turnover over 4 daysNA
Loewenstein et al., 2011Adult mouse auditory cortex>70% of spines changed size by >50% over 20 daysNA
Zuo et al., 2005Adult mouse (barrel, primary motor, frontal) cortex3–5% turnover over 2 weeks for all regions. 73.9 ± 2.8% of spines stable over 18 months (barrel cortex)NA
Nagaoka et al., 2016Adult mouse visual cortex8% turnover per 2 days in visually deprived environment. 15% in visually enriched environment. 7–8% in both environments under pharmacological suppression of spiking.50% (turnover)
Quinn et al., 2019Glutamatergic synapses, dissociated rat hippocampal culture28 ± 3.7% of synapses formed over 24 hr period. 28.6 ± 2.3% eliminated. Activity suppression through tetanus neurotoxin -light chain. Plasticity rate unmeasured.75% (turnover)
Yasumatsu et al., 2008CA1 pyramidal neurons, primary culture, rat hippocampusMeasured rates of synaptic turnover and spine-head volume change. Baseline conditions vs activity suppression (NMDAR inhibitors). Turnover rates: 32.8 ± 3.7% generation/elimination per day (control) vs 22.0 ± 3.6% (NMDAR inhibitor). Rate of spine-head volume change:67±17% (turnover). Size-dependent, but consistently >50% (spine-head volume)
Dvorkin and Ziv, 2016Glutamatergic synapses in cultured networks of mouse cortical neuronsPartitioned commonly innervated (CI) synapses sharing same axon and dendrite, and non-CI synapses. Quantified covariance in fluorescence change for CI vs non-CI synapses to estimate relative contribution of activity histories to synaptic remodelling62–64% (plasticity)
Minerbi et al., 2009Rat cortical neurons in primary cultureCreated ‘relative synaptic remodeling measure’ (RRM) based on frequency of changes in the rank ordering of synapses by fluorescence. Compared baseline RRM to when neural activity was suppressed by tetrodotoxin (TTX). RRM: 0.4 (control) vs 0.3 (TTX) after 30 hr.75% (plasticity)
Kasthuri et al., 2015Adult mouse neocortex (Three-dimensional post mortem reconstruction using electron microscopy).Data on 124 pairs of ‘redundant’ synapses sharing a pre/post-synaptic neuron was analysed in Dvorkin and Ziv, 2016. They calculated the correlation coefficient of spine volumes and post-synaptic density sizes between redundant pairs. This should be one if pre/post-synaptic activity history perfectly explains these variables.77% (post-synaptic density, r2=0.23).
66% (spine volume, r2=0.34)
Ziv and Brenner, 2018Literature review across multiple systemsCollectively these findings suggest that the contributions of spontaneous processes and specific activity histories to synaptic remodeling are of similar magnitudes50%

We focused on studies that measured ‘baseline’ synaptic changes that occur outside of any behavioural learning paradigm, and which controlled for stimuli that may induce widespread changes in synaptic strength. The approaches fall into two categories:

  1. Those that chemically suppress neural activity, and/or block known synaptic plasticity pathways, quantifying consequent changes in the rate of synaptic dynamics, in vitro (Yasumatsu et al., 2008; Minerbi et al., 2009; Quinn et al., 2019) and in vivo (Nagaoka et al., 2016). The latter study included a challenging experiment in which neural activity was pharmacologically suppressed in the visual cortex of mice raised in visually enriched conditions.

  2. Those that compare ‘redundant’ synapses sharing pre and post-synaptic neurons, and quantify the proportion of synaptic strength changes attributable to spontaneous processes independent of their shared activity history. These included in vitro studies that involved precise longitudinal imaging of dendritic spines in cultured cortical neurons (Dvorkin and Ziv, 2016). They also included in vivo studies, that used electron microscopy to reconstruct and compare the sizes of redundant synapses (Kasthuri et al., 2015) post mortem.

The studies in Table 1 consistently report that the the main component (more than 50%) of baseline synaptic dynamics is due to synaptic fluctuations that are independent of neural activity and/or easily identifiable plasticity signals. This is surprising because such a large contribution of fluctuations might be expected to disrupt circuit function. A key question that we address in this study is whether such a large relative magnitude of fluctuations can be accounted for from first principles, assuming that neural circuits need to protect overall function against perturbations.

The hypothesis we assumed is that some active plasticity mechanism compensates for the degradation of a learned memory trace or circuit function caused by ongoing synaptic fluctuations. We will thus express overall plasticity as a combination of synaptic fluctuations (task-independent processes that degrade memory quality) and compensatory plasticity, which counteracts this effect. There are various ways such a compensatory mechanism might access information on the integrity of overall circuit function, memory quality or ’task performance’. It could use external reinforcement signals (Kappel et al., 2018; Rule et al., 2020). Alternatively, such information could come from another brain region, as hypothesised in for example Acker et al., 2018, where cortical memories are stabilised by hippocampal replay events. Spontaneous, network-level reactivation events internal to the neural circuit itself could also plausibly induce performance-increasing plasticity (Fauth and van Rossum, 2019). Regardless, the decomposition of total ongoing plasticity into fluctuations and systematic plasticity allows us to derive relationships between both that are independent of the underlying mechanisms, which are not the focus of this study.

We must acknowledge that it is difficult, experimentally, to pin down and control for all physiological factors that regulate synaptic changes, or indeed to measure such changes accurately. However, even if one does not take the observations in Table 1 – or their interpretation – at face value, the conceptual question we ask remains relevant for any neural circuit that needs to retain information in the face of ongoing synaptic change.

Modelling setup

Suppose a neural circuit is maintaining previously learned information on a task. The circuit is subject to task-independent synaptic fluctuations which can degrade the quality of learned information. Meanwhile, some compensatory plasticity mechanism counteracts this degradation. Throughout this paper, we treat ‘memory’ and ‘task performance’ as interchangeable because our framework analyses the effect of synaptic weight change on overall circuit function. In this context, we ask:

if a network optimally maintains learned task performance, what rate of compensatory plasticity is required relative to the rate of synaptic fluctuations?

By ‘rate’ we mean magnitude of change in a given time interval. Our setup is depicted in Figure 1. We make the following assumptions, which are also stated mathematically in Box 1:

  1. The neural network has N adaptive elements that we call ‘synaptic weights’ for convenience, although they could include parameters controlling intrinsic neural excitability. We represent these elements through a vector 𝐰(t), which we call the neural network state. Changes to 𝐰(t) correspond to plasticity.

  2. Any state 𝐰(t) is associated with a quantifiable (scalar) level of task error, denoted F[𝐰(t)], and called the loss function. A higher value of F[𝐰(t)] implies greater corruption of previously learned information.

  3. The network state can be varied continuously. Task error varies smoothly with respect to changes in 𝐰(t).

  4. At any point of time, we can represent the rate of change (i.e. time-derivative) of the synaptic weights as

    𝐰˙(t)=𝐜˙(t)+ϵ˙(t).

    as discussed previously, which correspond to compensatory plasticity and synaptic fluctuations, respectively.

Motivating simulation results.

(a) We consider a network attempting to retain previously learned information that is subject to ongoing synaptic changes due to synaptic fluctuations and compensatory plasticity. (b) Simulations performed in this study use an abstract, rate based neural network (described in section Motivating example). The rate of synaptic fluctuations is constant over time. By iteratively increasing the compensatory plasticity rate in steps we observe a ‘sweet-spot’ compensatory plasticity rate, which is lower than that of the synaptic fluctuations, and which best controls task error. (c) A snapshot of the simulation described in b, at the point where the rates of synaptic fluctuations and compensatory plasticity are matched. Even as task error fluctuates around a mean value, individual weights experience systematic changes.

The magnitude and direction of plasticity may or may not change continually over time. Correspondingly, we may pick an appropriately small time interval, Δt, (which is not necessarily infinitesimally small) over which the directions of plasticity can be assumed constant, and write

(1) Δ𝐰(t)=Δ𝐜(t)+Δϵ(t),

where for any time-dependent variable x(t), we use the notation Δx(t):=x(t+Δt)-x(t). We regard Δ𝐜(t) and Δϵ(t) as coming from unknown probability distributions, which obey the following constraints:

  • Synaptic fluctuations Δϵ(t): We want to capture ‘task independent’ plasticity mechanisms. As such, we demand that the probability of the mechanism increasing or decreasing any particular synaptic weight over Δt is independent of whether such a change increases or decreases task error. A trivial example would be white noise, but systematic mechanisms, such as homeostatic plasticity, could also contribute (O’Leary, 2018; O'Leary and Wyllie, 2011).

  • Compensatory plasticity Δ𝐜(t): We demand that compensatory plasticity mechanisms change the network state in a direction of decreasing task error, on average. As such, they cause the network to preserve previously stored information, though not in general by restoring synaptic weights to their previous values following a perturbation.

Box 1.

Mathematical assumptions on plasticity.

To quantify memory quality/task performance we consider a loss function F[𝐰(t*)], which is twice differentiable in 𝐰(t). This loss function is simply an implicit measure of memory quality; we do not assume that the network explicitly represents F, or has direct access to it. Consider an infinitesimal weight-change Δ𝐰 over the infinitesimal time-interval Δt. We apply a second order Taylor expansion to express the consequent change in task error: ΔF=F[𝐰(t*)+Δ𝐰]-F[𝐰(t*)]:

(2) ΔF=ΔϵTF[w(t)]+ΔcTF[w(t)]+12ΔcT(2F[w(t)])Δc+12ΔϵT(2F[w(t)])Δϵ+ΔcT(2F[w(t)])Δϵ+𝒪(Δc+Δϵ23).

Here, F[𝐰(t*)] and 2F[𝐰(t)] represent the first two derivatives (gradient and hessian) of F[𝐰(t*)], with respect to a change in the weights 𝐰(t*). We assume that Δ𝐜 and Δϵ are sufficiently small (due to the short time interval) that the third-order term 𝒪(Δ𝐜+Δϵ23) can be ignored.

Next, we assume that Δ𝐜 and Δϵ are generated from unknown probability distributions. We place some constraints on these distributions. Firstly, synaptic fluctuations should be uncorrelated, in expectation, with the derivatives of F[𝐰], which govern learning. Accordingly,

(3a) 𝔼[ΔϵTF[𝐰(t*)]]=0,
(3b) 𝔼[Δ𝐜T(2F[𝐰(t*)])Δϵ]=0.

Secondly, we require that Δ𝐜 points in a direction of plasticity that decreases task error, for sufficiently small

(3c) Δc2ΔcTF[w(t)]<0.

Motivating example

Having described a generic modelling framework, we next uncover a key observation using a simple simulation.

Figure 1 depicts an abstract, artificial neural network trying to maintain a given input-output mapping over time, which is analogous to preservation of a memory trace or learned task. At every timestep, synaptic fluctuations corrupt the weights, and a compensatory plasticity mechanism acts to reduce any error in the input-output mapping (see Equation (1)). We fix the rate (i.e. magnitude per timestep) of synaptic fluctuations throughout. We increase the compensatory plasticity rate in stages, ranging from a level far below the synaptic fluctuation rate, to a level far above it. Each stage is maintained so that task error can settle to a steady state.

Two interesting phenomena emerge. The task error of the network is smallest when the compensatory plasticity rate is smaller than the synaptic fluctuation rate (Figure 1b). Meanwhile, individual weights in the network continually change even as overall task error remains stable due to redundancy in the weight configuration (Figure 1c), (see e.g. Rule et al., 2019 for a review).

In this simple simulation, we made a number of arbitrary and non-biologically motivated choices. In particular, we used an abstract, rate-based network, and synthesised compensatory plasticity directions using the biologically questionable backpropagation rule (see Materials and methods for full simulation details). Nevertheless, Figure 1 highlights a phenomenon that we claim is more general:

The ‘sweet-spot’ compensatory plasticity rate that leads to optimal, steady-state retention of previously learned information is at most equal to the rate of synaptic fluctuations, and often less.

In the remainder of the results section, we will build intuition as to when and why this claim holds. We will also explore factors influence the precise ‘sweet-spot’ compensatory plasticity rate.

The loss landscape

In order to analyse a general learning scenario that can accommodate biologically relevant assumptions about synaptic plasticity, we will develop a few general mathematical constructs that will allow us to draw conclusions about how synaptic weights affect the overall function of a network.

We first describe the ‘loss landscape’: a conceptually useful, geometrical visualisation of task error F[𝐰] (see also Figure 2). Every point on the landscape corresponds to a different network state 𝐰. Whereas any point on a standard three-dimensional landscape has two lateral (xy) co-ordinates, any point on the loss landscape has N co-ordinates representing each synaptic strength. Plasticity changes 𝐰, and thus corresponds to movement on the landscape. Any movement Δ𝐰 has both a direction Δ𝐰^ (where hats denote normalised vectors), and a magnitude Δ𝐰2. Meanwhile, the elevation of a point 𝐰 on the landscape represents the degree of task error, F[𝐰]. Compensatory plasticity improves task error, and thus moves downhill, regardless of the underlying plasticity mechanism.

Task error landscape and synaptic weight trajectories.

(a) Task error is visualised as the height of a ‘landscape’. Lateral co-ordinates represent the values of different synaptic strengths (only two are visualisable in 3D). Any point on the landscape defines a network state, and the height of the point is the associated task error. Both compensatory plasticity and synaptic fluctuations alter network state, and thus task error, by changing synaptic strengths. Compensatory plasticity reduces task error by moving ‘downwards’ on the landscape. (b) Eventually, an approximate steady state is reached where the effect of the two competing plasticity sources on task error cancel out. The synaptic weights wander over a rough level set of the landscape. (c) The effect of synaptic fluctuations on task error depends on local curvature in the landscape. Top: a flat landscape without curvature. Even though the landscape is sloped, synaptic fluctuations have no effect on task error in expectation: up/downhill directions are equally likely. Bottom: Although up/downhill synaptic fluctuations are still equally likely, most directions are upwardly curved. Thus, uphill directions increase task error more, and downhill directions decrease task error less. So in expectation, synaptic fluctuations wander uphill.

Understanding curvature in the loss landscape

Intuitively, one would expect task-independent synaptic fluctuations to increase task error. This is true even if fluctuations are unbiased in moving in an uphill or downhill direction on the loss landscape (see Equation (3a)) due to the curvature of the landscape (see Figure 2C). For instance, the slope (mathematically represented by the gradient F[𝐰]) at the bottom of a valley is zero. However, every direction is positively curved, and thus moves uphill. More generally, consider a fluctuation that is unbiased in selecting uphill or downhill directions, at a network state 𝐰. The fluctuation will increase task error in expectation if the total curvature of the upwardly curved directions at 𝐰 exceeds that of the downwardly curved directions, as illustrated in Figure 2c. We refer to such a state as partially trained. If all directions are upwardly curved, such as at/near the bottom of a valley, we refer to the state as highly trained. Mathematical definitions for these terms are provided in Box 2.

Box 2.

Curvature and the loss landscape.

Consider a fluctuation Δ𝐰 at a state 𝐰. The change in task error, to second order, can be written as

(4) F[w+Δw]F[w]ΔwTF[w]+12ΔwT2F[w]Δw

via a Taylor expansion. Suppose the fluctuation is task-independent. So it is unbiased with respect to selecting uphill/downhill, and more/less curved directions on the loss landscape. In this case

𝔼[F[𝐰]TΔ𝐰]=0

In expectation, Equation (4) thus becomes

E[F[w+Δw]F[w]]=Δw22Tr(2F[w])2N.

If Tr(2F[𝐰])>0, then the expected change in task error is positive, and we refer to the network state as ‘partially trained’. If additionally, 2F[𝐰]0, that is, Δ𝐰T2F[𝐰]Δ𝐰0 for any choice of Δ𝐰, then we refer to the network as highly trained. The ‘highly trained’ condition always holds in a neighbourhood of a local minimum of task error.

Comparison of the upward curvature of different plasticity directions plays an important role in the remainder of the section. Therefore, we introduce the following operator:

Qw[v]=v^T2F[w]v^.

Q𝐰[𝐯] is mathematical shorthand for the degree of curvature in the direction 𝐯, at point 𝐰 on the loss landscape, and is depicted in Figure 3a. Note that Q𝐰[𝐯] depends solely upon the direction, and not the magnitude, of 𝐯.

Quantifying effect of task error lanscape curvature on compensatory plasticity.

(a) Geometrical intuition behind the operator Q𝐰. The operator charts the degree to which a (normalised) direction is upwardly curved (i.e. lifts off the tangent plane depicted in grey). The red, shaded areas filling the region between the tangent plane and the upwardly curved directions are proportional to Q𝐰[v1], and Q𝐰[v2], respectively. (b) Compensatory plasticity points in a direction of locally decreasing task-error. Excessive plasticity in this direction can be detrimental, due to upward curvature (‘overshoot’). The optimal magnitude for a given direction is smaller if upward curvature (i.e. the Q-value) is large, as for cases (i) and (ii), and if the initial slope is shallow, as for case (ii). It is greater if the initial slope is steep, as for case (iii). This intuition underlies Equation (6) for the optimal magnitude of a given compensatory plasticity direction, which includes as a coefficient the ratio of slope to curvature. (c) Equation (11) depends upon the ratio of the upward curvatures in the two plasticity directions, Δ𝐜, and Δϵ. As illustrated, steep downhill directions often exhibit more upward curvature than arbitrary directions. In such cases, the optimal magnitude of compensatory plasticity should be outcompeted by synaptic fluctuations.

An expression for the optimal degree of compensatory plasticity during learning

The rates of compensatory plasticity and synaptic fluctuations, at time t, are 𝐜˙(t) and ϵ˙(t), respectively. These rates may change continually over time. Let’s temporarily assume they are fixed over a small time interval [t,t+Δt]. Thus,

(5) Δc=c˙(t)ΔtΔϵ=ϵ˙(t)Δt.

What magnitude of compensatory plasticity, Δ𝐜2, most decreases task error over Δt? The answer is

(6) Δc2=Δc^TF^[w]Qw[Δc]F[w]2.

A mathematical derivation is contained in Box 3, with geometric intuition in Figure 3b. Note that our answer turns out to be independent of the synaptic fluctuation rate ϵ˙(t). Here,

  • F[𝐰]2 represents the sensitivity of the task error to changes (i.e. the steepness of the loss landscape).

  • -Δ𝐜^TF^[𝐰] represents the accuracy of the compensatory plasticity direction in conforming to the steepest downhill direction on the loss landscape (in particular, their normalised correlation).

  • Q𝐰[Δ𝐜] represents the upward curvature of the compensatory plasticity direction. As shown in Figure 3b, excessive plasticity in an upwardly curved, but downhill, direction, can eventually increase task error. Thus, upward curvature limits the ideal magnitude of compensatory plasticity in the direction Δ𝐜^.

Box 3.

Optimal magnitude of compensatory plasticity.

Let us rewrite Equation (2), using the operator Q and omitting higher order terms, as justified in Box 1:

(7) ΔF=ΔϵTF[w(t)]+ΔcTF[w(t)]+12Δc22Qw(t)[Δc]+12Δϵ22Qw(t)[Δϵ]+ΔcT(2F[w(t)])Δϵ.

We can substitute our assumptions on synaptic fluctuations (Equations (3)) into Equation (7) to get

(8) 𝔼[ΔF]

Note that the requirement for assumption (3b) can be removed, but the alternative resulting derivation is more involved (see SI section two for this alternative).

We can differentiate Equation (8) in Δ𝐜2, to get:

dE[ΔF]dΔc2=Δc^TF[w(t)]+Δc2Qw(t)[Δc].

The root of this derivative gives a global minimum of the Equation (8) in Δ𝐜2, as long as Q𝐰(t*)[Δ𝐜]0 holds (justified in SI section 2.1). We get Equation (6), which defines the compensatory plasticity magnitude that minimises ΔF, and thus overall task error, at time t*+Δt.

For now, Equation (6) is valid only if the compensatory plasticity direction is fixed during Δt. If we want Equation (6) to also be compatible with continually changing compensatory plasticity directions, it needs to be valid for an arbitrarily small Δt. However, enacting a non-negligible magnitude Δ𝐜2* of plasticity over an arbitrarily small time interval Δt would require an unattainable, ‘infinitely-fast’ plasticity rate.

In fact, we show in the next section that our expression for Δ𝐜2* does become compatible with continuously changing plasticity at the end of learning, when task-error is stable.

Characterising the optimal rate of compensatory plasticity at steady state

Consider a scenario where task error is approximately stable. In this case, ΔF0 over Δt. In this scenario, Equation (6) simplifies to

(9a) Δc22Δϵ22=Qw[Δϵ]Qw[Δc],

as derived in Box 4 and illustrated geometrically in Figure 3c. We see that the magnitude Δ𝐜2* is proportional to Δϵ2, which is itself proportional to Δt from Equation (5), given some fixed rate of synaptic fluctuations. Thus, Δ𝐜2* is attainable even as Δt shrinks to zero, and is thus compatible with continually changing compensatory plasticity directions. In this case, Equation (9) can be rewritten as

(9b) c˙(t)2,2ϵ˙(t)22=Qw[ϵ˙(t)]Qw[c˙(t)].

Equation (9) is a key result of the paper. It applies regardless of the underlying plasticity mechanisms that induced Δ𝐜 and Δϵ. It is compatible with continually or occasionally changing directions of compensatory plasticity (i.e. infinitesimal or non-infinitesimal Δt). It says that the optimal compensatory plasticity rate, relative to the rate of synaptic fluctuations, depends on the relative upward curvature of these two plasticity directions on the loss landscape.

A corollary is that the optimal rate of compensatory plasticity is greater during learning than at steady state. If we substitute the steady-state requirement: 𝔼[ΔF]=0, with the condition for learning: 𝔼[ΔF]<0, in the derivation of Box 4, then we get

(10) Δc22Δϵ22Qw[Δϵ]Qw[Δc].

Indeed, the faster the optimal potential learning rate 𝔼[ΔF], the greater the optimal compensatory plasticity rate. Thus Δ𝐜2* decreases as learning slows to a halt, eventually reaching the level of Equation (9b).

Box 4.

Optimal compensatory plasticity magnitude at steady state error.

Let us substitute the special condition 𝔼[ΔF]=0 (steady-state task error) into Equation (8). This gives

0=ΔcTF[w(t)]+12Δc22Qw(t)[Δc]+12Δϵ|22Qw(t)[Δϵ].

Next, we substitute in our optimal reconsolidation magnitude (Equation (6)). This gives

0=12Δc22Q[Δc]+12Δϵ22Q[Δϵ],

which in turn implies the result (Equation (9)).

Note that Equation (9) is only valid when both the numerator and denominator of the right hand side are both positive. The converse is unlikely in a partially trained network, and impossible in a highly trained network (see SI section 2.1).

Main claim

We now claim that generically, the optimal compensatory plasticity rate should not outcompete the rate of synaptic fluctuations at steady state error. We will first provide geometric intuition for our claim, before bolstering with analytical arguments and making precise our notion of ‘generically’.

From Equation (9), our main claim holds if

(11) Qw[Δc]Qw[Δϵ],

that is, Δ𝐜 points in a more upwardly curved direction than Δϵ. When would this be true?

First consider Δϵ. Statistical independence from the task error means it should point in an ‘averagely’ curved direction. Mathematically (see SI secton 2.1), this means

(12) E[Qw[Δϵ]]=Tr(2F[w])N.

Our assumption of ‘average’ curvature fails if synaptic fluctuations are specialised to ‘unimportant’ synapses whose changes have little effect on task error. In this case Q𝐰[Δϵ] would be even smaller, since Δϵ would be constrained to consistently shallow, less-curved directions. Thus, this possibility does not interfere with our main claim.

For Equation (11) to hold, Δ𝐜 should point in directions of ‘more-than-average’ upward curvature. This follows intuitively because a steep downhill direction, which effectively reduces task error, will usually have higher upward curvature than an arbitrary direction (see Figure 3c for intuition). It remains to formalise this argument mathematically, and consider edge cases where it doesn’t hold.

Dependence of the optimal magnitude of steady-state, compensatory plasticity on the mechanism

Compensatory plasticity is analogous to learning, since it acts to reduce task error. We do not yet know the algorithms that neural circuits use to learn, although constructing biologically plausible learning algorithms is an active research area. Nevertheless, all the potential learning algorithms we are aware of fit into three broad categories. For each category, we shall show why and when our main claim holds. We will furthermore investigate quantitative differences in the optimal compensatory plasticity rate, across and within categories. A full mathematical justification of all the assertions we make is found in SI section 1.3.

We first highlight a few general points:

  • For any compensatory plasticity mechanism, Q𝐰[Δ𝐜] depends not only on the algorithm, but the point 𝐰 on the landscape. We cannot ever claim that Equation (11) holds for all network states.

  • We calculate the expected value of Q𝐰[Δ𝐜] for an ‘average’, trained, state 𝐰, across classes of algorithm. This corresponds to a plausible best-case tuning of compensatory plasticity that a neural circuit might be able to achieve. Any improvement would rely on online calculation of Q𝐰[Δ𝐜], which we do not believe would be plausible biologically.

Learning algorithms attempt to move to the bottom of the loss landscape. But they are blind. Spying a distant valley equates to ‘magically’ predicting that a very different network state will have very low task error. How do they find their way downhill? There are three broad strategies (Raman and O'Leary, 2021):

  • 0th order algorithms take small, exploratory steps in random directions. Information from the change in task error over these steps informs retained changes. For instance, steps that improve task error are retained. A notable 0-order algorithm is REINFORCE (Williams, 1992). Many computational models of biological learning in different circuits derive from this algorithm (Seung, 2003; Fee and Goldberg, 2011; Bouvier et al., 2018; Kornfeld et al., 2020).

  • 1st order algorithms explicitly approximate/calculate, and then step down the locally steepest direction (i.e. the gradient F[𝐰]). The backpropagation algorithm implements perfect gradient descent. Many approximate gradient descent methods with more biologically plausible assumptions have been developed in the recent literature (see e.g. Murray, 2019; Whittington and Bogacz, 2019; Bellec et al., 2020; Lillicrap et al., 2016; Guerguiev et al., 2017, and Lillicrap et al., 2020 for a review).

  • 2nd order algorithms additionally approximate/calculate the hessian 2F[𝐰], which provides information on local curvature. They look for descent directions that are both steep, and less upwardly curved. We doubt it is possible for biologically plausible learning rules to accurately approximate the hessian, which has N2 entries representing the interaction between every possible pair of synaptic weights.

Table 2 shows the categories for which our main claim holds.

Table 2
Table elements highlighted in teal correspond to scenarios in which our main claim holds, as Equation (11) is satisfied.
Quadratic F[𝐰]Nonlinear F[𝐰], low steady-state errorNonlinear F[𝐰], high steady-state error
0th order algorithmQ[Δ𝐜]Q[Δϵ]Q[Δ𝐜]Q[Δϵ]Q[Δ𝐜]Q[Δϵ]
0st order algorithmQ[Δ𝐜]Q[Δϵ]Q[Δ𝐜]Q[Δϵ]Q[Δ𝐜]Q[Δϵ]
0nd order algorithmQ[Δ𝐜]Q[Δϵ]Q[Δ𝐜]Q[Δϵ]Q[Δ𝐜]Q[Δϵ]

We first consider the simplest case of a quadratic loss function F[𝐰]. Here, directions of curvature in any direction are constant (mathematically, the hessian 2F[𝐰] does not vary with network state). Moreover, the gradient obeys a consistent relationship with the hessian:

(13) F[w]=2F[w](ww).

Components of (𝐰-𝐰*) with high upward curvature are magnified under the transformation 2F[𝐰*], since they correspond to eigenvectors of 2F[𝐰*] with high eigenvalue. Conversely, components with low upward curvature are shrunk. As the gradient F[𝐰] is the output of such a transformation from Equation (13), this suggests it is biased towards directions of high upward curvature. Indeed, we can quantify this bias. Let {λi} be the eigenvalues of 2F[𝐰*], and {ci} the projections of the corresponding eigenvectors onto 𝐰-𝐰*. Then

(14) Qw[F[w]]=i=1Nci2λi3i=1Nci2λi2.

The value of Equation (14) depends on the values {ci}. In the ‘average’ case, where they are equal, and 𝐰-𝐰* is thus a direction of ‘average’ curvature, Q𝐰[F[𝐰]]Q𝐰[Δϵ] holds. This inequality gap widens with increasing anisotropy in the curvature of different directions (i.e. with a wider spread of eigenvalues λi, corresponding to more elliptical/less circular level sets in the illustration of Figure 4b). Indeed, simulation results in Figure 5—figure supplement 1 (top row) show how the ratio Δ𝐜2:Δϵ2 that optimises steady-state task error is significantly less than one, in a quadratic error function where compensatory plasticity accurately follows the gradient, and for different synaptic fluctuation rates.

Geometric intuition for the optimal magnitude of different compensatory plasticity directions.

Colours depict level sets of the loss landscape. Elliptical level sets correspond to a quadratic loss function (which approximates any loss function in the neighbourhood of a local minimum). In c and d, we depict compensatory plasticity and synaptic fluctuations as sequential, alternating processes for illustrative purposes, although they are modelled as concurrent throughout the paper. (a) Compensatory plasticity directions locally decrease task error, so point from darker to lighter colours. Optimal magnitude is reached when the vectors ‘kiss’ a smaller level set, that is, intersect that level set while running parallel to its border. Increasing magnitude past this past this point increases task error, by moving network state to a higher-error level set. (b) If compensatory plasticity is parallel to the gradient (i.e. it enacts gradient descent), then it runs perpendicular to the border of the level set on which it lies (i.e. the tangent plane). This is shown explicitly for the ‘exact gradient’ direction of plasticity. The optimal magnitude of plasticity in this direction is smaller than that of a corrupted gradient descent direction, even though the former is more effective in reducing task error, because the exact gradient points in a more highly curved direction. (c) Synaptic fluctuations of a certain magnitude perturb the network state. The optimal magnitude of compensatory plasticity (in the exact gradient descent direction, for this example) is significantly smaller than that of the synaptic fluctuations, using the geometric heuristic explained in (a). If the magnitude of compensation increased to match the synaptic fluctuation magnitude there would be overshoot, and task error would converge to a higher steady state. (d) If compensatory plasticity mechanisms can perfectly calculate both the local gradient and hessian (curvature) of the loss landscape, then network state will move in the direction of the ‘Newton step’. In the quadratic case (elliptical level sets), this will directly ‘backtrack’ the synaptic fluctuations. Thus, the optimal magnitude of compensatory plasticity will be equal to that of the synaptic fluctuations. However, time delays in the sensing of synaptic fluctuations and limited precision of the compensatory plasticity mechanism will preclude this.

Figure 5 with 1 supplement see all
The dependence of steady state task performance in a nonlinear network on the magnitudes of compensatory plasticity and synaptic fluctuations, and on the learning rule quality.

Each (x,y) value on a given graph corresponds to an 8000 timepoint nonlinear neural network simulation (see ‘Methods’ for details). The y value gives the steady-state task error (average task error of the last 500 timepoints) of the simulation, while the x value gives the ratio of the magnitudes of the compensatory plasticity and synaptic fluctuations terms. Steady state error is averaged across 8 simulation repeats; shading depicts one standard deviation. Between graphs, we change simulation parameters. Down rows, we increase the proportionate noise corruption of the compensatory plasticity term (see Materials and methods section for details). Across columns, we increase the magnitude of synaptic fluctuations.

What about the case of a nonlinear loss function? Close to a minimum 𝐰*, the relationship of Equation (13) approximately holds (the loss function is locally quadratic). So if steady-state error is very low, we can directly transport the intuition of the quadratic case. However when steady state error increases, Equation (13) becomes increasingly approximate. In the limiting case, we could consider F[𝐰] as being completely uncorrelated from 2F[𝐰], in which case Q𝐰[F[𝐰]]Q𝐰[Δϵ] would hold. Numerical results in Figure 5 supports this assertion in nonlinear networks: the optimal ratio satisfies Δ𝐜2:Δϵ21 in conditions where steady-state task error is high, and Δ𝐜2:Δϵ21 in conditions where it is low.

Overall, we see that if Δ𝐜-F[𝐰] (i.e. compensatory plasticity enacts gradient descent), then we would expect compensatory plasticity to be outcompeted by synaptic fluctuations to maintain optimal steady-state error.

Even if compensatory plasticity does not move in the steepest direction of error decrease (i.e. the error gradient), it must move in an approximate downhill direction to improve task error (see e.g. Raman et al., 2019). Furthermore, the worse the quality of the gradient approximation, the larger the optimal level of compensatory plasticity (illustrated conceptually in Figure 4b–c, and numerically in Figure 5 and Figure 5—figure supplement 1). Why? We can rewrite such a learning rule as

Δ𝐜-F[𝐰]+ν,

where ν represents systematic error in the gradient approximation. The upward curvature in the direction Δ𝐜 becomes a (nonlinear) interpolation of the upward curvatures in the directions F[𝐰] and ν (see Equation (A6) of the SI). As long as ν is less biased towards high curvature directions than F[𝐰] itself, then this decreases the upward curvature in the direction Δ𝐜^, and thus increases the optimal compensatory plasticity rate. Indeed Figure 5 shows in simulation that this rate increases for more inaccurate compensatory plasticity mechanisms.

We now turn to zero-order learning algorithms, such as REINFORCE. These do not explicitly approximate a gradient, but generate random plasticity directions, which are retained/opposed based upon their observed effect on task error. We would expect randomly generated plasticity directions to have ‘average’ upward curvature, similarly to synaptic fluctuations. In this case, we would therefore get Q𝐰[Δ𝐜]Q𝐰[Δϵ], and compensatory plasticity should thus equal synaptic fluctuations in magnitude.

Finally, we consider second-order learning algorithms, and in particular the Newton update:

2F[𝐰]Δ𝐜=-F[𝐰].

As previously discussed, we assume that learning algorithms that require detailed information about the Hessian are biologically implausible. As such, our treatment is brief, and mainly contained in SI section 2.2.2.

In a linear network, the Newton update corresponds to compensatory plasticity making a direct ‘beeline’ for 𝐰* (see Figure 4d). As such Q𝐰[Δ𝐜]=Q𝐰[Δϵ] and the optimal magnitude of compensatory plasticity should match synaptic fluctuations. The same is true for a nonlinear network in a near-optimal state. However if steady-state task error is high in a nonlinear network, then compensatory plasticity should outcompete synaptic fluctuations. This case does not contradict our central claim however, since high task error at steady state implies that the task is not truly learned.

Together our results and analyses show that the magnitude of compensatory plasticity, at steady state task error, should be less or equal to that of synaptic fluctuations. This conclusion does not depend upon circuit architecture, or choice of biologically plausible learning algorithm.

Discussion

A long-standing question in neuroscience is how neural circuits maintain learned memories while being buffeted by synaptic fluctuations from noise and other task-independent processes (Fusi et al., 2005). There are several hypotheses that offer potential answers, none of which are mutually exclusive. One possibility is that fluctuations only occur in a subset of volatile connections that are relatively unimportant for learned behaviours (Moczulska et al., 2013; Chambers and Rumpel, 2017; Kasai et al., 2003). Following this line of thought, circuit models have been proposed that only require stability in a subset of synapses for stable function (Clopath et al., 2017; Mongillo et al., 2018; Susman et al., 2018).

Another hypothesis is that any memory degradation due to fluctuations is counteracted by restorative plasticity processes that allow circuits to continually ‘relearn’ stored associations. The information source directing this restorative plasticity could come from an external reinforcement signal (Kappel et al., 2018), from interactions with other circuits (Acker et al., 2018), or spontaneous, network-level reactivation events (Fauth and van Rossum, 2019). A final possibility is that ongoing synaptic fluctuations are accounted for by behavioural changes unrelated to learned task performance .

All these hypotheses share two core assumptions that we make, and several include a third that our results depend on:

  1. Not all synaptic changes are related to learning.

  2. Unchecked, these learning-independent plasticity sources generically hasten the degradation of previously stored information within a neural circuit.

  3. Some internal compensatory plasticity mechanism counteracts the degradation of previously stored information.

We extracted mathematical consequences of these three assumptions by building a general framework. We first modelled the the degree of degradation of previously learned information in terms of an abstract, scalar-valued, task error function or ‘loss landscape’. The brain may not have, and in any case does not require, explicit representation of such a function for a specific task. All that is required is error feedback from the environment and/or some internal prediction.

We then noted that compensatory plasticity should act to decrease task error, and thus point in a downhill direction on the ‘loss landscape’. We stress that we do not assume a gradient-based learning rule such as the backpropagation algorithm, the plausibility of which is an ongoing debate (Whittington and Bogacz, 2019).

Our results do not depend on whether synaptic changes during learning are gradual, or occur in large, abrupt steps. Although most theory work assumes plasticity to be gradual, there is evidence that plasticity can proceed in discrete jumps. For instance, abrupt potentiation of synaptic inputs that lead to the formation of place fields in mouse CA1 hippocampal neurons can occur within seconds as an animal explores a new environment (Bittner et al., 2017). Even classical plasticity paradigms that depend upon millisecond level precision in the relative timing of pre/post synaptic spikes follow a paradigm where there is a short ‘induction phase’ of a minute or so, following which there is a large and sustained change in synaptic efficacy (e.g. Markram et al., 1997; Bi and Poo, 1998). It is therefore an open question as to whether various forms of synaptic plasticity are best accounted for as an accumulation of small changes or a threshold phenomenon that results in a stepwise change. Our analysis is valid in either case. We quantify plasticity rate by picking a (large or small) time interval over which the net plasticity direction is approximately constant, and evaluate the optimal, steady-state magnitude of compensatory plasticity over this interval, relative to the magnitude of synaptic fluctuations.

A combination of learning-induced and learning-independent plasticity should lead to an eventual steady state level of task error, at which point the quality of stored information does not decay appreciably over time. The absolute quality of this steady state depends upon both the magnitude of the synaptic fluctuations, and the effectiveness of the compensatory plasticity.

Our main finding was that the quality of this steady state is optimal when the rate of compensatory plasticity does not outcompete that of the synaptic fluctuations. This result, which is purely mathematical in nature, is far from obvious. While it is intuitively clear that retention of circuit function will suffer when compensatory plasticity is absent or too weak, it is far less intuitive that the same is true generally when compensatory plasticity is too strong.

We also found that the precision of compensatory plasticity influenced its optimal rate. When ‘precision’ corresponds to the closeness of an approximation to a gradient-based compensatory plasticity rule, an increase in precision resulted in the optimal rate of compensatory plasticity being strictly less than that of fluctuations. In other words, sophisticated learning rules need to do less work to optimally overcome the damage done by learning-independent synaptic fluctuations. Indeed experimental estimates (see Table 1) suggest that activity-independent synaptic fluctuations can significantly outcompete systematic, activity-dependent changes in certain experimental contexts. Tentatively, this means that the high degree of synaptic turnover in these systems is in fact evidence for the operation of precise synaptic plasticity mechanisms as opposed to crude and imprecise mechanisms.

Our results are generic, in that they follow from fundamental mathematical relationships in optimisation theory, and hence are not dependent on particular circuit architectures or plasticity mechanisms. We considered cases in which synaptic fluctuations were distributed across an entire neural circuit. However, the basic framework easily extends, allowing for predictions in more specialised cases. For instance, recent theoretical work (Clopath et al., 2017; Mongillo et al., 2018; Susman et al., 2018) have hypothesised that synaptic fluctuations could be restricted to ‘unimportant’ synapses. These correspond to low curvature (globally insensitive) directions in the ‘loss landscape’. Our framework (Equation (9) in particular) immediately predicts that the optimal rate of compensatory plasticity will decrease proportionately with this curvature.

Precise experimental isolation/elimination of the plasticity sources attributable to learning and retention of memories remains challenging. Nevertheless, in conventional theories of learning (e.g. Hebbian learning), neural networks learn through plasticity induced by patterns of pre- and postsynaptic neural activity. A reasonable approximation, therefore, is to equate the ‘compensatory/learning-induced’ plasticity of our paper with ‘activity-dependent’ plasticity in experimental setups. With this assumption, our results provide several testable predictions.

Firstly, our results show that that the rate of compensatory (i.e. learning-dependent) plasticity is greater when a neural circuit is in a phase of active learning, as opposed to maintaining previously learned information (see Equation (10) and the surrounding discussion). Consequently, the relative contribution of synaptic fluctuations to the overall plasticity rate should be lower in this case. It would be interesting to test whether this were indeed the case, by comparing brain circuits in immature vs mature organisms, and in neural circuits thought to be actively learning vs those thought to be retaining previously learned information. One way to do this would be to measure the covariance of functional synaptic strengths at coinnervated synapses using EM reconstructions of neural tissue. A higher covariance implies a lower proportion of activity-dependent (i.e. compensatory) plasticity, since co-innervated synapses share presynaptic activity histories. Interestingly, two very similar experiments (Bartol et al., 2015) and (Dvorkin and Ziv, 2016) did indeed examine covariance in EM reconstructions of hippocampus and neocortex, respectively. This covariance appears to be much lower in hippocampus (compare Figure 1 of Bartol et al., 2015 to Figure 8 of Dvorkin and Ziv, 2016). Many cognitive theories characterise hippocampus as a continual learner and neocortex as a consolidator of previously learned information (e.g. O'Reilly and Rudy, 2001). Our analysis provides support for this hypothesis at a mechanistic level by linking low covariance in coinnervated hippocampal synapses to continual learning.

Secondly, a number of experimental studies (Nagaoka et al., 2016; Quinn et al., 2019; Yasumatsu et al., 2008; Minerbi et al., 2009; Dvorkin and Ziv, 2016) note a persistence of the bulk of synaptic plasticity in the absence of activity-dependent plasticity or other correlates of an explicit learning signal, as explained in our review of key experimental findings. However, there are two important caveats for relating our work to these experimental observations:

  • Experimentally isolating different plasticity mechanisms, measuring synaptic changes, and accounting for confounding behavioural/physiological changes is extremely challenging. The most compelling in vivo support comes from Nagaoka et al., 2016, where an analogue of compensatory plasticity in the mouse visual cortex was suppressed both chemically (by suppression of spiking activity) and behaviourally (by raising the mouse in visually impoverished conditions). Synaptic turnover was reduced by about half for both suppression protocols, and also when they were applied simultaneously. Further studies that quantified changes in synaptic strength in addition to spine turnover in an analogous setup would lend further credence to our results.

  • We do not know if observed synaptic plasticity in the experiments we cite truly reflect a neural circuit attempting to minimise steady-state error on a particular learning goal (as captured through an abstract, implicit, ‘loss function’). Our analysis simply shows that somewhat surprising levels of ongoing plasticity can be explained parsimoniously in such a framework. In particular, the concepts of ‘learning’ and behaviour have no clear relationship with neural circuit dynamics in vitro. Nevertheless, we might speculate that synapses could tune the extent to which they respond to ‘endogenous’ (task independent) signals versus external signals that could convey task information in the intact animal. Even if the information conveyed by activity-dependent signals were disrupted in vitro, the fact that activity-dependent signals determined such a small proportion of plasticity is notable, and seems to carry over to the in vivo case.

Thus, while our results offer a surprising agreement with a number of experimental observations, we believe it is important to further replicate measurements of synaptic modification in a variety of settings, both in vivo and in vitro. We hope our analysis provides an impetus for this difficult experimental work by offering a first-principles theory for the volatility of connections in neural circuits.

Materials and methods

Simulations

Request a detailed protocol

We simulated two types of network, which we refer to as linear (Figure 5—figure supplement 1) and nonlinear (Figures 1 and 5) respectively. We ran our simulations in the Julia programming language (version 1.3), and in particular used the Flux.jl software package (version 0.9) to construct and update networks. Source code is available at https://github.com/Dhruva2/OptimalPlasticityRatios (copy archived at swh:1:rev:fcb1717a822f90b733c49d62bfc2f970155b7364, Raman, 2021).

Nonlinear networks

Request a detailed protocol

Networks were rate-based, with the firing rate r(t) of a given neuron defined as

r(t)=σ(wT(t)u(t)),

where w is the vector of presynaptic strengths, u represents the firing rate of the associated presynaptic neurons, and σ(x):=11+exp(-x) is the sigmoid function. Initial weight values were generated randomly, according to the standard Xavier distribution (Glorot and Bengio, 2010). Networks were organised into three layers, containing 12, 20, and 10 neurons, respectively. Any given neuron was connected to all neurons in the previous layer. For the first layer, the firing rates of the ‘previous layer’ corresponded to the network inputs.

Linear networks

Request a detailed protocol

Networks were organised into an input layer of 12 neurons, and an output layer of 10 neurons. Each output neuron was connected to all input layer neurons. Networks were rate-based, with the firing rate r(t) of a given neuron defined as

r(t)=wT(t)u(t),

where ui(t) corresponds to the ith input (input-layer neuron) or the firing rate of the ith input-layer neuron (output-layer neuron). Initial weight values were generated randomly, according to the Xavier distribution (Glorot and Bengio, 2010).

Task error

Request a detailed protocol

For each network, we generated 1000 different, random, input vectors. Each component of the vector was generated from a unit Gaussian distribution. Task error, at the tth timestep, was taken as the mean squared error of the network in recreating the outputs of the initial (t=0) network, in response to the suite of inputs. Mathematically, this equates to

F[w(t)]=1|𝒰|u𝒰y(w(t),u)y(w(0),u)22,

where y(𝐰(t),u) denotes the output of the network given the synaptic strengths at time t, in response to input u𝒰. Note that this task error recreates the ‘student-teacher’ framework of e.g. (Levin et al., 1990; Seung et al., 1992), where a fixed copy of the initial network is the teacher.

Weight dynamics

Request a detailed protocol

At each simulation timestep, synaptic weights were updated as

Δwt+1=Δct+Δϵt.

We took the synaptic fluctuations term, Δϵt, as scaled white noise, that is,

Δϵt𝒩(0,𝕀)

The constant of proportionality was calculated so that the magnitude Δϵ2 conformed to a pre-specified value. This magnitude was 2 in the simulation of Figure 1, and was a graphed variable in the simulations of Figure 5 and Figure 5—figure supplement 1.

The compensatory plasticity term, Δ𝐜t, was calculated in two stages. First we applied the backpropagation algorithm, using y(𝐰(0),u) as the ideal network outputs to train against. This generated an ‘ideal’ direction of compensatory plasticity , proportional to the negative gradient F[𝐰(t)]. For Figure 5 and Figure 5—figure supplement 1 we then corrupted this gradient with a tunable proportion of white noise. Overall, this gives,

Δ𝐜t=-γ1F^[𝐰]t+γ2ν^t,

where νt𝒩(0,𝕀) is the noise corruption term, and γ1,γ2>0 are tunable hyperparameters. The higher the ratio γ2:γ1, the greater the noise corruption. Meanwhile, γ12+γ22 sets the overall magnitude of compensatory plasticity . By tuning γ1 and γ2, we can therefore independently modify the magnitude and precision of the compensatory plasticity term. In Figure 1, we set γ2=0.

Appendix 1

Alternative derivation of Equation (9a)

We provide an alternative derivation of Equation (9a) that removes the need for assumption (3b). We did not put this main derivation in the main text as we perceive it to have less clarity.

The derivation proceeds identically to that given in the main text until Equation (7). We can then use (3a) to simplify Equation (7). We get

E[ΔF]=ΔcTF[w(t)]+12Δc22Qw(t)[Δc]+12Δϵ22Qw(t)[Δϵ]+ΔcT(2F[w(t)])Δϵ.
E[ΔF]=ΔcTF[w(t)]+12ΔcT(2F[w(t)])Δc+12ΔϵT(2F[w(t)])Δϵ+ΔcT(2F[w(t)])Δϵ.

Recall that expectation is taken over an unknown probability distribution from which Δϵ is drawn, which satisfies Equation (3a).

We then assume that we are in a phase of stable memory retention, so that 𝔼[ΔF]=0. Now if the magnitude of compensatory plasticity Δ𝐜2 is tuned to minimise steady state error F, then any change to Δ𝐜2 will result in an increase in 𝔼[ΔF]. So 𝔼[ΔF] is locally minimal in Δ𝐜2. This implies

dE[ΔF]dΔc2=0.

We also claim that local minimality implies

(1) dE[ΔFΔc2]dΔc2=0.

Why? 𝔼[ΔF]=0 implies that 𝔼[ΔFΔ𝐜2]=0. If a small change to Δ𝐜2 results in 𝔼[ΔF]0, then it also results in 𝔼[ΔFΔ𝐜2]0, since Δ𝐜2 is non-negative.

Expanding the LHS of Equation (1), we get

ddΔc2{Δc^TF[w(t)]+12Δc2Δc^T(2F[w(t)])Δc^+12ΔϵT(2F[w(t)])ΔϵΔc2+Δc^T(2F[w(t)])Δϵ}=0.

Differentiating, we get

12Δc^T(2F[w(t)])Δc^=12ΔϵT(2F[w(t)])ΔϵΔc22.
Qw[Δc]=Δϵ22Δc22Qw[Δϵ],

from which (9a) follows.

Positivity of the numerator and denominator in Equation (9a)

Equation (9a) of the main text asserts that

Δc22Δϵ22=Qw[Δϵ]Qw[Δc]

holds as long as both the numerator and denominator of the RHS are positive. Here we describe sufficient conditions for positivity.

The inequality 2F[𝐰]0 must hold in some neighbourhood of any minimum 𝐰*. Recall that we referred to such a neighbourhood as a highly trained state of the network in the main text. In such a state, our assertion follows immediately, as Q𝐰[𝐯]:=1𝐯22𝐯T(2F[𝐰])𝐯0, for any vector 𝐯. Therefore, Q𝐰[Δϵ]0 and Q𝐰[Δ𝐜]0.

We now consider a partially trained network state, which we defined in the main text as any 𝐰 satisfying Tr(2F)0. Note that

F[w+Δϵ]=F[w]+F[w]TΔϵ+12ΔϵTF[w]Δϵ+𝒪(Δϵ23).

We assumed in the main text (Equation (3a)), that Δϵ is uncorrelated with the gradient F[𝐰] in expectation, since Δϵ is realised by memory-independent processes. Similarly we can assume that Δϵ is unbiased in how it projects onto the eigenvectors of 2F[𝐰]. In other words,

E[v^iTΔϵ]=E[v^jTΔϵ],

for any normalised eigenvectors v^i, v^j of 2F[w]. In expectation, we can therefore simplify to

E[F[w+Δϵ]]=F[w]+E[F[w]TΔϵ+12ΔϵTF[w]Δϵ]+𝒪(Δϵ23).
=0+Δϵ22Tr(2F[w])N+𝒪(Δϵ23),

where N is the dimensionality of the vector 𝐰. So a partially trained network is one for which small, memory-independent weight fluctuations (such as Δϵ, or white noise) are expected to decrease task performance.

Now recall that Q𝐰[Δϵ]=1Δϵ22ΔϵT2F[𝐰]Δϵ. So we have

E[Qw[Δϵ]]=Tr(2F[w])N>0,

where the positivity constraint comes from being in a partially trained network.

We now consider why Q𝐰[Δ𝐜] should be generically positive in a partially trained network. Suppose Q𝐰[Δ𝐜]<0 holds. We can rewrite this as Δ𝐜T2F[𝐰]Δ𝐜0. In this case, maintaining the same compensatory plasticity Δ𝐜 over the time interval [t*+Δt,t*+2Δt] would result in increased improvement in loss, as

F[w+Δc]TΔc=F[w]TΔc+ΔcT2F[w]Δc+𝒪(Δc22).

Effectively, memory improvement due to compensatory plasticity Δ𝐜 would be in an ‘accelerating’ direction, and maintaining the same direction Δ𝐜 of compensatory plasticity would lead to ever faster learning. However, by assumption, we are in a regime of steady state task performance, where

E[F(t+2Δt)F(t+Δt)]=E[F(t+Δt)F(t)]=0.

Optimal plasticity ratios in specific learning rules

Noise-free learning rules (first-order)

Let us first consider the case where Δ𝐜 can be computed with perfect access to the gradient F[𝐰], but without access to 2F[𝐰]. Such a Δ𝐜 is known as a first-order learning rule, as it has access only to the first derivative of F (Polyak, 1987). Imperfect access is considered subsequently. In this case, the optimal direction of compensatory plasticity is

Δ𝐜-F[𝐰].

In other words, Δ𝐜 would implement perfect gradient descent on F[𝐰]. The condition of Equation (11) for synaptic fluctuations to outcompete reconsolidation plasticity evaluates to

Q𝐰[F[𝐰]]Q𝐰[Δϵ].

To what extent can we quantify Q𝐰[F[𝐰]]? First let us relate the gradient and Hessian of F[𝐰]. Let 𝐰* be an optimal state of the network (i.e. one where F is minimised). Let us parameterise the straight line connecting 𝐰 with 𝐰*:

γ(s)=s𝐰*+(1-s)𝐰,s[0,1].

Then

F[w]=(ww)TM, where 
M=012F[γ(s)] ds.

This gives

Qw[F[w]]=(ww)TMT2F[w]M(ww)(ww)TMTM(ww).

First let us rewrite

(𝐰-𝐰*):=iNcivi,
M(𝐰-𝐰*):=iNdivi

where (λi,vi) is the ith eigenvalue/eigenvector pair of 2F (sorted in ascending order of λi), and ci, di are some scalar weights. Now

(2) Qw[F[w]]=i=1Ndi2λii=1Ndi2.

The value of Q𝐰[F[𝐰]] now depends upon the distribution of mass of the sequence {di}. If later elements of the sequence are larger (i.e. M(𝐰=𝐰*) projects more highly onto eigenvectors of 2F[w] with large eigenvalue), then Q𝐰[F[𝐰]] becomes larger, and the optimal magnitude of reconsolidation plasticity decreases, relative to the magnitude of synaptic fluctuations. The opposite is true if earlier elements of the sequence are larger.

Guaranteed bounds on the value of Equation (2) are vacuous. If we do not restrict M, then we can tailor the sequence {di} as we like, and we end up with λ1Q𝐰[F[𝐰]]λN. However, pragmatic bounds are much tighter. Let us now consider two plausibly extremal cases.

First consider the simplest case of a network that linearly transforms its outputs, and which has a quadratic loss function F[𝐰]. In this case 2F is a constant, (independent of 𝐰), positive-semidefinite matrix, and M=2F. This means that

di=ciλivi
Qw[F[w]]=i=1Nci2λi3i=1Nci2λi2.

Condition (11) then becomes

(3) Qw[F[w]]Qw[Δϵ]i=1Nci2λi3i=1Nci2λi2i=1NλiN.

A conservative sufficient condition for (Equation 3), using Chebyshev’s summation inequality, is that

(4) ci2λi2ci12λi12, for all i{1,,N}.

Under what conditions would a plausible reconsolidation mechanism choose to ‘outcompete’ synaptic fluctuations, in this linear example? For Q𝐰[F[𝐰]]<Q𝐰[Δϵ] to even hold, (26) would have to be broken, and significantly so due to conservatism in the inequality. In other words, 𝐰-𝐰* must project quite biasedly onto the eigenvectors of 2F with smaller-than-average eigenvalue. If the discrepancy between 𝐰 and 𝐰* were caused by fluctuations (which are independent of 2F), then this would not be the case, in expectation. Even if this were the case, the reconsolidation mechanism would have to know about the described bias. This requires knowledge of both 𝐰* and 2F, and is thus implausible.

Now let us consider the case of a generic nonlinear network. At one extreme, if 𝐰-𝐰*2 is small, then M2F[𝐰], and the discussion of the linear case is valid. This corresponds to the case where steady state error is close to the minimum achievable by the network. As 𝐰-𝐰*2 increases (i.e. steady state error gets worse), the correspondence between M and 2F[𝐰] will likely decrease. Thus the optimal magnitude of reconsolidation plasticity, relative to the level of synaptic fluctuations, will rise.

We could consider another ‘extreme’ case in which M and 2F[𝐰] were completely independent of each other. In this case,

(5) di21Ni=1Ndi2.

In other words, the projection of M(𝐰-𝐰*) onto the different eigenvectors of 2F[𝐰] is approximately even. Using (24), this gives

Qw[F[w]]i=1NλiN=Qw[Δϵ].

In summary, we have two plausible extremes. One occurs where M=2F[𝐰], and another occurs where M is completely independent of 2F[𝐰]. In either case, Q𝐰[F[𝐰]]Q𝐰[Δϵ], and so the magnitude of synaptic fluctuations should optimally outcompete/equal the magnitude of reconsolidation plasticity. Of course, there might be particular values of 𝐰 where the correspondence between M and 2F[𝐰] is ‘worse’ than chance. In other words, eigenvectors of M with large eigenvalue preferentially project onto eigenvectors of 2F[𝐰] with small eigenvalue. In such cases, we would have Q𝐰[F[𝐰]]Q𝐰[Δϵ]. However, we find it implausible that a reconsolidation mechanism would be able to gain sufficient information on M to determine this at particular points in time, and thereby increase its plasticity magnitude.

Noise-free learning rules (second-order)

Let us now suppose that Δ𝐜 can be computed with perfect access to both F[𝐰] and 2F[𝐰]. In this case, the reconsolidation mechanism would optimally apply plasticity in the direction of the Newton step: we would have

2F[w]Δc=F[w].

Note that the Newton step is often conceptualised as a weighted form of gradient descent, where movement on the loss landscape is biased towards direction of lower curvature. Thus we would expect Q𝐰[Δ𝐜] to be smaller, and the optimal proportion of reconsolidation plasticity to be larger. This is indeed the case. For mathematical tractability, we will restrict our discussion to the case in which 2F[𝐰]0, and M0. This would hold if F[𝐰] were convex, or if 𝐰 were sufficiently close to a unique local minimum 𝐰*. In this case we can rewrite

Δc=2F[w]1F[w],

which gives

(6a)Qw[Δc]=F[w]T(2F[w])1F[w]F[w]T(2F[w])2F[w](6b)=(ww)M(2F[w])1M(ww)(ww)M(2F[w])2M(ww)(6c)=i=1Ndi2λi1i=1Ndi2λi2.

Once again, we first consider the case of a linear network with quadratic loss function, and hence with constant Hessian 2F. This gives M=2F, and

Qw[Δc]=(ww)2F[w](ww)ww22
=i=1Nci2λii=1Nci2.

We again assume that the reconsolidation mechanism does not have knowledge of the relative projections of 𝐰-𝐰* onto the different eigenvectors of 2F, which requires knowledge of 𝐰*. Without such information, we can use an analogous argument to that preceding (Equation 5) to argue that the approximation ci21Ni=1Nci2 is reasonable. This gives Q𝐰[Δ𝐜]Q𝐰[Δϵ].

Note that the Newton step, in the linear-quadratic case just considered, corresponds to a direction 𝐰*-𝐰, that is, a direct path to a local minimum. So we could consider a compensatory plasticity mechanism implementing the Newton step as one directly undoing synaptic changes caused by Δϵ.

We now consider the case of a nonlinear network. As before, if 𝐰-𝐰*2 is small, then we have M2F[𝐰], and the arguments of the linear network hold. As 𝐰-𝐰*2 increases, the correspondence between M and 2F will decrease. We again consider the plausible extreme where M is completely uncorrelated with 2F[𝐰], and so the approximation (Equation 5) holds. In this case, Equation (6) can be simplified to give

Qw[Δc]i=1Nλi1i=1Nλi2.

We assumed that 2F[𝐰]0. Therefore, all eigenvalues are positive. This allows us to use Chebyshev’s summation inequality to arrive at

i=1Nλi1i=1Nλi2i=1NλiN=Qw[Δϵ].

So as 𝐰-𝐰*2 increases, the magnitude of reconsolidation plasticity will optimally outcompete that of synaptic fluctuations. This is the one case that contradicts our main claim.

Imperfect learning rules

The previous section applied in the implausible case where a reconsolidation mechanism had perfect access to F[𝐰] and/or 2F[𝐰]. Recall from the main text that at least some information on F[𝐰] is required, in order for compensatory plasticity to move in a direction of decreasing task error. What if Δ𝐜 contains a mean-zero noise term, corresponding to unbiased noise corruption of these quantities? We will now show how such noise pushes Q𝐰[Δ𝐜] towards equality with Q𝐰[Δϵ], and thus pushes the optimal magnitude of reconsolidation plasticity towards the magnitude of synaptic fluctuations. Let us use the model

(7) Δc=Δc~+ν,

where ν is some mean-zero random variable, and Δ𝐜~ is the ideal output of the reconsolidation mechanism, assuming perfect access to the derivatives of F[𝐰]. Here ν represents the portion of compensatory plasticity attributable to systematic error in the algorithm, due to imperfect information on F[𝐰]. This could arise due to imperfect sensory information or limited communication between synapses. We can therefore assume, as for Δϵ, that it does not contain information on 2F[𝐰]. We therefore get

Qw[ν]Tr(2F[w])N,

analogously to Equation (12). Now the operator Q𝐰 satisfies

(8) Qw[Δc]=Qw[Δc~](1+ν22Δc~22)1+Qw[ν](1+Δc~22ν22)1.

So depending upon the relative magnitudes of Δ𝐜~ and ν, Q𝐰[Δ𝐜] interpolates between Q𝐰[Δ𝐜~] and Q𝐰[ν]. In particular, as the crudeness of the learning rule (i.e. the ratio νΔ𝐜~ ) grows, Q𝐰[Δ𝐜] approaches equality (from below) with Q𝐰[ν], and thus Q𝐰[Δϵ], completing our argument.

Data availability

All code is publicly available on github at this URL: https://github.com/Dhruva2/OptimalPlasticityRatios (copy archived at https://archive.softwareheritage.org/swh:1:rev:fcb1717a822f90b733c49d62bfc2f970155b7364).

References

    1. Bi GQ
    2. Poo MM
    (1998)
    Synaptic modifications in cultured hippocampal neurons: dependence on spike timing, synaptic strength, and postsynaptic cell type
    The Journal of Neuroscience 18:10464–10472.
  1. Conference
    1. Glorot X
    2. Bengio Y
    (2010)
    Understanding the difficulty of training deep feedforward neural networks
    Proceedings of the Thirteenth International Conference on Artificial Intelligence and Statistics. pp. 249–256.
  2. Book
    1. Hebb DO
    (1949)
    The Organization of Behavior: A Neuropsychological Theory
    Wiley: Chapman & Hall.
  3. Book
    1. Polyak BT
    (1987)
    Introduction to Optimization
    New York: Optimization Software, Publications Division.

Article and author information

Author details

  1. Dhruva V Raman

    Department of Engineering, University of Cambridge, Cambridge, United Kingdom
    Contribution
    Conceptualization, Formal analysis, Investigation, Visualization, Methodology, Writing - original draft, Writing - review and editing
    For correspondence
    dvr23@cam.ac.uk
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-8992-1353
  2. Timothy O'Leary

    Department of Engineering, University of Cambridge, Cambridge, United Kingdom
    Contribution
    Conceptualization, Supervision, Funding acquisition, Visualization, Writing - original draft, Project administration, Writing - review and editing
    For correspondence
    tso24@cam.ac.uk
    Competing interests
    Reviewing editor, eLife
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-1029-0158

Funding

European Commission (StG 2016 716643 FLEXNEURO)

  • Dhruva V Raman
  • Timothy O'Leary

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

This work was supported by ERC grant StG 2016 716643 FLEXNEURO.

Copyright

© 2021, Raman and O'Leary

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 1,784
    views
  • 335
    downloads
  • 17
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Dhruva V Raman
  2. Timothy O'Leary
(2021)
Optimal plasticity for memory maintenance during ongoing synaptic change
eLife 10:e62912.
https://doi.org/10.7554/eLife.62912

Share this article

https://doi.org/10.7554/eLife.62912

Further reading

    1. Computational and Systems Biology
    2. Developmental Biology
    Juan Manuel Gomez, Hendrik Nolte ... Maria Leptin
    Research Article Updated

    The initially homogeneous epithelium of the early Drosophila embryo differentiates into regional subpopulations with different behaviours and physical properties that are needed for morphogenesis. The factors at top of the genetic hierarchy that control these behaviours are known, but many of their targets are not. To understand how proteins work together to mediate differential cellular activities, we studied in an unbiased manner the proteomes and phosphoproteomes of the three main cell populations along the dorso-ventral axis during gastrulation using mutant embryos that represent the different populations. We detected 6111 protein groups and 6259 phosphosites of which 3398 and 3433 were differentially regulated, respectively. The changes in phosphosite abundance did not correlate with changes in host protein abundance, showing phosphorylation to be a regulatory step during gastrulation. Hierarchical clustering of protein groups and phosphosites identified clusters that contain known fate determinants such as Doc1, Sog, Snail, and Twist. The recovery of the appropriate known marker proteins in each of the different mutants we used validated the approach, but also revealed that two mutations that both interfere with the dorsal fate pathway, Toll10B and serpin27aex do this in very different manners. Diffused network analyses within each cluster point to microtubule components as one of the main groups of regulated proteins. Functional studies on the role of microtubules provide the proof of principle that microtubules have different functions in different domains along the DV axis of the embryo.

    1. Computational and Systems Biology
    2. Physics of Living Systems
    Natanael Spisak, Gabriel Athènes ... Aleksandra M Walczak
    Tools and Resources Updated

    B-cell repertoires are characterized by a diverse set of receptors of distinct specificities generated through two processes of somatic diversification: V(D)J recombination and somatic hypermutations. B-cell clonal families stem from the same V(D)J recombination event, but differ in their hypermutations. Clonal families identification is key to understanding B-cell repertoire function, evolution, and dynamics. We present HILARy (high-precision inference of lineages in antibody repertoires), an efficient, fast, and precise method to identify clonal families from single- or paired-chain repertoire sequencing datasets. HILARy combines probabilistic models that capture the receptor generation and selection statistics with adapted clustering methods to achieve consistently high inference accuracy. It automatically leverages the phylogenetic signal of shared mutations in difficult repertoire subsets. Exploiting the high sensitivity of the method, we find the statistics of evolutionary properties such as the site frequency spectrum and dN/dS ratio do not depend on the junction length. We also identify a broad range of selection pressures spanning two orders of magnitude.