A partial velocity approach to subcycling structural dynamics

Subcycling, or the use of different timesteps at different nodes, can be an effective way of improving the computational efficiency of explicit transient dynamic structural solutions. The method that has been most widely adopted uses a nodal partition. extending the central difference method, in which small timestep updates are performed interpolating on the displacement at neighbouring large timestep nodes. This approach leads to narrow bands of unstable timesteps or statistical stability. It also can be in error due to lack of momentum conservation on the timestep interface. The author has previously proposed energy conserving algorithms that avoid the first problem of statistical stability. However, these sacrifice accuracy to achieve stability. An approach to conserve momentum on an element interface by adding partial velocities is considered here. Applied to extend the central difference method. this approach is simple. and has accuracy advantages. The method can be programmed by summing impulses of internal forces, evaluated using local element timesteps, in order to predict a velocity change at a node. However, it is still only statistically stable, so an adaptive timestep size is needed to monitor accuracy and to be adjusted if necessary. By replacing the central difference method with the explicit generalized alpha method. it is possible to gain stability by dissipating the high frequency response that leads to stability problems. However. coding the algorithm is less elegant, as the response depends on previous partial accelerations. Extension to implicit integration, is shown to be impractical due to the neglect of remote effects of internal forces acting across a timestep interface. (C) 2002 Elsevier Science B.V. All rights reserved.

Analysis of the mouse transcriptome for genes involved in the function of the nervous system

We analyzed the mouse Representative Transcript and Protein Set for molecules involved in brain function. We found full-length cDNAs of many known brain genes and discovered new members of known brain gene families, including Family 3 G-protein coupled receptors, voltage-gated channels, and connexins. We also identified previously unknown candidates for secreted neuroactive molecules. The existence of a large number of unique brain ESTs suggests an additional molecular complexity that remains to be explored. A list of genes containing CAG stretches in the coding region represents a first step in the potential identification of candidates for hereditary neurological disorders.

Screening of human primary melanocytes of defined melanocortin-1 receptor genotype: Pigmentation marker, ultrastructural and UV-survival studies

Recent population studies have demonstrated an association with the red-hair and fair-skin phenotype with variant alleles of the melanocortin-1 receptor (MC1R) which result in amino acid substitutions within the coding region leading to an altered receptor activity. In particular, Arg151Cys, Arg160Trp and Asp294His were the most commonly associated variants seen in the south-east Queensland population with at least one of these alleles found in 93% of those with red hair. In order to study the individual effects of these variants on melanocyte biology and melanocytic pigmentation, we established a series of human melanocyte strains genotyped for the MC1R receptor which included wild-type consensus, variant heterozygotes, compound heterozygotes and homozygotes for Arg151Cys, Arg160Trp, Val60Leu and Val92Met alleles. These strains ranged from darkly pigmented to amelanotic, with all strains of consensus sequence having dark pigmentation. UV sensitivity was found not to be associated with either MC1R genotype or the level of pigmentation with a range of sensitivities seen across all genotypes. Ultrastructural analysis demonstrated that while consensus strains contained stage IV melanosomes in their terminal dendrites, Arg151Cys and Arg160Trp homozygote strains contained only stage II melanosomes. This was despite being able to show expression of tyrosinase and tyrosinase-related protein-1 markers, although at reduced levels and an ability to convert exogenous 3,4-dihydroxyphenyl-alanine (DOPA) to melanin in these strains.

Rates of gene rearrangement and nucleotide substitution are correlated in the mitochondrial genomes of insects

A number of studies indicated that lineages of animals with high rates of mitochondrial (mt) gene rearrangement might have high rates of mt nucleotide substitution. We chose the hemipteroid assemblage and the Insecta to test the idea that rates of mt gene rearrangement and mt nucleotide substitution are correlated. For this purpose, we sequenced the mt genome of a lepidopsocid from the Psocoptera, the only order of hemipteroid insects for which an entire mtDNA sequence is not available. The mt genome of this lepidopsocid is circular, 16,924 bp long, and contains 37 genes and a putative control region; seven tRNA genes and a protein-coding gene in this genome have changed positions relative to the ancestral arrangement of mt genes of insects. We then compared the relative rates of nucleotide substitution among species from each of the four orders of hemipteroid insects and among the 20 insects whose mt genomes have been sequenced entirely. All comparisons among the hernipteroid insects showed that species with higher rates of gene rearrangement also had significantly higher rates of nucleotide substitution statistically than did species with lower rates of gene rearrangement. In comparisons among the 20 insects, where the mt genomes of the two species differed by more than five breakpoints, the more rearranged species always had a significantly higher rate of nucleotide substitution than the less rearranged species. However, in comparisons where the mt genomes of two species differed by five or less breakpoints, the more rearranged species did not always have a significantly higher rate of nucleotide substitution than the less rearranged species. We tested the statistical significance of the correlation between the rates of mt gene rearrangement and mt nucleotide substitution with nine pairs of insects that were phylogenetically independent from one 2 another. We found that the correlation was positive and statistically significant (R-2 = 0.73, P = 0.01; R-s = 0.67, P < 0.05). We propose that increased rates of nucleotide substitution may lead to increased rates of gene rearrangement in the mt genomes of insects.

Effects of additional sequences directly downstream from the AUG on the expression of GFP gene

We have studied the expression of the green fluorescent protein (GFP) gene to gain more understanding of the effects of additional nucleotide triplets (codons) downstream from the initiation codon on the translation of the GFP mRNA in CHO and Cos1 cells. A leader sequence of six consecutive identical codons (GUG, CUC, AGU or UCA) was introduced into a humanized GFP (hm gfp) gene downstream from the AUG to produce four GFP gene variants. Northern blot and RT-PCR analysis indicated that mRNA transcription from the GFP gene was not significantly affected by any of the additional sequences. However, immunoblotting and FACS analysis revealed that AGU and UCA GFP variants produced GFP at a mean level per cell 3.5-fold higher than the other two GFP variants and the hm gfp gene. [35S]-Methionine labeling and immunoprecipitation demonstrate that GFP synthesis was very active in UCA variant transfected-cells, but not in GUG variant and hm gfp transfected-cells. Moreover, proteasome inhibitor MG-132 treatment indicated that the GFPs encoded by each of the GFP variants and the hm gfp were equally stable, and this together with the comparable mRNA levels observed for each construct suggested that the different steady-state GFP concentrations observed reflected different translation efficiencies of the various GFP genes. In addition, the CUC GFP variant, when transiently transfected into CHO or COS-1 cells, did not produce any GFP expressing cells (fully green cells), and the GUG variant produced GFP expressing cells less than 10%, while AGU and UCA GFP variants up to 30–35% in a time course study from 8 to 36 h posttransfection. Analysis of the potential secondary structure of the GFP variant mRNAs especially in the translation initiation region suggested that the secondary structure of the GFP mRNAs was unlikely to explain the different translation efficiencies of the GFP variants. The present findings indicate that a change of the initiation context of the GFP gene by addition of extra coding sequence can alter the translation efficiency of GFP mRNA, providing a means of more efficient expression of GFP in eukaryotic cells.

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.