Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Evolutionary coordination”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 433 records · Page 24Linked to original sources

SpOct, a gene encoding the major octamer-binding protein in sea urchin embryos: expression profile, evolutionary relationships, and DNA binding of expressed protein.

We have characterized a sea urchin gene, SpOct, that encodes a 78-kDa POU-domain protein related to mammalian Oct-1 and Oct-2. The SpOct protein binds octamer elements in the promoters of the alpha H2B (Bell et al., 1992, Dev. Biol. 150, 363-371) and CyIIIa actin genes, and it closely resembles the major octamer-binding activity obtained from sea urchin blastula nuclear lysates in the size of its DNase I footprint on a canonical octamer element and in its relative binding affinity (Kr) for the octamer element versus poly(dAT) (1.4 x 10(4)). Moreover, partial protein sequences obtained from affinity-purified octamer-binding protein match sequences present in SpOct. These data suggest that SpOct is closely related to, if not identical with, the major octamer-binding activity in blastula nuclear extracts. RNA gel blots reveal four forms of SpOct mRNA, ranging in size from 4 to 12 kb. They are regulated coordinately in the embryo: all are present in the unfertilized egg, increase 28-fold in amount by the 8-hr blastula stage, and decline 6-fold by the 12-hr blastula stage. The same four size classes of SpOct mRNAs are present in several adult tissues, although their relative amounts vary. The temporal profile of SpOct mRNA expression in embryos closely resembles that of the alpha histone H2B gene. Our previous work (Bell et al., 1992) showed that expression of the alpha H2B gene in blastula-stage embryos was entirely dependent on an octamer element. Together, these data strongly suggest that SpOct may be the key regulator of the alpha H2B gene.

Amino Acid Sequence↗

Apoptosis in yeast: a new model system with applications in cell biology and medicine.

Apoptosis is a highly coordinated cellular suicide program crucial for metazoan health and diseases. Although its increasing importance in cancer, neurodegenerative disorders and AIDS led to intense research and a better understanding of apoptosis, many details of its regulation or the apoptotic phenotypes are poorly understood. The complex regulatory network and the often contradictory results obtained with human cell lines made application of an easier model system desirable. Apoptosis in yeast promises to provide a better understanding of the genetics of apoptosis. During the past 2 years, scientists were successful in identifying new cell-death regulators of humans, plants and fungi using Saccharomyces cerevisiae. The finding of apoptotic phenotypes, even in protists, suggests that apoptosis developed in unicellular organisms long before the evolutionary separation between fungi, plants and metazoan animals occurred.

Aging↗

The morph server: a standardized system for analyzing and visualizing macromolecular motions in a database framework.

The number of solved structures of macromolecules that have the same fold and thus exhibit some degree of conformational variability is rapidly increasing. It is consequently advantageous to develop a standardized terminology for describing this variability and automated systems for processing protein structures in different conformations. We have developed such a system as a 'front-end' server to our database of macromolecular motions. Our system attempts to describe a protein motion as a rigid-body rotation of a small 'core' relative to a larger one, using a set of hinges. The motion is placed in a standardized coordinate system so that all statistics between any two motions are directly comparable. We find that while this model can accommodate most protein motions, it cannot accommodate all; the degree to which a motion can be accommodated provides an aid in classifying it. Furthermore, we perform an adiabatic mapping (a restrained interpolation) between every two conformations. This gives some indication of the extent of the energetic barriers that need to be surmounted in the motion, and as a by-product results in a 'morph movie'. We make these movies available over the Web to aid in visualization. Many instances of conformational variability occur between proteins with somewhat different sequences. We can accommodate these differences in a rough fashion, generating an 'evolutionary morph'. Users have already submitted hundreds of examples of protein motions to our server, producing a comprehensive set of statistics. So far the statistics show that the median submitted motion has a rotation of approximately 10 degrees and a maximum Calpha displacement of 17 A. Almost all involve at least one large torsion angle change of >140 degrees. The server is accessible at http://bioinfo.mbb.yale. edu/MolMovDB

Computer Graphics↗

Evolutionary conservation of enzymatic catalysis: quantitative comparison of the effects of mutation of aligned residues in Saccharomyces cerevisiae and Escherichia coli inorganic pyrophosphatases on enzymatic activity.

Soluble inorganic pyrophosphatase (PPase) is one of the better understood phosphoryl-transfer enzymes and is distinctive in having four divalent metal ions at the active site. Here we determine pH profiles for wild-type Saccharomyces cerevisiae PPase (Y-PPase) and for 14 of its active site variants and consider the effects of active site mutation on the pH-independent parameters and acid dissociation constants that characterize these profiles against the framework of the proposed structure of the activated complex. The results obtained (a) support the current mechanistic model in which a hydroxide ion, stabilized by binding to two metal ions at the active site and by an extended system of hydrogen bonds within the active site, is the nucleophile that attacks enzyme-bound inorganic pyrophosphate and (b) provide evidence that the acid group that is necessary for maximal activity is a water molecule coordinated to a third metal ion, as shown by the general rise in the pKa of this group that is a consequence of almost all of the mutations. We further compare the present results to those previously observed for the corresponding mutations in Escherichia coli PPase [E-PPase; Salminen et al. (1995) Biochemistry 34, 782-791]. Such comparison provides a measure of the extent to which different portions of the active site are conserved. We find that some corresponding mutations have different effects on catalytic function, demonstrating that even in the context of very similar active sites, interactions of the mutated site with less well conserved portions of the enzyme, in this case outside the active site, can lead to different outcomes. On the other hand, one region of the active site is highly conserved, suggesting that it may represent a common feature of phosphoryl-transfer enzymes or a vestige of a primitive ur-PPase active site.

Amino Acid Sequence↗

Approximate protein structural alignment in polynomial time.

Alignment of protein structures is a fundamental task in computational molecular biology. Good structural alignments can help detect distant evolutionary relationships that are hard or impossible to discern from protein sequences alone. Here, we study the structural alignment problem as a family of optimization problems and develop an approximate polynomial-time algorithm to solve them. For a commonly used scoring function, the algorithm runs in O(n(10)/epsilon(6)) time, for globular protein of length n, and it detects alignments that score within an additive error of epsilon from all optima. Thus, we prove that this task is computationally feasible, although the method that we introduce is too slow to be a useful everyday tool. We argue that such approximate solutions are, in fact, of greater interest than exact ones because of the noisy nature of experimentally determined protein coordinates. The measurement of similarity between a pair of protein structures used by our algorithm involves the Euclidean distance between the structures (appropriately rigidly transformed). We show that an alternative approach, which relies on internal distance matrices, must incorporate sophisticated geometric ingredients if it is to guarantee optimality and run in polynomial time. We use these observations to visualize the scoring function for several real instances of the problem. Our investigations yield insights on the computational complexity of protein alignment under various scoring functions. These insights can be used in the design of scoring functions for which the optimum can be approximated efficiently and perhaps in the development of efficient algorithms for the multiple structural alignment problem.

Algorithms↗

Drosophila mutants of the kynurenine pathway as a model for ageing studies.

A search for Drosophila mutants with phenotypes similar to human diseases might help to unravel evolutionary conserved genes implicated in polygenic human disorders. Among these are neurodegenerative diseases, characterized by a late onset disturbance of memory, structural brain impairments and altered content of the intermediates of the kynurenine pathway. The ratio between kynurenate (KYNA) and 3-hydroxykynurenine (3-HOK) in the brain is a critical determinant of neuronal viability. Therefore, the Drosophila mutants cinnabar (KYNA excess) and cardinal (3-HOK excess) allow an evaluation of the specific roles of these metabolites which present in physiologic concentrations and mimic systemic administration. Previously we have demonstrated that the mutant cardinal can serve as a model for dementia and can help to unravel the earliest manifestations of brain dysfunction. Here we show that a state of the brain control of locomotor coordination characterized by the parameters of sound production in males results from the neuroprotective and neurotoxic effects of KYNA and 3-HOK accumulated in young and aged Drosophila mutants. The high instability of 1) cycle form and number in pulses; 2) of pulse amplitude and 3) rhythm in the courtship song of aged cardinal males are similar to the alterations in mutants with defective central complex of the brain. The cardinal mutants demonstrate apoptosis in the brain after stress treatment. This might reflect the misbalance in the content of excitatory amino acids' and the glycine site agonists revealed by HPLC-determination. The mutant cinnabar proved to be normal in respect of the parameters studied.

Aging↗

Is apolipoprotein D a mammalian bilin-binding protein?

Human apolipoprotein D (APO-D) is a serum glycoprotein that has no sequence similarity with other apolipoproteins but rather belongs to the alpha 2-microglobulin superfamily whose other members transport small hydrophobic ligands in a wide variety of biological contexts. To investigate the ligand specificity of APO-D, we analyzed its relationship with the other members of this superfamily and constructed a detailed molecular model using the atomic coordinates of its most closely related homolog--insecticyanin from the tobacco hornworm, Manduca sexta. We studied the geometry of the binding pocket of APO-D and the topology of characteristic patches of both hydrophobic and polar side chains that also occur in crystal structures of insecticyanin and bilin-binding protein from the butterfly Pieris brassicae. From the data obtained we hypothesize that heme-related compounds may be more favorable ligands for APO-D than either cholesterol or cholesteryl ester. Preliminary experiments showed that purified human APO-D binds bilirubin in an approximately one-to-one molar ratio. These results suggest a new biological role for APO-D that is more congruent with its tissue distribution and evolutionary history.

Amino Acid Sequence↗

Two distinct classes of keratin genes and their evolutionary significance.

Bacterial plasmids containing cDNA sequences specific for keratins were constructed from mRNA of cultured human epidermal cells. Two separate classes of cloned cDNAs were identified by positive hybrid selection: one class removed from total human epidermal mRNA a fraction that was translated into 56 and 58 kilodalton (kd) keratins, and the other class selected mRNAs that translated into a mixture of 50 kd and 46 kd keratins. When probes specific for these two keratin classes were hybridized with human DNA digested with a restriction endonuclease that does not cleave within the probe, two distinct patterns of about ten fragments each were observed. Most of the hybridizing genomic fragments corresponded to complete cDNA sequences, and it is estimated that each of the two classes is encoded by about 10 genes. When the probes were hybridized with DNA from different species, all vertebrates were found to contain discrete sequences homologous to both human keratin probes. Within each vertebrate species, the two probes always hybridized with approximately equal intensities to nonoverlapping sets of genomic sequences, suggesting a coordinate evolution between the two subfamilies of keratin genes. This finding has important functional implications for keratin filament assembly.

Animals↗

Repurposing anti-phage defenses to differentially arrest the viral lifecycle reveals the regulatory logic of a parasitic satellite.

Mobile genetic elements frequently encode defense mechanisms to protect their bacterial hosts from viral attack. In Vibrio cholerae, these defensive elements include phage-inducible chromosomal island-like elements (PLEs), which are phage satellites that act as highly specialized parasites of the lytic phage ICP1. While PLE transcriptional activation upon ICP1 infection is known to be temporally regulated, the underlying regulatory logic and dependencies on the progression of the phage's developmental program required for activation remain unclear. In this study, we took a novel approach to define these dependencies by introducing independent anti-phage defense systems, BREX and DarTG, as molecular roadblocks to impede the ICP1 lifecycle. We discovered that, for both ICP1 and PLE, late-stage gene expression is fundamentally uncoupled from genome replication, representing a striking departure from the standard paradigm for double-stranded DNA phages. While BREX restricts ICP1 to an immediate-early transcriptional state that stalls PLE activation, DarTG allows the phage to execute its full transcriptional cascade despite the total block in DNA replication. This permissive environment provides the necessary cue(s) for complete PLE induction, revealing that robust PLE activation is profoundly dependent on the transcriptional progression of its helper phage.IMPORTANCEBacteria and their viruses (phages) are locked in perpetual evolutionary conflict. Some bacteria harbor phage satellites, specialized parasites that are activated to hijack the phage's components to spread all the while inhibiting viral production. While some satellites respond to a single viral trigger, the regulation of many satellites, including clinically relevant phage-inducible chromosomal island-like elements (PLEs) in Vibrio cholerae, remains poorly understood. Here, we used bacterial defense systems as molecular roadblocks to probe how PLE activation depends on its helper phage. We found that severe disruptions to viral transcription stall PLE activation. Unexpectedly, both the virus and the satellite can execute their full transcriptional programs even when DNA replication is completely blocked, challenging a fundamental paradigm in virology. These insights reveal a sophisticated level of phage-satellite coordination, illustrating how satellite activation is tightly linked to the transcriptional state of its helper phage, a dependency that ultimately drives the dissemination of mobile genetic elements.

Vibrio cholerae↗

Hormones and history: the evolution and development of primate female sexuality.

Sexual behavior is required for reproduction in internally fertilizing species but poses significant social and physical risks. Females in many nonprimate species have evolved physical and behavioral mechanisms restricting sexual behavior to when females are fertile. The same hormones producing female fertility also control these mechanisms, assuring that sex only occurs when reproduction is possible. In contrast to nonprimate mammals, hormones do not regulate the capacity to engage in sex in female anthropoid primates, uncoupling fertility and the physical capacity to mate. Instead, in primates, sexual motivation has become the primary coordinator between sexual behavior and fertility. This dependence upon psychological mechanisms to coordinate physiology with behavior is possibly unique to primates, including humans, and allows a variety of nonphysiological influences, particularly social context, to regulate sexual behavior. The independence between hormonal state and sexual behavior allows sex to be used for social purposes. This complex regulation of primate sexuality develops during adolescence, where female monkeys show both hormonally influenced sexual motivation and socially modulated sexual behavior. We present findings from rhesus monkeys illustrating how social context and hormonal state interact to modulate adolescent and adult sexuality. It is argued that this flexibility in sexual behavior, combined with a tight regulation of sexual motivational systems by reproductive hormones, allows sexual behavior to be used for nonreproductive purposes while still assuring its occurrence during periods of female fertility. The evolutionary pressures that produced such flexibility in sexual behavior remain puzzling, but may reflect the importance of sexuality to primate social attraction and cohesion.

Adolescent↗

CoDIAC: A comprehensive approach for interaction analysis reveals novel insights into SH2 domain function and regulation.

Protein domains are conserved structural and functional units that serve as building blocks of proteins. Through evolutionary expansion, domain families are represented by multiple members in diverse configurations with other domains, evolving new specificities for their interacting partners. Here, we develop a structure-based interface analysis to comprehensively map domain interfaces from experimental and predicted structures, including interfaces with macromolecules and intraprotein interfaces. We hypothesized that comprehensive contact mapping of domains could yield new insights into domain selectivity, conservation of domain-domain interfaces across proteins, and identify conserved post-translational modifications (PTMs), relative to interaction interfaces, allowing for the inference of specific effects due to PTMs or mutations. We applied this approach to the human SH2 domain family, a modular unit central to phosphotyrosine-mediated signaling, identifying a novel approach to understanding binding selectivity and evidence of coordinated regulation of SH2 domain binding interfaces by tyrosine and serine/threonine phosphorylation and acetylation. These findings suggest multiple signaling systems can regulate protein activity and SH2 domain interactions in a coordinated manner. We provide the extensive features of the human SH2 domain family and this modular approach as an open source Python package for COmprehensive Domain Interface Analysis of Contacts (CoDIAC).

SH2 domains↗

The heme-copper oxidases of Thermus thermophilus catalyze the reduction of nitric oxide: evolutionary implications.

We show that the heme-copper terminal oxidases of Thermus thermophilus (called ba(3) and caa(3)) are able to catalyze the reduction of nitric oxide (NO) to nitrous oxide (N(2)O) under reducing anaerobic conditions. The rate of NO consumption and N(2)O production were found to be linearly dependent on enzyme concentration, and activity was abolished by enzyme denaturation. Thus, contrary to the eukaryotic enzyme, both T. thermophilus oxidases display a NO reductase activity (3.0 +/- 0.7 mol NO/mol ba(3) x min and 32 +/- 8 mol NO/mol caa(3) x min at [NO] approximately 50 microM and 20 degrees C) that, though considerably lower than that of bona fide NO reductases (300-4,500 mol NO/mol enzyme x min), is definitely significant. We also show that for ba(3) oxidase, NO reduction is associated to oxidation of cytochrome b at a rate compatible with turnover, suggesting a mechanism consistent with the stoichiometry of the overall reaction. We propose that the NO reductase activity of T. thermophilus oxidases may depend on a peculiar Cu(B)(+) coordination, which may be revealed by the forthcoming three-dimensional structure. These findings support the hypothesis of a common phylogeny of aerobic respiration and bacterial denitrification, which was proposed on the basis of structural similarities between the Pseudomonas stutzeri NO reductase and the cbb(3) terminal oxidases. Our findings represent functional evidence in support of this hypothesis.

Copper↗

Recommended nomenclature for the vertebrate alcohol dehydrogenase gene family.

The alcohol dehydrogenase (ADH) gene family encodes enzymes that metabolize a wide variety of substrates, including ethanol, retinol, other aliphatic alcohols, hydroxysteroids, and lipid peroxidation products. Studies on 19 vertebrate animals have identified ADH orthologs across several species, and this has now led to questions of how best to name ADH proteins and genes. Seven distinct classes of vertebrate ADH encoded by non-orthologous genes have been defined based upon sequence homology as well as unique catalytic properties or gene expression patterns. Each class of vertebrate ADH shares <70% sequence identity with other classes of ADH in the same species. Classes may be further divided into multiple closely related isoenzymes sharing >80% sequence identity such as the case for class I ADH where humans have three class I ADH genes, horses have two, and mice have only one. Presented here is a nomenclature that uses the widely accepted vertebrate ADH class system as its basis. It follows the guidelines of human and mouse gene nomenclature committees, which recommend coordinating names across species boundaries and eliminating Roman numerals and Greek symbols. We recommend that enzyme subunits be referred to by the symbol "ADH" (alcohol dehydrogenase) followed by an Arabic number denoting the class; i.e. ADH1 for class I ADH. For genes we recommend the italicized root symbol "ADH" for human and "Adh" for mouse, followed by the appropriate Arabic number for the class; i.e. ADH1 or Adh1 for class I ADH genes. For organisms where multiple species-specific isoenzymes exist within a class, we recommend adding a capital letter after the Arabic number; i.e. ADH1A, ADH1B, and ADH1C for human alpha, beta, and gamma class I ADHs, respectively. This nomenclature will accommodate newly discovered members of the vertebrate ADH family, and will facilitate functional and evolutionary studies.

Alcohol Dehydrogenase↗

An Integrated Sequence-Structure Database incorporating matching mRNA sequence, amino acid sequence and protein three-dimensional structure data.

We have constructed a non-homologous database, termed the Integrated Sequence-Structure Database (ISSD) which comprises the coding sequences of genes, amino acid sequences of the corresponding proteins, their secondary structure and straight phi,psi angles assignments, and polypeptide backbone coordinates. Each protein entry in the database holds the alignment of nucleotide sequence, amino acid sequence and the PDB three-dimensional structure data. The nucleotide and amino acid sequences for each entry are selected on the basis of exact matches of the source organism and cell environment. The current version 1.0 of ISSD is available on the WWW at http://www.protein.bio.msu.su/issd/ and includes 107 non-homologous mammalian proteins, of which 80 are human proteins. The database has been used by us for the analysis of synonymous codon usage patterns in mRNA sequences showing their correlation with the three-dimensional structure features in the encoded proteins. Possible ISSD applications include optimisation of protein expression, improvement of the protein structure prediction accuracy, and analysis of evolutionary aspects of the nucleotide sequence-protein structure relationship.

Algorithms↗

Molecular modeling of family GH16 glycoside hydrolases: potential roles for xyloglucan transglucosylases/hydrolases in cell wall modification in the poaceae.

Family GH16 glycoside hydrolases can be assigned to five subgroups according to their substrate specificities, including xyloglucan transglucosylases/hydrolases (XTHs), (1,3)-beta-galactanases, (1,4)-beta-galactanases/kappa-carrageenases, "nonspecific" (1,3/1,3;1,4)-beta-D-glucan endohydrolases, and (1,3;1,4)-beta-D-glucan endohydrolases. A structured family GH16 glycoside hydrolase database has been constructed (http://www.ghdb.uni-stuttgart.de) and provides multiple sequence alignments with functionally annotated amino acid residues and phylogenetic trees. The database has been used for homology modeling of seven glycoside hydrolases from the GH16 family with various substrate specificities, based on structural coordinates for (1,3;1,4)-beta-D-glucan endohydrolases and a kappa-carrageenase. In combination with multiple sequence alignments, the models predict the three-dimensional (3D) dispositions of amino acid residues in the substrate-binding and catalytic sites of XTHs and (1,3/1,3;1,4)-beta-d-glucan endohydrolases; there is no structural information available in the databases for the latter group of enzymes. Models of the XTHs, compared with the recently determined structure of a Populus tremulos x tremuloides XTH, reveal similarities with the active sites of family GH11 (1,4)-beta-D-xylan endohydrolases. From a biological viewpoint, the classification, molecular modeling and a new 3D structure of the P. tremulos x tremuloides XTH establish structural and evolutionary connections between XTHs, (1,3;1,4)-beta-D-glucan endohydrolases and xylan endohydrolases. These findings raise the possibility that XTHs from higher plants could be active not only on cell wall xyloglucans, but also on (1,3;1,4)-beta-D-glucans and arabinoxylans, which are major components of walls in grasses. A role for XTHs in (1,3;1,4)-beta-D-glucan and arabinoxylan modification would be consistent with the apparent overrepresentation of XTH sequences in cereal expressed sequence tags databases.

Amino Acid Sequence↗

[Structural, dynamic and interaction studies of proteins: the example of immunology].

The three-dimensional structure of a protein or of one of its functional domains cannot be predicted simply from the amino-acid sequence. It has to be determined either in the crystal, generally by X-ray diffraction, or in solution by high resolution NMR, or it can be modelled by comparison with a highly homologous protein. It cannot be represented simply by an average set of atomic coordinates. The dynamics of the molecule is part of its definition. The relative role of the structure and of the local and segmental mobility will be illustrated for recognition processes, mainly in immunological reactions. A unique conformational domain, member of the immunoglobulin superfamily, recognises differently its partners for an antigen-antibody reaction or for a cell adhesion process involving immunological determinants such as CD2 or CD4. Similar structural domains can be observed for completely different amino acid sequences, with the same type of dynamical properties, e.g. in neocarzinostatin. Restrained mobility of CDR loops and water molecules participation at the interface can explain the rapid binding of the ligands in AG/AB reactions, but cellular adhesion involves more hydrophobic and rigid interactions. This difference can be related to evolutionary aspects. Recognition mechanisms in cellular immunology, e.g. that of an intracellular foreign viral protein by a T lymphocyte cell, involve rather different biochemical processes in the so-called "peptidic self" model including proteolysis and cell surface presentation of the resulting peptides by MHC Class I proteins. These complexes are then recognised by the T cell receptors. The binding site of the TCR alpha and beta chains are structurally related to the antigen binding site of immunoglobulins. The peptide is anchored at both C- and N-terminals and should possess a restrained mobility. Recognition is thus a similar structural and dynamical process as in AG/AB reactions but is more secure due to the larger number of epitopes for a single protein.

Antigen-Antibody Reactions↗

Gamma carbonic anhydrases in plant mitochondria.

Three genes from Arabidopsis thaliana with high sequence similarity to gamma carbonic anhydrase (gammaCA), a Zn containing enzyme from Methanosarcina thermophila (CAM), were identified and characterized. Evolutionary and structural analyses predict that these genes code for active forms of gammaCA. Phylogenetic analyses reveal that these Arabidopsis gene products cluster together with CAM and related sequences from alpha and gamma proteobacteria, organisms proposed as the mitochondrial endosymbiont ancestor. Indeed, in vitro and in vivo experiments indicate that these gene products are transported into the mitochondria as occurs with several mitochondrial protein genes transferred, during evolution, from the endosymbiotic bacteria to the host genome. Moreover, putative CAM orthologous genes are detected in other plants and green algae and were predicted to be imported to mitochondria. Structural modeling and sequence analysis performed in more than a hundred homologous sequences show a high conservation of functionally important active site residues. Thus, the three histidine residues involved in Zn coordination (His 81, 117 and 122), Arg 59, Asp 61, Gin 75, and Asp 76 of CAM are conserved and properly arranged in the active site cavity of the models. Two other functionally important residues (Glu 62 and Glu 84 of CAM) are lacking, but alternative amino acids that might serve to their roles are postulated. Accordingly, we propose that photosynthetic eukaryotic organisms (green algae and plants) contain gammaCAs and that these enzymes codified by nuclear genes are imported into mitochondria to accomplish their biological function.

Amino Acid Sequence↗

Spatial, sequence-order-independent structural comparison of alpha/beta proteins: evolutionary implications.

We present a unique sequence-order independent approach which allows examination of three dimensional structures, searching for spatially similar substructural motifs. If the amino acids composing the motifs are contiguous in the primary chain, that is, they follow each other in the sequence, a common ancestor and a divergent evolutionary process may be implied. On the other hand, if the three-dimensional substructural motif consists of amino acids whose positions in the sequences vary between the different proteins, a convergent evolution might have taken place. Starting from different, ancient sequences, mutations may have occurred that brought about formation and conservation of a truly structural motif. Such a motif might be particularly suitable for fulfilling a specific function. Clearly, in order to be able to carry out such a task one needs a technique which allows comparisons of protein structures absolutely independent of their amino acid sequence-order. Our novel, efficient, computer vision based technique treats atoms (residues) as unconnected points in space, using strictly the atomic (either all atoms or only the C alpha atoms) coordinates. The order of the residues is completely disregarded. Detection, cataloging and analysis of "real" three-dimensional, sequence-order independent motifs in the crystallographic database is expected to be an invaluable tool for protein folding. Here we demonstrate the power of the technique by applying it to alpha/beta proteins. Our studies indicate that for some of the proteins, the "classical" structural alignments (conserving the amino acid order) are the optimal ones. Nevertheless, for others, truly spatial (out of sequential-order) amino acid equivalencing results in a better geometrical match.

Biological Evolution↗