Clues to reaction specificity in PLP-dependent fold type I aminotransferases of monosaccharide biosynthesis.

Novel functions can emerge in an enzyme family while conserving catalytic mechanism, motif or fold. Pyridoxal 5'-phosphate-dependent enzymes have evolved into seven fold-types and catalyze diverse reactions using the same mechanism for the formation of external aldimine. Nucleotide sugar aminotransferases (which will be henceforth referred to as aminotransferases) belong to fold type I and mediate the biosynthesis of several monosaccharides. They use diverse substrates but are highly selective to the C3 or C4 carbon to which amine group is transferred. Profile hidden Markov models (HMMs) were able to identify aminotransferases but could not capture reaction specificity. A search for discriminating features led to the discovery of sequence motifs that are located near the pyranose binding site suggesting their role in imparting reaction specificity. Using a position weight matrix for this motif, we were able to assign reaction specificity to a large number of aminotransferases. Inferences from this analysis set way for future experiments that can shed light on mechanisms of functional diversification in nucleotide sugar aminotransferases of fold type I.

Monosaccharide biosynthesis pathways database.

A distinctive feature of glycans vis-à-vis proteins and nucleic acids is its structural complexity, which arises from the huge repertoire of monosaccharides, isomeric linkages and branching. A very large number of monosaccharides have so far been discovered in natural glycans. Experimentally, pathways for the biosynthesis have been characterized completely for 55 monosaccharides and partially for a few more. However, there is no single platform, which provides information about monosaccharide biosynthesis pathways and associated enzymes We have gathered 572 experimentally characterized enzymes of 66 biosynthesis pathways from literature and set up a first of its kind database called the Monosaccharide Biosynthesis Pathways Database http://www.bio.iitb.ac.in/mbpd/). Annotations such as the reaction catalyzed, substrate specificity, biosynthesis pathway and PubMed IDs are provided for all the enzymes in the database. Sequence homologs of the experimentally characterized enzymes found in nearly 13,000 completely sequenced genomes from Bacteria and Archaea have also been included in the database. This platform will help in the deduction of evolutionary relationships among enzymes such as aminotransferases, nucleotidyltransferases, acetyltransferases and SDR family enzymes. It can also facilitate experimental studies such as direct enzyme assays to validate putative annotations, establish structure-function relationship, expression profiling to determine the function, determine the phenotypic consequences of gene knock-out/knock-in and complementation studies.

Characterization of left-handed beta helix-domains, and identification and functional annotation of proteins containing such domains.

Only about 0.3% of the entries in UniProt database have manually curated annotation. Annotation at the molecular level often relies on low-throughput one-protein-at-a-time approach. Computational methods bridge this gap by assigning function based on sequence and/or fold similarity. Left-handed beta helix (LbH) consists of three repeating six-stranded beta-strands forming an 18-mer turn of the helix. Analysis of LbH-domains showed that variations are found in the number of residues in a beta-strand (5-7, 6 being the most common), number of turns (4-10) of the helix, insertions of one or more loops of variable length (0-36 residues), and the location of loop insertion. An 18-mer HMM profile was created which identifies LbH-domain containing proteins using sequence as the only input; the number of false positives is zero when proteins tested were those with known 3D structures. 136 474 entries of TrEMBL database were found to contain LbH-domain. Rules developed by analyzing LbH-domain containing acyltransferases, gamma-class carbonic anhydrases, and nucleotidyltransferases have led to the annotation of 17 389 TrEMBL entries which currently have no functional tag.

EpsN from Bacillus subtilis 168 has UDP-2,6-dideoxy 2-acetamido 4-keto glucose aminotransferase activity in vitro.

The gene epsN of Bacillus subtilis 168 was cloned and overexpressed in Escherichia coli. Purified recombinant EpsN is shown to be a pyridoxal 5'-phosphate (PLP)-dependent aminotransferase by absorption spectroscopy, l-cycloserine inhibition and reverse phase HPLC studies. EpsN catalyzes the conversion of UDP-2,6-dideoxy 2-acetamido 4-keto glucose to UDP-2,6-dideoxy 2-acetamido 4-amino glucose. Lys190 was found by sequence comparison and site-directed mutagenesis to form Schiff base with PLP. Mutagenesis studies showed that, in addition to Lys190, Ser185, Glu164, Gly58 and Thr59 are essential for aminotransferase activity.

EpsM from Bacillus subtilis 168 has UDP-2,4,6-trideoxy-2-acetamido-4-amino glucose acetyltransferase activity in vitro.

Bacillus subtilis 168 EpsM (UniProt id P71063) has been electronically annotated as putative acetyltransferase in the UniProt database. The gene epsM was cloned and overexpressed in E. coli with an N-terminal GST tag. The purified fusion protein was shown by absorption spectroscopy, autoradiography and reverse phase HPLC to catalyse the conversion of UDP-2,4,6-trideoxy-2-acetamido-4-amino glucose to UDP-2,4,6-trideoxy-2,4-diacetamido glucose, commonly known as N,N'-diacetylbacillosamine, using acetyl coenzyme A as the donor substrate. His146 was shown by site-directed mutagenesis to be essential for acetyltransferase activity. It is hypothesized that EpsC (NAD+ dependent UDP GlcNAc 4,6-dehydratase), EpsN (PLP dependent aminotransferase) and EpsM, all of which are part of the eps operon, are involved in the biosynthesis of N,N'-diacetylbacillosamine.

In vitro characterization of N-terminal truncated EpsC from Bacillus subtilis 168, a UDP-N-acetylglucosamine 4,6-dehydratase.

Bacillus subtilis 168 EpsC is annotated as "Probable polysaccharide biosynthesis protein" in the SwissProt database. epsC is part of the eps operon, thought to be involved in the biosynthesis of exopolymeric substances (EPS). The present study was undertaken to determine the molecular function of EpsC. Sequence analysis of EpsC suggested the presence of a transmembrane domain. Two N-terminal deletion mutants in which residues 1-89 (EpsC89) and 1-115 (EpsC115) are deleted were cloned and overexpressed. Enzyme activity and substrate preferences were investigated by reverse phase HPLC, surface plasmon resonance (SPR) spectroscopy and absorption spectroscopy. These data show that EpsC has UDP-GlcNAc 4,6-dehydratase activity in vitro. Purified recombinant proteins were found to utilise UDP-Glc and TDP-Glc also as substrates. In addition, EpsC115 could utilise UDP-Gal and UDP-GalNAc as substrates whereas EpsC89 could only bind these two sugar nucleotides. These results show that deletion of a longer N-terminal region broadens substrate specificity. These broadened specificity is perhaps an outcome of the deletion of the putative transmembrane domain and may not be present in vivo. EpsC, together with the aminotransferase EpsN (Kaundinya CR et al., Glycobiology, 2018) and acetyltransferase EpsM (unpublished data), appears to be involved in the biosynthesis of N,N'-diacetylbacillosamine.

Molecular events during the early stages of aggregation of GNNQQNY: An all atom MD simulation study of randomly dispersed peptides.

This study probes the early events during lag phase of aggregation of GNNQQNY using all atom MD simulations in explicit solvent. Simulations were performed by varying system size, temperature and starting configuration. Peptides dispersed randomly in the simulation box come together early on in the simulation and form aggregates. These aggregates are dynamic implying the absence of stabilizing interactions. This facilitates the exploration of alternate arrangements. The constituent peptides sample a variety of conformations, frequently re-orient and re-arrange with respect to each other and dissociate from/re-associate with the aggregate. The size and lifetime of aggregates vary depending upon the number of inter-peptide backbone H-bonds. Most of the aggregates formed are amorphous but crystalline aggregates of smaller size (mainly 2-mers) do appear and sustain for varying durations of time. The peptides in crystalline 2-mers are mostly anti-parallel. The largest crystalline aggregate that appears is a 4-mer in a single sheet and a 4-, 5-, or 6-mer in double layered arrangement. Crystalline aggregates grow either by the sequential addition of peptides, or by the head-on or lateral collision-adhesion of 2-mers. The formation of various smaller aggregates suggests the polymorphic nature of oligomers and heterogeneity in the lag phase.

Size, orientation and organization of oligomers that nucleate amyloid fibrils: clues from MD simulations of pre-formed aggregates.

All-atom MD simulations of pre-formed aggregates of GNNQQNY with variable size (5 to 16 peptides), orientation (parallel or anti-parallel), organization (single or double sheet, with or without twist), charge status of termini and temperature (300 and 330K) have been performed for 50ns each (68 simulations; total time=3.4µs). Double-layer systems are stable irrespective of whether the peptides within the sheet are oriented parallel or anti-parallel. The lifetime of single sheet systems is determined by the protonation status, nature of association of peptides and the size of the aggregates. For example, single sheet 8-mers are stable with parallel arrangement and neutral termini, or with anti-parallel arrangement and charged termini. This suggests that the residues flanking the amyloidogenic sequence also play an important role in determining the organization of peptides in an aggregate. Twist of the cross-beta sheets is found to be intrinsic to the aggregates. Main chain H-bonds are key determinants of stability and loss of these H-bonds is followed by disorder and/or dissociation of the peptide despite the presence of side chain hydrogen bonds. Aggregates are inherently asymmetric along the fiber axis and dissociation from the C-edge is observed more often. An aggregate can disintegrate into smaller-sized oligomers or the edge peptides can dissociate sequentially. A variety of dissociation and disintegration events are observed pointing to the existence of multiple pathways for association during nucleation. It appears that a heterogeneous mixture of oligomers of different sizes exist prior to the formation of the critical nucleus.

Conformational mapping and energetics of saccharide-aromatic residue interactions: implications for the discrimination of anomers and epimers and in protein engineering.

Aromatic residues play a key role in saccharide-binding sites. Experimental studies have given an estimate of the energetics of saccharide-aromatic residue interactions. In this study, dependence of the energetics on the mutual position-orientation (PO) of saccharide and aromatic residue has been investigated by geometry optimization of a very large number (164) of complexes at MP2/6-31G(d,p) level of theory. The complexes are of Tyr and Phe analogs with α/ß-D-Glc, ß-D-Gal, α-D-Man and α/ß-L-Fuc. A number of iso-energy POs are found for the complexes of all six saccharides. Stacking and non-stacking modes of binding are found to be of comparable strengths. In general, complexes of p-OHTol are stronger than those of Tol, and those dominated by OH···O interactions are more stable than ones dominated by CH···π interactions. The strengths of OH···O/π interactions, but not those of CH···π, show large variations. Even though an aromatic residue has a large variety of POs to interact with a saccharide, distinct preferences are found due to anomeric and epimeric differences. An aromatic residue can interact from either the a- or b-face of Glc, but only through the b-face with Gal, its C4-epimer. In contrast, stacking interaction with Man (C2-epimer of Glc) requires the participation of the -CH(2)OH group and free rotation of this group, as is observed in solution, precludes all modes of stacking interactions. It is also found that an aromatic residue can be strategically placed either to discriminate or to accommodate (i) anomers of Glc and of Fuc and (ii) Gal/Fuc. Thus, analysis of the optimized geometries of by far the largest number of complexes, and with six different saccharides, at this level of theory has given insights into how Nature cleverly uses aromatic residues to fine tune saccharide specificities of proteins. These are of immense utility for protein engineering and protein design studies.

E93R substitution of Escherichia coli FtsZ induces bundling of protofilaments, reduces GTPase activity, and impairs bacterial cytokinesis.

Recently, we found that divalent calcium has no detectable effect on the assembly of Mycobacterium tuberculosis FtsZ (MtbFtsZ), whereas it strongly promoted the assembly of Escherichia coli FtsZ (EcFtsZ). While looking for potential calcium binding residues in EcFtsZ, we found a mutation (E93R) that strongly promoted the assembly of EcFtsZ. The mutation increased the stability and bundling of the FtsZ protofilaments and produced a dominating effect on the assembly of the wild type FtsZ (WT-FtsZ). Although E93R-FtsZ was found to bind to GTP similarly to the WT-FtsZ, it displayed lower GTPase activity than the WT-FtsZ. E93R-FtsZ complemented for its wild type counterpart as observed by a complementation test using JKD7-1/pKD3 cells. However, the bacterial cells became elongated upon overexpression of the mutant allele. We modeled the structure of E93R-FtsZ using the structures of MtbFtsZ/Methanococcus jannaschi FtsZ (MjFtsZ) dimers as templates. The MtbFtsZ-based structure suggests that the Arg(93)-Glu(138) salt bridge provides the additional stability, whereas the effect of mutation appears to be indirect (allosteric) if the EcFtsZ dimer is similar to that of MjFtsZ. The data presented in this study suggest that an increase in the stability of the FtsZ protofilaments is detrimental for the bacterial cytokinesis.

Quantification of binding affinities of essential sugars with a tryptophan analogue and the ubiquitous role of C-H···π interactions.

The role of noncovalent interactions in carbohydrate recognition by aromatic amino acids has long been reported. To develop a molecular understanding of noncovalent interactions in the recognition process, we have examined a series of binary complexes between 3-methylindole (3-MeIn) and sugars. In particular, the geometries and binding affinities of 3-MeIn with α/ß-D-glucose, ß-D-galactose, α-D-mannose and α/ß-L-fucose are obtained using the MP2(full)/6-31G(d,p) and the M06/TZV2D//MP2/6-31G(d,p) level of theories. The conventional hydrogen bonding such as N-H···O and C-H···O as well as nonconventional O-H···π and C-H···π type of interactions is, in general, identified as responsible for the moderately strong interaction energies. Large variations in the position-orientations of 3-MeIn with respect to saccharide are noticed, within the same sugar family, as well as across different sugar series. Furthermore, complexes with large differences in their geometries are recognized as capable of exhibiting very similar interaction energies, underscoring the significance of exhaustive conformation sampling, as carried out in the present study. These observations are readily attributed to the differences in the efficiency of the type of interactions enlisted above. The highest and lowest interaction energies, upon inclusion of 50% BSSE correction, are found to be -16.02 and -6.22 kcal mol(-1), respectively, for α-D-glucose (1a) and α-L-fucose (5j). While more number of prominent conventional hydrogen bonding contacts remains as a characteristic feature of the strongly bound complexes, the lower end of the interaction energy spectrum is dominated by multiple C-H···π interactions. The complexes exhibiting as many as four C-H···π contacts are identified in the case of α/ß-D-glucose, ß-D-galactose, and α/ß-L-fucose with an interaction energy hovering around -8 kcal mol(-1). The presence of effective C-H···π interactions is found to be dependent on the saccharide configuration as well as the area of the apolar patch constituted by the C-H groups. The study offers a comprehensive set of binary complexes, across different saccharides, which serves as an illustration of the significance and ubiquitous nature of C-H···π interactions in carbohydrate binding in saccharide-protein complexes.

Griseofulvin stabilizes microtubule dynamics, activates p53 and inhibits the proliferation of MCF-7 cells synergistically with vinblastine.

BACKGROUND: Griseofulvin, an antifungal drug, has recently been shown to inhibit proliferation of various types of cancer cells and to inhibit tumor growth in athymic mice. Due to its low toxicity, griseofulvin has drawn considerable attention for its potential use in cancer chemotherapy. This work aims to understand how griseofulvin suppresses microtubule dynamics in living cells and sought to elucidate the antimitotic and antiproliferative action of the drug. METHODS: The effects of griseofulvin on the dynamics of individual microtubules in live MCF-7 cells were measured by confocal microscopy. Immunofluorescence microscopy, western blotting and flow cytometry were used to analyze the effects of griseofulvin on spindle microtubule organization, cell cycle progression and apoptosis. Further, interactions of purified tubulin with griseofulvin were studied in vitro by spectrophotometry and spectrofluorimetry. Docking analysis was performed using autodock4 and LigandFit module of Discovery Studio 2.1. RESULTS: Griseofulvin strongly suppressed the dynamic instability of individual microtubules in live MCF-7 cells by reducing the rate and extent of the growing and shortening phases. At or near half-maximal proliferation inhibitory concentration, griseofulvin dampened the dynamicity of microtubules in MCF-7 cells without significantly disrupting the microtubule network. Griseofulvin-induced mitotic arrest was associated with several mitotic abnormalities like misaligned chromosomes, multipolar spindles, misegregated chromosomes resulting in cells containing fragmented nuclei. These fragmented nuclei were found to contain increased concentration of p53. Using both computational and experimental approaches, we provided evidence suggesting that griseofulvin binds to tubulin in two different sites; one site overlaps with the paclitaxel binding site while the second site is located at the alphabeta intra-dimer interface. In combination studies, griseofulvin and vinblastine were found to exert synergistic effects against MCF-7 cell proliferation. CONCLUSIONS: The study provided evidence suggesting that griseofulvin shares its binding site in tubulin with paclitaxel and kinetically suppresses microtubule dynamics in a similar manner. The results revealed the antimitotic mechanism of action of griseofulvin and provided evidence suggesting that griseofulvin alone and/or in combination with vinblastine may have promising role in breast cancer chemotherapy.

The glycan alphabet is not universal: a hypothesis.

Several monosaccharides constitute naturally occurring glycans, but it is uncertain whether they constitute a universal set like the alphabets of proteins and DNA. Based on the available experimental observations, it is hypothesized herein that the glycan alphabet is not universal. Data on the presence/absence of pathways for the biosynthesis of 55 monosaccharides in 12 939 completely sequenced archaeal and bacterial genomes are presented in support of this hypothesis. Pathways were identified by searching for homologues of biosynthesis pathway enzymes. Substantial variations were observed in the set of monosaccharides used by organisms belonging to the same phylum, genera and even species. Monosaccharides were grouped as common, less common and rare based on their prevalence in Archaea and Bacteria. It was observed that fewer enzymes are sufficient to biosynthesize monosaccharides in the common group. It appears that the common group originated before the formation of the three domains of life. In contrast, the rare group is confined to a few species in a few phyla, suggesting that these monosaccharides evolved much later. Fold conservation, as observed in aminotransferases and SDR (short-chain dehydrogenase reductase) superfamily members involved in monosaccharide biosynthesis, suggests neo- and sub-functionalization of genes led to the formation of the rare group monosaccharides. The non-universality of the glycan alphabet begets questions about the role of different monosaccharides in determining an organism's fitness.

Characterization of the conformational and orientational dynamics of ganglioside GM1 in a dipalmitoylphosphatidylcholine bilayer by molecular dynamics simulations.

The structure and dynamics of a single GM1 (Gal5-beta1,3-GalNAc4-beta1,4-(NeuAc3-alpha2,3)-Gal2-beta1,4-Glc1-beta1,1-Cer) embedded in a DPPC bilayer have been studied by MD simulations. Eleven simulations, each of 10 ns productive run, were performed with different initial conformations of GM1. Simulations of GM1-Os in water and of a DPPC bilayer were also performed to delineate the effects of the bilayer and GM1 on the conformational and orientational dynamics of each other. The conformation of the GM1 headgroup observed in the simulations is in agreement with those reported in literature; but the headgroup is restricted when embedded in the bilayer. NeuAc3 is the outermost saccharide towards the water phase. Glc1 and Gal2 prefer a parallel, and NeuAc3, GalNac4 and Gal5 prefer a perpendicular, orientation with respect to the bilayer normal. The overall characteristics of the bilayer are not affected by the presence of GM1; however, GM1 does influence the DPPC molecules in its immediate vicinity. The implications of these observations on the specific recognition and binding of GM1 embedded in a lipid bilayer by exogenous proteins as well as proteins embedded in lipids have been discussed.

Characterization of symmetric and asymmetric lipid bilayers composed of varying concentrations of ganglioside GM1 and DPPC.

Gangliosides are a group of structurally diverse, sialic acid containing glycosphingolipids embedded into the membrane via their hydrophobic ceramide moiety. To gain atomic level insights into the structural perturbations caused by Galbeta3GalNAcbeta4(NeuAcalpha3)Galbeta4Glc1Cer (GM1), molecular dynamics (MD) simulations of a 1,2-dipalmitoyl-sn-glycero-3-phosphatidylcholine (DPPC) bilayer containing GM1 at five different concentrations have been performed. Biological membranes contain GM1 only on the exoplasmic leaflet. However, vesicles prepared in the laboratory contain GM1 in both the leaflets albeit unequally. Hence, simulations were performed with GM1 present in only one (asymmetric bilayers) or in both of the leaflets (symmetric bilayers) of the bilayer. In symmetric bilayers, there is a decrease in surface area, an increase in deuterium order parameter, and an increase in peak-to-peak distance of DPPC with increasing concentration of GM1. Thus, the overall area of the lipid bilayer decreases (condensation effect) and the thickness increases with increasing concentrations of GM1. Even in asymmetric systems, decrease in surface area and increase in deuterium order parameter of hydrocarbon chains of DPPC are observed. However, the decrease in bilayer area and the increase in bilayer thickness are not as much as in the symmetric bilayer.

Fold-recognition and comparative modeling of human beta3GalT I, II, IV, V and VI and beta3GalNAcT I: prediction of residues conferring acceptor substrate specificity.

beta3GalTs are type II transmembrane proteins that transfer galactose from UDP-Gal donor substrate to acceptor GlcNAc, GalNAc or Gal in beta1-->3-linkage. beta1-->3-linked galactose have been found to be a part of many glycans like glycosphingolipids, core tetrasaccharide of proteoglycans, type 1 chains. The 3-D structure of none of the beta3GalTs is known to date. In this study, the 3-D structures of human beta3GalT I, II, IV, V, VI and beta3GalNAcT I have been modeled using fold-recognition and comparative modeling methods. Residues that constitute the UDP-Gal binding site have been predicted. The models are able to qualitatively rationalize data from the site-directed mutagenesis experiments reported in the literature. Residues likely to be involved in conferring differential acceptor substrate specificity have been predicted by a combination of specificity determining positions prediction (SDPs) and subsequent mapping on the generated 3-D models.

Length and composition analysis of the cytoplasmic, transmembrane and stem regions of human Golgi glycosyltransferases.

A dataset of experimentally characterized, human Golgi GlyTs with type II membrane topology was created. Based on the experimentally observed acceptor substrate preferences, the GlyTs were classified into five functional categories: biosynthesis of blood group antigens, glycolipids, N-glycans, O-glycans and glycosaminoglycans. The cytoplasmic, transmembrane and stem (CTS) regions were predicted and their length and composition were analyzed. The stem region of GlyTs involved in the biosynthesis of glycolipids and blood group antigens appear to have a shorter stem region compared to those GlyTs which participate in the biosynthesis of N- and O-linked glycans and glycosaminoglycans. The stem regions of all the GlyTs, irrespective of the functional category to which they belong, were found to be rich in disorder-promoting amino acid residues. Thus, the stem region is largely devoid of any regular secondary structure thereby facilitating its tethering role. A higher frequency of occurrence of basic amino acids is observed towards the N-terminus of the transmembrane domain and this is suggested to be important for topogenesis of these enzymes.

Rv3634c from Mycobacterium tuberculosis H37Rv encodes an enzyme with UDP-Gal/Glc and UDP-GalNAc 4-epimerase activities.

A bioinformatics study revealed that Mycobacterium tuberculosis H37Rv (Mtb) contains sequence homologs of Campylobacter jejuni protein glycosylation enzymes. The ORF Rv3634c from Mtb was identified as a sequence homolog of C. jejuni UDP-Gal/GalNAc 4-epimerase. This study reports the cloning of Rv3634c and its expression as an N-terminal His-tagged protein. The recombinant protein was shown to have UDP-Gal/Glc 4-epimerase activity by GOD-POD assay and by reverse phase HPLC. This enzyme was shown to have UDP-GalNAc 4-epimerase activity also. Residues Ser121, Tyr146 and Lys150 were shown by site-directed mutagenesis to be important for enzyme activity. Mutation of Ser121 and Tyr146 to Ala and Phe, respectively, led to complete loss of activity whereas mutation of Lys150 to Arg led to partial loss of activity. There were no gross changes in the secondary structures of any of these three mutants. These results suggest that Ser121 and Tyr146 are essential for epimerase activity of Rv3634c. UDP-Gal/Glc 4-epimerases from other organisms also have a catalytic triad consisting of Ser, Tyr and Lys. The triad carries out proton transfer from nucleotide sugar to NAD+ and back, thus effecting the epimerization of the substrate. Addition of NAD+ to Lys150 significantly abrogates the loss of activity, suggesting that, as in other epimerases, NAD+ is associated with Rv3634c.

Fold-recognition and comparative modeling of human alpha2,3-sialyltransferases reveal their sequence and structural similarities to CstII from Campylobacter jejuni.

BACKGROUND: The 3-D structure of none of the eukaryotic sialyltransferases (SiaTs) has been determined so far. Sequence alignment algorithms such as BLAST and PSI-BLAST could not detect a homolog of these enzymes from the protein databank. SiaTs, thus, belong to the hard/medium target category in the CASP experiments. The objective of the current work is to model the 3-D structures of human SiaTs which transfer the sialic acid in alpha2,3-linkage viz., ST3Gal I, II, III, IV, V, and VI, using fold-recognition and comparative modeling methods. The pair-wise sequence similarity among these six enzymes ranges from 41 to 63%. RESULTS: Unlike the sequence similarity servers, fold-recognition servers identified CstII, a alpha2,3/8 dual-activity SiaT from Campylobacter jejuni as the homolog of all the six ST3Gals; the level of sequence similarity between CstII and ST3Gals is only 15-20% and the similarity is restricted to well-characterized motif regions of ST3Gals. Deriving template-target sequence alignments for the entire ST3Gal sequence was not straightforward: the fold-recognition servers could not find a template for the region preceding the L-motif and that between the L- and S-motifs. Multiple structural templates were identified to model these regions and template identification-modeling-evaluation had to be performed iteratively to choose the most appropriate templates. The modeled structures have acceptable stereochemical properties and are also able to provide qualitative rationalizations for some of the site-directed mutagenesis results reported in literature. Apart from the predicted models, an unexpected but valuable finding from this study is the sequential and structural relatedness of family GT42 and family GT29 SiaTs. CONCLUSION: The modeled 3-D structures can be used for docking and other modeling studies and for the rational identification of residues to be mutated to impart desired properties such as altered stability, substrate specificity, etc. Several studies in literature have focused on the development of tools and/or servers for the large-scale/automated modeling of 3-D structures of proteins. In contrast, the present study focuses on modeling the 3-D structure of a specific protein of interest to a biochemist and illustrates the associated difficulties. It is also able to establish a sequence/structure relationship between sialyltransferases of two distinct families.

Beta-hairpins with native-like and non-native hydrogen bonding patterns could form during the refolding of staphylococcal nuclease.

Refolding of staphylococcal nuclease has been studied recently by hydrogen-deuterium exchange and NMR spectroscopy. These studies infer that beta-hairpin formed by strand 2 and strand 3 connected by reverse turn forms early during the refolding of nuclease. Typically, hydrogen-deuterium exchange NMR techniques are usually carried out on a time scale of milliseconds whereas beta-hairpins are known to fold on a much shorter time scale. It follows that in the experiments, the hydrogen-deuterium exchange protection patterns could be arising from a significant population of fully formed hairpins. In order to demonstrate it is the fully formed hairpins which gives rise to the hydrogen-deuterium exchange protection patterns, we have considered molecular dynamics simulation of the peptide (21)DTVKLMYKGQPMTFR(35) from staphylococcal nuclease corresponding to the beta-hairpin region, using GROMOS96 force field under NVT conditions. Starting from unfolded conformational states, the peptide folds into hairpin conformations with native-like and non-native hydrogen bonding patterns. Subsequent to folding, equilibrium conditions prevail. The computed protection factors and atom depth values, at equilibrium, of the various amide protons agree qualitatively with experimental observations. A collection of molecules following the trajectories observed in the simulations can account for experimental observations. These simulations provide a molecular picture of the formed hairpins and their conformational features during the refolding experiments on nuclease, monitored by hydrogen-deuterium exchange.