Stress correlation function evolution in lattice solid elasto-dynamic models of shear and fracture zones and earthquake prediction

It has been argued that power-law time-to-failure fits for cumulative Benioff strain and an evolution in size-frequency statistics in the lead-up to large earthquakes are evidence that the crust behaves as a Critical Point (CP) system. If so, intermediate-term earthquake prediction is possible. However, this hypothesis has not been proven. If the crust does behave as a CP system, stress correlation lengths should grow in the lead-up to large events through the action of small to moderate ruptures and drop sharply once a large event occurs. However this evolution in stress correlation lengths cannot be observed directly. Here we show, using the lattice solid model to describe discontinuous elasto-dynamic systems subjected to shear and compression, that it is for possible correlation lengths to exhibit CP-type evolution. In the case of a granular system subjected to shear, this evolution occurs in the lead-up to the largest event and is accompanied by an increasing rate of moderate-sized events and power-law acceleration of Benioff strain release. In the case of an intact sample system subjected to compression, the evolution occurs only after a mature fracture system has developed. The results support the existence of a physical mechanism for intermediate-term earthquake forecasting and suggest this mechanism is fault-system dependent. This offers an explanation of why accelerating Benioff strain release is not observed prior to all large earthquakes. The results prove the existence of an underlying evolution in discontinuous elasto-dynamic, systems which is capable of providing a basis for forecasting catastrophic failure and earthquakes.

Long-range automaton models of earthquakes: Power-law accelerations, correlation evolution, and mode-switching

We introduce a conceptual model for the in-plane physics of an earthquake fault. The model employs cellular automaton techniques to simulate tectonic loading, earthquake rupture, and strain redistribution. The impact of a hypothetical crustal elastodynamic Green's function is approximated by a long-range strain redistribution law with a r(-p) dependance. We investigate the influence of the effective elastodynamic interaction range upon the dynamical behaviour of the model by conducting experiments with different values of the exponent (p). The results indicate that this model has two distinct, stable modes of behaviour. The first mode produces a characteristic earthquake distribution with moderate to large events preceeded by an interval of time in which the rate of energy release accelerates. A correlation function analysis reveals that accelerating sequences are associated with a systematic, global evolution of strain energy correlations within the system. The second stable mode produces Gutenberg-Richter statistics, with near-linear energy release and no significant global correlation evolution. A model with effectively short-range interactions preferentially displays Gutenberg-Richter behaviour. However, models with long-range interactions appear to switch between the characteristic and GR modes. As the range of elastodynamic interactions is increased, characteristic behaviour begins to dominate GR behaviour. These models demonstrate that evolution of strain energy correlations may occur within systems with a fixed elastodynamic interaction range. Supposing that similar mode-switching dynamical behaviour occurs within earthquake faults then intermediate-term forecasting of large earthquakes may be feasible for some earthquakes but not for others, in alignment with certain empirical seismological observations. Further numerical investigation of dynamical models of this type may lead to advances in earthquake forecasting research and theoretical seismology.

Circular proteins - no end in sight

Circular proteins are a recently discovered phenomenon. They presumably evolved to confer advantages over ancestral linear proteins while maintaining the intrinsic biological functions of those proteins. In general, these advantages include a reduced sensitivity to proteolytic cleavage and enhanced stability. In one remarkable family of circular proteins, the cyclotides, the cyclic backbone is additionally braced by a knotted arrangement of disulfide bonds that confers additional stability and topological complexity upon the family. This article describes the discovery, structure, function and biosynthesis of the currently known circular proteins. The discovery of naturally occurring circular proteins in the past few years has been complemented by new chemical and biochemical methods to make synthetic circular proteins; these are also briefly described.

Prediction of Golgi Type II membrane proteins based on their transmembrane domains

Motivation: A major issue in cell biology today is how distinct intracellular regions of the cell, like the Golgi Apparatus, maintain their unique composition of proteins and lipids. The cell differentially separates Golgi resident proteins from proteins that move through the organelle to other subcellular destinations. We set out to determine if we could distinguish these two types of transmembrane proteins using computational approaches. Results: A new method has been developed to predict Golgi membrane proteins based on their transmembrane domains. To establish the prediction procedure, we took the hydrophobicity values and frequencies of different residues within the transmembrane domains into consideration. A simple linear discriminant function was developed with a small number of parameters derived from a dataset of Type II transmembrane proteins of known localization. This can discriminate between proteins destined for Golgi apparatus or other locations (post-Golgi) with a success rate of 89.3% or 85.2%, respectively on our redundancy-reduced data sets.

The role of conserved CXXXRXG motif in the expression and function of the human norepinephrine transporter

Highly conserved motifs in the monoamine transporters, e.g. the human norepinephrine transporter (hNET) GXXXRXG motif which was the focus of the present study, are likely to be important structural features in determining function. This motif was investigated by mutating the glycines to glutamate (causing loss of function) and alanine, and the arginine to glycine. The effects of hG117A, hR121G and hG123A mutations on function were examined in COS-7 cells and compared to hNET. Substrate K-m values were decreased for hG117A and hG123A, and their K values for inhibition of [3 H]nisoxetine binding were decreased 3-4-fold and 4-6-fold, respectively. Transporter turnover was reduced to 65% of hNET for hG117A and hR121G and to 28% for hG123A, suggesting that substrate translocation is impaired. K values of nisoxetine and desipramine for inhibition of [H-3]norepinephrine uptake were increased by 5-fold for hG117A, with no change for cocaine. The K-i value of cocaine was increased by 3-fold for hG123A, with no change for nisoxetine and desipramine. However, there were no effects of the mutations on the K-d of [H-3]nisoxetine binding or K-i values of desipramine or cocaine for inhibition of [H-3]nisoxetine binding. Hence, glycine residues of the GXXXRXG motif are important determinants of NET expression and function, while the arginine residue does not have a major role. This study also showed that antidepressants and psychostimulants have different NET binding sites and provided the first evidence that different sites on the NET are involved in the binding of inhibitors and their competitive inhibition of substrate uptake. (C) 2002 Elsevier Science B.V. All rights reserved.

Morphology, ultrastructure, and implied function of ciliated sensory structures on the developmental stages of Merizocotyle icopae (Monogenea : Monocotylidae)

Experimental infections were used to track the fate of the dorsal sensilla of Merizocotyle icopae (Monogenea: Monocotylidae) from nasal tissue of the shovelnose ray, Rhinobatos typus (Rhinobatidae). Scanning and transmission electron microscopy revealed that 3 types of uniciliate dorsal sensilla exist at different times in the development of the monogenean. Type 1 sensilla have little or no invagination where the cilium exits the distal end of the dendrite and possess a ring of epidermis surrounding the cilium distal to the invagination. Type 2 sensilla have a deep invagination where the cilium exits the dendrite. Type 3 sensilla can be distinguished from the other types by the shape of the dendrite. The larvae have predominantly Type I dorsal sensilla, most of which are lost approximately 24 h after infection and a few Type 2 sensilla, which are retained. Additional Type 2 sensilla (termed Adult Type 2 sensilla), which are slightly different morphologically from the Type 2 sensilla of the larvae, form in later stages of development. Numerous Type 3 sensilla are unique to the dorsal surface of adults. Loss of all Type I sensilla upon attachment to the host, R. typus, suggests that these may be chemo- or mechanoreceptors responsible for host location by the swimming infective larvae. Type 2 sensilla appear to be important in the larvae, juveniles, and adults whereas the modality mediated by Type 3 is specific to adults. (C) 2003 Wiley-Liss, Inc.

Identification of residues and domains of Raf important for function in vivo and in vitro

Random mutagenesis and genetic screens for impaired Raf function in Caenorhabditis elegans were used to identify six loss-of-function alleles of lin-45 raf that result in a substitution of a single amino acid. The mutations were classified as weak, intermediate, and strong based on phenotypic severity. We engineered these mutations into the homologous residues of vertebrate Raf-1 and analyzed the mutant proteins for their underlying biochemical defects. Surprisingly, phenotype strength did not correlate with the catalytic activity of the mutant proteins. Amino acid substitutions Val-589 and Ser-619 severely compromised Raf kinase activity, yet these mutants displayed weak phenotypes in the genetic screen. Interestingly, this is because these mutant Raf proteins efficiently activate the MAPK (mitogen-activated protein kinase) cascade in living cells, a result that may inform the analysis of knockout mice. Equally intriguing was the observation that mutant proteins with non-functional Ras-binding domains, and thereby deficient in Ras-mediated membrane recruitment, displayed only intermediate strength phenotypes. This confirms that secondary mechanisms exist to couple Ras to Raf in vivo. The strongest phenotype in the genetic screens was displayed by a S508N mutation that again did not correlate with a significant loss of kinase activity or membrane recruitment by oncogenic Ras in biochemical assays. Ser-508 lies within the Raf-1 activation loop, and mutation of this residue in Raf-1 and the equivalent Ser-615 in B-Raf revealed that this residue regulates Raf binding to MEK. Further characterization revealed that in response to activation by epidermal growth factor, the Raf-S508N mutant protein displayed both reduced catalytic activity and aberrant activation kinetics: characteristics that may explain the C. elegans phenotype.

Impulse-response function of splanchnic circulation with model-independent constraints: theory and experimental validation

Modeling physiological processes using tracer kinetic methods requires knowledge of the time course of the tracer concentration in blood supplying the organ. For liver studies, however, inaccessibility of the portal vein makes direct measurement of the hepatic dual-input function impossible in humans. We want to develop a method to predict the portal venous time-activity curve from measurements of an arterial time-activity curve. An impulse-response function based on a continuous distribution of washout constants is developed and validated for the gut. Experiments with simultaneous blood sampling in aorta and portal vein were made in 13 anesthetized pigs following inhalation of intravascular [O-15] CO or injections of diffusible 3-O[ C-11] methylglucose (MG). The parameters of the impulse-response function have a physiological interpretation in terms of the distribution of washout constants and are mathematically equivalent to the mean transit time ( T) and standard deviation of transit times. The results include estimates of mean transit times from the aorta to the portal vein in pigs: (T) over bar = 0.35 +/- 0.05 min for CO and 1.7 +/- 0.1 min for MG. The prediction of the portal venous time-activity curve benefits from constraining the regression fits by parameters estimated independently. This is strong evidence for the physiological relevance of the impulse-response function, which includes asymptotically, and thereby justifies kinetically, a useful and simple power law. Similarity between our parameter estimates in pigs and parameter estimates in normal humans suggests that the proposed model can be adapted for use in humans.

Product form approximations for highly linear loss networks with trunk reservation

Admission controls, such as trunk reservation, are often used in loss networks to optimise their performance. Since the numerical evaluation of performance measures is complex, much attention has been given to finding approximation methods. The Erlang Fixed-Point (EFP) approximation, which is based on an independent blocking assumption, has been used for networks both with and without controls. Several more elaborate approximation methods which account for dependencies in blocking behaviour have been developed for the uncontrolled setting. This paper is an exploratory investigation of extensions and synthesis of these methods to systems with controls, in particular, trunk reservation. In order to isolate the dependency factor, we restrict our attention to a highly linear network. We will compare the performance of the resulting approximations against the benchmark of the EFP approximation extended to the trunk reservation setting. By doing this, we seek to gain insight into the critical factors in constructing an effective approximation. (C) 2003 Elsevier Ltd. All rights reserved.

Mean anisotropy of homogeneous Gaussian random fields and anisotropic norms of linear translation-invariant operators on multidimensional integer lattices

Sensitivity of output of a linear operator to its input can be quantified in various ways. In Control Theory, the input is usually interpreted as disturbance and the output is to be minimized in some sense. In stochastic worst-case design settings, the disturbance is considered random with imprecisely known probability distribution. The prior set of probability measures can be chosen so as to quantify how far the disturbance deviates from the white-noise hypothesis of Linear Quadratic Gaussian control. Such deviation can be measured by the minimal Kullback-Leibler informational divergence from the Gaussian distributions with zero mean and scalar covariance matrices. The resulting anisotropy functional is defined for finite power random vectors. Originally, anisotropy was introduced for directionally generic random vectors as the relative entropy of the normalized vector with respect to the uniform distribution on the unit sphere. The associated a-anisotropic norm of a matrix is then its maximum root mean square or average energy gain with respect to finite power or directionally generic inputs whose anisotropy is bounded above by a≥0. We give a systematic comparison of the anisotropy functionals and the associated norms. These are considered for unboundedly growing fragments of homogeneous Gaussian random fields on multidimensional integer lattice to yield mean anisotropy. Correspondingly, the anisotropic norms of finite matrices are extended to bounded linear translation invariant operators over such fields.

An augmented Lanczos algorithm for the efficient computation of a dot-product of a function of a large sparse symmetric matrix

The Lanczos algorithm is appreciated in many situations due to its speed. and economy of storage. However, the advantage that the Lanczos basis vectors need not be kept is lost when the algorithm is used to compute the action of a matrix function on a vector. Either the basis vectors need to be kept, or the Lanczos process needs to be applied twice. In this study we describe an augmented Lanczos algorithm to compute a dot product relative to a function of a large sparse symmetric matrix, without keeping the basis vectors.