Resolving protein structure‐function‐binding site relationships from a binding site similarity network perspective

Functional annotation is seldom straightforward with complexities arising due to functional divergence in protein families or functional convergence between non‐homologous protein families, leading to mis‐annotations. An enzyme may contain multiple domains and not all domains may be involved in a given function, adding to the complexity in function annotation. To address this, we use binding site information from bound cognate ligands and catalytic residues, since it can help in resolving fold‐function relationships at a finer level and with higher confidence. A comprehensive database of 2,020 fold‐function‐binding site relationships has been systematically generated. A network‐based approach is employed to capture the complexity in these relationships, from which different types of associations are deciphered, that identify versatile protein folds performing diverse functions, same function associated with multiple folds and one‐to‐one relationships. Binding site similarity networks integrated with fold, function, and ligand similarity information are generated to understand the depth of these relationships. Apart from the observed continuity in the functional site space, network properties of these revealed versatile families with topologically different or dissimilar binding sites and structural families that perform very similar functions. As a case study, subtle changes in the active site of a set of evolutionarily related superfamilies are studied using these networks. Tracing of such similarities in evolutionarily related proteins provide clues into the transition and evolution of protein functions. Insights from this study will be helpful in accurate and reliable functional annotations of uncharacterized proteins, poly‐pharmacology, and designing enzymes with new functional capabilities. Proteins 2017; 85:1319–1335. © 2017 Wiley Periodicals, Inc.     

Assigning genomic sequences to CATH

We report the latest release (version 1.6) of the CATH protein domains database (http://www.biochem.ucl. ac.uk/bsm/cath ). This is a hierarchical classification of 18 577 domains into evolutionary families and structural groupings. We have identified 1028 homo-logous superfamilies in which the proteins have both structural, and sequence or functional similarity. These can be further clustered into 672 fold groups and 35 distinct architectures. Recent developments of the database include the generation of 3D templates for recognising structural relatives in each fold group, which has led to significant improvements in the speed and accuracy of updating the database and also means that less manual validation is required. We also report the establishment of the CATH-PFDB (Protein Family Database), which associates 1D sequences with the 3D homologous superfamilies. Sequences showing identifiable homology to entries in CATH have been extracted from GenBank using PSI-BLAST. A CATH-PSIBLAST server has been established, which allows you to scan a new sequence against the database. The CATH Dictionary of Homologous Superfamilies (DHS), which contains validated multiple structural alignments annotated with consensus functional information for evolutionary protein superfamilies, has been updated to include annotations associated with sequence relatives identified in GenBank. The DHS is a powerful tool for considering the variation of functional properties within a given CATH superfamily and in deciding what functional properties may be reliably inherited by a newly identified relative.     

Evolution of function in protein superfamilies, from a structural perspective

The recent growth in protein databases has revealed the functional diversity of many protein superfamilies. We have assessed the functional variation of homologous enzyme superfamilies containing two or more enzymes, as defined by the CATH protein structure classification, by way of the Enzyme Commission (EC) scheme. Combining sequence and structure information to identify relatives, the majority of superfamilies display variation in enzyme function, with 25 % of superfamilies in the PDB having members of different enzyme types. We determined the extent of functional similarity at different levels of sequence identity for 486,000 homologous pairs (enzyme/enzyme and enzyme/non-enzyme), with structural and sequence relatives included. For single and multi-domain proteins, variation in EC number is rare above 40 % sequence identity, and above 30 %, the first three digits may be predicted with an accuracy of at least 90 %. For more distantly related proteins sharing less than 30 % sequence identity, functional variation is significant, and below this threshold, structural data are essential for understanding the molecular basis of observed functional differences. To explore the mechanisms for generating functional diversity during evolution, we have studied in detail 31 diverse structural enzyme superfamilies for which structural data are available. A large number of variations and peculiarities are observed, at the atomic level through to gross structural rearrangements. Almost all superfamilies exhibit functional diversity generated by local sequence variation and domain shuffling. Commonly, substrate specificity is diverse across a superfamily, whilst the reaction chemistry is maintained. In many superfamilies, the position of catalytic residues may vary despite playing equivalent functional roles in related proteins. The implications of functional diversity within supefamilies for the structural genomics projects are discussed. More detailed information on these superfamilies is available at http://www.biochem.ucl.ac.uk/bsm/FAM-EC/.     

The structure of a doripenem‐bound OXA‐51 class D β‐lactamase variant with enhanced carbapenemase activity

OXA‐51 is a class D β‐lactamase that is thought to be the native carbapenemase of Acinetobacter baumannii. Many variants of OXA‐51 containing active site substitutions have been identified from A. baumannii isolates, and some of these substitutions increase hydrolytic activity toward carbapenem antibiotics. We have determined the high‐resolution structures of apo OXA‐51 and OXA‐51 with one such substitution (I129L) with the carbapenem doripenem trapped in the active site as an acyl‐intermediate. The structure shows that acyl‐doripenem adopts an orientation very similar to carbapenem ligands observed in the active site of OXA‐24/40 (doripenem) and OXA‐23 (meropenem). In the OXA‐51 variant/doripenem complex, the indole ring of W222 is oriented away from the doripenem binding site, thereby eliminating a clash that is predicted to occur in wildtype OXA‐51. Similarly, in the OXA‐51 variant complex, L129 adopts a different rotamer compared to I129 in wildtype OXA‐51. This alternative position moves its side chain away from the hydroxyethyl moiety of doripenem and relieves another potential clash between the enzyme and carbapenem substrates. Molecular dynamics simulations of OXA‐51 and OXA‐51 I129L demonstrate that compared to isoleucine, a leucine at this position greatly favors a rotamer that accommodates the ligand. These results provide a molecular justification for how this substitution generates enhanced binding affinity for carbapenems, and therefore helps explain the prevalence of this substitution in clinical OXA‐51 variants.     

Using a neural network and spatial clustering to predict the location of active sites in enzymes

Structural genomics projects aim to provide a sharp increase in the number of structures of functionally unannotated, and largely unstudied, proteins. Algorithms and tools capable of deriving information about the nature, and location, of functional sites within a structure are increasingly useful therefore. Here, a neural network is trained to identify the catalytic residues found in enzymes, based on an analysis of the structure and sequence. The neural network output, and spatial clustering of the highly scoring residues are then used to predict the location of the active site.A comparison of the performance of differently trained neural networks is presented that shows how information from sequence and structure come together to improve the prediction accuracy of the network. Spatial clustering of the network results provides a reliable way of finding likely active sites. In over 69% of the test cases the active site is correctly predicted, and a further 25% are partially correctly predicted. The failures are generally due to the poor quality of the automatically generated sequence alignments.We also present predictions identifying the active site, and potential functional residues in five recently solved enzyme structures, not used in developing the method. The method correctly identifies the putative active site in each case. In most cases the likely functional residues are identified correctly, as well as some potentially novel functional groups.     

Role of active‐site residues Tyr55 and Tyr114 in catalysis and substrate specificity of Corynebacterium diphtheriae C‐S lyase

In recent years, there has been increased interest in bacterial methionine biosynthesis enzymes as antimicrobial targets because of their pivotal role in cell metabolism. C‐S lyase from Corynebacterium diphtheriae is a pyridoxal 5′‐phosphate‐dependent enzyme in the transsulfuration pathway that catalyzes the α,β‐elimination of sulfur‐containing amino acids, such as l ‐cystathionine, to generate ammonia, pyruvate, and homocysteine, the immediate precursor of L ‐methionine. In order to gain deeper insight into the functional and dynamic properties of the enzyme, mutants of two highly conserved active‐site residues, Y55F and Y114F, were characterized by UV‐visible absorbance, fluorescence, and CD spectroscopy in the absence and presence of substrates and substrate analogs, as well as by steady‐state kinetic studies. Substitution of Tyr55 with Phe apparently causes a 130‐fold decrease in at pH 8.5 providing evidence that Tyr55 plays a role in cofactor binding. Moreover, spectral data show that the mutant accumulates the external aldimine intermediate suggesting that the absence of interaction between the hydroxyl moiety and PLP‐binding residue Lys222 causes a decrease in the rate of substrate deprotonation. Mutation of Tyr114 with Phe slightly influences hydrolysis of l ‐cystathionine, and causes a change in substrate specificity towards l ‐serine and O‐acetyl‐l ‐serine compared to the wild type enzyme. These findings, together with computational data, provide useful insights in the substrate specificity of C‐S lyase, which seems to be regulated by active‐site architecture and by the specific conformation in which substrates are bound, and will aid in development of inhibitors. Proteins 2015; 83:78–90. © 2014 Wiley Periodicals, Inc.     

Functional evolution of PLP‐dependent enzymes based on active‐site structural similarities

Families of distantly related proteins typically have very low sequence identity, which hinders evolutionary analysis and functional annotation. Slowly evolving features of proteins, such as an active site, are therefore valuable for annotating putative and distantly related proteins. To date, a complete evolutionary analysis of the functional relationship of an entire enzyme family based on active‐site structural similarities has not yet been undertaken. Pyridoxal‐5′‐phosphate (PLP) dependent enzymes are primordial enzymes that diversified in the last universal ancestor. Using the comparison of protein active site structures (CPASS) software and database, we show that the active site structures of PLP‐dependent enzymes can be used to infer evolutionary relationships based on functional similarity. The enzymes successfully clustered together based on substrate specificity, function, and three‐dimensional‐fold. This study demonstrates the value of using active site structures for functional evolutionary analysis and the effectiveness of CPASS. Proteins 2014; 82:2597–2608. © 2014 Wiley Periodicals, Inc.     

Heterologous expression and periplasmic secretion of an antifungal Bacillus amyloliquefaciensBLB369 endo‐β‐1,3‐1,4‐glucanase in Escherichia coli

The endo‐β‐1,3‐1,4‐glucanases are glycoside hydrolases involved in the enzymatic depolymerization of 1,3‐1,4 β‐glucans and showed an antifungal activity against some fungi. Bacillus amyloliquefaciensBLB369 has a high antagonistic activity against phytopathogenic fungi. Its glu369 full‐coding sequence of the endo‐β‐1,3‐1,4‐glucanase gene (732 bp) was sequenced, cloned and successfully expressed in Escherichia coli Top10. The encoded protein (243 amino acids) has a calculated molecular mass of 27.3 kDa. To simplify the purification procedure, the glu369 coding sequence was cloned into the vector pKJD4. The produced OmpA‐His‐Glu369 harboured OmpA signal sequence for E. coli periplasmic localization and followed by a 6His residues for its purification. The purified His‐tagged proteins revealed two bands on SDS‐PAGE analysis with molecular masses of about 30.5 (His‐Glu369) and 32.5 kDa (OmpA‐His‐Glu369). They had the ability to inhibit the growth of phytopathogenic fungus Alternaria alternata. These favourable properties make the endo‐β‐1,3‐1,4‐glucanase a good candidate for biotechnological applications.     

Distinctive structural motifs co‐ordinate the catalytic nucleophile and the residues of the oxyanion hole in the alpha/beta‐hydrolase fold enzymes

The alpha/beta‐hydrolases (ABH) are among the largest structural families of proteins that are found in nature. Although they vary in their sequence and function, the ABH enzymes use a similar acid–base‐nucleophile catalytic mechanism to catalyze reactions on different substrates. Because ABH enzymes are biocatalysts with a wide range of potential applications, protein engineering has taken advantage of their catalytic versatility to develop enzymes with industrial applications. This study is a comprehensive analysis of 40 ABH enzyme families focusing on two identified substructures: the nucleophile zone and the oxyanion zone, which co‐ordinate the catalytic nucleophile and the residues of the oxyanion hole, and independently reported as critical for the enzymatic activity. We also frequently observed an aromatic cluster near the nucleophile and oxyanion zones, and opposite the ligand‐binding site. The nucleophile zone, the oxyanion zone and the residue cluster enriched in aromatic side chains comprise a three‐dimensional structural organization that shapes the active site of ABH enzymes and plays an important role in the enzymatic function by structurally stabilizing the catalytic nucleophile and the residues of the oxyanion hole. The structural data support the notion that the aromatic cluster can participate in co‐ordination of the catalytic histidine loop, and properly place the catalytic histidine next to the catalytic nucleophile.     

Structures of Arg‐ and Gln‐type bacterial cysteine dioxygenase homologs

In some bacteria, cysteine is converted to cysteine sulfinic acid by cysteine dioxygenases (CDO) that are only ～15–30% identical in sequence to mammalian CDOs. Among bacterial proteins having this range of sequence similarity to mammalian CDO are some that conserve an active site Arg residue (“Arg‐type” enzymes) and some having a Gln substituted for this Arg (“Gln‐type” enzymes). Here, we describe a structure from each of these enzyme types by analyzing structures originally solved by structural genomics groups but not published: a Bacillus subtilis “Arg‐type” enzyme that has cysteine dioxygenase activity (BsCDO), and a Ralstonia eutropha “Gln‐type” CDO homolog of uncharacterized activity (ReCDOhom). The BsCDO active site is well conserved with mammalian CDO, and a cysteine complex captured in the active site confirms that the cysteine binding mode is also similar. The ReCDOhom structure reveals a new active site Arg residue that is hydrogen bonding to an iron‐bound diatomic molecule we have interpreted as dioxygen. Notably, the Arg position is not compatible with the mode of Cys binding seen in both rat CDO and BsCDO. As sequence alignments show that this newly discovered active site Arg is well conserved among “Gln‐type” CDO enzymes, we conclude that the “Gln‐type” CDO homologs are not authentic CDOs but will have substrate specificity more similar to 3‐mercaptopropionate dioxygenases.     

High tolerance to simultaneous active‐site mutations in TEM‐1 β‐lactamase: Distinct mutational paths provide more generalized β‐lactam recognition

The diversity in substrate recognition spectra exhibited by various β‐lactamases can result from one or a few mutations in the active‐site area. Using Escherichia coli TEM‐1 β‐lactamase as a template that efficiently hydrolyses penicillins, we performed site‐saturation mutagenesis simultaneously on two opposite faces of the active‐site cavity. Residues 104 and 105 as well as 238, 240, and 244 were targeted to verify their combinatorial effects on substrate specificity and enzyme activity and to probe for cooperativity between these residues. Selection for hydrolysis of an extended‐spectrum cephalosporin, cefotaxime (CTX), led to the identification of a variety of novel mutational combinations. In vivo survival assays and in vitro characterization demonstrated a general tendency toward increased CTX and decreased penicillin resistance. Although selection was undertaken with CTX, productive binding (K_M) was improved for all substrates tested, including benzylpenicillin for which catalytic turnover (k_cat) was reduced. This indicates broadened substrate specificity, resulting in more generalized (or less specialized) variants. In most variants, the G238S mutation largely accounted for the observed properties, with additional mutations acting in an additive fashion to enhance these properties. However, the most efficient variant did not harbor the mutation G238S but combined two neighboring mutations that acted synergistically, also providing a catalytic generalization. Our exploration of concurrent mutations illustrates the high tolerance of the TEM‐1 active site to multiple simultaneous mutations and reveals two distinct mutational paths to substrate spectrum diversification.     

The extracellular β‐1,3‐endoglucanase EngA is involved in autolysis of Aspergillus nidulans

Aims: To elucidate the roles of the β‐1,3‐endoglucanase EngA in autolysis of the filamentous fungus Aspergillus nidulans and to identify the common regulatory elements of autolytic hydrolases. Methods and Results: A β‐1,3‐endoglucanase was purified from carbon‐starving cultures of A. nidulans. This enzyme is found to be encoded by the engA gene (locus ID: AN0472.3). Functional and gene‐expression studies demonstrated that EngA is involved in the autolytic cell wall degradation resulting from carbon starvation of the fungus. Moreover, regulation of engA is found to be dependent on the FluG/BrlA asexual sporulation signalling pathway in submerged culture. The deletion of either engA or chiB (encoding an endochitinase) caused highly reduced production of hydrolases in general. Conclusions: The β‐1,3‐endoglucanase EngA plays a pivotal role in fungal autolysis, and activities of both EngA and ChiB are necessary to orchestrate the expression of autolytic hydrolases. The production of cell wall–degrading enzymes was coordinately controlled in a highly sophisticated and complex manner. Significance and Impact of the Study: No information was available on the autolytic glucanase(s) of the euascomycete A. nidulans. This study demonstrates that EngA is a key element in fungal autolysis, and normal activities of both EngA and ChiB are crucial for balanced production of hydrolases.     

Polyglycine hydrolases: Fungal β‐lactamase‐like endoproteases that cleave polyglycine regions within plant class IV chitinases

Polyglycine hydrolases are secreted fungal proteases that cleave glycine–glycine peptide bonds in the inter‐domain linker region of specific plant defense chitinases. Previously, we reported the catalytic activity of polyglycine hydrolases from the phytopathogens Epicoccum sorghi (Es‐cmp) and Cochliobolus carbonum (Bz‐cmp). Here we report the identity of their encoding genes and the primary amino acid sequences of the proteins responsible for these activities. Peptides from a tryptic digest of Es‐cmp were analyzed by LC‐MS/MS and the spectra obtained were matched to a draft genome sequence of E. sorghi. From this analysis, a 642 amino acid protein containing a predicted β‐lactamase catalytic region of 280 amino acids was identified. Heterologous strains of the yeast Pichia pastoris were created to express this protein and its homolog from C. carbonum from their cDNAs. Both strains produced recombinant proteins with polyglycine hydrolase activity as shown by SDS‐PAGE and MALDI‐MS based assays. Site directed mutagenesis was used to mutate the predicted catalytic serine of Es‐cmp to glycine, resulting in loss of catalytic activity. BLAST searching of publicly available fungal genomes identified full‐length homologous proteins in 11 other fungi of the class Dothideomycetes, and in three fungi of the related class Sordariomycetes while significant BLAST hits extended into the phylum Basidiomycota. Multiple sequence alignment led to the identification of a network of seven conserved tryptophans that surround the β‐lactamase‐like region. This is the first report of a predicted β‐lactamase that is an endoprotease.     

A TLC bioautographic method for the detection of α‐ and β‐glucosidase inhibitors in plant extracts

Introduction – Bioautographic assays using TLC play an important role in the search for active compounds from plants. A TLC assay has previously been established for the detection of β‐glucosidase inhibitors but not for α‐glucosidase. Nonetheless, α‐glucosidase inhibition is an important target for therapeutic agents against of type 2 diabetes and anti‐viral infections. Objective – To develop a TLC bioautographic method to detect α‐ and β‐glucosidase inhibitors in plant extracts. Methodology – The enzymes α‐ and β‐d ‐glucosidase were dissolved in sodium acetate buffer. After migration of the samples, the TLC plate was sprayed with enzyme solution and incubated at room temperature for 60 min in the case of α‐d ‐glucosidase, and 37°C for 20 min in the case of β‐d ‐glucosidase. For detection of the active enzyme, solutions of 2‐naphthyl‐α‐D‐glucopyranoside or 2‐naphthyl‐β‐D‐glucopyranoside and Fast Blue Salt were mixed at a ratio of 1 : 1 (for α‐d ‐glucosidase) or 1 : 4 (for β‐d ‐glucosidase) and sprayed onto the plate to give a purple background colouration after 2–5 min. Results – Enzyme inhibitors were visualised as white spots on the TLC plates. Conduritol B epoxide inhibited α‐d ‐glucosidase and β‐d ‐glucosidase down to 0.1 µg. Methanol extracts of Tussilago farfara and Urtica dioica after migration on TLC gave enzymatic inhibition when applied in amounts of 100 µg for α‐glucosidase and 50 µg for β‐glucosidase. Conclusion – The screening test was able to detect inhibition of α‐ and β‐glucosidases by pure reference substances and by compounds present in complex matrices, such as plant extracts. Copyright © 2009 John Wiley & Sons, Ltd.     

Chromosomal integration of a synthetic expression control sequence achieves poly‐γ‐glutamate production in a Bacillus subtilis strain

Poly‐γ‐glutamate (γ‐PGA) has applications in food, medical, cosmetic, animal feed, and wastewater industries. Bacillus subtilis DB430, which possesses the γ‐PGA synthesis ywsC‐ywtAB genes in its chromosome, cannot produce γ‐PGA. An efficient synthetic expression control sequence (SECS) was introduced into the upstream region of the ywtABC genes, and this resulted in γ‐PGA‐producing B. subtilis mutant strains. Mutant B. subtilis PGA6‐2 stably produces high levels of γ‐PGA in medium A without supplementation of extra glutamic acid or ammonium chloride. The mutant B. subtilis PGA 6‐2 is not only a γ‐PGA producer, but it is also a candidate for the genetic and metabolic engineering of γ‐PGA production. © 2010 American Institute of Chemical Engineers Biotechnol. Prog., 2010     

Site‐specific inhibition of integrin αvβ3‐vitronectin association by a ser‐asp‐val sequence through an Arg‐Gly‐Asp‐binding site of the integrin

A functional proteomic technology using protein chip and molecular simulation was used to demonstrate a novel biomolecular interaction between P11, a peptide containing the Ser‐Asp‐Val (SDV) sequence and integrin αvβ3. P11 (HSDVHK) is a novel antagonistic peptide of integrin αvβ3 screened from hexapeptide library through protein chip system. An in silico docking study and competitive protein chip assay revealed that the SDV sequence of P11 is able to create a stable inhibitory complex onto the vitronectin‐binding site of integrin αvβ3. The Arg‐Gly‐Asp (RGD)‐binding site recognition by P11 was site specific because the P11 was inactive for the complex formation of a denatured form of integrin–vitronectin. P11 showed a strong antagonism against αvβ3‐GRGDSP interaction with an IC₅₀ value of 25.72±3.34 nM, whereas the value of GRGDSP peptide was 1968.73±444.32 nM. The binding‐free energies calculated from the docking simulations for each P11 and RGD peptide were ?3.99 and ?3.10 kcal/mol, respectively. The free energy difference between P11 and RGD corresponds to approximately a 4.5‐fold lower K_i value for the P11 than the RGD peptide. The binding orientation of the docked P11 was similar to the crystal structure of the RGD in αvβ3. The analyzed docked poses suggest that a divalent metal–ion coordination was a common driving force for the formation of both SDV/αvβ3 and RGD/αvβ3 complexes. This is the first report on the specific recognition of the RGD‐binding site of αvβ3 by a non‐RGD containing peptide using a computer‐assisted proteomic approach.     

Direct and indirect roles of His‐418 in metal binding and in the activity of β‐galactosidase (E. coli)

The active site of ß‐galactosidase (E. coli) contains a Mg²⁺ ion ligated by Glu‐416, His‐418 and Glu‐461 plus three water molecules. A Na⁺ ion binds nearby. To better understand the role of the active site Mg²⁺ and its ligands, His‐418 was substituted with Asn, Glu and Phe. The Asn‐418 and Glu‐418 variants could be crystallized and the structures were shown to be very similar to native enzyme. The Glu‐418 variant showed increased mobility of some residues in the active site, which explains why the substitutions at the Mg²⁺ site also reduce Na⁺ binding affinity. The Phe variant had reduced stability, bound Mg²⁺ weakly and could not be crystallized. All three variants have low catalytic activity due to large decreases in the degalactosylation rate. Large decreases in substrate binding affinity were also observed but transition state analogs bound as well or better than to native. The results indicate that His‐418, together with the Mg²⁺, modulate the central role of Glu‐461 in binding and as a general acid/base catalyst in the overall catalytic mechanism. Glucose binding as an acceptor was also dramatically decreased, indicating that His‐418 is very important for the formation of allolactose (the natural inducer of the lac operon).     

NMR characterization of an engineered domain fusion between maltose binding protein and TEM1 β‐lactamase provides insight into its structure and allosteric mechanism

RG13 is a 72 kDa engineered allosteric enzyme comprised of a fusion between maltose binding protein (MBP) and TEM1 β‐lactamase (BLA) for which maltose is a positive effector of BLA activity. We have used NMR spectroscopy to acquire [¹⁵N, ¹H]‐TROSY‐HSQC spectra of RG13 in the presence and absence of maltose. The RG13 chemical shift data was compared to the published chemical shift data of MBP and BLA. The spectra are consistent with the expectation that the individual domain structures of RG13 are substantially conserved from MBP and BLA. Differences in the spectra are consistent with the fusion geometry of MBP and BLA and the maltose‐dependent differences in the kinetics of RG13 enzyme activity. In particular, the spectra provide evidence for a maltose‐dependent conformational change of a key active site glutamate involved in deacylation of the enzyme‐substrate intermediate. Proteins 2010. © 2009 Wiley‐Liss, Inc.