Predicting the columnar-to-equiaxed transition for a distribution of nucleation undercoolings

A deterministic mathematical model for steady-state unidirectional solidification is proposed to predict the columnar-to-equiaxed transition. In the model, which is an extension to the classic model proposed by Hunt [Hunt JD. Mater Sci Eng 1984;65:75], equiaxed grains nucleate according to either a normal or a log-normal distribution of nucleation undercoolings. Growth maps are constructed, indicating either columnar or equiaxed solidification as a function of the velocity of isotherms and temperature gradient. The fields A columnar and equiaxed growth change significantly with the spread of the nucleation undercooling distribution. Increasing the spread Favors columnar solidification if the dimensionless velocity of the isotherms is larger than 1. For a velocity less than 1, however, equiaxed solidification is initially favored, but columnar solidification is enhanced for a larger increase in the spread. This behavior was confirmed by a stochastic model, which showed that an increase in the distribution spread Could change the grain structure from completely columnar to 50% columnar grains. (c) 2008 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Formulating Mixed Models for Experiments, Including Longitudinal Experiments

Mixed models have become important in analyzing the results of experiments, particularly those that require more complicated models (e.g., those that involve longitudinal data). This article describes a method for deriving the terms in a mixed model. Our approach extends an earlier method by Brien and Bailey to explicitly identify terms for which autocorrelation and smooth trend arising from longitudinal observations need to be incorporated in the model. At the same time we retain the principle that the model used should include, at least, all the terms that are justified by the randomization. This is done by dividing the factors into sets, called tiers, based on the randomization and determining the crossing and nesting relationships between factors. The method is applied to formulate mixed models for a wide range of examples. We also describe the mixed model analysis of data from a three-phase experiment to investigate the effect of time of refinement on Eucalyptus pulp from four different sources. Cubic smoothing splines are used to describe differences in the trend over time and unstructured covariance matrices between times are found to be necessary.

Root Water Extraction under Combined Water and Osmotic Stress

Using a numerical implicit model for root water extraction by a single root in a symmetric radial flow problem, based on the Richards equation and the combined convection-dispersion equation, we investigated some aspects of the response of root water uptake to combined water and osmotic stress. The model implicitly incorporates the effect of simultaneous pressure head and osmotic head on root water uptake, and does not require additional assumptions (additive or multiplicative) to derive the combined effect of water and salt stress. Simulation results showed that relative transpiration equals relative matric flux potential, which is defined as the matric flux potential calculated with an osmotic pressure head-dependent lower bound of integration, divided by the matric flux potential at the onset of limiting hydraulic conditions. In the falling rate phase, the osmotic head near the root surface was shown to increase in time due to decreasing root water extraction rates, causing a more gradual decline of relative transpiration than with water stress alone. Results furthermore show that osmotic stress effects on uptake depend on pressure head or water content, allowing a refinement of the approach in which fixed reduction factors based on the electrical conductivity of the saturated soil solution extract are used. One of the consequences is that osmotic stress is predicted to occur in situations not predicted by the saturation extract analysis approach. It is also shown that this way of combining salinity and water as stressors yields results that are different from a purely multiplicative approach. An analytical steady state solution is presented to calculate the solute content at the root surface, and compared with the outputs of the numerical model. Using the analytical solution, a method has been developed to estimate relative transpiration as a function of system parameters, which are often already used in vadose zone models: potential transpiration rate, root length density, minimum root surface pressure head, and soil theta-h and K-h functions.

Structural complexes of human adenine phosphoribosyltransferase reveal novel features of the APRT catalytic mechanism

Adenine phosphoribosyltransferase (APRT) is an important enzyme component of the purine recycling pathway. Parasitic protozoa of the order Kinetoplastida are unable to synthesize purines de novo and use the salvage pathway for the synthesis of purine bases rendering this biosynthetic pathway an attractive target for antiparasitic drug design. The recombinant human adenine phosphoribosyltransferase (hAPRT) structure was resolved in the presence of AMP in the active site to 1.76 angstrom resolution and with the substrates PRPP and adenine simultaneously bound to the catalytic site to 1.83 angstrom resolution. An additional structure was solved containing one subunit of the dimer in the apo-form to 2.10 angstrom resolution. Comparisons of these three hAPRT structures with other `type I` PRTases revealed several important features of this class of enzymes. Our data indicate that the flexible loop structure adopts an open conformation before and after binding of both substrates adenine and PRPR Comparative analyses presented here provide structural evidence to propose the role of Glu 104 as the residue that abstracts the proton of adenine N9 atom before its nucleophilic attack on the PRPP anomeric carbon. This work leads to new insights to the understanding of the APRT catalytic mechanism.

Probabilistic Action System Trace Semantics

Action systems are a construct for reasoning about concurrent, reactive systems, in which concurrent behaviour is described by interleaving atomic actions. Sere and Troubitsyna have proposed an extension to action systems in which actions may be expressed and composed using discrete probabilistic choice as well as demonic nondeterministic choice. In this paper we develop a trace-based semantics for probabilistic action systems. This semantics provides a simple theoretical base on which practical refinement rules for probabilistic action systems may be justified.

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.

Structural basis of autoregulation of phenylalanine hydroxylase

Phenylalanine hydroxylase converts phenylalanine to tyrosine, a rate-limiting step in phenylalanine catabolism and protein and neurotransmitter biosynthesis. It is tightly regulated by the substrates phenylalanine and tetrahydrobiopterin and by phosphorylation. We present the crystal structures of dephosphorylated and phosphorylated forms of a dimeric enzyme with catalytic and regulatory properties of the wild-type protein. The structures reveal a catalytic domain flexibly linked to a regulatory domain. The latter consists of an N-terminal autoregulatory sequence (containing Ser 16, which is the site of phosphorylation) that extends over the active site pocket, and an alpha-beta sandwich core that is, unexpectedly, structurally related to both pterin dehydratase and the regulatory domains of metabolic enzymes. Phosphorylation has no major structural effects in the absence of phenylalanine, suggesting that phenylalanine and phosphorylation act in concert to activate the enzyme through a combination of intrasteric and possibly allosteric mechanisms.

Structure of a putative ancestral protein encoded by a single sequence repeat from a multidomain proteinase inhibitor gene from Nicotiana alata

Background: The ornamental tobacco Nicotiana alata produces a series of proteinase inhibitors (Pls) that are derived from a 43 kDa precursor protein, NaProPl. NaProPl contains six highly homologous repeats that fold to generate six separate structural domains, each corresponding to one of the native Pls. An unusual feature of NaProPl is that the structural domains lie across adjacent repeats and that the sixth Pl domain is generated from fragments of the first and sixth repeats. Although the homology of the repeats suggests that they may have arisen from gene duplication, the observed folding does not appear to support this. This study of the solution structure of a single NaProPl repeat (aPl1) forms a basis for unravelling the mechanism by which this protein may have evolved, Results: The three-dimensional structure of aPl1 closely resembles the triple-stranded antiparallel beta sheet observed in each of the native Pls. The five-residue sequence Glu-Glu-Lys-Lys-Asn, which forms the linker between the six structural domains in NaProPl, exists as a disordered loop in aPl1. The presence of this loop in aPl1 results in a loss of the characteristically flat and disc-like topography of the native inhibitors. Conclusions: A single repeat from NaProPl is capable of folding into a compact globular domain that displays native-like Pl activity. Consequently, it is possible that a similar single-domain inhibitor represents the ancestral protein from which NaProPl evolved.

Aluminium alloys - their physical and mechanical properties: Grain formation in AlSi7Mg0.35 foundry alloy at low superheat

Hypoeutectic AI-Si alloys represent the most widely used alloy system for cast aluminium applications. This system has a unique behaviour with respect to grain formation where an increase in silicon content results in a transition to larger grain sizes after a minimum at an intermediate concentration. As a result of the already large solute content, grain refinement by solute additions is inefficient and nucleant particles from the common aluminium grain refiners are not as effective as in wrought alloys. However, casting conditions, such as a low pouring temperature, that promote the formation of wall crystals tie. crystals nucleated in the thermally undercooled layer at or next to mould walls) are very effective in yielding a small grain size. This paper presents results of an investigation of the effect of low superheat and mould preheat temperature on grain size. It was found that pouring temperature controls the effectiveness of the wall mechanism while mould preheat has little effect until high preheat temperatures at which a large increase in grain size occurs. The observed changes in grain size are explained in terms of the balance between nucleation rate and survival rate of crystal nuclei resulting from changes in superheat and mould temperature.

Determination of the stability constants of mercury(II) complexes of mixed donor macrobicyclic encapsulating ligands

The reactions of mercury(II) with the mixed donor encapsulating ligands 3,6,16-trithia-6,11,19-triazabicyclo[6.6.6]icosane (AMN(3)S(3)sar) and 1-amino-8-methyl-6,19-dithia-3,10,13,16-tetraazabicyclo[6.6.6]icosane (AMN(4)S(2)sar) have been studied. NMR ligand-ligand competition experiments with the ligands 1,4,8,11-tetraazaeyclotetradecane ([14]aneN(4)), 1-thia-4,7,10-triazacyclododecane ([12]aneN(3)S) and ethylenediaminetetraacetic acid (EDTA) with AMN(3)S(3)sar and Hg(II) indicated that [14]aneN(4) would be an appropriate competing ligand for the, determination of the Hg(II) stability constant. Calculations indicated the ratio of concentrations of AMN3S3sar, [14]aneN(4) and Hg(II) required for the determination of the stability constant ranged from 1:1:1 to 1:5:1. Refinement of the titration curves yielded log(10)K[Hg(AMN(3)S(3)sar)](2+) = 17.7. A similar competition titration resulted in the determination of the stability constant for the AMN(4)S(2)sar system as log(10)K[Hg(AMN(4)S(2)sar)](2+) = 19.5. The observed binding constants for the mixed N/S donor systems and the hexaaza analogues sar (3,6,10,13,16,19-hexaazabicyclo [6.6.6]icosane) and diamsar (1,8-diamino-3,6,10,13,16,19 -hexazabicyclo [6.6.6] icosane (log(10)K-[Hg(diamsar)](2+) = 26.4; log(10)K[Hg(sar)](2+) = 28.1) differ by approximately ten orders of magnitude. The difference is ascribed not to a cryptate effect but to a mismatch in the Hg-N and Hg-S bond lengths in the N/S systems.

Structuring real-time Object-Z specifications

This paper presents a means of structuring specifications in real-time Object-Z: an integration of Object-Z with the timed refinement calculus. Incremental modification of classes using inheritance and composition of classes to form multi-component systems are examined. Two approaches to the latter are considered: using Object-Z's notion of object instantiation and introducing a parallel composition operator similar to those found in process algebras. The parallel composition operator approach is both more concise and allows more general modelling of concurrency. Its incorporation into the existing semantics of real-time Object-Z is presented.

Flexibility revealed by the 1.85 angstrom crystal structure of the beta sliding-clamp subunit of Escherichia coli DNA polymerase III

The beta subunit of the Escherichia coli replicative DNA polymerase III holoenzyme is the sliding clamp that interacts with the alpha (polymerase) subunit to maintain the high processivity of the enzyme. The beta protein is a ring-shaped dimer of 40.6 kDa subunits whose structure has previously been determined at a resolution of 2.5 Angstrom [Kong et al. (1992), Cell, 69, 425-437]. Here, the construction of a new plasmid that directs overproduction of beta to very high levels and a simple procedure for large-scale purification of the protein are described. Crystals grown under slightly modified conditions diffracted to beyond 1.9 Angstrom at 100 K at a synchrotron source. The structure of the beta dimer solved at 1.85 Angstrom resolution shows some differences from that reported previously. In particular, it was possible at this resolution to identify residues that differed in position between the two subunits in the unit cell; side chains of these and some other residues were found to occupy alternate conformations. This suggests that these residues are likely to be relatively mobile in solution. Some implications of this flexibility for the function of beta are discussed.

A framework for specification-based testing

Test templates and a test template framework are introduced as useful concepts in specification-based testing. The framework can be defined using any model-based specification notation and used to derive tests from model-based specifications-in this paper, it is demonstrated using the Z notation. The framework formally defines test data sets and their relation to the operations in a specification and to other test data sets, providing structure to the testing process. Flexibility is preserved, so that many testing strategies can be used. Important application areas of the framework are discussed, including refinement of test data, regression testing, and test oracles.

Protein crystallography and examples of its applications in medicinal chemistry

The field of protein crystallography inspires and enthrals, whether it be for the beauty and symmetry of a perfectly formed protein crystal, the unlocked secrets of a novel protein fold, or the precise atomic-level detail yielded from a protein-ligand complex. Since 1958, when the first protein structure was solved, there have been tremendous advances in all aspects of protein crystallography, from protein preparation and crystallisation through to diffraction data measurement and structure refinement. These advances have significantly reduced the time required to solve protein crystal structures, while at the same time substantially improving the quality and resolution of the resulting structures. Moreover, the technological developments have induced researchers to tackle ever more complex systems, including ribosomes and intact membrane-bound proteins, with a reasonable expectation of success. In this review, the steps involved in determining a protein crystal structure are described and the impact of recent methodological advances identified. Protein crystal structures have proved to be extraordinarily useful in medicinal chemistry research, particularly with respect to inhibitor design. The precise interaction between a drug and its receptor can be visualised at the molecular level using protein crystal structures, and this information then used to improve the complementarity and thus increase the potency and selectivity of an inhibitor. The use of protein crystal structures in receptor-based drug design is highlighted by (i) HIV protease, (ii) influenza virus neuraminidase and (iii) prostaglandin H-2-synthetase. These represent, respectively, examples of protein crystal structures that (i) influenced the design of drugs currently approved for use in the treatment of HIV infection, (ii) led to the design of compounds currently in clinical trials for the treatment of influenza infection and (iii) could enable the design of highly specific non-steroidal anti-inflammatory drugs that lack the common side-effects of this drug class.