Getting to the inside of cells using metabolic control analysis

Metabolic control analysis can relate control properties of an intact system to kinetic properties (elasticity coefficients) of the enzymes within that system. The method formulating the former as matrix inverse of the latter is elaborated here for the general case and founded in standard metabolic control theory. Then a method is developed that accomplishes the reverse: it is shown that a matrix containing all elasticity coefficients and information concerning the pathway structure equals the inverse of a matrix containing flux and concentration control coefficients. As a consequence, by measuring the control properties of an intact system, one is able to deduce its in situ pathway structure and enzyme kinetic properties: This solves the ever-present question of whether the kinetic properties of enzymes in their isolated state differ from those under the conditions prevailing in the cell.     

The matrix method of metabolic control analysis: its validity for complex pathway structures

The sensitivities of the variables of a metabolic system (such as fluxes and concentrations) to variations in enzyme concentration are expressed in metabolic control analysis as control coefficients. The matrix method is a system of writing matrix equations that generate expressions for the control coefficients in terms of the characteristics of the components (principally the enzymes). Previously, the matrix method has been considered in terms of simple pathway structures; here we justify its applicability to complex pathways, such as those with multiple branches. It is shown that this requires modification of the branch point relationship to take account of changes of flux along the limbs of the branch and of stoichiometric factors. The method of deriving the flux control coefficients with respect to different fluxes in the system is extended to cope with these circumstances.     

High-Performance Tunable Multichannel Absorbers Coupled with Graphene-Based Grating and Dual-Tamm Plasmonic Structures

Plasmonics - We present a hybrid Tamm system targeting the tunable multichannel absorber. The proposed optical absorber is analyzed and investigated by using the transfer matrix method (TMM). The...     

Asymptotic stability of delay-difference system of hopfield neural networks via matrix inequalities and application

In this paper, we derive a sufficient condition for asymptotic stability of the zero solution of delay-difference system of Hopfield neural networks in terms of certain matrix inequalities by using a discrete version of the Lyapunov second method. The result is applied to obtain new asymptotic stability condition for some class of delay-difference system such as delay-difference system of Hopfield neural networks with multiple delays in terms of certain matrix inequalities. Our results can be well suited for computational purposes.     

High-sensitivity protein detection by a new "contact-copy" method using a protein A-neomycin phosphotransferase II fusion protein

A new system for high-sensitivity protein detection by an immunoenzymatic "contact-copy" procedure is described. It is based on two components: (i) a microbiologically produced bifunctional fusion protein of protein A and neomycin phosphotransferase II (protein A-NPT II) in which the protein A moiety acts as a second immunological reagent while NPT II catalyzes the detection reaction and (ii) a novel kanamycin-loaded substrate matrix (kanamycin-cyanuric chloride-activated and sulfanilic acid-derivatized paper) brought into direct contact with a protein-carrying matrix after blot or dot application and initial immunoreaction--the NPT II enzyme reaction with [gamma-32P]ATP as cosubstrate leads to phosphorylation of the substrate kanamycin on the substrate matrix, which is used for further analysis. The contact-copy method has at least the same detection sensitivity as procedures employing 125I-protein A, but allows extremely short exposure times and avoids probe prelabeling. Twenty-five picograms of specific protein blotted from sodium dodecyl sulfate-polyacrylamide gels onto nitrocellulose is detected after 15 min of autoradiography. The limit of detection in dot tests was found to be 10 pg per dot (3 mm2). The method is suitable for quantitative determination of antigens in the range down to 100 pg. Several contact copies of the same original protein-carrying matrix can be produced and used for detection or quantitative analysis without destroying the original matrix.     

Toward a solution of allocation in life cycle inventories: the use of least-squares techniques

Purpose  

The matrix method for the solution of the so-called inventory problem in LCA generally determines the inventory vector related to a specific system of processes by solving a system of linear equations. The paper proposes a new approach to deal with systems characterized by a rectangular (and thus non-invertible) coefficients matrix. The approach, based on the application of regression techniques, allows solving the system without using computational expedients such as the allocation procedure.     

A method for testing the stability of a steady-state system during the calculation of a response to large changes in regulator concentrations

A modification of the ‘finite decomposition’ method (Crabtree and Newsholme (1985) Curr. Top. Cell. Regul. 25, 21–76) for calculating physiological responses from sensitivities is described, to enable the system to be tested for stability at each step of the procedure. Instability is indicated by a change of sign of the determinant of the square matrix (N) in the governing equation for the system. The method cannot be used to predict responses beyond any step at which instability occurs.     

Using priced timed automaton to analyse the energy consumption in cloud computing environment

According to the fact that cloud servers have different energy consumption on different running states, as well as the energy waste problem caused by the mismatching between cloud servers and cloud tasks, we carry out researches on the energy optimal method achieved by a priced timed automaton for the cloud computing center in this paper. The priced timed automaton is used to model the running behaviors of the cloud computing system. After introducing the matching matrix of cloud tasks and cloud resources as well as the power matrix of the running states of cloud servers, we design a generation algorithm for the cloud system automaton based on the generation rules and reduction rules given ahead. Then, we propose another algorithm to settle the minimum path energy consumption problem in the cloud system automaton, therefore obtaining an energy optimal solution and an energy optimal value for the cloud system. A case study and repeated experimental analyses manifest that our method is effective and feasible.     

Metabolic network simulation using logical loop algorithm and Jacobian matrix

A novel method to accomplish efficient numerical simulation of metabolic networks for flux analysis was developed. The only inputs required are the set of stoichiometric balances and the atom mapping matrices of all components of the reaction network. The latter are used to automatically calculate isotopomer mapping matrices. Using the symbolic toolbox of MATLAB the analytical solution of the stoichiometric balance equation system, isotopomer balances and the analytical Jacobian matrix of the total set of stoichiometric and isotopomer balances are created automatically. The number of variables in the isotopomer distribution equation system is significantly reduced applying modified isotopomer mapping matrices. These allow lumping of several consecutive isotopomer reactions into a single one. The solution of the complete system of equations is improved by implementing an iterative logical loop algorithm and using the analytical Jacobian matrix. This new method provided quick and robust convergence to the root of such equation systems in all cases tested. The method was applied to a network of lysine producing Corynebacterium glutamicum. The resulting equation system with the dimension of 546 x 546 was directly derived from 12 isotopomer balance equations. The results obtained yielded identical labeling patterns for metabolites as compared to the relaxation method.     

Development of vacuum-ultraviolet circular dichroism measurement system using a polarizing undulator

We have developed an improved circular dichroism (CD) and linear dichroism (LD) simultaneous measurement system for the vacuum ultraviolet (VUV) region by polarization modulation techniques using a four-period Onuki-type crossed undulator as a polarized light source. The system has been constructed at the VUV beamline BL-5 in the electron storage ring TERAS, at AIST. Our improvements, in particular the adoption of an optical chopper as the detection method of incident light, have resulted in a flat baseline and a consequent simplification of the Mueller matrix calculation for our optical system. Based on the Mueller matrix calculation, we have successfully measured real VUV-CD and LD spectra of leucine films for wavelengths down to 160 nm with absolute optical constants. The obtained spectra show good consistency with spectra measured by conventional methods.     

Osteocyte calcium signaling response to bone matrix deformation

Osteocytes embedded in calcified bone matrix have been widely believed to play important roles in mechanosensing to achieve adaptive bone remodeling in a changing mechanical environment. In vitro studies have clarified several types of mechanical stimuli such as hydrostatic pressure, fluid shear stress, and direct deformation influence osteocyte functions. However, osteocyte response to mechanical stimuli in the bone matrix has not been clearly understood. In this study, we observed the osteocyte calcium signaling response to the quantitatively applied deformation in the bone matrix. A novel experimental system was developed to apply deformation to cultured bone tissue with osteocytes on a microscope stage. As a mechanical stimulus to the osteocytes in bone matrix, in-plane shear deformation was applied using a pair of glass microneedles to bone fragments, obtained from 13-day-old embryonic chick calvariae. Deformation of bone matrix and cells was quantitatively evaluated using an image correlation method by applying for differential interference contrast images of the matrix and fluorescent images of immunolabeled osteocytes, together with imaging of the cellular calcium transient using a ratiometric method. As a result, it was confirmed that the newly developed system enables us to apply deformation to bone matrix and osteocytes successfully under the microscope without significant focal plane shift or deviation from the observation view field. The system could be a basis for further development to investigate the mechanosensing mechanism of osteocytes in bone matrix through examination of various types of rapid biochemical signaling responses and intercellular communication induced by matrix deformation.     

Control involving metabolism and gene expression

Control of DNA supercoiling by the free-energy of hydrolysis of ATP that involves gene expression is analyzed in terms of three levels of unconnected metabolic pathways. These are synthesis and breakdown of topoisomerase mRNAs, synthesis and breakdown of topoisomerase proteins and supercoiling and relaxation of DNA. The so-called square-matrix method previously developed for the control of metabolic pathways, is extended to deal with this hierarchical control system. It turns out that also in this case, the matrix of control coefficients is equal to the inverse of the so-called elasticity matrix, which contains all relevant elasticity coefficients as well as information about the structure and connectedness of the pathways involved. For a simpler case of a hierarchy of two systems, we demonstrate that the explicit matrix inversion method may be replaced by an implicit method in which the regulatory effects that run through the other level are described by an additional elasticy coefficient which may then be treated as if local.     

Molecular confinement influences protein structure and enhances thermal protein stability

The sol-gel method of encapsulating proteins in a silica matrix was investigated as a potential experimental system for testing the effects of molecular confinement on the structure and stability of proteins. We demonstrate that silica entrapment (1) is fully compatible with structure analysis by circular dichroism, (2) allows conformational studies in contact with solvents that would otherwise promote aggregation in solution, and (3) generally enhances thermal protein stability. Lysozyme, alpha-lactalbumin, and metmyoglobin retained native-like solution structures following sol-gel encapsulation, but apomyoglobin was found to be largely unfolded within the silica matrix under control buffer conditions. The secondary structure of encapsulated apomyoglobin was unaltered by changes in pH and ionic strength of KCl. Intriguingly, the addition of other neutral salts resulted in an increase in the alpha-helical content of encapsulated apomyoglobin in accordance with the Hofmeister ion series. We hypothesize that protein conformation is influenced directly by the properties of confined water in the pores of the silica. Further work is needed to differentiate the steric effects of the silica matrix from the solvent effects of confined water on protein structure and to determine the extent to which this experimental system mimics the effects of crowding and confinement on the function of macromolecules in vivo.     

An efficient method of computing the numerator relationship matrix and its inverse matrix with inbreeding for large sets of animals

Summary The numerator relationship matrix describes the genetic relationships between individuals of a population. Its inverse is used for the prediction of breeding values, as outlined by Henderson (1975a).For large populations, the recursive method commonly used is difficult to apply because of the size of the relationship matrix. Recently Henderson (1975b) derived a method which allows computing the inverse of the numerator relationship matrix itself for a large number of animals, provided the population is non-inbred. The method presented here is an extension of Henderson's method to allow for inbreeding with large number of animals. It takes inbreeding into account and computes the numerator relationship matrix as well as its inverse. The method is particularly efficient in computer storage in that it allows handling of sets of animals larger than 5000 animals, and is almost as fast as the recursive method.     

Inclusion of genetically identical animals to a numerator relationship matrix and modification of its inverse

In the field of animal breeding, estimation of genetic parameters and prediction of breeding values are routinely conducted by analyzing quantitative traits. Using an animal model and including the direct inverse of a numerator relationship matrix (NRM) into a mixed model has made these analyses possible. However, a method including a genetically identical animal (GIA) in NRM if genetic relationships between pairs of GIAs are not perfect, is still lacking. Here, we describe a method to incorporate GIAs into NRM using a K matrix in which diagonal elements are set to 1.0, off-diagonal elements between pairs of GIAs to (1-x) and the other elements to 0, where x is a constant less than 0.05. The inverse of the K matrix is then calculated directly by a simple formula. Thus, the inverse of the NRM is calculated by the products of the lower triangular matrix that identifies the parents of each individual, its transpose matrix, the inverse of the K matrix and the inverse of diagonal matrix D, in which the diagonal elements comprise a number of known parents and their inbreeding coefficients. The computing method is adaptable to the analysis of a data set including pairs of GIAs with imperfect relationships.     

Enhancing interacting residue prediction with integrated contact matrix prediction in protein-protein interaction

Identifying the residues in a protein that are involved in protein-protein interaction and identifying the contact matrix for a pair of interacting proteins are two computational tasks at different levels of an in-depth analysis of protein-protein interaction. Various methods for solving these two problems have been reported in the literature. However, the interacting residue prediction and contact matrix prediction were handled by and large independently in those existing methods, though intuitively good prediction of interacting residues will help with predicting the contact matrix. In this work, we developed a novel protein interacting residue prediction system, contact matrix-interaction profile hidden Markov model (CM-ipHMM), with the integration of contact matrix prediction and the ipHMM interaction residue prediction. We propose to leverage what is learned from the contact matrix prediction and utilize the predicted contact matrix as “feedback” to enhance the interaction residue prediction. The CM-ipHMM model showed significant improvement over the previous method that uses the ipHMM for predicting interaction residues only. It indicates that the downstream contact matrix prediction could help the interaction site prediction.     

A Class-Information-Based Penalized Matrix Decomposition for Identifying Plants Core Genes Responding to Abiotic Stresses

In terms of making genes expression data more interpretable and comprehensible, there exists a significant superiority on sparse methods. Many sparse methods, such as penalized matrix decomposition (PMD) and sparse principal component analysis (SPCA), have been applied to extract plants core genes. Supervised algorithms, especially the support vector machine-recursive feature elimination (SVM-RFE) method, always have good performance in gene selection. In this paper, we draw into class information via the total scatter matrix and put forward a class-information-based penalized matrix decomposition (CIPMD) method to improve the gene identification performance of PMD-based method. Firstly, the total scatter matrix is obtained based on different samples of the gene expression data. Secondly, a new data matrix is constructed by decomposing the total scatter matrix. Thirdly, the new data matrix is decomposed by PMD to obtain the sparse eigensamples. Finally, the core genes are identified according to the nonzero entries in eigensamples. The results on simulation data show that CIPMD method can reach higher identification accuracies than the conventional gene identification methods. Moreover, the results on real gene expression data demonstrate that CIPMD method can identify more core genes closely related to the abiotic stresses than the other methods.     

Confocal assay for invasion: use of propidium iodide fluorescence and laser reflectance to quantify the rate of migration of cells through a matrix

BACKGROUND: Most assays used to measure invasion are based on manual counting of the number of cells that have migrated completely through commercial coated filters. We describe here a confocal fluorescence-imaging method that can assess the relative rates of invasion of cells into a matrix. METHODS: After being seeded on the matrix and a period of incubation, the cells are fixed and treated with RNase. Propidium iodide is then added to stain the double-stranded DNA. A confocal microscope system is used to obtain high-resolution images of the red propidium iodide fluorescence and laser reflectance from optical sections at increasing depths in the matrix. The section with high laser reflectance marks the top of the matrix. RESULTS: Data were calculated as the total area of red fluorescence above background in each section and were plotted as a percentage of the summed fluorescent areas in all sections. CONCLUSIONS: Because the distance into the matrix of the nuclei can be calculated by measuring from the reflective upper surface of the matrix, the method is useful for assessing the rate of cell migration and for comparing the ability of different cells to invade through different matrices under varying conditions.