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 initiation sites for T4 lysozyme folding using CD and NMR spectroscopy of peptide fragments

Using CD and 2D H-1 NMR spectroscopy, we have identified potential initiation sites for the folding of T4 lysozyme by examining the conformational preferences of peptide fragments corresponding to regions of secondary structure. CD spectropolarimetry showed most peptides were unstructured in water, but adopted partial helical conformations in TFE and SDS solution. This was also consistent with the H-1 NMR data which showed that the peptides were predominantly disordered in water, although in some cases, nascent or small populations of partially folded conformations could be detected. NOE patterns, coupling constants, and deviations from random coil Her chemical shift values complemented the CD data and confirmed that many of the peptides were helical in TFE and SDS micelles. In particular, the peptide corresponding to helix E in the native enzyme formed a well-defined helix in both TFE and SDS, indicating that helix E potentially forms an initiation site for T4 lysozyme folding. The data for the other peptides indicated that helices D, F, G, and H are dependent on tertiary interactions for their folding and/or stability. Overall, the results from this study, and those of our earlier studies, are in agreement with modeling and IID-deuterium exchange experiments, and support an hierarchical model of folding for T4 lysozyme.

Functional implications of the human T-lymphotropic virus type 1 transmembrane glycoprotein helical hairpin structure

Retrovirus entry into cells follows receptor binding by the surface exposed envelope glycoprotein (Env) subunit (SU), which triggers the membrane fusion activity of the transmembrane (TM) protein. TM protein fragments expressed in the absence of SU adopt helical hairpin structures comprising a central coiled coil, a region of chain reversal containing a disulfide-bonded loop, and a C-terminal segment that packs onto the exterior of the coiled coil in an antiparallel manner. Here we used in vitro mutagenesis to test the functional role of structural elements observed in a model helical hairpin, gp21 of human T-lymphotropic virus type 1. Membrane fusion activity requires the stabilization of the N and C termini of the central coiled coil by a hydrophobic N cap and a small hydrophobic core, respectively. A conserved Gly-Gly hinge motif preceding the disulfide-bonded loop, a salt bridge that stabilizes the chain reversal region, and interactions between the C-terminal segment and the coiled coil are also critical for fusion activity. Our data support a model whereby the chain reversal region transmits a conformational signal from receptor-bound SU to induce the fusion-activated helical hairpin conformation of the TM protein.

Conotoxin TVIIA, a novel peptide from the venom of Conus tulipa - 2. Three-dimensional solution structure

The three-dimensional solution structure of conotoxin TVIIA, a 30-residue polypeptide from the venom of the piscivorous cone snail Conus tulipa, has been determined using 2D H-1 NMR spectroscopy. TVIIA contains six cysteine residues which form a 'four-loop' structural framework common to many peptides from Conus venoms including the omega-, delta-, kappa-, and mu O-conotoxins. However, TVIIA does not belong to these well-characterized pharmacological classes of conotoxins, but displays high sequence identity with conotoxin GS, a muscle sodium channel blocker from Conus geographus. Structure calculations were based on 562 interproton distance restraints inferred from NOE data, together with 18 backbone and nine side-chain torsion angle restraints derived from spin-spin coupling constants. The final family of 20 structures had mean pairwise rms differences over residues 2-27 of 0.18 +/- 0.05 Angstrom for the backbone atoms and 1.39 +/- 0.33 Angstrom for all heavy atoms. The structure consists of a triple-stranded, antiparallel beta sheet with +2x, -1 topology (residues 7-9, 16-20 and 23-27) and several beta turns. The core of the molecule is formed by three disulfide bonds which form a cystine knot motif common to many toxic and inhibitory polypeptides. The global fold, molecular shape and distribution of amino-acid sidechains in TVIIA is similar to that previously reported for conotoxin GS, and comparison with other four-loop conotoxin structures provides further indication that TVIIA and GS represent a new and distinct subgroup of this structural family. The structure of TVIIA determined in this study provides the basis for determining a structure-activity relationship for these molecules and their interaction with target receptors.

Structural interpretation of mutations in phenylalanine hydroxylase protein aids in identifying genotype-phenotype correlations in phenylketonuria

Phenylalanine hydroxylase (PAH) is the enzyme that converts phenylalanine to tyrosine as a rate-limiting step in phenylalanine catabolism and protein and neurotransmitter biosynthesis. Over 300 mutations have been identified in the gene encoding PAH that result in a deficient enzyme activity and lead to the disorders hyperphenylalaninaemia and phenylketonuria. The determination of the crystal structure of PAH now allows the determination of the structural basis of mutations resulting in PAH deficiency. We present an analysis of the structural basis of 120 mutations with a 'classified' biochemical phenotype and/or available in vitro expression data. We find that the mutations can be grouped into five structural categories, based on the distinct expected structural and functional effects of the mutations in each category. Missense mutations and small amino acid deletions are found in three categories:'active site mutations', 'dimer interface mutations', and 'domain structure mutations'. Nonsense mutations and splicing mutations form the category of 'proteins with truncations and large deletions'. The final category, 'fusion proteins', is caused by frameshift mutations. We show that the structural information helps formulate some rules that will help predict the likely effects of unclassified and newly discovered mutations: proteins with truncations and large deletions, fusion proteins and active site mutations generally cause severe phenotypes; domain structure mutations and dimer interface mutations spread over a range of phenotypes, but domain structure mutations in the catalytic domain are more likely to be severe than domain structure mutations in the regulatory domain or dimer interface mutations.

Th structure of the P-II-ATP complex

P-II is a signal transduction protein that is part of the cellular machinery used by many bacteria to regulate the activity of glutamine synthetase and the transcription of its gene. The structure of P-II was solved using a hexagonal crystal form (form I). The more physiologically relevant form of P-II is a complex with small molecule effecters. We describe the structure of P-II with ATP obtained by analysis of two different crystal forms (forms II and III) that were obtained by co-crystallization of P-II with ATP. Both structures have a disordered recognition (T) loop and show differences at their C termini. Comparison of these structures with the form I protein reveals changes that occur on binding ATP. Surprisingly, the structure of the P-II/ATP complex differs with that of GlnK, a functional homologue. The two proteins bind the base and sugar of ATP in a similar manner but show differences in the way that they interact with the phosphates. The differences in structure could account for the differences in their activities, and these have been attributed to a difference in sequence at position 82. It has been demonstrated recently that P-II and GlnK form functional heterotrimers in vivo. We construct models of the heterotrimers and examine the junction between the subunits.

In vivo protein cyclization promoted by a circularly permuted Synechocystis sp PCC6803 DnaB mini-intein

A synthetic Synechocystis sp. PCC6803 DnaB split mini-intein gene was constructed for the in vivo cyclization of recombinant proteins expressed in Escherichia coli. The system was used to cyclize the NH2-terminal domain of E. coli DnaB, the structure of which had been determined previously by NMR spectroscopy. Cyclization was found to proceed efficiently, with little accumulation of precursor, and the product was purified in high yield. The solution structure of cyclic DnaB-N is not significantly different from that of linear DnaB-N and it unfolds reversibly at temperatures similar to14 degreesC higher. Improved hydrogen bonding was observed in the first and last helices, and the length of the last helix was increased, while the 9-amino acid linker used to join the NH2 and COOH termini was found to be highly mobile. The measured thermodynamic stabilization of the structure (DeltaDeltaG approximate to 2 kcal/mol) agrees well with the value estimated from the reduced conformational entropy in the unfolded form. Simple polymer theory can be used to predict likely free energy changes resulting from protein cyclization and how the stabilization depends on the size of the protein and the length of the linker used to connect the termini.

A computational analysis of substrate binding strength by phosphorylase kinase and protein kinase A

Protein kinases exhibit various degrees of substrate specificity. The large number of different protein kinases in the eukaryotic proteomes makes it impractical to determine the specificity of each enzyme experimentally. To test if it were possible to discriminate potential substrates from non-substrates by simple computational techniques, we analysed the binding enthalpies of modelled enzyme-substrate complexes and attempted to correlate it with experimental enzyme kinetics measurements. The crystal structures of phosphorylase kinase and cAMP-dependent protein kinase were used to generate models of the enzyme with a series of known peptide substrates and non-substrates, and the approximate enthalpy of binding assessed following energy minimization. We show that the computed enthalpies do not correlate closely with kinetic measurements, but the method can distinguish good substrates from weak substrates and non-substrates. Copyright (C) 2002 John Wiley Sons, Ltd.

Sequence analysis and expression of a virus-like particle protein,VLP2, from the parasitic wasp Venturia canescens

Endoparasitoid wasps produce maternal protein secretions, which are transported into the body of insect hosts at oviposition to regulate host physiology for successful development of their offspring. Venturia canescens calyx fluid contains so-called virus-like particles (VLPs) that are essential for immune evasion of the developing parasitoid inside the host. VLPs consist of four major proteins. In this paper, we describe the isolation and molecular cloning of a gene (vlp2) that is a constituent of VLPs and discuss its possible role in VLP structure and function.

Homomeric ring assemblies of eukaryotic Sm proteins have affinity for both RNA and DNA - Crystal structure of an oligomeric complex of yeast SmF

Sm and Sm-like proteins are key components of small ribonucleoproteins involved in many RNA and DNA processing pathways. In eukaryotes, these complexes contain seven unique Sm or Sm-like (Lsm) proteins assembled as hetero-heptameric rings, whereas in Archaea and bacteria six or seven-membered rings are made from only a single polypeptide chain. Here we show that single Sm and Lsm proteins from yeast also have the capacity to assemble into homo-oligomeric rings. Formation of homo-oligomers by the spliceosomal small nuclear ribonucleoprotein components SmE and SmF preclude hetero-interactions vital to formation of functional small nuclear RNP complexes in vivo. To better understand these unusual complexes, we have determined the crystal structure of the homomeric assembly of the spliceosomal protein SmF. Like its archaeal/bacterial homologs, the SmF complex forms a homomeric ring but in an entirely novel arrangement whereby two heptameric rings form a co-axially stacked dimer via interactions mediated by the variable loops of the individual SmF protein chains. Furthermore, we demonstrate that the homomeric assemblies of yeast Sm and Lsm proteins are capable of binding not only to oligo(U) RNA but, in the case of SmF, also to oligo(dT) single-stranded DNA.

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.

The crystal structure of free human hypoxanthineguanine phosphoribosyltransferase reveals extensive conformational plasticity throughout the catalytic cycle

Human hypoxanthine-guanine phosphoribosyltransferase (HGPRT) catalyses the synthesis of the purine nucleoside monophosphates, IMP and GMP, by the addition of a 6-oxopurine base, either hypoxanthine or guanine, to the 1-beta-position of 5-phospho-U-D-ribosyl-1-pyrophosphate (PRib-PP). The mechanism is sequential, with PRib-PP binding to the free enzyme prior to the base. After the covalent reaction, pyrophosphate is released followed by the nucleoside monophosphate. A number of snapshots of the structure of this enzyme along the reaction pathway have been captured. These include the structure in the presence of the inactive purine base analogue, 7-hydroxy [4,3-d] pyrazolo pyrimidine (HPP) and PRib-PP. Mg2+, and in complex with IMP or GMP. The third structure is that of the immucillinHP.Mg2+.PPi complex, a transition-state analogue. Here, the first crystal structure of free human HGPRT is reported to 1.9 angstrom resolution, showing that significant conformational changes have to occur for the substrate(s) to bind and for catalysis to proceed. Included in these changes are relative movement of subunits within the tetramer, rotation and extension of an active-site alpha-helix (D137-D153), reorientation of key active-site residues K68, D137 and K165, and the rearrangement of three active-site loops (100-128, 165-173 and 186-196). Toxoplasina gondii HGXPRT is the only other 6-oxopurine phosphoribosyltransferase structure solved in the absence of ligands. Comparison of this structure with human HGPRT reveals significant differences in the two active sites, including the structure of the flexible loop containing K68 (human) or K79 (T gondii). (c) 2005 Elsevier Ltd. All rights reserved.

The 2.0 angstrom crystal structure of a pocilloporin at pH 3.5: The structural basis for the linkage between color transition and halide binding

The pocilloporin Rtms5 and an engineered variant Rtms5(H146S) undergo distinct color transitions (from blue to red to yellow to colorless) in a pH-dependent manner. pK(a) values of 4.1 and 3.2 were determined for the blue (absorption lambda(max), 590 nm) to yellow (absorption lambda(max), similar to 453 nm) transitions of Rtms5 and Rtms5H(146). The pK(a) for the blue-yellow transition of Rtms5H(146S) increased by 1.4 U in the presence of 0.1 M KI, whereas the pK(a) for the same transition of Rtms5 was relatively insensitive to added halides. To understand the structural basis for these observations, we have determined to 2.0 A resolution the crystal structure of a yellow form of Rtms5(H146S) at pH 3.5 in the presence of iodide. Iodide was found occupying a pocket in the structure with a pH of 3.5, forming van der Waals contacts with the tyrosyl moiety of the chromophore. Elsewhere, it was determined that this pocket is occupied by a water molecule in the Rtms5(H141S) structure (pH 8.0) and by the side chain of histidine 146 in the wild-type Rtms5 structure. Collectively, our data provide an explanation for the observed linkage between color transitions for Rtms5(H146S) and binding to halides.

Giant protein kinases: Domain interactions and structural basis of autoregulation

The myosin-associated giant protein kinases twitchin and titin are composed predominantly of fibronectin- and immunoglobulin-like modules, We report the crystal structures of two autoinhibited twitchin kinase fragments, one from Aplysia and a larger fragment from Caenorhabditis elegans containing an additional C-terminal immunoglobulin-like domain, The structure of the longer fragment shoes that the immunoglobulin domain contacts the protein kinase domain on the opposite side from the catalytic cleft, laterally exposing potential myosin binding residues, Together, the structures reveal the cooperative interactions between the autoregulatory region and the residues from the catalytic domain involved in protein substrate binding, ATP binding, catalysis and the activation loop, and explain the differences between the observed autoinhibitory mechanism and the one found in the structure of calmodulin-dependent kinase I.