PocketMatch: A new algorithm to compare binding sites in protein structures

Recognizing similarities and deriving relationships among protein molecules is a fundamental requirement in present-day biology. Similarities can be present at various levels which can be detected through comparison of protein sequences or their structural folds. In some cases similarities obscure at these levels could be present merely in the substructures at their binding sites. Inferring functional similarities between protein molecules by comparing their binding sites is still largely exploratory and not as yet a routine protocol. One of the main reasons for this is the limitation in the choice of appropriate analytical tools that can compare binding sites with high sensitivity. To benefit from the enormous amount of structural data that is being rapidly accumulated, it is essential to have high throughput tools that enable large scale binding site comparison. Results: Here we present a new algorithm PocketMatch for comparison of binding sites in a frame invariant manner. Each binding site is represented by 90 lists of sorted distances capturing shape and chemical nature of the site. The sorted arrays are then aligned using an incremental alignment method and scored to obtain PMScores for pairs of sites. A comprehensive sensitivity analysis and an extensive validation of the algorithm have been carried out. A comparison with other site matching algorithms is also presented. Perturbation studies where the geometry of a given site was retained but the residue types were changed randomly, indicated that chance similarities were virtually non-existent. Our analysis also demonstrates that shape information alone is insufficient to discriminate between diverse binding sites, unless combined with chemical nature of amino acids. Conclusion: A new algorithm has been developed to compare binding sites in accurate, efficient and high-throughput manner. Though the representation used is conceptually simplistic, we demonstrate that along with the new alignment strategy used, it is sufficient to enable binding comparison with high sensitivity. Novel methodology has also been presented for validating the algorithm for accuracy and sensitivity with respect to geometry and chemical nature of the site. The method is also fast and takes about 1/250(th) second for one comparison on a single processor. A parallel version on BlueGene has also been implemented.

Structural and functional studies of Bacillus stearothermophilus serine hydroxymethyltransferase: the role of Asn(341), Tyr(60) and Phe(351) in tetrahydrofolate binding

SHMT (serine hydoxymethyltransferase), a type I pyridoxal 5'-phosphate-dependent enzyme, catalyses the conversion of L-serine and THF (tetrahydrofolate) into glycine and 5,10-methylene THE SHMT also catalyses several THF-independent side reactions such as cleavage of P-hydroxy amino acids, trans-amination, racemization and decarboxylation. In the present study, the residues Asn(341), Tyr(60) and Phe(351), which are likely to influence THF binding, were mutated to alanine, alanine and glycine respectively, to elucidate the role of these residues in THF-dependent and -independent reactions catalysed by SHMT. The N341A and Y60A bsSHMT (Bacillus stearothermophilus SHMT) mutants were inactive for the THF-dependent activity, while the mutations had no effect on THF-independent activity. However, mutation of Phe(351) to glycine did not have any effect oil either of the activities. The crystal structures of the glycine binary complexes of the mutants showed that N341A bsSHMT forms an external aldimine as in bsSHMT, whereas Y60A and F351G bsSHMTs exist as a Mixture of internal/external aldimine and gem-diamine forms. Crystal structures of all of the three Mutants obtained in the presence of L-allo-threonine were similar to the respective glycine binary complexes. The structure of the ternary complex of F351G bsSHMT with glycine and FTHF (5-formyl THF) showed that the monoglutamate side chain of FTHF is ordered in both the subunits of the asymmetric unit, unlike in the wild-type bsSHMT. The present studies demonstrate that the residues Asn(341) and Tyr(60) are pivotal for the binding of THF/FTHF, whereas Phe(351) is responsible for the asymmetric binding of FTHF in the two subunits of the dimer.

Orchestration of Haemophilus influenzae RecJ Exonuclease by Interaction with Single-Stranded DNA-Binding Protein

RecJ exonuclease plays crucial roles in several DNA repair and recombination pathways, and its ubiquity in bacterial species points to its ancient origin and vital cellular function. RecJ exonuclease from Haemophilus influenzae is a 575-amino-acid protein that harbors the characteristic motifs conserved among RecJ homologs. The purified protein exhibits a process 5'-3' single-stranded-DNA-specific exonuclease activity. The exonuclease activity of H. influenzae RecJ (HiRecJ) was supported by Mg2+ or Mn2+ and inhibited by Cd2+ suggesting a different mode of metal binding in HiRecJ as compared to Escherichia coli RecJ (EcoRecJ). Site-directed mutagenesis of highly conserved residues in HiRecJ abolished enzymatic activity. Interestingly, substitution of alanine for aspartate 77 resulted in a catalytically inactive enzyme that bound to DNA with a significantly higher affinity as compared to the wild-type enzyme. Noticeably, steady-state kinetic studies showed that H. influenzae single-stranded DNA-binding protein (HiSSB) increased the affinity of HiRecJ for single-stranded DNA and stimulated its exonuclease activity. HiSSB, whose C-terminal tail had been deleted, failed to enhance RecJ exonuclease activity. More importantly, HiRecJ was found to directly associate with its cognate single-stranded DNA-binding protein (SSB), as demonstrated by various in vitro assays, Interaction studies carried out with the truncated variants of HiRecJ and HiSSB revealed that the two proteins interact via the C-terminus of SSB protein and the core-catalytic domain of RecJ. Taken together, these results emphasize direct interactio between RecJ and SSB, which confers functional cooperativity to these two proteins. In addition, these results implicate SSB as being involved in the recruitment of RecJ to DNA and provide insights into the interplay between these proteins in repair and recombination pathways.

Stability of Dimeric Interface in Banana Lectin: Insight from Molecular Dynamics Simulations

Banana lectin (Banlec) is a homodimeric non-glycosylated protein. It exhibits the b-prism I structure. High-temperature molecular dynamics simulations have been utilized to monitor and understand early stages of thermally induced unfolding of Banlec. The present study elucidates the behavior of the dimeric protein at four different temperatures and compares the structural and conformational changes to that of the minimized crystal structure. The process of unfolding was monitored by following the radius of gyration, the rms deviation of each residue, change in relative solvent accessibility and the pattern of inter- and intra-subunit interactions. The overall study demonstrates that the Banlec dimer is a highly stable structure, and the stability is mostly contributed by interfacial interactions. It maintains its overall conformation during high-temperature (400–500 K) simulations, with only the unstructured loop regions acquiring greater momentum under such condition. Nevertheless, at still higher temperatures (600 K) the tertiary structure is gradually lost which later extends to loss of secondary structural elements. The pattern of hydrogen bonding within the subunit and at the interface across different stages has been analyzed and has provided rationale for its intrinsic high stability.

Theoretical-Studies on the Modes of Binding of Some of the Beta-Glycosidically Linked Glucobioses to Concanavalin-A

The probable modes of binding for methyl-α-d-sophoroside, methyl-β-d-sophoroside, laminariboise and cellobiose to concanavalin A have been determined using theoretical methods. Methyl-d-sophorosides can bind to concanavalin A in two modes, i.e. by placing their reducing as well as non-reducing sugar units in the carbohydrate specific binding site, whereas laminaribiose and cellobiose can reach the binding site only with their non-reducing glucose units. However, the probability for methyl-α-d-sophoroside to bind to concanavalin A with its reducing sugar residue as the occupant of the binding site is much higher than it is with its non-reducing sugar residue as the occupant of the sugar binding site. A few of the probable conformers of methyl-β-d-sophoroside can bind to concanavalin A with either the reducing or non-reducing glucose unit. Higher energy conformers of cellobiose or laminaribiose can reach the binding site with their non-reducing residues alone. The relative differences in the binding affinities of these disaccharides are mainly due to the differences in the availability of proper conformers which can reach the binding site and to non-covalent interactions between the sugar and the protein. This study also suggests that though the sugar binding site of concanavalin A accommodates a single sugar residue, the residue outwards from the binding site also interacts with concanavalin A, indicating the existence of extended concanavalin A carbohydrate interactions.

Cholinergic Muscarinic Receptors in Human Fetal Brain: Ontogeny of [3H]Quinuclidinyl Benzilate Binding Sites in Corpus Striatum, Brainstem, and Cerebellum

The ontogeny of muscarinic receptors was studied in human fetal striatum, brainstem, and cerebellum to investigate general principles of synaptogenesis as well as the physiological balance between various chemical synapses during development in a given region of the brain. [3H]Quinuclidinyl benzilate ([-'H]QNB) binding was assayed in total particulate fraction (TPF) from various parts of brain. In the corpus striatum, QNB binding sites are present at 16 weeks of gestation (average concentration 180 fmol/mg protein of TPF), slowly increase up to 24 weeks (average concentration 217 fmol/mg protein), and rapidly increase during the third trimester to 480 fmol/mg protein of TPF. In contrast, dopaminergic receptors exist as two subpopulations. one with low affinity and the other with high affinity up to the 24th week of gestation; all of them acquire the highaffinity characteristic during the third trimester. In brainstem, the muscarinic receptors show maximum concentration by 16 weeks of age (360 fmolimg protein of TPF). Subsequently the muscarinic receptor concentration shows a gradual decline in the brainstem. In cerebellum, except for a slight increase at 24 weeks (average concentration 90 fmol/mg protein of TPF), the receptor concentration remained nearly constant at about 60-70 fmolimg protein of TPF throughout fetal life. This study demonstrates that the ontogeny of muscarinic receptors varies among the different regions, and the patterns observed suggest that receptor formation occurs principally in the third trimester. Also noteworthy is the finding that the QNB binding sites decreased in all regions of the human brain during adult life. Key Words: Cholinergic muscarinic receptors-Quinuclidinyl benzilate-Corpus striaturn-Brainstem-Cerebellum. Ravikumar B. V. and Sastry P. S. Cholinergic muscarinic receptors in human fetal brain: Ontogeny of [3H]quinuclidinyl benzilate binding sites in corpus striatum, brainstem, and cerebellum. J. Neurochem. 45, 1948- 1950 (1985).

Base specific binding of deoxyguanylate and deoxycytidylate antibodies to double stranded DNA

Antibodies raised in rabbits against daoxyguanylate and daoxycytidylate bind to 3H-(lambda) double stranded DNA and the binding is base specific. The concentrations of antibody populations that bind to double stranded DNA are much less than those binding to denatured DNA. Due to their low concentrations, these antibodies ware not detected in earlier studies. These antibodies are expected to be useful to probe the conformational flexibilities of double stranded DNAs.

Fluorescent labelling of strychnine: A novel approach for recognition of strychnine binding sites on neuronal membrane

trychnine was coupled to fluorescein isothiocyanate to mark strychnine binding sites in spinal cord of rat. Specific binding of strychnine could be demonstrated in synaptosomal fraction. Addition of glycine to the strychninised membrane led to a decrease in fluorescence indicating same receptor loci.

Isolation and characterization of binding proteins for retinol from the cytosol, nucleosol and chromatin of the oviduct magnum of laying hens

Protein fractions that bind retinol were isolated from the cytosol, nucleosol and chromatin of the oviduct magnum of laying hens. The proteins isolated from the three sources showed similar elution profiles on chromatography through Sephadex G-75 and G-50 columns, and comparable mobility during electrophoresis on sodium dodecyl sulphate/polyacrylamide gels. Their molecular weights were calculated to be around 14500. When oviducts from vitamin A-depleted and vitamin A-repleted immature chicks given oestrogen injections for 6 consecutive days were incubated with [3H]retinyl acetate, uptake of the radioactivity in the nuclei of the vitamin A-depleted tissue was severalfold higher than that in the nuclei from the vitamin A-repleted tissue.

Structural Polymorphism of the Nucleosomal DNA and Implications for Protein Binding

A nucleosome forms a basic unit of the chromosome structure. A biologically relevant question is how much of the nucleosomal conformational space is accessible to protein-free DNA, and what proportion of the nucleosomal conformations are induced by bound histones. To investigate this, we have analysed high resolution xray crystal structure datasets of DNA in protein-free as well as protein-bound forms, and compared the dinucleotide step parameters for the two datasets with those for high resolution nucleosome structures. Our analysis shows that most of the dinucleotide step parameter values for the nucleosome structures lie within the range accessible to protein-free DNA, indirectly indicating that the histone core plays more of a stabilizing role. The nucleosome structures are observed to assume smooth and nearly planar curvature, implying that ‘normal’ B-DNA like parameters can give rise to a curved geometry at the gross structural level. Different nucleosome

Small local variations in B-form DNA lead to a large variety of global geometries which can accommodate most DNA-binding protein motifs

An important question of biological relevance is the polymorphism of the double-helical DNA structure in its free form, and the changes that it undergoes upon protein-binding. We have analysed a database of free DNA crystal structures to assess the inherent variability of the free DNA structure and have compared it with a database of protein-bound DNA crystal structures to ascertain the protein-induced variations.

Crystallographic and modelling studies on Mycobacterium tuberculosis RuvA: Additional role of RuvB-binding domain and inter species variability

RuvA, along with RuvB, is involved in branch migration of heteroduplex DNA in homologous recombination. The structures of three new crystal forms of RuvA from Mycobacterium tuberculosis (MtRuvA) have been determined. The RuvB-binding domain is cleaved off in one of them. Detailed models of the complexes of octameric RuvA from different species with the Holliday junction have also been constructed. A thorough examination of the structures presented here and those reported earlier brings to light the hitherto unappreciated role of the RuvB-binding domain in determining inter-domain orientation and oligomerization. These structures also permit an exploration of the interspecies variability of structural features such as oligomerization and the conformation of the loop that carries the acidic pin, in terms of amino acid substitutions. These models emphasize the additional role of the RuvB-binding domain in Holliday junction binding. This role along with its role in oligomerization could have important biological implications.

Estrogen modulation of retinol-binding protein in immature chicks: Comparison with riboflavin carrier protein

The kinetics of estrogen (E) modulation of retinol-binding protein (RBP) production in the liver of immature chicks were compared with those governing de novo induction of riboflavin carrier protein (RCP) in the same tissue. A single dose of E markedly enhanced the plasma levels of RBP without any detectable lag period to reach peak value by 24 h and this was followed by a decline to attain the baseline by 4 days. There was no amplification of the response during secondary stimulation unlike the case with RCP induction. With multiple E administration, the 4-fold increased plasma RBP concentrations were sustained at a steady state during both primary and secondary stimulations, whereas concomitant RCP concentration progressively increased with each hormone administration and this response was further amplified during secondary stimulation. Unlike RCP induction, enhanced RBP accumulation was not strictly E dose dependent although a minimal threshold level of the steroid was required to elicit measurable response. Progesterone (P) could neither modulate nor substitute for E in enhancing plasma levels of either of the 2 proteins while the anti-estrogens, en- and zuclomifene citrate severely suppressed the production of both the proteins. RCP induction was completely inhibited by both α-amanitin and cycloheximide for prolonged periods while E-stimulated RBP production was affected only partially by α-amanitin. Likewise, cycloheximide inhibition of RBP accumulation followed a pattern similar to that of hepatic general protein synthesis.

Effect of Crowding Agents, Signal Peptide, and Chaperone SecB on the Folding and Aggregation of E. coli Maltose Binding Protein

SecB, a soluble cytosolic chaperone component of the Secexport pathway, binds to newly synthesized precursor proteins and prevents their premature aggregation and folding and subsequently targets them to the translocation machinery on the membrane. PreMBP, the precursor form of maltose binding protein, has a 26-residue signal sequence attached to the N-terminus of MBP and is a physiological substrate of SecB. We examine the effect of macromolecular crowding and SecB on the stability and refolding of denatured preMBP and MBP. PreMBP was less stable than MBP (ΔTm =7( 0.5 K) in both crowded and uncrowded solutions. Crowding did not cause any substantial changes in the thermal stability ofMBP(ΔTm=1(0.4 K) or preMBP (ΔTm=0(0.6 K), as observed in spectroscopically monitored thermal unfolding experiments. However, both MBP and preMBP were prone to aggregation while refolding under crowded conditions. In contrast to MBP aggregates, which were amorphous, preMBP aggregates form amyloid fibrils.Under uncrowded conditions, a molar excess of SecB was able to completely prevent aggregation and promote disaggregation of preformed aggregates of MBP. When a complex of the denatured protein and SecB was preformed, SecB could completely prevent aggregation and promote folding of MBP and preMBP even in crowded solution. Thus, in addition to maintaining substrates in an unfolded, export-competent conformation, SecB also suppresses the aggregation of its substrates in the crowded intracellular environment. SecB is also able to promote passive disaggregation of macroscopic aggregates of MBP in the absence of an energy source such as ATP or additional cofactors. These experiments also demonstrate that signal peptide can reatly influence protein stability and aggregation propensity.

Copper binding by Pseudomonas aeruginosa

A copper-binding complex formed in the exopolysaccharide fraction of Image was isolated and characterized using a variety of techniques. By comparison with model Cu(II) complexes of uronic acids, it is shown that the Image forms a square-planer, cupric complex similar to cupric glucuronates.