ATP binding proteins with different folds share a common ATP-binding structural motif

Expanding FTMap for Fragment-Based Identification of Pharmacophore Regions in Ligand Binding Sites

The knowledge of ligand binding hot spots and of the important interactions within such hot spots is crucial for the design of lead compounds in the early stages of structure-based drug discovery. The computational solvent mapping server FTMap can reliably identify binding hot spots as consensus clusters, free energy minima that bind a variety of organic probe molecules. However, in its current implementation, FTMap provides limited information on regions within the hot spots that tend to interact with specific pharmacophoric features of potential ligands. E-FTMap is a new server that expands on the original FTMap protocol. E-FTMap uses 119 organic probes, rather than the 16 in the original FTMap, to exhaustively map binding sites, and identifies pharmacophore features as atomic consensus sites where similar chemical groups bind. We validate E-FTMap against a set of 109 experimentally derived structures of fragment-lead pairs, finding that highly ranked pharmacophore features overlap with the corresponding atoms in both fragments and lead compounds. Additionally, comparisons of mapping results to ensembles of bound ligands reveal that pharmacophores generated with E-FTMap tend to sample highly conserved protein-ligand interactions. E-FTMap is available as a web server at https://eftmap.bu.edu.

How Ligands Interact with the Kinase Hinge

ATP-competitive kinase inhibitors form hydrogen bond interactions with the kinase hinge region at the adenine binding site. Thus, it is crucial to explore hinge-ligand recognition as part of a rational drug design strategy. Here, harnessing known ligand-bound kinase structures and experimental assay resources, we first created a kinase structure-assay database (KSAD) containing 2705 nM ligand-bound kinase complexes. Then, using KSAD, we systematically investigate hinge-ligand binding patterns using interaction fingerprints, thereby delineating 15 different hydrogen-bond interaction modes. We believe these results will be valuable for de novo drug design and/or scaffold hopping of kinase-targeted drugs.

On the evolution of protein-adenine binding

Proteins' interactions with ancient ligands may reveal how molecular recognition emerged and evolved. We explore how proteins recognize adenine: a planar rigid fragment found in the most common and ancient ligands. We have developed a computational pipeline that extracts protein-adenine complexes from the Protein Data Bank, structurally superimposes their adenine fragments, and detects the hydrogen bonds mediating the interaction. Our analysis extends the known motifs of protein-adenine interactions in the Watson-Crick edge of adenine and shows that all of adenine's edges may contribute to molecular recognition. We further show that, on the proteins' side, binding is often mediated by specific amino acid segments ("themes") that recur across different proteins, such that different proteins use the same themes when binding the same adenine-containing ligands. We identify numerous proteins that feature these themes and are thus likely to bind adenine-containing ligands. Our analysis suggests that adenine binding has emerged multiple times in evolution.

Anion induced conformational preference of Cα NN motif residues in functional proteins

Among different ligand binding motifs, anion binding C^α NN motif consisting of peptide backbone atoms of three consecutive residues are observed to be important for recognition of free anions, like sulphate or biphosphate and participate in different key functions. Here we study the interaction of sulphate and biphosphate with C^α NN motif present in different proteins. Instead of total protein, a peptide fragment has been studied keeping C^α NN motif flanked in between other residues. We use classical force field based molecular dynamics simulations to understand the stability of this motif. Our data indicate fluctuations in conformational preferences of the motif residues in absence of the anion. The anion gives stability to one of these conformations. However, the anion induced conformational preferences are highly sequence dependent and specific to the type of anion. In particular, the polar residues are more favourable compared to the other residues for recognising the anion.

Deciphering common recognition principles of nucleoside mono/di and tri-phosphates binding in diverse proteins via structural matching of their binding sites

Nucleoside triphosphate (NTP) ligands are of high biological importance and are essential for all life forms. A pre-requisite for them to participate in diverse biochemical processes is their recognition by diverse proteins. It is thus of great interest to understand the basis for such recognition in different proteins. Towards this, we have used a structural bioinformatics approach and analyze structures of 4677 NTP complexes available in Protein Data Bank (PDB). Binding sites were extracted and compared exhaustively using PocketMatch, a sensitive in-house site comparison algorithm, which resulted in grouping the entire dataset into 27 site-types. Each of these site-types represent a structural motif comprised of two or more residue conservations, derived using another in-house tool for superposing binding sites, PocketAlign. The 27 site-types could be grouped further into 9 super-types by considering partial similarities in the sites, which indicated that the individual site-types comprise different combinations of one or more site features. A scan across PDB using the 27 structural motifs determined the motifs to be specific to NTP binding sites, and a computational alanine mutagenesis indicated that residues identified to be highly conserved in the motifs are also most contributing to binding. Alternate orientations of the ligand in several site-types were observed and rationalized, indicating the possibility of some residues serving as anchors for NTP recognition. The presence of multiple site-types and the grouping of multiple folds into each site-type is strongly suggestive of convergent evolution. Knowledge of determinants obtained from this study will be useful for detecting function in unknown proteins. Proteins 2017; 85:1699-1712. © 2017 Wiley Periodicals, Inc.

Systematic investigation of sequence and structural motifs that recognize ATP

Interaction between ATP, a multifunctional and ubiquitous nucleotide, and proteins initializes phosphorylation, polypeptide synthesis and ATP hydrolysis which supplies energy for metabolism. However, current knowledge concerning the mechanisms through which ATP is recognized by proteins is incomplete, scattered, and inaccurate. We systemically investigate sequence and structural motifs of proteins that recognize ATP. We identified three novel motifs and refined the known p-loop and class II aminoacyl-tRNA synthetase motifs. The five motifs define five distinct ATP-protein interaction modes which concern over 5% of known protein structures. We demonstrate that although these motifs share a common GXG tripeptide they recognize ATP through different functional groups. The p-loop motif recognizes ATP through phosphates, class II aminoacyl-tRNA synthetase motif targets adenosine and the other three motifs recognize both phosphates and adenosine. We show that some motifs are shared by different enzyme types. Statistical tests demonstrate that the five sequence motifs are significantly associated with the nucleotide binding proteins. Large-scale test on PDB reveals that about 98% of proteins that include one of the structural motifs are confirmed to bind ATP.

Visualization of Protein 3D Structures in ‘Double-Centroid’ Reduced Representation

This article is of two parts: (a) the development of a protein reduced representation and its implementation in a Web server; and (b) the use of the reduced protein representation in the modeling of the binding site of a given ligand and the screening for the model in other protein 3D structures. Current methods of reduced protein 3D structure representation such as the Ca trace method not only lack essential molecular detail, but also ignore the chemical properties of the component amino acid side chains. This chapter describes a reduced protein 3D structure representation called “double-centroid reduced representation” and presents a visualization tool called the “DCRR Web Server” that graphically displays a protein 3D structure in DCRR along with non-covalent intra- and intermolecular hydrogen bonding and van der Waals interactions. In the DCRR model, each amino acid residue is represented as two points: the centroid of the backbone atoms and that of the side chain atoms; in the visualization Web server, they and the non-bonded interactions are color-coded for easy identification. The visualization tool in this chapter is implemented in MATLAB and is the first for a reduced protein representation as well as one that simultaneously displays non-covalent interactions in the molecule. The DCRR model reduces the atomicity of the protein structure by ~75% while capturing the essential chemical properties of the component amino acids. The second half of this chapter describes the application of this reduced representation to the modeling and screening of ligand binding sites using a data model termed the “tetrahedral motif.” This type of ligand binding site modeling and screening presents a novel type of pharmacophore modeling and screening, one that depends on a reduced protein representation.

GIRAF: a method for fast search and flexible alignment of ligand binding interfaces in proteins at atomic resolution

Comparison and classification of protein structures are fundamental means to understand protein functions. Due to the computational difficulty and the ever-increasing amount of structural data, however, it is in general not feasible to perform exhaustive all-against-all structure comparisons necessary for comprehensive classifications. To efficiently handle such situations, we have previously proposed a method, now called GIRAF. We herein describe further improvements in the GIRAF protein structure search and alignment method. The GIRAF method achieves extremely efficient search of similar structures of ligand binding sites of proteins by exploiting database indexing of structural features of local coordinate frames. In addition, it produces refined atom-wise alignments by iterative applications of the Hungarian method to the bipartite graph defined for a pair of superimposed structures. By combining the refined alignments based on different local coordinate frames, it is made possible to align structures involving domain movements. We provide detailed accounts for the database design, the search and alignment algorithms as well as some benchmark results.

LigAlign: flexible ligand-based active site alignment and analysis

Ligand-based active site alignment is a widely adopted technique for the structural analysis of protein-ligand complexes. However, existing tools for ligand alignment treat the ligands as rigid objects even though most biological ligands are flexible. We present LigAlign, an automated system for flexible ligand alignment and analysis. When performing rigid alignments, LigAlign produces results consistent with manually annotated structural motifs. In performing flexible alignments, LigAlign automatically produces biochemically reasonable ligand fragmentations and subsequently identifies conserved structural motifs that are not detected by rigid alignment.

Ligand-binding site prediction of proteins based on known fragment-fragment interactions

Motivation: The identification of putative ligand-binding sites on proteins is important for the prediction of protein function. Knowledge-based approaches using structure databases have become interesting, because of the recent increase in structural information. Approaches using binding motif information are particularly effective. However, they can only be applied to well-known ligands that frequently appear in the structure databases.

Results: We have developed a new method for predicting the binding sites of chemically diverse ligands, by using information about the interactions between fragments. The selection of the fragment size is important. If the fragments are too small, then the patterns derived from the binding motifs cannot be used, since they are many-body interactions, while using larger fragments limits the application to well-known ligands. In our method, we used the main and side chains for proteins, and three successive atoms for ligands, as fragments. After superposition of the fragments, our method builds the conformations of ligands and predicts the binding sites. As a result, our method could accurately predict the binding sites of chemically diverse ligands, even though the Protein Data Bank currently contains a large number of nucleotides. Moreover, a further evaluation for the unbound forms of proteins revealed that our building up procedure was robust to conformational changes induced by ligand binding.

Availability: Our method, named ‘BUMBLE’, is available at http://bumble.hgc.jp/

Contact:kasahara@cb.k.u-tokyo.ac.jp

Supplementary information:Supplementary Material is available at Bioinformatics online.

Comprehensive structural classification of ligand-binding motifs in proteins

Comprehensive knowledge of protein-ligand interactions should provide a useful basis for annotating protein functions, studying protein evolution, engineering enzymatic activity, and designing drugs. To investigate the diversity and universality of ligand-binding sites in protein structures, we conducted the all-against-all atomic-level structural comparison of over 180,000 ligand-binding sites found in all the known structures in the Protein Data Bank by using a recently developed database search and alignment algorithm. By applying a hybrid top-down-bottom-up clustering analysis to the comparison results, we determined approximately 3000 well-defined structural motifs of ligand-binding sites. Apart from a handful of exceptions, most structural motifs were found to be confined within single families or superfamilies, and to be associated with particular ligands. Furthermore, we analyzed the components of the similarity network and enumerated more than 4000 pairs of structural motifs that were shared across different protein folds.

Cytoplasmic steps of peptidoglycan biosynthesis

The biosynthesis of bacterial cell wall peptidoglycan is a complex process that involves enzyme reactions that take place in the cytoplasm (synthesis of the nucleotide precursors) and on the inner side (synthesis of lipid-linked intermediates) and outer side (polymerization reactions) of the cytoplasmic membrane. This review deals with the cytoplasmic steps of peptidoglycan biosynthesis, which can be divided into four sets of reactions that lead to the syntheses of (1) UDP-N-acetylglucosamine from fructose 6-phosphate, (2) UDP-N-acetylmuramic acid from UDP-N-acetylglucosamine, (3) UDP-N-acetylmuramyl-pentapeptide from UDP-N-acetylmuramic acid and (4) D-glutamic acid and dipeptide D-alanyl-D-alanine. Recent data concerning the different enzymes involved are presented. Moreover, special attention is given to (1) the chemical and enzymatic synthesis of the nucleotide precursor substrates that are not commercially available and (2) the search for specific inhibitors that could act as antibacterial compounds.

Cognate ligand domain mapping for enzymes

Here, we present an automatic assignment of potential cognate ligands to domains of enzymes in the CATH and SCOP protein domain classifications on the basis of structural data available in the wwPDB. This procedure involves two steps; firstly, we assign the binding of particular ligands to particular domains; secondly, we compare the chemical similarity of the PDB ligands to ligands in KEGG in order to assign cognate ligands. We find that use of the Enzyme Commission (EC) numbers is necessary to enable efficient and accurate cognate ligand assignment. The PROCOGNATE database currently has cognate ligand mapping for 3277 (4118) protein structures and 351 (302) superfamilies, as described by the CATH and (SCOP) databases, respectively. We find that just under half of all ligands are only and always bound by a single domain, with 16% bound by more than one domain and the remainder of the ligands showing a variety of binding modes. This finding has implications for domain recombination and the evolution of new protein functions. Domain architecture or context is also found to affect substrate specificity of particular domains, and we discuss example cases. The most popular PDB ligands are all found to be generic components of crystallisation buffers, highlighting the non-cognate ligand problem inherent in the PDB. In contrast, the most popular cognate ligands are all found to be universal cellular currencies of reducing power and energy such as NADH, FADH2 and ATP, respectively, reflecting the fact that the vast majority of enzymatic reactions utilise one of these popular co-factors. These ligands all share a common adenine ribonucleotide moiety, suggesting that many different domain superfamilies have converged to bind this chemical framework.

An empirical approach for detecting nucleotide-binding sites on proteins

Protein structure data in the PDB (Protein Data Bank) were used to construct empirical scores of nucleotide-protein interactions. A simple strategy to evaluate the spatial distribution of protein atoms around the base moieties of nucleotides was applied to categorize adenine, guanine, nicotinamide and flavin nucleotide-binding sites. In addition to the known nucleotide-binding motifs, the empirical scores detected several other features that were shared among proteins with different folds. The empirical scores were also used to predict the binding sites on protein molecules and a comprehensive test of the prediction system was performed. As a result, adenine, guanine, nicotinamide and flavin sites were detected with efficiencies of 31, 29, 32 and 40%, respectively. The predictions were judged to be successful if the predicted base with the best score was located within a 3.0 A r.m.s.d. from the known ligand positions.

Functional annotation by identification of local surface similarities: a novel tool for structural genomics

Background

Protein function is often dependent on subsets of solvent-exposed residues that may exist in a similar three-dimensional configuration in non homologous proteins thus having different order and/or spacing in the sequence. Hence, functional annotation by means of sequence or fold similarity is not adequate for such cases.

Results

We describe a method for the function-related annotation of protein structures by means of the detection of local structural similarity with a library of annotated functional sites. An automatic procedure was used to annotate the function of local surface regions. Next, we employed a sequence-independent algorithm to compare exhaustively these functional patches with a larger collection of protein surface cavities. After tuning and validating the algorithm on a dataset of well annotated structures, we applied it to a list of protein structures that are classified as being of unknown function in the Protein Data Bank. By this strategy, we were able to provide functional clues to proteins that do not show any significant sequence or global structural similarity with proteins in the current databases.

Conclusion

This method is able to spot structural similarities associated to function-related similarities, independently on sequence or fold resemblance, therefore is a valuable tool for the functional analysis of uncharacterized proteins. Results are available at

FSSA: a novel method for identifying functional signatures from structural alignments

MOTIVATION

It is commonly believed that sequence determines structure, which in turn determines function. However, the presence of many proteins with the same structural fold but different functions suggests that global structure and function do not always correlate well.

RESULTS

We propose a method for accurate functional annotation, based on identification of functional signatures from structural alignments (FSSA) using the Structural Classification of Proteins (SCOP) database. The FSSA method is superior at function discrimination and classification compared with several methods that directly inherit functional annotation information from homology inference, such as Smith-Waterman, PSI-BLAST, hidden Markov models and structure comparison methods, for a large number of structural fold families. Our results indicate that the contributions of amino acid residue types and positions to structure and function are largely separable for proteins in multi-functional fold families.

Novel CalphaNN structural motif for protein recognition of phosphate ions

Phosphate is one of the most frequently exploited chemical moieties in nature, present in a wide range of naturally occurring and critically important small molecules. Several phosphate group recognition motifs have been found for a few narrow groups of proteins, but for many protein families and folds the mode of phosphate recognition remains unclear. Here, we have analyzed the structures of all fold-representative protein-ligand complexes listed in the FSSP database, regardless of whether the bound ligand included a phosphate group. Based on a phosphate-binding motif that we identified in pyridoxal phosphate binding proteins, we have identified a new anion-binding structural motif, CalphaNN, common to 104 fold-representative protein structures that belong to 62 different folds, of which 86% of the fold-representative structures (51 folds) bind phosphate or lone sulfate ions. This motif leads to a precise mode for phosphate group recognition forming a structure where atoms of the phosphate group occupy the most favorable stabilizing positions. The anion-binding CalphaNN motif is based only on main-chain atoms from three adjacent residues, has a conservative betaalphaalpha or betaalphabeta geometry, and recognizes the free phosphate (sulfate) ion as well as one or more phosphate groups in nucleotides and in a variety of cofactors. Moreover, the CalphaNN motif is positioned in functionally important regions of protein structures and often residues of the motif directly participate in the function of the protein.

Activation of phosphorothionate pesticides based on a cytochrome P450 BM-3 (CYP102 A1) mutant for expanded neurotoxin detection in food using acetylcholinesterase biosensors

A novel enzymatic in vitro activation method for phosphorothionates has been developed to allow their detection with acetylcholinesterase (AChE) biosensors. Activation is necessary because this group of insecticides shows nearly no inhibitory effect toward AChE in their pure nonmetabolized form. In contrast, they exert a strong inhibitory effect on AChE after oxidation as it takes place by metabolic activation in higher organisms. Standard chemical methods to oxidize phosphorothionates showed inherent disadvantages that impede their direct use in food analysis. In contrast, a genetically engineered triple mutant of P450 BM-3 (CYP102 A1) could convert the two frequently used insecticides parathion and chlorpyrifos into their oxo variants as was confirmed by GC/MS measurements. The wild-type protein was unable to do so. In the case of chlorpyrifos, the enzymatic activation was as good as the chemical oxidation. In the case of parathion, the P450 activation was more efficient than the oxidation by NBS but neither activation method yielded an AChE inhibition that was as high as with paraoxon. The application of the method to infant food in combination with a disposable AChE biosensor enabled detection of chlorpyrifos and parathion at concentrations down to 20 microg/kg within an overall assay time of 95 min.

Electrostatic Potential of Nucleotide-free Protein is Sufficient for Discrimination Between Adenine and Guanine-specific Binding Sites

Despite sharing many common features, adenine-binding and guanine-binding sites in proteins often show a clear preference for the cognate over the non-cognate ligand. We have analyzed electrostatic potential (ESP) patterns at adenine and guanine-binding sites of a large number of non-redundant proteins where each binding site was first annotated as adenine/guanine-specific or non-specific from a survey of primary literature. We show that more than 90% of ESP variance at the binding sites is accounted for by only two principal component ESP vectors, each aligned to molecular dipoles of adenine and guanine. Projected on these principal component vectors, the adenine/guanine-specific and non-specific binding sites, including adenine-containing dinucleotides, show non-overlapping distributions. Adenine or guanine specificities of the binding sites also show high correlation with the corresponding electrostatic replacement (cognate by non-cognate ligand) energies. High correlation coefficients (0.94 for 35 adenine-binding sites and 1.0 for 20 guanine-binding sites) were obtained when adenine/guanine specificities were predicted using the replacement energies. Our results demonstrate that ligand-free protein ESP is an excellent indicator for discrimination between adenine and guanine-specific binding sites and that ESP of ligand-free protein can be used as a tool to annotate known and putative purine-binding sites in proteins as adenine or guanine-specific.

"Acceptor-donor-acceptor" motifs recognize the Watson-Crick, Hoogsteen and Sugar "donor-acceptor-donor" edges of adenine and adenosine-containing ligands

Nucleotides are among the most extensively exploited chemical moieties in nature and, as a part of a handful of different protein ligands, nucleotides play key roles in energy transduction, enzymatic catalysis and regulation of protein function. We have previously reported that in many proteins with different folds and functions a distinctive adenine-binding motif is involved in the recognition of the Watson-Crick edge of adenine. Here, we show that many proteins do have clear structural motifs that recognize adenosine (and some other nucleotides and nucleotide analogs) not only through the Watson-Crick edge, but also through the sugar and Hoogsteen edges. Each of the three edges of adenosine has a donor-acceptor-donor (DAD) pattern that is often recognized by proteins via a complementary acceptor-donor-acceptor (ADA) motif, whereby three distinct hydrogen bonds are formed: two conventional N-H...O and N-H...N hydrogen bonds, and one weak C-H...O hydrogen bond. The local conformation of the adenine-binding loop is betabetabeta or betabetaalpha and reflects the mode of nucleotide binding. Additionally, we report 21 proteins from five different folds that simultaneously recognize both the sugar edge and the Watson-Crick edge of adenine. In these proteins a unique beta-loop-beta supersecondary structure grasps an adenine-containing ligand between two identical adenine-binding motifs as part of the betaalphabeta-loop-beta fold.

Nucleotide binding to nucleoside diphosphate kinases: X-ray structure of human NDPK-A in complex with ADP and comparison to protein kinases

NDPK-A, product of the nm23-H1 gene, is one of the two major isoforms of human nucleoside diphosphate kinase. We analyzed the binding of its nucleotide substrates by fluorometric methods. The binding of nucleoside triphosphate (NTP) substrates was detected by following changes of the intrinsic fluorescence of the H118G/F60W variant, a mutant protein engineered for that purpose. Nucleoside diphosphate (NDP) substrate binding was measured by competition with a fluorescent derivative of ADP, following the fluorescence anisotropy of the derivative. We also determined an X-ray structure at 2.0A resolution of the variant NDPK-A in complex with ADP, Ca(2+) and inorganic phosphate, products of ATP hydrolysis. We compared the conformation of the bound nucleotide seen in this complex and the interactions it makes with the protein, with those of the nucleotide substrates, substrate analogues or inhibitors present in other NDP kinase structures. We also compared NDP kinase-bound nucleotides to ATP bound to protein kinases, and showed that the nucleoside monophosphate moieties have nearly identical conformations in spite of the very different protein environments. However, the beta and gamma-phosphate groups are differently positioned and oriented in the two types of kinases, and they bind metal ions with opposite chiralities. Thus, it should be possible to design nucleotide analogues that are good substrates of one type of kinase, and poor substrates or inhibitors of the other kind.

A consensus-binding structure for adenine at the atomic level permits searching for the ligand site in a wide spectrum of adenine-containing complexes

Attempts to derive structural features of ligand-binding sites have traditionally involved seeking commonalities at the residue level. Recently, structural studies have turned to atomic interactions of small molecular fragments to extract common binding-site properties. Here, we explore the use of larger ligand elements to derive a consensus binding structure for the ligand as a whole. We superimposed multiple molecular structures from a nonredundant set of adenosine-5'-triphosphate (ATP) protein complexes, using the adenine moiety as template. Clustered binding-site atoms of compatible atomic classes forming attractive contacts with the adenine probe were extracted. A set of atomic clusters characterizing the adenine binding pocket was then derived. Among the clusters are three vertices representing the interactions of adenine atom N6 with its protein-binding niche. These vertices, together with atom C6 of the purine ring system, complete the set of four vertices for the pyramid-like structure of the N6 anchor atom. Also, the sequence relationship for the adenine-binding loop interacting with the C2-N6 end of the conjugated ring system is expanded to include a third hydrophilic cluster interacting with atom N1. A search procedure involving interatomic distances between cluster centers was formulated and applied to seek putative binding sites in test cases. The results show that a consensus network of clusters, based on an adenine probe and an ATP-complexed training set of proteins, is sufficient to recognize the experimental cavity for adenine in a wide spectrum of ligand-protein complexes.

Conserved positions for ribose recognition: importance of water bridging interactions among ATP, ADP and FAD-protein complexes

Analysis of the spatial arrangement of protein and water atoms that form polar interactions with ribose has been performed for a structurally non-redundant dataset of ATP, ADP and FAD-protein complexes. The 26 ligand-protein structures were separated into two groups corresponding to the most populated furanose ring conformations (N and S-domains). Four conserved positions were found for S-domain protein-ligand complexes and five for N-domain complexes. Multiple protein folds and secondary structural elements were represented at a single conserved position. The following novel points were revealed: (i) Two complementary positions sometimes combine to describe a putative atomic spatial location for a specific conserved binding spot. (ii) More than one third of the interactions scored were water-mediated. Thus, conserved spatial positions rich in water atoms are a significant feature of ribose-protein complexes.

Structure of triphosphoryl nucleotide bound at the active site of yeast hexokinase: 1H-nuclear magnetic resonance study

Conformation of a nonhydrolyzable adenosine triphosphate (ATP) analogue, adenylyl-(beta,gamma-methylene)-diphosphonate (AMPPCP) bound at the active site of yeast hexokinase-PII was determined by proton two-dimensional transferred nuclear Overhauser effect spectroscopy (TRNOESY) and molecular dynamics simulations. The effect of the glucose-induced domain closure on the conformation of the nucleotide was evaluated by making measurements on two different complexes: PII AMPPCPMg(II) and PII-Glc.AMPPCPMg(Il). TRNOE measurements were made at 500 MHz, 10 degress C, as a function of several mixing times varying in the range of 40 to 200 ms. Interproton distances derived from the analysis of NOE buildup curves were used as restraints in molecular dynamics simulations to determine the conformation of the enzyme bound nucleotide. The adenosine moiety was found to bind in high anti conformation with a glycosidic torsion angle chi = 48 +/- 5 degrees in both complexes. However, significant differences in the conformations of the ribose and triphosphoryl chain of the nucleotide are observed between the two complexes. The phase angles of pseudorotation P in PII.AMPPCPMg(II) and PII.Glc.AMPPCPMg(II) are 87 degrees and 77 degrees, describing a OE and OT4 sugar pucker and the amplitudes of the sugar pucker (tau) are 37 degrees and 61 degrees, respectively.

Adenine recognition: a motif present in ATP-, CoA-, NAD-, NADP-, and FAD-dependent proteins

Adenosine triphosphate (ATP) plays an essential role in energy transfer within the cell. In the form of NAD, adenine participates in multiple redox reactions. Phosphorylation and ATP-hydrolysis reactions have key roles in signal transduction and regulation of many proteins, especially enzymes. In each cell, proteins with many different functions use adenine and its derivatives as ligands; adenine, of course, is present in DNA and RNA. We show that an adenine binding motif, which differs according to the backbone chain direction of a loop that binds adenine (and in one variant by the participation of an aspartate side-chain), is common to many proteins; it was found from an analysis of all adenylate-containing protein structures from the Protein Data Bank. Indeed, 224 protein-ligand complexes (86 different proteins) from a total of 645 protein structure files bind ATP, CoA, NAD, NADP, FAD, or other adenine-containing ligands, and use the same structural elements to recognize adenine, regardless of whether the ligand is a coenzyme, cofactor, substrate, or an allosteric effector. The common adenine-binding motif shown in this study is simple to construct. It uses only (1) backbone polar interactions that are not dependent on the protein sequence or particular properties of amino acid side-chains, and (2) nonspecific hydrophobic interactions. This is probably why so many different proteins with different functions use this motif to bind an adenylate-containing ligand. The adenylate-binding motif reported is present in "ancient proteins" common to all living organisms, suggesting that adenine-containing ligands and the common motif for binding them were exploited very early in evolution. The geometry of adenine binding by this motif mimics almost exactly the geometry of adenine base-pairing seen in DNA and RNA.

Structural homology between DNA binding sites of DNA polymerase beta and DNA topoisomerase II

Unsaturated long-chain fatty acids selectively bind to the DNA binding sites of DNA polymerase beta and DNA topoisomerase II, and inhibit their activities, although the amino acid sequences of these enzymes are markedly different from each other. Computer modeling analysis revealed that the fatty acid interaction interface in both enzymes has a group of four amino acid residues in common, forming a pocket which binds to the fatty acid molecule. The four amino acid residues were Thr596, His735, Leu741 and Lys983 for yeast DNA topoisomerase II, corresponding to Thr79, His51, Leu11 and Lys35 for rat DNA polymerase beta. Using three-dimensional structure model analysis, we determined the spatial positioning of specific amino acid residues binding to the fatty acids in DNA topoisomerase II, and subsequently obtained supplementary information to build the structural model.

The molecular basis of vancomycin resistance in clinically relevant Enterococci: crystal structure of D-alanyl-D-lactate ligase (VanA)

d-alanine-d-lactate ligase from Enterococcus faecium BM4147 is directly responsible for the biosynthesis of alternate cell-wall precursors in bacteria, which are resistant to the glycopeptide antibiotic vancomycin. The crystal structure has been determined with data extending to 2.5-A resolution. This structure shows that the active site has unexpected interactions and is distinct from previous models for d-alanyl-d-lactate ligase mechanistic studies. It appears that the preference of the enzyme for lactate as a ligand over d-alanine could be mediated by electrostatic effects and/or a hydrogen-bonding network, which principally involve His-244. The structure of d-alanyl-d-lactate ligase provides a revised interpretation of the molecular events that lead to vancomycin resistance.

Phosphatidylinositol phosphate kinase: a link between protein kinase and glutathione synthase folds

Comparisons of serine/threonine protein kinase (PK) and type IIbeta phosphatidylinositol phosphate kinase (PIPK) structures with each other and also with other proteins reveal structural and functional similarity between the two kinases and proteins of the glutathione synthase fold (ATP-grasp). This suggests that these enzymes are evolutionarily related. The structure of PIPK, which clearly resembles both PK and ATP-grasp, provides a link between the two proteins and establishes that the C-terminal domains of PK, PIPK and ATP-grasp share the same fold. The functional implications of the proposed homology are discussed.

Binding of nucleotide triphosphates to cardiotoxin analogue II from the Taiwan cobra venom (Naja naja atra). Elucidation of the structural interactions in the dATP-cardiotoxin analogue ii complex

Snake venom cardiotoxins have been recently shown to block the enzymatic activity of phospholipid protein kinase and Na+,K+-ATPase. To understand the molecular basis for the inhibitory effects of cardiotoxin on the action of these enzymes, the nucleotide triphosphate binding ability of cardiotoxin analogue II (CTX II) from the Taiwan cobra (Naja naja atra) venom is investigated using a variety of spectroscopic techniques such as fluorescence, circular dichroism, and two-dimensional NMR. CTX II is found to bind to all the four nucleotide triphosphates (ATP, UTP, GTP, and CTP) with similar affinity. Detailed studies of the binding of dATP to CTX II indicated that the toxin molecule is significantly stabilized in the presence of the nucleotide. Molecular modeling, based on the NOEs observed for the dATP.CTX II complex, reveals that dATP binds to the CTX II molecule at the groove enclosed between the N- and C-terminal ends of the toxin molecule. Based on the results obtained in the present study, a molecular mechanism to account for the inhibition of the enzymatic activity of the phospholipid-sensitive protein kinase and Na+,K+-ATPase is also proposed.

From protein structure to function

Several databases of protein structural families now exist-organised according to both evolutionary relationships and common folding arrangements. Although these lag behind sequence databases in size, the prospect of structural genomics initiatives means that they may soon include representatives of many of the sequence families. To some extent, functional information can be derived from structural similarity. For some structural families, their function is highly conserved, whereas, for others, it can only be inherited or derived on the basis of additional information (e.g. sequence patterns, common residue clusters and characteristic surface properties).

Coding sequence mutations in the alpha subunit of propionyl-CoA carboxylase in patients with propionic acidemia

Propionic acidemia is a rare autosomal recessive disorder of intermediary metabolism. It is caused by a deficiency of the mitochondrial enzyme propionyl-CoA carboxylase (PCC, EC 6.4.1.3), a heteropolymeric protein composed of two subunits, alpha and beta. PCC requires ATP and biotin as cofactors for the reaction, the latter enzymatically added onto the alpha subunit. We investigated coding sequence mutations in the alpha subunit of PCC by analyzing fibroblast RNA from propionic acidemia patients deficient in alpha subunit function by single-strand conformation polymorphism and direct sequencing. Five missense mutations and one short in-frame deletion were found among different patients. Four mutations were located in the putative biotin carboxylase domain, whereas the two others were within the 67-amino-acid C-terminal domain previously shown to be required to obtain biotinylation of the alpha subunit. We analyzed fibroblast extracts for the presence of a biotinylated alpha subunit by Western blot analysis using streptavidin coupled to alkaline phosphatase. Four of five cell lines failed to show a biotinylated alpha subunit, regardless of the position of the mutations within the coding sequence. Two mutations located in the biotinylation domain were expressed in an Escherichia coli-based system and shown to abolish biotinylation of the domain. The results suggest that most mutations have a severe impact on the stability or the functionality of the alpha subunit.

Finding local structural similarities among families of unrelated protein structures: a generic non-linear alignment algorithm

We have developed a generic tool for the automatic identification of regions of local structural similarity in unrelated proteins having different folds, as well as for defining more global similarities that result from homologous protein structures. The computer program GENFIT has evolved from the genetic algorithm-based three-dimensional protein structure comparison program GA_FIT. GENFIT, however, can locate and superimpose regions of local structural homology regardless of their position in a pair of structures, the fold topology, or the chain direction. Furthermore, it is possible to restrict the search to a volume centered about a region of interest (e.g., catalytic site, ligand-binding site) in two protein structures. We present a number of examples to illustrate the function of the program, which is a parallel processing implementation designed for distribution to multiple machines over a local network or to run on a single multiprocessor computer.

Structural interpretation of site-directed mutagenesis and specificity of the catalytic subunit of protein kinase CK2 using comparative modelling

The catalytic subunit of protein kinase casein kinase 2 (CK2alpha), which has specificity for both ATP and GTP, shows significant amino acid sequence similarity to the cyclin-dependent kinase 2 (CDK2). We constructed site-directed mutants of CK2alpha and used a three-dimensional model to investigate the basis for the dual specificity. Introduction of Phe and Gly at positions 50 and 51, in order to restore the pattern of the glycine-rich motif, did not seriously affect the specificity for ATP or GTP. We show that the dual specificity probably originates from the loop situated around the position His115 to Asp120 (HVNNTD). The insertion of a residue in this loop in CK2 alpha subunits, compared with CDK2 and other kinases, might orient the backbone to interact with the base A and G; this insertion is conserved in all known CK2alpha. The mutant deltaN118, the design of which was based on the modelling, showed reduced affinity for GTP as predicted from the model. Other mutants were intended to probe the integrity of the catalytic loop, alter the polarity of a buried residue and explore the importance of the carboxy terminus. Introduction of Arg to replace Asn189, which is mapped on the activation loop, results in a mutant with decreased k(cat), possibly as a result of disruption of the interaction between this residue and basic residues in the vicinity. Truncation at position 331 eliminates the last 60 residues of the alpha subunit and this mutant has a reduced catalytic efficiency compared with the wild-type. Catalytic efficiency is restored in the truncation mutant by the replacement of a potentially buried Glu at position 252 by Lys, probably owing to a higher stability resulting from the formation of a salt bridge between Lys252 and Asp208.

Structural motif of phosphate-binding site common to various protein superfamilies: all-against-all structural comparison of protein-mononucleotide complexes

In order to search for a common structural motif in the phosphate-binding sites of protein-mononucleotide complexes, we investigated the structural variety of phosphate-binding schemes by an all-against-all comparison of 491 binding sites found in the Protein Data Bank. We found four frequently occurring structural motifs composed of protein atoms interacting with phosphate groups, each of which appears in different protein superfamilies with different folds. The most frequently occurring motif, which we call the structural P-loop, is shared by 13 superfamilies and is characterized by a four-residue fragment, GXXX, interacting with a phosphate group through the backbone atoms. Various sequence motifs, including Walker's A motif or the P-loop, turn out to be a structural P-loop found in a few specific superfamilies. The other three motifs are found in pairs of superfamilies: protein kinase and glutathione synthetase ATPase domain like, actin-like ATPase domain and nucleotidyltransferase, and FMN-linked oxidoreductase and PRTase.

Enzyme-mononucleotide interactions: three different folds share common structural elements for ATP recognition

Three ATP-dependent enzymes with different folds, cAMP-dependent protein kinase, D-Ala:D-Ala ligase and the alpha-subunit of the alpha2beta2 ribonucleotide reductase, have a similar organization of their ATP-binding sites. The most meaningful similarity was found over 23 structurally equivalent residues in each protein and includes three strands each from their beta-sheets, in addition to a connecting loop. The equivalent secondary structure elements in each of these enzymes donate four amino acids forming key hydrogen bonds responsible for the common orientation of the "AMP" moieties of their ATP-ligands. One lysine residue conserved throughout the three families binds the alpha-phosphate in each protein. The common fragments of structure also position some, but not all, of the equivalent residues involved in hydrophobic contacts with the adenine ring. These examples of convergent evolution reinforce the view that different proteins can fold in different ways to produce similar structures locally, and nature can take advantage of these features when structure and function demand it, as shown here for the common mode of ATP-binding by three unrelated proteins.

Two "unrelated" families of ATP-dependent enzymes share extensive structural similarities about their cofactor binding sites

Two proteins, D-alanine:D-alanine ligase and cAMP-dependent protein kinase, share a remarkable degree of structural convergence despite having different three-dimensional folds and different enzymatic functions. Here we report that as many as 103 residues from 10 segments form two identical super-secondary structures between which the cofactor ATP is bound. The cofactor, two bound metal cations, and several water molecules form a large network of electrostatic and hydrophobic interactions common to both enzymes, and these are mediated by the similar placement of equivalent amino acids within the common supersecondary structures.

How do kinases transfer phosphoryl groups?

Understanding how phosphoryl transfer is accomplished by kinases, a ubiquitous group of enzymes, is central to many biochemical processes. Qualitative analysis of the crystal structures of enzyme-substrate complexes of kinases reveals structural features of these enzymes important to phosphoryl transfer. Recently determined crystal structures which mimic the transition state complex have added new insight into the debate as to whether kinases use associative or dissociative mechanisms of catalysis.

The structure of SAICAR synthase: an enzyme in the de novo pathway of purine nucleotide biosynthesis

BACKGROUND

The biosynthesis of key metabolic components is of major interest to biologists. Studies of de novo purine synthesis are aimed at obtaining a deeper understanding of this central pathway and the development of effective chemotherapeutic agents. Phosphoribosylaminoimidazolesuccinocarboxamide (SAICAR) synthase catalyses the seventh step out of ten in the biosynthesis of purine nucleotides. To date, only one structure of an enzyme involved in purine biosynthesis has been reported: adenylosuccinate synthetase, which catalyses the first committed step in the synthesis of AMP from IMP.

RESULTS

We report the first three-dimensional structure of a SAICAR synthase, from Saccharomyces cerevisiae. It is a monomer with three domains. The first two domains consist of antiparallel beta sheets and the third is composed of two alpha helices. There is a long deep cleft made up of residues from all three domains. Comparison of SAICAR synthases by alignment of their sequences reveals a number of conserved residues, mostly located in the cleft. The presence of two sulphate ions bound in the cleft, the structure of SAICAR synthase in complex with ATP and a comparison of this structure with that of other ATP-dependent proteins point to the interdomain cleft as the location of the active site.

CONCLUSIONS

The topology of the first domain of SAICAR synthase resembles that of the N-terminal domain of proteins belonging to the cyclic AMP-dependent protein kinase family. The fold of the second domain is similar to that of members of the D-alanine:D-alanine ligase family. Together these enzymes form a new superfamily of mononucleotide-binding domains. There appears to be no other enzyme, however, which is composed of the same combination of three domains, with the individual topologies found in SAICAR synthase.