An approach to verifying and debugging simulation models governed by ordinary differential equations: Part 1. Methodology for residual generation

For dynamic simulations to be credible, verification of the computer code must be an integral part of the modelling process. This two-part paper describes a novel approach to verification through program testing and debugging. In Part 1, a methodology is presented for detecting and isolating coding errors using back-to-back testing. Residuals are generated by comparing the output of two independent implementations, in response to identical inputs. The key feature of the methodology is that a specially modified observer is created using one of the implementations, so as to impose an error-dependent structure on these residuals. Each error can be associated with a fixed and known subspace, permitting errors to be isolated to specific equations in the code. It is shown that the geometric properties extend to multiple errors in either one of the two implementations. Copyright (C) 2003 John Wiley Sons, Ltd.

An approach to verifying and debugging simulation models governed by ordinary differential equations: Part 2. Residuals analysis and a case study

In Part 1 of this paper a methodology for back-to-back testing of simulation software was described. Residuals with error-dependent geometric properties were generated. A set of potential coding errors was enumerated, along with a corresponding set of feature matrices, which describe the geometric properties imposed on the residuals by each of the errors. In this part of the paper, an algorithm is developed to isolate the coding errors present by analysing the residuals. A set of errors is isolated when the subspace spanned by their combined feature matrices corresponds to that of the residuals. Individual feature matrices are compared to the residuals and classified as 'definite', 'possible' or 'impossible'. The status of 'possible' errors is resolved using a dynamic subset testing algorithm. To demonstrate and validate the testing methodology presented in Part 1 and the isolation algorithm presented in Part 2, a case study is presented using a model for biological wastewater treatment. Both single and simultaneous errors that are deliberately introduced into the simulation code are correctly detected and isolated. Copyright (C) 2003 John Wiley Sons, Ltd.

Isolation and characterization of a cone snail protease with homology to CRISP proteins of the pathogenesis-related protein superfamily

The pathogenesis-related (PR) protein superfamily is widely distributed in the animal, plant, and fungal kingdoms and is implicated in human brain tumor growth and plant pathogenesis. The precise biological activity of PR proteins, however, has remained elusive. Here we report the characterization, cloning and structural homology modeling of Tex31 from the venom duct of Conus textile. Tex31 was isolated to >95% purity by activity-guided fractionation using a para-nitroanilide substrate based on the putative cleavage site residues found in the propeptide precursor of conotoxin TxVIA. Tex31 requires four residues including a leucine N-terminal of the cleavage site for efficient substrate processing. The sequence of Tex31 was determined using two degenerate PCR primers designed from N-terminal and tryptic digest Edman sequences. A BLAST search revealed that Tex31 was a member of the PR protein superfamily and most closely related to the CRISP family of mammalian proteins that have a cysteine-rich C-terminal tail. A homology model constructed from two PR proteins revealed that the likely catalytic residues in Tex31 fall within a structurally conserved domain found in PR proteins. Thus, it is possible that other PR proteins may also be substrate-specific proteases.

On the evaluation of some differential formulations for the pom-pom constitutive model

In recent years, progress has been made in modelling long chain branched polymers by the introduction of the so-called pompom model. Initially developed by McLeish and Larson (1998), the model has undergone several improvements or alterations, leading to the development of new formulations. Some of these formulations however suffer from certain mathematical defects. The purpose of the present paper is to review some of the formulations of the pom-pom constitutive model, and to investigate their possible mathematical defects. Next, an alternative formulation is proposed, which does not appear to exhibit mathematical defects, and we explore its modelling performance by comparing the predictions with experiments in non-trivial rheometric flows of an LDPE melt. The selected rheometric flows are the double step strain, as well as the large amplitude oscillatory shear experiments. For LAOS experiments, the comparison involves the use of Fourier-transform analysis.

Surface diffusion of strongly adsorbing vapors in activated carbon by a differential permeation method

Conventional methods to determine surface diffusion of adsorbed molecules are proven to be inadequate for strongly adsorbing vapors on activated carbon. Knudsen diffusion permeability (B-k) for strongly adsorbing vapors cannot be directly estimated from that of inert gases such as helium. In this paper three models are considered to elucidate the mechanism of surface diffusion in activated carbon. The transport mechanism in all three models is a combination of Knudsen diffusion, viscous flow and surface diffusion. The collision reflection factor f (which is the fraction of molecules undergoing collision to the solid surface over reflection from the surface) of the Knudsen diffusivity is assumed to be a function of loading. It was found to be 1.79 in the limit of zero loading, and decreases as loading increases. The surface diffusion permeability increases sharply at very low pressures and then starts to decrease after it has reached a maximum (B(mum)s) at a threshold pressure. The initial rapid increase in the total permeability is mainly attributed to surface diffusion. Interestingly the B(mum)s for all adsorbates appear at the same volumetric adsorbed phase concentration, suggesting that the volume of adsorbed molecules may play an important role in the surface diffusion mechanism in activated carbon. (C) 2003 Elsevier Ltd. All rights reserved.

Purification of post-translationally modified proteins from bacteria: homologous expression and purification of histidine-tagged pilin from Neisseria meningitidis

Until recently, glycosylation of proteins in prokaryotes was regarded as uncommon and thought to be limited to special cases such as S-layer proteins and some archeal outer membrane proteins. Now, there are an increasing number of reports of bacterial proteins that are glycosylated. Pilin of pathogenic Neisseria is one of the best characterised post-translation ally modified bacterial proteins, with four different types of modifications reported, including a novel glycosylation. Pilin monomers assemble to form pilus fibres, which are long protein filaments that protrude from the surface of bacterial cells and are key virulence factors. To aid in the investigation of these modifications, pure pilin is required. A number of pilin purification methods have been published, but none are appropriate for the routine purification of pilin from many different isolates. This study describes a novel, rapid, and simple method of pilin purification from Neisseria meningitidis C311#3, which facilitates the production of consistent quantities of pure, native pilin. A 6 x histidine tag was fused to the C-terminus of the pilin subunit structural gene, pilE, via homologous recombination placing the 6 x histidine-tagged allele in the chromosome of N. meningitidis C311#3. Pilin was purified under non-denaturing conditions via a two-step process using immobilised metal affinity chromatography (IMAC), followed by dye affinity chromatography. Analysis of the purified pilin confirmed that it retained both of the post-translational modifications examined. This novel approach may prove to be a generally applicable method for purification and analysis of post-translationally modified proteins in bacteria. (C) 2003 Elsevier Science (USA). All rights reserved.

Wavelength Selective a-SiC:H p-i-n/p-i-n Heterostructure for Fluorescent Proteins Detection

In this paper we present results on the optimization of multilayered a-SiC:H heterostructures that can be used as optical transducers for fluorescent proteins detection using the Fluorescence Resonance Energy Transfer approach. Double structures composed by pin based aSiC:H cells are analyzed. The color discrimination is achieved by ac photocurrent measurement under different externally applied bias. Experimental data on spectral response analysis, current-voltage characteristics and color and transmission rate discrimination are reported. An electrical model, supported by a numerical simulation gives insight into the device operation. Results show that the optimized a-SiC:H heterostructures act as voltage controlled optical filters in the visible spectrum. When the applied voltages are chosen appropriately those optical transducers can detect not only the selective excitation of specimen fluorophores, but also the subsequent weak acceptor fluorescent channel emission.

The folding of knotted proteins: insights from lattice simulations

We carry out systematic Monte Carlo simulations of Go lattice proteins to investigate and compare the folding processes of two model proteins whose native structures differ from each other due to the presence of a trefoil knot located near the terminus of one of the protein chains. We show that the folding time of the knotted fold is larger than that of the unknotted protein and that this difference in folding time is particularly striking in the temperature region below the optimal folding temperature. Both proteins display similar folding transition temperatures, which is indicative of similar thermal stabilities. By using the folding probability reaction coordinate as an estimator of folding progression we have found out that the formation of the knot is mainly a late folding event in our shallow knot system.

The expression of tubulin cofactor A (TBCA) is regulated by a noncoding antisense Tbca RNA during testis maturation

Abstract - Recently, long noncoding RNAs have emerged as pivotal molecules for the regulation of coding genes' expression. These molecules might result from antisense transcription of functional genes originating natural antisense transcripts (NATs) or from transcriptional active pseudogenes. TBCA interacts with β-tubulin and is involved in the folding and dimerization of new tubulin heterodimers, the building blocks of microtubules. Methodology/Principal findings: We found that the mouse genome contains two structurally distinct Tbca genes located in chromosomes 13 (Tbca13) and 16 (Tbca16). Interestingly, the two Tbca genes albeit ubiquitously expressed, present differential expression during mouse testis maturation. In fact, as testis maturation progresses Tbca13 mRNA levels increase progressively, while Tbca16 mRNA levels decrease. This suggests a regulatory mechanism between the two genes and prompted us to investigate the presence of the two proteins. However, using tandem mass spectrometry we were unable to identify the TBCA16 protein in testis extracts even in those corresponding to the maturation step with the highest levels of Tbca16 transcripts. These puzzling results led us to re-analyze the expression of Tbca16. We then detected that Tbca16 transcription produces sense and natural antisense transcripts. Strikingly, the specific depletion by RNAi of these transcripts leads to an increase of Tbca13 transcript levels in a mouse spermatocyte cell line. Conclusions/Significance: Our results demonstrate that Tbca13 mRNA levels are post-transcriptionally regulated by the sense and natural antisense Tbca16 mRNA levels. We propose that this regulatory mechanism operates during spermatogenesis, a process that involves microtubule rearrangements, the assembly of specific microtubule structures and requires critical TBCA levels.

Molecular imprinted nanoelectrodes for ultra sensitive detection of ovarian cancer marker

The relentless discovery of cancer biomarkers demands improved methods for their detection. In this work, we developed protein imprinted polymer on three-dimensional gold nanoelectrode ensemble (GNEE) to detect epithelial ovarian cancer antigen-125 (CA 125), a protein biomarker associated with ovarian cancer. CA 125 is the standard tumor marker used to follow women during or after treatment for epithelial ovarian cancer. The template protein CA 125 was initially incorporated into the thin-film coating and, upon extraction of protein from the accessible surfaces on the thin film, imprints for CA 125 were formed. The fabrication and analysis of the CA 125 imprinted GNEE was done by using cyclic voltammetry (CV), differential pulse voltammetry (DPV) and electrochemical impedance spectroscopy (EIS) techniques. The surfaces of the very thin, protein imprinted sites on GNEE are utilized for immunospecific capture of CA 125 molecules, and the mass of bound on the electrode surface can be detected as a reduction in the faradic current from the redox marker. Under optimal conditions, the developed sensor showed good increments at the studied concentration range of 0.5–400 U mL−1. The lowest detection limit was found to be 0.5 U mL−1. Spiked human blood serum and unknown real serum samples were analyzed. The presence of non-specific proteins in the serum did not significantly affect the sensitivity of our assay. Molecular imprinting using synthetic polymers and nanomaterials provides an alternative approach to the trace detection of biomarker proteins.

Differential expression of CDC25 phosphatases splice variants in human breast cancer cells

Background: CDC25 phosphatases control cell cycle progression by activating cyclin dependent kinases. The three CDC25 isoforms encoding genes are submitted to alternative splicing events which generate at least two variants for CDC25A and five for both CDC25B and CDC25C. An over-expression of CDC25 was reported in several types of cancer, including breast cancer, and is often associated with a poor prognosis. Nevertheless, most of the previous studies did not address the expression of CDC25 splice variants. Here, we evaluated CDC25 spliced transcripts expression in anti-cancerous drug-sensitive and resistant breast cancer cell lines in order to identify potential breast cancer biomarkers. Methods: CDC25 splice variants mRNA levels were evaluated by semi-quantitative RT-PCR and by an original real-time RT-PCR assay. Results: CDC25 spliced transcripts are differentially expres-sed in the breast cancer cell lines studied. An up-regulation of CDC25A2 variant and an increase of the CDC25C5/C1 ratio are associated to the multidrug-resistance in VCREMS and DOXOR breast cancer cells, compared to their sensitive counterpart cell line MCF-7. Additionally, CDC25B2 tran-script is exclusively over-expressed in VCREMS resistant cells and could therefore be involved in the development of certain type of drug resistance. Conclusions: CDC25 splice variants could represent interesting potential breast cancer prognostic biomarkers.

Differential expression of GSPT1 GGCn alleles in cancer

The human eukaryotic release factor 3a (eRF3a), encoded by the G1 to S phase transition 1 gene (GSPT1; alias eRF3a), is upregulated in various human cancers. GSPT1 contains a GGCn polymorphism in exon 1, encoding a polyglycine expansion in the N-terminal of the protein. The longer allele, GGC12, was previously shown to be associated to cancer. The GGC12 allele was present in 2.2% of colorectal cancer patients but was absent in Crohn disease patients and in the control group. Real-time quantitative RT-PCR analysis showed that the GGC12 allele was present at up to 10-fold higher transcription levels than the GGC10 allele (P < 0.001). No GSPT1 amplifications were detected, and there was no correlation between the length of the alleles and methylation levels of the CpG sites inside the GGC expansion. Using flow cytometry, we compared the levels of apoptosis and proliferation rates between cell lines with different genotypes, but detected no significant differences. Finally, we used a cytokinesis-block micronucleus assay to evaluate the frequency of micronuclei in the same cell lines. Cell lines with the longer alleles had higher frequencies of micronuclei in binucleated cells, which is probably a result of defects in mitotic spindle formation. Altogether, these findings indicate that GSPT1 should be considered a potential proto-oncogene.

A class of singular first order differential equations with applications in reaction-diffusion

Agências Financiadoras: FCT e MIUR

Cholesterol 24S-Hydroxylase overexpression inhibits the liver X receptor (LXR) pathway by activating small guanosine triphosphate-binding proteins (sGTPases) in neuronal cells

The neuronal-specific cholesterol 24S-hydroxylase (CYP46A1) is important for brain cholesterol elimination. Cyp46a1 null mice exhibit severe deficiencies in learning and hippocampal long-term potentiation, suggested to be caused by a decrease in isoprenoid intermediates of the mevalonate pathway. Conversely, transgenic mice overexpressing CYP46A1 show an improved cognitive function. These results raised the question of whether CYP46A1 expression can modulate the activity of proteins that are crucial for neuronal function, namely of isoprenylated small guanosine triphosphate-binding proteins (sGTPases). Our results show that CYP46A1 overexpression in SH-SY5Y neuroblastoma cells and in primary cultures of rat cortical neurons leads to an increase in 3-hydroxy-3-methyl-glutaryl-CoA reductase activity and to an overall increase in membrane levels of RhoA, Rac1, Cdc42 and Rab8. This increase is accompanied by a specific increase in RhoA activation. Interestingly, treatment with lovastatin or a geranylgeranyltransferase-I inhibitor abolished the CYP46A1 effect. The CYP46A1-mediated increase in sGTPases membrane abundance was confirmed in vivo, in membrane fractions obtained from transgenic mice overexpressing this enzyme. Moreover, CYP46A1 overexpression leads to a decrease in the liver X receptor (LXR) transcriptional activity and in the mRNA levels of ATP-binding cassette transporter 1, sub-family A, member 1 and apolipoprotein E. This effect was abolished by inhibition of prenylation or by co-transfection of a RhoA dominant-negative mutant. Our results suggest a novel regulatory axis in neurons; under conditions of membrane cholesterol reduction by increased CYP46A1 expression, neurons increase isoprenoid synthesis and sGTPase prenylation. This leads to a reduction in LXR activity, and consequently to a decrease in the expression of LXR target genes.

Controllability results for impulsive mixed type functional integro-differential evolution equations with nonlocal conditions

In this paper, we establish the controllability for a class of abstract impulsive mixed-type functional integro-differential equations with finite delay in a Banach space. Some sufficient conditions for controllability are obtained by using the Mönch fixed point theorem via measures of noncompactness and semigroup theory. Particularly, we do not assume the compactness of the evolution system. An example is given to illustrate the effectiveness of our results.