A numerical study of large sparse matrix exponentials arising in Markov chains

Krylov subspace techniques have been shown to yield robust methods for the numerical computation of large sparse matrix exponentials and especially the transient solutions of Markov Chains. The attractiveness of these methods results from the fact that they allow us to compute the action of a matrix exponential operator on an operand vector without having to compute, explicitly, the matrix exponential in isolation. In this paper we compare a Krylov-based method with some of the current approaches used for computing transient solutions of Markov chains. After a brief synthesis of the features of the methods used, wide-ranging numerical comparisons are performed on a power challenge array supercomputer on three different models. (C) 1999 Elsevier Science B.V. All rights reserved.AMS Classification: 65F99; 65L05; 65U05.

Interaction of zinc(II) with the cyclic octapeptides, cyclo[Ile(Oxn)-D-Val(Thz)](2) and ascidiacyclamide, a cyclic peptide from Lissoclinum patella

The interactions between zinc salts and the naturally occurring cyclic octapeptide ascidiacyclamide in methanol, as well as a synthetic analogue cyclo[Ile(Oxn)-D-Val(Thz)](2), were monitored by H-1 NMR and CD spectroscopy. Three zinc complexes were identified, their relative amounts depending on the nature of the anion (perchlorate, triflate or chloride) and the presence or absence of base. Binding constants for two of the zinc species were calculated from CD or H-1 NMR spectra, [Zn(L - H)](+) (KZn(L-H) = [Zn(L - H)(+)]/[Zn2+][(L - H)(-)] = 10(7 +/- 2) M-1; 95% methanol/5% water, 298.0 K, NEt3/HClO4 buffer 0.04 M) and [ZnLCl](+) (K-ZnCIL = [ZnCIL+]/[Zn2+][Cl-][L] = 10(7.2) (+/-) (0.1) M-2; d(3)-methanol, 301 K).

Strategies of accommodation: Development of a coding system for conversational interaction

This study describes a coding system developed to operationalize the sociolinguistic strategies proposed by communication accommodation theory (CAT) in an academic context. Fifty interactions between two students (of Australian or Chinese ethnic background) or a student and faculty member were videotaped. A turn- and episode-based coding system was developed, focusing on verbal and nonverbal behavior. The development of this system is described in detail, before results are presented. Results indicated that status was the main influence on choice of strategies, particularly the extent and type of discourse management and interpersonal control. Participants' sew and ethnicity also played a role: Male participants made more use of interpretability (largely questions), whereas female participants used discourse management to develop a shared perspective. The results make clear that there is no automatic correspondence between behaviors and the strategies they constitute, and they point to the appropriateness of conceptualizing behavior and strategies separately in CAT.

Interaction between Helicoverpa armigera and Colletotrichum gloeosporioides on the tropical pasture legume Stylosanthes scabra

Glasshouse experiments determined effects of a moth, Helicoverpa armigera (Lepidoptera: Noctuidae), and the anthracnose pathogen, Colletotrichum gloeosporioides (Penz.) Penz. and Sacc., on each other when attacking the same host plant, Stylosanthes scabra (Vog.) (Leguminosae) cv. Fitzroy. The host was treated with both organisms in 2 ways of succession and at 2 different life stages each. Larvae of the moth preferred to feed on healthy plants rather than plants recently infected with C. gloeosporioides, and preferred such newly infected plants to severely diseased ones. Adult female moths laid more eggs on healthy and recently infected plants than on diseased plants, when given a choice of all 3 plant types. Severity of anthracnose disease was neither promoted nor retarded by damage to leaves caused by larvae of the moth.

Identification of regions within the third FnIII-like domain of the IL-5R alpha involved in IL-5 interaction

Previously, two binding sites for interleukin 5 (IL-5) were identified on the IL-5 receptor alpha chain (IL-5R alpha). They are located within the CD loop of the first fibronectin type III (FnIII)-like domain and the EF loop of the second FnIII-like domain. The first binding site was identified by exploiting the different abilities of human IL-5R alpha (hIL-5R alpha) and mouse IL-5R alpha (mIL-5R alpha) to bind hIL-5. Here we show that ovine IL-5 (oIL-5) has the ability to activate the hIL-5R alpha but not the mIL-5R alpha. By using chimeras of the mIL-5R alpha and hIL-5R alpha we demonstrate that residues within the first and third FnIII-like domains of mIL-5R alpha are responsible for this lack of activity. Furthermore, mutation of residues on hIL-5R alpha to mIL-5R alpha within the predicted DE and FG loop regions of the third FnIII domain reduces oIL-5 activity, These results show that regions of the third FnIII domain of IL-5R alpha are involved in binding, in addition to the regions in domains one and two of the IL-5R alpha that were identified in an earlier study. (C) 2000 Academic Press.

Conotoxin TVIIA, a novel peptide from the venom of Conus tulipa - 2. Three-dimensional solution structure

The three-dimensional solution structure of conotoxin TVIIA, a 30-residue polypeptide from the venom of the piscivorous cone snail Conus tulipa, has been determined using 2D H-1 NMR spectroscopy. TVIIA contains six cysteine residues which form a 'four-loop' structural framework common to many peptides from Conus venoms including the omega-, delta-, kappa-, and mu O-conotoxins. However, TVIIA does not belong to these well-characterized pharmacological classes of conotoxins, but displays high sequence identity with conotoxin GS, a muscle sodium channel blocker from Conus geographus. Structure calculations were based on 562 interproton distance restraints inferred from NOE data, together with 18 backbone and nine side-chain torsion angle restraints derived from spin-spin coupling constants. The final family of 20 structures had mean pairwise rms differences over residues 2-27 of 0.18 +/- 0.05 Angstrom for the backbone atoms and 1.39 +/- 0.33 Angstrom for all heavy atoms. The structure consists of a triple-stranded, antiparallel beta sheet with +2x, -1 topology (residues 7-9, 16-20 and 23-27) and several beta turns. The core of the molecule is formed by three disulfide bonds which form a cystine knot motif common to many toxic and inhibitory polypeptides. The global fold, molecular shape and distribution of amino-acid sidechains in TVIIA is similar to that previously reported for conotoxin GS, and comparison with other four-loop conotoxin structures provides further indication that TVIIA and GS represent a new and distinct subgroup of this structural family. The structure of TVIIA determined in this study provides the basis for determining a structure-activity relationship for these molecules and their interaction with target receptors.

The interaction between the equity premium and the risk-free rate

The equity premium arises from the interaction between the atemporal risk premium for equity, the risk-free rate of intertemporal substitution and the impact of risk on the precautionary motive for saving. Depending on parameter values, the equity premium may either be increased or reduced by the presence of undiversifiable background risk. (C) 2000 Elsevier Science S.A. All rights reserved.

Finite element modeling of fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins

In order to use the finite element method for solving fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins effectively and efficiently, we have presented, in this paper, the new concept and numerical algorithms to deal with the fundamental issues associated with the fluid-rock interaction problems. These fundamental issues are often overlooked by some purely numerical modelers. (1) Since the fluid-rock interaction problem involves heterogeneous chemical reactions between reactive aqueous chemical species in the pore-fluid and solid minerals in the rock masses, it is necessary to develop the new concept of the generalized concentration of a solid mineral, so that two types of reactive mass transport equations, namely, the conventional mass transport equation for the aqueous chemical species in the pore-fluid and the degenerated mass transport equation for the solid minerals in the rock mass, can be solved simultaneously in computation. (2) Since the reaction area between the pore-fluid and mineral surfaces is basically a function of the generalized concentration of the solid mineral, there is a definite need to appropriately consider the dependence of the dissolution rate of a dissolving mineral on its generalized concentration in the numerical analysis. (3) Considering the direct consequence of the porosity evolution with time in the transient analysis of fluid-rock interaction problems; we have proposed the term splitting algorithm and the concept of the equivalent source/sink terms in mass transport equations so that the problem of variable mesh Peclet number and Courant number has been successfully converted into the problem of constant mesh Peclet and Courant numbers. The numerical results from an application example have demonstrated the usefulness of the proposed concepts and the robustness of the proposed numerical algorithms in dealing with fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins. (C) 2001 Elsevier Science B.V. All rights reserved.

Transfer matrix eigenvalues of the anisotropic multiparametric U model

A multiparametric extension of the anisotropic U model is discussed which maintains integrability. The R-matrix solving the Yang-Baxter equation is obtained through a twisting construction applied to the underlying U-q(sl (2/1)) superalgebraic structure which introduces the additional free parameters that arise in the model. Three forms of Bethe ansatz solution for the transfer matrix eigenvalues are given which we show to be equivalent.

Measuring the interaction forces between protein inclusion bodies and an air bubble using an atomic force microscope

Interaction forces between protein inclusion bodies and an air bubble have been quantified using an atomic force microscope (AFM). The inclusion bodies were attached to the AFM tip by covalent bonds. Interaction forces measured in various buffer concentrations varied from 9.7 nN to 25.3 nN (+/- 4-11%) depending on pH. Hydrophobic forces provide a stronger contribution to overall interaction force than electrostatic double layer forces. It also appears that the ionic strength affects the interaction force in a complex way that cannot be directly predicted by DLVO theory. The effects of pH are significantly stronger for the inclusion body compared to the air bubble. This study provides fundamental information that will subsequently facilitate the rational design of flotation recovery system for inclusion bodies. It has also demonstrated the potential of AFM to facilitate the design of such processes from a practical viewpoint.

Interaction of RTD-1, a cyclic antimicrobial defensin from Rhesus macaque leukocytes, and its open chain analogue with membrane mimics

In contrast to other mammalian defensins, rhesus theta defensin-1 (RTD-1) is composed of just 18 amino acids with the backbone cyclized through peptide bonds. Antibacterial activities of both the native cyclic peptide and a linear form were examined, showing that the cyclic form was 3-fold more active than the open chain analogue, oRTD-1, although both peptides adopt very similar structures in water. It was suggested that the additional charges at the termini of oRTD-1 are the cause for its lower antimicrobial activity. Therefore, we studied the interaction of both peptides with membrane mimics composed of zwitterionic (PC) and negatively charged (PG) phospholipids, major lipid components of erythrocyte and bacterial cell membranes, respectively. Microcalorimetry showed that RTD-1 and oRTD-1 did not affect the phase behavior of PC liposomes, while in PG liposomes both peptides induced new phase transitions above the chain melting transition of the lipid. The shape and fraction differed between both peptides, depending also on their concentration, which will be discussed in terms of their antimicrobial activity.

Expression, purification, crystallization, and NMR studies of the helicase interaction domain of Escherichia coli DnaG primase

in Escherichia coli, the DnaG primase is the RNA polymerase that synthesizes RNA primers at replication forks. It is composed of three domains, a small N-terminal zinc-binding domain, a larger central domain responsible for RNA synthesis, and a C-terminal domain comprising residues 434-581 [DnaG(434-581)] that interact with the hexameric DnaB helicase. Presumably because of this interaction, it had not been possible previously to express the C-terminal domain in a stably transformed E coli strain. This problem was overcome by expression of DnaG(434-581) under control of tandem bacteriophage gimel-promoters, and the protein was purified in yields of 4-6 mg/L of culture and studied by NMR. A TOCSY spectrum of a 2 mM solution of the protein at pH 7.0, indicated that its structured core comprises residues 444-579. This was consistent with sequence conservation among most-closely related primases. Linewidths in a NOESY spectrum of a 0.5 mM sample in 10 mM phosphate, pH 6.05, 0.1 M NaCl, recorded at 36 degreesC, indicated the protein to be monomeric. Crystals of selenomethionine-substituted DnaG(434-581) obtained by the hanging-drop vapor-diffusion method were body-centered tetragonal, space group I4(1)22, with unit cell parameters a = b 142.2 Angstrom, c = 192.1 Angstrom, and diffracted beyond 2.7 Angstrom resolution with synchrotron radiation. (C) 2003 Elsevier Inc. All rights reserved.