An L-system-based plant modeling language

Cpfg is a program for simulating and visualizing plant development, based on the theory of L-systems. A special-purpose programming language, used to specify plant models, is an essential feature of cpfg. We review postulates of L-system theory that have influenced the design of this language. We then present the main constructs of this language, and evaluate it from a user's perspective.

Physiologically based synthetic models of hepatic disposition

Current Physiologically based pharmacokinetic (PBPK) models are inductive. We present an additional, different approach that is based on the synthetic rather than the inductive approach to modeling and simulation. It relies on object-oriented programming A model of the referent system in its experimental context is synthesized by assembling objects that represent components such as molecules, cells, aspects of tissue architecture, catheters, etc. The single pass perfused rat liver has been well described in evaluating hepatic drug pharmacokinetics (PK) and is the system on which we focus. In silico experiments begin with administration of objects representing actual compounds. Data are collected in a manner analogous to that in the referent PK experiments. The synthetic modeling method allows for recognition and representation of discrete event and discrete time processes, as well as heterogeneity in organization, function, and spatial effects. An application is developed for sucrose and antipyrine, administered separately and together PBPK modeling has made extensive progress in characterizing abstracted PK properties but this has also been its limitation. Now, other important questions and possible extensions emerge. How are these PK properties and the observed behaviors generated? The inherent heuristic limitations of traditional models have hindered getting meaningful, detailed answers to such questions. Synthetic models of the type described here are specifically intended to help answer such questions. Analogous to wet-lab experimental models, they retain their applicability even when broken apart into sub-components. Having and applying this new class of models along with traditional PK modeling methods is expected to increase the productivity of pharmaceutical research at all levels that make use of modeling and simulation.

Development of disabled, replication-defective gene transfer vectors from the Jembrana disease virus, a new infectious agent of cattle

Jembrana disease virus (JDV) is a newly isolated and characterised bovine lentivirus. It causes an acute disease in Ball cattle (Bos javanicus). which can be readily transmitted to susceptible cattle with 17% mortality. There is as yet no treatment or preventive vaccine. We have developed a gene transfer vector system based on JDV that has three components. The first of the components is a bicistronic transfer vector plasmid that was constructed to contain cis-sequences from the JDV genome, including 5 '- and 3 ' -long terminal repeats (LTRs), 0.4 kb of truncated gag and 1.1 kb of 3 ' -env, a multiple cloning site to accommodate the gene(s) of interest for transfer, and an internal ribosome entry site plus the neomycin phosphotransferase (Neo) gene cassette for antibiotic selection. The second element is a packaging plasmid that contains trans-sequences. including gag, pol. vif, tar and rev: but without the env and packaging signals. The third is a plasmid encoding the G glycoprotein of vesicular stomatitis virus (VSV-G) to supply the vector an envelope for pseudotyping. Cotransfection of 293T cells with these three plasmid components produced VSV-G pseudotyped. disabled, replication defective, bicistronic JDV vectors encoding the green fluorescent protein (EGFP) and the Neo resistance selection maker simultaneously with a titre range of (0.4-1.2) x 10(6) CFU/ml. Transduction of several replicating primary and transformed cells from cattle, primate and human sources and importantly growth-arrested cells with the JDV vectors showed high efficiency of EGFP gene transfer at 35-75%, which was stable and the expression of EGFP was long term. Furthermore, these JDV vectors were designed to suit the inclusion and expression of genes corresponding to JDV specific proteins, such as gag or env, for the development of vaccines for Jembrana disease. This strategy should also be applicable to other bovine diseases as wall. The design and construction of the JDV vector system should facilitate the study of the lentivirology and pathogenesis of the diseases associated with JDV or other bovine virus infections. To our knowledge, this is the first such vector system developed from a cattle virus. (C) 2001 Elsevier Science B.V. All rights reserved.

A framework for table driven testing of Java classes

With the advent of object-oriented languages and the portability of Java, the development and use of class libraries has become widespread. Effective class reuse depends on class reliability which in turn depends on thorough testing. This paper describes a class testing approach based on modeling each test case with a tuple and then generating large numbers of tuples to thoroughly cover an input space with many interesting combinations of values. The testing approach is supported by the Roast framework for the testing of Java classes. Roast provides automated tuple generation based on boundary values, unit operations that support driver standardization, and test case templates used for code generation. Roast produces thorough, compact test drivers with low development and maintenance cost. The framework and tool support are illustrated on a number of non-trivial classes, including a graphical user interface policy manager. Quantitative results are presented to substantiate the practicality and effectiveness of the approach. Copyright (C) 2002 John Wiley Sons, Ltd.

Insights into the structural basis for zinc inhibition of the glycine receptor

Histidines 107 and 109 in the glycine receptor ( GlyR) alpha(1) subunit have previously been identified as determinants of the inhibitory zinc-binding site. Based on modeling of the GlyR alpha(1) subunit extracellular domain by homology to the acetylcholine-binding protein crystal structure, we hypothesized that inhibitory zinc is bound within the vestibule lumen at subunit interfaces, where it is ligated by His(107) from one subunit and His(109) from an adjacent subunit. This was tested by co-expressing alpha(1) subunits containing the H107A mutation with alpha(1) subunits containing the H109A mutation. Although sensitivity to zinc inhibition is markedly reduced when either mutation is individually incorporated into all five subunits, the GlyRs formed by the co-expression of H107A mutant subunits with H109A mutant subunits exhibited an inhibitory zinc sensitivity similar to that of the wild type alpha(1) homomeric GlyR. This constitutes strong evidence that inhibitory zinc is coordinated at the interface between adjacent alpha(1) subunits. No evidence was found for beta subunit involvement in the coordination of inhibitory zinc, indicating that a maximum of two zinc-binding sites per alpha(1)beta receptor is sufficient for maximal zinc inhibition. Our data also show that two zinc-binding sites are sufficient for significant inhibition of alpha(1) homomers. The binding of zinc at the interface between adjacent alpha(1) subunits could restrict intersubunit movements, providing a feasible mechanism for the inhibition of channel activation by zinc.

Active site mutations and substrate inhibition in human sulfotransferase 1A1 and 1A3

Human SULT1A1 is primarily responsible for sulfonation of xenobiotics, including the activation of promutagens, and it has been implicated in several forms of cancer. Human SULT1A3 has been shown to be the major sulfotransferase that sulfonates dopamine. These two enzymes shares 93% amino acid sequence identity and have distinct but overlapping substrate preferences. The resolution of the crystal structures of these two enzymes has enabled us to elucidate the mechanisms controlling their substrate preferences and inhibition. The presence of two p-nitrophenol (pNP) molecules in the crystal structure of SULT1A1 was postulated to explain cooperativity at low and inhibition at high substrate concentrations, respectively. In SULT1A1, substrate inhibition occurs with pNP as the substrate but not with dopamine. For SULT1A3, substrate inhibition is found for dopamine but not with pNP. We investigated how substrate inhibition occurs in these two enzymes using molecular modeling, site-directed mutagenesis, and kinetic analysis. The results show that residue Phe-247 of SULT1A1, which interacts with both p-nitrophenol molecules in the active site, is important for substrate inhibition. Mutation of phenylalanine to leucine at this position in SULT1A1 results in substrate inhibition by dopamine. We also propose, based on modeling and kinetic studies, that substrate inhibition by dopamine in SULT1A3 is caused by binding of two dopamine molecules in the active site. © 2004 by The American Society for Biochemistry and Molecular Biology, Inc.

3D protein structure matching by patch signatures

For determining functionality dependencies between two proteins, both represented as 3D structures, it is an essential condition that they have one or more matching structural regions called patches. As 3D structures for proteins are large, complex and constantly evolving, it is computationally expensive and very time-consuming to identify possible locations and sizes of patches for a given protein against a large protein database. In this paper, we address a vector space based representation for protein structures, where a patch is formed by the vectors within the region. Based on our previews work, a compact representation of the patch named patch signature is applied here. A similarity measure of two patches is then derived based on their signatures. To achieve fast patch matching in large protein databases, a match-and-expand strategy is proposed. Given a query patch, a set of small k-sized matching patches, called candidate patches, is generated in match stage. The candidate patches are further filtered by enlarging k in expand stage. Our extensive experimental results demonstrate encouraging performances with respect to this biologically critical but previously computationally prohibitive problem.

Off-line signature verification and forgery detection system based on fuzzy modeling

This paper presents an innovative approach for signature verification and forgery detection based on fuzzy modeling. The signature image is binarized and resized to a fixed size window and is then thinned. The thinned image is then partitioned into a fixed number of eight sub-images called boxes. This partition is done using the horizontal density approximation approach. Each sub-image is then further resized and again partitioned into twelve further sub-images using the uniform partitioning approach. The features of consideration are normalized vector angle (α) from each box. Each feature extracted from sample signatures gives rise to a fuzzy set. Since the choice of a proper fuzzification function is crucial for verification, we have devised a new fuzzification function with structural parameters, which is able to adapt to the variations in fuzzy sets. This function is employed to develop a complete forgery detection and verification system.

Fuzzy modeling based recognition of multi-font numerals

In this paper, we present a new scheme for off-line recognition of multi-font numerals using the Takagi-Sugeno (TS) model. In this scheme, the binary image of a character is partitioned into a fixed number of sub-images called boxes. The features consist of normalized vector distances (gamma) from each box. Each feature extracted from different fonts gives rise to a fuzzy set. However, when we have a small number of fonts as in the case of multi-font numerals, the choice of a proper fuzzification function is crucial. Hence, we have devised a new fuzzification function involving parameters, which take account of the variations in the fuzzy sets. The new fuzzification function is employed in the TS model for the recognition of multi-font numerals.

Gene transfer and expression of a non-viral polycation-based vector in CD4(+) cells

CD4-selective targeting of an antibody-polycation-DNA complex was investigated The complex was synthesized with the anti-CD4 monoclonal antibody B-F5, polylysine(268) (pLL) and either the pGL3 control vector containing the luciferase reporter gene or the pGeneGrip vector containing the green fluorescent protein (GFP) gene. B-F5-pLL-DNA complexes inhibited the binding of I-125-B-F5 to CD4(+) Jurkat cells, while complexes synthesised either without B-F5 or using a non-specific mouse IgG1 antibody had little or no effect Expression of the luciferase reporter gene was achieved in Jurkat cells using the B-F5-pLL-pGL3 complex and was enhanced in the presence of PMA. Negligible luciferase activity was defected with the non-specific antibody complex in Jurkat cells or with the B-F5-pLL-pGL3 complex in the CD4(-) K-562 cells. Using complexes synthesised with the pGeneGrip vector, the transfection efficiency in Jurkat and K-562 cells was examined using confocal microscopy. More than 95% of Jurkat cells expressed GFP and the level of this expression was markedly enhanced by PMA. Negligible GFP expression was seen in K-562 cells or when B-F5 was replaced by a nonspecific antibody. Using flow cytometry, fluorescein-labelled complex showed specific targeting to CD4(+) cells in a mixed cell population from human peripheral blood. These studies demonstrate the selective transfection of CD4(+) T-lymphoid cells using a polycation-based gene delivery system. The complex may provide a means of delivering anti-HIV gene therapies to CD4(+) cells in vivo.

Modeling of hepatic elimination and organ distribution kinetics with the extended convection-dispersion model

The conventional convection-dispersion (also called axial dispersion) model is widely used to interrelate hepatic availability (F) and clearance (Cl) with the morphology and physiology of the liver and to predict effects such as changes in liver blood flow on F and Cl. An extended form of the convection-dispersion model has been developed to adequately describe the outflow concentration-time profiles for vascular markers at both short and long times after bolus injections into perfused livers. The model, based on flux concentration and a convolution of catheters and large vessels, assumes that solute elimination in hepatocytes follows either fast distribution into or radial diffusion in hepatocytes. The model includes a secondary vascular compartment, postulated to be interconnecting sinusoids. Analysis of the mean hepatic transit time (MTT) and normalized variance (CV2) of solutes with extraction showed that the discrepancy between the predictions of MTT and CV2 for the extended and conventional models are essentially identical irrespective of the magnitude of rate constants representing permeability, volume, and clearance parameters, providing that there is significant hepatic extraction. In conclusion, the application of a newly developed extended convection-dispersion model has shown that the unweighted conventional convection-dispersion model can be used to describe the disposition of extracted solutes and, in particular, to estimate hepatic availability and clearance in booth experimental and clinical situations.

Simulated-annealing-based genetic algorithm for modeling the optical constants of solids

We propose a simulated-annealing-based genetic algorithm for solving model parameter estimation problems. The algorithm incorporates advantages of both genetic algorithms and simulated annealing. Tests on computer-generated synthetic data that closely resemble optical constants of a metal were performed to compare the efficiency of plain genetic algorithms against the simulated-annealing-based genetic algorithms. These tests assess the ability of the algorithms to and the global minimum and the accuracy of values obtained for model parameters. Finally, the algorithm with the best performance is used to fit the model dielectric function to data for platinum and aluminum. (C) 1997 Optical Society of America.

Modeling of phase equilibrium and vapor adsorption on carbon black based on a combination of a lattice theory and equation of state

A thermodynamic approach is developed in this paper to describe the behavior of a subcritical fluid in the neighborhood of vapor-liquid interface and close to a graphite surface. The fluid is modeled as a system of parallel molecular layers. The Helmholtz free energy of the fluid is expressed as the sum of the intrinsic Helmholtz free energies of separate layers and the potential energy of their mutual interactions calculated by the 10-4 potential. This Helmholtz free energy is described by an equation of state (such as the Bender or Peng-Robinson equation), which allows us a convenient means to obtain the intrinsic Helmholtz free energy of each molecular layer as a function of its two-dimensional density. All molecular layers of the bulk fluid are in mechanical equilibrium corresponding to the minimum of the total potential energy. In the case of adsorption the external potential exerted by the graphite layers is added to the free energy. The state of the interface zone between the liquid and the vapor phases or the state of the adsorbed phase is determined by the minimum of the grand potential. In the case of phase equilibrium the approach leads to the distribution of density and pressure over the transition zone. The interrelation between the collision diameter and the potential well depth was determined by the surface tension. It was shown that the distance between neighboring molecular layers substantially changes in the vapor-liquid transition zone and in the adsorbed phase with loading. The approach is considered in this paper for the case of adsorption of argon and nitrogen on carbon black. In both cases an excellent agreement with the experimental data was achieved without additional assumptions and fitting parameters, except for the fluid-solid potential well depth. The approach has far-reaching consequences and can be readily extended to the model of adsorption in slit pores of carbonaceous materials and to the analysis of multicomponent adsorption systems. (C) 2002 Elsevier Science (USA).

Modeling of gas adsorption equilibrium over a wide range of pressure: A thermodynamic approach based on equation of state

A thermodynamic approach based on the Bender equation of state is suggested for the analysis of supercritical gas adsorption on activated carbons at high pressure. The approach accounts for the equality of the chemical potential in the adsorbed phase and that in the corresponding bulk phase and the distribution of elements of the adsorption volume (EAV) over the potential energy for gas-solid interaction. This scheme is extended to subcritical fluid adsorption and takes into account the phase transition in EAV The method is adapted to gravimetric measurements of mass excess adsorption and has been applied to the adsorption of argon, nitrogen, methane, ethane, carbon dioxide, and helium on activated carbon Norit R I in the temperature range from 25 to 70 C. The distribution function of adsorption volume elements over potentials exhibits overlapping peaks and is consistently reproduced for different gases. It was found that the distribution function changes weakly with temperature, which was confirmed by its comparison with the distribution function obtained by the same method using nitrogen adsorption isotherm at 77 K. It was shown that parameters such as pore volume and skeleton density can be determined directly from adsorption measurements, while the conventional approach of helium expansion at room temperature can lead to erroneous results due to the adsorption of helium in small pores of activated carbon. The approach is a convenient tool for analysis and correlation of excess adsorption isotherms over a wide range of pressure and temperature. This approach can be readily extended to the analysis of multicomponent adsorption systems. (C) 2002 Elsevier Science (USA).

The modeling of turbulent reactive flows based on multiple mapping conditioning

A new modeling approach-multiple mapping conditioning (MMC)-is introduced to treat mixing and reaction in turbulent flows. The model combines the advantages of the probability density function and the conditional moment closure methods and is based on a certain generalization of the mapping closure concept. An equivalent stochastic formulation of the MMC model is given. The validity of the closuring hypothesis of the model is demonstrated by a comparison with direct numerical simulation results for the three-stream mixing problem. (C) 2003 American Institute of Physics.