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.

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.

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.

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.

Structure of TcpG, the DsbA protein folding catalyst from Vibrio cholerae

The efficient and correct folding of bacterial disulfide bonded proteins in vivo is dependent upon a class of periplasmic oxidoreductase proteins called DsbA, after the Escherichia coli enzyme. In the pathogenic bacterium Vibrio cholerae, the DsbA homolog (TcpG) is responsible for the folding, maturation and secretion of virulence factors. Mutants in which the tcpg gene has been inactivated are avirulent; they no longer produce functional colonisation pill and they no longer secrete cholera toxin. TcpG is thus a suitable target for inhibitors that could counteract the virulence of this organism, thereby preventing the symptoms of cholera. The crystal structure of oxidized TcpG (refined at a resolution of 2.1 Angstrom) serves as a starting point for the rational design of such inhibitors. As expected, TcpG has the same fold as E. coli DsbA, with which it shares similar to 40% sequence identity. Ln addition, the characteristic surface features of DsbA are present in TcpG, supporting the notion that these features play a functional role. While the overall architecture of TcpG and DsbA is similar and the surface features are retained in TcpG, there are significant differences. For example, the kinked active site helix results from a three-residue loop in DsbA, but is caused by a proline in TcpG (making TcpG more similar to thioredoxin in this respect). Furthermore, the proposed peptide binding groove of TcpG is substantially shortened compared with that of DsbA due to a six-residue deletion. Also, the hydrophobic pocket of TcpG is more shallow and the acidic patch is much less extensive than that of E. coli DsbA. The identification of the structural and surface features that are retained or are divergent in TcpG provides a useful assessment of their functional importance in these protein folding catalysts and is an important prerequisite for the design of TcpG inhibitors. (C) 1997 Academic Press Limited.

Estimation of body water compartments in cirrhosis by multiple-frequency bioelectrical-impedance analysis

Estimation of total body water by measuring bioelectrical impedance at a fixed frequency of 50 kHz is useful in assessing body composition in healthy populations. However, in cirrhosis, the distribution of total body water between the extracellular and intracellular compartments is of greater clinical importance. We report an evaluation of a new multiple-frequency bioelectrical-impedance analysis technique (MFBIA) that may quantify the distribution of total body water in cirrhosis. In 21 cirrhotic patients and 21 healthy control subjects, impedance to the Row of current was measured at frequencies ranging from 4 to 1012 kHz. These measurements were used to estimate body water compartments and then compared with total body water and extracellular water determined by isotope methodology. In cirrhotic patients, extracellular water and total body water (as determined by isotope methods) were well predicted by MFBIA (r = 0.73 and 0.89, respectively).;However, the 95% confidence intervals of the limits of agreement between MFBIA and the isotope methods were +/- 14% and +/-9% for cirrhotics (extracellular water and total body water, respectively) and +/-9% and +/-9% for cirrhotics without ascites. The 95% confidence intervals estimated from the control group were +/-10% and +/-5% for extracellular water and total body water, respectively. Thus, despite strong correlations between MFBIA and isotope measurements, the relatively large limits of agreement with accepted techniques suggest that the MFBIA technique requires further refinement before it can be routinely used to determine the nutritional assessment of individual cirrhotic patients. Nutrition 2001,17.31-34. (C)Elsevier Science Inc. 2001.

Three-dimensional structure of RTD-1, a cyclic antimicrobial defensin from rhesus macaque leukocytes

Most mammalian defensins are cationic peptides of 29-42 amino acids long, stabilized by three disulfide bonds. However, recently Tang et al. (1999, Science 286, 498-502) reported the isolation of a new defensin type found in the leukocytes of rhesus macaques. In contrast to all the other defensins found so far, rhesus theta defensin-1 (RTD-1) is composed of just 18 amino acids with the backbone cyclized through peptide bonds. Antibacterial activities of both the native cyclic peptide and a linear form were examined, showing that the cyclic form was 3-fold more active than the open chain analogue [Tang et al. (1999) Science 286, 498-502]. To elucidate the three-dimensional structure of RTD-1 and its open chain analogue, both peptides were synthesized using solid-phase peptide synthesis and tert-butyloxycarbonyl chemistry. The structures of both peptides in aqueous solution were determined from two-dimensional H-1 NMR data recorded at 500 and 750 MHz. Structural constraints consisting of interproton distances and dihedral angles were used as input for simulated-annealing calculations and water refinement with the program CNS. RTD-1 and its open chain analogue oRTD-1 adopt very similar structures in water. Both comprise an extended beta -hairpin structure with turns at one or both ends. The turns are well defined within themselves and seem to be flexible with respect to the extended regions of the molecules. Although the two strands of the beta -sheet are connected by three disulfide bonds, this region displays a degree of flexibility. The structural similarity of RTD-1 and its open chain analogue oRTD-1, as well as their comparable degree of flexibility, support the theory that the additional charges at the termini of the open chain analogue rather than overall differences in structure or flexibility are the cause for oRTD-1's lower antimicrobial activity. In contrast to numerous other antimicrobial peptides, RTD-1 does not display any amphiphilic character, even though surface models of RTD-1 exhibit a certain clustering of positive charges. Some amide protons of RTD-1 that should be solvent-exposed in monomeric beta -sheet structures show low-temperature coefficients, suggesting the possible presence of weak intermolecular hydrogen bonds.

Getting the adrenaline going: Crystal structure of the adrenaline-synthesizing enzyme PNMT

Background: Adrenaline is localized to specific regions of the central nervous system (CNS), but its role therein is unclear because of a lack of suitable pharmacologic agents. Ideally, a chemical is required that crosses the blood-brain barrier, potently inhibits the adrenaline-synthesizing enzyme PNMT, and does not affect other catecholamine processes. Currently available PNMT inhibitors do not meet these criteria. We aim to produce potent, selective, and CNS-active PNMT inhibitors by structure-based design methods. The first step is the structure determination of PNMT. Results: We have solved the crystal structure of human PNMT complexed with a cofactor product and a submicromolar inhibitor at a resolution of 2.4 Angstrom. The structure reveals a highly decorated methyltransferase fold, with an active site protected from solvent by an extensive cover formed from several discrete structural motifs. The structure of PNMT shows that the inhibitor interacts with the enzyme in a different mode from the (modeled) substrate noradrenaline. Specifically, the position and orientation of the amines is not equivalent. Conclusions: An unexpected finding is that the structure of PNMT provides independent evidence of both backward evolution and fold recruitment in the evolution of a complex enzyme from a simple fold. The proposed evolutionary pathway implies that adrenaline, the product of PNMT catalysis, is a relative newcomer in the catecholamine family. The PNMT structure reported here enables the design of potent and selective inhibitors with which to characterize the role of adrenaline in the CNS. Such chemical probes could potentially be useful as novel therapeutics.

A refinement calculus for logic programs

Existing refinement calculi provide frameworks for the stepwise development of imperative programs from specifications. This paper presents a refinement calculus for deriving logic programs. The calculus contains a wide-spectrum logic programming language, including executable constructs such as sequential conjunction, disjunction, and existential quantification, as well as specification constructs such as general predicates, assumptions and universal quantification. A declarative semantics is defined for this wide-spectrum language based on executions. Executions are partial functions from states to states, where a state is represented as a set of bindings. The semantics is used to define the meaning of programs and specifications, including parameters and recursion. To complete the calculus, a notion of correctness-preserving refinement over programs in the wide-spectrum language is defined and refinement laws for developing programs are introduced. The refinement calculus is illustrated using example derivations and prototype tool support is discussed.

Mapping glider habitat in dry eucalypt forests for Montreal Process indicator 1.1: Fragmentation of forest types

Accurate habitat mapping is critical to landscape ecological studies such as required for developing and testing Montreal Process indicator 1.1e, fragmentation of forest types. This task poses a major challenge to remote sensing, especially in mixedspecies, variable-age forests such as dry eucalypt forests of subtropical eastern Australia. In this paper, we apply an innovative approach that uses a small section of one-metre resolution airborne data to calibrate a moderate spatial resolution model (30 m resolution; scale 1:50 000) based on Landsat Thematic Mapper data to estimate canopy structural properties in St Marys State Forest, near Maryborough, south-eastern Queensland. The approach applies an image-processing model that assumes each image pixel is significantly larger than individual tree crowns and gaps to estimate crown-cover percentage, stem density and mean crown diameter. These parameters were classified into three discrete habitat classes to match the ecology of four exudivorous arboreal species (yellowbellied glider Petaurus australis, sugar glider P. breviceps, squirrel glider P. norfolcensis , and feathertail glider Acrobates pygmaeus), and one folivorous arboreal marsupial, the greater glider Petauroides volans. These species were targeted due to the known ecological preference for old trees with hollows, and differences in their home range requirements. The overall mapping accuracy, visually assessed against transects (n = 93) interpreted from a digital orthophoto and validated in the field, was 79% (KHAT statistic = 0.72). The KHAT statistic serves as an indicator of the extent that the percentage correct values of the error matrix are due to ‘true’ agreement verses ‘chance’ agreement. This means that we are able to reliably report on the effect of habitat loss on target species, especially those with a large home range size (e.g. yellow-bellied glider). However, the classified habitat map failed to accurately capture the spatial patterning (e.g. patch size and shape) of stands with a trace or sub-dominance of senescent trees. This outcome makes the reporting of the effects of habitat fragmentation more problematic, especially for species with a small home range size (e.g. feathertail glider). With further model refinement and validation, however, this moderateresolution approach offers an important, cost eff e c t i v e advancement in mapping the age of dry eucalypt forests in the region.

A formal model of real-time program compilation

Program compilation can be formally defined as a sequence of equivalence-preserving transformations, or refinements, from high-level language programs to assembler code, Recent models also incorporate timing properties, but the resulting formalisms are intimidatingly complex. Here we take advantage of a new, simple model of real-time refinement, based on predicate transformer semantics, to present a straightforward compilation formalism that incorporates real-time constraints. (C) 2002 Elsevier Science B.V. All rights reserved.