Critical Behaviour in DOPC/DPPC/Cholesterol Mixtures: Static 2H NMR Line Shapes Near the Critical Point

Static ²H NMR spectroscopy is used to study the critical behavior of mixtures of 1,2-dioleoyl-phosphatidylcholine/1,2-dipalmitoyl-phosphatidylcholine (DPPC)/cholesterol in molar proportion 37.5:37.5:25 using either chain perdeuterated DPPC-d₆₂ or chain methyl deuterated DPPC-d₆. The temperature dependence of the first moment of the ²H spectrum of the sample made with DPPC-d₆₂ and of the quadrupolar splittings of the chain-methyl-labeled DPPC-d₆ sample are directly related to the temperature dependence of the critical order parameter η, which scales as [(T_c?T)/T_c]^β_c near the critical temperature. Analysis of the data reveals that for the chain perdeuterated sample, the value of T_c is 301.51 ± 0.1 K, and that of the critical exponent, β_c = 0.391 ± 0.02. The line shape analysis of the methyl labeled (d₆) sample gives T_c = 303.74 ± 0.07 K and β_c = 0.338 ± 0.009. These values obtained for β_c are in good agreement with the predictions of a three-dimensional Ising model. The difference in critical temperature between the two samples having nominally the same molar composition arises because of the lowering of the phase transition temperature that occurs due to the perdeuteration of the DPPC.     

Rectal,telemetry pill and tympanic membrane thermometry during exercise heat stress

We compared the accuracy of an ingestible telemetry pill method of core temperature (T_c) measurement and an infrared tympanic membrane thermometer to values from a rectal thermistor during exercise-induced heat stress. Ten well-trained subjects completed four exercise trials consisting of 40 min constant-load exercise at 63% of maximum work rate followed by a 16.1 km time trial at 30 °C and 70% relative humidity. Temperature at rest was not different between the three methods of T_c measurement (T_re: 37.2±0.3 °C; T_p: 37.2±0.2 °C; T_ty: 37.1±0.3 °C; P=0.40$P=0.40$). Temperature rose continuously during the exercise period (ΔT_re: 2.2±0.5 °C; ΔT_p: 2.2±0.5 °C; ΔT_ty: 1.9±0.5 ±°C and there were no differences between T_re and T_p measurements at any time throughout exercise (P=0.32$P=0.32$). While there were no differences between T_re and T_ty after 10 min (P=0.11$P=0.11$) and 20 min (P=0.06$P=0.06$) of exercise, T_ty was lower than T_re after 30 min of exercise (P<0.01$P<0.01$) and remained significantly lower throughout the remainder of the exercise period. These results demonstrate that the telemetry pill system provides a valid measurement of trunk temperature during rest and exercise-induced thermal strain. T_ty was significantly lower than T_re when temperature exceeded 37.5 °C. However, whether these differences are due to selective brain cooling or imperfections in the tympanic membrane thermometer methodology remains to be determined.     

Ultrashort echo imaging of cyclically loaded rabbit patellar tendon

Tendinopathy affects individuals who perform repetitive joint motion. Magnetic resonance imaging (MRI) is frequently used to qualitatively assess tendon health, but quantitative evaluation of inherent MRI properties of loaded tendon has been limited. This study evaluated the effect of cyclic loading on _T2^?${T_2}^{?}$ values of fresh and frozen rabbit patellar tendons using ultra short echo (UTE) MRI. Eight fresh and 8 frozen rabbit lower extremities had MR scans acquired for tendon _T2^?${T_2}^{?}$ evaluation. The tendons were then manually cyclically loaded for 100 cycles to 45 N at approximately 1 Hz. The MR scanning was repeated to reassess the _T2^?${T_2}^{?}$ values. Analyses were performed to detect differences of tendon _T2^?${T_2}^{?}$ values between fresh and frozen samples prior to and after loading, and to detect changes of tendon _T2^?${T_2}^{?}$ values between the unloaded and loaded configurations. No difference of _T2^?${T_2}^{?}$ was found between the fresh and frozen samples prior to or after loading, p=0.8 and p  =0.1, respectively. The tendons had significantly shorter _T2^?${T_2}^{?}$ values, p  =0.023, and reduced _T2^?${T_2}^{?}$ variability, p  =0.04, after cyclic loading. Histologic evaluation confirmed no induced tendon damage from loading. Shorter _T2^?${T_2}^{?}$, from stronger spin–spin interactions, may be attributed to greater tissue organization from uncrimping of collagen fibrils and lateral contraction of the tendon during loading. Cyclic tensile loading of tissue reduces patellar tendon _T2^?${T_2}^{?}$ values and may provide a quantitative metric to assess tissue organization.     

The gating charge should not be estimated by fitting a two-state model to a Q-V curve

The voltage dependence of charges in voltage-sensitive proteins, typically displayed as charge versus voltage (Q-V) curves, is often quantified by fitting it to a simple two-state Boltzmann function. This procedure overlooks the fact that the fitted parameters, including the total charge, may be incorrect if the charge is moving in multiple steps. We present here the derivation of a general formulation for Q-V curves from multistate sequential models, including the case of infinite number of states. We demonstrate that the commonly used method to estimate the charge per molecule using a simple Boltzmann fit is not only inadequate, but in most cases, it underestimates the moving charge times the fraction of the field.Many ion channels, transporters, enzymes, receptors, and pumps are voltage dependent. This voltage dependence is the result of voltage-induced translocation of intrinsic charges that, in some way, affects the conformation of the molecule. The movement of such charges is manifested as a current that can be recorded under voltage clamp. The best-known examples of these currents are “gating” currents in voltage-gated channels and “sensing” currents in voltage-sensitive phosphatases. The time integral of the gating or sensing current as a function of voltage (V) is the displaced charge Q(V), normally called the Q-V curve.It is important to estimate how much is the total amount of net charge per molecule (Q_max) that relocates within the electric field because it determines whether a small or a large change in voltage is necessary to affect the function of the protein. Most importantly, knowing Q_max is critical if one wishes to correlate charge movement with structural changes in the protein. The charge is the time integral of the current, and it corresponds to the product of the actual moving charge times the fraction of the field it traverses. Therefore, correlating charge movement with structure requires knowledge of where the charged groups are located and the electric field profile. In recent papers by Chowdhury and Chanda (2012) and Sigg (2013), it was demonstrated that the total energy of activating the voltage sensor is equal to Q_max V_M, where V_M is the median voltage of charge transfer, a value that is only equal to the half-point of activation V_1/2 for symmetrical Q-V curves. V_M is easily estimated from the Q-V curve, but Q_max must be obtained with other methods because, as we will show here, it is not directly derived from the Q-V curve in the general case.The typical methods used to estimate charge per molecule Q_max include measurements of limiting slope (Almers, 1978) and the ratio of total charge divided by the number of molecules (Schoppa et al., 1992). The discussion on implementation, accuracy, and reliability of these methodologies has been addressed many times in the literature, and it will not be discussed here (see Sigg and Bezanilla, 1997). However, it is worth mentioning that these approaches tend to be technically demanding, thus driving researchers to seek alternative avenues toward estimating the total charge per molecule. Particularly, we will discuss here the use of a two-state Boltzmann distribution for this purpose. Our intention is to demonstrate that this commonly used method to estimate the charge per molecule is generally incorrect and likely to give a lower bound of the moving charge times the fraction of the field.The two-state Boltzmann distribution describes a charged particle that can only be in one of two positions or states that we could call S₁ and S₂. When the particle with charge Q_max (in units of electronic charge) moves from S₁ to S₂, or vice versa, it does it in a single step. The average charge found in position S₂, Q(V), will depend on the energy difference between S₁ and S₂, and the charge of the particle. The equation that describes Q(V) is:Q(V)=Qmax1+exp[−Qmax(V−V1/2)kT],(1)where V_1/2 is the potential at which the charge is equally distributed between S₁ and S₂, and k and T are the Boltzmann constant and absolute temperature, respectively. The Q(V) is typically normalized by dividing Eq. 1 by the total charge Q_max. The resulting function is frequently called a “single Boltzmann” in the literature and is used to fit normalized, experimentally obtained Q-V curves. The fit yields an apparent V_1/2 (V∼1/2) and an apparent Q_MAX (Q−max), and this last value is then attributed to be the total charge moving Q_max. Indeed, this is correct but only for the case of a charge moving between two positions in a single step. However, the value of Q−max thus obtained does not represent the charge per molecule for the more general (and frequent) case when the charge moves in more than one step.To demonstrate the above statement and also estimate the possible error in using the fitted Q−max from Eq. 1, let us consider the case when the gating charge moves in a series of n steps between n + 1 states, each step with a fractional charge z_i (in units of electronic charge e₀) that will add up to the total charge Q_max.S1↔μ1S2↔μ2…Si↔μiSi+1…Sn↔μnSn+1The probability of being in each of the states S_i is labeled as P_i, and the equilibrium constant of each step is given byμi=exp[zi(V−Vi)kT], i=1…n,where z_i is the charge (in units of e₀) of step i, and V_i is the membrane potential that makes the equilibrium constant equal 1. In steady state, the solution of P_i can be obtained by combiningPi+1Pi=μi, i=1…n and ∑i=1i=n+1Pi=1,givingPi+1=∏m=1iμm1+∑j=1n∏k=1jμk, i=1…n and P1=11+∑j=1n∏k=1jμk.We define the reaction coordinate along the moved charged q asqi=∑j=1izj, i=1…n.The Q-V curve is defined asQ(V)=∑i=1nqiPi+1.Then, replacing P_i yieldsQ(V)=∑i=1n[∑j=1izj][∏m=1iμm]1+∑j=1n∏k=1jμk,or written explicitly as a function of V:Q(V)=∑i=1n[∑j=1izj][∏m=1iexp[zm(V−Vm)kT]]1+∑j=1n∏k=1jexp[zk(V−Vk)kT].(2)Eq. 2 is a general solution of a sequential model with n + 1 states with arbitrary valences and V_i’s for each transition. We can easily see that Eq. 2 has a very different form than Eq. 1, except when there is only a single transition (n = 1). In this latter case, Eq. 2 reduces to Eq. 1 because z₁ and V₁ are equal to Q_max and V_1/2, respectively. For the more general situation where n > 1, if one fits the Q(V) relation obeying Eq. 2 with Eq. 1, the fitted Q−max value will not correspond to the sum of the z_i values (see examples below and Fig. 1). A simple way to visualize the discrepancy between the predicted value of Eqs. 1 and 2 is to compute the maximum slope of the Q-V curve. This can be done analytically assuming that V_i = V_o for all transitions and that the total charge Q_max is evenly divided among those transitions. The limit of the first derivative of the Q(V) with respect to V evaluated at V = V_o is given by this equation:dQ(V)dV|V=V0=Qmax(n+2)12nkT.(3)From Eq. 3, it can be seen that the slope of the Q-V curve decreases with the number of transitions being maximum and equal to Q_max /(4kT) when n = 1 (two states) and a minimum equal to Q_max /(12kT) when n goes to infinity, which is the continuous case (see next paragraph).Open in a separate windowFigure 1.Examples of normalized Q-V curves for a Q_max = 4 computed with Eq. 2 for the cases of one, two, three, four, and six transitions and the continuous case using Eq. 5 (squares). All the Q-V curves were fitted with Eq. 1 (lines). The insets show the fitted valence (Q−max) and half-point (V∼1/2).

Infinite number of steps

Eq. 2 can be generalized to the case where the charge moves continuously, corresponding to an infinite number of steps. If we makez_i = Q_max/n, ?i = 1…n, ??V_i = V_o, ?i = 1…n, then all µ_i = µ, and we can write Eq. 2 as the normalized Q(V) in the limit when n goes to infinity:Qnor(V)=limn→∞∑i=1n[∑j=1iQmaxn]∏m=1iexp[Qmax(V−Vo)nkT]Qmax[1+∑i=1n∏j=1iexp[Qmax(V−Vo)nkT]]=[Qmax(V−Vo)−kT]exp[Qmax(V−Vo)kT]+kTQmax(V−Vo)[exp[Qmax(V−Vo)kT]−1].(4)Eq. 4 can also be written asQnor(V)=12[1+coth[Qmax(V−Vo)2kT]−2kTQmax(V−V0)],(5)which is of the same form of the classical equation of paramagnetism (see Kittel, 2005).

Examples

We will illustrate now that data generated by Eq. 2 can be fitted quite well by Eq. 1, thus leading to an incorrect estimate of the total charge moved. Typically, the experimental value of the charge plotted is normalized to its maximum because there is no knowledge of the absolute amount of charge per molecule and the number of molecules. The normalized Q-V curve, Q_nor, is obtained by dividing Q(V) by the sum of all the partial charges.Fig. 1 shows Q_nor computed using Eq. 2 for one, two, three, four, and six transitions and for the continuous case using Eq. 5 (squares) with superimposed fits to a two-state Boltzmann distribution (Eq. 1, lines). The computations were done with equal charge in each step (for a total charge Q_max = 4e₀) and also the same V_i = −25 mV value for all the steps. It is clear that fits are quite acceptable for cases up to four transitions, but the fit significantly deviates in the continuous case.Considering that experimental data normally have significant scatter, it is then quite likely that the experimenter will accept the single-transition fit even for cases where there are six or more transitions (see Fig. 1). In general, the case up to four transitions will look as a very good fit, and the fitted Q−max value may be inaccurately taken and the total charge transported might be underestimated. To illustrate how bad the estimate can be for these cases, we have included as insets the fitted value of Q−max for the cases presented in Fig. 1. It is clear that the estimated value can be as low as a fourth of the real total charge. The estimated value of V∼1/2 is very close to the correct value for all cases, but we have only considered cases in which all V_i’s are the same.It should be noted that if µ_i of the rightmost transition is heavily biased to the last state (V_i is very negative), then the Q−max estimated by fitting a two-state model is much closer to the total gating charge. In a three-state model, it can be shown that the fitted value is exact when V₁→∞ and V₂→−∞ because in that case, it converts into a two-state model. Although these values of V are unrealistic, the fitted value of Q−max can be very close to the total charge when V₂ is much more negative than V₁ (that is, V₁ >> V₂). On the other hand, If V₁ << V₂, the Q-V curve will exhibit a plateau region and, as the difference between V₁ and V₂ decreases, the plateau becomes less obvious and the curve looks monotonic. These cases have been discussed in detail for the two-transition model in Lacroix et al. (2012).We conclude that it is not possible to estimate unequivocally the gating charge per sensor from a “single-Boltzmann” fit to a Q-V curve of a charge moving in multiple transitions. The estimated Q−max value will be a low estimate of the gating charge Q_max, except in the case of the two-state model or the case of a heavily biased late step, which are rare occurrences. It is then safer to call “apparent gating charge” the fitted Q−max value of the single-Boltzmann fit.

Addendum

The most general case in which transitions between states include loops, branches, and steps can be derived directly from the partition function and follows the general thermodynamic treatment by Sigg and Bezanilla (1997), Chowdhury and Chanda (2012), and Sigg (2013). The reaction coordinate is the charge moving in the general case where it evolves from q = 0 to q = Q_max by means of steps, loops, or branches. In that case, the partition function is given byZ=∑iexp(qi(V−Vi)kT).(6)We can compute the mean gating charge, also called the Q-V curve, asQ(V)=〈q〉=kTZ’Z=kTdlnZdV=∑iqiexp(qi(V−Vi)kT)∑iexp(qi(V−Vi)kT).(7)The slope of the Q-V is obtained by taking the derivative of 〈q〉 with respect to V:dQ(V)dV=(kT)2d2lnZdV2.(8)Let us now consider the gating charge fluctuation. The charge fluctuation will depend on the number of possible conformations of the charge and is expected to be a maximum when there are only two possible charged states to dwell. As the number of intermediate states increases, the charge fluctuation decreases. Now, a measure of the charge fluctuation is given by the variance of the gating charge, which can be computed from the partition function as:〈Δq2〉=〈q2〉−〈q〉2=(kT)2(Z’’Z−(Z’Z)2)=(kT)2d2lnZdV2.(9)But the variance (Eq. 9) is identical to the slope of Q(V) (Eq. 8). This implies that the slope of the Q-V is maximum when there are only two states.     

Ionic microenvironmental effects on triplex DNA stabilization: Cationic counterion effects on poly(dT)·poly(dA)·poly(dT)

The structure and conformation of nucleic acids are influenced by metal ions, polyamines, and the microenvironment. In poly(purine) · poly(pyrimidine) sequences, triplex DNA formation is facilitated by metal ions, polyamines and other ligands. We studied the effects of mono- and di-valent metal ions, and ammonium salts on the stability of triple- and double-stranded structures formed from poly(dA) and poly(dT) by measuring their respective melting temperatures. In the presence of metal ions, the absorbance versus temperature profile showed two transitions: T_m1 for triplex to duplex and single stranded DNA, and T_m2 for duplex DNA melting to single stranded DNA. Monovalent cations (Li⁺, Na⁺, K⁺, Rb⁺, Cs⁺ and ₄NH⁺${{\text{NH}}_4}^{+}$) promoted triplex DNA at concentrations ≥150 mM. T_m1 varied from 49.8 °C in the presence of 150 mM Li⁺ to 30.6 °C in the presence of 150 mM K⁺. ₄NH⁺${{\text{NH}}_4}^{+}$ was very effective in stabilizing triplex DNA and its efficacy decreased with increasing substitution of the hydrogen atoms with methyl, ethyl, propyl and butyl groups. As in the case of monovalent cations, a concentration-dependent increase in T_m1 was observed with divalent ions and triplex DNA stabilization decreased in the order: Mg²⁺ > Ca²⁺ > Sr²⁺ > Ba²⁺. All positively charged cations increased the melting temperature of duplex DNA. Values of Δn (number of ions released) on triplex DNA melting were 0.46 ± 0.06 and 0.18 ± 0.02, respectively, for mono- and di-valent cations, as calculated from 1/T_m1 versus ln[M^+,2+] plots. The corresponding values for duplex DNA were 0.25 ± 0.02 and 0.12 ± 0.02, respectively, for mono- and di-valent cations. Circular dichroism spectroscopic studies showed distinct conformational changes in triplex DNA stabilized by alkali metal and ammonium ions. Our results might be useful in developing triplex forming oligonucleotide based gene silencing techniques.     

The effects of fasting on swimming performance in juvenile qingbo (Spinibarbus sinensis) at two temperatures

We measured the following variables to investigate the effects of fasting and temperature on swimming performance in juvenile qingbo (Spinibarbus sinensis): the critical swimming speed (U_crit), resting metabolic rate ${(\dot{\mathrm{M}}O}_{\mathrm{2rest}})$ and active metabolic rate ${(\dot{\mathrm{M}}O}_{2active})$ of fish fasting for 0 (control), 1, 2 and 4 weeks at low and high acclimation temperatures (15 and 25 °C). Both fasting treatment and temperature acclimation had significant effects on all parameters measured (P<0.05). Fasting at the higher temperature had a negative effect on all measured parameters after 1 week (P<0.05). However, when acclimated to the lower temperature, fasting had a negative effect on U_crit until week 2 and on ${(\dot{\mathrm{M}}O}_{\mathrm{2rest}})$, ${(\dot{\mathrm{M}}O}_{\mathrm{2active}})$ and metabolic scope (MS, ${(\dot{\mathrm{M}}O}_{\mathrm{2active}})$–${(\dot{\mathrm{M}}O}_{\mathrm{2rest}})$) until week 4 (P<0.05). The values of all parameters at the lower temperature were significantly lower than those at the higher temperature in the identical fasting period groups except for ${(\dot{\mathrm{M}}O}_{\mathrm{2rest}})$ of the fish that fasted for 2 weeks. The relationship between fasting time (T) and U_crit was described as U_crit(15)=−0.302T²−0.800T+35.877 (r=0.781, n=32, P<0.001) and U_crit(25)=0.471T²−3.781T+50.097 (r=0.766, n=32, P<0.001₎ at 15 and 25 °C, respectively. The swimming performance showed less decrease in the early stage of fasting but more decrease in the later stage at the low temperature compared to the high temperature, which might be related to thermal acclimation time, resting metabolism, respiratory capacity, energy stores, enzyme activity in muscle tissue and energy substrate utilization changes with fasting between low and high temperatures. The divergent response of the swimming performance to fasting in qingbo at different temperatures might be an adaptive strategy to seasonal temperature and food resource variation in their habitat.     

The limitations,dangers, and benefits of simple methods for testing identifiability

In their Commentary paper, Villaverde and Massonis (On testing structural identifiability by a simple scaling method: relying on scaling symmetries can be misleading) have commented on our paper in which we proposed a simple scaling method to test structural identifiability. Our scaling invariance method (SIM) tests for scaling symmetries only, and Villaverde and Massonis correctly show the SIM may fail to detect identifiability problems when a model has other types of symmetries. We agree with the limitations raised by these authors but, also, we emphasize that the method is still valuable for its applicability to a wide variety of models, its simplicity, and even as a tool to introduce the problem of identifiability to investigators with little training in mathematics.

In their Commentary paper, Villaverde and Massonis (On testing structural identifiability by a simple scaling method: relying on scaling symmetries can be misleading [1]) have commented on our paper in which we proposed a simple scaling method to test structural identifiability [2]. Our scaling invariance method (SIM) tests for scaling symmetries only, and Villaverde and Massonis correctly show the SIM may fail to detect identifiability problems when a model has other types of symmetries (we indeed indicated but not investigated the importance of generalizing the method to other symmetries). Thus, we agree that our simple method provides a necessary but not sufficient condition for identifiability, and we appreciate their careful analysis and constructive criticism.We nevertheless think that the simple method remains useful because it is so simple. Even for investigators with little training in mathematics, the method provides a necessary condition for structural identifiability that can be derived in a few minutes with pen and paper. Similarly, we have found its pedagogic strength by teaching the method to our own graduate students and colleagues. More advanced methods (such as STRIKE-GOLDD [3,4], COMBOS [5], or SIAN [6]) are typically intimidating for researchers with a background in Biology or Bioinformatics. This simple method can help those practitioners to familiarize themselves with the identifiability problem and better understand their models.Finally, it is worth noting that if scaling invariance is the only symmetry (as it was in all the cases we analyzed), our SIM remains valuable (albeit uncontrolled), and surprisingly effective for a wide variety of problems (as the extensive list collected in the Supplementary Material our paper [2]). We guess that the SIM especially fails when applied to linear models (as more potential rotations of the variables leave the system invariant), and in non-linear scenarios where some parameters are identical. For instance, the FitzHugh-Nagumo model raised by Villaverde and Massonis, x˙1(t)=c(x1(t)−x13(t)3−x2(t)+d),x˙2(t)=1c(x1(t)+a−b·x2(t)),y(t)=x1(t), could have been written as x˙1(t)=λ1x1(t)−λ2x13(t)3−λ3x2(t)+d,x˙2(t)=λ4x1(t)+a−b·x2(t),y(t)=x1(t) where λ₁ = λ₂ = λ₃ = 1/λ₄ = c. One of the reasons why our method fails, in this case, might be these additional symmetries introduced in this more elaborate notation of the model.Hence, it is worth understanding generic conditions under which the SIM method is expected to be fragile, possibly using STRIKE-GOLDD to test large families of nonlinear models.As a final remark, we appreciate that Villaverde and Massonis have shared their source code, so researchers might have a gold standard to test identifiability.     

Fast Label-Free Cytoskeletal Network Imaging in Living Mammalian Cells

We present a full-field technique that allows label-free cytoskeletal network imaging inside living cells. This noninvasive technique allows monitoring of the cytoskeleton dynamics as well as interactions between the latter and organelles on any timescale. It is based on high-resolution quantitative phase imaging (modified Quadriwave lateral shearing interferometry) and can be directly implemented using any optical microscope without modification. We demonstrate the capability of our setup on fixed and living Chinese hamster ovary cells, showing the cytoskeleton dynamics in lamellipodia during protrusion and mitochondria displacement along the cytoskeletal network. In addition, using the quantitative function of the technique, along with simulation tools, we determined the refractive index of a single tubulin microtubule to be n_tubu=2.36±0.6$n_{tubu}=2.36\pm 0.6$ at λ=527$\lambda =527$ nm.     

Biaxial mechanical testing of human sclera

The biomechanical environment of the optic nerve head (ONH), of interest in glaucoma, is strongly affected by the biomechanical properties of sclera. However, there is a paucity of information about the variation of scleral mechanical properties within eyes and between individuals. We thus used biaxial testing to measure scleral stiffness in human eyes. Ten eyes from 5 human donors (age 55.4±3.5 years; mean±SD) were obtained within 24 h of death. Square scleral samples (6 mm on a side) were cut from each ocular quadrant 3–9 mm from the ONH centre and were mechanically tested using a biaxial extensional tissue tester (BioTester 5000, CellScale Biomaterials Testing, Waterloo). Stress–strain data in the latitudinal (toward the poles) and longitudinal (circumferential) directions, here referred to as directions 1 and 2, were fit to the four-parameter Fung constitutive equation W=c(e^Q?1), where $Q=c_1{E}_{11}^2+c_2{E}_{22}^2+2c_3{E}_{11}{E}_{22}$ and W, c’s and E_ij are the strain energy function, material parameters and Green strains, respectively. Fitted material parameters were compared between samples. The parameter c₃ ranged from 10^?7 to 10^?8, but did not contribute significantly to the accuracy of the fitting and was thus fixed at 10^?7. The products c?c₁ and c?c₂, measures of stiffness in the 1 and 2 directions, were 2.9±2.0 and 2.8±1.9 MPa, respectively, and were not significantly different (two-sided t-test; p=0.795). The level of anisotropy (ratio of stiffness in orthogonal directions) was 1.065±0.33. No statistically significant correlations between sample thickness and stiffness were found (correlation coefficients=?0.026 and ?0.058 in directions 1 and 2, respectively). Human sclera showed heterogeneous, near-isotropic, nonlinear mechanical properties over the scale of our samples.     

Graphic method for determination of allosteric constants

The maximum slope of the plot, appearing in the paper of Watari & Isogai (1976), was derived algebraically as a function of allosteric constants c and α_mor β_m (= cα_m), and the relation between L, c, and α_mor β_m, was also obtained, where $\text{L =}\text{T}{}_{\mathrm{o}}\text{R}{}_{\mathrm{o}}\text{, c =}\text{K}{}_{\mathrm{R}}\text{K}{}_{\mathrm{T}}\text{, α}{}_{\mathrm{m}}\text{=}\text{F}{}_{\mathrm{m}}\text{K}{}_{\mathrm{R}}\text{, β}{}_{\mathrm{m}}\text{=}\text{F}{}_{\mathrm{m}}\text{K}{}_{\mathrm{T}}\text{, R}{}_{\mathrm{o}}\text{and}\text{T}{}_{\mathrm{o}}$ are concentrations of unligated R and T states respectively, K_Rand K_T are microscopic dissociation constants, and F_m is the ligand concentration at the maximum slope of the plot. When the maximum slope is increased by one, the value becomes Hill constant, n. Nomographs which enable easier estimation of allosteric constants, L and c, were constructed from the two given values, the maximum slope of the plot, n ? 1, and α_mor β_m, in the cases where the maximum number of ligands, N, was 2 and 4. In the nomograph, log c is plotted against log L²c^N keeping the value of the maximum slope of the plot and that of α_mor β_m constant. These nomographs show that the representation is symmetrical in the cases of L²c^N > 1 and L²c^N < 1.     

Kinetics and Energetics of Biomolecular Folding and Binding

The ability of biomolecules to fold and to bind to other molecules is fundamental to virtually every living process. Advanced experimental techniques can now reveal how single biomolecules fold or bind against mechanical force, with the force serving as both the regulator and the probe of folding and binding transitions. Here, we present analytical expressions suitable for fitting the major experimental outputs from such experiments to enable their analysis and interpretation. The fit yields the key determinants of the folding and binding processes: the intrinsic on-rate and the location and height of the activation barrier.Dynamic processes in living cells are regulated through conformational changes in biomolecules—their folding into a particular shape or binding to selected partners. The ability of biomolecules to fold and to bind enables them to act as switches, assembly factors, pumps, or force- and displacement-generating motors (1). Folding and binding transitions are often hindered by a free energy barrier. Overcoming the barrier requires energy-demanding rearrangements such as displacing water from the sites of native contacts and breaking nonnative electrostatic contacts, as well as loss of configurational entropy. Once the barrier is crossed, the folded and bound states are stabilized by short-range interactions: hydrogen bonds, favorable hydrophobic effects, and electrostatic and van der Waals attractions (2).Mechanistic information about folding and binding processes is detailed in the folding and binding trajectories of individual molecules: observing an ensemble of molecules may obscure the inherent heterogeneity of these processes. Single-molecule trajectories can be induced, and monitored, by applying force to unfold/unbind a molecule and then relaxing the force until folding or binding is observed (3–5) (Fig. 1). Varying the force relaxation rate shifts the range of forces at which folding or binding occurs, thus broadening the explorable spectrum of molecular responses to force and revealing conformational changes that are otherwise too fast to detect. The measured force-dependent kinetics elucidates the role of force in physiological processes (6) and provides ways to control the timescales, and even the fate, of these processes. The force-dependent data also provides a route to understanding folding and binding in the absence of force—by extrapolating the data to zero force via a fit to a theory.Open in a separate windowFigure 1Schematic of the output from a force-relaxation experiment. The applied force is continuously relaxed from the initial value F₀ until the biomolecule folds or binds, as signified by a sharp increase in the measured force. From multiple repeats of this experiment, distributions of the folding or binding forces are collected (inset). Fitting the force distributions with the derived analytical expression yields the key parameters that determine the kinetics and energetics of folding or binding.In this letter, we derive an analytical expression for the distribution of transition forces, the major output of force-relaxation experiments that probe folding and binding processes. The expression extracts the key determinants of these processes: the on-rate and activation barrier in the absence of force. The theory is first developed in the context of biomolecular folding, and is then extended to cover the binding of a ligand tethered to a receptor. In contrast to unfolding and unbinding, the reverse processes of folding and binding require a theory that accounts for the compliance of the unfolded state, as well as the effect of the tether, to recover the true kinetic parameters of the biomolecule of interest.In a force-relaxation experiment, an unfolded biomolecule or unbound ligand-receptor complex is subject to a stretching force, which is decreased from the initial value F₀ as the pulling device approaches the sample at speed V until a folding or binding transition is observed (Fig. 1) (3–5). Define S(t) as the probability that the molecule has not yet escaped from the unfolded (implied: or unbound) state at time t. When escape is limited by one dominant barrier, S(t) follows the first-order rate equationS˙(t)≡dS(t)dt=−k←(F(t))S(t),where k_←(F(t)) is the on-rate at force F at time t. Because, prior to the transition, the applied force decreases monotonically with time, the distribution of transition forces, p(F), is related to S(t) through p(F)dF=S˙(t)dt, yieldingp(F)=−k←(F)F˙(F)e−∫F0Fk←(F′)F˙(F′)dF′.(1)Here F˙(F)≡dF(t)/dt<0 is the force relaxation rate. The proper normalization of p(F) is readily confirmed by integrating Eq. 1 from the initial force F₀ to negative infinity, the latter accounting for transitions that do not occur by the end of the experiment. Note that the expression for the distribution of folding/binding forces in Eq. 1 differs from its analog for the unfolding process (7) by the limits of integration and a negative sign, reflecting the property of a relaxation experiment to decrease the survival probability S(t) by decreasing the force. Converting the formal expression in Eq. 1 into a form suitable for fitting experimental data requires establishing functional forms for k_←(F) and F˙(F) and analytically solving the integral. These steps are accomplished below.The on-rate k_←(F) is computed by treating the conformational dynamics of the molecule as a random walk on the combined free energy profile G(x,t) = G₀(x) + G_pull(x,t) along the molecular extension x. Here G₀(x) is the intrinsic molecular potential and G_pull(x,t) is the potential of the pulling device. When G(x,t) features a high barrier on the scale of k_BT (k_B is the Boltzmann constant and T the temperature), the dynamics can be treated as diffusive. The unfolded region of the intrinsic potential for a folding process, unlike that for a barrierless process (8), can be captured by the functionG0(x)=ΔG‡ν1−ν(xx‡)11−ν−ΔG‡ν(xx‡),which has a sharp (if ν = 1/2, Fig. 2, inset) or smooth (if ν = 2/3) barrier of height ΔG^‡ and location x^‡. The potential of a pulling device of stiffness κ_S is G_pull(x,t) = κ_S/2(X₀ – Vt – x)² with an initial minimum at X₀ (corresponding to F₀). Applying Kramers formalism (9) to the combined potential G(x,t), we establish the analytical form of the on-rate at force F(t),k←(F)=k0(1+κSκU(F))1ν−12(1+νFx‡ΔG‡)1ν−1×eβΔG‡[1−(1+κSκU(F))2ν1−ν−1(1+νFx‡ΔG‡)1ν],where k₀ is the intrinsic on-rate, β ≡ (k_BT)⁻¹, andκU(F)=ν(1−ν)2ΔG‡x‡2(1+νFx‡ΔG‡)2−1νis the stiffness of the unfolded biomolecule under force F (see the Supporting Material for details on all derivations). The full nonlinear form of G_pull(x,t) was necessary in the derivation because, in contrast to the typically stiff folded state, the unfolded state may be soft (to be exact, 1/2κ_S x^‡2(F) << k_BT may not be satisfied) and thus easily deformed by the pulling device. Because of this deformation, the folding transition faces an extra contribution (regulated by the ratio κ_S/κ_U(F)) to the barrier height, typically negligible for unfolding, that decreases the on-rate in addition to the applied force F.Open in a separate windowFigure 2Contributions to the free energy profile for folding (inset) and binding (main figure). The derived expression (Eq. 2) extracts the on-rate and the location and height of the activation barrier to folding. When applied to binding data, the expression extracts the parameters of the ligand-tether-receptor (LTR) potential G˜0 (x); the proposed algorithm (Eqs. 3 and 4) removes the contribution of the tether potential G_teth(x) to recover the parameters of the intrinsic ligand-receptor (LR) potential G₀(x).The last piece required for Eq. 1, the loading rate F˙(F), is computed as the time derivative of the force F(t) on the unfolded molecule at its most probable extension at time t:F˙(F)=−κSV1+κS/κU(F).Finally, we realize that the integral in Eq. 1 can be solved analytically exactly, both for ν = 1/2 and ν = 2/3, resulting in the analytical expression for the distribution of folding forces:p(F)=k←(F)|F˙(F)|e−k←(F)β|F˙(F)|x‡(1+κSκU(F))νν−1(1+νFx‡ΔG‡)1−1ν.(2)Equation 2 can be readily applied to (normalized) histograms from force-relaxation experiments to extract the parameters of the intrinsic kinetics and energetics of folding. Being exact for ν = 1/2 and ν = 2/3, Eq. 2 is also an accurate approximation for any ν in the interval 1/2 < ν < 2/3 as long as κ_S ≲ κ_U (F) (see Fig. S1 in the Supporting Material). For simplicity, in Eq. 2 we have omitted the term containing F₀ as negligible if F₀ is large enough to prevent folding events.The solution in Eq. 2 reveals properties of the distribution of folding forces that distinguish it from its unfolding counterpart (7):

• 1.The distribution has a positive skew (Fig. 3), as intuitively expected: the rare folding events occur at high forces when the barrier is still high.Open in a separate windowFigure 3Force histograms from folding (left) and binding (right) simulations at several values of the force-relaxation speed (in nanometers per second, indicated at each histogram). Fitting the histograms with the analytical expression in Eq. 2 (lines) recovers the on-rate and activation barrier for folding or binding (2.Increasing the relaxation speed shifts the distribution to lower forces (Fig. 3): faster force relaxation leaves less time for thermal fluctuations to push the system over a high barrier, causing transitions to occur later (i.e., at lower forces), when the barrier is lower.
• 3.The stiffness κ_S and speed V enter Eq. 2 separately, providing independent routes to control the range of folding forces and thus enhance the robustness of a fit.

The application of the above framework to binding experiments on a ligand and receptor connected by a tether (3) involves an additional step—decoupling the effect of the tether—to reconstruct the parameters of ligand-receptor binding. Indeed, the parameters extracted from a fit of experimental histograms to Eq. 2 characterize the ligand-tether-receptor (LTR) potential (k˜0, x˜‡, ΔG˜‡, ν) (Fig. 2). The parameters of the natural ligand-receptor (LR) potential (k₀, x^‡, ΔG^‡) can be recovered using three characteristics of the tether: contour length L; persistence length p; and extension Δℓ of the tether along the direction of the force in the LTR transition state. The values of L and p can be determined from the force-extension curve of the tether (10); these define the tether potential G_teth(x) (Fig. 2). The value of Δℓ can be found from an unbinding experiment (7) on LTR and the geometry of the tether attachment points (see Fig. S3). Approximating the region of the LR potential between the transition and unbound states as harmonic, with no assumptions about the shape of the potential beyond x^‡, the ligand-receptor barrier parameters are thenx‡=α−1α−2x˜‡,ΔG‡=(α−1)22(α−2)x˜‡Fteth(Δℓ+x˜‡),(3)and the intrinsic unimolecular association rate isk0≈k˜0(βΔG‡)32(βΔG˜‡)1ν−12(x˜‡x‡)2eβ(ΔG˜‡−ΔG‡).(4)Here, the force value Fteth(Δℓ+x˜‡) is extracted from the force-extension curve of the tether at extension Δℓ+x˜‡ andα=2(ΔG˜‡−Gteth(Δℓ)+Gteth(Δℓ+x˜‡))x˜‡Fteth(Δℓ+x˜‡),where G_teth(x) is the wormlike-chain potential (see Eq. S13 in the Supporting Material). Equations 3–4 confirm that a tether decreases the height and width of the barrier (see Fig. 2), thus increasing the on-rate.In Fig. 3, the developed analytical framework is applied to folding and binding force histograms from Brownian dynamics simulations at parameters similar to those in the analogous experimental and computational studies (3,5,11) (for details on simulations and fitting procedure, see the Supporting Material). For the stringency of the test, the simulations account for the wormlike-chain nature of the molecular unfolded and LTR unbound states that is not explicitly accounted for in the theory. With optimized binning (12) of the histograms and a least-squares fit, Eqs. 2–4 recover the on-rate, the location and the height of the activation barrier, and the value of ν that best captures how the kinetics scale with force (

1.Multiple relaxation speeds,

2.Folding/binding events at low forces, and

3.A large number of events at each speed.

Table 1

On-rate and the location and height of the activation barrier from the fit of simulated data to the theory in Eq. 2

Kinetics of H2O2-driven catalysis by a lytic polysaccharide monooxygenase from the fungus Trichoderma reesei

Owing to their ability to break glycosidic bonds in recalcitrant crystalline polysaccharides such as cellulose, the catalysis effected by lytic polysaccharide monooxygenases (LPMOs) is of major interest. Kinetics of these reductant-dependent, monocopper enzymes is complicated by the insoluble nature of the cellulose substrate and parallel, enzyme-dependent, and enzyme-independent side reactions between the reductant and oxygen-containing cosubstrates. Here, we provide kinetic characterization of cellulose peroxygenase (oxidative cleavage of glycosidic bonds in cellulose) and reductant peroxidase (oxidation of the reductant) activities of the LPMO TrAA9A of the cellulose-degrading model fungus Trichoderma reesei. The catalytic efficiency (kcat/Km(H2O2)) of the cellulose peroxygenase reaction (k_cat = 8.5 s⁻¹, and Km(H2O2)=30μM) was an order of magnitude higher than that of the reductant (ascorbic acid) peroxidase reaction. The turnover of H₂O₂ in the ascorbic acid peroxidase reaction followed the ping-pong mechanism and led to irreversible inactivation of the enzyme with a probability of 0.0072. Using theoretical analysis, we suggest a relationship between the half-life of LPMO, the values of kinetic parameters, and the concentrations of the reactants.