Conformational analysis of alternative protein structures

MOTIVATION: Alternative structural models determined experimentally are available for an increasing number of proteins. Structural and functional studies of these proteins need to take these models into consideration as they can present considerable structural differences. The characterization of the structural differences and similarities between these models is a fundamental task in structural biology requiring appropriate methods. RESULTS: We propose a method for characterizing sets of alternative structural models. Three types of analysis are performed: grouping according to structural similarity, visualization and detection of structural variation and comparison of subsets for identifying and locating distinct conformational states. The alpha carbon atoms are used in order to analyse the backbone conformations. Alternatively, side-chain atoms are used for detailed conformational analysis of specific sites. The method takes into account estimates of atom coordinate uncertainty. The invariant regions are used to generate optimal superpositions of these models. We present the results obtained for three proteins showing different degrees of conformational variability: relative motion of two structurally conserved subdomains, a disordered subdomain and flexibility in the functional site associated with ligand binding. The method has been applied in the analysis of the alternative models available in SCOP. Considerable structural variability can be observed for most proteins. AVAILABILITY: The results of the analysis of the SCOP alternative models, the estimates of coordinate uncertainty as well as the source code of the implementation are available in the STRuster web site: http://struster.bioinf.mpi-inf.mpg.de.     

Tripeptide analysis of collagen synthesized by silica-exposed fibroblast cultures

The effect of silica on collagen biosynthesis by confluent monolayers of WI-38 fibroblast cultures was examined by a more comprehensive method of analysis. The presence of the particulates had no direct effect on protein (collagen) synthesis, proline incorporation or prolyl hydroxylase activity; the latter is determined by the degree of hydroxylation. Silica, however, was highly toxic to the cells.     

Conformational analysis of protein structures derived from NMR data

A study is presented of the conformational characteristics of NMR-derived protein structures in the Protein Data Bank compared to X-ray structures. Both ensemble and energy-minimized average structures are analyzed. We have addressed the problem using the methods developed for crystal structures by examining the distribution of ?, Ψ, and χ angles as indicators of global conformational irregularity. All these features in NMR structures occur to varying degrees in multiple conformational states. Some measures of local geometry are very tightly constrained by the methods used to generate the structure, e.g., proline ? angles, α-helix ?, Ψ angles, ω angles, and C_α chirality. The more lightly restrained torsion angles do show increasead clustering as the number of overall experimental observations increases. ?, Ψ, and χ₁ angle conformational heterogeneity is strongly correlated with accessibility but shows additional differences which reflect the differing number of observations possible in NMR for the various side chains (e.g., many for Trp, few for Ser). In general, we find that the core is defined to a notional resolution of 2.0 to 2.3 Å. Of real interest is the behavior of surface residues and in particular the side chains where multiple rotameric states in different structures can vary from 10% to 88%. Later generation structures show a much tighter definition which correlates with increasing use of J-coupling information, stereospecific assignments, and heteronumclear techniques. A suite of programs is being developed to address the special needs of NMR-derived structures which will take into account the existence of increased mobility in solution. © 1993 Wiley-Liss, Inc.     

Fourier transform infrared spectroscopic analysis of protein secondary structures

Infrared spectroscopy is one of the oldest and well established experimental techniques for the analysis of secondary structure of polypeptides and proteins. It is convenient, non-destructive, requires less sample preparation, and can be used under a wide variety of conditions. This review introduces the recent developments in Fourier transform infrared （FTIR） spectroscopy technique and its applications to protein structural studies. The experimental skills, data analysis, and correlations between the FTIR spectroscopic bands and protein secondary structure components are discussed. The applications of FTIR to the second- ary structure analysis, conformational changes, structural dynamics and stability studies of proteins are also discussed.     

Network analysis of the protein chain tertiary structures of heterocomplexes

In this paper, the tertiary structures of protein chains of heterocomplexes were mapped to 2D networks; based on the mapping approach, statistical properties of these networks were systematically studied. Firstly, our experimental results confirmed that the networks derived from protein structures possess small-world properties. Secondly, an interesting relationship between network average degree and the network size was discovered, which was quantified as an empirical function enabling us to estimate the number of residue contacts of the protein chains accurately. Thirdly, by analyzing the average clustering coefficient for nodes having the same degree in the network, it was found that the architectures of the networks and protein structures analyzed are hierarchically organized. Finally, network motifs were detected in the networks which are believed to determine the family or superfamily the networks belong to. The study of protein structures with the new perspective might shed some light on understanding the underlying laws of evolution, function and structures of proteins, and therefore would be complementary to other currently existing methods.     

A comprehensive analysis of non-sequential alignments between all protein structures

Background  

The majority of relations between proteins can be represented as a conventional sequential alignment. Nevertheless, unusual non-sequential alignments with different connectivity of the aligned fragments in compared proteins have been reported by many researchers. It is interesting to understand those non-sequential alignments; are they unique, sporadic cases or they occur frequently; do they belong to a few specific folds or spread among many different folds, as a common feature of protein structure. We present here a comprehensive large-scale study of non-sequential alignments between available protein structures in Protein Data Bank.     

Statistical analysis of unstructured amino acid residues in protein structures

We have performed a statistical analysis of unstructured amino acid residues in protein structures available in the databank of protein structures. Data on the occurrence of disordered regions at the ends and in the middle part of protein chains have been obtained: in the regions near the ends (at distance less than 30 residues from the N- or C-terminus), there are 66% of unstructured residues (38% are near the N-terminus and 28% are near the C-terminus), although these terminal regions include only 23% of the amino acid residues. The frequencies of occurrence of unstructured residues have been calculated for each of 20 types in different positions in the protein chain. It has been shown that relative frequencies of occurrence of unstructured residues of 20 types at the termini of protein chains differ from the ones in the middle part of the protein chain; amino acid residues of the same type have different probabilities to be unstructured in the terminal regions and in the middle part of the protein chain. The obtained frequencies of occurrence of unstructured residues in the middle part of the protein chain have been used as a scale for predicting disordered regions from amino acid sequence using the method (FoldUnfold) previously developed by us. This scale of frequencies of occurrence of unstructured residues correlates with the contact scale (previously developed by us and used for the same purpose) at a level of 95%. Testing the new scale on a database of 427 unstructured proteins and 559 completely structured proteins has shown that this scale can be successfully used for the prediction of disordered regions in protein chains.     

Reactivity-based analysis of domain structures in native replication protein A

Replication protein A (RPA) is an essential heterotrimeric ssDNA binding protein that participates in DNA repair, replication, and recombination. Though X-ray and NMR experiments have been used to determine three-dimensional structure models of the protein's domain fragments, a complete RPA structural model has not been reported. To test whether the fragment structures faithfully represent the same portions in the native solution-state protein, we have examined the structure of RPA under biologically relevant conditions. We have probed the location of multiple amino acids within the native RPA three-dimensional structure using reactivity of these amino acids toward proteolytic and chemical modification reagents. In turn, we evaluated different structural models by comparing the observed native RPA reactivities with anticipated reactivities based on candidate structural models. Our results show that our reactivity analysis approach is capable of critically assessing structure models and can be a basis for selecting the most relevant from among alternate models of a protein structure. Using this analytical approach, we verified the relevance of RPA fragment models to the native protein structure. Our results further indicate several important features of native RPA's structure in solution, such as flexibility at specific locations in RPA, particularly in the C-terminal region of RPA70. Our findings are consistent with reported DNA-free structural models and support the role of conformational change in the ssDNA binding mechanism of RPA.     

A comparative analysis of 23 structures of the amyloidogenic protein transthyretin

Self-assembly of the human plasma protein transthyretin (TTR) into unbranched insoluble amyloid fibrils occurs as a result of point mutations that destabilize the molecule, leading to conformational changes. The tertiary structure of native soluble TTR and many of its disease-causing mutants have been determined. Several independent studies by X-ray crystallography have suggested structural differences between TTR variants which are claimed to be of significance for amyloid formation. As these changes are minor and not consistent between the studies, we have compared all TTR structures available at the protein data bank including three wild-types, three non-amyloidogenic mutants, seven amyloidogenic mutants and nine complexes. The reference for this study is a new 1.5 A resolution structure of human wild-type TTR refined to an R-factor/R-free of 18.6 %/21.6 %. The present findings are discussed in the light of the previous structural studies of TTR variants, and show the reported structural differences to be non-significant.     

Shotgun cross-linking analysis for studying quaternary and tertiary protein structures

We developed a new approach that employs a novel computer algorithm for the sensitive and high-throughput analysis of tertiary and quaternary interaction sites from chemically cross-linked proteins or multi-protein complexes. First, we directly analyze the digests of the chemically cross-linked proteins using only high-accuracy LC-MS/MS data. We analyze these data using a computer algorithm, we term X!Link, to find cross-links between two peptides. Our algorithm is rapid, taking only a few seconds to analyze approximately 5000 MS/MS spectra. We applied this algorithm to analyze cross-linked sites generated chemically using the amino specific reagent, BS3, in both cytochrome c and the mitochondrial division dynamin mutant, Dnm1G385D, which exists as a stable homodimer. From cytochrome c, a well-established test protein, we identified a total of 31 cross-links, 21 interpeptide and 10 intrapeptide cross-links, in 257 MS/MS spectra from a single LC-MS/MS data set. The high sensitivity of this technique is indicated by the fact that all 19 lysines in cytochrome c were detected as a cross-link product and 33% of all the Lys pairs within 20 A were also observed as a cross-link. Analysis of the cross-linked dimeric form of Dnm1G385D identified a total of 46 cross-links, 38 interpeptide and 8 intrapeptide cross-links, in 98 MS/MS spectra in a single LC-MS/MS data set. These results represent the most abundant cross-links identified in a single protein or protein dimer to date. Statistical analysis suggests a 1% false discovery rate after optimization of filtering parameters. Further analysis of the cross-links identified using our approach indicates that careful manual inspection is important for the correct assignment of cross-linking sites when multiple cross-linkable sites or several similar sequences exist. In summary, we have developed a sensitive MS-based approach to identify peptide-peptide cross-links that does not require isotopic labeling or comparison with non-cross-linked controls, making it faster and simpler than current methodologies.     

Functional structures of the recA protein found by chimera analysis.

We developed a novel genetic method for finding functional regions of a protein by the analysis of chimeras formed between homologous proteins. Sets of chimeric genes were made by intramolecular homologous recombination in a linearized plasmid DNA carrying both recA genes of Escherichia coli and Pseudomonas aeruginosa. A recBCsbcA strain of E. coli was used for isolation of plasmids carrying recombinants between these genes. Examination of properties of E. coli strains deleting the recA gene and carrying a plasmid with a chimeric gene shows that chimera formation at certain positions inactivates a RecA function. Frequently, all chimeras with a junction in a certain region of the protein inactivate a function. Rather than a direct effect of the presence of the junction at a particular position, mismatching of the regions both sides of the junction that are derived from the different species is responsible for the inactivation. For a chimeric protein to be functional, certain pairs of sequences in different regions of the protein must derive from the same parent. Four pairs of such sequences were found: two are involved in activities for genetic recombination and for resistance to ultraviolet light irradiation and the others in formation of active oligomers. Regions defined by these sequences are located in the looped regions of the protein. A pair of regions may co-operate to form a functional folded structure.     

Topological analysis and interactive visualization of biological networks and protein structures

Computational analysis and interactive visualization of biological networks and protein structures are common tasks for gaining insight into biological processes. This protocol describes three workflows based on the NetworkAnalyzer and RINalyzer plug-ins for Cytoscape, a popular software platform for networks. NetworkAnalyzer has become a standard Cytoscape tool for comprehensive network topology analysis. In addition, RINalyzer provides methods for exploring residue interaction networks derived from protein structures. The first workflow uses NetworkAnalyzer to perform a topological analysis of biological networks. The second workflow applies RINalyzer to study protein structure and function and to compute network centrality measures. The third workflow combines NetworkAnalyzer and RINalyzer to compare residue networks. The full protocol can be completed in ～2 h.     

Tripeptide hydroxamate from Corynebacterium kutscheri.

A tripeptide hydroxamate was isolated from cultures of Corynebacterium kutscheri grown on an iron-limiting medium. The metabolite was characterized by spectral measurements and by chemical degradation as L-alpha-aspartyl-L-alpha-N-hydroxy-aspartyl-D-cycloserine.     

Tripeptide analogues of MG132 as protease inhibitors

The 26S proteasome and calpain are linked to a number of important human diseases. Here, we report a series of analogues of the prototypical tripeptide aldehyde inhibitor MG132 that show a unique combination of high activity and selectivity for calpains over proteasome. Tripeptide aldehydes (1–3) with an aromatic P3 substituent show enhanced activity and selectivity against ovine calpain 2 relative to chymotrypsin-like activity of proteasome. Docking studies reveal the key contacts between inhibitors and calpain to confirm the importance of the S3 pocket with respect to selectivity between calpains 1 and 2 and the proteasome.     

Flexibility of alpha-helices: results of a statistical analysis of database protein structures

Alpha-helices stand out as common and relatively invariant secondary structural elements of proteins. However, alpha-helices are not rigid bodies and their deformations can be significant in protein function (e.g. coiled coils). To quantify the flexibility of alpha-helices we have performed a structural principal-component analysis of helices of different lengths from a representative set of protein folds in the Protein Data Bank. We find three dominant modes of flexibility: two degenerate bend modes and one twist mode. The data are consistent with independent Gaussian distributions for each mode. The mode eigenvalues, which measure flexibility, follow simple scaling forms as a function of helix length. The dominant bend and twist modes and their harmonics are reproduced by a simple spring model, which incorporates hydrogen-bonding and excluded volume. As an application, we examine the amount of bend and twist in helices making up all coiled-coil proteins in SCOP. Incorporation of alpha-helix flexibility into structure refinement and design is discussed.     

Secondary structures without backbone: an analysis of backbone mimicry by polar side chains in protein structures.

Backbone mimicry by the formation of closed-loop C7, C10 and C13 (mimics of gamma-, beta- and alpha-turns) conformations through side chain-main chain hydrogen bonds by polar groups is a frequent observation in protein structures. A data set of 250 non-homologous and high-resolution protein crystal structures was used to analyze these conformations for their characteristic features. Seven out of the nine polar residues (Ser, Thr, Asn, Asp, Gln, Glu and His) have hydrogen bonding groups in their side chains which can participate in such mimicry and as many as 15% of all these polar residues engage in such conformations. The distributions of dihedral angles of these mimics indicate that only certain combinations of the dihedral angles involved aid the formation of these mimics. The observed examples were categorized into various classes based on these combinations, resulting in well defined motifs. Asn and Asp residues show a very high capability to perform such backbone secondary structural mimicry. The most highly mimicked backbone structure is of the C10 conformation by the Asx residues. The mimics formed by His, Ser, Thr and Glx residues are also discussed. The role of such conformations in initiating the formation of regular secondary structures during the course of protein folding seems significant.     

Incorporation of crystallographic temperature factors in the statistical analysis of protein tertiary structures

A method to identify statistically significant differences between equivalent atoms in two closely related protein X-ray crystallographic structures is described. This method uses the linear relationship found between the logarithm of the distance between equivalent atoms and their mean temperature factor to determine, by linear regression, the expected difference and variance.     

Functional analysis of substrate and cofactor complex structures of a thymidylate synthase-complementing protein

Like thymidylate synthase (TS) in eukaryotes, the thymidylate synthase-complementing proteins (TSCPs) are mandatory for cell survival of many prokaryotes in the absence of external sources of thymidylate. Details of the mechanism of this novel family of enzymes are unknown. Here, we report the structural and functional analysis of a TSCP from Thermotoga maritima and its complexes with substrate, analogs, and cofactor. The structures presented here provide a basis for rationalizing the TSCP catalysis and reveal the possibility of the design of an inhibitor. We have identified a new helix-loop-strand FAD binding motif characteristic of the enzymes in the TSCP family. The presence of a hydrophobic core with residues conserved among the TSCP family suggests a common overall fold.