extfieldsofstudy: Chemistry Biology Physics Medicine

text:
Origins of concentration dependence of waiting times for single-molecule fluorescence binding

Binary fluorescence time series obtained from single-molecule imaging experiments can be used to infer protein binding kinetics, in particular, association and dissociation rate constants from waiting time statistics of fluorescence intensity changes. In many cases, rate constants inferred from fluorescence time series exhibit nonintuitive dependence on ligand concentration. Here we examine several possible mechanistic and technical origins that may induce ligand dependence of rate constants. Using aggregated Markov models, we show under the condition of detailed balance that non-fluorescent bindings and missed events due to transient interactions, instead of conformation fluctuations, may underly the dependence of waiting times and thus apparent rate constants on ligand concentrations. In general, waiting times are rational functions of ligand concentration. The shape of concentration dependence is qualitatively affected by the number of binding sites in the single molecule and is quantitatively tuned by model parameters. We also show that ligand dependence can be caused by non-equilibrium conditions which result in violations of detailed balance and require an energy source. As to a different but significant mechanism, we examine the effect of ambient buffers that can substantially reduce the effective concentration of ligands that interact with the single molecules. To demonstrate the effects by these mechanisms, we applied our results to analyze the concentration dependence in a single-molecule experiment EGFR binding to fluorophore-labeled adaptor protein Grb2 by Morimatsu et al. (PNAS,104:18013,2007).

I. INTRODUCTION
Single-molecule fluorescence techniques measure realtime kinetics of chemical reactions, such as RNA folding 2 , enzymatic reactions 3,4 , protein-protein, proteinoligonucleotides and protein-DNA bindings 1,[5][6][7][8][9] . An important application of single-molecule fluorescence techniques is to probe the binding kinetics between interaction partners and to infer binding rate constants from the fluorescence time series. In contrast to conventional ensemble-averaged measurements, single-molecule fluorescence techniques directly observes the binding stochasticity at the molecular level, allowing investigations of conformation fluctuations of individual molecules under various conditions.
Single-molecule fluorescence time series often exhibit transitions alternating between two observable states: fluorescent (on) and non-fluorescent (off) states. One can summarize information in such "binary" time series using one-dimensional waiting time distributions and kinetic rate constants for association (k on ) and dissociation (k off ) can therefore be derived from the mean waiting times. Common experimental procedures usually involve measuring single-molecule binding fluorescence time series under varied ligand concentrations.
These waiting time distributions for protein binding in many cases are fit by sums of multiple exponentials a) Electronic mail: jinyang2004@gmail.com or empirically by stretched exponentials 1,4,10 , suggesting that conformation of the single molecule fluctuates during the course of interaction. Therefore, transitions between the two macroscopic states ("on" and "off") proceed through diverse conformation channels connecting microscopic states which cannot be directly distinguished by fluorescence time series. Such temporal variants in protein conformation are referred to as dynamic disorder if conformations fluctuate on a time scale comparable to that of protein binding 3,11 .
To analyze a binary fluorescence time series, the kinetics of binding between a single molecule and its ligand is usually described by the simple phenomenological twostate model: where state "off" indicates that the molecule is free and state "on" indicates that the molecule is bound to a ligand, and the transition rate from "off" to "on" is proportional to the ligand concentration [L] due to the law of mass action. One can obtain the apparent rate constants k on and k off from the mean waiting times, τ on and τ off , respectively as: where τ off is the mean waiting time for association (or, the mean dwell time at the macroscopic "off" state) and τ on is the mean waiting time for dissociation (or, the mean dwell time at the macroscopic "on" state). The apparent dissociation constant is then given as K d = k off /k on . In this two-state model, when forward and backward transitions are Markovian, in other words, if the two-state model is biochemically elementary, the rate constants k on and k off are independent of [L]. However, measured single-molecule kinetics may potentially be affected by various mechanistic and technical factors including protein concentrations, protein conformations and cellular environments, etc. These factors may cause significant ligand concentration dependence of kinetic rate constants inferred by the mean waiting times as calculated by Eq. [2].
In particular, transitions between observed "on" and "off" states can be complex due to conformation fluctuations, such that the two-state minimalistic model of Eq. [1] is non-Markovian and thus is not adequate to describe the true transition mechanism. For such cases, mechanistic models are required to analyze the data. Here, we examine several possible mechanisms that might cause ligand dependence of the kinetic rate constants k on and k off . Conformation fluctuations in single molecules are described by aggregated Markov models as usually being treated in analysis of single-ion channel recordings [12][13][14][15] . In addition, we show that non-equilibrium models that violate the detailed balance constraint can also generate strong ligand concentration dependence of rate constants.
We further examine the ligand dependence caused by two technical sources that are not directly related to the binding mechanism: (1) missed events, and (2) background buffer. The former is concerned with transient events that are not captured by observations, which causes overestimation and distortion of ligand dependence in waiting times. We show that missed events have an effect on waiting times similar to that by models of single protein binding site with non-fluorescent binding. The latter mechanism is concerned with ligand interactions with ambient buffer molecules, which causes a substantial reduction in the effective concentration of ligands that bind to the single molecules. Such a background buffering mechanism if unaccounted will cause a strong dependence of the apparent association rate constant on the total ligand concentration [L].
We apply our models to analyze an in vitro experiment of single epidermal growth factor receptor (EGFR) binding to adaptor protein Grb2 by Morimatsu et al.s 1,7 . The experiment showed that waiting time distributions had multiple exponential decay, suggesting that EGFR molecule may have conformational changes on time scales comparable to Grb2 binding. Moreover, the apparent association rate constant k on had a counter-intuitive dependence on Grb2 concentrations.

A. Aggregated Markov models
The theory of aggregated Markov models was developed for analyzing ion channel gating mechanisms from single ion channel recordings 12,13 . Similar Markov models have been developed more recently in analyzing times series obtained from single-molecule fluorescence experiments 16 . An aggregated Markov model is a special case of hidden Markov chain, in which states in a particular category (an aggregate) correspond to a signal of an identical fluorescence intensity. In a binary single-molecule model, the system fluoresces in states in the "on" aggregate and does not fluoresce in states in the "off" aggregate. An aggregated Markov model for a single molecule can be fully characterized in terms of the "generator matrix", Q, which has an off-diagonal structure that encodes the reaction scheme. Entry q ij is the transition rate from state i to state j. The diagonal entries are defined as q ii = − j q ij . In a matrix form, one can write Qu = 0, where u is the right null vector of all ones. The steady-state occupancies w is the normalized left null vector of Q, i.e., wQ = 0 and i w i = 1. For systems with aggregates "on" and "off", Q can be organized and partitioned as where diagonal blocks contain intra-aggregate transition rates and off-diagonal blocks contain inter-aggregate transition rates. o and c in Eq. [3] denote "on" and "off" aggregates, respectively. We can correspondingly partition the null vectors: w = [w o w c ] and u = [u o u c ] T . The mean waiting times are given as each of which is the ratio of steady-state aggregate occupancy, P on ≡ w o u o or P off ≡ w c u c , to the total interaggregate probability flux (e.g., J oc ≡ w o Q oc u c ). The inter-aggregate probability fluxes are balanced at the steady state, i.e., J oc = J co .

B. Detailed balance
The principle of microscopic reversibility (or the law of detailed balance) 17 states that at thermodynamic equilibrium for any reversible transition between two neighboring states i and j the probability flux from state i to state j is balanced by that from j to i, i.e., w i q ij = w j q ji . In general, the occupancy w i has a complex relationship with the ligand concentration [L]. However, under the detailed balance, one can derive a relative occupancy w, in which each entry has a monomial dependence on [L]. This treatment eases analysis of ligand dependence of k on and k off . By the law of mass action, the rate of a ligand FIG. 1. Template I consists of an "on" aggregate in which the receptor has a ligand bound, and an "off" aggregate with no ligand bound. Models constructed from Template I may have any number of "on" and "off" states with any connectivity. Rates for transitions from "off" states (empty circles) to "on" states (filled circles) are proportional to the ligand concentration [L]. All reversible rates for transitions from "on" states to "off" states are constant, and all intra-aggregate transitions are spontaneous and have constant transition rates. Template II extends Template I with unresolved ligand-bound "off" aggregate (gray circles). Template III describes a single molecule with two binding sites. One site is fluorescent (dark dot) when binding to a ligand, whereas the other is nonfluorescent (gray dot). Again there can be any number of states in each of the four aggregates of Template III.
binding is proportional to [L], whereas the reversible dissociation rate is a constant. With the detailed balance condition, for two neighboring states i and j, we have where b = −1 if the state transition from i to j is induced by a ligand binding, b = 1 if the state transition from j to i is induced by a ligand binding, and b = 0 if both forward and backward transitions do not involve ligand binding. The coefficient K ij is a constant. Designate a reference state r at which the molecule does not bind to a ligand. We define the relative occupancy as the ratio w i ≡ w i /w r , and note that state i and r are connected by a path involving one or more transitions. Applying Eq. [5] successively along a path in the reaction scheme from state i to state r, we can show w i = k i [L] ni , where the non-negative integer n i is the number of ligands bound at state i and k i is the product of equilibrium constants for all the reversible reactions along the path from state i to state r. The numerical value of k i depends on the choice of the reference state r but does not depend on the choice of a path that connects state i and r due to detailed balance. The mean waiting times are now given as From the above derivation, it is evident that in general both τ on and τ off are ligand dependent as rational functions of [L].

C. Model templates
Instead of using models with specific mechanisms, we discuss the apparent rate constants under different model classes depicted in Fig. 1 by "template" reaction schemes. A template consists of several different aggregates. Template I has two aggregates, an "on" aggregate which has ligand bound and an "off" aggregate with no ligand bound. Models constructed from Templates I and II can have any number of on and off states and any connectivity. For Template III, the connectivity is arbitrary except that no state in the unliganded aggregate can be directly linked to states in the doubly liganded aggregate. We consider models with two experimentally-distinguishable fluorescent ("on") and non-fluorescent ("off") aggregates with three categories of states for the single molecule; (i) a dark state ("off") with no ligand bound, (ii) a ligandbound dark state ("off"), and (iii) a bright state ("on") which has ligand bound. The advantage of this approach is that our conclusions do not depend on model details such as the number of states and the connectivity between states but only on key biochemical aspects such as the number of binding sites and whether the binding protein fluoresces while associating with the receptor. Although a waiting time distribution that is a sum of multiple exponents depends on the number of states, the ligand dependence of the apparent association rate does not.

A. Conformation fluctuations do not cause concentration-dependence of rate constants
We first consider a single molecule that has a single ligand binding site. We assume that every ligand binding event is experimentally observed, which therefore switches the molecule from an "off" state to an "on" state. Each ligand departure switches the molecule from an "on" state to an "off" state. Intra-aggregate state transitions do not involve ligand arrival or departure. Clearly, any such kinetic model can be constructed from a two-state base scheme as shown in Template I (Fig. 1). Conformation fluctuations could be modeled by extending the base scheme with multiple "on" and "off" states with an arbitrary connections between states. We can write w o = k o [L] and w c = k c , where k o and k c are constant vectors. The mean waiting times are given as Since Q oc only contains ligand dissociation rate constants τ on is a constant and τ off is proportional to the inverse of [L]. Thus, we have shown that conformation fluctuations do not generate ligand concentration dependence of the apparent rate constants, k on and k off as defined in Eq. [2]. This result holds for a reaction scheme with arbitrary number of states and arbitrary connections between states as extended from Template I, as long as each binding event was directly observed by the experiment.
Notice that in such case k on and k off calculated from the mean waiting times are in fact identical to those obtained in ensemble-averaged measurements.

B. Effect of non-fluorescent binding
In some cases, ligand binding may not be resolved experimentally and become unnoticed, which could cause ligand dependence as we show below. We consider a model that incorporates non-fluorescent ligand binding to the single molecule. A model (Template II, Fig. 1) of a molecule that has a single binding site allows a "dark" conformation (a "c" state) in which a bound ligand does not fluoresce (such as due to transient interactions). From Eq. [6], models from Template II give an [L]-independent τ on and a τ off that has a linear dependence on [L]. In particular, one special acyclic scheme from Template II describes a following two-step binding model.
where the first step from the "off" state to the nonfluorescent "dark" state represents a ligand-receptor contact due to ligand diffusion, and the second step to the "on" state models a reaction-limited transition. An alternative mechanism is the 2D lateral diffusion of a ligand on the cell surface to search for a binding molecule after a 3D diffusion in the bulk solution onto the cell surface, which may also contribute complications in analyzing the mean waiting times. From the model in Eq. [8], the apparent association rate constant is given as: where k 1 can be considered as the diffusion-limited rate constant. Note that Eq. [9] contains only two free parameters. Morimatsu et al. 1 hypothesized that frequent and short-term Grb2-EGFR interactions that escaped the instrumental resolution may induce conformation memory in EGFR molecules and thus account for the concentration dependence of the mean off-time. In their experiments, single EGFR molecules were monitored by a total internal reflection fluorescence microscope (TIR-FM) for binding and dissociation with fluorophore (Cy3)labeled adaptor protein Grb2 that reversibly binds to specific phospho-tyrosine residues on EGFR. Statistics of the fluorescence time series showed that k on decreased from 220 µM −1 s −1 to 7.1 µM −1 s −1 when Grb2 concentration increased from 0.1 to 100 nM, whereas the apparent dissociation rate constant k off was found to be constant about 3.4 s −1 . To apply this model to the data reported by Morimatsu et al 1 , we assume that the second-order association rate constant is diffusionlimited, k 1 = 4πDs = 1.51 × 10 3 µM −1 s −1 , where the diffusion coefficient D = 100µm 2 s −1 and the spherical contact radius s = 2 nm as in Ref. 1 . The best fit to the data gives k −1 = 4.54s −1 and k 2 = 0.42s −1 . k −2 is measured by the mean on time: k −2 = 1/τ on = 3.4s −1 . As shown by the fitting quality in Fig. 2, even though k on qualitatively decreases with increasing [L], this scheme is yet to fully quantify the observed dependence on [L], suggesting that an alternative mechanism might better account for the intriguing results reported by Morimatsu et al. 1 .

C. Molecule with multiple binding sites
A single molecule with multiple binding sites to a ligand may potentially cause a ligand dependence of the apparent rate constants that have different forms from a molecule with a single binding site. One cause of ligand dependence is that transient ligand binding escapes observation as discussed below in the missed events section (Section III E). Here, we assume the existence of a site which does not fluoresce for an unknown reason which does result in different ligand dependence of τ off than the above model with a single non-fluorescent binding site.
With Template III, we consider a molecule that has two binding sites as shown in Fig. 1. Although each reversible reaction in Template III involves ligand binding, only one site is monitored for binding. Models constructed from Template III have four state classes: (i) both sites unbound ("off"), (ii) ligand bound to the dark site ("off"), (iii) ligand bound to the bright site ("on"), (iv) ligand bound to both sites ("on"). Calculate the relative occupancy based on the rules described in the previous section (Sec. II B). Then, by Eq. [6], the mean on-time is given by which ranges from 1/k −11 to 1/k −12 as [L] increases from zero to infinity. When the ligand concentration is small the kinetics of the model is biased to the transitions between the upper two states in Template III, whereas the kinetics is shifted to the transitions between the lower two states when the ligand concentration becomes large. If the off rates k −11 and k −12 are nearly identical for fluorescent ligand to dissociate from the single molecule bound or unbound to the non-fluorescent ligand, τ on is independent of [L] and an experiment in this case can only resolve the off rate for the fluorescent binding site. The mean off-time is given by: Using the above equation and the condition of detailed balance, the apparent association rate constant is given as .
As [L] increases from zero to infinity, k on is bounded between k 11 and k 12 , respectively. If ligand binding to the fluorescent site is independent of binding to the nonfluorescent site (see Fig. 1) i.e. k 11 = k 12 , then k on does not have ligand dependence. We use Eq. [12] to fit the data from Morimatsu et al. 1 . In fact, more than one sites on EGFR molecule including phospho-tyrosine sites Y1068 and Y1086 were identified as Grb2 binding sites 18 . As shown in Fig. 2 (dotted curve), models from Template III with the best fit are able to generate a closer agreement to the data. The fitting gives k 11 = 2.16×10 2 µM −1 s −1 , k 12 = 5.84µM −1 s −1 and k 21 /k −21 = 2.02 × 10 3 µM −1 . Since the off-rates k −11 = k −12 = 3.4s −1 , these results indicate a near two orders of magnitude reduction in ligand affinity to the second ligand site after the first ligand binding, suggesting a negative cooperativity of the two binding sites. The parameters also indicate a high affinity non-fluorescent binding with a dissociation constant k −21 /k 21 ≈ 0.5nM .

D. Effect of detailed balance violation
The above analysis breaks down when the assumption that a system obeys detailed balance becomes invalid, which refers to the situation that a system reaches a nonequilibrium steady state because of reactions driven by an implicit external energy source such as a sustained chemical or electrical potential. The results of Eq. [6] are not applicable when detailed balance does not hold. The steady-state occupancy w must be obtained alternatively (see Appendix).
Here, we first use a minimalistic three-state model (Fig. 3) that contains a reaction loop to show that violation of detailed balance causes ligand dependence of rate constants. The model has two "off" states without ligand binding and one "on" state bound to a ligand. To isolate the effect of violation of detailed balance, the model does not invoke any non-fluorescent "dark" state. One can derive the mean waiting times as (see Appendix): We note that τ on is constant regardless the condition of detailed balance. This is a special case for this particular model. In general, violation of detailed balance causes ligand dependence in both mean "on" and "off" times. The association rate constant k on has the same structure as Eq. 12 from Template III and would achieve identical best fit to a time series data set as a model constructed from Template III. Under the condition of detailed balance (k 13 k 32 k 21 = k 31 k 12 k 21 ), the mean off-time τ off is reduced to: which is inversely proportional to [L] and is identical to the result obtained using Eq. [6].
For an arbitrary reaction scheme, we show (see Appendix for detailed derivation using a graph-theoretic method, Fig. 4) that the apparent rate constants are given as rational functions of [L]: where P on and P off are unnormalized steady-state occupancies of "on" and "off" aggregates, respectively. J is an unnormalized inter-aggregate probability flux. These three terms are all polynomials of ligand concentration L. Note that for both k on and k off the denominator and nominator polynomials have corresponding terms of [L] to same powers (See Appendix for mathematical derivation). The exact form of the polynomials is specific to a model topology and the coefficients of the polynomials are in terms of model parameters. If parameters in a model satisfies detailed balance, in the above equations the ligand-dependent terms factor out from both the nominator and denominator and cancel out, leaving k on and k off ligand independent. The in vitro experiments by Morimatsu et al. appear to be done under equilibrium conditions so detailed balance violation is an unlikely explanation for the observed ligand dependence of τ off .

E. Effect of missed events
A missed event is a short-lived binding that escapes the instrumental resolution because it cannot be distinguished from the background noise or because the detector has an intrinsic dead time. Unaccounted missed events distort waiting time distributions and increase mean waiting times. This issue was studied extensively in the field of single ion channel recordings 19,20 . Here, we show that missed events may cause the dependence of k on and k off on [L].
Here we analyze the effect of missed events using the two-state model in Eq. [1]. Assume that the measurement has a fixed dead time σ and that an event is missed if its waiting time is shorter than σ. The apparent mean off-time, τ off , is given as (See Appendix): where the σ accounts for the dead time skipped before the onset of the next detectable on-time interval. The effect of missed events on τ on is ignorable for small ligand concentrations that only induce less frequent binding. This condition is often satisfied in TIR-FM experiments because binding events should be made rare enough to reduce spatial crowding and the background noise and thus to allow detection in changes in the level of fluorescence signals. This technical requirement limits the concentration at the order of 10nM 21 . The association rate constant is obtained as This result is mathematically equivalent to that from the single site protein with non-fluorescent interactions (Eq. [9]). The model can be derived from Template II without the transitions between "on" and "dark" states: where the ligand binding from "off" to "dark" is used to model the missed events.
We apply the above result to estimate the k on and the relative dead time σ/τ on in the experiment by Morimatsu et al. 1 . From the best fitting shown in Fig. 2 (dashed curve), we identify that k on = 1.15 nM −1 s −1 and σ/τ on = 2.19. The dead time σ is more than 2 times the mean on-time and the dissociation constant K d is about 3 nM, suggesting that the affinity between the phosphorylated EGFR and Grb2 is somewhat overestimated, compared to experimental measurements at 700 nM 22 and 30nM 18 . This result reflects the same structural limitation by models from Template II, which did not generate a good fit to the measured data of k on .

F. Effect of an external buffer
Here we examine another possible mechanism in which an ambient buffer may sequester ligands (specifically or non-specifically) and consequently reduce the concentration of free ligands available to bind the single molecule of interest. Unaccounted background buffering may cause ligand dependence of k on . In the absence of a buffer, the effective ligand concentration [L] that interacts with the single molecule equals the total ligand concentration [L] tot . Otherwise, ambient buffers may offset the effective ligand concentration available for binding. We consider a following simple buffering mechanism.
where α ≥ 1 measures the (average) degree of binding cooperativity between the ligand and the buffer group. In the presence of the buffer, the free ligand concentration [L] is the effective concentration of ligands that interact with the single-molecule of interest. We note that the specified ambient buffer B can be a mixture of several kinds of molecules that may interact with the ligand pool. With a phenomenological equilibrium dissociation constant K B and the total ligand concentration [L] tot , we have Considering that buffer B is in excess ( The fitting results suggest that a strong cooperativity (α ≈ 2) existed for Grb2 binding to the buffer. Although the fitting did not directly resolve k + and β, we can make a crude estimation to β ≈ 1.8 × 10 −6 µM (i.e., K B /[B] = 3.6 × 10 −6 µM ) by assuming a diffusion-limited second-order association rate constant k + = 4πDs = 1.51 × 10 3 µM −1 s −1 (with α ≈ 2). We note that although the above external-buffer model produced the closest agreement (partially due to the mathematical properties of the fitting function) to the data by Morimatsu et al. 1 in comparison to the previous ones (Fig. 2), it remains unclear whether the experiment setup introduced a chemical or physical environment that might serve as ambient buffers for Grb2.

IV. DISCUSSION
Molecular binding is an essential biochemical interaction, which can now be probed at the single-molecule level with fluorescence techniques such as Förster (Fluorescence) resonance energy transfer and TIR-FM 21 . These techniques unveil interaction details that are often unavailable in data obtained from ensemble-averaged experiments. Proper interpretation of the fluorescence time series for single molecule binding by its partner protein (or ligand) requires caution. Especially, phenomenological binding constants k on and k off as well as the dissociation constant K d extracted from the fluorescence time series may change as the ligand concentration varies, which carries important information about the binding biochemistry and its experimental environment. Modelbased analysis of the ligand-dependence of kinetic parameters can help to uncover the underlying mechanisms.
In this paper we explore influences by various mechanistic and technical factors, specifically, single-site and multisite non-fluorescent binding, non-equilibrium steady-states, missed events, and ambient buffers, which could potentially introduce dependence of mean waiting times and thus apparent kinetic rate constants on ligand concentration. A combination of these factors can further obscure the analysis of single-molecule kinetics, requiring assistance of appropriate kinetic models.
We have shown that molecular conformation fluctuation (or dynamic disorder) alone does not cause concentration dependence under the condition of detailed balance in models that reach equilibrium steady states as long as each ligand-induced state transition is experimentally resolved (Template I, Fig. 1). In this case, kinetic rate constants inferred from mean waiting times reconcile with those measured by ensemble-averaged experiments.
Unobserved ligand binding, due to unknown biochemical reasons, are the essential sources of ligand dependence of the waiting times, which we analyzed using kinetic models that invoke non-fluorescent liganded states. Different models generate different mathematical structures of ligand dependence. In general, a kinetic rate constant, k on or k off , is a rational function of ligand concentration. Models with non-fluorescent liganded states for a molecule that has a single ligand binding site (Template II, Fig. 1) predict that k on has an inverse linear relationship with ligand concentration [L] (Eq. [9]), whereas k off remains unmodulated by [L]. Models of a molecule with two ligand binding sites with one site non-fluorescent when bound to ligand (Template III, Fig. 1) predict that both k on and k off have sigmoidal shaped relationship with the ligand concentration.
Unmonitored binding can also be caused by short transitions called missed events whose time durations fall within the length of the dead time of the experimental instrument, which has a similar form of concentration dependence by the single site non-fluorescent binding models (Template II, Fig. 1). Our results coincide with a similar three-state model proposed by Crouzy and Sigworth 23 to account for missed events in single-ion channel recordings, in which transient transitions between a closed state to a short-lived state were used to capture events that was off the scope of the instrumental resolution. It was known in analysis of single-ion channel recordings that unaccounted missed events due to fixed dead time can distort the waiting time distributions and cause overestimation of waiting times. Such limitation may be carried over to cause ligand concentration dependence in single-protein fluorescent binding experiments.
The aforementioned models were studied under the condition of detailed balance. Another source that likely causes ligand dependence is the violation of detailed balance in model parameters, which can be studied using non-equilibrium models. This mechanism is ignored in most studies. The typical assumption of a singlemolecule analysis is that the system relaxed to its thermodynamic equilibrium at the steady state. The equilibrium assumption is rather strong and requires the system to meet stringent conditions (the thermodynamics requires the system being isolated without energy and material exchange with its external environment), and it may not be always justified in particular for in vivo systems that entail many energy-driven reactions 24 or for in vitro systems that are sustained by energy sources. Detailed balance violation can be tested by analyzing time series data. For example, two-dimensional joint waiting time distributions that account for two consecutive events, waiting time for binding event followed by that of a dissociation event, can be used to test whether detailed balance holds by checking the time reversibility (also see Ref. 25 for other methods). As a consequence of adopting non-equilibrium model, model parameters might not be constrained by detailed balance. A non-equilibrium model achieves a steady state with net fluxes around reaction loops, which gives rise to ligand dependence of rate constants as rational functions of ligand concentration.
We note that our analysis of the effect of detailed balance is closely related to recent works that study the substrate-dependence of enzymatic turnover rate v in fluctuating enzymes with multiple conformation channels 26,27 . Under the detailed balance condition, the dependence of production formation velocity on substrate concentration [S] maintains the classic Michaelis-Menten form, v = k[S]/( K M + [S]), where effective catalytic rate constant k and apparent Michaelis-Menten constant K m can be derived from kinetic parameters of the model for the enzyme system. When the condition of detailed balance does not hold, v becomes in general a rational function of [S]. To demonstrate that the results from our work can also be applied to analyze the turnover rate of multi-conformational enzymes, consider a general scheme of enzymatic network, where an enzyme fluctuates among several (m) conformations, forming parallel and interconnected catalytic channels. Through each channel, the enzyme engages the substrate and then undergos multiple (n) reversible intermediate steps before finally converting the substrate into a product. The turnover rate can be expressed as the summation of turnover rates in all individual channels: v = m i=1 η ni k i , where η ni is the steadystate residence probability at the last (nth) substratebound step of the ith channel and k i is the catalytic rate constant of that channel. We can write η ni = η ni /Z, where the unnormalized residence probability η ni and the partition function Z = n j=1 m i=1 η ji . As shown in Section II B, under the detailed balance η ni is proportional to the substrate concentration [S] and Z is a linear function of [S]. Thus, the conventional Michaelis-Menten form is preserved in v. Without detailed balance, the turnover rate v assumes a rational functional form of [S], which can be obtained systematically using the graphic method as shown in the Appendix C.
Finally we studied the effect by an external buffer group that sequesters ligands, which if unaccounted could cause strong ligand dependence in rate constants. The extent of the buffering effect depends on biochemical nature of ligand-buffer interaction and the relative availability of the buffer group. It is natural to consider buffering in cellular environment of living cells where molecules are subject to ubiquitous binding reactions in a crowded molecular surrounding by specific and/or non-specific interactions.
We applied our results to analyze the experiment data of labeled Grb2 binding to EGFR molecules by Morimatsu et al. 1 . We examined the possibility that missed events due to transient binding (and showed that this is equivalent to Template II) were the source of the ligand dependence of the apparent association rate constant and found that the best fit could not accurately reproduce the data. The mathematically simplest and best fit resulted from assuming there were background Grb2 buffers characterized by two parameters accounting for cooperativity and affinity. Non-equilibrium models with detailed balance violation were not applied to analyze the data because the in vitro experiments by Morimatsu et al. 1 were apparently performed under the equilibrium condition. Elucidation of the most likely mechanism requires further experimental investigation.
where entry p ij in matrix P is the probability of being in state j at time t when the system was in state i at t = 0. Matrix Q is called "generator matrix". For systems with aggregates "on" and "off", Q can be organized and partitioned as where diagonal blocks contain intra-aggregate transition rates and off-diagonal blocks contain inter-aggregate transition rates. Letters o and c denote "on" and "off" aggregates, respectively. The on-time distribution is given by 12,13 : where vector π o is the steady-state distribution of "on" aggregate entry probabilities over the "on" states, and it is given as the steady-state probability flux into individual "on" states from "off" states normalized by the total probability flux into the "on" aggregate: π o = w c Q co /(w c Q co u o ). The mean on-waiting time is calculated as: The off-time τ off is similarly obtained.

B. Violation of detailed balance in the three-state model
Here we derive the mean "on" and "off" waiting times for the three-state model shown in Fig. 3 in the main text. If we arrange the states in the order as labeled in the figure, the generator matrix for this model is given by: One can find the steady-state occupancy (the left null space of Q) as: where Z is the partition function that normalizes w, i.e., Z = 3 i=1 w i . The magnitude of the net flux (regardless of the direction) around the reaction loop can be calculated as: (A7) If the model satisfies detailed balance, then J n = 0 and the model parameters obey the following constraint: Thus, the equilibrium state probability can be reduced to: The mean "on" and "off" times can be calculated from Eq. [4] in the main text: , τ off = k 23 k 31 + k 13 k 32 k 13 k 23 (k 31 + k 32 ) [L] . (A10) According to k on = 1/τ off [L], k off = 1/τ on , both apparent association and dissociation rate constants k a and k off are ligand independent. This is consistent with the results we shown in the main text (Detailed balance). We note that in the three-state model the mean "off" time is also independent of the intra-aggregate transition rates k 12 and k 21 . When detailed balance does not hold in the model the mean "on" time remains unchanged while the mean "off" time will assume a form as follows: and the apparent association rate constant is: which has a ligand dependence similar in form to that by models from Template III (Eq. [12]) and may potentially achieve fitting to date with the same quality. In this specific model, the mean "on" time is a constant and does not have a dependence on ligand concentration. In general, as we shown below that in non-equilibrium models both τ on and τ off are ligand dependent. C. Ligand dependence in a general scheme for single-site binding Obtaining an analytical solution to the steady-state probability distribution w for a general reaction scheme is unwieldy by directly finding the left null space of the generator matrix Q. As an alternative, one can obtain w by a known graph-theoretical approach used in nonequilibrium statistical mechanics 28 , which solves for the steady-state distribution for a non-equilibrium system, as we show below. Here, note that we only consider single-site fluorescent binding and assume that connections between any two states consist of both forward and backward transitions. Below, we first introduce how to use the method to systematically obtain the steady-state probabilities, and then derive the general formula for the ligand dependence of rate constants.
The method involves enumerating all distinct spanning trees of the topology of a given reaction scheme. A spanning tree of a (undirected) graph is a tree with edges from the original graph that connects all the nodes from the graph. For a topology that has N nodes (states), the maximum number of distinct spanning trees possible is N N −2 for the fully connected topology with every pair of nodes directly connected. Fig. 4A shows all distinct spanning trees for an example four-state model.
For a state k, any given undirected spanning tree has a corresponding directed spanning tree s with all unidirectional edges (transitions) leading toward k (see Fig. 4A for examples). One can view state k as a root of the tree and any directed edge has a direction pointing from an offspring node toward the root. Let V ks be the product of all transition rates associated with the edges in s. It is an established result that the steady-state probability for the system to reside in state k is given by 28 : where N s is the number of distinct spanning trees and Z ≡ N k Ns s V ks is the partition function for normalization purpose. We define the un-normalized steady-state probability vector as: which we will show has each entry as a polynomial function of ligand concentration [L]. We consider an aggregated Markov model of a single molecule binding by a ligand with two aggregates of states, liganded (fluorescent) aggregates and unliganded (non-fluorescent) aggregates.
With the above preparation, we now can show that for the ith state of "on" aggregate in a given directed spanning tree s, V is is a monomial function of [L] with a form V is = α is [L] cs , where integer c s is the number of disjointed subtrees of c states in s. α is is the product of rate constants of the transitions in s. A spanning tree s partitions all c states into c s (1 ≤ c s ≤ N c ) disjointed c [A13] the corresponding directed spanning tree provides an additive contribution to the un-normalized steady-state probability w oi , which consists one and only one term proportional to [L] by a gateway transition from each DCS such that V is = α is [L] cs . The claim of only one [L]-dependent transition from a DCS leading to the ith o state is based on the observation that if more than one such transitions exist there will be loop(s) in the spanning tree, which is an obvious contradiction. The result holds for any arbitrary o state in the spanning tree s. The result V js = β js [L] cs−1 can be derived using similar arguments above for the jth c state.
Summing up contributions from all distinct spanning trees for a topology, we then obtain the un-normalized steady-state probability for state i in the o aggregate and state j in the c aggregate: Therefore, the un-normalized steady-state "on" and "off" probabilities are given as: Ns s=1 α is [L] cs , P off = Nc j=1 Ns s=1 β js [L] cs−1 .
(A16) The steady-state inter-aggregate flux calculated using the unnormalized "on" probability is: where γ i is the ith entry in the vector Q oc u c . Thus, the mean "on" and "off" times are given as: Let q σ = 1 − p σ . The mean dead time is calculated as: t τon e −t/τon dt p σ = τ on − σ q σ p σ . (A24) Assuming binding and dissociation events are independent, the probability that one misses k consecutive short events is p k σ q σ . The apparent mean off-time, τ off , is given as: where the σ accounts for the dead time skipped before the onset of the next detectable on-time interval.
Similarly, the apparent mean on-time is given as τ on = (τ on + p δ τ off )/q δ , where p δ = 1 − e −k+[L]δ is the probability that an off-time is shorter than the dead time δ for detecting a dissociation event. When [L] is very small (i.e., p δ ≈ 0), it is unlikely that a waiting time for a dissociation event falls within the dead time δ and therefore τ on is not significantly affected by missed events.


extfieldsofstudy: Chemistry Medicine

text:
4,4,5,5-Tetramethyl-2-[1,3,6,8-tetrabromo-7-(4,4,5,5-tetramethyl-1,3,2-dioxaborolan-2-yl)pyren-2-yl]-1,3,2-dioxaborolane

The complete molecule of the title compound, C28H28B2Br4O4, is generated by the application of a centre of inversion. In the molecule, the BO2 plane is perpendicular to that through the pyrene ring [dihedral angle = 86.27 (13)°]. In the crystal, molecules stack into columns along the b axis, the closest contact between these being of the type C—Br⋯π.

The complete molecule of the title compound, C 28 H 28 B 2 -Br 4 O 4 , is generated by the application of a centre of inversion. In the molecule, the BO 2 plane is perpendicular to that through the pyrene ring [dihedral angle = 86.27 (13) ]. In the crystal, molecules stack into columns along the b axis, the closest contact between these being of the type C-BrÁ Á Á.

Experimental
Cg1 is the centroid of the C3-C5/C8/C9/C14 benzene ring.  The chemistry of pyrene is strongly position-dependent. For example, in the bromination reaction of 2,7-di-t-butylpyrene, the bromide atoms were connected at the 4,5,9,10-positions of pyrene (Miura & Yamano, 1995). However, to our surprise, when the t-butyl group is changed to pinacol boronate, the bromination reaction resulted in bromination at the 1,3,6,8-positions as confirmed by the crystal structure described herein.
The molecules pack into columns along the b axis, Fig. 2. The most prominent contacts in the structure appear to be of the type C-Br···π, Table 1.

Experimental
The title compound was synthesized via a one-step bromination reaction. The precursor compound, 2,7-di-Bpinpyrene (pin = O 2 C 2 Me 4 ) was prepared using the method of Ir-catalyzed borylation (Coventry et al., 2005). To a stirred mixture of 1.36 g (3.0 mmol) of 2,7-di-Bpinpyrene and a small amount of Fe powder (ca 0.10 g) in 80 ml of CCl 4 was added dropwise a solution of 2.88 g (0.78 mL, 18 mmol) of bromine in 20 ml of CCl 4 at room temperature. After stirring for 5 h, the mixture was slowly poured into ice water. Then the organic layer was collected and washed with aqueous Na 2 S 2 O 3 and dried over MgSO 4 . After evaporation, the residue was crystallized from hexane, giving 1.80 g (78% yield) of gray powdered product. Crystals were grown by slow evaporation from its hexane/dichloromethane solution.

Refinement
Carbon-bound H-atoms were placed in calculated positions [C-H 0.93-0.96 Å, U iso (H) 1.2-1.5U eq (C)] and were included in the refinement in the riding model approximation.   Refinement. Refinement of F 2 against ALL reflections. The weighted R-factor wR and goodness of fit S are based on F 2 , conventional R-factors R are based on F, with F set to zero for negative F 2 . The threshold expression of F 2 > σ(F 2 ) is used only for calculating R-factors(gt) etc. and is not relevant to the choice of reflections for refinement. R-factors based on F 2 are statistically about twice as large as those based on F, and R-factors based on ALL data will be even larger.

Fractional atomic coordinates and isotropic or equivalent isotropic displacement parameters (Å 2 )
x y z U iso */U eq


Synthesis of Cobalt-Zinc Phosphates Templated by Polyamines

Open-framework structures of metal phosphates, especially those zinc and cobalt phosphates, have exhibited many fascinating structural features and potential applications in catalysis, separation processes, and as photoluminantphosphors (Rao et al., 2000; Chang et al., 2004). The feasibility of zinc and cobalt to tetrahedrally coordinate their phosphates makes them to be an important group of transition metal phosphates (Mandal & Natarajan, 2002; Harrison et al., 1991; Neeraj et al., 1999; Fan et al., 2000). Usually, these materials have been synthesized by employing hydrothermal or solvothermal conditions in the presence of organic amines as the structure-directing agents (SDA). Studies have shown that the important role of host–guest charge matching between the inorganic frameworks and the organic amines makes linear polyamines the preferred candidates for the structure-directing agents to form the higher negative charge framework of zinc and cobalt phosphates (Bu et al., 1997, 1998; Feng et al., 1997). Moreover, diethylenetriamine (DETA), triethylenetetramine (TETA) and tetraethylenepentamine (TEPA) become typical amines and have directed the formation of a large number of open-framework structures of zinc (or cobalt) phosphates by varying their concentrations accompanying with the changes in inorganic compositions (Choudhury et al., 2000). Recently, a new synthesis route involving alkylformamide as a template precursors has been reported (Lakiss et al., 2005; Vidal et al., 2000). It seems that organic amines generated in situ by the decomposition of alkylformamide show a special effect on the formation of metal phosphates with a novel structure. In our studies aimed to synthesize transition metal phosphates, more than one cobalt-zinc phosphates with different framework structures have been obtained from a single gel using diethylenetriamine as the organic basis (Lu et al., 2008; Li et al., 2009). It is noticed that the guest species encapsulated in their structures are not diethylenetriamine, but some smaller amines decomposed from it. Controlling the in situ decomposition of these linear polyamines plays the key role for directing the synthesis of pure-phase cobalt-zinc phosphates. According to these matters, large amines can also be used as a source for smaller amines as they will provide hydrolysis decomposition in the reaction mixture during crystallization. In the synthesis a temperature of crystallization is an important factor for controlling the hydrolysis of linear polyamines, and higher temperature makes decomposition of polyamines easier. On the other hand, adding small amines, such as propylamine, butylamine, can control the decomposition of the chain-type polyamines.

Introduction
Open-framework structures of metal phosphates, especially those zinc and cobalt phosphates, have exhibited many fascinating structural features and potential applications in catalysis, separation processes, and as photoluminantphosphors Chang et al., 2004). The feasibility of zinc and cobalt to tetrahedrally coordinate their phosphates makes them to be an important group of transition metal phosphates (Mandal & Natarajan, 2002;Harrison et al., 1991;Neeraj et al., 1999;Fan et al., 2000). Usually, these materials have been synthesized by employing hydrothermal or solvothermal conditions in the presence of organic amines as the structure-directing agents (SDA). Studies have shown that the important role of host-guest charge matching between the inorganic frameworks and the organic amines makes linear polyamines the preferred candidates for the structure-directing agents to form the higher negative charge framework of zinc and cobalt phosphates (Bu et al., 1997(Bu et al., , 1998Feng et al., 1997). Moreover, diethylenetriamine (DETA), triethylenetetramine (TETA) and tetraethylenepentamine (TEPA) become typical amines and have directed the formation of a large number of open-framework structures of zinc (or cobalt) phosphates by varying their concentrations accompanying with the changes in inorganic compositions (Choudhury et al., 2000). Recently, a new synthesis route involving alkylformamide as a template precursors has been reported (Lakiss et al., 2005;Vidal et al., 2000). It seems that organic amines generated in situ by the decomposition of alkylformamide show a special effect on the formation of metal phosphates with a novel structure. In our studies aimed to synthesize transition metal phosphates, more than one cobalt-zinc phosphates with different framework structures have been obtained from a single gel using diethylenetriamine as the organic basis (Lu et al., 2008;Li et al., 2009). It is noticed that the guest species encapsulated in their structures are not diethylenetriamine, but some smaller amines decomposed from it. Controlling the in situ decomposition of these linear polyamines plays the key role for directing the synthesis of pure-phase cobalt-zinc phosphates. According to these matters, large amines can also be used as a source for smaller amines as they will provide hydrolysis decomposition in the reaction mixture during crystallization. In the synthesis a temperature of crystallization is an important factor for controlling the hydrolysis of linear polyamines, and higher temperature makes decomposition of polyamines easier. On the other hand, adding small amines, such as propylamine, butylamine, can control the decomposition of the chain-type polyamines.

Synthesis and initial characterization
The main reactants are Zn(CH 3 COO) 2 (denoted as ZnAc 2 ), Co(CH 3 COO) 2 (denoted as CoAc 2 ), H 3 PO 4 , linear polyamines, and ethylene glycol. In a typical synthetic procedure, solution A was prepared by mixing ZnAc 2 ·2H 2 O with 85% H 3 PO 4 and ethylene glycol (EG), the mixture was stirred at room temperature until ZnAc 2 ·2H 2 O was dissolved. Solution B was prepared by mixing CoAc 2 ·4H 2 O, distilled water, and 85% H 3 PO 4 . Solution A and a kind of linear polyamine, such as diethylenetriamine were added to solution B. Then the pH value of the mixture was adjusted to 8.0 by adding an assistant amine. Finally the mixture was transferred to a stainless steel autoclave and heated at 140˚C~160˚C for 3-10 days. Specific gel compositions and crystallization conditions of three cobalt-zinc phosphates are shown in Table. 1). XRD patterns of the three compounds also show their phases being distinct from each other (Fig.1). www.intechopen.com

Crystal structure of CoZnPO 4 -V
Single crystal analysis and Inductively Coupled Plasma atomic emission spectroscopy (ICP) analysis reveal that CoZnPO 4 -V is composed of C 12 N 9 H 52 Zn 4.5 Co 1.5 P 8 O 32 ; it crystallizes in the monoclinic system, space group P2 1 /c (No. 14), with a = 31.936(3) Å, b = 8.3775(7) Å, c = 15.7874(13) Å, = = 90°, = 97.0530(10)°, V = 4191.8(6) Å 3 , and z = 4. Three edge-shared 4-ring units of Zn-O and P-O tetrahedral groups could be the building unit, which was named as 4 3 unit (Fig. 3a). The neighboring units link with each other through 4-ring to form an infinite zigzag chain along the b-axis direction (Fig. 3b). 2D (4, 10) nets are conformably formed of the chains of 4-ring units from the corner-shared 4 3 units ( Fig. 3b). Viewed along the plane parallel to b-axis, the neighboring 4 3 units link with each other in a corner-shared way to form another infinite chain (Fig. 3c). The neighboring chains link together from 4-ring every two 4 3 units to form a layer. Therefore, 16-ring channels are constructed (Fig. 3c). On the other hand, the wall of each 16-ring channel can be described as the coiling of a net composed of 10-rings and infinite 4-ring zigzag chain (Fig. 3b).

Crystal structure of CoZnPO 4 -VI
Single crystal analysis and Inductively Coupled Plasma atomic emission spectroscopy (ICP) analysis reveal that CoZnPO 4 -VI is composed of C 6 N 4 H 22 Co 0.16 Zn 5.84 (PO 4 ) 4 (HPO 4 ) 2 ; it crystallizes in the triclinic system, space group P 1 (No. 2), with a = 5.1919(10) Å, b = 8.7263 (17)  The three-dimensional framework of CoZnPO 4 -VI can be described by with two secondary structural units (SBU-a and SBU-b). SBU-a (Fig.4a) is a special structure composed of two 6rings ladder linked face to face; SBU-b (Fig.4b) is a common 4-ring ladder in the molecular sieves. The linkage of SBU-a conducted by sharing 4-rings form the layers along ab-plane. SBU-a and SBU-b connect with one sharing oxygen atom in the structure, forming the threedimensional framework of CoZnPO 4 -VI and resulting in 16-ring large channels along the direction of c-axis (Fig.4c). The arms of channels are constituted mainly by SBU-a and SBU-b.

Controlling the in situ decomposition of linear polyamines to synthesize open-framework cobalt-zinc phosphates
CoZnPO 4 -I, CoZnPO 4 -II, and CoZnPO 4 -III, have been synthesized from a single gel at 200 °C. Organic amine in the starting mixture is diethylenetriamine. Structural analysis reveals that species encapsulated in the structures of CoZnPO 4 -I, CoZnPO 4 -II, and CoZnPO 4 -III are ammonium, ethylenediaminum, and diethylenetriaminum, respectively, which can be attributed to the partly in situ decomposition of diethylenetriamine molecules during the crystallization procedure ( Fig.13a, 13b, and 13c). When a second small organic amine, such as propylamine, or butylamine, was added to reacting mixture, single phase of CoZnPO 4 -I, CoZnPO 4 -II, or CoZnPO 4 -III can be obtained respectively, according to the kind or concentration of small organic amine, which plays an important role on controlling the in situ decomposition of linear polyamines.  Three open-framework cobalt-zinc phosphates have been synthesized hydrothermally from a single gel at 200°C by using DETA as the SDA with n-propylamine as the assistant organic amine. They are all large vivid blue crystals with sizes more than 400 μm. Scanning electron micrograph (SEM) images show their morphologies being distinct from each other, shaped as elongated hexagonal-dipyramid, square bifrustum, and cubical, respectively (Fig.5). XRD patterns of CoZnPO 4 -I, -II, and -III (Fig.6) have shown that they are three distinguishing phases with different diffraction peaks.

Synthesis and initial characterization
www.intechopen.com

Crystal structures of three compounds
It is clear that there are only diethylenetriamine and n-propylamine as the SDAs present in the starting materials. However, organic species encapsulated in the structure of CoZnPO 4 -I, -II, and -III are ammonium, ethylenediaminum, and diethylenetriaminum cations, respectively (Fig.13a, 13b, 13c). The appearance of ethylenediaminum and ammonium may be attributed to the in situ decomposition of diethylenetriamine during the crystallization. Single crystal analysis reveals that CoZnPO 4 -I is a novel cobalt-zinc phosphate. It is composed of (NH 4 ) 2 Co 0.34 Zn 1.66 (PO 4 ) 2, which was modeled as (NH 4 ) 2 Co 2 (PO 4 ) 2 in the single crystal data and crystallizes in the hexagonal system, space group P6 3 (No. 173), with a =b =10.7207(5)Å, c =8.7241(8) Å, = =90°, λ =120°, V=868.36(10) Å 3 , and z =4. Viewing along the plane parallel to c-axis, a net consisting of 6-ring channels has been found. The 6ring channels with regularly alternating Zn(Co)-and P-centered tetrahedral can be considered as the building units of CoZnPO 4 -I (Fig. 7a). J. V. Smith determined a simple hexagonal net, which has only one type of 3-connected node and is described as the Schläfli symbol 6 3 because each node lies between three circuits of 6 nodes (Smith, J. V. Am. Mineral. 1977).For any node in a simple hexagonal net in horizational position, an additional perpendicular linkage pointing either upward (U) or downward (D) may result in the linkage joined to a node of another simple hexagonal net lying either above or below the first hexagonal net. Eight ways of the sequence of U and D linkages around each 6-ring has been used to enumerate 4-connected, 3-dimensional nets and classification of framework of silicates and aluminophosphates. The linking sequences of UDUDUD around each 6-ring produces a dense tridymite framework, whereas the UUUDDD results in an open framework ABW containing 8-ring channels. However, infinity of frameworks can be produced if more than one sequence occurs in the hexagons of a 6 3 net. This has just happened in the structure of CoZnPO 4 -I. In its framework, the 6-rings of tetrahedral are present in two different conformations based on the more than one type circuit as the following sequence of up and down tetrahedral: UDUDUD and UUUDDD in the proportion 1:3. A second building unit, which is called a 6 9 4 6 cage, has been formed (Fig. 13a).

www.intechopen.com
On the other hand, it can be found along the (010) plane that double 4-rings and one 6-ring linked in turn by sharing edges to form an infinite chain (Fig. 7b). These chains are crosslinked along the c-axis direction to the similar chains from the linkage Zn-O-P or P-O-Zn to form the sheet structure with edgesharing infinite 10-ring chains (Fig. 7b). Adjacent sheets are cross-linked with each other along the vertical direction of (010) plane via the Zn-O-P bonds displaced by +1/2 in the [001] direction to form the three-dimensional framework of CoZnPO 4 -I (Fig.7c). Each 10-ring is blocked by double 4-and 6-rings around it. So the tendency to construct 10-ring channels is interdicted with only clathrate 6 9 4 6 cages to be formed.   6 ]. Its structure is made of CoO 6 , ZnO 4 (or CoO 4 ), and PO 4 polyhedra. The cobalt octahedron is surrounded by six zinc tetrahedra to form the CoZn(or Co) 6 O 6 cluster which are linked via six PO 4 tetrahedra to other clusters to form its framework. The clusters are arranged i n s u c h a m a n n e r t h a t e a c h c l u s t e r i s displaced by half its a-distance from its neighbors, forming a honeycomb layer. The next layer of clusters is identical to the first but is displaced along the a-axis by half the unit cell so that the honeycomb channels are capped. The framework of AB-type stacking results in the formation of 8-ring channels at an angle to a-axis (Fig. 9).

The in situ framework transformation of cobalt-zinc phosphates in the presence of linear polyamines
Although linear polyamines cannot decompose at low crystallization temperature, experiments have shown the crystallization time also affects the crystal structure transformation. Another typical result that we have found is the transformation of two cobalt-zinc phosphates, CoZnPO 4 -VI and CoZnPO 4 -VII. CoZnPO 4 -VI was synthesized by using TETA as the template, which can be transformed into CoZnPO 4 -VII with the increase of crystallization time in 140°C.  .10) shows the morphology of CoZnPO 4 -VII is as a relatively large irregular block, it is distinct from CoZnPO 4 -VI (Fig.1c). XRD patterns of CoZnPO 4 -VI and CoZnPO 4 -VII have also shown they are two distinguishing phases with different diffraction peaks ( Fig.10 and Fig.1c).

Crystal structure of CoZnPO 4 -VII
Single crystal analysis and Inductively Coupled Plasma atomic emission spectroscopy (ICP) analysis reveal that CoZnPO 4 -VII is composed of C 6 H 22 Co 0 N 4 O 16 P 4 Zn 4 ; it crystallizes in the triclinic system, space group P 1 (No. 2), with a = 8.0871 (16)   The biggest feature in the structure of CoZnPO 4 -VII is that there are two cross-linked 8-ring channels in the structure. The three-dimensional framework of CoZnPO 4 -VII can be also described as the stack of a secondary structural unit (SBU-c) (Fig.11a). In the direction of aaxis, the connection of 8-rings by three 4-rings forms a ladder (SBU-c) (Fig.11b); in the direction of b-axis ladders are joined by sharing edges, spreading the layers running abplane. Then adjacent layers are connected with 4-rings. All described above spread the three-dimensional framework of CoZnPO 4 -VII and exhibit a 8-ring channel system running along the direction of a-axis (Fig.11c). The arms of channels are constituted mainly by SBU-c. At the same time they form another 8-ring channel system running the direction of c-axis (Fig.11d).

Relationship between CoZnPO 4 -VI and CoZnPO 4 -VII
The similarities between CoZnPO 4 -VI and CoZnPO 4 -VII provide the possibility of transformation. On one hand, both of them are composed of Co,Zn,P,O elements, and all Zn (Co) atoms and P atoms in their structures are tetrahedrally coordinated by oxygens, which have been seen from their asymmetric unit structures (Fig.12). On the other hand, the templates filled in their inorganic frameworks are TETA, and both of them crystallize in the triclinic system, space group P 1 (No. 2). All of these result in the transformation of CoZnPO 4 -VI and CoZnPO 4 -VII with the increase of crystallization time.

Investigation of the in situ decomposition of linear polyamines during hydrothermal synthesis of cobalt-zinc phosphates
In this section, polyamines such as diethylenetriamine, triethylenetetramine, and tetraethylenepentamine were used as structure directing agents to investigate the in situ decomposition of linear polyamines during the crystallization of microporous zinc-cobalt phosphate materials. Mixtures of water and ethylene glycol were used as solvents.

Experimental methods
The important basis for determination of polyamines' decomposition is consideration of the template molecules in the framework of all crystals. From the summarized above, we know that the template molecules in the framework of CoZnPO 4 -I, CoZnPO 4 -II, CoZnPO 4 -III, CoZnPO 4 -IV, CoZnPO 4 -V and CoZnPO 4 -VI are NH 4 + , en, DETA, DETA, DETA and TETA, respectively (Fig.13).
Previous studies have shown preliminary that temperature is an important factor for controlling the hydrolysis of linear polyamines. Other, the addition of different assistant templates has different effects on decomposition of propylamines while using the straightchain polyamines as a precursor of the template in the synthesis of cobalt-zinc phosphates. In order to investigate the decomposition of polyamines, the same ratio of gel should react for a certain time at different temperatures, and the pH value of gel should be adjusted to pH (7.0-9.0) by the addition of different assistant templates. The selected molar ratio of gel is 2ZnAc 2 •2H 2 O:1CoAc 2 •4H 2 O:9.6EG:10.5H 3 PO 4 :95.4H 2 O:1.8R (R is DETA, TETA, TEPA respectively).

Discussions of factors affecting decomposition of polyamine
The results are shown in Table. 2 when crystallizing at different temperatures for 6 days.  However, when an assistant amine, such as n-propylamine, was added to the mixture to increase the pH value to 7-9, even if the mixture is crystallized at 200 °C, three metal phosphates (CoZnPO 4 -I, -II, -III) have been obtained with diethylenetriamine still encapsulated in the pores of CoZnPO 4 -III. It indicates that the presence of n-propylamine can restrain the hydrolysis of diethylenetriamine during the crystallization. It can be found from Table 2 that, apart from n-propylamine, other assistant amines, such as nbutylamine, dipropylamine, and tripropylamine, can also restrain the decomposition of diethylenetriamine to some extent. The influence of second organic amine may be attributed to the effect of the pH on hydrolysis of the linear polyamines. In the case of triethylenetetramine or tetraethylenepentamine, the controlling effects of assistant amine are less prominent than for diethylenetriamine. CoZnPO 4 -II and -III have only been obtained as the pure phases by varying assistant amine or crystallizing temperature. The reacting mixture with the molar ratio of 2:1:9.6:10.5:95.4:1.8 ZnAc2:CoAc2:EG:H3PO4:H2O:R was crystallized for 6 days; y = yes, compound has been obtained; n = no,compound has not been obtained. Table 2. Products obtained by adding various assistant amines to adjust the pH value of the mixture to 7-9 a

Conclusions
Synthesis of zinc-cobalt phosphate materials have been investigated by using linear polyamines, such as tetraethylene-pentamine, triethylenetetramine or diethylenetriamine as SDA. During the investigation of decomposition of polyamines seven cobalt-substituted zinc phosphates were obtained: CoZnPO 4 -I-VII. CoZnPO 4 -I, CoZnPO 4 -II, and CoZnPO 4 -III, have been synthesized from a single gel at 200 °C in the presence of diethylenetriamine as the structure-directing agent. Structural analysis reveals that species encapsulated in the structures of CoZnPO 4 -I, CoZnPO 4 -II, and CoZnPO 4 -III are ammonium, ethylenediaminum, and diethylenetriaminum, respectively, which can be attributed to the partly in situ decomposition of diethylenetriamine molecules during the crystallization procedure. CoZnPO 4 -IV，CoZnPO 4 -V and CoZnPO 4 -VI with differently unique 3D framework topologies have been synthesized under hydrothermal or solvothermal conditions in the presence of linear polyamines as structure directing agents (SDA). CoZnPO 4 -VI and CoZnPO 4 -VII, which were synthesized by using TETA as the template, are symbiotic and CoZnPO 4 -VI can be transformed into CoZnPO 4 -VII with the increase of crystallization time in 140°C.
In the synthesis of zinc-cobalt phosphates, temperature is an important factor, favoring the hydrolysis of linear polyamines at higher temperature. Investigations have also found that the in situ decomposition of tetraethylene-pentamine, triethylenetetramine, or diethylenetriamine can be controlled by using second organic amines, such as propylamine as the assistant organic additive. The presence of an assistant amine can restrain the hydrolysis of diethylenetriamine during the crystallization. If synthesis conditions can be controlled, different structures can be transformed into each other.

Acknowledgment
This work was funded by the Natural Science Foundation of China (20873069) and State Basic Research Project (2009CB623502).