Protein structural alignments and functional genomics

Structural genomics-the systematic solution of structures of the proteins of an organism-will increasingly often produce molecules of unknown function with no close relative of known function. Prediction of protein function from structure has thereby become a challenging problem of computational molecular biology. The strong conservation of active site conformations in homologous proteins suggests a method for identifying them. This depends on the relationship between size and goodness-of-fit of aligned substructures in homologous proteins. For all pairs of proteins studied, the root-mean-square deviation (RMSD) as a function of the number of residues aligned varies exponentially for large common substructures and linearly for small common substructures. The exponent of the dependence at large common substructures is well correlated with the RMSD of the core as originally calculated by Chothia and Lesk (EMBO J 1986;5:823-826), affording the possibility of reconciling different structural alignment procedures. In the region of small common substructures, reduced aligned subsets define active sites and can be used to suggest the locations of active sites in homologous proteins.     

Regular grid-like substructures in the midgut epithelial basement membrane of some Coleoptera

Summary The midgut epithelial basement membranes in 13 species of Coleoptera belonging to 11 families have been examined ultrastructurally and are described in the present work. Regular grid-like substructures are present in 6 species. One of the basement membranes possessing a regular structure is roughly characterized histochemically on the light microscopic level; it contains tyrosine-rich protein and PAS-positive carbohydrate but apparently no acid mucosubstances.The function of these peculiar substructures is entirely unknown. No correlation between the feeding biology of the beetles and the occurrence and morphology of the substructures has been found. It may, however, be possible to link the substructures with the systematic position of the animals. This is supported by the fact that exactly the same characteristic pattern of substructures has been found in the four investigated species belonging to Polyphaga-Haplogastra and not anywhere else.     

The Origin of Modern 5S rRNA: A Case of Relating Models of Structural History to Phylogenetic Data

Evolutionary models of molecular structures must incorporate molecular information at different levels of structural complexity and must be phrased within a phylogenetic perspective. In this regard, phylogenetic trees of substructures that are reconstructed from molecular features that contribute to order and thermodynamic stability show that a gradual model of evolution of 5S rRNA structure is more parsimonious than models that invoke large segmental duplications of the molecule. The search for trees of substructures that are most parsimonious, by their very nature, defines an objective strategy to select models of molecular change that best fit structural data. When combined with additional data, such as the age of protein domains that interact with RNA substructures, these trees can be used to falsify unlikely hypotheses.     

Supramolecular structure of tripeptidyl peptidase II from human erythrocytes as studied by electron microscopy, and its correlation to enzyme activity.

Tripeptidyl peptidase II is an extralysosomal serine peptidase of an unusually large size, i.e. Mr greater than 10(6) for the native enzyme and Mr 135000 for the subunit. The enzyme from human erythrocytes was studied by electron microscopy on samples negatively stained by ammonium molybdate. Two different structural representations of the purified enzyme were obtained, both with a length of about 50 nm, and consisting of repetitive substructures. Upon dialysis of the enzyme against a Tris/HCl buffer, the activity was gradually decreased. This decrease was shown to parallel the dissociation of the large enzyme structures into smaller ones, the smallest measuring 3 nm by 10 nm and apparently corresponding to the repetitive substructures. The results indicate that a large polymeric form of the enzyme is a prerequisite for full activity.     

A new type of semi-empirical molecular orbital method for large molecules.

A new type of molecular orbital method is proposed. It is applicable to large molecules containing large conjugated substructures. Only π-electrons in the conjugated part, but all-valence electrons in the non-conjugated part of a molecule, are taken into account explicitly. The Fock matrix elements are evaluated from the semi-empirical values employed in the existing all-valence-electron methods. The examples presented here suggest that the new type of MO method predicts electronic structures which are quite similar to those obtained by complete semi-empirical MO calculations. This new method may make it possible to reasonably well describe the electronic structure of, and interaction between, large molecules using considerably less computation time and core storage than the complete calculation analogs.     

A data-mining approach for multiple structural alignment of proteins

Comparing the 3D structures of proteins is an important but computationally hard problem in bioinformatics. In this paper, we propose studying the problem when much less information or assumptions are available. We model the structural alignment of proteins as a combinatorial problem. In the problem, each protein is simply a set of points in the 3D space, without sequence order information, and the objective is to discover all large enough alignments for any subset of the input. We propose a data-mining approach for this problem. We first perform geometric hashing of the structures such that points with similar locations in the 3D space are hashed into the same bin in the hash table. The novelty is that we consider each bin as a coincidence group and mine for frequent patterns, which is a well-studied technique in data mining. We observe that these frequent patterns are already potentially large alignments. Then a simple heuristic is used to extend the alignments if possible. We implemented the algorithm and tested it using real protein structures. The results were compared with existing tools. They showed that the algorithm is capable of finding conserved substructures that do not preserve sequence order, especially those existing in protein interfaces. The algorithm can also identify conserved substructures of functionally similar structures within a mixture with dissimilar ones. The running time of the program was smaller or comparable to that of the existing tools.     

Mechanical unfoldons as building blocks of maltose-binding protein

Identifying independently folding cores or substructures is important for understanding and assaying the structure, function and assembly of large proteins. Here, we suggest mechanical stability as a criterion to identify building blocks of the 366 amino acid maltose-binding protein (MBP). We find that MBP, when pulled at its termini, unfolds via three (meta-) stable unfolding intermediates. Consequently, the MBP structure consists of four structural blocks (unfoldons) that detach sequentially from the folded structure upon force application. We used cysteine cross-link mutations to characterize the four unfoldons structurally. We showed that many MBP constructs composed of those building blocks indeed form stably folded structures in solution. Mechanical unfoldons may provide a new tool for a systematic search for stable substructures of large proteins.     

Localization and Organization of the Translational Apparatus in Mitochondrial Membranes of Astasia longa

In isolated mitochondrial membranes of Astasia longa, large and small circular substructures were observed with diameters of 1.30 and 0.35 m and thickness of 0.08 and 0.03 m, respectively. Such substructures were isolated by membrane treatment with proteolytic enzymes (proteinase K, trypsin) or by lipid solubilization with Triton X-100. After the removal of surface protein layer, we uncovered circular polyribosomes with similar diameters as those of original substructures. Polyribosomes were identified on the basis of their morphology, positive staining with uranyl acetate, a capacity for chloramphenicol-sensitive incorporation of ¹⁴C-amino acids into polypeptides, and from their buoyant density as estimated by equilibrium centrifugation in the CsCl density gradient. The conclusion is that, in mitochondrial membranes of A. longa, the translational apparatus is organized similarly to that in the membranes of chloroplasts and cyanobacteria, e.g., large and small circular polyribosomes situated within the membrane ring-shaped substructures are the basics for the formation of the latter.     

Chemical substructures that enrich for biological activity

MOTIVATION: Certain chemical substructures are present in many drugs. This has led to the claim of 'privileged' substructures which are predisposed to bioactivity. Because bias in screening library construction could explain this phenomenon, the existence of privilege has been controversial. RESULTS: Using diverse phenotypic assays, we defined bioactivity for multiple compound libraries. Many substructures were associated with bioactivity even after accounting for substructure prevalence in the library, thus validating the privileged substructure concept. Determinations of privilege were confirmed in independent assays and libraries. Our analysis also revealed 'underprivileged' substructures and 'conditional privilege'-rules relating combinations of substructure to bioactivity. Most previously reported substructures have been flat aromatic ring systems. Although we validated such substructures, we also identified three-dimensional privileged substructures. Most privileged substructures display a wide variety of substituents suggesting an entropic mechanism of privilege. Compounds containing privileged substructures had a doubled rate of bioactivity, suggesting practical consequences for pharmaceutical discovery.     

Investigation of Phycobilisome Subunit Interaction Interfaces by Coupled Cross-linking and Mass Spectrometry

The phycobilisome (PBS) is an extremely large light-harvesting complex, common in cyanobacteria and red algae, composed of rods and core substructures. These substructures are assembled from chromophore-bearing phycocyanin and allophycocyanin subunits, nonpigmented linker proteins and in some cases additional subunits. To date, despite the determination of crystal structures of isolated PBS components, critical questions regarding the interaction and energy flow between rods and core are still unresolved. Additionally, the arrangement of minor PBS components located inside the core cylinders is unknown. Different models of the general architecture of the PBS have been proposed, based on low resolution images from electron microscopy or high resolution crystal structures of isolated components. This work presents a model of the assembly of the rods onto the core arrangement and for the positions of inner core components, based on cross-linking and mass spectrometry analysis of isolated, functional intact Thermosynechococcus vulcanus PBS, as well as functional cross-linked adducts. The experimental results were utilized to predict potential docking interactions of different protein pairs. Combining modeling and cross-linking results, we identify specific interactions within the PBS subcomponents that enable us to suggest possible functional interactions between the chromophores of the rods and the core and improve our understanding of the assembly, structure, and function of PBS.     

Fast detection of common geometric substructure in proteins.

We consider the problem of identifying common three-dimensional substructures between proteins. Our method is based on comparing the shape of the alpha-carbon backbone structures of the proteins in order to find three-dimensional (3D) rigid motions that bring portions of the geometric structures into correspondence. We propose a geometric representation of protein backbone chains that is compact yet allows for similarity measures that are robust against noise and outliers. This representation encodes the structure of the backbone as a sequence of unit vectors, defined by each adjacent pair of alpha-carbons. We then define a measure of the similarity of two protein structures based on the root mean squared (RMS) distance between corresponding orientation vectors of the two proteins. Our measure has several advantages over measures that are commonly used for comparing protein shapes, such as the minimum RMS distance between the 3D positions of corresponding atoms in two proteins. A key advantage is that this new measure behaves well for identifying common substructures, in contrast with position-based measures where the nonmatching portions of the structure dominate the measure. At the same time, it avoids the quadratic space and computational difficulties associated with methods based on distance matrices and contact maps. We show applications of our approach to detecting common contiguous substructures in pairs of proteins, as well as the more difficult problem of identifying common protein domains (i.e., larger substructures that are not necessarily contiguous along the protein chain).     

Duplication, divergence and formation of novel protein topologies

The rearrangement or permutation of protein substructures is an important mode of divergence. Recent work explored one possible underlying mechanism called permutation-by-duplication, which produces special forms of motif rearrangements called circular permutations. Permutation-by-duplication, involving gene duplication, fusion and truncation, can produce fully functional intermediate proteins and thus represents a feasible mechanism of protein evolution. In spite of this, circular permutations are relatively rare and we discuss possible reasons for their existence.     

A simple autoradiographic method for investigating long range chromosome substructure: size and number of DNA molecules in isolated nucleoids of Escherichia coli.

A method is described for gently dissociating large DNA-protein complexes and for visualizing and quantitating the substructures by autoradiography. Using this technique, it is shown that nucleoids isolated from exponentially growing Escherichia coli (mean generation time = 35 min) contain on average 2.8 genome equivalents of DNA and that this nucleoid can be dissociated by deproteinization into two substructures having on average 1.4 genome equivalents. This result is correlated with previous sedimentation studies on the unfolded nucleoid DNA to explain prior inconsistencies. Scanning electron microscopy studies demonstrate that the shape and size of the isolated nucleoid is consistent with the proposed subunit structure of the in vivo nucleoid.     

Automatic domain decomposition of proteins by a Gaussian Network Model

Proteins are often comprised of domains of apparently independent folding units. These domains can be defined in various ways, but one useful definition divides the protein into substructures that seem to move more or less independently. The same methods that allow fairly accurate calculation of motion can be used to help classify these substructures. We show how the Gaussian Network Model (GNM), commonly used for determining motion, can also be adapted to automatically classify domains in proteins. Parallels between this physical network model and graph theory implementation are apparent. The method is applied to a nonredundant set of 55 proteins, and the results are compared to the visual assignments by crystallographers. Apart from decomposing proteins into structural domains, the algorithm can generally be applied to any large macromolecular system to decompose it into motionally decoupled sub-systems.     

Structural characterization of lignin during Pinus taeda wood treatment with Ceriporiopsis subvermispora

Pinus taeda wood chips were biotreated with Ceriporiopsis subvermispora under solid-state fermentation for periods varying from 15 to 90 days. Milled wood lignins extracted from sound and biotreated wood samples were characterized by wet-chemical and spectroscopic techniques. Treatment of the lignins by derivatization followed by reductive cleavage (DFRC) made it possible to detect DFRC monomers and dimers that are diagnostic of the occurrence of arylglycerol-beta-O-aryl and beta-beta, beta-5, beta-1, and 4-O-5 units in the lignin structure. Quantification of these DFRC products indicated that beta-O-aryl cleavage was a significant route for lignin biodegradation but that beta-beta, beta-5, beta-1, and 4-O-5 linkages were more resistant to the biological attack. The amount of aromatic hydroxyls did not increase with the split of beta-O-4 linkages, suggesting that the beta-O-4 cleavage products remain as quinone-type structures as detected by UV and visible spectroscopy. Nuclear magnetic resonance techniques also indicated the formation of new substructures containing nonoxygenated, saturated aliphatic carbons (CH(2) and CH(3)) in the side chains of lignins extracted from biotreated wood samples.     

Tableau-based protein substructure search using quadratic programming

Background  

Searching for proteins that contain similar substructures is an important task in structural biology. The exact solution of most formulations of this problem, including a recently published method based on tableaux, is too slow for practical use in scanning a large database.     

TESTLoc: protein subcellular localization prediction from EST data

Background  

The eukaryotic cell has an intricate architecture with compartments and substructures dedicated to particular biological processes. Knowing the subcellular location of proteins not only indicates how bio-processes are organized in different cellular compartments, but also contributes to unravelling the function of individual proteins. Computational localization prediction is possible based on sequence information alone, and has been successfully applied to proteins from virtually all subcellular compartments and all domains of life. However, we realized that current prediction tools do not perform well on partial protein sequences such as those inferred from Expressed Sequence Tag (EST) data, limiting the exploitation of the large and taxonomically most comprehensive body of sequence information from eukaryotes.     

Functional interaction of nuclear transport-defective simian virus 40 large T antigen with chromatin and nuclear matrix.

We analyzed the subcellular distribution of nuclear transport-defective simian virus 40 Lys-128-mutant (cT-3 [R. E. Lanford and J. S. Butel, Cell 37:801-813, 1984] and d10 [D. Kalderon, W. D. Richardson, A. F. Markham, and A. E. Smith, Nature (London) 311:33-38, 1984]) large T antigens in various Lys-128-mutant-transformed rodent cells and in Lys-128-mutant d10-infected TC7 cells. Small but significant amounts of the mutant large T antigens were found in association with nuclear substructures, both in mutant-transformed and in mutant-infected cells. Experiments with TC7 cells made incompetent for cell division by 60Co irradiation supported the assumption that Lys-128-mutant large T antigen did not associate with nuclear components during mitosis but most likely was transported into the nucleus because the Lys-128 mutation was leaky for nuclear transport. Low-level simian virus 40 DNA replication and production of infectious mutant virus progeny in TC7 cells indicated that the association of Lys-128-mutant large T antigen with nuclear substructures is functional.     

Fractionation of the epithelial apical junctional complex: reassessment of protein distributions in different substructures

The epithelial apical junctional complex (AJC) is an important regulator of cell structure and function. The AJC is compartmentalized into substructures comprising the tight and adherens junctions, and other membrane complexes containing the membrane proteins nectin, junctional adhesion molecule, and crumbs. In addition, many peripheral membrane proteins localize to the AJC. Studies of isolated proteins indicate a complex map of potential binding partners in which there is extensive overlap in the interactions between proteins in different AJC substructures. As an alternative to a direct search for specific protein-protein interactions, we sought to separate membrane substructures of the AJC in iodixanol density gradients and define their protein constituents. Results show that the AJC can be fractured into membrane substructures that contain specific membrane and peripheral membrane proteins. The composition of each substructure reveals a more limited overlap in common proteins than predicted from the inventory of potential interactions; some of the overlapping proteins may be involved in stepwise recruitment and assembly of AJC substructures.     

Secondary structure of streptokinase in aqueous solution: a Fourier transform infrared spectroscopic study.

The secondary structure of streptokinase (Sk) in aqueous solution was quantitatively examined by using Fourier transform infrared (FT-IR) spectroscopy. Resolution enhancement techniques, including Fourier deconvolution and derivative spectroscopy, were combined with band curve-fitting procedures to quantitate the spectral information from the amide I bands. Nine component bands were found under the broad, nearly featureless amide I bands which reflect the presence of various substructures. The relative areas of these component bands indicate an amount of beta-sheet between 30 and 37% and an alpha-helix content of only 12-13% in Sk. Further conformational substructures are assigned to turns (25-26%) and to "random" structures (15-16%). Additionally, the correlation of a pronounced component band near 1640 cm-1 (10-16% fractional area) with the possible presence of 3(10)-helices is discussed.