User guide for shock and blast simulation with the OctVCE code (version 3.5+), Department of Mechanical Engineering Report 2007/13

OctVCE is a cartesian cell CFD code produced especially for numerical simulations of shock and blast wave interactions with complex geometries. Virtual Cell Embedding (VCE) was chosen as its cartesian cell kernel as it is simple to code and sufficient for practical engineering design problems. This also makes the code much more ‘user-friendly’ than structured grid approaches as the gridding process is done automatically. The CFD methodology relies on a finite-volume formulation of the unsteady Euler equations and is solved using a standard explicit Godonov (MUSCL) scheme. Both octree-based adaptive mesh refinement and shared-memory parallel processing capability have also been incorporated. For further details on the theory behind the code, see the companion report 2007/12.

Improved methods for predicting individual leaf area and leaf senescence in maize (Zea mays)

The ability to predict leaf area and leaf area index is crucial in crop simulation models that predict crop growth and yield. Previous studies have shown existing methods of predicting leaf area to be inadequate when applied to a broad range of cultivars with different numbers of leaves. The objectives of the study were to (i) develop generalised methods of modelling individual and total plant leaf area, and leaf senescence, that do not require constants that are specific to environments and/or genotypes, (ii) re-examine the base, optimum, and maximum temperatures for calculation of thermal time for leaf senescence, and (iii) assess the method of calculation of individual leaf area from leaf length and leaf width in experimental work. Five cultivars of maize differing widely in maturity and adaptation were planted in October 1994 in south-eastern Queensland, and grown under non-limiting conditions of water and plant nutrient supplies. Additional data for maize plants with low total leaf number (12-17) grown at Katumani Research Centre, Kenya, were included to extend the range in the total leaf number per plant. The equation for the modified (slightly skewed) bell curve could be generalised for modelling individual leaf area, as all coefficients in it were related to total leaf number. Use of coefficients for individual genotypes can be avoided, and individual and total plant leaf area can be calculated from total leaf number. A single, logistic equation, relying on maximum plant leaf area and thermal time from emergence, was developed to predict leaf senescence. The base, optimum, and maximum temperatures for calculation of thermal time for leaf senescence were 8, 34, and 40 degrees C, and apply for the whole crop-cycle when used in modelling of leaf senescence. Thus, the modelling of leaf production and senescence is simplified, improved, and generalised. Consequently, the modelling of leaf area index (LAI) and variables that rely on LAI will be improved. For experimental purposes, we found that the calculation of leaf area from leaf length and leaf width remains appropriate, though the relationship differed slightly from previously published equations.

Structural analysis and molecular model of a self-incompatibility RNase from wild tomato

Self-incompatibility RNases (S-RNases) are an allelic series of style glycoproteins associated with rejection of self-pollen in solanaceous plants. The nucleotide sequences of S-RNase alleles from several genera have been determined, but the structure of the gene products has only been described for those from Nicotiana alata. We report on the N-glycan structures and the disulfide bonding of the S-3-RNase from wild tomato (Lycopersicon peruvianum) and use this and other information to construct a model of this molecule. The S-3-RNase has a single N-glycosylation site (Asn-28) to which one of three N-glycans is attached. S-3-RNase has seven Cys residues; six are involved in disulfide linkages (Cys-16-Cys-21, Cys-46-Cys-91, and Cys-166-Cys-177), and one has a free thiol group (Cys-150). The disulfide-bonding pattern is consistent with that observed in RNase Rh, a related RNase for which radiographic-crystallographic information is available. A molecular model of the S-3-RNase shows that four of the most variable regions of the S-RNases are clustered on one surface of the molecule. This is discussed in the context of recent experiments that set out to determine the regions of the S-RNase important for recognition during the self-incompatibility response.

Calcium is essential for the structural integrity of the cysteine-rich, ligand-binding repeat of the low-density lipoprotein receptor

Seven cysteine-rich repeats form the ligand-binding region of the low-density lipoprotein (LDL) receptor. Each of these repeats is assumed to bind a calcium ion, which is needed for association of the receptor with its ligands, LDL and beta-VLDL. The effects of metal ions on the folding of the reduced N-terminal cysteine-rich repeat have been examined by using reverse-phase high-performance liquid chromatography to follow the formation of fully oxidized isomers with different disulfide connectivities. in the absence of calcium many of the 15 possible isomers formed on oxidation, whereas in its presence the predominant product at equilibrium had the native disulfide bond connectivities. Other metals were far less effective at directing disulfide bond formation: Mn2+ partly mimicked the action of Ca2+, but Ba2+, Sr2+, and Mg2+ had little effect. This metal-ion specificity was also observed in two-dimensional H-1 NMR spectral studies: only Ca2+ induced the native three-dimensional fold. The two paramagnetic ions, Gd3+ and Mn2+, and Cd2+ did not promote adoption of a well-defined structure, and the two paramagnetic ions did not displace calcium ions. The location of calcium ion binding sites in the repeat was also explored by NMR spectroscopy. The absence of chemical shift changes for the side chain proton resonances of Asp26, Asp36, and Glu37 from pH 3.9 to 6.8 in the presence of calcium ions and their proximal location in the NMR structures implicated these side chains as calcium ligands. Deuterium exchange NMR experiments also revealed a network of hydrogen bonds that stabilizes the putative calcium-binding loop.

Structural and functional characterisation of human sulfotransferases

The human aryl sulfotransferases HAST4 and HAST4v vary by only two amino acids but exhibit markedly different affinity towards the sulfonate acceptor p-nitrophenol and the sulfonate donor 3'-phosphoadenosine-5'-phosphosulfate (PAPS). To determine the importance of each of these amino acid differences, chimeric constructs were made of HAST4 and HAST4v. By attaching the last 120 amino acids of HAST-4v to HAST4 (changing Thr235 to Asn235) we have been able to produce a protein that has a K-m for PAPS similar to HAST4v. The reverse construct, HAST4v/4 produces a protein with a K-m for PAPS similar to HAST4. These data suggests that the COOH-terminal of sulfotransferases is involved in co-factor binding. (C) 1998 Elsevier Science Ireland Ltd. All rights reserved.

Functional and structural characterization of isolated perfused stingray liver including effects of ischaemia/reperfusion

The morphological and functional characteristics of stingray liver were studied, including the effect of ischaemia/reperfusion. With an isolated perfused model, it was shown that the stingray liver was more resistant than the rat liver to ischaemia/reperfusion injury; this was consistent with the differing partial oxygen tensions usually present in the two species. This study confirmed that whereas stingray hepatocytes form tubules with central bile canaliculi as in other fish, the stingray liver has portal triads and a lobular architecture as in mammals. Apoptosis of hepatocytes, demonstrated in the normal liver, was only marginally enhanced by ischaemia/reperfusion. Resulting apoptotic bodies were phagocytized by macrophage-like cells in hepatocyte tubules. In contrast to rat liver, the stingray liver showed no necrosis after ischaemia-reperfusion. (C) 1998 W.B. Saunders Company Limited.

A study of the stability of subcycling algorithms in structural dynamics

Algorithms for explicit integration of structural dynamics problems with multiple time steps (subcycling) are investigated. Only one such algorithm, due to Smolinski and Sleith has proved to be stable in a classical sense. A simplified version of this algorithm that retains its stability is presented. However, as with the original version, it can be shown to sacrifice accuracy to achieve stability. Another algorithm in use is shown to be only statistically stable, in that a probability of stability can be assigned if appropriate time step limits are observed. This probability improves rapidly with the number of degrees of freedom in a finite element model. The stability problems are shown to be a property of the central difference method itself, which is modified to give the subcycling algorithm. A related problem is shown to arise when a constraint equation in time is introduced into a time-continuous space-time finite element model. (C) 1998 Elsevier Science S.A.

On the graded quantum Yang-Baxter and reflection equations

We clarify the extra signs appearing in the graded quantum Yang-Baxter reflection equations, when they are written in a matrix form. We find the boundary K-matrix for the Perk-Schultz six-vertex model, thus give a general solution to the graded reflection equation associated with it.

Structure determination of the three disulfide bond isomers of alpha-conotoxin GI: A model for the role of disulfide bonds in structural stability

The three possible disulfide bonded isomers of alpha-conotoxin GI have been selectively synthesised and their structures determined by H-1 NMR spectroscopy. alpha-Conotoxin GI derives from the venom of Conus geographus and is a useful neuropharmacological tool as it selectively binds to the nicotinic acetylcholine receptor (nAChR), a ligand-gated ion channel involved in nerve signal transmission. The peptide has the sequence ECCNPACGRHYSC-NH2, and the three disulfide bonded isomers are referred to as GI(2-7;3-13), GI(2-13;3-7) and GI(2-3;7-13). The NMR structure for the native isomer GI(2-7;3-13) is of excellent quality, with a backbone pairwise RMSD of 0.16 Angstrom for a family of 35 structures, and comprises primarily a distorted 3(10),, helix between residues 5 to 11. The two non-native isomers exhibit multiple conformers in solution, with the major populated forms being different in structure both from each other and from the native form. Structure-activity relationships for the native GI(2-7;3-13) as well as the role of the disulfide bonds on folding and stability of the three isomers are examined. It is concluded that the disulfide bonds in alpha-conotoxin GI play a crucial part in determining both the structure and stability of the peptide. A trend for increased conformational heterogeneity was observed in the order of GI(2-7;3-13) < GI(2-13;3-7) < GI(2-3;7-13). It was found that the peptide bond joining Cys2 to Cys3 in GI(2-3;7-13) is predominantly trans, rather than cis as theoretically predicted. These structural data are used to interpret the varying nAChR binding of the non-native forms. A model for the binding of native GI(2-7;3-13) to the mammalian nAChR is proposed, with an alpha-subunit binding face made up of Cys2, Asn4, Pro5, Ala6 and Cys7 and a selectivity face, comprised of Arg9 and His10. These two faces orient the molecule between the alpha and delta subunits of the receptor. The structure of the CCNPAC sequence of the native GI(2-7;3-13) is compared to the structure of the identical sequence from the toxic domain of heat-stable enterotoxins, which forms part of the receptor binding region of the enterotoxins, but which has a different disulfide connectivity. (C) 1998 Academic Press Limited.

Stage-specific proteophosphoglycan from Leishmania mexicana amastigotes - Structural characterization of novel mono-, di-, and triphosphorylated phosphodiester-linked oligosaccharides

Intracellular amastigotes of the protozoan parasite Leishmania mexicana secrete a macromolecular proteophosphoglycan (aPPG) into the phagolysosome of their host cell, the mammalian macrophage. The structures of aPPG glycans were analyzed by a combination of high pH anion exchange high pressure liquid chromatography, gas chromatography-mass spectrometry, enzymatic digestions, electrospray-mass spectrometry as well as H-1 and P-31 NMR spectroscopy. Some glycans are identical to oligosaccharides known from Leishmania mexicana promastigote lipophosphoglycan and secreted acid phosphatase, However, the majority of the aPPG glycans represent amastigote stage-specific and novel structures. These include neutral glycans ([Glc beta(1-3)](1-2)Gal beta 1-4Man, Gal beta 1-3Gal beta 1-4Man, Gal beta 1-3Glc beta 1-3Gal beta 1-4Man), several monophosphorylated glycans containing the conserved phosphodisaccharide backbone (R-3-[PO4-6-Gal]beta 1-4Man) but carrying stage-specific modifications (R = Gal beta 1-, [Glc beta 1-3](1-2)Glc beta 1-), and monophosphorylated aPPG tri- and tetrasaccharides that are uniquely phosphorylated on the terminal hexose (PO4-6-Glc beta 1-3Gal beta 1-4Man, PO4-6-Glc beta 1-3Glc beta 1-3Gal beta 1-4Man, PO4-6-Gal beta 1-3Glc beta 1-3Gal beta 1-4Man), In addition aPPG contains highly unusual di- and triphosphorylated glycans whose major species are PO4-6-Glc beta 1-3Glc beta 1-3[PO4-6-Gal]beta 1-4Man, PO4-6-Gal beta 1-3Glc beta 1-3 [PO4-6-Gal]beta 1-4Man, PO4-6-GaL beta 1-3Glc beta 1-3Glc beta 1-3[PO4-6-Gal]beta 1-4Man, PO4-6-Glc beta 1-3[PO4-6-Glc]beta 1-3[PO4-6-Gal]beta 1-4Man, PO4-6Gal beta 1-3[PO4-6-Glc]beta 1-3Glc beta 1-3[PO4-6-Gal]beta 1-4Man, and PO4-6-Glc beta 1-3[PO4-6-Glc]beta 1-3Glc beta 1-3[PO4-6-Gal]beta 1-4Man. These glycans are linked together by the conserved phosphodiester R-Man alpha 1-PO4-6-Gal-R or the novel phosphodiester R-Man alpha 1-PO4-6-Glc-R and are connected to Ser(P) of the protein backbone most likely via the linkage R-Man alpha 1-PO4-Ser. The variety of stage-specific glycan structures in Leishmania mexicana aPPG suggests the presence of developmentally regulated amastigote glycosyltransferases which may be potential anti-parasite drug targets.

Modelling fungal growth on surfaces

Fungal growth in time and space at the substrate surface was modelled for a simple system mimicking solid-state fermentation, using a polycarbonate Nucleopore membrane laid over a glucose solution. Biomass production depends on both tip density and the diffusion of glucose within the fungal hyphae. The model predicts early increases in both height and concentration, followed by a period in which the biomass profile moves with a constant wavefront. The rate of increase in height increases as tip diffusivity increases or as the Monod saturation constant for glucose decreases.

Subcycling first- and second-order generalizations of the trapezoidal rule

Subcycling algorithms which employ multiple timesteps have been previously proposed for explicit direct integration of first- and second-order systems of equations arising in finite element analysis, as well as for integration using explicit/implicit partitions of a model. The author has recently extended this work to implicit/implicit multi-timestep partitions of both first- and second-order systems. In this paper, improved algorithms for multi-timestep implicit integration are introduced, that overcome some weaknesses of those proposed previously. In particular, in the second-order case, improved stability is obtained. Some of the energy conservation properties of the Newmark family of algorithms are shown to be preserved in the new multi-timestep extensions of the Newmark method. In the first-order case, the generalized trapezoidal rule is extended to multiple timesteps, in a simple way that permits an implicit/implicit partition. Explicit special cases of the present algorithms exist. These are compared to algorithms proposed previously. (C) 1998 John Wiley & Sons, Ltd.

New integrable boundary conditions for the q-deformed supersymmetric U model and Bethe ansatz equations

New classes of integrable boundary conditions for the q-deformed (or two-parameter) supersymmetric U model are presented. The boundary systems are solved by using the coordinate space Bethe ansatz technique and Bethe ansatz equations are derived. (C) 1998 Elsevier Science B.V.