Sausage: protein threading with flexible force fields

Sausage is a protein sequence threading program, but with remarkable run-time flexibility. Using different scripts, it can calculate protein sequence-structure alignments, search structure libraries, swap force fields, create models form alignments, convert file formats and analyse results. There are several different force fields which might be classed as knowledge-based, although they do not rely on Boltzmann statistics. Different force fields are used for alignment calculations and subsequent ranking of calculated models.

Identification of novel non-autonomous CemaT transposable elements and evidence of their mobility within the C. elegans genome

We describe here two new transposable elements, CemaT4 and CemaT5, that were identified within the sequenced genome of Caenorhabditis elegans using homology based searches. Five variants of CemaT4 were found, all non-autonomous and sharing 26 bp inverted terminal repeats (ITRs) and segments (152-367 bp) of sequence with similarity to the CemaT1 transposon of C. elegans. Sixteen copies of a short, 30 bp repetitive sequence, comprised entirely of an inverted repeat of the first 15 bp of CemaT4's ITR, were also found, each flanked by TA dinucleotide duplications, which are hallmarks of target site duplications of mariner-Tc transposon transpositions. The CemaT5 transposable element had no similarity to maT elements, except for sharing identical ITR sequences with CemaT3. We provide evidence that CemaT5 and CemaT3 are capable of excising from the C. elegans genome, despite neither transposon being capable of encoding a functional transposase enzyme. Presumably, these two transposons are cross-mobilised by an autonomous transposon that recognises their shared ITRs. The excisions of these and other non-autonomous elements may provide opportunities for abortive gap repair to create internal deletions and/or insert novel sequence within these transposons. The influence of non-autonomous element mobility and structural diversity on genome variation is discussed.

Casimir invariants and characteristic identities for gl(infinity)

A full set of (higher-order) Casimir invariants for the Lie algebra gl(infinity) is constructed and shown to be well defined in the category O-FS generated by the highest weight (unitarizable) irreducible representations with only a finite number of nonzero weight components. Moreover, the eigenvalues of these Casimir invariants are determined explicitly in terms of the highest weight. Characteristic identities satisfied by certain (infinite) matrices with entries from gl(infinity) are also determined and generalize those previously obtained for gl(n) by Bracken and Green [A. J. Bracken and H. S. Green, J. Math. Phys. 12, 2099 (1971); H. S. Green, ibid. 12, 2106 (1971)]. (C) 1997 American Institute of Physics.

Three-dimensional orthotropic viscoelastic finite element model of a human ligament

Ligaments undergo finite strain displaying hyperelastic behaviour as the initially tangled fibrils present straighten out, combined with viscoelastic behaviour (strain rate sensitivity). In the present study the anterior cruciate ligament of the human knee joint is modelled in three dimensions to gain an understanding of the stress distribution over the ligament due to motion imposed on the ends, determined from experimental studies. A three dimensional, finite strain material model of ligaments has recently been proposed by Pioletti in Ref. [2]. It is attractive as it separates out elastic stress from that due to the present strain rate and that due to the past history of deformation. However, it treats the ligament as isotropic and incompressible. While the second assumption is reasonable, the first is clearly untrue. In the present study an alternative model of the elastic behaviour due to Bonet and Burton (Ref. [4]) is generalized. Bonet and Burton consider finite strain with constant modulii for the fibres and for the matrix of a transversely isotropic composite. In the present work, the fibre modulus is first made to increase exponentially from zero with an invariant that provides a measure of the stretch in the fibre direction. At 12% strain in the fibre direction, a new reference state is then adopted, after which the material modulus is made constant, as in Bonet and Burton's model. The strain rate dependence can be added, either using Pioletti's isotropic approximation, or by making the effect depend on the strain rate in the fibre direction only. A solid model of a ligament is constructed, based on experimentally measured sections, and the deformation predicted using explicit integration in time. This approach simplifies the coding of the material model, but has a limitation due to the detrimental effect on stability of integration of the substantial damping implied by the nonlinear dependence of stress on strain rate. At present, an artificially high density is being used to provide stability, while the dynamics are being removed from the solution using artificial viscosity. The result is a quasi-static solution incorporating the effect of strain rate. Alternate approaches to material modelling and integration are discussed, that may result in a better model.

Lanczos subspace filter diagonalization: Homogeneous recursive filtering and a low-storage method for the calculation of matrix elements

We develop a new iterative filter diagonalization (FD) scheme based on Lanczos subspaces and demonstrate its application to the calculation of bound-state and resonance eigenvalues. The new scheme combines the Lanczos three-term vector recursion for the generation of a tridiagonal representation of the Hamiltonian with a three-term scalar recursion to generate filtered states within the Lanczos representation. Eigenstates in the energy windows of interest can then be obtained by solving a small generalized eigenvalue problem in the subspace spanned by the filtered states. The scalar filtering recursion is based on the homogeneous eigenvalue equation of the tridiagonal representation of the Hamiltonian, and is simpler and more efficient than our previous quasi-minimum-residual filter diagonalization (QMRFD) scheme (H. G. Yu and S. C. Smith, Chem. Phys. Lett., 1998, 283, 69), which was based on solving for the action of the Green operator via an inhomogeneous equation. A low-storage method for the construction of Hamiltonian and overlap matrix elements in the filtered-basis representation is devised, in which contributions to the matrix elements are computed simultaneously as the recursion proceeds, allowing coefficients of the filtered states to be discarded once their contribution has been evaluated. Application to the HO2 system shows that the new scheme is highly efficient and can generate eigenvalues with the same numerical accuracy as the basic Lanczos algorithm.

Wildlife-based tourism & increased support for nature conservation financially and otherwise: evidence from sea-turtle ecotourism at Mon Repos

The arguments of most conservationists supporting ecotourism have been based on the view that it is environmentally friendly as a resource-use and that receipts from it can counter demands to use the natural resources involved for more extractive economic purposes. But wildlife-based ecotourism can also have positive impacts in itself on the willingness of tourists to pay for wildlife conservation, strengthen the pro-conservation attitudes of tourists, and foster personal actions by them that contribute to wildlife conservation. These aspects are explored in this article on the basis of a survey of tourists visiting Mon Repos Beach near Bundaberg, Queensland, for the purpose of watching marine turtles. The results enable several of the conservation impacts of this experience on tourists to be quantified, and highlight important relationships between specific socio-economic variables and the willingness of tourists to pay for the protection of sea turtles. Furthermore, it is shown that the on-site experiences of ecotourists have positive impacts on the willingness of tourists to pay for the conservation of wildlife, and that willingness to pay is sensitive to whether or not wildlife is seen. It is suggested that in situ ecotourism is likely to be a more powerful force for fostering pro-conservation attitudes and actions among visitors than ex situ wildlife-based tourism in aquaria and zoos.

Repeatability of maximal voluntary force and of surface EMG variables during voluntary isometric contraction of quadriceps muscles in healthy subjects

The repeatability of initial values and rate of change of EMG signal mean spectral frequency (MNF), average rectified values (ARV), muscle fiber conduction velocity (CV) and maximal voluntary contraction (MVC) was investigated in the vastus medialis obliquus (VMO) and vastus lateralis (VL) muscles of both legs of nine healthy male subjects during voluntary, isometric contractions sustained for 50 s at 50% MVC. The values of MVC were recorded for both legs three times on each day and for three subsequent days, while the EMG signals have been recorded twice a day for three subsequent days. The degree of repeatability was investigated using the Fisher test based upon the ANalysis Of VAriance (ANOVA), the Standard Error of the Mean (SEM) and the Intraclass Correlation Coefficient (ICC). Data collected showed a high level of repeatability of MVC measurement (normalized SEM from 1.1% to 6.4% of the mean). MNF and ARV initial values also showed a high level of repeatability (ICC > 70% for all muscles and legs except right VMO). At 50% MVC level no relevant pattern of fatigue was observed for the VMO and VL muscles, suggesting that other portions of the quadriceps might have contributed to the generated effort. These observations seem to suggest that in the investigation of muscles belonging to a multi-muscular group at submaximal level, the more selective electrically elicited contractions should be preferred to voluntary contractions. (C) 2001 Elsevier Science Ltd. All rights reserved.

Amino acid residues in conodont elements

Thermally unaltered conodont elements, brachiopods. and vertebrates were analyzed with reverse phase high profile liquid chromatography to locate and quantify amino acid remnants of the original organic matrix in the fossils. No consistent similarities in amino acid content were found in conodont taxa. and criteria based on organic residues appear to have no taxonomic significance in the fossils tested from these localities. However, hydroxyproline. an amino acid that is found in the collagen molecules of animals. as well as in the glycoproteins in the cell walls and reproductive tissues of certain plants, is represented in most taxa. The organic matter retained in the impermeable crowns of conodont elements might have been derived originally from a form of collagen. Biochemical analyses. correlated with histochemical tests, demonstrate that organic matter is an integral part of the hyaline tissue of the element crown and not the result of surface contamination. Tests of a range of vertebrate and invertebrate fossil hard tissues produced similar results. The analyses indicate that hyaline tissue in the conodont element crown is not a form of vertebrate enamel. which contains no collagen. Albid tissue. with little or no organic content. is not a form of vertebrate bone or dentine, both based on collagen and low in mineral. Although these results do not help to determine the phylogenetic affinities of conodont animals, they indicate teat conodont elements do not contain hard tissues characteristic of vertebrate animals.

Whither our Art? Clinical Wisdom & Evidence-based Medicine

The relationship between evidence-based medicine (EBM) and clinical judgement is the subject of conceptual and practical dispute. For example, EBM and clinical guidelines are seen to increasingly dominate medical decision-making at the expense of other, human elements, and to threaten the art of medicine. Clinical wisdom always remains open to question. We want to know why particular beliefs are held, and the epistemological status of claims based in wisdom or experience. The paper critically appraises a number of claims and distinctions, and attempts to clarify the connections between EBM, clinical experience and judgement, and the objective and evaluative categories of medicine. I conclude that to demystify clinical wisdom is not to devalue it. EBM ought not be conceived as needing to be limited or balanced by clinical wisdom, since if its language is translatable into terms comprehensible and applicable to individuals, it helps constitute clinical wisdom. Failure to appreciate this constitutive relation will help perpetuate medical paternalism and delay the adoption of properly evidence-based practice, which would be both unethical and unwise.

Finite difference time domain (FDTD) method for modeling the effect of switched gradients on the human body in MRI

Numerical modeling of the eddy currents induced in the human body by the pulsed field gradients in MRI presents a difficult computational problem. It requires an efficient and accurate computational method for high spatial resolution analyses with a relatively low input frequency. In this article, a new technique is described which allows the finite difference time domain (FDTD) method to be efficiently applied over a very large frequency range, including low frequencies. This is not the case in conventional FDTD-based methods. A method of implementing streamline gradients in FDTD is presented, as well as comparative analyses which show that the correct source injection in the FDTD simulation plays a crucial rule in obtaining accurate solutions. In particular, making use of the derivative of the input source waveform is shown to provide distinct benefits in accuracy over direct source injection. In the method, no alterations to the properties of either the source or the transmission media are required. The method is essentially frequency independent and the source injection method has been verified against examples with analytical solutions. Results are presented showing the spatial distribution of gradient-induced electric fields and eddy currents in a complete body model.

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).

Analysis of the mouse transcriptome based on functional annotation of 60,770 full-length cDNAs

Only a small proportion of the mouse genome is transcribed into mature messenger RNA transcripts. There is an international collaborative effort to identify all full-length mRNA transcripts from the mouse, and to ensure that each is represented in a physical collection of clones. Here we report the manual annotation of 60,770 full-length mouse complementary DNA sequences. These are clustered into 33,409 'transcriptional units', contributing 90.1% of a newly established mouse transcriptome database. Of these transcriptional units, 4,258 are new protein-coding and 11,665 are new non-coding messages, indicating that non-coding RNA is a major component of the transcriptome. 41% of all transcriptional units showed evidence of alternative splicing. In protein-coding transcripts, 79% of splice variations altered the protein product. Whole-transcriptome analyses resulted in the identification of 2,431 sense-antisense pairs. The present work, completely supported by physical clones, provides the most comprehensive survey of a mammalian transcriptome so far, and is a valuable resource for functional genomics.

Knowledge-based system on optimum design of liquid retaining structures with genetic algorithms

This paper delineates the development of a prototype hybrid knowledge-based system for the optimum design of liquid retaining structures by coupling the blackboard architecture, an expert system shell VISUAL RULE STUDIO and genetic algorithm (GA). Through custom-built interactive graphical user interfaces under a user-friendly environment, the user is directed throughout the design process, which includes preliminary design, load specification, model generation, finite element analysis, code compliance checking, and member sizing optimization. For structural optimization, GA is applied to the minimum cost design of structural systems with discrete reinforced concrete sections. The design of a typical example of the liquid retaining structure is illustrated. The results demonstrate extraordinarily converging speed as near-optimal solutions are acquired after merely exploration of a small portion of the search space. This system can act as a consultant to assist novice designers in the design of liquid retaining structures.

A coupled knowledge-based expert system for design of liquid-retaining structures

This paper describes a coupled knowledge-based system (KBS) for the design of liquid-retaining structures, which can handle both the symbolic knowledge processing based on engineering heuristics in the preliminary synthesis stage and the extensive numerical crunching involved in the detailed analysis stage. The prototype system is developed by employing blackboard architecture and a commercial shell VISUAL RULE STUDIO. Its present scope covers design of three types of liquid-retaining structures, namely, a rectangular shape with one compartment, a rectangular shape with two compartments and a circular shape. Through custom-built interactive graphical user interfaces, the user is directed throughout the design process, which includes preliminary design, load specification, model generation, finite element analysis, code compliance checking and member sizing optimization. It is also integrated with various relational databases that provide the system with sectional properties, moment and shear coefficients and final member details. This system can act as a consultant to assist novice designers in the design of liquid-retaining structures with increase in efficiency and optimization of design output and automated record keeping. The design of a typical example of the liquid-retaining structure is also illustrated. (C) 2003 Elsevier B.V All rights reserved.