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 comparative analysis of parechovirus protein structures with other picornaviruses

Parechoviruses belong to the genus Parechovirus within the family Picornaviridae and are non-enveloped icosahedral viruses with a single-stranded RNA genome. Parechoviruses include human and animal pathogens classified into six species. Those that infect humans belong to the Parechovirus A species and can cause infections ranging from mild gastrointestinal or respiratory illness to severe neonatal sepsis. There are no approved antivirals available to treat parechovirus (nor any other picornavirus) infections. In this parechovirus review, we focus on the cleaved protein products resulting from the polyprotein processing after translation comparing and contrasting their known or predicted structures and functions to those of other picornaviruses. The review also includes our original analysis from sequence and structure prediction. This review highlights significant structural differences between parechoviral and other picornaviral proteins, suggesting that parechovirus drug development should specifically be directed to parechoviral targets.     

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.     

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.     

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.     

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.     

Tripeptide inhibitors of Yersinia protein-tyrosine phosphatase

The protein-tyrosine phosphatase (PTP) 'YopH' is a virulence factor of Yersinia pestis, the causative agent of plague. Potential use of Yersinia as a bioterrorism agent renders YopH inhibitors of therapeutic importance. Previously, we had examined the inhibitory potencies of a variety of phosphotyrosyl (pTyr) mimetics against the human PTP1B enzyme by displaying them in the EGFR-derived hexapeptide sequence, 'Ac-Asp-Ala-Asp-Glu-Xxx-Leu-amide', where Xxx=pTyr mimetic. The poor inhibitory potencies of certain of these pTyr mimetics were attributed to restricted orientation within the PTP1B catalytic pocket incurred by extensive peripheral interaction of the hexapeptide platform. Utilizing the smaller tripeptide platform, 'Fmoc-Glu-Xxx-Leu-amide' we demonstrate herein that several of the low affinity hexapeptide-expressed pTyr mimetics exhibit high PTP1B affinity within the context of the tripeptide platform. Of particular note, the mono-anionic 4-(carboxydifluoromethyl)Phe residue exhibits affinity equivalent to the di-anionic F(2)Pmp residue, which had previously been among the most potent PTP-binding motifs. Against YopH, it was found that all tripeptides having Glu residues with an unprotected side chain carboxyl were inactive. Alternatively, in their Glu-OBn ester forms, several of the tripeptides exhibited good YopH affinity with the mono-anionic peptide, Fmoc-Glu(OBn)-Xxx-Leu-amide, where Xxx=4-(carboxymethyloxy)Phe providing an IC(50) value of 2.8 microM. One concern with such inhibitors is that they may potentially function by non-specific mechanisms. Studies with representative inhibitors, while failing to provide evidence of a non-specific promiscuous mode of inhibition, did indicate that non-classical inhibition may be involved.     

Flexibility of beta-sheets: principal component analysis of database protein structures

Protein folds are built primarily from the packing together of two types of structures: alpha-helices and beta-sheets. Neither structure is rigid, and the flexibility of helices and sheets is often important in determining the final fold (e.g., coiled coils and beta-barrels). Recent work has quantified the flexibility of alpha-helices using a principal component analysis (PCA) of database helical structures (J. Mol. Bio. 2003, 327, pp. 229-237). Here, we extend the analysis to beta-sheet flexibility using PCA on a database of beta-sheet structures. For sheets of varying dimension and geometry, we find two dominant modes of flexibility: twist and bend. The distributions of amplitudes for these modes are found to be Gaussian and independent, suggesting that the PCA twist and bend modes can be identified as the soft elastic normal modes of sheets. We consider the scaling of mode eigenvalues with sheet size and find that parallel beta-sheets are more rigid than antiparallel sheets over the entire range studied. Finally, we discuss the application of our PCA results to modeling and design of beta-sheet proteins.     

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.     

Genome analysis: Assigning protein coding regions to three-dimensional structures

We describe the results of a procedure for maximizing the number of sequences that can be reliably linked to a protein of known three-dimensional structure. Unlike other methods, which try to increase sensitivity through the use of fold recognition software, we only use conventional sequence alignment tools, but apply them in a manner that significantly increases the number of relationships detected. We analyzed 11 genomes and found that, depending on the genome, between 23 and 32% of the ORFs had significant matches to proteins of known structure. In all cases, the aligned region consisted of either >100 residues or >50% of the smaller sequence. Slightly higher percentages could be attained if smaller motifs were also included. This is significantly higher than most previously reported methods, even those that have a fold-recognition component. We survey the biochemical and structural characteristics of the most frequently occurring proteins, and discuss the extent to which alignment methods can realistically assign function to gene products.     

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.     

Amphipathic alpha-helices in proteins: results from analysis of protein structures

Amphipathic alpha-helices play a crucial role in mediating the interaction of peptides and proteins with membranes. We have analyzed protein structures for the occurrence of 18-residue amphipathic helices. We find several of these alpha-helices having average hydrophobic moments and average hydrophobicities that would favor their interaction with membranes. We have analyzed the distribution of net charge, helix length, normalized frequency of occurrence, and propensities of the 20 amino acids in the delineated 18-residue helices. We have observed distinct differences in the frequencies of occurrence of polar and hydrophobic amino acids at positions 1-18 in amphipathic and nonamphipathic helices. There are also differences in propensities of the 20 amino acids to occur at positions 1-18 of amphipathic and nonamphipathic helices. Synthetic peptides corresponding to some of these surface-seeking helices do possess antibacterial and/or hemolytic activities. Knowledge of the distribution of charges in 18-residue surface-seeking amphipathic alpha-helices, as well as propensity of occurrence of amino acids at various positions, would be useful inputs in the de novo design of amphipathic peptides.     

StoneHinge: Hinge prediction by network analysis of individual protein structures

Hinge motions are important for molecular recognition, and knowledge of their location can guide the sampling of protein conformations for docking. Predicting domains and intervening hinges is also important for identifying structurally self‐determinate units and anticipating the influence of mutations on protein flexibility and stability. Here we present StoneHinge, a novel approach for predicting hinges between domains using input from two complementary analyses of noncovalent bond networks: StoneHingeP, which identifies domain‐hinge‐domain signatures in ProFlex constraint counting results, and StoneHingeD, which does the same for DomDecomp Gaussian network analyses. Predictions for the two methods are compared to hinges defined in the literature and by visual inspection of interpolated motions between conformations in a series of proteins. For StoneHingeP, all the predicted hinges agree with hinge sites reported in the literature or observed visually, although some predictions include extra residues. Furthermore, no hinges are predicted in six hinge‐free proteins. On the other hand, StoneHingeD tends to overpredict the number of hinges, while accurately pinpointing hinge locations. By determining the consensus of their results, StoneHinge improves the specificity, predicting 11 of 13 hinges found both visually and in the literature for nine different open protein structures, and making no false‐positive predictions. By comparison, a popular hinge detection method that requires knowledge of both the open and closed conformations finds 10 of the 13 known hinges, while predicting four additional, false hinges.     

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.