PubMed HealthSearch

Biomedical subjects

R Sowdhamini

Publications and source records attributed to R Sowdhamini.

12 recordsLinked to original sources

A database of globular protein structural domains: clustering of representative family members into similar folds.

BACKGROUND: A database of globular domains, derived from a non-redundant set of proteins, is useful for the sequence analysis of aligned domains, for structural comparisons, for understanding domain stability and flexibility and for fold recognition procedures. Domains are defined by the program DIAL and classified structurally using the procedure SEA. RESULTS: The DIAL-derived domain database (DDBASE) consists of 436 protein chains involving 695 protein domains. Of these, 206 are alpha-class, 191 are beta-class and 294 alpha and beta class. The domains, 63% from multidomain proteins and 73% less than 150 residues in length, were clustered automatically using both single-link cluster analysis and hierarchical clustering to give a quantitative estimate of similarity in the domain-fold space. CONCLUSIONS: Highly populated and well described folds (doubly wound alpha/beta, singly wound alpha/beta barrels, globins alpha, large Greek-key beta and flavin-binding alpha/beta) are recognized at a SEA cut-off score of 0.55 in single-link clustering and at 0.65 in hierarchical clustering, although functionally related families are usually clearly distinguished at more stringent values.

Cluster Analysis

An automatic method involving cluster analysis of secondary structures for the identification of domains in proteins.

With a growing number of structures available in the Brookhaven Protein Data Bank, automatic methods for domain identification are required for the construction of databases. Domains are considered to be clusters of secondary structure elements. Thus, helices and strands are first clustered using intersecondary structural distances between C alpha positions, and dendrograms based on this distance measure are used to identify domains. Individual domains are recognized by a disjoint factor, which enables the automatic identification and classification into disjoint, interacting, and conjoint domains. Application to a database of 83 protein families and 18 unique structures shows that the approach provides an effective delineation of boundaries and identifies those proteins that can be considered as a single domain. A quantitative estimate of the interaction between domains has been proposed. The database of protein domains is a useful tool for understanding protein folding, for recognizing protein folds, and for understanding structure-activity relationships.

Algorithms

The recognition of protein structure and function from sequence: adding value to genome data.

The explosion of DNA sequence data from genome projects presents many challenges. For instance, we must extend our current knowledge of protein structure and function so that it can be applied to these new sequences. The derivation of rules for the relationships between sequence and structure allow us to recognize a common fold by the use of tertiary templates. New techniques enable us to begin to meet the challenge of rule-based modelling of distantly related proteins. This paper describes an integrated and knowledge-based approach to the prediction of protein structure and function which can maximize the value of sequence information.

Amino Acid Sequence

Molecular evolution and domain structure of plasminogen-related growth factors (HGF/SF and HGF1/MSP).

Plasminogen-related growth factors, a new family of polypeptide growth factors with the basic domain organization and mechanism of activation of the blood proteinase plasminogen, include hepatocyte growth factor/scatter factor (HGF/SF), a potent effector of the growth, movement, and differentiation of epithelia and endothelia, and hepatocyte growth factor-like/macrophage stimulating protein (HGF1/MSP), an effector of macrophage chemotaxis and phagocytosis. Phylogeny of the serine proteinase domains and analysis of intron-exon boundaries and kringle sequences indicate that HGF/SF, HGF1/MSP, plasminogen, and apolipoprotein (a) have evolved from a common ancestral gene that consisted of an N-terminal domain corresponding to plasminogen activation peptide (PAP), 3 copies of the kringle domain, and a serine proteinase domain. Models of the N domains of HGF/SF, HGF1/MSP, and plasminogen, characterized by the presence of 4 conserved Cys residues forming a loop in a loop, have been modeled based on disulfide-bond constraints. There is a distinct pattern of charged and hydrophobic residues in the helix-strand-helix motif proposed for the PAP domain of HGF/SF; these may be important for receptor interaction. Three-dimensional structures of the 4 kringle and the serine proteinase domains of HGF/SF were constructed by comparative modeling using the suite of programs COMPOSER and were energy minimized. Docking of a lysine analogue indicates a putative lysine-binding pocket within kringle 2 (and possibly another in kringle 4). The models suggest a mechanism for the formation of a noncovalent HGF/SF homodimer that may be responsible for the activation of the Met receptor. These data provide evidence for the divergent evolution and structural similarity of plasminogen, HGF/SF, and HGF1/MSP, and highlight a new strategy for growth factor evolution, namely the adaptation of a proteolytic enzyme to a role in receptor activation.

Amino Acid Sequence

Conformational characteristics of asparaginyl residues in proteins.

Backbone conformations at 1064 asparaginyl residues in 123 non-homologous, high-resolution X-ray structures of proteins were analysed. Asn adopts conformations in left-handed alpha-helical region and other partially allowed regions in the Ramachandran map more readily than any other non-glycyl residue. Asn conformational clusters in the (phi, psi) regions of left-handed alpha-helix right-handed alpha-helix and extended (beta) strands were investigated in detail for their occurrence in various secondary structures, especially in beta-turn regions. Preferences were observed for Asn conformations in different positions in various beta-turn types, including the first and fourth positions of the turn. Asparaginyl residues with extended conformations are found to occur frequently in irregular regions, although they are expected to occur predominantly in extended strands or in the third position of type II beta-turns. Asn conformations at the N-cap positions of helices strongly prefer extended conformation than alpha L, which seems to be characteristic of non-glycyl residues at that position. In the liners connecting two extended strands and those connecting an alpha-helix and an extended strand, Asn with alpha L or alpha R conformation is more favoured than Asn with the beta-conformation. Analysis of Asn-Asn doublets and Asn-X-Asn triplets permitted identification of conformational families in such sequences. Results of this investigation provide useful hints in modelling Asn-rich regions in proteins such as malaria parasite coat protein.

Asparagine

Knowledge-based protein modeling.

Knowledge, both from the three-dimensional structures of homologous proteins and from the general analysis of protein structure, is of value in modeling a protein of known sequence but unknown structure. While many models are still constructed at least in part by manual methods on graphics devices, automated procedures have come into greater use. These procedures include those that assemble fragments of structure from other known structures and those that derive coordinates for the model from the satisfaction of restraints placed on atomic positions.

Amino Acid Sequence

Termination of right handed helices in proteins by residues in left handed helical conformations.

An analysis of 636 helical segments, ranging in length from 4 to 32 residues, from 123 independent protein crystal structures reveals that helix termination by residues in left handed (alpha 1) helical conformations is a common occurrence. Gly and Asn residues are the most frequent alpha L helix terminators, with the former having a very high propensity to adopt such conformations. The alpha R-alpha R-alpha R-alpha L segment at the C termini of protein helices often possesses a 6--> 1 (pi-type) hydrogen bond between the CO of residue i and the NH of residue i + 5 with residue i + 4 occurring in the alpha L conformation. A stereochemical analysis of 216 examples shows that in 62 cases the 6-->1 hydrogen bond is absent. The present analysis provides a quantitative measure of the propensity of the 20 amino acids to adopt alpha L helix terminating conformations.

Hydrogen Bonding

Modelling multiple disulphide loop containing polypeptides by random conformation generation. The test cases of alpha-conotoxin GI and endothelin I.

A general procedure for arriving at 3-D models of disulphide-rich polypeptide systems based on the covalent cross-link constraints has been developed. The procedure, which has been coded as a computer program, RANMOD, assigns a large number of random, permitted backbone conformations to the polypeptide and identifies stereochemically acceptable structures as plausible models based on strainless disulphide bridge modelling. Disulphide bond modelling is performed using the procedure MODIP developed earlier, in connection with the choice of suitable sites where disulphide bonds could be engineered in proteins (Sowdhamini, R., Srinivasan, N., Shoichet, B., Santi, D.V., Ramakrishnan, C. and Balaram, P. (1989) Protein Engng, 3, 95-103). The method RANMOD has been tested on small disulphide loops and the structures compared against preferred backbone conformations derived from an analysis of putative disulphide subdatabase and model calculations. RANMOD has been applied to disulphide-rich peptides and found to give rise to several stereochemically acceptable structures. The results obtained on the modelling of two test cases, alpha-conotoxin GI and endothelin I, are presented. Available NMR data suggest that such small systems exhibit conformational heterogeneity in solution. Hence, this approach for obtaining several distinct models is particularly attractive for the study of conformational excursions.

Amino Acid Sequence

Orthogonal beta beta motifs in proteins.

A super-secondary structural motif comprising two orthogonally oriented beta-strands connected by short linking segments of less than or equal to 5 residues has been identified from a data set of 65 independent protein crystal structures. Of the 42 examples from 14 proteins, a vast majority have only a single residue as the linking element. Analysis of the conformational angles at the junction reveals that the recently described type VIII beta-turn occurs frequently at the connecting hinge, while the type II beta-turn is also fairly common.

Animals

Conformations of disulfide bridges in proteins.

The conformational characteristics of disulfide bridges in proteins have been analyzed using a dataset of 22 protein structures, available at a resolution of less than or equal to 2.0 A, containing a total of 72 disulfide crosslinks. The parameters used in the analysis include (phi, psi) values at Cys residues, bridge dihedral angles chi ss, chi i1, chi j1, chi i2, and chi j2, the distances C alpha i-C alpha j and C beta i-C beta j between the C alpha and C beta atoms of Cys(i) and Cys(j). Eight families of bridge conformations with three or more occurrences have been identified on the basis of these stereochemical parameters. The most populated family corresponds to the "left handed spiral" identified earlier by Richardson [1981) Adv. Protein Chem. 34, 167-330). Disulfide bridging across antiparallel extended strands is observed in alpha-lytic protease, crambin, and beta-trypsin and this structure is shown to be very similar to those obtained in small cystine peptides. Solvent accessible surface area calculations show that the overwhelming majority of disulfide bridges are inaccessible to solvent.

Cystine

Stereochemical modeling of disulfide bridges. Criteria for introduction into proteins by site-directed mutagenesis.

A computer modeling procedure for assessing the stereochemical suitability of pairs of residues in proteins as potential sites for introduction of cystine disulfide crosslinks has been developed. Residue pairs with C alpha-C alpha distances of less than or equal to 6.5 A and C beta-C beta distances of less than or equal to 4.5 A are chosen for geometrical fixation of S atoms using the program MODIP. The stereochemistry of the modeled disulfides is evaluated using limits for the structural parameters of the various torsion angles and S-S bond length in the disulfide bridge. The ability of the procedure to correctly model disulfides has been checked with examples of cystine peptides of known crystal structures and 103 disulfide bridges from 25 available protein crystal structures determined at less than or equal to 2 A resolution. An analysis of results on three proteins with engineered disulfides, T4 lysozyme, dihydrofolate reductase and subtilisin, is presented. Two positions for the introduction of 'stereochemically optimal' disulfides are identified in subtilisin.

Disulfides