Finite volume evolution Galerkin method for hyperbolic conservation laws with spatially varying flux functions

We present a generalization of the finite volume evolution Galerkin scheme [M. Lukacova-Medvid'ova,J. Saibertov'a, G. Warnecke, Finite volume evolution Galerkin methods for nonlinear hyperbolic systems, J. Comp. Phys. (2002) 183 533-562; M. Luacova-Medvid'ova, K.W. Morton, G. Warnecke, Finite volume evolution Galerkin (FVEG) methods for hyperbolic problems, SIAM J. Sci. Comput. (2004) 26 1-30] for hyperbolic systems with spatially varying flux functions. Our goal is to develop a genuinely multi-dimensional numerical scheme for wave propagation problems in a heterogeneous media. We illustrate our methodology for acoustic waves in a heterogeneous medium but the results can be generalized to more complex systems. The finite volume evolution Galerkin (FVEG) method is a predictor-corrector method combining the finite volume corrector step with the evolutionary predictor step. In order to evolve fluxes along the cell interfaces we use multi-dimensional approximate evolution operator. The latter is constructed using the theory of bicharacteristics under the assumption of spatially dependent wave speeds. To approximate heterogeneous medium a staggered grid approach is used. Several numerical experiments for wave propagation with continuous as well as discontinuous wave speeds confirm the robustness and reliability of the new FVEG scheme.

targetTB: A target identification pipeline for Mycobacterium tuberculosis through an interactome, reactome and genome-scale structural analysis

Background: Tuberculosis still remains one of the largest killer infectious diseases, warranting the identification of newer targets and drugs. Identification and validation of appropriate targets for designing drugs are critical steps in drug discovery, which are at present major bottle-necks. A majority of drugs in current clinical use for many diseases have been designed without the knowledge of the targets, perhaps because standard methodologies to identify such targets in a high-throughput fashion do not really exist. With different kinds of 'omics' data that are now available, computational approaches can be powerful means of obtaining short-lists of possible targets for further experimental validation. Results: We report a comprehensive in silico target identification pipeline, targetTB, for Mycobacterium tuberculosis. The pipeline incorporates a network analysis of the protein-protein interactome, a flux balance analysis of the reactome, experimentally derived phenotype essentiality data, sequence analyses and a structural assessment of targetability, using novel algorithms recently developed by us. Using flux balance analysis and network analysis, proteins critical for survival of M. tuberculosis are first identified, followed by comparative genomics with the host, finally incorporating a novel structural analysis of the binding sites to assess the feasibility of a protein as a target. Further analyses include correlation with expression data and non-similarity to gut flora proteins as well as 'anti-targets' in the host, leading to the identification of 451 high-confidence targets. Through phylogenetic profiling against 228 pathogen genomes, shortlisted targets have been further explored to identify broad-spectrum antibiotic targets, while also identifying those specific to tuberculosis. Targets that address mycobacterial persistence and drug resistance mechanisms are also analysed. Conclusion: The pipeline developed provides rational schema for drug target identification that are likely to have high rates of success, which is expected to save enormous amounts of money, resources and time in the drug discovery process. A thorough comparison with previously suggested targets in the literature demonstrates the usefulness of the integrated approach used in our study, highlighting the importance of systems-level analyses in particular. The method has the potential to be used as a general strategy for target identification and validation and hence significantly impact most drug discovery programmes.

Scheme of squares: Part I. Systems formalism for potentiostatic studies

The so-called “Scheme of Squares”, displaying an interconnectivity of heterogeneous electron transfer and homogeneous (e.g., proton transfer) reactions, is analysed. Explicit expressions for the various partial currents under potentiostatic conditions are given. The formalism is applicable to several electrode geometries and models (e.g., semi-infinite linear diffusion, rotating disk electrodes, spherical or cylindrical systems) and the analysis is exact. The steady-state (t→∞) expressions for the current are directly given in terms of constant matrices whereas the transients are obtained as Laplace transforms that need to be inverted by approximation of numerical methods. The methodology employs a systems approach which replaces a system of partial differential equations (governing the concentrations of the several electroactive species) by an equivalent set of difference equations obeyed by the various partial currents.

The credit risk(+) model with general sector correlations

We consider an enhancement of the credit risk+ model to incorporate correlations between sectors. We model the sector default rates as linear combinations of a common set of independent variables that represent macro-economic variables or risk factors. We also derive the formula for exact VaR contributions at the obligor level.

Finite volume TVD formulation of lattice Boltzmann simulation on unstructured mesh

A numerical scheme is presented for accurate simulation of fluid flow using the lattice Boltzmann equation (LBE) on unstructured mesh. A finite volume approach is adopted to discretize the LBE on a cell-centered, arbitrary shaped, triangular tessellation. The formulation includes a formal, second order discretization using a Total Variation Diminishing (TVD) scheme for the terms representing advection of the distribution function in physical space, due to microscopic particle motion. The advantage of the LBE approach is exploited by implementing the scheme in a new computer code to run on a parallel computing system. Performance of the new formulation is systematically investigated by simulating four benchmark flows of increasing complexity, namely (1) flow in a plane channel, (2) unsteady Couette flow, (3) flow caused by a moving lid over a 2D square cavity and (4) flow over a circular cylinder. For each of these flows, the present scheme is validated with the results from Navier-Stokes computations as well as lattice Boltzmann simulations on regular mesh. It is shown that the scheme is robust and accurate for the different test problems studied.

Thermodynamic properties of Pb2PtO4 and PbPt2O4 and phase equilibria in the system Pb–Pt–O

The compounds Pb2PtO4 and PbPt2O4 were synthesized from an intimate mixture of yellow PbO and Pt metal powders by heating under pure oxygen gas at 973 K for periods up to 600 ks with intermediate grinding and recompacting. Both compounds were found to decompose on heating in pure oxygen to PbO and Pt, apparently in conflict with the requirements for equilibrium phase relations in the ternary system Pb–Pt–O. The oxygen chemical potential corresponding to the three-phase mixtures, Pb2PtO4 + PbO + Pt and PbPt2O4 + PbO + Pt were measured as a function of temperature using solid-state electrochemical cells incorporating yttria-stabilized zirconia as the solid electrolyte and pure oxygen gas at 0.1 MPa pressure as the reference electrode. The standard Gibbs free energies of formation of the ternary oxides were derived from the measurements. Analysis of the results indicated that the equilibrium involving three condensed phases Pb2PtO4 + PbO + Pt is metastable. Under equilibrium conditions, Pb2PtO4 should have decomposed to a mixture of PbO and PbPt2O4. Measurement of the oxygen potential corresponding to this equilibrium decomposition as a function of temperature indicated that decomposition temperature in pure oxygen is 1014(±2) K. This was further confirmed by direct determination of phase relations in the ternary Pb–Pt–O by equilibrating several compositions at 1023 K for periods up to 850 ks and phase identification of quenched samples using X-ray diffraction (XRD), scanning electron microscopy (SEM) and energy dispersive X-ray spectroscopy (EDS). Only one ternary oxide PbPt2O4 was stable at 1023 K under equilibrium conditions. Alloys and intermetallic compounds along the Pb–Pt binary were in equilibrium with PbO.

ESR evidence for 2 coexisting liquid phases in deeply supercooled bulk water

Using electron spin resonance spectroscopy (ESR), we measure the rotational mobility of probe molecules highly diluted in deeply supercooled bulk water and negligibly constrained by the possible ice fraction. The mobility increases above the putative glass transition temperature of water, T-g = 136 K, and smoothly connects to the thermodynamically stable region by traversing the so called "no man's land" (the range 150-235 K), where it is believed that the homogeneous nucleation of ice suppresses the liquid water. Two coexisting fractions of the probe molecules are evidenced. The 2 fractions exhibit different mobility and fragility; the slower one is thermally activated (low fragility) and is larger at low temperatures below a fragile-to-strong dynamic cross-over at approximate to 225 K. The reorientation of the probe molecules decouples from the viscosity below approximate to 225 K. The translational diffusion of water exhibits a corresponding decoupling at the same temperature [Chen S-H, et al. (2006) The violation of the Stokes-Einstein relation in supercooled water. Proc Natl Acad Sci USA 103:12974-12978]. The present findings are consistent with key issues concerning both the statics and the dynamics of supercooled water, namely the large structural fluctuations [Poole PH, Sciortino F, Essmann U, Stanley HE (1992) Phase behavior of metastable water. Nature 360: 324-328] and the fragile-to-strong dynamic cross-over at approximate to 228 K [Ito K, Moynihan CT, Angell CA (1999) Thermodynamic determination of fragility in liquids and a fragile-tostrong liquid transition in water. Nature 398: 492-494].

Static RAM cell for ternary logic

A novel CMOS static RAM cell for ternary logic systems is described. This cell is based on the lambda diode. The operation of the cell has been simulated using the SPICE 2G program. The results of the simulation are given.

Transcriptional events during early phase of germination of rice embryos

GERMINATION transfers a metabolically inert embryo into an active state of growth and development. The presence of conserved mRNAs has been demonstrated in different species of eggs and seeds1–4. In rice embryos, germination was shown to be independent of the synthesis of RNA up to 18–24 h after the start of imbibition5, although RNA synthesis was detected as early as 9 h after the start of imbibition. In this report, the sequence of the transcriptional events taking place during the early phase of the germination of rice embryos are presented.

A Finite Variable Difference Relaxation Scheme for hyperbolic-parabolic equations

Using the framework of a new relaxation system, which converts a nonlinear viscous conservation law into a system of linear convection-diffusion equations with nonlinear source terms, a finite variable difference method is developed for nonlinear hyperbolic-parabolic equations. The basic idea is to formulate a finite volume method with an optimum spatial difference, using the Locally Exact Numerical Scheme (LENS), leading to a Finite Variable Difference Method as introduced by Sakai [Katsuhiro Sakai, A new finite variable difference method with application to locally exact numerical scheme, journal of Computational Physics, 124 (1996) pp. 301-308.], for the linear convection-diffusion equations obtained by using a relaxation system. Source terms are treated with the well-balanced scheme of Jin [Shi Jin, A steady-state capturing method for hyperbolic systems with geometrical source terms, Mathematical Modeling Numerical Analysis, 35 (4) (2001) pp. 631-645]. Bench-mark test problems for scalar and vector conservation laws in one and two dimensions are solved using this new algorithm and the results demonstrate the efficiency of the scheme in capturing the flow features accurately.

Graphical simulation of Petri Nets

An interactive graphics package for modeling with Petri Nets has been implemented. It uses the VT-11 graphics terminal supported on the PDP-11/35 computer to draw, execute, analyze, edit and redraw a Petri Net. Each of the above mentioned tasks can be performed by selecting appropriate items from a menu displayed on the screen. Petri Nets with a reasonably large number of nodes can be created and analyzed using this package. The number of nodes supported may be increased by making simple changes in the program. Being interactive, the program seeks information from the user after displaying appropriate messages on the terminal. After completing the Petri Net, it may be executed step by step and the changes in the number of tokens may be observed on the screen, at each place. Some properties of Petri Nets like safety, boundedness, conservation and redundancy can be checked using this package. This package can be used very effectively for modeling asynchronous (concurrent) systems with Petri Nets and simulating the model by “graphical execution.”

The fully implicit stochastic-alpha method for stiff stochastic differential equations

A fully implicit integration method for stochastic differential equations with significant multiplicative noise and stiffness in both the drift and diffusion coefficients has been constructed, analyzed and illustrated with numerical examples in this work. The method has strong order 1.0 consistency and has user-selectable parameters that allow the user to expand the stability region of the method to cover almost the entire drift-diffusion stability plane. The large stability region enables the method to take computationally efficient time steps. A system of chemical Langevin equations simulated with the method illustrates its computational efficiency.

Pi-turns in proteins and peptides: Classification, conformation, occurrence, hydration and sequence.

The i + 5-->i hydrogen bonded turn conformation (pi-turn) with the fifth residue adopting alpha L conformation is frequently found at the C-terminus of helices in proteins and hence is speculated to be a "helix termination signal." An analysis of the occurrence of i + 5-->i hydrogen bonded turn conformation at any general position in proteins (not specifically at the helix C-terminus), using coordinates of 228 protein crystal structures determined by X-ray crystallography to better than 2.5 A resolution is reported in this paper. Of 486 detected pi-turn conformations, 367 have the (i + 4)th residue in alpha L conformation, generally occurring at the C-terminus of alpha-helices, consistent with previous observations. However, a significant number (111) of pi-turn conformations occur with (i + 4)th residue in alpha R conformation also, generally occurring in alpha-helices as distortions either at the terminii or at the middle, a novel finding. These two sets of pi-turn conformations are referred to by the names pi alpha L and pi alpha R-turns, respectively, depending upon whether the (i + 4)th residue adopts alpha L or alpha R conformations. Four pi-turns, named pi alpha L'-turns, were noticed to be mirror images of pi alpha L-turns, and four more pi-turns, which have the (i + 4)th residue in beta conformation and denoted as pi beta-turns, occur as a part of hairpin bend connecting twisted beta-strands. Consecutive pi-turns occur, but only with pi alpha R-turns. The preference for amino acid residues is different in pi alpha L and pi alpha R-turns. However, both show a preference for Pro after the C-termini. Hydrophilic residues are preferred at positions i + 1, i + 2, and i + 3 of pi alpha L-turns, whereas positions i and i + 5 prefer hydrophobic residues. Residue i + 4 in pi alpha L-turns is mainly Gly and less often Asn. Although pi alpha R-turns generally occur as distortions in helices, their amino acid preference is different from that of helices. Poor helix formers, such as His, Tyr, and Asn, also were found to be preferred for pi alpha R-turns, whereas good helix former Ala is not preferred. pi-Turns in peptides provide a picture of the pi-turn at atomic resolution. Only nine peptide-based pi-turns are reported so far, and all of them belong to pi alpha L-turn type with an achiral residue in position i + 4. The results are of importance for structure prediction, modeling, and de novo design of proteins.

Negative nuclear Overhauser effects as probes of macromolecular structure

Nuclear Overhauser effects (NOE) and circular dichroism (CD) techniques have been used to probe @-turn conformations in acyclic and cyclic peptides containingPro-Xsequences. The model peptides studied are of the type Piv-Pro-X-NHMe (X = Aib, D-Ala, Gly, Val, and Leu) and Boc-Cys-Pro-X-C s NHMe (X = Aib, L-Ala, D-Ala, Gly, and Leu). In the acyclic series, observation of NOES between Pro C"H and X-NH, together with solvent and temperature dependence of NH chemical shifts, establishes a 4 - 1 hydrogen bond stabilized type I1 @-turn in the Gly, D-Ala, and Aib peptides, in CDC13 and (CD3)2S0. A positive n-r* CD band at -225-230 nm appears to be characteristic of this structure. For the acyclic Pro-Leu peptide the observation of NOE's for both Pro and Leu C"H resonances on saturation of Leu NH is compatible with a type V bend or consecutive y-turn conformation. In the cyclic disulfide series the Pro-Aib and Pro-D-Ala peptides favor type I1 @-turns, whereas all other peptides adopt type I (111) conformations. All the cyclic disulfides exhibit an intense negative CD band at -228-230 nm. The results suggest thatgeneralcorrelations between CD spectral type and specific 0-turn conformations may not be obtained. Evidence for solvent-dependent structural changes in the Pro-Aib sequence in both cyclic and acyclic peptides is presented.

Transduction in mycobacterium smegmatis

SEVERAL people have studied lysogeny and lysogenic conversion in the advanced prokaryote, Mycobacterium1-3, and presented evidence for reciprocal genetic exchange from the lysogenic to the non-lysogenic state. Biochemical and physiological changes induced by lysogenization of M. phlei with mycobacteriophage B2 have been demonstrated by Juhasz4. We wish to report the isolation of a mycobacteriophage capable of mediating transduction in Mycobacterium smegmatis strain SN2, thus providing an opportunity to study the genetics of mycobacteria.