Criteria that discriminate between native proteins and incorrectly folded models

Various theoretical concepts, such as free energy potentials, electrostatic interaction potentials, atomic packing, solvent-exposed surface, and surface charge distribution, were tested for their ability to discriminate between native proteins and misfolded protein models. Misfolded models were constructed by introducing incorrect side chains onto polypeptide backbones: side chains of the alpha-helical hemerythrin were modeled on the beta-sheeted backbone of immunoglobulin VL domain, whereas those of the VL domain were similarly modeled on the hemerythrin backbone. CONGEN, a conformational space sampling program, was used to construct the side chains, in contrast to the previous work, where incorrect side chains were modeled in all trans conformations. Capability of the conformational search procedure to reproduce native conformations was gauged first by rebuilding (the correct) side chains in hemerythrin and the VL domain: constructs with r.m.s. differences from the x-ray side chains 2.2-2.4 A were produced, and many calculated conformations matched the native ones quite well. Incorrectly folded models were then constructed by the same conformational protocol applied to incorrect amino acid sequences. All CONGEN constructs, both correctly and incorrectly folded, were characterized by exceptionally small molecular surfaces and low potential energies. Surface charge density, atomic packing, and Coulomb formula-based electrostatic interactions of the misfolded structures and the correctly folded proteins were similar, and therefore of little interest for diagnosing incorrect folds. The following criteria clearly favored the native structures over the misfolded ones: 1) solvent-exposed side-chain nonpolar surface, 2) number of buried ionizable groups, and 3) empirical free energy functions that incorporate solvent effects.     

A new method for mapping discontinuous antibody epitopes to reveal structural features of proteins.

Antibodies that bind to protein surfaces of interest can be used to report the three-dimensional structure of the protein as follows: Proteins are composed of linear polypeptide chains that fold together in complex spatial patterns to create the native protein structure. These folded structures form binding sites for antibodies. Antibody binding sites are typically "assembled" on the protein surface from segments that are far apart in the primary amino acid sequence of the target proteins. Short amino acid probe sequences that bind to the active region of each antibody can be used as witnesses to the antibody epitope surface and these probes can be efficiently selected from random sequence peptide libraries. This paper presents a new method to align these antibody epitopes to discontinuous regions of the one-dimensional amino acid sequence of a target protein. Such alignments of the epitopes indicate how segments of the protein sequence must be folded together in space and thus provide long-range constraints for solving the 3-D protein structure. This new antibody-based approach is applicable to the large fraction of proteins that are refractory to current approaches for structure determination and has the additional advantage of requiring very small amounts of the target protein. The binding site of an antibody is a surface, not just a continuous linear sequence, so the epitope mapping alignment problem is outside the scope of classical string alignment algorithms, such as Smith-Waterman. We formalize the alignment problem that is at the heart of this new approach, prove that the epitope mapping alignment problem is NP-complete, and give some initial results using a branch-and-bound algorithm to map two real-life cases. Initial results for two validation cases are presented for a graph-based protein surface neighbor mapping procedure that promises to provide additional spatial proximity information for the amino acid residues on the protein surface.     

Conformations of folded proteins in restricted spaces

A new method is presented to examine the complete range of folded topologies accessible in the compact state of globular proteins. The procedure is to generate all conformations, with volume exclusion, upon a lattice in a space restricted to the individual protein's known compact conformational space. Using one lattice point per residue, we find 10(2)-10(4) possible compact conformations for the five small globular proteins studied. Subsequently, these conformations are evaluated in terms of residue-specific, pairwise contact energies that favor nonbonded, hydrophobic interactions. Native structures for the five proteins are always found within the best 2% of all conformers generated. This novel method is simple and general and can be used to determine a small group of most favorable overall arrangements for the folding of specific amino acid sequences within a restricted space.     

Connection topology of proteins

One-dimensional amino acid sequences and three-dimensional foldedpolypeptide chains were modelled as non-directed graphs in whichnodes corresponded to amino acids and arcs repre sented connectionsbetween them. In the case of folded chains, non-backbone connectionswere assigned to amino acid pairs if their distance was lowerthan a threshold. Two topological indices, the connectednessnumber and the effective chain length were devised to comparefolding topologies. Loops created by non-backbone connectionsin the structure graphs were studied by simple graphical representations,revealing the hierarchy in native protein structures.     

Folding of polypeptide chains induced by the amino acid side-chains

Conformational calculations with the use of semi-empirical potential functions have been applied to the analysis of the folding of peptide chains. In particular, the part played by the amino acid side-chains in the adoption of folded conformations has been investigated.The results show that the preferred conformations of short peptides are mostly extended ones. However, from a given peptide chain-length, the side-chain to backbone and side-chain to side-chain interactions become strong enough so that definite sequences of amino acids can induce a transition from extended to folded conformations. We propose to call these folded structures “conformational nuclei”. The type of “nucleus” formed is dependent on both the amino acid composition and the sequence.Our results strongly support the hypothesis that folding of polypeptide chains can occur through a nucleation process that could be induced by the side-chains.     

Denatured proteins and early folding intermediates simulated in a reduced conformational space

Conformations of globular proteins in the denatured state were studied using a high-resolution lattice model of proteins and Monte Carlo dynamics. The model assumes a united-atom and high-coordination lattice representation of the polypeptide conformational space. The force field of the model mimics the short-range protein-like conformational stiffness, hydrophobic interactions of the side chains and the main-chain hydrogen bonds. Two types of approximations for the short-range interactions were compared: simple statistical potentials and knowledge-based protein-specific potentials derived from the sequence-structure compatibility of short fragments of protein chains. Model proteins in the denatured state are relatively compact, although the majority of the sampled conformations are globally different from the native fold. At the same time short protein fragments are mostly native-like. Thus, the denatured state of the model proteins has several features of the molten globule state observed experimentally. Statistical potentials induce native-like conformational propensities in the denatured state, especially for the fragments located in the core of folded proteins. Knowledge-based protein-specific potentials increase only slightly the level of similarity to the native conformations, in spite of their qualitatively higher specificity in the native structures. For a few cases, where fairly accurate experimental data exist, the simulation results are in semiquantitative agreement with the physical picture revealed by the experiments. This shows that the model studied in this work could be used efficiently in computational studies of protein dynamics in the denatured state, and consequently for studies of protein folding pathways, i.e. not only for the modeling of folded structures, as it was shown in previous studies. The results of the present studies also provide a new insight into the explanation of the Levinthal's paradox.     

Conformational stability of P22 tailspike proteins carrying temperature-sensitive folding mutations

The thermostable tailspike endorhamnosidase of Salmonella phage P22 provides a model system for comparing the role of amino acid sequences in determining the intracellular folding pathway with their role in stabilizing the mature structural protein. Complete Raman band assignments are given here for the native form of the tailspike trimer in aqueous solution. Once correctly folded and assembled, the wild-type and two well-characterized mutant proteins, tsfIle258----Leu and tsfGly323----Asp, exhibit the same secondary structure in solution, consisting predominantly of beta-strand (56 +/- 5%) and turns (17 +/- 2%). Raman bands that are sensitive indicators of hydrogen-bonding interactions of tyrosine (phenolic OH) and tryptophan (indole NH) are unchanged between 30 and 80 degrees C in both wild type and tsf mutants. Similarly, Raman bands that are sensitive to changes in the hydrophobic environment of nonpolar side chains exhibit no significant temperature dependence in wild type and tsf mutants. In contrast, these conformational features are greatly altered by chemical denaturation of the tailspike with lithium halide and guanidine hydrochloride. In the chemically denatured tailspike, the beta-strand structure is substantially converted to irregular or "random coil" conformation. These findings confirm conclusions from physiological studies that the three-dimensional structures of the tsf mutants, once stabilized at permissive temperatures, are equivalent to the native structure of the wild type, and this structure is maintained at temperatures far above those that block the folding of the chain into the final native conformation.(ABSTRACT TRUNCATED AT 250 WORDS)     

Computer simulations of protein folding with a small number of distance restraints

A high coordination lattice model was used to represent the protein chain. Lattice points correspond to amino-acid side groups. A complicated force field was designed in order to reproduce a protein-like behavior of the chain. Long-distance tertiary restraints were also introduced into the model. The Replica Exchange Monte Carlo method was applied to find the lowest energy states of the folded chain and to solve the problem of multiple minima. In this method, a set of replicas of the model chain was simulated independently in different temperatures with the exchanges of replicas allowed. The model chains, which consisted of up to 100 residues, were folded to structures whose root-mean-square deviation (RMSD) from their native state was between 2.5 and 5 A. Introduction of restrain based on the positions of the backbone hydrogen atoms led to an improvement in the number of successful simulation runs. A small improvement (about 0.5 A) was also achieved in the RMSD of the folds. The proposed method can be used for the refinement of structures determined experimentally from NMR data.     

The dominant role of side-chain backbone interactions in structural realization of amino acid code. ChiRotor: a side-chain prediction algorithm based on side-chain backbone interactions

The basic differences between the 20 natural amino acid residues are due to differences in their side-chain structures. This characteristic design of protein building blocks implies that side-chain-side-chain interactions play an important, even dominant role in 3D-structural realization of amino acid codes. Here we present the results of a comparative analysis of the contributions of side-chain-side-chain (s-s) and side-chain-backbone (s-b) interactions to the stabilization of folded protein structures within the framework of the CHARMm molecular data model. Contrary to intuition, our results suggest that side-chain-backbone interactions play the major role in side-chain packing, in stabilizing the folded structures, and in differentiating the folded structures from the unfolded or misfolded structures, while the interactions between side chains have a secondary effect. An additional analysis of electrostatic energies suggests that combinatorial dominance of the interactions between opposite charges makes the electrostatic interactions act as an unspecific folding force that stabilizes not only native structure, but also compact random conformations. This observation is in agreement with experimental findings that, in the denatured state, the charge-charge interactions stabilize more compact conformations. Taking advantage of the dominant role of side-chain-backbone interactions in side-chain packing to reduce the combinatorial problem, we developed a new algorithm, ChiRotor, for rapid prediction of side-chain conformations. We present the results of a validation study of the method based on a set of high resolution X-ray structures.     

Protein refolding in silico with atom-based statistical potentials and conformational search using a simple genetic algorithm

A distance-dependent atom-pair potential that treats long range and local interactions separately has been developed and optimized to distinguish native protein structures from sets of incorrect or decoy structures. Atoms are divided into 30 types based on chemical properties and relative position in the amino acid side-chains. Several parameters affecting the calculation and evaluation of this statistical potential, such as the reference state, the bin width, cutoff distances between pairs, and the number of residues separating the atom pairs, are adjusted to achieve the best discrimination. The native structure has the lowest energy for 39 of the 40 sets of original ROSETTA decoys (1000 structures per set) and 23 of the 25 improved decoys (approximately 1900 structures per set). Combined with the orientation-dependent backbone hydrogen bonding potential used by ROSETTA and a statistical solvation potential based on the solvent exclusion model of Lazaridis & Karplus, this potential is used as a scoring function for conformational search based on a genetic algorithm method. After unfolding the native structure by changing every phi and psi angle by either +/-3, +/-5 or +/-7 degrees, five small proteins can be efficiently refolded, in some cases to within 0.5 A C(alpha) distance matrix error (DME) to the native state. Although no significant correlation is found between the total energy and structural similarity to the native state, a surprisingly strong correlation exists between the radius of gyration and the DME for low energy structures.     

Unique oligomeric intermediates of bovine liver catalase

Catalases, although synthesized from single genes and built up from only one type of subunit, exist in heterogeneous form with respect to their conformations and association states in biological systems. This heterogeneity is not of genetic origin, but rather reflects the instability of this oligomeric heme enzyme. To understand better the factors that stabilize the various association states of catalase, we performed studies on the multimeric intermediates that are stabilized during guanidine-hydrochloride- and urea-induced unfolding of bovine liver catalase (BLC). For the first time, we have observed an enzymatically active, folded dimer of native BLC. This dimer has slightly higher enzymatic activity and altered structural properties compared to the native tetramer. Comparative studies of the effect of NaCl, GdmCl, and urea on BLC show that cation binding to negatively charged groups present in amino acid side chains of the enzyme leads to stabilization of an enzymatically active, folded dimer of BLC. Besides the folded dimer, an enzymatically active expanded tetramer and a partially unfolded, enzymatically inactive dimer of BLC were also observed. A complete recovery of native enzyme was observed on refolding of expanded tetramers and folded dimers; however, a very low recovery (maximum of approximately 5%) of native enzyme was observed on refolding of partially unfolded dimers and fully unfolded monomers.     

Surface amino acids as sites of temperature-sensitive folding mutations in the P22 tailspike protein

Temperature-sensitive folding (tsf) mutations in the gene for the thermostable P22 tailspike interfere with the polypeptide chain folding and association pathway at restrictive temperature without altering the thermostability of the protein once correctly folded and assembled at permissive temperature. Though the native proteins matured at permissive temperature are biologically active, many of them display alterations in electrophoretic mobility. The native forms of 15 of these tsf mutant proteins have been purified and characterized. The purified proteins differed in electrophoretic mobility and isoelectric point from wild type but did not show evidence of major conformational alterations. The results suggest that the electrophoretic variations conferred by the 15 tsf amino acid substitutions are due to changes in the net charge at solvent-accessible sites in the native form of the mutant protein. During the maturation of the chains at restrictive temperature, these sites influence the conformation of intermediates in chain folding and association. The amino acid sequences at these sites resemble those found at turns in polypeptide chains. The isolation of tsf mutations requires that the mature structure of the tailspike accommodates the mutant amino acid substitution without loss of function. The solvent-accessible sites are probably at the surface of this structural protein. This would explain how bulky mutant substitutions, such as arginines for glycines, are accommodated in the native tailspike structure. Such sites, stabilizing intermediates in the folding pathway and located on the surface of the mature protein, probably represent a general class of conformational substrates for tsf mutations.     

Automated search of natively folded protein fragments for high-throughput structure determination in structural genomics

Structural genomic projects envision almost routine protein structure determinations, which are currently imaginable only for small proteins with molecular weights below 25,000 Da. For larger proteins, structural insight can be obtained by breaking them into small segments of amino acid sequences that can fold into native structures, even when isolated from the rest of the protein. Such segments are autonomously folding units (AFU) and have sizes suitable for fast structural analyses. Here, we propose to expand an intuitive procedure often employed for identifying biologically important domains to an automatic method for detecting putative folded protein fragments. The procedure is based on the recognition that large proteins can be regarded as a combination of independent domains conserved among diverse organisms. We thus have developed a program that reorganizes the output of BLAST searches and detects regions with a large number of similar sequences. To automate the detection process, it is reduced to a simple geometrical problem of recognizing rectangular shaped elevations in a graph that plots the number of similar sequences at each residue of a query sequence. We used our program to quantitatively corroborate the premise that segments with conserved sequences correspond to domains that fold into native structures. We applied our program to a test data set composed of 99 amino acid sequences containing 150 segments with structures listed in the Protein Data Bank, and thus known to fold into native structures. Overall, the fragments identified by our program have an almost 50% probability of forming a native structure, and comparable results are observed with sequences containing domain linkers classified in SCOP. Furthermore, we verified that our program identifies AFU in libraries from various organisms, and we found a significant number of AFU candidates for structural analysis, covering an estimated 5 to 20% of the genomic databases. Altogether, these results argue that methods based on sequence similarity can be useful for dissecting large proteins into small autonomously folding domains, and such methods may provide an efficient support to structural genomics projects.     

Exploring critical determinants of protein amyloidogenesis: a review

The transformation of polypeptide chains from their globular native structure to fibrillar aggregates has been a matter of great concern because of the involvement of these aggregates in the onset of several degenerative diseases. These aggregates exhibit highly ordered cross β sheet structures and are known as ‘amyloids’. Formation of amyloids in the body is associated with cytotoxicity due to direct interaction of the aggregated species with the cell membrane leading to cellular permeability or due to loss of functionality of the proteins involved in the amyloid formation. The preference of polypeptide chains to remain in their native conformation or to aggregate into amyloids is guided by several factors such as its conformation at specific condition, concentration, physicochemical properties of the amino acid sequence and so on. In the current review, we have reviewed the different factors that guide the transition of proteins from their natively folded state to the amyloidogenic state. Understanding the critical determinants of amyloidogenesis is vital towards deciphering the molecular mechanism of amyloidogenesis and for the development of effective therapeutics against amyloidosis. Copyright © 2013 European Peptide Society and John Wiley & Sons, Ltd.     

Studies on protein folding,unfolding, and fluctuations by computer simulation. II. A. Three-dimensional lattice model of lysozyme

A three-dimensional lattice model of protein designed to assimilate lysozyme is introduced. An attractive interaction is assumed to work between preassigned specific pairs of units, when they occupy the nearest-nighbor lattice points. The behavior of this lattice lysozyme is studied by a Monte Carlo simulation method. Because of the specific interunit interactions,“native state” of the lattice lysozyme is stable at low temperatures. Conformational fluctuations in the native state are observed to occur at both termini and loop regions of the main chain existing on the surface. The process of unfolding and denatured states of this model are discussed. Complete refolding from a denatured state was not observed. However, by starting from partially folded structures, the native conformation could be attained. From these observation it is concluded that, in the process of folding of proteins as simplified in a lattice model, nulceation is a rate-limiting factor. The artificial character of this model and possible improvement are discussed.     

Assembly of protein tertiary structures from secondary structures using optimized potentials

We present a simulated annealing-based method for the prediction of the tertiary structures of proteins given knowledge of the secondary structure associated with each amino acid in the sequence. The backbone is represented in a detailed fashion whereas the sidechains and pairwise interactions are modeled in a simplified way, following the LINUS model of Srinivasan and Rose. A perceptron-based technique is used to optimize the interaction potentials for a training set of three proteins. For these proteins, the procedure is able to reproduce the tertiary structures to below 3 A in root mean square deviation (rmsd) from the PDB targets. We present the results of tests on twelve other proteins. For half of these, the lowest energy decoy has a rmsd from the native state below 6 A and, in 9 out of 12 cases, we obtain decoys whose rmsd from the native states are also well below 5 A.     

Acid-induced partly folded conformation resembling a molten globule state of xylanase from an alkalothermophilic Bacillus sp.

Nonnative protein structures having a compact secondary, but not rigid tertiary structure, have been increasingly observed as intermediate states in protein folding. We have shown for the first time during acid-induced unfolding of xylanase (Xyl II) the presence of a partially structured intermediate form resembling a molten globule state. The conformation and stability of Xyl II at acidic pH was investigated by equilibrium unfolding methods. Using intrinsic fluorescence and CD spectroscopic studies, we have established that Xyl II at pH 1.8 (A-state) retains the helical secondary structure of the native protein at pH 7.0, while the tertiary interactions are much weaker. At variance, from the native species (N-state), Xyl II in the A-state binds 1-anilino-8-sulfonic acid (ANS) indicating a considerable exposure of aromatic side chains. Lower concentration of Gdn HCl are required to unfold the A-state. For denaturation by Gdn HCl, the midpoint of the cooperative unfolding transition measured by fluorescence for the N-state is 3.5 +/- 0.1 M, which is higher than the value (2.2 +/- 0.1 M) observed for the A-state at pH 1.8. This alternatively folded state exhibits certain characteristics of the molten globule but differs distinctly from it by its structural stability that is characteristic for native proteins.     

Detection of local structures in reduced unfolded bovine pancreatic trypsin inhibitor.

The structure of BPTI and reduced BPTI in concentrated guanidinium HCl (GUHCl) in the presence of glycerol has been probed by measurements of dynamic nonradiative excitation energy transfer between probes attached to its amino groups. Interprobe distance distributions were obtained from analysis of donor fluorescence decay curves and used to characterize local structures in unordered states of the protein. Site specifically fluorescently labeled BPTI derivatives (1-n)BPTI (n = 15, 20, 41, 46) were used, each carrying a 2-methoxy-naphthyl-1-methylenyl group (MNA) at the N-terminal amino group of arg1 and 7-(dimethylamino)-coumarin-4-yl-acetyl residue (DA-coum) at one of its epsilon-NH2 groups of the lysine side chains. Analysis of donor fluorescence decay kinetics gave the interprobe distance distributions in the native and denatured states. The N-terminal-segment, residues 1-15, is in an extended conformation (with an average interprobe distance of 34 +/- 2 A) in the native state. Upon unfolding by reduction with DTT or beta-mercapto ethanol in 6 M GUHCl/glycerol mixture, the conformation of this segment relaxed to a state characterized by a reduced average interprobe distance and a larger width of the distances distribution. The average distance between residues 1 and 26, i.e., between the N-terminus and the turn of the twisted beta sheet element (residues 18-35), increased upon unfolding. At -30 degrees C in the above solvent, the distribution between these two sites was probably composed of two conformational subpopulations. About 45 +/- 20% of the molecules were characterized by a short interprobe distance (like the native state) representing a compact conformation, and 55 +/- 20% of the molecules showed large interprobe distances representing an expanded (unfolded) conformation. Thus local structures seem to exist in reduced denatured BPTI even under denaturing conditions in 6 M GUHCl/glycerol mixtures. Some of those structures are unstable in guanidinium isothiocyanate (GUSCN). The method introduced here is suitable for probing local structures and very long range interactions in unfolded proteins and for search for folding initiation sites (FISs) and early folding intermediates.     

A genetic algorithm that seeks native states of peptides and proteins.

We describe a computer algorithm to predict native structures of proteins and peptides from their primary sequences, their known native radii of gyration, and their known disulfide bonding patterns, starting from random conformations. Proteins are represented as simplified real-space main chains with single-bead side chains. Nonlocal interactions are taken from structural database-derived statistical potentials, as in an earlier treatment. Local interactions are taken from simulations of (phi, psi) energy surfaces for each amino acid generated using the Biosym Discover program. Conformational searching is done by a genetic algorithm-based method. Reasonable structures are obtained for melittin (a 26-mer), avian pancreatic polypeptide inhibitor (a 36-mer), crambin (a 46-mer), apamin (an 18-mer), tachyplesin (a 17-mer), C-peptide of ribonuclease A (a 13-mer), and four different designed helical peptides. A hydrogen bond interaction was tested and found to be generally unnecessary for helical peptides, but it helps fold some sheet regions in these structures. For the few longer chains we tested, the method appears not to converge. In those cases, it appears to recover native-like secondary structures, but gets incorrect tertiary folds.     

Prediction of natively unfolded regions in protein chain

We have shown that the ability of a protein to be in globular or in natively unfolded state (under native conditions) may be determined (besides low overall hydrophobicity and a large net charge) by such a property as the average environment density, the average number of residues enclosed at the given distance. A statistical scale of the average number of residues enclosed at the given distance for 20 types of amino acid residues in globular state has been created on the basis of 6626 protein structures. Using this scale for separation of 80 globular and 90 natively unfolded proteins we fail only in 11% of proteins (compared with 17% of errors which are observed if to use hydrophobicity scale). The present scale may be used both for prediction of form (folded or unfolded) of the native state of protein and for prediction of natively unfolded regions in protein chains. The results of comparison of our method of predicting natively unfolded regions with the other known methods show that our method has the highest fraction of correctly predicted natively unfolded regions (that is 87% and 77% if to make averaging over residues and over proteins correspondingly).