Search PubMed⌕ Search

Biomedical subjects

S Subbiah

Publications and source records attributed to S Subbiah.

At least 19 recordsLinked to original sources

SEGE: A database on 'intron less/single exonic' genes from eukaryotes.

UNLABELLED: Eukaryotes have both 'intron containing' and 'intron less' genes. Several databases are available for 'intron containing' genes in eukaryotes. In this note, we describe a database for 'intron less' genes from eukaryotes. 'Intron less' eukaryotic genes having prokaryotic architecture will help to understand gene evolution in a much simpler way unlike 'intron containing' genes. AVAILABILITY: SEGE is available at http://intron.bic.nus.edu.sg/seg/ CONTACT: mmeena@ntu.edu.sg

Animals↗

APBioNet: the Asia-Pacific regional consortium for bioinformatics.

Bioinformatics and computational biology, along with the related fields of genomics, proteomics, functional genomics and systems biology are new wave scientific disciplines that harness composite computational power across networks to advance biological knowledge at the most basic level and to direct traditional laboratory-based research efforts in the biomedical sciences. 'Fostering the growth of bioinformatics and allied disciplines in the Asia-Pacific region' is the motto of the first regional bioinformatics society, the Asia-Pacific Bioinformatics Network (APBioNet). APBioNet addresses the issues of hardware, software, databases and networks pertaining to bioinformatics, with the additional layer of pertinent education, training and research. Recent milestones achieved include hosting an international bioinformatics symposium in Asia and setting up large-scale regional grid-computing projects.

Asia↗

Endothelial function in patients with chronic heart failure.

Congestive heart failure is characteristically associated with systemic vasoconstriction and impaired vasodilatory capacity, leading to decreased peripheral perfusion. Factors identified as possible causes of reduced vasodilatory capacity include activation of the sympathetic nervous system and the renin-angiotensin system, vascular stiffness due to increased sodium and fluid retention, and structural vascular changes. More recently, the role of the endothelium as a mediator of vasoregulation and tissue perfusion has been recognized.

Endothelium, Vascular↗

Validation of NMR side-chain conformations by packing calculations.

The precision and accuracy of protein structures determined by nuclear magnetic resonance (NMR) spectroscopy depend on the completeness of input experimental data set. Typically, rather than a single structure, an ensemble of up to 20 equally representative conformers is generated and routinely deposited in the Protein Database. There are substantially more experimentally derived restraints available to define the main-chain coordinates than those of the side chains. Consequently, the side-chain conformations among the conformers are more variable and less well defined than those of the backbone. Even when a side chain is determined with high precision and is found to adopt very similar orientations among all the conformers in the ensemble, it is possible that its orientation might still be incorrect. Thus, it would be helpful if there were a method to assess independently the side-chain orientations determined by NMR. Recently, homology modeling by side-chain packing algorithms has been shown to be successful in predicting the side-chain conformations of the buried residues for a protein when the main-chain coordinates and sequence information are given. Since the main-chain coordinates determined by NMR are consistently more reliable than those of the side-chains, we have applied the side-chain packing algorithms to predict side-chain conformations that are compatible with the NMR-derived backbone. Using four test cases where the NMR solution structures and the X-ray crystal structure of the same protein are available, we demonstrate that the side-chain packing method can provide independent validation for the side-chain conformations of NMR structures. Comparison of the side-chain conformations derived by side-chain packing prediction and by NMR spectroscopy demonstrates that when there is agreement between the NMR model and the predicted model, on average 78% of the time the X-ray structure also concurs. While the side-chain packing method can confirm the reliable residue conformations in NMR models, more importantly, it can also identify the questionable residue conformations with an accuracy of 60%. This validation method can serve to increase the confidence level for potential users of structural models determined by NMR.

Animals↗

BioJAKE: a tool for the creation, visualization and manipulation of metabolic pathways.

The BioJAKE program has been created for the visualization, creation and manipulation of metabolic pathways. It has been designed to provide a familiar and easy-to-use interface while still allowing for the input and manipulation of complex and detailed metabolic data. In recognition of the detailed and diverse sources of data available across the Internet, it also provides a mechanism by which remote database queries can be stored and performed with respect to individual molecules within a pathway. This remote database access functionality is offered in addition to comprehensive local database creation, management and querying capability. The program has been developed in Java so as to provide for platform independence and maximum extendibility.

Biochemistry↗

In vivo expression of an alternatively spliced human tumor message that encodes a truncated form of cathepsin B. Subcellular distribution of the truncated enzyme in COS cells.

Cathepsin B is a lysosomal cysteine protease whose increased expression is believed to be linked to the malignant progression of tumors. Alternative splicing and the use of alternative transcription initiation sites in humans produce cathepsin B mRNAs that differ in their 5'- and 3'-untranslated ends. Some human tumors also contain cathepsin B-related transcripts that lack exon 3 which encodes the N-terminal signal peptide and 34 of the 62-amino acid inhibitory propeptide. In this study we show that one such transcript, CB(-2,3), which is missing exons 2 and 3, is likely to be a functional message in tumors. Thus, CB(-2,3) was found to be otherwise complete, containing the remainder of the cathepsin B coding sequence and the part of the 3'-untranslated region that is common to all previously characterized cathepsin B mRNAs in humans. Its in vitro translation product can be folded to produce enzymatic activity against the cathepsin B-specific substrate, Nalpha-benzyloxycarbonyl-L-Arg-L-Arg-4-methylcoumaryl-7-amide. Endogenous CB(-2,3) from the metastatic human melanoma cell line, A375M, co-sediments with polysomes, indicating that it engages the eukaryotic translation machinery in these cells. Epitope-tagged forms of the truncated cathepsin B from CB(-2,3) are produced in amounts comparable to the normal protein after transient transfection into COS cells. Immunofluorescence microscopy and subcellular fractionation show this novel tumor form of cathepsin B to be associated with nuclei and other membranous organelles, where it is likely to be bound to the cytoplasmic face of the membranes. This subcellular distribution was different from the lysosomal pattern shown by the epitope-tagged, full-length cathepsin B in COS cells. These results indicate that the message missing exons 2 and 3 is likely to be translated into a catalytically active enzyme, and that alternative splicing (exon skipping) could contribute to the aberrant intracellular trafficking of cathepsin B that is observed in some human cancers.

Alternative Splicing↗

Integration of bioInformatics tools at the National University of Singapore (NUS).

In the past decade "Big Science" such as the Genome Project has generated an enormous amount of data in the life sciences. Concurrently, the synergy of this project with existing research has quickened the pace of biological discovery. But the major drawback that is beginning to be felt worldwide is the primitive level of organisation in the data accumulated. Without a proper framework or knowledge scaffold to hang and interconnect the various bits of data and information, the national knowledge-to-data ratio is declining rapidly. We are trying to serve a solution to this enigma by providing a World Wide Web (WWW) interface to Biosoftware and at the same time have come up with a database integration tool that can query heterogeneous, geographically scattered and disparate databases simultaneously. In this report we will talk about BioInformatics in general with specific reference to BioInformatics Centre (BIC) at the National University of Singapore.

Computational Biology↗

Development of software tools at BioInformatics Centre (BIC) at the National University of Singapore (NUS).

There is burgeoning volume of information and data arising from the rapid research and unprecedented progress in molecular biology. This has been particularly affected by the Human Genome Project which is trying to completely sequence three billion nucleotides of the human genome (1),(1a). Other genome sequencing projects are also contributing substantially to this exponential growth in the number of DNA nucleotides and proteins sequenced. The number of journals, reports and research papers and tools required for the analysis of these sequences has also increased. For this the life sciences today needs tools in information technology and computation to prevent degeneration of this data into an inchoate accretion of unconnected facts and figures. The recently formed BioInformatics Centre (BIC) at the National University of Singapore (NUS) provides access to various commonly used computational tools available over the World Wide Web (WWW)--using a uniform interface and easy access. We have also come up with a new database tool. BioKleisli, which allows you to interact with various geographically scattered, heterogeneous, structurally complex and constantly evolving data sources. This paper summarises the importance of network access and database integration to biomedical research and gives a glimpse of current research conducted at BIC.

Animals↗

From sequence to structure to literature: the protocol approach to bioinformation.

Until the recent advent of high-throughput experimental data-acquisition in biology, the computational analysis of the biological data was predominantly on an ad hoc basis--i.e., the application of a given piece of software on the biological data depended on the need of the moment. This "functional approach" often resulted in piecemeal computational analysis with large amount of intervening "dead-time". The present high-throughput availability of experimental biological data requires a more streamlined and integrated "protocol approach". In this work, we illustrate such a user-friendly protocol using a common question frequently faced by a wet-lab bench-biologist--"Now that I have a DNA or protein sequence, what can I do with it using a computer?" As phrased, this question is steeped in the functional approach. In contrast, the protocol approach would re-phrase the same question as "Now that I have a DNA or protein sequence, what can a computer do for me?" Our integrating tool can start with a sequence and build a substantial custom data-warehouse of computationally derived sequence information, structure information and relevant published literature, that is continually updated.

Amino Acid Sequence↗

Protein folding: the endgame.

The last stage of protein folding, the "endgame," involves the ordering of amino acid side-chains into a well defined and closely packed configuration. We review a number of topics related to this process. We first describe how the observed packing in protein crystal structures is measured. Such measurements show that the protein interior is packed exceptionally tightly, more so than the protein surface or surrounding solvent and even more efficiently than crystals of simple organic molecules. In vitro protein folding experiments also show that the protein is close-packed in solution and that the tight packing and intercalation of side-chains is a final and essential step in the folding pathway. These experimental observations, in turn, suggest that a folded protein structure can be described as a kind of three-dimensional jigsaw puzzle and that predicting side-chain packing is possible in the sense of solving this puzzle. The major difficulty that must be overcome in predicting side-chain packing is a combinatorial "explosion" in the number of possible configurations. There has been much recent progress towards overcoming this problem, and we survey a variety of the approaches. These approaches differ principally in whether they use ab initio (physical) or more knowledge-based methods, how they divide up and search conformational space, and how they evaluate candidate configurations (using scoring functions). The accuracy of side-chain prediction depends crucially on the (assumed) positioning of the main-chain. Methods for predicting main-chain conformation are, in a sense, not as developed as that for side-chains. We conclude by surveying these methods. As with side-chain prediction, there are a great variety of approaches, which differ in how they divide up and search space and in how they score candidate conformations.

Animals↗

A structural explanation for the twilight zone of protein sequence homology.

Homology modeling of protein structures as a function of sequence breaks down at the twilight zone limit of sequence identity between the template and target proteins. Our results suggest that protein sequences that have diverged from a common ancestor beyond the twilight zone may adopt side-chain interactions that are very different from those endowed by the ancestral sequence.

Computer Simulation↗

Using a hydrophobic contact potential to evaluate native and near-native folds generated by molecular dynamics simulations.

There are several knowledge-based energy functions that can distinguish the native fold from a pool of grossly misfolded decoys for a given sequence of amino acids. These decoys, which are typically generated by mounting, or "threading", the sequence onto the backbones of unrelated protein structures, tend to be non-compact and quite different from the native structure: the root-mean-squared (RMS) deviations from the native are commonly in the range of 15 to 20 angstroms. Effective energy functions should also demonstrate a similar recognition capability when presented with compact decoys that depart only slightly in conformation from the correct structure (i.e. those with RMS deviations of approximately 5 angstroms or less). Recently, we developed a simple yet powerful method for native fold recognition based on the tendency for native folds to form hydrophobic cores. Our energy measure, which we call the hydrophobic fitness score, is challenged to recognize the native fold from 2000 near-native structures generated for each of five small monomeric proteins. First, 1000 conformations for each protein were generated by molecular dynamics simulation at room temperature. The average RMS deviation of this set of 5000 was 1.5 angstroms. A total of 323 decoys had energies lower than native; however, none of these had RMS deviations greater than 2 angstroms. Another 1000 structures were generated for each at high temperature, in which a greater range of conformational space was explored (4.3 angstroms RMS deviation). Out of this set, only seven decoys were misrecognized. The hydrophobic fitness energy of a conformation is strongly dependent upon the RMS deviation. On average our potential yields energy values which are lowest for the population of structures generated at room temperature, intermediate for those produced at high temperature and highest for those constructed by threading methods. In general, the lowest energy decoy conformations have backbones very close to native structure. The possible utility of our method for screening backbone candidates for the purpose of modelling by side-chain packing optimization is discussed.

Protein Folding↗

How similar must a template protein be for homology modeling by side-chain packing methods?

Given the correct backbone coordinates of a globular protein, side-chain packing methods can be generally expected to predict the side-chain coordinates of the buried core residues accurately. In the context of a study in modeling a family of bacteriophage DNA-binding proteins, we observed that when the coordinates of the actual perfect backbone are not available, the side-chain packing methods are still of predictive value using homologous but imperfect backbones. This is the situation in practical homology modeling where a target protein sequence is modeled from template structures of known protein homologs. In order to assess the quality and degree of accuracy of such predictions and their dependence on the extent of homology, we have now extended these studies to a well characterized family of globin structures that span a much wider range of sequence-structure similarity. The collective results show a clear relationship that is independent of protein family between side-chain prediction accuracy and the level of similarity between the template and target proteins. We judge this similarity in terms of sequence identity and the backbone r.m.s. deviation of the template structure used for modeling and the actual target structure in cases where the target structures are available. In summary, as sequence identity drops from 100% to about 50%, or when the backbone r.m.s. deviation between template and target structures increases from 0 A to about 1 A, the overall average r.m.s. error for the buried-core residues rises from 1.2 A to 1.5 A while the chi 1 prediction accuracy drops from 85% to 70-75% and the chi 2 prediction accuracy drops from 80% to 60-65%. When the sequence identity drops below 50% or the backbone r.m.s. deviation rises above 1 A, all 3 measures of prediction accuracy decrease rapidly. When the sequence identity edges to the so-called twilight zone of sequence similarity at around 22%, or when the backbone r.m.s. deviation exceeds 2 A, the prediction accuracy approaches the values to be expected for random predictions, namely, 3.1 A for average r.m.s. error, 22% and 29% for accuracy of chi 1 and chi 2 prediction. These observations provide a practical evaluation of the side-chain packing methods and are of value to the homology-modeler. The extent and degree to which the backbone topology of a protein fold can constrain internal side-chain orientation gives insight into the plasticity of the sequence-structure relationship found in the architecture of proteins.

Amino Acid Sequence↗

Recognizing native folds by the arrangement of hydrophobic and polar residues.

Central to the ab initio protein folding problem is the development of an energy function for which the correct native structure has a lower energy than all other conformations. Existing potentials of mean force typically rely extensively on database-derived contact frequencies or knowledge of three-dimensional structural information in order to be successful in the problem of recognizing the native fold for a given sequence from a set of decoy backbone conformations. Is the detailed statistical information or sophisticated analysis used by these knowledge-based potentials needed to achieve the observed degree of success in fold recognition? Here we introduce a novel pairwise energy function that enumerates contacts between hydrophobic residues while weighting their sum by the total number of residues surrounding these hydrophobic residues. Thus it effectively selects compact folds with the desired structural feature of a buried, intact core. This approach represents an advance over using pairwise terms whose energies of interaction that are independent of the position in the protein and greatly improves the discrimination capability of an energy function. Our results show that 85% of a set of 195 representative native folds were recognized correctly. The 29 exceptions were lipophilic proteins, small proteins with prosthetic groups or disulfide bonds, and oligomeric proteins. Overall, our method separates the native fold from incorrect folds by a larger margin (measured in standard deviation units) than has been previously demonstrated by more sophisticated methods. The arrangement of hydrophobic and polar residues alone as evaluated by our novel scoring scheme, is unexpectedly effective at recognizing native folds in general. It is surprising that a simple binary pattern of hydrophobic and polar residues apparently selects a give unique fold topology.

Amino Acids↗

The use of side-chain packing methods in modeling bacteriophage repressor and cro proteins.

In recent years, it has been repeatedly demonstrated that the coordinates of the main-chain atoms alone are sufficient to determine the side-chain conformations of buried residues of compact proteins. Given a perfect backbone, the side-chain packing method can predict the side-chain conformations to an accuracy as high as 1.2 A RMS deviation (RMSD) with greater than 80% of the chi angles correct. However, similarly rigorous studies have not been conducted to determine how well these apply, if at all, to the more important problem of homology modeling per se. Specifically, if the available backbone is imperfect, as expected for practical application of homology modeling, can packing constraints alone achieve sufficiently accurate predictions to be useful? Here, by systematically applying such methods to the pairwise modeling of two repressor and two cro proteins from the closely related bacteriophages 434 and P22, we find that when the backbone RMSD is 0.8 A, the prediction on buried side chain is accurate with an RMS error of 1.8 A and approximately 70% of the chi angles correctly predicted. When the backbone RMSD is larger, in the range of 1.6-1.8 A, the prediction quality is still significantly better than random, with RMS error at 2.2 A on the buried side chains and 60% accuracy on chi angles. Together these results suggest the following rules-of-thumb for homology modeling of buried side chains. When the sequence identity between the modeled sequence and the template sequence is > 50% (or, equivalently, the expected backbone RMSD is < 1 A), side-chain packing methods work well. When sequence identity is between 30-50%, reflecting a backbone RMS error of 1-2 A, it is still valid to use side-chain packing methods to predict the buried residues, albeit with care. When sequence identity is below 30% (or backbone RMS error greater than 2 A), the backbone constraint alone is unlikely to produce useful models. Other methods, such as those involving the use of database fragments to reconstruct a template backbone, may be necessary as a complementary guide for modeling.

Amino Acid Sequence↗

Low-resolution real-space envelopes: the application of the condensing-protocol approach to the ab initio macromolecular phase problem of a variety of examples.

A recently reported method - the improved condensing protocol [Subbiah (1991). Science, 252, 128-133; (1993). Acta Cryst. D49, 108-119] - for obtaining low-resolution macromolecular envelopes is applied to five varied test cases. These examples were chosen to illustrate the general applicability of the method to the wide range that typical macromolecular crystals adopt. The cases include small and large asymmetric unit volumes (4.7 x 10(4) to 1.17 x 10(6) A(3)), low and high symmetry (2 to 12 symmetry elements), small and large proteins (1570 to 12 216 non-H atoms), orthogonal and non-orthogonal unit cells, a wide variety of space groups (P2(1) to P6(3)22), small and large solvent contents (33-80%), and a case of non-crystallographic symmetry (threefold). In all five cases the inherent ambiguity of the condensing protocol in differentiating between bulk matter and bulk solvent is then resolved by use of the recently reported sign-fixing method (Subbiah, 1993).

Journal Article↗

Different protein sequences can give rise to highly similar folds through different stabilizing interactions.

We report an interesting case of structural similarity between 2 small, nonhomologous proteins, the third domain of ovomucoid (ovomucoid) and the C-terminal fragment of ribosomal L7/L12 protein (CTF). The region of similarity consists of a 3-stranded beta-sheet and an alpha-helix. This region is highly similar; the corresponding elements of secondary structure share a common topology, and the RMS difference for "equivalent" C alpha atoms is 1.6 A. Surprisingly, this common structure arises from completely different sequences. For the common core, the sequence identity is less than 3%, and there is neither significant sequence similarity nor similarity in the position or orientation of conserved hydrophobic residues. This superposition raises the question of how 2 entirely different sequences can produce an identical structure. Analyzing this common region in ovomucoid revealed that it is stabilized by disulfide bonds. In contrast, the corresponding structure in CTF is stabilized in the alpha-helix by a composition of residues with high helix-forming propensities. This result suggests that different sequences and different stabilizing interactions can produce an identical structure.

Amino Acid Sequence↗