Focusing in on structural genomics: The University of Queensland structural biology pipeline

The flood of new genomic sequence information together with technological innovations in protein structure determination have led to worldwide structural genomics (SG) initiatives. The goals of SG initiatives are to accelerate the process of protein structure determination, to fill in protein fold space and to provide information about the function of uncharacterized proteins. In the long-term, these outcomes are likely to impact on medical biotechnology and drug discovery, leading to a better understanding of disease as well as the development of new therapeutics. Here we describe the high throughput pipeline established at the University of Queensland in Australia. In this focused pipeline, the targets for structure determination are proteins that are expressed in mouse macrophage cells and that are inferred to have a role in innate immunity. The aim is to characterize the molecular structure and the biochemical and cellular function of these targets by using a parallel processing pipeline. The pipeline is designed to work with tens to hundreds of target gene products and comprises target selection, cloning, expression, purification, crystallization and structure determination. The structures from this pipeline will provide insights into the function of previously uncharacterized macrophage proteins and could lead to the validation of new drug targets for chronic obstructive pulmonary disease and arthritis. (c) 2006 Elsevier B.V. All rights reserved.

Phosphate forms an unusual tripodal complex with the Fe-Mn center of sweet potato purple acid phosphatase

Purple acid phosphatases (PAPs) are a family of binuclear metalloenzymes that catalyze the hydrolysis of phosphoric acid esters and anhydrides. A PAP in sweet potato has a unique, strongly antiferromagnetically coupled Fe(III)-Mn(II) center and is distinguished from other PAPs by its increased catalytic efficiency for a range of activated and unactivated phosphate esters, its strict requirement for Mn(II), and the presence of a mu-oxo bridge at pH 4.90. This enzyme displays maximum catalytic efficiency (k(cat)/K-m) at pH 4.5, whereas its catalytic rate constant (k(cat)) is maximal at near-neutral pH, and, in contrast to other PAPs, its catalytic parameters are not dependent on the pK(a) of the leaving group. The crystal structure of the phosphate-bound Fe(III)-Mn(II) PAP has been determined to 2.5-Angstrom resolution (final R-free value of 0.256). Structural comparisons of the active site of sweet potato, red kidney bean, and mammalian PAPs show several amino acid substitutions in the sweet potato enzyme that can account for its increased catalytic efficiency. The phosphate molecule binds in an unusual tripodal mode to the two metal ions, with two of the phosphate oxygen atoms binding to Fe(III) and Mn(II), a third oxygen atom bridging the two metal ions, and the fourth oxygen pointing toward the substrate binding pocket. This binding mode is unique among the known structures in this family but is reminiscent of phosphate binding to urease and of sulfate binding to A protein phosphatase. The structure and kinetics support the hypothesis that the bridging oxygen atom initiates hydrolysis.

Structure-activity relationships for a new family of sulfonylurea herbicides

Acetohydroxyacid synthase (AHAS; EC 2.2.1.6) catalyzes the first common step in branched-chain amino acid biosynthesis. The enzyme is inhibited by several chemical classes of compounds and this inhibition is the basis of action of the sulfonylurea and imidazolinone herbicides. The commercial sulfonylureas contain a pyrimidine or a triazine ring that is substituted at both meta positions, thus obeying the initial rules proposed by Levitt. Here we assess the activity of 69 monosubstituted sulfonylurea analogs and related compounds as inhibitors of pure recombinant Arabidopsis thaliana AHAS and show that disubstitution is not absolutely essential as exemplified by our novel herbicide, monosulfuron (2-nitro-N-(4'-methyl-pyrimidin-2'-yl) phenyl-sulfonylurea), which has a pyrimidine ring with a single meta substituent. A subset of these compounds was tested for herbicidal activity and it was shown that their effect in vivo correlates well with their potency in vitro as AHAS inhibitors. Three-dimensional quantitative structure-activity relationships were developed using comparative molecular field analysis and comparative molecular similarity indices analysis. For the latter, the best result was obtained when steric, electrostatic, hydrophobic and H-bond acceptor factors were taken into consideration. The resulting fields were mapped on to the published crystal structure of the yeast enzyme and it was shown that the steric and hydrophobic fields are in good agreement with sulfonylurea-AHAS interaction geometry.

Herbicide-binding Sites Revealed in the Structure of Plant Acetohydroxyacid Synthase

The sulfonylureas and imidazolinones are potent commercial herbicide families. They are among the most popular choices for farmers worldwide, because they are nontoxic to animals and highly selective. These herbicides inhibit branched-chain amino acid biosynthesis in plants by targeting acetohydroxyacid synthase (AHAS, EC 2.2.1.6). This report describes the 3D structure of Arabidopsis thaliana AHAS in complex with five sulfonylureas (to 2.5 angstrom resolution) and with the imidazolinone, imazaquin (IQ; 2.8 angstrom). Neither class of molecule has a structure that mimics the substrates for the enzyme, but both inhibit by blocking a channel through which access to the active site is gained. The sulfonylureas approach within 5 angstrom of the catalytic center, which is the C2 atom of the cofactor thiamin diphosphate, whereas IQ is at least 7 angstrom from this atom. Ten of the amino acid residues that bind the sulfonylureas also bind IQ. Six additional residues interact only with the sulfonylureas, whereas there are two residues that bind IQ but not the sulfonylureas. Thus, the two classes of inhibitor occupy partially overlapping sites but adopt different modes of binding. The increasing emergence of resistant weeds due to the appearance of mutations that interfere with the inhibition of AHAS is now a worldwide problem. The structures described here provide a rational molecular basis for understanding these mutations, thus allowing more sophisticated AHAS inhibitors to be developed. There is no previously described structure for any plant protein in complex with a commercial herbicide.

Antimicrobial peptides in action

Molecular dynamics simulations of the magainin MG-H2 peptide interacting with a model phospholipid membrane have been used to investigate the mechanism by which antimicrobial peptides act. Multiple copies of the peptide were randomly placed in solution close to the membrane. The peptide readily bound to the membrane, and above a certain concentration, the peptide was observed to cooperatively induce the formation of a nanometer- sized, toroidally shaped pore in the bilayer. In sharp contrast with the commonly accepted model of a toroidal pore, only one peptide was typically found near the center of the pore. The remaining peptides lay close to the edge of the pore, maintaining a predominantly parallel orientation with respect to the membrane.

Crystallization and preliminary X-ray analysis of a Kunitz-type inhibitor, textilinin-1 from Pseudonaja textilis textilis

Textilinin-1 (Txln-1), a Kunitz-type serine protease inhibitor, is a 59-amino-acid polypeptide isolated from the venom of the Australian Common Brown snake Pseudonaja textilis textilis. This molecule has been suggested as an alternative to aprotinin, also a Kunitz-type serine protease inhibitor, for use as an anti-bleeding agent in surgical procedures. Txln-1 shares only 47% amino-acid identity to aprotinin; however, six cysteine residues in the two peptides are in conserved locations. It is therefore expected that the overall fold of these molecules is similar but that they have contrasting surface features. Here, the crystallization of recombinant textilinin-1 (rTxln-1) as the free molecule and in complex with bovine trypsin (229 amino acids) is reported. Two organic solvents, phenol and 1,4-butanediol, were used as additives to facilitate the crystallization of free rTxln-1. Crystals of the rTxln-1-bovine trypsin complex diffracted to 2.0 angstrom resolution, while crystals of free rTxln-1 diffracted to 1.63 angstrom resolution.

Identification and molecular modeling of a novel, plant-like, human purple acid phosphatase

Purple acid phosphatases are a family of binuclear metallohydrolases that have been identified in plants, animals and fungi. Only one isoform of similar to 35 kDa has been isolated from animals, where it is associated with bone resorption and microbial killing through its phosphatase activity, and hydroxyl radical production, respectively. Using the sensitive PSI-BLAST search method, sequences representing new purple acid phosphatase-like proteins have been identified in mammals, insects and nematodes. These new putative isoforms are closely related to the similar to 55 kDa purple acid phosphatase characterized from plants. Secondary structure prediction of the new human isoform further confirms its similarity to a purple acid phosphatase from the red kidney bean. A structural model for the human enzyme was constructed based on the red kidney bean purple acid phosphatase structure. This model shows that the catalytic centre observed in other purple acid phosphatases is also present in this new isoform. These observations suggest that the sequences identified in this study represent a novel subfamily of plant-like purple acid phosphatases in animals and humans. (c) 2006 Elsevier B.V. All rights reserved.

Molecular dynamics simulations of the hydrophobin SC3 at a hydrophobic/hydrophilic interface

Hydrophobins are small (similar to 100 aa) proteins that have an important role in the growth and development of mycelial fungi. They are surface active and, after secretion by the fungi, self-assemble into amphipathic membranes at hydrophobic/hydrophilic interfaces, reversing the hydrophobicity of the surface. In this study, molecular dynamics simulation techniques have been used to model the process by which a specific class I hydrophobin, SC3, binds to a range of hydrophobic/ hydrophilic interfaces. The structure of SC3 used in this investigation was modeled based on the crystal structure of the class II hydrophobin HFBII using the assumption that the disulfide pairings of the eight conserved cysteine residues are maintained. The proposed model for SC3 in aqueous solution is compact and globular containing primarily P-strand and coil structures. The behavior of this model of SC3 was investigated at an air/water, an oil/water, and a hydrophobic solid/water interface. It was found that SC3 preferentially binds to the interfaces via the loop region between the third and fourth cysteine residues and that binding is associated with an increase in a-helix formation in qualitative agreement with experiment. Based on a combination of the available experiment data and the current simulation studies, we propose a possible model for SC3 self-assembly on a hydrophobic solid/water interface.

Determination of the orientational distribution and orientation factor for transfer between membrane-bound fluorophores using a confocal microscope

Orientational fluorophores have been a useful tool in physical chemistry, biochemistry, and more recently structural biology due to the polarized nature of the light they emit and that fact that energy can be transferred between them. We present a practical scheme in which measurements of the intensity of emitted fluorescence can be used to determine limits on the mean and distribution of orientation of the absorption transition moment of membrane-bound. uorophores. We demonstrate how information about the orientation of. uorophores can be used to calculate the orientation factor k(2) required for use in FRET spectroscopy. We illustrate the method using images of AlexaFluor probes bound to MscL mechanosensitive transmembrane channel proteins in spherical liposomes.

Direct protein interaction underlies gene-for-gene specificity and coevolution of the flax resistance genes and flax rust avirulence genes

Plant resistance proteins (R proteins) recognize corresponding pathogen avirulence (Avr) proteins either indirectly through detection of changes in their host protein targets or through direct R-Avr protein interaction. Although indirect recognition imposes selection against Avr effector function, pathogen effector molecules recognized through direct interaction may overcome resistance through sequence diversification rather than loss of function. Here we show that the flax rust fungus AvrLS67 genes, whose products are recognized by the L5, L6, and L7 R proteins of flax, are highly diverse, with 12 sequence variants identified from six rust strains. Seven AvrL567 variants derived from Avr alleles induce necrotic responses when expressed in flax plants containing corresponding resistance genes (R genes), whereas five variants from avr alleles do not. Differences in recognition specificity between AvA567 variants and evidence for diversifying selection acting on these genes suggest they have been involved in a gene-specific arms race with the corresponding flax R genes. Yeast two-hybrid assays indicate that recognition is based on direct R-Avr protein interaction and recapitulate the interaction specificity observed in planta. Biochemical analysis of Escherichia coli-produced AvrL567 proteins shows that variants that escape recognition nevertheless maintain a conserved structure and stability, suggesting that the amino acid sequence differences directly affect the R-Avr protein interaction. We suggest that direct recognition associated with high genetic diversity at corresponding R and Avr gene loci represents an alternative outcome of plant-pathogen coevolution to indirect recognition associated with simple balanced polymorphisms for functional and nonfunctional R and Avr genes.

Functional recycling of C2 domains throughout evolution: A comparative study of synaptotagmin, protein kinase C and phospholipase C by sequence, structural and modelling approaches

The C2 domain is one of the most frequent and widely distributed calcium-binding motifs. Its structure comprises an eight-stranded beta-sandwich with two structural types as if the result of a circular permutation. Combining sequence, structural and modelling information, we have explored, at different levels of granularity, the functional characteristics of several families of C2 domains. At the coarsest level,the similarity correlates with key structural determinants of the C2 domain fold and, at the finest level, with the domain architecture of the proteins containing them, highlighting the functional diversity between the various subfamilies. The functional diversity appears as different conserved surface patches throughout this common fold. In some cases, these patches are related to substrate-binding sites whereas in others they correspond to interfaces of presumably permanent interaction between other domains within the same polypeptide chain. For those related to substrate-binding sites, the predictions overlap with biochemical data in addition to providing some novel observations. For those acting as protein-protein interfaces' our modelling analysis suggests that slight variations between families are a result of not only complementary adaptations in the interfaces involved but also different domain architecture. In the light of the sequence and structural genomic projects, the work presented here shows that modelling approaches along with careful sub-typing of protein families will be a powerful combination for a broader coverage in proteomics. (C) 2003 Elsevier Ltd. All rights reserved.

Structural basis and prediction of substrate specificity in protein serine/threonine kinases

The large number of protein kinases makes it impractical to determine their specificities and substrates experimentally. Using the available crystal structures, molecular modeling, and sequence analyses of kinases and substrates, we developed a set of rules governing the binding of a heptapeptide substrate motif (surrounding the phosphorylation site) to the kinase and implemented these rules in a web-interfaced program for automated prediction of optimal substrate peptides, taking only the amino acid sequence of a protein kinase as input. We show the utility of the method by analyzing yeast cell cycle control and DNA damage checkpoint pathways. Our method is the only available predictive method generally applicable for identifying possible substrate proteins for protein serine/threonine kinases and helps in silico construction of signaling pathways. The accuracy of prediction is comparable to the accuracy of data from systematic large-scale experimental approaches.

Structural basis for the specificity of bipartite nuclear localization sequence binding by importin-α

Importin-alpha is the nuclear import receptor that recognizes cargo proteins carrying conventional basic monopartite and bipartite nuclear localization sequences (NLSs) and facilitates their transport into the nucleus. Bipartite NLSs contain two clusters of basic residues, connected by linkers of variable lengths. To determine the structural basis of the recognition of diverse bipartite NLSs by mammalian importin-alpha, we co-crystallized a non-autoinhibited mouse receptor protein with peptides corresponding to the NLSs from human retinoblastoma protein and Xenopus laevis phosphoprotein N1N2, containing diverse sequences and lengths of the linker. We show that the basic clusters interact analogously in both NLSs, but the linker sequences adopt different conformations, whereas both make specific contacts with the receptor. The available data allow us to draw general conclusions about the specificity of NLS binding by importin-alpha and facilitate an improved definition of the consensus sequence of a conventional basic/bipartite NLS (KRX10-12KRRK) that can be used to identify novel nuclear proteins.

Twists, knots, and rings in proteins: Structural definition of the cyclotide framework

In recent years an increasing number of miniproteins containing an amide-cyclized backbone have been discovered. The cyclotide family is the largest group of such proteins and is characterized by a circular protein backbone and six conserved cysteine residues linked by disulfide bonds in a tight core of the molecule. These form a cystine knot in which an embedded ring formed by two of the disulfide bonds and the connecting backbone segment is threaded by a third disulfide bond. In the current study we have undertaken high resolution structural analysis of two prototypic cyclotides, kalata B1 and cycloviolacin O1, to define the role of the conserved residues in the sequence. We provide the first comprehensive analysis of the topological features in this unique family of proteins, namely rings (a circular backbone), twists (a cis-peptide bond in the Mobius cyclotides) and knots (a knotted arrangement of the disulfide bonds).