C\bar lr and C\bar ls subcomponents of human complement: Two serine proteinases lacking the ‘histidine-loop’ disulphide bridge

The N-terminal amino acid sequence of human C \(\bar l\) s b chain has been extended to 52 residues. The histidine residue involved in the charge-relay system is located at position 38, whereas the ‘histidine-loop’ disulphide bridge is missing. So far, human complement subcomponents C \(\bar l\) r and C \(\bar l\) s are the only known mammalian serine proteinases lacking this disulphide bridge.     

Modulation of the cAMP Response by G\alpha _i and G\beta \gamma : A Computational Study of G Protein Signaling in Immune Cells

Cyclic AMP is important for the resolution of inflammation, as it promotes anti-inflammatory signaling in several immune cell lines. In this paper, we present an immune cell specific model of the cAMP signaling cascade, paying close attention to the specific isoforms of adenylyl cyclase (AC) and phosphodiesterase that control cAMP production and degradation, respectively, in these cells. The model describes the role that G protein subunits, including G \(\alpha _s\) , G \(\alpha _i\) , and G \(\beta \gamma \) , have in regulating cAMP production. Previously, G \(\alpha _i\) activation has been shown to increase the level of cAMP in certain immune cell types. This increase in cAMP is thought to be mediated by \(\beta \gamma \) subunits which are released upon G \(\alpha \) activation and can directly stimulate specific isoforms of AC. We conduct numerical experiments in order to explore the mechanisms through which G \(\alpha _i\) activation can increase cAMP production. An important conclusion of our analysis is that the relative abundance of different G protein subunits is an essential determinant of the cAMP profile in immune cells. In particular, our model predicts that limited availability of \(\beta \gamma \) subunits may both \((i)\) enable immune cells to link inflammatory G \(\alpha _i\) signaling to anti-inflammatory cAMP production thereby creating a balanced immune response to stimulation with low concentrations of PGE2, and \((ii)\) prohibit robust anti-inflammatory cAMP signaling in response to stimulation with high concentrations of PGE2.     

Kinetic models of coupling between H+ and Na+-translocation and ATP synthesis/hydrolysis by F0F1-ATPases: Can a cell utilize both\Delta \bar \mu _{H^ + } and\Delta \bar \mu _{Na^ + } for ATP synthesis underin vivo conditions using the same enzyme?

Kinetic models of the F₀F₁-ATPase able to transport H⁺ or/and Na⁺ ions are proposed. It is assumed that (i) H⁺ and Na⁺ compete for the same binding sites, (ii) ion translocation through F₀ is coupled to the rate-limiting step of the F₁-catalyzed reaction. The main characteristics of the dependences of ATP synthesis and hydrolysis rates on Δφ, ΔpH, and ΔpNa are predicted for various versions of the coupling model. The mechanism of the switchover from \(\Delta \bar \mu _{H^ + } \) -dependent synthesis to the \(\Delta \bar \mu _{Na^ + } \) -dependent one is demonstrated. It is shown that even with a drastic drop in \(\Delta \bar \mu _{H^ + } \) , ATP hydrolysis by the proton mode of catalysis can be effectively inhibited by Δφ and ΔpNa. The results obtained strongly support the possibility that the same F₀F₁-ATPase in bacterial cells can utilize both \(\Delta \bar \mu _{H^ + } \) and \(\Delta \bar \mu _{Na^ + } \) for ATP synthesis underin vivo conditions.     

A thermomechanical framework for reconciling the effects of ultraviolet radiation exposure time and wavelength on connective tissue elasticity

Augmentation of the mechanical properties of connective tissue using ultraviolet (UV) radiation—by targeting collagen cross-linking in the tissue at predetermined UV exposure time \((t)\) and wavelength \((\lambda )\) —has been proposed as a therapeutic method for supporting the treatment for structural-related injuries and pathologies. However, the effects of \(\lambda \) and \(t\) on the tissue elasticity, namely elastic modulus \((E)\) and modulus of resilience \((u_\mathrm{Y})\) , are not entirely clear. We present a thermomechanical framework to reconcile the \(t\) - and \(\lambda \) -related effects on \(E\) and \(u_\mathrm{Y}\) . The framework addresses (1) an energy transfer model to describe the dependence of the absorbed UV photon energy, \(\xi \) , per unit mass of the tissue on \(t\) and \(\lambda \) , (2) an intervening thermodynamic shear-related parameter, \(G\) , to quantify the extent of UV-induced cross-linking in the tissue, (3) a threshold model for the \(G\) versus \(\xi \) relationship, characterized by   \(t_\mathrm{C}\) —the critical \(t\) underpinning the association of \(\xi \) with \(G\) —and (4) the role of \(G\) in the tissue elasticity. We hypothesized that \(G\) regulates \(E\) (UV-stiffening hypothesis) and \(u_\mathrm{Y}\) (UV-resilience hypothesis). The framework was evaluated with the support from data derived from tensile testing on isolated ligament fascicles, treated with two levels of \(\lambda \) (365 and 254 nm) and three levels of \(t\) (15, 30 and 60 min). Predictions from the energy transfer model corroborated the findings from a two-factor analysis of variance of the effects of \(t\) and \(\lambda \) treatments. Student’s t test revealed positive change in \(E\) and \(u_\mathrm{Y}\) with increases in \(G\) —the findings lend support to the hypotheses, implicating the implicit dependence of UV-induced cross-links on \(t\) and \(\lambda \) for directing tissue stiffness and resilience. From a practical perspective, the study is a step in the direction to establish a UV irradiation treatment protocol for effective control of exogenous cross-linking in connective tissues.     

Quantifying interspecific spatial overlap in aquatic macrophyte communities

Levins’s asymmetrical α index quantifies between species overlap over resources more realistically than similar-purpose single-value indices. The associated community-wide \(\bar \alpha\) index expresses the degree of “species packing”. Both indices were formulated upon competing animal (i.e., mobile) organisms and are independent of population densities. However, overlap over resources for nonmobile organisms such as plants may have an impact even below carrying capacity. The proposed \(\hat \alpha\) index, based on Levins’s α index, quantifies spatial overlap for plants integrating information on species spatial distribution and crowding conditions. The \(\hat \alpha\) index is specifically designed for plant distribution data collected in discrete plots with density expressed as percent coverage (%cover) of substratum. We also propose a community-wide \({\hat \alpha_{\text{c}}}\) index, conceptually analogous to \(\bar \alpha\) , but furnished with a measure of dispersion (se \({\hat \alpha_{\text{c}}}\) ). Species importance within the community is inferred from comparisons of pairwise \(\hat \alpha\) ’s with \({\hat \alpha_{\text{c}}}\) . The \(\hat \alpha\) and \({\hat \alpha_{\text{c}}}\) indices correlate closely and exponentially with plant density, and correct apparent over- and underestimations of interaction intensity at low and very high crowding by Levins’s α and \(\bar \alpha\) , respectively. Index application to aquatic plant communities gave results consistent with within-community and general ecological patterns, suggesting a high potential of the proposed \(\hat \alpha\) and \({\hat \alpha_{\text{c}}}\) indices in basic and applied macrophyte ecological studies and management.     

Sur l'espace topologique lié à une nouvelle théorie de l'apprentissage

There have been two contrasting doctrines concerning learning, more generally about acquisition of knowledge: empiricism and rationalism. The theory of learning in such a field as artificial intelligence seems to fall within the empiricist framework. On the hand, N. Chomsky and his followers have discussed, during the last decade, concerning learning, especially about language learning, from the rationalist point of view (Chomsky, 1965). The main feature in the rationalist approach toward a theory of learning lies in the speculation that in order to acquire knowledge it is indispensable for a learner to be endowed with “innate ideas”, and that “experience” in the external world are merely subsidiary types of information for the learner. If this is acceptable, we can inquire: Under what kind of innate ideas can the learner understand the structure of the external world? In our previous paper (Uesaka, Aizawa, Ebara, and Ozeki, 1973), we formalized this by introducing the mathematical notion of “learnability”, and gave a partial answer to the above inquiry. In this formalization we assumed that the set F of objects to be learned consists of mappings of N to itself, where N is the set of positive integers. Then, constructing a topological space (F, \(\mathcal{O}\) ) by an appropriate family \(\mathcal{O}\) of open sets, we observed that the notion of learnability can be well described in terms of topological properties of the learning space (F, \(\mathcal{O}\) ). Many problems must be solved, however, before we raise the theory to a complete model of the rationalist theory of learning. The topological study of the space (F, \(\mathcal{O}\) ) is, we believe, the first step toward this approach. In this context, we discuss the topological aspects of this space. Now we define \(\mathcal{O}\) as follows: By N ² we mean the direct product of two N's. Let s be a subset of N ². If, for any (x, y), (x′, y′) in s, x=x′ implies y=y′, then we say that s is single-valued. Let f ∈ F, If, for any (x, y) in s, y=f(x), then f is said to be on s, denoted as \(f\underline \supseteq s\) . Let \(\pi \left( s \right) = \left\{ {g;g \in F,g\underline \supseteq s} \right\}\) . A single-valued finite subset of N ² is called datum. Let D denote the family of all data. Let \(\mathcal{O}* = \left\{ \phi \right\} \cup \left\{ {\pi \left( d \right);d \in D} \right\}\) , and \(\mathcal{O}\) denote the family of all subsets of F, each of which is written as \(\mathop \cup \limits_\alpha W_{\alpha }\) , where W _α is in \(\mathcal{O}*\) . Then, it is easily seen that \(\mathcal{O}\) satisfies the axiom of the open system of a topological space. It is shown that the learning space (F, \(\mathcal{O}\) ) has the following properties:

1. It satisfies the first and the second countability axioms.
2. It is separable and is totally disconnected.
3. It is a Hausdorff space and, further, is regular and normal.
4. It is neither compact nor locally compact.
5. It is metrizable, or more precisely there exists a complete but not totally bounded metric space which is homeomorphic to learning space.
6. Any of its subspace can be embedded into its special subspace.

     

Investigation of the possibility of exceeding the Greenwald density limit during ECRH in T-10

In T-10 experiments, attempts were made to significantly exceed the Greenwald limit $\bar n_{Gr} $ during high-power (P _ab=750 kW) electron-cyclotron resonance heating (ECRH) and gas puffing. Formally, the density limit $(\bar n_e )_{\lim } $ exceeding the Greenwald limit $({{(\bar n_e )_{\lim } } \mathord{\left/ {\vphantom {{(\bar n_e )_{\lim } } {\bar n_{Gr} }}} \right. \kern-0em} {\bar n_{Gr} }} = 1.8)$ was achieved for q _L=8.2. However, as q _L decreased, the ratio ${{(\bar n_e )_{\lim } } \mathord{\left/ {\vphantom {{(\bar n_e )_{\lim } } {\bar n_{Gr} }}} \right. \kern-0em} {\bar n_{Gr} }}$ also decreased, approaching unity at q _L≈3. It was suggested that the “current radius” (i.e., the radius of the magnetic surface enclosing the bulk of the plasma current I _p), rather than the limiter radius, was the parameter governing the value of $(\bar n_e )_{\lim } $ . In the ECRH experiments, no substantial degradation of plasma confinement was observed up to $\bar n_e \sim 0.9(\bar n_e )_{\lim } $ regardless of the ratio ${{(\bar n_e )_{\lim } } \mathord{\left/ {\vphantom {{(\bar n_e )_{\lim } } {\bar n_{Gr} }}} \right. \kern-0em} {\bar n_{Gr} }}$ . In different scenarios of the density growth up to $(\bar n_e )_{\lim } $ , two types of disruptions related to the density limit were observed.     

Approximation of sojourn-times via maximal couplings: motif frequency distributions

Sojourn-times provide a versatile framework to assess the statistical significance of motifs in genome-wide searches even under non-Markovian background models. However, the large state spaces encountered in genomic sequence analyses make the exact calculation of sojourn-time distributions computationally intractable in long sequences. Here, we use coupling and analytic combinatoric techniques to approximate these distributions in the general setting of Polish state spaces, which encompass discrete state spaces. Our approximations are accompanied with explicit, easy to compute, error bounds for total variation distance. Broadly speaking, if \({\mathsf{T}}_n\) is the random number of times a Markov chain visits a certain subset \({\mathsf{T}}\) of states in its first \(n\) transitions, then we can usually approximate the distribution of \({\mathsf{T}}_n\) for \(n\) of order \((1-\alpha )^{-m}\) , where \(m\) is the largest integer for which the exact distribution of \({\mathsf{T}}_m\) is accessible and \(0\le \alpha \le 1\) is an ergodicity coefficient associated with the probability transition kernel of the chain. This gives access to approximations of sojourn-times in the intermediate regime where \(n\) is perhaps too large for exact calculations, but too small to rely on Normal approximations or stationarity assumptions underlying Poisson and compound Poisson approximations. As proof of concept, we approximate the distribution of the number of matches with a motif in promoter regions of C. elegans. Mathematical properties of the proposed ergodicity coefficients and connections with additive functionals of homogeneous Markov chains as well as ergodicity of non-homogeneous Markov chains are also explored.     

N-site Phosphorylation Systems with 2N-1 Steady States

Multisite protein phosphorylation plays a prominent role in intracellular processes like signal transduction, cell-cycle control and nuclear signal integration. Many proteins are phosphorylated in a sequential and distributive way at more than one phosphorylation site. Mathematical models of \(n\) -site sequential distributive phosphorylation are therefore studied frequently. In particular, in Wang and Sontag (J Math Biol 57:29–52, 2008), it is shown that models of \(n\) -site sequential distributive phosphorylation admit at most \(2n-1\) steady states. Wang and Sontag furthermore conjecture that for odd \(n\) , there are at most \(n\) and that, for even \(n\) , there are at most \(n+1\) steady states. This, however, is not true: building on earlier work in Holstein et al. (Bull Math Biol 75(11):2028–2058, 2013), we present a scalar determining equation for multistationarity which will lead to parameter values where a \(3\) -site system has \(5\) steady states and parameter values where a \(4\) -site system has \(7\) steady states. Our results therefore are counterexamples to the conjecture of Wang and Sontag. We furthermore study the inherent geometric properties of multistationarity in \(n\) -site sequential distributive phosphorylation: the complete vector of steady state ratios is determined by the steady state ratios of free enzymes and unphosphorylated protein and there exists a linear relationship between steady state ratios of phosphorylated protein.     

A locally convergent rotationally invariant particle swarm optimization algorithm

Several well-studied issues in the particle swarm optimization algorithm are outlined and some earlier methods that address these issues are investigated from the theoretical and experimental points of view. These issues are the: stagnation of particles in some points in the search space, inability to change the value of one or more decision variables, poor performance when the swarm size is small, lack of guarantee to converge even to a local optimum (local optimizer), poor performance when the number of dimensions grows, and sensitivity of the algorithm to the rotation of the search space. The significance of each of these issues is discussed and it is argued that none of the particle swarm optimizers we are aware of can address all of these issues at the same time. To address all of these issues at the same time, a new general form of velocity update rule for the particle swarm optimization algorithm that contains a user-definable function \(f\) is proposed. It is proven that the proposed velocity update rule guarantees to address all of these issues if the function \(f\) satisfies the following two conditions: (i) the function \(f\) is designed in such a way that for any input vector \(\vec {y}\) in the search space, there exists a region \(A\) which contains \(\vec {y}\) and \( f\!\left( {\vec {y}} \right) \) can be located anywhere in \(A\) , and (ii) \(f\) is invariant under any affine transformation. An example of function \(f\) is provided that satisfies these conditions and its performance is examined through some experiments. The experiments confirm that the proposed algorithm (with an appropriate function \(f)\) can effectively address all of these issues at the same time. Also, comparisons with earlier methods show that the overall ability of the proposed method for solving benchmark functions is significantly better.     

Individual-Based Model for Quorum Sensing with Background Flow

Quorum sensing is a wide-spread mode of cell–cell communication among bacteria in which cells release a signalling substance at a low rate. The concentration of this substance allows the bacteria to gain information about population size or spatial confinement. We consider a model for \(N\) cells which communicate with each other via a signalling substance in a diffusive medium with a background flow. The model consists of an initial boundary value problem for a parabolic PDE describing the exterior concentration \(u\) of the signalling substance, coupled with \(N\) ODEs for the masses \(a_i\) of the substance within each cell. The cells are balls of radius \(R\) in \(\mathbb {R} ^3\) , and under some scaling assumptions we formally derive an effective system of \(N\) ODEs describing the behaviour of the cells. The reduced system is then used to study the effect of flow on communication in general, and in particular for a number of geometric configurations.     

Age,smoking habits,heat stress,and their interactive effects with carbon monoxide and peroxyacetylnitrate on man's aerobic power

Metabolic, body temperature, and cardiorespiratory responses of 16 healthy middle-aged (40–57 years) men, 9 nonsmokers and 7 smokers, were obtained during tests of maximal aerobic power at ambient environmental temperatures of 25 ± 0.5 and 35 ± 0.5°C and 20% relative humidity under four conditions: (a) filtered air, FA; (b) 50 ppm carbon monoxide in filtered air, CO; (c) 0.27 ppm peroxyacetylnitrate in filtered air, PAN; and (d) a combination of all three mixtures, PANCO. There was no significant change in maximum aerobic power \(\left( {\dot VO2max} \right)\) related to the presence of air pollutants, although total working time was lowered in the 25°C environment while breathing CO. Older nonsmokers did have a decrement in \(\left( {\dot VO2max} \right)\) while breathing 50 ppm CO, while older smokers failed to show any change. This difference was related to the initial COHb levels of the smokers, who, when breathing this level of ambient CO, had only a 14% increase in COHb over their initial levels in contrast to the 200% increase in the nonsmokers. Smoking habits were the most influential factor affecting the cardiorespiratory responses of these older men to maximal exercise. Regardless of ambient conditions, smokers had a significantly lower (27%) aerobic power than nonsmokers, were breathing closer to their maximal breathing capacities throughout the walk, and had a higher respiratory exchange ratio. While the \(\left( {\dot VO2max} \right)\) of nonsmokers was only 6% less than that of younger nonsmoking males ( \(\bar x\) age = 25 years) working under similar conditions, the aerobic power of the older smokers was 26% lower than that of young smokers ( \(\bar x\) age = 24 years).     

Synchronization of non-linear biochemical oscillators coupled by diffusion

The behaviour of similar coupled non-linear oscillators of the type \(\dot x\) =f(x, y, µ \(\dot y\) =g(x, y, µ is to be investigated. The oscillators are assumed to be coupled by diffusion gradients. If some conditions on the magnitude of the diffusion coefficients are satisfied, it is proved that: 1) if the oscillators have the same period (identical value of the parameter μ) and different phases before coupling, after coupling they tend to synchronize the phases; 2) if the periods of the oscillators are not too different (in terms of the values of the parameter μ) before coupling, after coupling they tend to oscillate with the same period. It is suggested the possible role of diffusion as a synchronizing mechanism in some biological phenomena.     

On linear birth-and-death processes in a random environment

We study the probability of extinction for single-type and multi-type continuous-time linear birth-and-death processes in a finite Markovian environment. The probability of extinction is equal to 1 almost surely if and only if the basic reproduction number \(R_0\) is \(\le 1\) , the key point being to identify a suitable definition of \(R_0\) for such a result to hold.     

Elastic behavior of a red blood cell with the membrane’s nonuniform natural state: equilibrium shape,motion transition under shear flow,and elongation during tank-treading motion

Direct numerical simulations of the mechanics of a single red blood cell (RBC) were performed by considering the nonuniform natural state of the elastic membrane. A RBC was modeled as an incompressible viscous fluid encapsulated by an elastic membrane. The in-plane shear and area dilatation deformations of the membrane were modeled by Skalak constitutive equation, while out-of-plane bending deformation was formulated by the spring model. The natural state of the membrane with respect to in-plane shear deformation was modeled as a sphere ( \(\alpha =0\) ), biconcave disk shape ( \(\alpha =1\) ) and their intermediate shapes ( \(0<\alpha <1\) ) with the nonuniformity parameter \(\alpha \) , while the natural state with respect to out-of-plane bending deformation was modeled as a flat plane. According to the numerical simulations, at an experimentally measured in-plane shear modulus of \(2.5\times 10^{-6}\,\hbox {N}/\hbox {m}\) and an out-of-plane bending rigidity of \(2.0\times 10^{-19}\,\hbox {N}\cdot \hbox {m}\) of the cell membrane, the following results were obtained. (i) The RBC shape at equilibrium was biconcave discoid for \(\alpha >0.22\) and cupped otherwise; (ii) the experimentally measured fluid shear stress at the transition between tumbling and tank-treading motions under shear flow was reproduced for \(0.05<\alpha <0.34\) ; (iii) the elongation deformation of the RBC during tank-treading motion from the simulation was consistent with that from in vitro experiments, irrespective of the \(\alpha \) value. Based on our RBC modeling, the three phenomena (i), (ii), and (iii) were mechanically consistent for \(0.22<\alpha <0.34\) . The condition \(0.05<\alpha <0.22\) precludes a biconcave discoid shape at equilibrium (i); however, it gives appropriate fluid shear stress at the motion transition under shear flow (ii), suggesting that a combined effect of \(\alpha \) and the natural state with respect to out-of-plane bending deformation is necessary for understanding details of the RBC mechanics at equilibrium. Our numerical results demonstrate that moderate nonuniformity in a membrane’s natural state with respect to in-plane shear deformation plays a key role in RBC mechanics.     

Counting glycans revisited

We present an algorithm for counting glycan topologies of order \(n\) that improves on previously described algorithms by a factor \(n\) in both time and space. More generally, we provide such an algorithm for counting rooted or unrooted \(d\) -ary trees with labels or masses assigned to the vertices, and we give a “recipe” to estimate the asymptotic growth of the resulting sequences. We provide constants for the asymptotic growth of \(d\) -ary trees and labeled quaternary trees (glycan topologies). Finally, we show how a classical result from enumeration theory can be used to count glycan structures where edges are labeled by bond types. Our method also improves time bounds for counting alkanes.