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.

On derivatives and norms of generalized matrix functions and respective symmetric powers

In recent papers, the authors obtained formulas for directional derivatives of all orders, of the immanant and of the m-th xi-symmetric tensor power of an operator and a matrix, when xi is a character of the full symmetric group. The operator norm of these derivatives was also calculated. In this paper, similar results are established for generalized matrix functions and for every symmetric tensor power.

An Iterative Method for the Matrix Principal n-th Root

In this paper we give an iterative method to compute the principal n-th root and the principal inverse n-th root of a given matrix. As we shall show this method is locally convergent. This method is analyzed and its numerical stability is investigated.

Observability of nonlinear fractional dynamical systems

We study the observability of linear and nonlinear fractional differential systems of order 0 < α < 1 by using the Mittag-Leffler matrix function and the application of Banach’s contraction mapping theorem. Several examples illustrate the concepts.

Not the outside but the inside matters : functional analysis of Capricious and Tartan during Drosophila tracheal morphogenesis

The tubular structures, which transport essential gases, liquids, or cells from one site to another, are shared among various divergent organisms. These highly organized tubular networks include lung, kidney, vasculature and mammary gland in mammals as well as trachea and salivary gland in Drosophila melanogaster. Many questions regarding the tubular morphogenesis cannot be addressed sufficiently by investigating the mammalian organs because their structures are extremely complex and therefore, systematic analyses of genetic and cellular programs guiding the development is not possible. In contrast, the Drosophila tracheal development provides an excellent model system since many molecular markers and powerful tools for genetic manipulations are available. Two mechanisms were shown to be important for the outgrowth of tracheal cells: the FGF signaling pathway and the interaction between the tracheal cells and the surrounding mesodermal cells. The Drosophila FGF ligand encoded by branchless (bnl) is localized in groups of cells near tracheal metameres. The tracheal cells expressing the FGF receptor breathless (btl) respond to these sources of FGF ligand and extend towards them. However, this FGF signaling pathway is not sufficient for the formation of continuous dorsal trunk, the only muticellular tube in tracheal system. Recently, it was found out that single mesodermal cells called bridge-cells are essential for the formation of continuous dorsal trunk as they direct the outgrowth of dorsal trunk cells towards the correct targets. The results in this PhD thesis demonstrate that a cell adhesion molecule Capricious (Caps), which is specifically localized on the surface of bridge-cells, plays an essential role in guiding the outgrowing dorsal trunk cells towards their correct targets. When caps is lacking, some bridge-cells cannot stretch properly towards the adjacent posterior tracheal metameres and thus fail to interconnect the juxtaposing dorsal trunk cells. Consequently, discontinuous dorsal trunks containing interruptions at several positions are formed. On the other hand, when caps is ectopically expressed in the mesodermal cells through a twi-GAL4 driver, these mesodermal cells acquire a guidance function through ectopic caps and misguide the outgrowing dorsal trunk cells in abnormal directions. As a result, disconnected dorsal trunks are formed. These loss- and gain-of-function studies suggest that Caps presumably establishes the cell-to-cell contact between the bridge-cells and the tracheal cells and thereby mediates directly the guidance function of bridge-cells. The most similar protein known to Caps is another cell adhesion molecule called Tartan (Trn). Interestingly, trn is expressed in the mesodermal cells but not in the bridge-cells. When trn is lacking, the outgrowth of not only the dorsal trunks but also the lateral trunks are disrupted. However, in contrast to the ectopic expression of caps, the misexpression of trn does not affect tracheal development. Whereas Trn requires only its extracellular domain to mediate the matrix function, Caps requires both its extracellular and intracellular domains to function as a guidance molecule in the bridge-cells. These observations suggest that Trn functions differently from Caps during tracheal morphogenesis. Presumably, Trn mediates a matrix function of mesodermal cells, which support the tracheal cells to extend efficiently through the surrounding mesodermal tissue. In order to determine which domains dictate the functional specificity of Caps, two hybrid proteins CapsEdTrnId, which contains the Caps extracellular domain and the Trn intracellular domain, and TrnEdCapsId, which consists of the Trn extracellular domain and the Caps intracellular domain, were constructed. Gain of function and rescue experiments with these hybrid proteins suggest on one hand that the extracellular domains of Caps and Trn are functionally redundant and on the other hand that the intracellular domain dictates the functional specificity of Caps. In order to identify putative interactors of Caps, yeast two-hybrid screening was performed. An in vivo interaction assay in yeast suggests that Ras64B interacts specifically with the Caps intracellular domain. In addition, an in vitro binding assay reveals a direct interaction between an inactive form of Ras64B and the Caps intracellular domain. ras64B, which encodes a small GTPase, is expressed in the mesodermal cells concurrently as caps. Finally, a gain-of-function study with the constitutively active Ras64B suggests that Ras64B presumably functions downstream of Caps. All these results suggest consistently that the small GTPase Ras64B binds specifically to the Caps intracellular domain and may thereby mediate the guidance function of Caps.

Weighted geometric mean of large-scale matrices: numerical analysis and algorithms

Computing the weighted geometric mean of large sparse matrices is an operation that tends to become rapidly intractable, when the size of the matrices involved grows. However, if we are not interested in the computation of the matrix function itself, but just in that of its product times a vector, the problem turns simpler and there is a chance to solve it even when the matrix mean would actually be impossible to compute. Our interest is motivated by the fact that this calculation has some practical applications, related to the preconditioning of some operators arising in domain decomposition of elliptic problems. In this thesis, we explore how such a computation can be efficiently performed. First, we exploit the properties of the weighted geometric mean and find several equivalent ways to express it through real powers of a matrix. Hence, we focus our attention on matrix powers and examine how well-known techniques can be adapted to the solution of the problem at hand. In particular, we consider two broad families of approaches for the computation of f(A) v, namely quadrature formulae and Krylov subspace methods, and generalize them to the pencil case f(A\B) v. Finally, we provide an extensive experimental evaluation of the proposed algorithms and also try to assess how convergence speed and execution time are influenced by some characteristics of the input matrices. Our results suggest that a few elements have some bearing on the performance and that, although there is no best choice in general, knowing the conditioning and the sparsity of the arguments beforehand can considerably help in choosing the best strategy to tackle the problem.

The block numerical range of analytic operator functions

We introduce the block numerical range Wn(L) of an operator function L with respect to a decomposition H = H1⊕. . .⊕Hn of the underlying Hilbert space. Our main results include the spectral inclusion property and estimates of the norm of the resolvent for analytic L . They generalise, and improve, the corresponding results for the numerical range (which is the case n = 1) since the block numerical range is contained in, and may be much smaller than, the usual numerical range. We show that refinements of the decomposition entail inclusions between the corresponding block numerical ranges and that the block numerical range of the operator matrix function L contains those of its principal subminors. For the special case of operator polynomials, we investigate the boundedness of Wn(L) and we prove a Perron-Frobenius type result for the block numerical radius of monic operator polynomials with coefficients that are positive in Hilbert lattice sense.

Matrix Metalloproteinases And Left Ventricular Function And Structure In Spinal Cord Injured Subjects.

Subjects with spinal cord injury (SCI) exhibit impaired left ventricular (LV) diastolic function, which has been reported to be attenuated by regular physical activity. This study investigated the relationship between circulating matrix metalloproteinases (MMPs) and tissue inhibitors of MMPs (TIMPs) and echocardiographic parameters in SCI subjects and the role of physical activity in this regard. Forty-two men with SCI [19 sedentary (S-SCI) and 23 physically-active (PA-SCI)] were evaluated by clinical, anthropometric, laboratory, and echocardiographic analysis. Plasmatic pro-MMP-2, MMP-2, MMP-8, pro-MMP-9, MMP-9, TIMP-1 and TIMP-2 levels were determined by enzyme-linked immunosorbent assay and zymography. PA-SCI subjects presented lower pro-MMP-2 and pro-MMP-2/TIMP-2 levels and improved markers of LV diastolic function (lower E/Em and higher Em and E/A values) than S-SCI ones. Bivariate analysis showed that pro-MMP-2 correlated inversely with Em and directly with E/Em, while MMP-9 correlated directly with LV mass index and LV end-diastolic diameter in the whole sample. Following multiple regression analysis, pro-MMP-2, but not physical activity, remained associated with Em, while MMP-9 was associated with LV mass index in the whole sample. These findings suggest differing roles for MMPs in LV structure and function regulation and an interaction among pro-MMP-2, diastolic function and physical activity in SCI subjects.

Expression of smooth muscle and extracellular matrix proteins in relation to airway function in asthma

Background: Smooth muscle content is increased within the airway wall in patients with asthma and is likely to play a role in airway hyperresponsiveness. However, smooth muscle cells express several contractile and structural proteins, and each of these proteins may influence airway function distinctly. Objective: We examined the expression of contractile and structural proteins of smooth muscle cells, as well as extracellular matrix proteins, in bronchial biopsies of patients with asthma, and related these to lung function, airway hyperresponsiveness, and responses to deep inspiration. Methods: Thirteen patients with asthma (mild persistent, atopic, nonsmoking) participated in this cross-sectional study. FEV1 % predicted, PC20 methacholine, and resistance of the respiratory system by the forced oscillation technique during tidal breathing and deep breath were measured. Within 1 week, a bronchoscopy was performed to obtain 6 bronchial biopsies that were immunuhistochemically stained for alpha-SM-actin, desmin, myosin light chain kinase (MLCK), myosin, calponin, vimentin, elastin, type III collagen, and fibronectin. The level of expression was determined by automated densitometry. Results: PC20 methacholine was inversely related to the expression of alpha-smooth muscle actin (r = -0.62), desmin (r = -0.56), and elastin (r = -0.78). In addition, FEV1% predicted was positively related and deep inspiration-induced bronchodilation inversely related to desmin (r = -0.60), MLCK (r = -0.60), and calponin (r = -0.54) expression. Conclusion: Airway hyperresponsiveness, FEV1% predicted, and airway responses to deep inspiration are associated with selective expression of airway smooth muscle proteins and components of the extracellular matrix.

Extracellular matrix and platelet function in patients with musculocontractural Ehlers-Danlos syndrome caused by mutations in the CHST14 gene.

We report on a consanguineous, Afghani family with two sisters affected with characteristic facial features, multiple contractures, progressive joint and skin laxity, hemorrhagic diathesis following minor trauma and multisystem fragility-related manifestations suggestive of a diagnosis of musculocontractural Ehlers-Danlos syndrome (EDS). This novel form of connective tissue disorder was recently reported in patients of Japanese, Turkish, and Indian descent who were formerly classified as having EDS type VIB and has now been recognized to be a part of spectrum including patients previously classified as having adducted thumb-clubfoot syndrome. We identified a previously unreported mutation in the CHST14 gene, which codes for the enzyme dermatan 4-O-sulfotransferase. We discuss the prenatal presentation, detailed clinical manifestations, and neurological findings in two sisters with this newly described musculocontractural EDS-CHST14 type. We demonstrate that fibroblasts from one of our patients produce more chondroitin sulfate than normal and show lower than normal deposition of collagens I and II and fibrillin 1-containing microfibrills. These findings suggest that the imbalance in the glycosaminoglycan content in developing tissues might interfere with normal deposition of other extracellular matrix components and ultimately contribute to the development of the phenotype observed in these patients. Furthermore, we ruled out the contribution of intrinsic platelet factors to the bleeding diathesis observed in some affected individuals. © 2012 Wiley Periodicals, Inc.

Numerical computation of an analytic singular value decomposition of a matrix valued function

This paper extends the singular value decomposition to a path of matricesE(t). An analytic singular value decomposition of a path of matricesE(t) is an analytic path of factorizationsE(t)=X(t)S(t)Y(t) T whereX(t) andY(t) are orthogonal andS(t) is diagonal. To maintain differentiability the diagonal entries ofS(t) are allowed to be either positive or negative and to appear in any order. This paper investigates existence and uniqueness of analytic SVD's and develops an algorithm for computing them. We show that a real analytic pathE(t) always admits a real analytic SVD, a full-rank, smooth pathE(t) with distinct singular values admits a smooth SVD. We derive a differential equation for the left factor, develop Euler-like and extrapolated Euler-like numerical methods for approximating an analytic SVD and prove that the Euler-like method converges.

Pole placement method using the system matrix transfer function and sparsity

This paper presents a pole placement method using both the augmented Jacobian and the corresponding system transfer function matrices. From the manipulation of these matrices a straightforward approach results to get the coefficients of a non-linear system, whose solution gives the parameters of the stabilizers that can provide a pre-specified minimum damping to the system. (C) 2001 Elsevier B.V. Ltd. All rights reserved.

Extracellular matrix in airway smooth muscle is associated with dynamics of airway function in asthma

Background: Altered deposition of extracellular matrix (ECM) in the airway smooth muscle (ASM) layer as observed in asthma may influence ASM mechanical properties. We hypothesized that ECM in ASM is associated with airway function in asthma. First, we investigated the difference in ECM expression in ASM between asthma and controls. Second, we examined whether ECM expression is associated with bronchoconstriction and bronchodilation in vivo. Methods: Our cross-sectional study comprised 19 atopic mild asthma patients, 15 atopic and 12 nonatopic healthy subjects. Spirometry, methacholine responsiveness, deep-breath-induced bronchodilation (Delta R-rs) and bronchoscopy with endobronchial biopsies were performed. Positive staining of elastin, collagen I, III and IV, decorin, versican, fibronectin, laminin and tenascin in ASM was quantified as fractional area and mean density. Data were analysed using Pearson's or Spearman's correlation coefficient. Results: Extracellular matrix expression in ASM was not different between asthma and controls. In asthmatics, fractional area and mean density of collagen I and III were correlated with methacholine dose-response slope and DRrs, respectively (r = 0.71, P < 0.01; r = 0.60, P = 0.02). Furthermore, ASM collagen III and laminin in asthma were correlated with FEV1 reversibility (r = -0.65, P = 0.01; r = -0.54, P = 0.04). Conclusion: In asthma, ECM in ASM is related to the dynamics of airway function in the absence of differences in ECM expression between asthma and controls. This indicates that the ASM layer in its full composition is a major structural component in determining variable airways obstruction in asthma.

Molecular characterization of the immune modulatory function of matrix metalloproteinase-7, a factor in the tumour micromilieu

Matrix metalloproteinases are the components of the tumour microenvironment which play a crucial role in tumour progression. Matrix metalloproteinase-7 (MMP-7) is expressed in a variety of tumours and the expression is associated with an aggressive malignant phenotype and poor prognosis. A role for MMP-7 in the immune escape of tumours has been postulated, but the mechanisms are not clearly understood. The present study was focused on identifying physiological inactivators of MMP-7 and also to unravel the mechanisms involved in MMP-7 mediated immune escape. This study shows that human leukocyte elastase (HLE), secreted by polymorphonuclear leukocytes cleaves MMP-7 in the catalytic domain as revealed by N-terminal sequencing. Further analysis demonstrates that the activity of MMP-7 was drastically decreased after HLE treatment in a time and dose dependent manner. MMP-7 induces apoptosis resistance in tumour cells by cleaving CD95 and CD95L. The effect of HLE on MMP-7 mediated apoptosis resistance was analysed. In vitro stimulation of apoptosis by anti-Apo-1 (anti-CD95 antibody) and the chemotherapeutic drug doxorubicin is reduced by MMP-7. Also tumour specific cytotoxic T cells do not effectively kill tumour cells in the presence of MMP-7. This study revealed that HLE abrogates the negative effect of MMP-7 on apoptosis induced by CD95 stimulation, doxorubicin or cytotoxic T cells and restores apoptosis sensitivity of tumour cells. To gain insight into the possible immune modulatory functions of MMP-7, experiments were performed to identify new immune relevant substrates. The human T cell line, Jurkat, was selected for these studies. Hsc70 which is involved in uncoating of clathrin vesicles was found in the supernatants of the MMP-7 treated cells indicating a modulatory role of MMP-7 on endocytosis. Further studies demonstrated that MMP-7 leads to decreased clathrin staining in HEK293, HepG2, Jurkat, CD4+ T cells and dendritic cells. Results also show MMP-7 treatment increased surface expression of cytotoxic T lymphocyte associated protein-4 (CTLA-4) which accumulated due to inhibition of the clathrin mediated internalization in CD4+CD25+ cells.

Should 3K zoom function be used for detection of pneumothorax in cesium iodide/amorphous silicon flat-panel detector radiographs presented on 1K-matrix soft copies?

The purpose of the study was to evaluate observer performance in the detection of pneumothorax with cesium iodide and amorphous silicon flat-panel detector radiography (CsI/a-Si FDR) presented as 1K and 3K soft-copy images. Forty patients with and 40 patients without pneumothorax diagnosed on previous and subsequent digital storage phosphor radiography (SPR, gold standard) had follow-up chest radiographs with CsI/a-Si FDR. Four observers confirmed or excluded the diagnosis of pneumothorax according to a five-point scale first on the 1K soft-copy image and then with help of 3K zoom function (1K monitor). Receiver operating characteristic (ROC) analysis was performed for each modality (1K and 3K). The area under the curve (AUC) values for each observer were 0.7815, 0.7779, 0.7946 and 0.7066 with 1K-matrix soft copies and 0.8123, 0.7997, 0.8078 and 0.7522 with 3K zoom. Overall detection of pneumothorax was better with 3K zoom. Differences between the two display methods were not statistically significant in 3 of 4 observers (p-values between 0.13 and 0.44; observer 4: p = 0.02). The detection of pneumothorax with 3K zoom is better than with 1K soft copy but not at a statistically significant level. Differences between both display methods may be subtle. Still, our results indicate that 3K zoom should be employed in clinical practice.