Comparative protein structure modeling of genes and genomes

Cytoskeletal protein 4.1G binds to the third intracellular loop of the A1 adenosine receptor and inhibits receptor action

To identify binding partners of the A1AR (A1 adenosine receptor), yeast two-hybrid screening of a rat embryonic cDNA library was performed. This procedure led to the identification of erythrocyte membrane cytoskeletal protein (represented as 4.1G) as an A1AR-binding partner. Truncation studies revealed that the C-terminal domain of 4.1G was essential for binding to A1ARs and that the C-terminal domain of 4.1G and the third intracellular loop of A1ARs interacted. A1AR-4.1G interaction was also confirmed in studies using brain tissue. Studies in HEK-293 (human embryonic kidney 293) cells and Chinese-hamster ovary cells showed that 4.1G interfered with A1AR signal transduction, as 4.1G reduced A1AR-mediated inhibition of cAMP accumulation and intracellular calcium release. 4.1G also altered cell-surface A1AR expression. These observations identify 4.1G as a novel A1AR-binding partner that can regulate adenosine action.

Structural insights into molecular function of the metastasis-associated phosphatase PRL-3

Phosphatases and kinases are the cellular signal transduction enzymes that control protein phosphorylation. PRL phosphatases constitute a novel class of small (20 kDa), prenylated phosphatases with oncogenic activity. In particular, PRL-3 is consistently overexpressed in liver metastasis in colorectal cancer cells and represents a new therapeutic target. Here, we present the solution structure of PRL-3, the first structure of a PRL phosphatase. The structure places PRL phosphatases in the class of dual specificity phosphatases with closest structural homology to the VHR phosphatase. The structure, coupled with kinetic studies of site-directed mutants, identifies functionally important residues and reveals unique features, differentiating PRLs from other phosphatases. These differences include an unusually hydrophobic active site without the catalytically important serine/threonine found in most other phosphatases. The position of the general acid loop indicates the presence of conformational change upon catalysis. The studies also identify a potential regulatory role of Cys(49) that forms an intramolecular disulfide bond with the catalytic Cys(104) even under mildly reducing conditions. Molecular modeling of the highly homologous PRL-1 and PRL-2 phosphatases revealed unique surface elements that are potentially important for specificity.

Structural and functional analysis of integrin alpha2I domain interaction with echovirus 1

Integrins are cell surface receptors for several microbial pathogens including echovirus 1 (EV1), a picornavirus. Cryo-electron microscopy revealed that the functional domain (alpha(2)I) of human alpha(2)beta(1) integrin binds to a surface depression on the EV1 capsid. This three-dimensional structure of EV1 bound to alpha(2)I domain provides the first structural details of an integrin interacting with a picornavirus. The model indicates that alpha(2)beta(1) integrin cannot simultaneously bind both EV1 and the physiological ligand collagen. Compared with collagen binding to the alpha(2)I domain, the virus binds with a 10-fold higher affinity but in vitro uncoating of EV1 was not observed as a result of attachment of alpha(2)I. A molecular model, constructed on the basis of the EV1-integrin complex, shows that multiple alpha(2)beta(1) heterodimers can bind at adjacent sites around the virus 5-fold symmetry axes without steric hindrance. In agreement with this, virus attachment to alpha(2)beta(1) integrin on the cell surface was found to result in integrin clustering, which can give rise to signaling and facilitate the initiation of the viral entry process that takes place via caveolae-mediated endocytosis.

Structural basis for PAS domain heterodimerization in the basic helix--loop--helix-PAS transcription factor hypoxia-inducible factor

Biological responses to oxygen availability play important roles in development, physiological homeostasis, and many disease processes. In mammalian cells, this adaptation is mediated in part by a conserved pathway centered on the hypoxia-inducible factor (HIF). HIF is a heterodimeric protein complex composed of two members of the basic helix-loop-helix Per-ARNT-Sim (PAS) (ARNT, aryl hydrocarbon receptor nuclear translocator) domain family of transcriptional activators, HIFalpha and ARNT. Although this complex involves protein-protein interactions mediated by basic helix-loop-helix and PAS domains in both proteins, the role played by the PAS domains is poorly understood. To address this issue, we have studied the structure and interactions of the C-terminal PAS domain of human HIF-2alpha by NMR spectroscopy. We demonstrate that HIF-2alpha PAS-B binds the analogous ARNT domain in vitro, showing that residues involved in this interaction are located on the solvent-exposed side of the HIF-2alpha central beta-sheet. Mutating residues at this surface not only disrupts the interaction between isolated PAS domains in vitro but also interferes with the ability of full-length HIF to respond to hypoxia in living cells. Extending our findings to other PAS domains, we find that this beta-sheet interface is widely used for both intra- and intermolecular interactions, suggesting a basis of specificity and regulation of many types of PAS-containing signaling proteins.

Architecture and selectivity in aquaporins: 2.5 a X-ray structure of aquaporin Z

Aquaporins are a family of water and small molecule channels found in organisms ranging from bacteria to animals. One of these channels, the E. coli protein aquaporin Z (AqpZ), has been shown to selectively conduct only water at high rates. We have expressed, purified, crystallized, and solved the X-ray structure of AqpZ. The 2.5 Å resolution structure of AqpZ suggests aquaporin selectivity results both from a steric mechanism due to pore size and from specific amino acid substitutions that regulate the preference for a hydrophobic or hydrophilic substrate. This structure provides direct evidence on the molecular mechanisms of specificity between water and glycerol in this family of channels from a single species. It is to our knowledge the first atomic resolution structure of a recombinant aquaporin and so provides a platform for combined genetic, mutational, functional, and structural determinations of the mechanisms of aquaporins and, more generally, the assembly of multimeric membrane proteins.

The 2.5 Å resolution structure of the E. coli water- specific aquaporin AqpZ provides direct evidence for the molecular mechanisms of specificity for water, as opposed to glycerol, in this family of channels

Ligand-supported homology modelling of protein binding-sites using knowledge-based potentials

A new approach, MOBILE, is presented that models protein binding-sites including bound ligand molecules as restraints. Initially generated, homology models of the target protein are refined iteratively by including information about bioactive ligands as spatial restraints and optimising the mutual interactions between the ligands and the binding-sites. Thus optimised models can be used for structure-based drug design and virtual screening. In a first step, ligands are docked into an averaged ensemble of crude homology models of the target protein. In the next step, improved homology models are generated, considering explicitly the previously placed ligands by defining restraints between protein and ligand atoms. These restraints are expressed in terms of knowledge-based distance-dependent pair potentials, which were compiled from crystallographically determined protein-ligand complexes. Subsequently, the most favourable models are selected by ranking the interactions between the ligands and the generated pockets using these potentials. Final models are obtained by selecting the best-ranked side-chain conformers from various models, followed by an energy optimisation of the entire complex using a common force-field. Application of the knowledge-based pair potentials proved efficient to restrain the homology modelling process and to score and optimise the modelled protein-ligand complexes. For a test set of 46 protein-ligand complexes, taken from the Protein Data Bank (PDB), the success rate of producing near-native binding-site geometries (rmsd<2.0A) with MODELLER is 70% when the ligand restrains the homology modelling process in its native orientation. Scoring these complexes with the knowledge-based potentials, in 66% of the cases a pose with rmsd <2.0A is found on rank 1. Finally, MOBILE has been applied to two case studies modelling factor Xa based on trypsin and aldose reductase based on aldehyde reductase.

The quaternary structure of the Saccharomyces cerevisiae succinate dehydrogenase. Homology modeling, cofactor docking, and molecular dynamics simulation studies

Succinate dehydrogenases and fumarate reductases are complex mitochondrial or bacterial respiratory chain proteins with remarkably similar structures and functions. Succinate dehydrogenase oxidizes succinate and reduces ubiquinone using a flavin adenine dinucleotide cofactor and iron-sulfur clusters to transport electrons. A model of the quaternary structure of the tetrameric Saccharomyces cerevisiae succinate dehydrogenase was constructed based on the crystal structures of the Escherichia coli succinate dehydrogenase, the E. coli fumarate reductase, and the Wolinella succinogenes fumarate reductase. One FAD and three iron-sulfur clusters were docked into the Sdh1p and Sdh2p catalytic dimer. One b-type heme and two ubiquinone or inhibitor analog molecules were docked into the Sdh3p and Sdh4p membrane dimer. The model is consistent with numerous experimental observations. The calculated free energies of inhibitor binding are in excellent agreement with the experimentally determined inhibitory constants. Functionally important residues identified by mutagenesis of the SDH3 and SDH4 genes are located near the two proposed quinone-binding sites, which are separated by the heme. The proximal quinone-binding site, located nearest the catalytic dimer, has a considerably more polar environment than the distal site. Alternative low energy conformations of the membrane subunits were explored in a molecular dynamics simulation of the dimer embedded in a phospholipid bilayer. The simulation offers insight into why Sdh4p Cys-78 may be serving as the second axial ligand for the heme instead of a histidine residue. We discuss the possible roles of heme and of the two quinone-binding sites in electron transport.

The export of coat protein from enteroaggregative Escherichia coli by a specific ATP-binding cassette transporter system

Enteroaggregative Escherichia coli (EAEC) is an emerging enteric pathogen characterized by aggregative adherence (AA) to cultured human mucosal epithelium cells. We have recently characterized a 10.2-kDa protein, called dispersin, which is exported from the bacteria and which promotes dispersal of EAEC across the intestinal mucosa. Here, we present evidence that dispersin is exported by a putative ABC transporter complex, which is encoded by a genetic locus of the EAEC virulence plasmid pAA2. We demonstrate that the locus comprises a cluster of five genes (designated aat-PABCD), including homologs of an inner-membrane permease (AatP), an ATP-binding cassette protein (AatC) and the outer membrane protein TolC (AatA). We show that, like TolC, AatA localizes to the outer membrane independently of its ABC partner. Dispersin appears to require the Aat complex for outer membrane translocation but not for secretion across the inner membrane. We also show that, like the dispersin gene, transcription of the aat cluster is dependent on AggR, a regulator of virulence genes in EAEC. We propose that the aat cluster encodes a specialized ABC transporter, which plays a role in the pathogenesis of EAEC by transporting dispersin out of the bacterial cell.

OXBench: a benchmark for evaluation of protein multiple sequence alignment accuracy

Background

The alignment of two or more protein sequences provides a powerful guide in the prediction of the protein structure and in identifying key functional residues, however, the utility of any prediction is completely dependent on the accuracy of the alignment. In this paper we describe a suite of reference alignments derived from the comparison of protein three-dimensional structures together with evaluation measures and software that allow automatically generated alignments to be benchmarked. We test the OXBench benchmark suite on alignments generated by the AMPS multiple alignment method, then apply the suite to compare eight different multiple alignment algorithms. The benchmark shows the current state-of-the art for alignment accuracy and provides a baseline against which new alignment algorithms may be judged.

Results

The simple hierarchical multiple alignment algorithm, AMPS, performed as well as or better than more modern methods such as CLUSTALW once the PAM250 pair-score matrix was replaced by a BLOSUM series matrix. AMPS gave an accuracy in Structurally Conserved Regions (SCRs) of 89.9% over a set of 672 alignments. The T-COFFEE method on a data set of families with <8 sequences gave 91.4% accuracy, significantly better than CLUSTALW (88.9%) and all other methods considered here. The complete suite is available from .

Conclusions

The OXBench suite of reference alignments, evaluation software and results database provide a convenient method to assess progress in sequence alignment techniques. Evaluation measures that were dependent on comparison to a reference alignment were found to give good discrimination between methods. The STAMP S_c Score which is independent of a reference alignment also gave good discrimination. Application of OXBench in this paper shows that with the exception of T-COFFEE, the majority of the improvement in alignment accuracy seen since 1985 stems from improved pair-score matrices rather than algorithmic refinements. The maximum theoretical alignment accuracy obtained by pooling results over all methods was 94.5% with 52.5% accuracy for alignments in the 0–10 percentage identity range. This suggests that further improvements in accuracy will be possible in the future.

Cosolvent-Induced Aggregation Inhibits Myosin ATPase Activity by Stabilizing the Predominant Transition Intermediate

High concentration of the cosolvent poly(ethylene glycol) (PEG) induces reversible aggregation of skeletal myosin subfragment 1 (S1) and inhibition of its Mg-ATPase activity [Highsmith et al. (1998) Biophys. J. 74, 1465-1472]. In the present work the effect of aggregation on the various steps of the ATPase cycle was studied. The isomerization and hydrolysis steps of the cycle were not affected by S1 aggregation since the formation of the "trapped" S1.MgADP.phosphate analogue complexes, which mimic the prehydrolysis M*ATP and posthydrolysis M**ADP.P(i) transition states, proceeded without any hindrance. Similar conclusions could be reached from the chemical modification of Lys-83 and Cys-707 in the presence of MgATP and MgATPgammaS, which indicated that the most populated intermediate of the cycle in solubilized and aggregated S1 is M**ADP.P(i). The dissociation of the trapped S1.MgADP.phosphate analogue complexes resembling the M**ADP.P(i) state was strongly inhibited by PEG-6000, showing that the transition from this intermediate is prevented by the aggregation. This step is presumably inhibited because the coupled swinging of the lever arm from the closed to the open position is constrained by the close packing of aggregated S1.

Intramolecular interactions of the regulatory region with the catalytic core in the plasma membrane calcium pump

The access of three proteases to their sites of cleavage was used as a measure of regulatory interactions in the plasma membrane Ca2+ pump isoform 4b (PMCA4b). When the proteases could not cut at their sites in the C-terminal regulatory region, the interaction was judged to be tight. This was the case in the absence of Ca2+, when chymotrypsin and caspase cut PMCA only very slowly. Ca2+ accelerated the fragmentation, but the digestion remained incomplete. In the presence of Ca2+ plus calmodulin, the digestion became nearly complete in all cases, indicating a more flexible conformation of the carboxyl terminus in the fully activated state. The acceleration of proteolysis by Ca2+ or Ca2+ plus calmodulin occurred equally at the caspase site upstream of the calmodulin-binding domain and the chymotrypsin and calpain sites downstream of that domain. Replacing Trp1093 (a key residue within the calmodulin-binding domain) with alanine had a much more specific effect, because it exposed only proteolytic sites within the calmodulin-binding domain that had previously been shielded in the native protein. At these sites, both calpain and chymotrypsin cut the Trp1093 --> Ala mutant in the absence of calmodulin. These data indicate that, in the auto-inhibited conformation, the calmodulin-binding/auto-inhibitory sequence and the regions both upstream and downstream are in close contact with the catalytic core. Trp1093 plays an essential role not only in stabilizing the Ca2+-calmodulin/calmodulin-binding domain complex but also in the formation or stability of the inhibitory conformation of that domain when it interacts with the catalytic core of PMCA4b.

Comparative modeling of GABA(A) receptors: limits, insights, future developments

GABA(A) receptors are chloride ion channels that mediate fast synaptic transmission and belong to a superfamily of pentameric ligand-gated ion channels. The recently published crystal structure of the acetylcholine binding protein can be used as a template for comparative modeling of the extracellular domain of GABA(A) receptors. In this commentary, difficulties with comparative modeling at low sequence identity are discussed, the degree of structural conservation to be expected within the superfamily is analyzed and numerical estimates of model uncertainties in functional regions are provided. Topography of the binding sites at subunit-interfaces is examined and possible targets for rational mutagenesis studies are suggested. Allosteric motions are considered and a mechanism for mediation of positive cooperativity at the benzodiazepine site is proposed.

Potential for dramatic improvement in sequence alignment against structures of remote homologous proteins by extracting structural information from multiple structure alignment

A novel method has been developed for acquiring the correct alignment of a query sequence against remotely homologous proteins by extracting structural information from profiles of multiple structure alignment. A systematic search algorithm combined with a group of score functions based on sequence information and structural information has been introduced in this procedure. A limited number of top solutions (15,000) with high scores were selected as candidates for further examination. On a test-set comprising 301 proteins from 75 protein families with sequence identity less than 30%, the proportion of proteins with completely correct alignment as first candidate was improved to 39.8% by our method, whereas the typical performance of existing sequence-based alignment methods was only between 16.1% and 22.7%. Furthermore, multiple candidates for possible alignment were provided in our approach, which dramatically increased the possibility of finding correct alignment, such that completely correct alignments were found amongst the top-ranked 1000 candidates in 88.3% of the proteins. With the assistance of a sequence database, completely correct alignment solutions were achieved amongst the top 1000 candidates in 94.3% of the proteins. From such a limited number of candidates, it would become possible to identify more correct alignment using a more time-consuming but more powerful method with more detailed structural information, such as side-chain packing and energy minimization, etc. The results indicate that the novel alignment strategy could be helpful for extending the application of highly reliable methods for fold identification and homology modeling to a huge number of homologous proteins of low sequence similarity. Details of the methods, together with the results and implications for future development are presented.

Structural basis of the lisinopril-binding specificity in N- and C-domains of human somatic ACE

Angiotensin I-converting enzyme (ACE) is a dipeptidyl carboxypeptidase which converts angiotensin I into the vasopressor peptide angiotensin II and also inactivates the hypotensive peptide bradykinin, playing an important role in blood pressure regulation. The present work describes the molecular modeling of the N-terminal human somatic ACE in complex with the inhibitor lisinopril, identifying the residues involved in the inhibitor-binding pocket. The obtained results identify differences in the lisinopril lysine moiety-binding residues for N- and C-terminals of sACE domains and an important carboxy-terminal proline hydrophobic accommodations mediated by the aromatic ring of Tyr532 and Tyr1128 residues, respectively. The present model will be useful for the development of a new inhibitor family based on the natural BPP peptides and derivatives, or even to improve the binding capacities and the domain specificity of the already known inhibitors.

A structural model for the open conformation of the mdr1 P-glycoprotein based on the MsbA crystal structure

The validity of the structure of the Escherichia coli MsbA lipid transporter as a model from the mdr1 P-glycoprotein has been evaluated. Comparative sequence analyses, motif search and secondary structure prediction indicated that each of the two P-glycoprotein halves is structurally similar to the MsbA monomer and also suggested that the open dimer structure is valid for P-glycoprotein. Homology modeling was used to predict the structure of P-glycoprotein using MsbA as a template. The resulting modeled structure allowed a detailed study of the interactions between the intracellular domain and the nucleotide binding domain and suggested that these contacts are involved in mediating the coupling between nucleotide binding domain conformational changes and transmembrane helices reorientation during transport. In P-glycoprotein, the internal chamber open to the inner leaflet and the inner medium is significantly different in size and charge than in MsbA. These differences can be related to those of the transported substrates. Moreover an ensemble of 20 conserved aromatic residues appears to border the periphery of each side of the chamber in P-glycoprotein. These may be important for size selection and proper positioning of drugs for transport. The relevance of the modeled conformation to P-gp function is discussed.

Detecting potential IgE-reactive sites on food proteins using a sequence and structure database, SDAP-food

The high incidence of food allergies, including oral allergy syndrome, represent major considerations when introducing new crops and foods. A new structural database of allergenic proteins, SDAP-Food, http://fermi.utmb.edu/SDAP/, has been developed to aid in predicting the IgE-binding potential of novel food proteins and cross-reactivities among known allergens. The site is designed to facilitate the first steps of a decision tree approach to determine the allergenicity of a given protein, based on the sequence and structural similarity to known allergens and their IgE binding sites. Immunological tests can then be used to confirm the predictions. A hierarchical procedure for identifying potential allergens, using a physical property-based sequence similarity index, has been designed to identify regions that resemble known IgE binding sites. As an example, SDAP tools were used to find food allergen sequences similar to an IgE binding site of the Jun a 3 allergen from mountain cedar pollen. The SDAP sequence similarity search matched the Jun a 3 epitope to regions in several food allergens, including cherry (Pru av 2), apple (Mal d 2) and pepper (Cap a 1), which are, like Jun a 3, members of the plant pathogenesis-related (PR-5) protein family. Homology modeling, using our EXDIS/DIAMOD/FANTOM program suite, indicated a similar surface location and structure for the potential epitope region on all of these allergens. The quantitative approach presented here can be used as part of a screening process for potential allergenicity of recombinant food products.

Comparative protein structure modeling by iterative alignment, model building and model assessment

Comparative or homology protein structure modeling is severely limited by errors in the alignment of a modeled sequence with related proteins of known three-dimensional structure. To ameliorate this problem, we have developed an automated method that optimizes both the alignment and the model implied by it. This task is achieved by a genetic algorithm protocol that starts with a set of initial alignments and then iterates through re-alignment, model building and model assessment to optimize a model assessment score. During this iterative process: (i) new alignments are constructed by application of a number of operators, such as alignment mutations and cross-overs; (ii) comparative models corresponding to these alignments are built by satisfaction of spatial restraints, as implemented in our program MODELLER; (iii) the models are assessed by a variety of criteria, partly depending on an atomic statistical potential. When testing the procedure on a very difficult set of 19 modeling targets sharing only 4-27% sequence identity with their template structures, the average final alignment accuracy increased from 37 to 45% relative to the initial alignment (the alignment accuracy was measured as the percentage of positions in the tested alignment that were identical to the reference structure-based alignment). Correspondingly, the average model accuracy increased from 43 to 54% (the model accuracy was measured as the percentage of the C(alpha) atoms of the model that were within 5 A of the corresponding C(alpha) atoms in the superposed native structure). The present method also compares favorably with two of the most successful previously described methods, PSI-BLAST and SAM. The accuracy of the final models would be increased further if a better method for ranking of the models were available.

Mutational and structural analysis of aglycone specificity in maize and sorghum beta-glucosidases

Plant beta-glucosidases display varying substrate specificities. The maize beta-glucosidase isozyme Glu1 (ZmGlu1) hydrolyzes a broad spectrum of substrates in addition to its natural substrate DIMBOA-Glc (2-O-beta-d-glucopyranosyl-4-hydroxy-7-methoxy-1,4-benzoxaxin-3-one), whereas the sorghum beta-glucosidase isozyme Dhr1 (SbDhr1) hydrolyzes exclusively its natural substrate dhurrin (p-hydroxy-(S)-mandelonitrile-beta-d-glucoside). Structural data from cocrystals of enzyme-substrate and enzyme-aglycone complexes have shown that five amino acid residues (Phe198, Phe205, Trp378, Phe466, and Ala467) are located in the aglycone-binding site of ZmGlu1 and form the basis of aglycone recognition and binding, hence substrate specificity. To study the mechanism of substrate specificity further, mutant beta-glucosidases were generated by replacing Phe198, Phe205, Asp261, Met263, Phe377, Phe466, Ala467, and Phe473 of Glu1 by Dhr1 counterparts. The effects of mutations on enzyme activity and substrate specificity were studied using both natural and artificial substrates. The simple mutant replacing Phe198 by a valine had the most drastic effect on activity, because the capacity of this enzyme to hydrolyze beta-glucosides was almost completely abolished. The analysis of this mutation was completed by a structural study of the double mutant ZmGlu1-E191D,F198V in complex with the natural substrate. The structure reveals that the single mutation F198V causes a cascade of conformational changes, which are unpredictable by standard molecular modeling techniques. Some other mutations led to drastic effects: replacing Asp261 by an asparagine decreases the catalytic efficiency of this simple mutant by 75% although replacing Tyr473 by a phenylalanine increase its efficiency by 300% and also provides a new substrate specificity by hydrolyzing dhurrin.

Development of natural product-derived receptor tyrosine kinase inhibitors based on conservation of protein domain fold

Receptor tyrosine kinases (RTKs) such as Tie-2, IGF1R, Her-2/Neu, EGFR, and VEGFR1-3 play crucial roles in the control of cell growth and differentiation. Inhibition of such RTKs has become a major focus of current anticancer drug development, and therefore the discovery of new classes of inhibitors for these signal-transducing proteins is of prime importance. We have recently proposed a novel concept for improving the hit-finding process by employing natural products as biologically validated starting points in structural space for compound library development. In this concept, natural products are regarded as evolutionary chosen ligands for protein domains which are structurally conserved yet genetically mobile. Here we report on the discovery of novel and highly selective VEGFR-2 and -3, Tie-2, and IGF1R inhibitors derived from the naturally occurring Her-2/Neu kinase inhibitor nakijiquinone C and developed on the basis of this concept. Based on the structure of the natural product, a small library (74 members) was synthesized and investigated for inhibition of kinases with highly similar ATP-binding domains. The library yielded inhibitors with IC(50)s in the low micromolar range with high frequency (7 out of 74). In particular, four inhibitors of Tie-2 were found, a kinase critically involved in the formation of new blood vessels from preexisting ones (angiogenesis) and believed to be a new promising target in antitumor therapy. These results support the "domain concept". To advance the development of improved inhibitors, extensive molecular modeling studies were undertaken, including the construction of new homology models for VEGFR-2 and Tie-2. These studies revealed residues in the kinase structure which are crucial to the development of tailor-made receptor tyrosine kinase inhibitors.

Mechanisms of thermal adaptation revealed from the genomes of the Antarctic Archaea Methanogenium frigidum and Methanococcoides burtonii

We generated draft genome sequences for two cold-adapted Archaea, Methanogenium frigidum and Methanococcoides burtonii, to identify genotypic characteristics that distinguish them from Archaea with a higher optimal growth temperature (OGT). Comparative genomics revealed trends in amino acid and tRNA composition, and structural features of proteins. Proteins from the cold-adapted Archaea are characterized by a higher content of noncharged polar amino acids, particularly Gln and Thr and a lower content of hydrophobic amino acids, particularly Leu. Sequence data from nine methanogen genomes (OGT 15 degrees -98 degrees C) were used to generate 1111 modeled protein structures. Analysis of the models from the cold-adapted Archaea showed a strong tendency in the solvent-accessible area for more Gln, Thr, and hydrophobic residues and fewer charged residues. A cold shock domain (CSD) protein (CspA homolog) was identified in M. frigidum, two hypothetical proteins with CSD-folds in M. burtonii, and a unique winged helix DNA-binding domain protein in M. burtonii. This suggests that these types of nucleic acid binding proteins have a critical role in cold-adapted Archaea. Structural analysis of tRNA sequences from the Archaea indicated that GC content is the major factor influencing tRNA stability in hyperthermophiles, but not in the psychrophiles, mesophiles or moderate thermophiles. Below an OGT of 60 degrees C, the GC content in tRNA was largely unchanged, indicating that any requirement for flexibility of tRNA in psychrophiles is mediated by other means. This is the first time that comparisons have been performed with genome data from Archaea spanning the growth temperature extremes from psychrophiles to hyperthermophiles.

SWISS-MODEL: An automated protein homology-modeling server

SWISS-MODEL (http://swissmodel.expasy.org) is a server for automated comparative modeling of three-dimensional (3D) protein structures. It pioneered the field of automated modeling starting in 1993 and is the most widely-used free web-based automated modeling facility today. In 2002 the server computed 120 000 user requests for 3D protein models. SWISS-MODEL provides several levels of user interaction through its World Wide Web interface: in the 'first approach mode' only an amino acid sequence of a protein is submitted to build a 3D model. Template selection, alignment and model building are done completely automated by the server. In the 'alignment mode', the modeling process is based on a user-defined target-template alignment. Complex modeling tasks can be handled with the 'project mode' using DeepView (Swiss-PdbViewer), an integrated sequence-to-structure workbench. All models are sent back via email with a detailed modeling report. WhatCheck analyses and ANOLEA evaluations are provided optionally. The reliability of SWISS-MODEL is continuously evaluated in the EVA-CM project. The SWISS-MODEL server is under constant development to improve the successful implementation of expert knowledge into an easy-to-use server.

Tools for comparative protein structure modeling and analysis

The following resources for comparative protein structure modeling and analysis are described (http://salilab.org): MODELLER, a program for comparative modeling by satisfaction of spatial restraints; MODWEB, a web server for automated comparative modeling that relies on PSI-BLAST, IMPALA and MODELLER; MODLOOP, a web server for automated loop modeling that relies on MODELLER; MOULDER, a CPU intensive protocol of MODWEB for building comparative models based on distant known structures; MODBASE, a comprehensive database of annotated comparative models for all sequences detectably related to a known structure; MODVIEW, a Netscape plugin for Linux that integrates viewing of multiple sequences and structures; and SNPWEB, a web server for structure-based prediction of the functional impact of a single amino acid substitution.

Study of the structural dynamics of the E coli 70S ribosome using real-space refinement

Cryo-EM density maps showing the 70S ribosome of E. coli in two different functional states related by a ratchet-like motion were analyzed using real-space refinement. Comparison of the two resulting atomic models shows that the ribosome changes from a compact structure to a looser one, coupled with the rearrangement of many of the proteins. Furthermore, in contrast to the unchanged inter-subunit bridges formed wholly by RNA, the bridges involving proteins undergo large conformational changes following the ratchet-like motion, suggesting an important role of ribosomal proteins in facilitating the dynamics of translation.

Tetratricopeptide-like repeats in type-III-secretion chaperones and regulators

Efficient type-III secretion depends on cytosolic molecular chaperones, which bind specifically to the translocators and effectors. In the past there has been a tendency to shoe-horn all type-III-secretion chaperones into a single structural and functional class. However, we have shown that the LcrH/SycD-like chaperones consist of three central tetratricopeptide-like repeats that are predicted to fold into an all-alpha-helical array that is quite distinct from the known structure of the SycE class of chaperones. Furthermore, we predict that this array creates a peptide-binding groove that is utterly different from the helix-binding groove in SycE. We present a homology model of LcrH/SycD that is consistent with existing mutagenesis data. We also report the existence of tetratricopeptide-like repeats in regulators of type-III secretion, such as HilA from Salmonella enterica and HrpB from Ralstonia solanacearum. The discovery of tetratricopeptide-like repeats in type-III-secretion regulators and chaperones provides a new conceptual framework for structural and mutagenesis studies and signals a potential unification of prokaryotic and eukaryotic chaperone biology.

Generation of accurate protein loop conformations through low-barrier molecular dynamics

Prediction and refinement of protein loop structures are important and challenging tasks for which no general solution has been found. In addition to the accuracy of scoring functions, the main problems reside in (1) insufficient statistical sampling and (2) crossing energy barriers that impede conformational rearrangements of the loop. We approach these two issues by using "low-barrier molecular dynamics," a combination of energy smoothing techniques. To address statistical sampling, locally enhanced sampling (LES) is used to produce multiple copies of the loop, thus improving statistics and reducing energy barriers. We introduce a novel extension of LES that can improve local sampling even further through hierarchical subdivision of copies. Even though LES reduces energy barriers, it cannot provide for crossing infinite barriers, which can be problematic when substantial rearrangement of residues is necessary. To permit this kind of loop residue repacking, a "soft-core" potential energy function is introduced, so that atomic overlaps are temporarily allowed. We tested this new combined methodology to a loop in anti-influenza antibody Fab 17/9 (7 residues long) and to another loop in the antiprogesterone antibody DB3 (8 residues). In both cases, starting from random conformations, we were able to locate correct loop structures (including sidechain orientations) with heavy-atom root-mean-square deviation (fit to the nonloop region) of approximately 1.1 A in Fab 17/9 and approximately 1.8 A in DB3. We show that the combination of LES and soft-core potential substantially improves sampling compared to regular molecular dynamics. Moreover, the sampling improvement obtained with this combined approach is significantly better than that provided by either of the two methods alone.

A structural perspective on genome evolution

Protein translations of over 100 complete genomes are now available. About half of these sequences can be provided with structural annotation, thereby enabling some profound insights into protein and pathway evolution. Whereas the major domain structure families are common to all kingdoms of life, these are combined in different ways in multidomain proteins to give various domain architectures that are specific to kingdoms or individual genomes, and contribute to the diverse phenotypes observed. These data argue for more targets in structural genomics initiatives and particularly for the selection of different domain architectures to gain better insights into protein functions.

Genetic strategies for antibacterial drug discovery

The availability of genome sequences is revolutionizing the field of microbiology. Genetic methods are being modified to facilitate rapid analysis at a genome-wide level and are blossoming for human pathogens that were previously considered intractable. This revolution coincided with a growing concern about the emergence of microbial drug resistance, compelling the pharmaceutical industry to search for new antimicrobial agents. The availability of the new technologies, combined with many genetic strategies, has changed the way that researchers approach antibacterial drug discovery.

Molecular determinants of proton-sensitive N-methyl-D-aspartate receptor gating

Extracellular protons inhibit N-methyl-D-aspartate (NMDA) receptors with an IC50 value in the physiological pH range. To identify the molecular determinants of proton sensitivity, we used scanning mutagenesis of the NR1 subunit to search for residues that control proton inhibition of NMDA receptors. Homology modeling of the extracellular domains suggested that residues at which mutations perturbed pH sensitivity were localized in discrete regions. The majority of mutations that strongly affected proton sensitivity were clustered in the extracellular end of the second transmembrane domain (M3) and adjacent linker leading to the S2 portion of the glycine-binding domain of NR1. Mutations in NR2A confirmed that the analogous region controls the pH sensitivity of this subunit and also identified the linker region between the third transmembrane domain (M4) and the S2 portion of the NR2 glutamate binding domain as an additional determinant of proton sensitivity. One mutant receptor, NR1(A649C)/NR2A(A651T), showed a 145-fold reduction in the IC50 for protons (IC50, 17.3 microM corresponding to pH 4.9). The M3-S2 linker region has been suggested to control NMDA receptor gating, leading to the hypothesis that the proton sensor and receptor gate may be structurally and functionally integrated.

In silico protein recombination: enhancing template and sequence alignment selection for comparative protein modelling

Comparative modelling of proteins is a predictive technique to build an atomic model for a given amino acid sequence, on the basis of the structures of other proteins (templates) that have been determined experimentally. Critical problems arise in this procedure: selecting the correct templates, aligning the query sequence with them and building the non-conserved surface loops. In this work, we apply a genetic algorithm, with crossover and mutation, as a new tool to overcome the first two. In silico protein recombination proves to be an effective way to exploit the variability of templates and sequence alignments to produce populations of optimized models by artificial selection. Despite some limitations, the procedure is shown to be robust to alignment errors, while simplifying the task of selecting templates, making it a good candidate for automatic building of reliable protein models.

Chemometrical classification of ephrin ligands and Eph kinases using GRID/CPCA approach

Eph receptor tyrosine kinases are divided on two subfamilies based on their affinity for ephrin ligands and play a crucial role in the intercellular processes such as angiogenesis, neurogenesis, and carcinogenesis. As such, Eph kinases represent potential targets for drug design, which requires the knowledge of structural features responsible for their specific interactions. To overcome the existing gap between available sequence and structure information we have built 3D models of eight ephrins and 13 Eph kinase ligand-binding domains using homology modeling techniques. The interaction energies for several molecular probes with binding sites of these models were calculated using GRID and subjected to chemometrical classification based on consensus principal component analysis (CPCA). Despite inherent limitations of the homology models, CPCA was able to successfully distinguish between ephrins and Eph kinases, between Eph kinase subfamilies, and between ephrin subfamilies. As a result we have identified several amino acids that may account for selectivity in ephrin-Eph kinase interactions. In general, although the difference in charge between ephrin and Eph kinase binding domains creates an attractive long-range electrostatic force, the hydrophobic and steric interactions are highly important for the short-range interactions between two proteins. The chemometrical analysis also provides the pharmacophore model, which could be used for virtual screening and de novo ligand design.

Mechanism of proton pumping by plant plasma membrane H+-ATPase: role of residues in transmembrane segments 5 and 6

The mechanism of proton pumping by P-type plasma membrane H(+)-ATPases is not well clarified. Site-directed mutagenesis studies suggest that Asp684, situated in transmembrane segment M6, is involved in coordination of proton(s) in plant plasma membrane H(+)-ATPase. This hypothesis is supported by atomic models of H(+)-ATPases built on the basis of the crystal structure of the related SERCA1a Ca(2+)-ATPase. However, more biochemical, genetic, and structural studies are required before we will be able to understand the nature of the proton binding site(s) in P-type H(+)-ATPases and the mechanism of action of these pumps.

Bioinorganic chemistry in the postgenomic era

Genome sequencing has revolutionized all fields of life sciences. Bioinorganic chemistry is certainly not immune to this influence, which is presenting unprecedented challenges. A new goal for bioinorganic chemistry is the investigation of the linkages between inorganic elements and genomic information. This requires new advancements andor the development of new expertise in fields such as bioinformatics and genetics but also provides a driving force to push forward the exploitation of traditional analytical techniques and spectroscopic tools. The "case study" of metal homeostasis in cells is discussed to provide a flavor of the current evolution of the field.

Crystal structure of Drosophila angiotensin I-converting enzyme bound to captopril and lisinopril

Angiotensin I-converting enzymes (ACEs) are zinc metallopeptidases that cleave carboxy-terminal dipeptides from short peptide hormones. We have determined the crystal structures of AnCE, a Drosophila homolog of ACE, with and without bound inhibitors to 2.4 A resolution. AnCE contains a large internal channel encompassing the entire protein molecule. This substrate-binding channel is composed of two chambers, reminiscent of a peanut shell. The inhibitor and zinc-binding sites are located in the narrow bottleneck connecting the two chambers. The substrate and inhibitor specificity of AnCE appears to be determined by extensive hydrogen-bonding networks and ionic interactions in the active site channel. The catalytically important zinc ion is coordinated by the conserved Glu395 and histidine residues from a HExxH motif.

From words to literature in structural proteomics

Technical advances on several frontiers have expanded the applicability of existing methods in structural biology and helped close the resolution gaps between them. As a result, we are now poised to integrate structural information gathered at multiple levels of the biological hierarchy - from atoms to cells - into a common framework. The goal is a comprehensive description of the multitude of interactions between molecular entities, which in turn is a prerequisite for the discovery of general structural principles that underlie all cellular processes.

Structural genomics of proteins involved in copper homeostasis

Genome sequencing projects have provided a wealth of data, most notably the primary sequences of all the proteins that a given organism can produce. The understanding of this information at the functional level is still in the beginning stages. Three-dimensional structural information is necessary to unravel at the atomic level the mechanisms by which a protein carries out its function, and such information can often be very useful to predict at least gross functional features, even in the absence of biochemical data. An exhaustive structural characterization of the proteins encoded in the genomes is thus highly desirable. To enhance the functional insights provided by genome-scale structural determination, we have prioritized our research to target specific processes of the cell, i.e., those responsible for controlling metal homeostasis. In this Account, we present the results obtained by the Magnetic Resonance Center of the University of Florence on proteins involved in the homeostasis of copper. The general research strategy is presented, followed by a discussion focused on different key experimental aspects. An overview of the initial results and of their relevance to the understanding of molecular function and cellular processes is also given.

Interactions of substrate with calreticulin, an endoplasmic reticulum chaperone

Calreticulin is a molecular chaperone found in the endoplasmic reticulum in eukaryotes, and its interaction with N-glycosylated polypeptides is mediated by the glycan Glc(1)Man(7-9)GlcNAc(2) present on the target glycoproteins. Here, we report the thermodynamic parameters of its interaction with di-, tri-, and tetrasaccharide, which are truncated versions of the glucosylated arm of Glc(1)Man(7-9)GlcNAc(2), determined by the quantitative technique of isothermal titration calorimetry. This method provides a direct estimate of the binding constants (K(b)) and changes in enthalpy of binding (Delta H(b) degrees ) as well as the stoichiometry of the reaction. Unlike past speculations, these studies demonstrate unambiguously that calreticulin has only one site per molecule for binding its complementary glucosylated ligands. Although the binding of glucose by itself is not detectable, a binding constant of 4.19 x 10(4) m(-1) at 279 K is obtained when glucose occurs in alpha-1,3 linkage to Man alpha Me as in Glc alpha 1-3Man alpha Me. The binding constant increases by 25-fold from di- to trisaccharide and doubles from tri- to tetrasaccharide, demonstrating that the entire Glc alpha 1-3Man alpha 1-2Man alpha 1-2Man alpha Me structure of the oligosaccharide is recognized by calreticulin. The thermodynamic parameters thus obtained were supported by modeling studies, which showed that increased number of hydrogen bonds and van der Waals interactions occur as the size of the oligosaccharide is increased. Also, several novel findings about the recognition of saccharide ligands by calreticulin vis á vis legume lectins, which have the same fold as this chaperone, are discussed.

Conservation of structure and function among tyrosine recombinases: homology-based modeling of the lambda integrase core-binding domain

Tyrosine recombinases participate in diverse biological processes by catalyzing recombination between specific DNA sites. Although a conserved protein fold has been described for the catalytic (CAT) domains of five recombinases, structural relationships between their core-binding (CB) domains remain unclear. Despite differences in the specificity and affinity of core-type DNA recognition, a conserved binding mechanism is suggested by the shared two-domain motif in crystal structure models of the recombinases Cre, XerD and Flp. We have found additional evidence for conservation of the CB domain fold. Comparison of XerD and Cre crystal structures showed that their CB domains are closely related; the three central alpha-helices of these domains are superposable to within 1.44 A. A structure-based multiple sequence alignment containing 25 diverse CB domain sequences provided evidence for widespread conservation of both structural and functional elements in this fold. Based upon the Cre and XerD crystal structures, we employed homology modeling to construct a three-dimensional structure for the lambda integrase CB domain. The model provides a conceptual framework within which many previously identified, functionally important amino acid residues were investigated. In addition, the model predicts new residues that may participate in core-type DNA binding or dimerization, thereby providing hypotheses for future genetic and biochemical experiments.

Toxoplasma gondii: analysis of the active site insertion of its ferredoxin-NADP(+)-reductase by peptide-specific antibodies and homology-based modeling

Apicomplexan parasites possess an apicoplast-localized redox system consisting of a plant-type ferredoxin-NADP(+)-reductase (FNR) and its redox partner ferredoxin, a small [2Fe-2S] protein. We show here that several apicomplexan FNRs contain unique amino acid insertions of various lengths which are located in close proximity to the enzymatically important FAD and ferredoxin-binding sites of these proteins. Using the insertion of the Toxoplasma gondii reductase as an example we raised epitope-specific antibodies against an 11 amino acids long peptide predicted to be surface-exposed within this insertion. This peptide was found to be immunogenic when presented to the immune system as part of a carrier protein, but also in its natural structural context in the whole recombinant protein, implying that the epitope is surface-exposed. Three-dimensional modeling of T. gondii FNR based on the known 3D-structure of maize root FNR predicts that the overall structure of plant and apicomplexan FNRs are very similar and that the 11 amino acids are part of an alpha-helix, looping out of the molecule. Collectively, these data suggest that the insertion in T. gondii FNR does not affect the overall structure of the protein but may have an effect on the binding dynamics of FAD, NADP(+), and/or ferredoxin to FNR.

Species-specific protein sequence and fold optimizations

BACKGROUND

An organism's ability to adapt to its particular environmental niche is of fundamental importance to its survival and proliferation. In the largest study of its kind, we sought to identify and exploit the amino-acid signatures that make species-specific protein adaptation possible across 100 complete genomes.

RESULTS

Environmental niche was determined to be a significant factor in variability from correspondence analysis using the amino acid composition of over 360,000 predicted open reading frames (ORFs) from 17 archaea, 76 bacteria and 7 eukaryote complete genomes. Additionally, we found clusters of phylogenetically unrelated archaea and bacteria that share similar environments by amino acid composition clustering. Composition analyses of conservative, domain-based homology modeling suggested an enrichment of small hydrophobic residues Ala, Gly, Val and charged residues Asp, Glu, His and Arg across all genomes. However, larger aromatic residues Phe, Trp and Tyr are reduced in folds, and these results were not affected by low complexity biases. We derived two simple log-odds scoring functions from ORFs (CG) and folds (CF) for each of the complete genomes. CF achieved an average cross-validation success rate of 85 +/- 8% whereas the CG detected 73 +/- 9% species-specific sequences when competing against all other non-redundant CG. Continuously updated results are available at http://genome.mshri.on.ca.

CONCLUSION

Our analysis of amino acid compositions from the complete genomes provides stronger evidence for species-specific and environmental residue preferences in genomic sequences as well as in folds. Scoring functions derived from this work will be useful in future protein engineering experiments and possibly in identifying horizontal transfer events.

Structuring the universe of proteins

High-throughput sequencing of human genomes and those of important model organisms (mouse, Drosophila melanogaster, Caenorhabditis elegans, fungi, archaea) and bacterial pathogens has laid the foundation for another "big science" initiative in biology. Together, X-ray crystallographers, nuclear magnetic resonance (NMR) spectroscopists, and computational biologists are pursuing high-throughput structural studies aimed at developing a comprehensive three-dimensional view of the protein structure universe. The new science of structural genomics promises more than 10,000 experimental protein structures and millions of calculated homology models of related proteins. The evolutionary underpinnings and technological challenges of automating target selection, protein expression and purification, sample preparation, NMR and X-ray data measurement/analysis, homology modeling, and structure/function annotation are discussed in detail. An informative case study from one of the structural genomics centers funded by the National Institutes of Health and the National Institute of General Medical Sciences (NIH/NIGMS) demonstrates how this experimental/computational pipeline will reveal important links between form and function in biology and provide new insights into evolution and human health and disease.

A single amino acid difference between human and monkey interleukin (IL)-1beta dictates effective binding to soluble type II IL-1 receptor

Soluble type II interleukin (IL)-1 receptor (sIL1R-II) binds human IL-1beta with high affinity and neutralizes its activity. Recombinant sIL1R-II is considered a potentially useful anti-IL-1 therapeutic, and preclinical studies have been undertaken with this molecule in primates. To better understand the cytokine-receptor interactions occurring in this nonhuman context, monkey IL-1 and IL1R-II were cloned, and their binding abilities were examined in vitro. IL-1beta from cynomolgus monkey was capable of binding and activating the human type I IL-1 receptor. However, unlike human IL-1beta, it was unable to effectively bind and become neutralized by sIL1R-II. Human and cynomolgus IL-1beta proteins are 96% identical, differing by only six amino acids. Structural and mutational analysis revealed that the unique sIL1R-II binding ability of human IL-1beta is due to a single amino acid difference compared with monkey IL-1beta.

Structural analysis of the inhibition of Cdk4 and Cdk6 by p16(INK4a) through molecular dynamics simulations

Cyclin-dependent kinases 4, 6 and 2 (Cdk4/6/2), are proteins that lead progression through the G1-S transition, a step strictly regulated in the process of cell proliferation. The p16(INK4a) tumor suppressor, whose expression is inhibited in a high number of cancers, binds to Cdk4/6 and inhibits phosphorylation of the retinoblastoma protein, forcing cells to remain in the G1 phase and therefore, arresting cell division. Accordingly, the design of small compounds mimicking the inhibition of p16(INK4a) appears to be a promising way to treat cancer. In order to get some insight into the key interactions governing recognition between different cyclin-dependent kinases and the p16(INK4a) tumor suppressor, the present work reports the results of molecular dynamics simulations of both, the Cdk6-p16(INK4a) complex and the Cdk4-p16(INK4a) complex, respectively at 300 K. Most of the key interactions observed, were already anticipated in the analysis of the crystal structure of Cdk6-p16(INK4a). However, a few different features found out from the analysis of these calculations provide a better understanding of the role of the T-loop conformation, a fragment of Cdks, and the way the ATP binding-site is distorted upon binding of p16(INK4a).

Alternative splicing of a beta4 subunit proline-rich motif regulates voltage-dependent gating and toxin block of Cav2.1 Ca2+ channels

Ca2+ channel beta subunits modify alpha1 subunit gating properties through direct interactions with intracellular linker domains. In a previous report (Helton and Horne, 2002), we showed that alternative splicing of the beta4 subunit had alpha1 subunit subtype-specific effects on Ca2+ channel activation and fast inactivation. We extend these findings in the present report to include effects on slow inactivation and block by the peptide toxin omega-conotoxin (CTx)-MVIIC. N-terminal deletion and site-directed mutagenesis experiments revealed that the effects of alternative splicing on toxin block and all aspects of gating could be attributed to a proline-rich motif found within N-terminal beta4b amino acids 10-20. Interestingly, this motif is conserved within the third postsynaptic density-95 (PSD-95)/Discs large/zona occludens-1 domain of the distantly related membrane-associated guanylate kinase homolog, PSD-95. Sequence identity of approximately 30% made possible the building of beta4a and beta4b three-dimensional structural models using PSD-95 as the target sequence. The models (1) reveal that alternative splicing of the beta4 N terminus results in dramatic differences in surface charge distribution and (2) localize the proline-rich motif of beta4b to an extended arm structure that flanks what would be the equivalent of a highly modified PSD-95 carboxylate binding loop. Northern blot analysis revealed a markedly different pattern of distribution for beta4a versus beta4b in the human CNS. Whereas beta4a is distributed throughout evolutionarily older regions of the CNS, beta4b is concentrated heavily in the forebrain. These results raise interesting questions about the functional role that alternative splicing of the beta4 subunit has played in the evolution of complex neural networks.

Solution structure of human Mts1 (S100A4) as determined by NMR spectroscopy

Mts1 is a member of the S100 family of Ca2+-binding proteins and is implicated in promoting tumor progression and metastasis. To better understand the structure-function relationships of this protein and to begin characterizing its Ca2+-dependent interaction with protein binding targets, the three-dimensional structure of mts1 was determined in the apo state by NMR spectroscopy. As with other S100 protein family members, mts1 is a symmetric homodimer held together by noncovalent interactions between two helices from each subunit (helices 1, 4, 1', and 4') to form an X-type four-helix bundle. Each subunit of mts1 has two EF-hand Ca2+-binding domains: a pseudo-EF-hand (or S100-hand) and a typical EF-hand that are brought into proximity by a small two-stranded antiparallel beta-sheet. The S100-hand is formed by helices 1 and 2, and is similar in conformation to other members of the S100 family. In the typical EF-hand, the position of helix 3 is similar to that of another member of the S100 protein family, calcyclin (S100A6), and less like that of other S100 family members for which three-dimensional structures are available in the calcium-free state (e.g., S100B and S100A1). The differences in the position of helix 3 in the apo state of these four S100 proteins are likely due to variations in the amino acid sequence in the C-terminus of helix 4 and in loop 2 (the hinge region) and could potentially be used to subclassify the S100 protein family.

X-ray Structures of Threonine Aldolase Complexes: Structural Basis of Substrate Recognition†,‡

L-Threonine acetaldehyde-lyase (threonine aldolase, TA) is a pyridoxal-5'-phosphate-dependent (PLP) enzyme that catalyzes conversion of L-threonine or L-allo-threonine to glycine and acetaldehyde in a secondary glycine biosynthetic pathway. X-ray structures of Thermatoga maritima TA have been determined as the apo-enzyme at 1.8 A resolution and bound to substrate L-allo-threonine and product glycine at 1.9 and 2.0 A resolution, respectively. Despite low pairwise sequence identities, TA is a member of aspartate aminotransferase (AATase) fold family of PLP enzymes. The enzyme forms a 222 homotetramer with the PLP cofactor bound via a Schiff-base linkage to Lys199 within a domain interface. The structure reveals bound calcium and chloride ions that appear to contribute to catalysis and oligomerization, respectively. Although L-threonine and L-allo-threonine are substrates for T. maritima TA, enzymatic assays revealed a strong preference for L-allo-threonine. Structures of the external aldimines with substrate/product reveal a pair of histidines that may provide flexibility in substrate recognition. Variation in the threonine binding pocket may explain preferences for L-allo-threonine versus L-threonine among TA family members.