Biomimetic in silico devices

We introduce biomimetic in silico devices, and means for validation along with methods for testing and refining them. The devices are constructed from adaptable software components designed to map logically to biological components at multiple levels of resolution. In this report we focus on the liver; the goal is to validate components that mimic features of the lobule (the hepatic primary functional unit) and dynamic aspects of liver behavior, structure, and function. An assembly of lobule-mimetic devices represents an in silico liver. We validate against outflow profiles for sucrose administered as a bolus to isolated, perfused rat livers. Acceptable in silico profiles are experimentally indistinguishable from those of the in situ referent. This new technology is intended to provide powerful Dew tools for challenging our understanding of how biological functional units function in vivo.

Disruption of NaS1 sulfate transport function in mice leads to enhanced acetaminophen-induced hepatotoxicity

Sulfate is required for detoxification of xenobiotics such as acetaminophen (APAP), a leading cause of liver failure in humans. The NaS1 sulfate transporter maintains blood sulfate levels sufficiently high for sulforiation reactions to work effectively for drug detoxification. In the present study, we identified two loss-of-function polymorphisms in the human NaS1 gene and showed the Nas1-null mouse to be hypersensitive to APAP hepatotoxicity. APAP treatment led to increased liver damage and decreased hepatic glutathione levels in the hyposulfatemic Nas1-null mice compared with that in normosulfatemic wild-type mice. Analysis of urinary APAP metabolites revealed a significantly lower ratio of APAP-sulfate to APAP-glucuronide in the Nas1-null mice. These results suggest hyposulfatemia increases sensitivity to APAP-induced hepatotoxicity by decreasing the sulfonation capacity to metabolize APAP. In conclusion, the results of this study highlight the importance of plasma sulfate level as a key modulator of acetaminophen metabolism and suggest that individuals with reduced NaS1 sulfate transporter function would be more sensitive to hepatotoxic agents.

A stochastic frontier production function with flexible risk properties

This paper considers a stochastic frontier production function which has additive, heteroscedastic error structure. The model allows for negative or positive marginal production risks of inputs, as originally proposed by Just and Pope (1978). The technical efficiencies of individual firms in the sample are a function of the levels of the input variables in the stochastic frontier, in addition to the technical inefficiency effects. These are two features of the model which are not exhibited by the commonly used stochastic frontiers with multiplicative error structures, An empirical application is presented using cross-sectional data on Ethiopian peasant farmers. The null hypothesis of no technical inefficiencies of production among these farmers is accepted. Further, the flexible risk models do not fit the data on peasant farmers as well as the traditional stochastic frontier model with multiplicative error structure.

The structure of versutoxin (delta-atracotoxin-Hv1) provides insights into the binding of site 3 neurotoxins to the voltage-gated sodium channel

Background: Versutoxin (delta-ACTX-Hv1) is the major component of the venom of the Australian Blue Mountains funnel web spider, Hadronyche versuta. delta-ACTX-Hv1 produces potentially fatal neurotoxic symptoms in primates by slowing the inactivation of voltage-gated sodium channels; delta-ACTX-Hv1 is therefore a useful tool for studying sodium channel function. We have determined the three-dimensional structure of delta ACTX-Hv1 as the first step towards understanding the molecular basis of its interaction with these channels. Results: The solution structure of delta-ACTX-Hv1, determined using NMR spectroscopy, comprises a core beta region containing a triple-stranded antiparallel beta sheet, a thumb-like extension protruding from the beta region and a C-terminal 3(10) helix that is appended to the beta domain by virtue of a disulphide bond. The beta region contains a cystine knot motif similar to that seen in other neurotoxic polypeptides. The structure shows homology with mu-agatoxin-l, a spider toxin that also modifies the inactivation kinetics of vertebrate voltage-gated sodium channels. More surprisingly, delta-ACTX-Hv1 shows both sequence and structural homology with gurmarin, a plant polypeptide. This similarity leads us to suggest that the sweet-taste suppression elicited by gurmarin may result from an interaction with one of the downstream ion channels involved in sweet-taste transduction. Conclusions: delta-ACTX-Hv1 shows no structural homology with either sea anemone or alpha-scorpion toxins, both of which also modify the inactivation kinetics of voltage-gated sodium channels by interacting with channel recognition site 3. However, we have shown that delta-ACTX-Hv1 contains charged residues that are topologically related to those implicated in the binding of sea anemone and alpha-scorpion toxins to mammalian voltage-gated sodium channels, suggesting similarities in their mode of interaction with these channels.

Tyrosines in the Carboxyl Terminus Regulate Syk Kinase Activity and Function

The Syk tyrosine kinase family plays an essential role in immunoreceptor tyrosine-based activation motif (ITAM) signaling. The binding of Syk to tyrosine-phosphorylated ITAM subunits of immunoreceptors, such as Fc epsilon RI on mast cells, results in a conformational change, with an increase of enzymatic activity of Syk. This conformational change exposes the COOH-terminal tail of Syk, which has three conserved Tyr residues (Tyr-623, Tyr-624, and Tyr-625 of rat Syk). To understand the role of these residues in signaling, wild-type and mutant Syk with these three Tyr mutated to Phe was expressed in Syk-deficient mast cells. There was decreased Fc epsilon RI-induced degranulation, nuclear factor for T cell activation and NF kappa B activation with the mutated Syk together with reduced phosphorylation of MAP kinases p38 and p42/44 ERK. In non-stimulated cells, the mutated Syk was more tyrosine phosphorylated predominantly as a result of autophosphorylation. In vitro, there was reduced binding of mutated Syk to phosphorylated ITAM due to this increased phosphorylation. This mutated Syk from non-stimulated cells had significantly reduced kinase activity toward an exogenous substrate, whereas its autophosphorylation capacity was not affected. However, the kinase activity and the autophosphorylation capacity of this mutated Syk were dramatically decreased when the protein was dephosphorylated before the in vitro kinase reaction. Furthermore, mutation of these tyrosines in the COOH-terminal region of Syk transforms it to an enzyme, similar to its homolog ZAP-70, which depends on other tyrosine kinases for optimal activation. In testing Syk mutated singly at each one of the tyrosines, Tyr-624 but especially Tyr-625 had the major role in these reactions. Therefore, these results indicate that these tyrosines in the tail region play a critical role in regulating the kinase activity and function of Syk.

Gangliosides Are Important for the Preservation of the Structure and Organization of RBL-2H3 Mast Cells

Gangliosides are known to be important in many biological processes. However, details concerning the exact function of these glycosphingolipids in cell physiology are poorly understood. in this study, the role of gangliosides present on the surface of rodent mast cells in maintaining cell structure was examined using RBL-2H3 mast cells and two mutant cell lines (E5 and D1) deficient in the gangliosides, GM(1) and the alpha-galactosyl derivatives of the ganglioside GD(1b). The two deficient cell lines were morphologically different from each other as well as from the parental RBL-2H3 cells. Actin filaments in RBL-2H3 and E5 cells were under the plasma membrane following the spindle shape of the cells, whereas in D1 cells, they were concentrated in large membrane ruffles. Microtubules in RBL-2H3 and E5 cells radiated from the centrosome and were organized into long, straight bundles. The bundles in D1 cells were thicker and organized circumferentially under the plasma membrane. The endoplasmic reticulum, the Golgi complex, and the secretory granule matrix were also altered in the mutant cell lines. These results suggest that the mast cell-specific alpha-galactosyl derivatives of ganglioside GD(1b) and GM(1) are important in maintaining normal cell morphology. (J Histochern Cytochem 58:83-93, 2010)

Neonatal endotoxin exposure changes neuroendocrine, cardiovascular function and mortality during polymicrobial sepsis in adult rats

Our aim was to investigate whether neonatal LPS challenge may improve hormonal, cardiovascular response and mortality, this being a beneficial adaptation when adult rats are submitted to polymicrobial sepsis by cecal ligation and puncture (CLP). Fourteen days after birth, pups received an intraperitoneal injection of lipopolysaccharide (LPS; 100 mu g/kg) or saline. After 8-12 weeks, they were submitted to CLP, decapitated 4,6 or 24 h after surgery and blood was collected for vasopressin (AVP), corticosterone and nitrate measurement, while AVP contents were measured in neurohypophysis, supra-optic (SON) and paraventricular (PVN) nuclei. Moreover, rats had their mean arterial pressure (MAP) and heart rate (HR) evaluated, and mortality and bacteremia were determined at 24 h. Septic animals with neonatal LPS exposure had higher plasma AVP and corticosterone levels, and higher c-Fos expression in SON and PVN at 24 h after surgery when compared to saline treated rats. The LPS pretreated group showed increased AVP content in SON and PVN at 6 h, while we did not observe any change in neurohypophyseal AVP content. The nitrate levels were significantly reduced in plasma at 6 and 24 h after surgery, and in both hypothalamic nuclei only at 6 h. Septic animals with neonatal LPS exposure showed increase in MAP during the initial phase of sepsis, but HR was not different from the neonatal saline group. Furthermore, neonatally LPS exposed rats showed a significant decrease in mortality rate as well as in bacteremia. These data suggest that neonatal LPS challenge is able to promote beneficial effects on neuroendocrine and cardiovascular responses to polymicrobial sepsis in adulthood. (C) 2011 Elsevier B.V. All rights reserved.

Integrating high and moderate spatial resolution image data to estimate forest age structure

The collection of spatial information to quantify changes to the state and condition of the environment is a fundamental component of conservation or sustainable utilization of tropical and subtropical forests, Age is an important structural attribute of old-growth forests influencing biological diversity in Australia eucalypt forests. Aerial photograph interpretation has traditionally been used for mapping the age and structure of forest stands. However this method is subjective and is not able to accurately capture fine to landscape scale variation necessary for ecological studies. Identification and mapping of fine to landscape scale vegetative structural attributes will allow the compilation of information associated with Montreal Process indicators lb and ld, which seek to determine linkages between age structure and the diversity and abundance of forest fauna populations. This project integrated measurements of structural attributes derived from a canopy-height elevation model with results from a geometrical-optical/spectral mixture analysis model to map forest age structure at a landscape scale. The availability of multiple-scale data allows the transfer of high-resolution attributes to landscape scale monitoring. Multispectral image data were obtained from a DMSV (Digital Multi-Spectral Video) sensor over St Mary's State Forest in Southeast Queensland, Australia. Local scene variance levels for different forest tapes calculated from the DMSV data were used to optimize the tree density and canopy size output in a geometric-optical model applied to a Landsat Thematic Mapper (TU) data set. Airborne laser scanner data obtained over the project area were used to calibrate a digital filter to extract tree heights from a digital elevation model that was derived from scanned colour stereopairs. The modelled estimates of tree height, crown size, and tree density were used to produce a decision-tree classification of forest successional stage at a landscape scale. The results obtained (72% accuracy), were limited in validation, but demonstrate potential for using the multi-scale methodology to provide spatial information for forestry policy objectives (ie., monitoring forest age structure).

Structure-hepatic disposition relationships for cationic drugs in isolated perfused rat livers: Transmembrane exchange and cytoplasmic binding process

This work studied the structure-hepatic disposition relationships for cationic drugs of varying lipophilicity using a single-pass, in situ rat liver preparation. The lipophilicity among the cationic drugs studied in this work is in the following order: diltiazem. propranolol. labetalol. prazosin. antipyrine. atenolol. Parameters characterizing the hepatic distribution and elimination kinetics of the drugs were estimated using the multiple indicator dilution method. The kinetic model used to describe drug transport (the two-phase stochastic model) integrated cytoplasmic binding kinetics and belongs to the class of barrier-limited and space-distributed liver models. Hepatic extraction ratio (E) (0.30-0.92) increased with lipophilicity. The intracellular binding rate constant (k(on)) and the equilibrium amount ratios characterizing the slowly and rapidly equilibrating binding sites (K-S and K-R) increase with the lipophilicity of drug (k(on) : 0.05-0.35 s(-1); K-S : 0.61-16.67; K-R : 0.36-0.95), whereas the intracellular unbinding rate constant (k(off)) decreases with the lipophilicity of drug (0.081-0.021 s(-1)). The partition ratio of influx (k(in)) and efflux rate constant (k(out)), k(in)/k(out), increases with increasing pK(a) value of the drug [from 1.72 for antipyrine (pK(a) = 1.45) to 9.76 for propranolol (pK(a) = 9.45)], the differences in k(in/kout) for the different drugs mainly arising from ion trapping in the mitochondria and lysosomes. The value of intrinsic elimination clearance (CLint), permeation clearance (CLpT), and permeability-surface area product (PS) all increase with the lipophilicity of drug [CLint (ml . min(-1) . g(-1) of liver): 10.08-67.41; CLpT (ml . min(-1) . g(-1) of liver): 10.80-5.35; PS (ml . min(-1) . g(-1) of liver): 14.59-90.54]. It is concluded that cationic drug kinetics in the liver can be modeled using models that integrate the presence of cytoplasmic binding, a hepatocyte barrier, and a vascular transit density function.

The solution structure of C1-T1, a two-domain proteinase inhibitor derived from a circular precursor protein from Nicotiana alata

A two-domain portion of the proteinase inhibitor precursor from Nicotiana alata (NaProPI) has been expressed and its structure determined by NMR spectroscopy. NaProPI contains six almost identical 53 amino acid repeats that fold into six highly similar domains; however, the sequence repeats do nut coincide with the structural domains. Five of the structural domains comprise the C-terminal portion of one repeat and the N-terminal portion of the next. The sixth domain contains the C-terminal portion of the sixth repeat and the N-terminal portion of the first repeat. Disulphide bonds link these C and N-terminal fragments to generate the clasped-bracelet fold of NaProPI. The three-dimensional structure of NaProPI is not known, but it is conceivable that adjacent domains in NaProPI interact to generate the circular bracelet with the N and C termini in close enough proximity to facilitate formation of the disulphide bonds that form the clasp The expressed protein, examined in the current study, comprises residues 25-135 of NaProPI and encompasses the first two contiguous structural domains, namely the chymotrypsin inhibitor C1 and the trypsin inhibitor T1, joined by a five-residue linker, and is referred to as C1-T1. The tertiary structure of each domain in C1-T1 is identical to that found in the isolated inhibitors. However, no nuclear Overhauser effect contacts are observed between the two domains and the five-residue linker adopts an extended conformation. The absence of interactions between the domains indicates that adjacent domains do not specifically interact to drive the circularisation of NaProPI. These results are in agreement with recent data which describe similar PI precursors from other members of the Solanaceae having two, three, or four repeats. The lack of strong interdomain association is likely to be important for the function of individual inhibitors by ensuring that there is no masking of reactive sites upon release from the precursor. (C) 2001 Academic Press.

Growth regime map for liquid-bound granules: further development and experimental validation

An attempt was made to quantify the boundaries and validate the granule growth regime map for liquid-bound granules recently proposed by Iveson and Litster (AlChE J. 44 (1998) 1510). This regime map postulates that the type of granule growth behaviour is a function of only two dimensionless groups: the amount of granule deformation during collision (characterised by a Stokes deformation number, St(def)) and the maximum granule pore saturation, s(max). The results of experiments performed with a range of materials (glass ballotini, iron ore fines, copper chalcopyrite powder and a sodium sulphate and cellulose mixture) using both drum and high shear mixer granulators were examined. The drum granulation results gave good agreement with the proposed regime map. The boundary between crumb and steady growth occurs at St(def) of order 0.1 and the boundary between steady and induction growth occurs at St(def) of order 0.001. The nucleation only boundary occurs at pore saturations that increase from 70% to 80% with decreasing St(def). However, the high shear mixer results all had St(def) numbers which were too large. This is most likely to be because the chopper tip-speed is an over-estimate of the average impact velocity granules experience and possibly also due to the dynamic yield strength of the materials being significantly greater than the yield strengths measured at low strain rates. Hence, the map is only a useful tool for comparing the granulation behaviour of different materials in the same device. Until we have a better understanding of the flow patterns and impact velocities in granulators, it cannot be used to compare different types of equipment. Theoretical considerations also revealed that several of the regime boundaries are also functions of additional parameters not explicitly contained on the map, such as binder viscosity. (C) 2001 Elsevier Science B.V. All rights reserved.

Identification of residues and domains of Raf important for function in vivo and in vitro

Random mutagenesis and genetic screens for impaired Raf function in Caenorhabditis elegans were used to identify six loss-of-function alleles of lin-45 raf that result in a substitution of a single amino acid. The mutations were classified as weak, intermediate, and strong based on phenotypic severity. We engineered these mutations into the homologous residues of vertebrate Raf-1 and analyzed the mutant proteins for their underlying biochemical defects. Surprisingly, phenotype strength did not correlate with the catalytic activity of the mutant proteins. Amino acid substitutions Val-589 and Ser-619 severely compromised Raf kinase activity, yet these mutants displayed weak phenotypes in the genetic screen. Interestingly, this is because these mutant Raf proteins efficiently activate the MAPK (mitogen-activated protein kinase) cascade in living cells, a result that may inform the analysis of knockout mice. Equally intriguing was the observation that mutant proteins with non-functional Ras-binding domains, and thereby deficient in Ras-mediated membrane recruitment, displayed only intermediate strength phenotypes. This confirms that secondary mechanisms exist to couple Ras to Raf in vivo. The strongest phenotype in the genetic screens was displayed by a S508N mutation that again did not correlate with a significant loss of kinase activity or membrane recruitment by oncogenic Ras in biochemical assays. Ser-508 lies within the Raf-1 activation loop, and mutation of this residue in Raf-1 and the equivalent Ser-615 in B-Raf revealed that this residue regulates Raf binding to MEK. Further characterization revealed that in response to activation by epidermal growth factor, the Raf-S508N mutant protein displayed both reduced catalytic activity and aberrant activation kinetics: characteristics that may explain the C. elegans phenotype.

The ultrastructure and function of the silk-producing basitarsus in the Hilarini (Diptera : Ernpididae)

The tribe Hilarini (Diptera: Empididae), commonly known as dance flies, can be recognised by their swollen silk-producing prothoracic basitarsus, a male secondary sexual characteristic. The ultrastructure and function of the silk-producing basitarsus from one undescribed morphospecies of Hilarini, 'Hilarempis 20', is presented. Male H. 20 collect small parcels of diatomaceous algae from the surface of freshwater creeks that they bind with silk produced by the gland in the basitarsus. The gift is then presented to females in a nearby swarm, composed predominately of females. The basitarsus houses approximately 12 pairs of class III dermal glandular units that congregate on the ventral side of the cavity. Each gland cell has a large extracellular lumen where secretion accumulates. The lumen drains to the outside via a conducting canal encompassed by a canal cell and a duct extending through the shaft of a specialised secretory spine. The secretory spines lie in pairs in a ventral groove that runs the length of the basitarsus. A comparison of the basitarsal secretory spines with sensilla on the basitarsi of non gland-bearing legs of males, and with non gland-bearing prothoracic. basitarsi of females, suggests that the glandular units are derived from contact chemosensory sensilla. (C) 2003 Elsevier Ltd. All rights reserved.

Epitope map for a growth hormone receptor agonist monoclonal antibody, MAb 263

Monoclonal antibody (MAb) 263 is a widely used monoclonal antibody that recognizes the extracellular domain (ECD) of the GH receptor. It has been shown to act as a GH agonist both in vitro and in vivo, and we report here that it must be divalent to exert its effect on the full-length receptor. To understand the mechanism of its agonist action, we have determined the precise epitope for this antibody using a novel random PCR mutagenesis approach together with expression screening in yeast. A library of 5200 clones of rabbit GH receptor ECD mutants were screened both with MAb 263 and with an anticarboxy-tag antibody to verify complete ECD expression. Sequencing for clones that expressed complete ECD but were not MAb 263 positive identified 20 epitope residues distributed in a discontinuous manner throughout the ECD. The major part of the epitope, as revealed after mapping onto the crystal structure model of the ECD molecule, was located on the side and upper portion of domain 1, particularly within the D - E strand disulfide loop 79 - 96. Molecular dynamics docking of an antibody of the same isotype as MAb 263 was used to dock the bivalent antibody to the 1528-Angstrom(2) epitope and to visualize the likely consequences of MAb binding. The minimized model enables the antibody to grasp two receptors in a pincer-like movement from opposite sides, facilitating alignment of the receptor dimerization domains in a manner similar to, but not identical with, GH.

Conformational changes in SP-B as a function of surface pressure

X-ray reflectivity of bovine and sheep surfactant-associated protein B (SP-B) monolayers is used in conjunction with pressure-area isotherms and protein models to suggest that the protein undergoes changes in its tertiary structure at the air/water interface under the influence of surface pressure, indicating the likely importance of such changes to the phenomena of protein squeeze out as well as lipid exchange between the air-water interface and subphase structures. We describe an algorithm based on the well-established box- or layer-models that greatly assists the fitting of such unknown scattering-length density profiles, and which takes the available instrumental resolution into account. Scattering-length density profiles from neutron reflectivity of bovine SP-B monolayers on aqueous subphases are shown to be consistent with the exchange of a large number of labile protons as well as the inclusion of a significant amount of water, which is partly squeezed out of the protein monolayer at elevated surface pressures.