Fold Space Research Articles

The role of protein homochirality in shaping the energy landscape of folding

The homochirality, or isotacticity, of the natural amino acids facilitates the formation of regular secondary structures such as alpha-helices and beta-sheets. However, many examples exist in nature where novel polypeptide topologies use both l- and d-amino acids. In this study, we explore how stereochemistry of the polypeptide backbone influences basic properties such as compactness and the size of fold space by simulating both lattice and all-atom polypeptide chains. We formulate a rectangular lattice chain model in both two and three dimensions, where monomers are chiral, having the effect of restricting local conformation. Syndiotactic chains with alternating chirality of adjacent monomers have a very large ensemble of accessible conformations characterized predominantly by extended structures. Isotactic chains on the other hand, have far fewer possible conformations and a significant fraction of these are compact. Syndiotactic chains are often unable to access maximally compact states available to their isotactic counterparts of the same length. Similar features are observed in all-atom models of isotactic versus syndiotactic polyalanine. Our results suggest that protein isotacticity has evolved to increase the enthalpy of chain collapse by facilitating compact helical states and to reduce the entropic cost of folding by restricting the size of the unfolded ensemble of competing states.

Connecting protein structure with predictions of regulatory sites

A common task posed by microarray experiments is to infer the binding site preferences for a known transcription factor from a collection of genes that it regulates and to ascertain whether the factor acts alone or in a complex. The converse problem can also be posed: Given a collection of binding sites, can the regulatory factor or complex of factors be inferred? Both tasks are substantially facilitated by using relatively simple homology models for protein-DNA interactions, as well as the rapidly expanding protein structure database. For budding yeast, we are able to construct reliable structural models for 67 transcription factors and with them redetermine factor binding sites by using a Bayesian Gibbs sampling algorithm and an extensive protein localization data set. For 49 factors in common with a prior analysis of this data set (based largely on phylogenetic conservation), we find that half of the previously predicted binding motifs are in need of some revision. We also solve the inverse problem of ascertaining the factors from the binding sites by assigning a correct protein fold to 25 of the 49 cases from a previous study. Our approach is easily extended to other organisms, including higher eukaryotes. Our study highlights the utility of enlarging current structural genomics projects that exhaustively sample fold structure space to include all factors with significantly different DNA-binding specificities.

Protein folding in the HP model on grid lattices with diagonals

The protein folding problem, i.e., the computational prediction of the three-dimensional structure of a protein from its amino acid sequence, is one of the most important and challenging problems in computational biology. Since a complete simulation of the folding process of a protein is far too complex to handle, one tries to find an approximate solution by using a simplified, abstract model. One of the most popular models is the so-called HP model, where the hydrophobic interactions between the amino acids are considered to be the main force in the folding process, and furthermore the folding space is modeled by a two- or three-dimensional grid lattice. In this paper, we will present some approximation algorithms for the protein folding problem in the HP model on an extended grid lattice with plane diagonals. The choice of this kind of lattice removes one of the major drawbacks of the original HP model, namely the bipartiteness of the grid which severely restricts the set of possible foldings. Our algorithms achieve an approximation ratio of 26 15 ≈ 1.733 for the two-dimensional and of 8 5 = 1.6 for the three-dimensional lattice. This improves significantly over the best previously known approximation ratios for the protein folding problem in the HP model on any lattice.

A reference database for circular dichroism spectroscopy covering fold and secondary structure space

Circular Dichroism (CD) spectroscopy is a long-established technique for studying protein secondary structures in solution. Empirical analyses of CD data rely on the availability of reference datasets comprised of far-UV CD spectra of proteins whose crystal structures have been determined. This article reports on the creation of a new reference dataset which effectively covers both secondary structure and fold space, and uses the higher information content available in synchrotron radiation circular dichroism (SRCD) spectra to more accurately predict secondary structure than has been possible with existing reference datasets. It also examines the effects of wavelength range, structural redundancy and different means of categorizing secondary structures on the accuracy of the analyses. In addition, it describes a novel use of hierarchical cluster analyses to identify protein relatedness based on spectral properties alone. The databases are shown to be applicable in both conventional CD and SRCD spectroscopic analyses of proteins. Hence, by combining new bioinformatics and biophysical methods, a database has been produced that should have wide applicability as a tool for structural molecular biology.

Multipolar representation of protein structure.

BackgroundThat the structure determines the function of proteins is a central paradigm in biology. However, protein functions are more directly related to cooperative effects at the residue and multi-residue scales. As such, current representations based on atomic coordinates can be considered inadequate. Bridging the gap between atomic-level structure and overall protein-level functionality requires parameterizations of the protein structure (and other physicochemical properties) in a quasi-continuous range, from a simple collection of unrelated amino acids coordinates to the highly synergistic organization of the whole protein entity, from a microscopic view in which each atom is completely resolved to a "macroscopic" description such as the one encoded in the three-dimensional protein shape.ResultsHere we propose such a parameterization and study its relationship to the standard Euclidian description based on amino acid representative coordinates. The representation uses multipoles associated with residue Cα coordinates as shape descriptors. We demonstrate that the multipoles can be used for the quantitative description of the protein shape and for the comparison of protein structures at various levels of detail. Specifically, we construct a (dis)similarity measure in multipolar configuration space, and show how such a function can be used for the comparison of a pair of proteins. We then test the parameterization on a benchmark set of the protein kinase-like superfamily. We prove that, when the biologically relevant portions of the proteins are retained, it can robustly discriminate between the various families in the set in a way not possible through sequence or conventional structural representations alone. We then compare our representation with the Cartesian coordinate description and show that, as expected, the correlation with that representation increases as the level of detail, measured by the highest rank of multipoles used in the representation, approaches the dimensionality of the fold space.ConclusionThe results described here demonstrate how a granular description of the protein structure can be achieved using multipolar coefficients. The description has the additional advantage of being immediately generalizable for any residue-specific property therefore providing a unitary framework for the study and comparison of the spatial profile of various protein properties.

Repeat protein architectures predicted by a continuum representation of fold space

It is an open question whether nature has utilized all possible protein folds. For a simple protein architecture, the helical repeats, we report a method to address this question based on a mapping between the set of repetitive curves and a space of parameters specifying the curve. The exploration of the parameter space for a particular architecture enables a systematic exploration of the fold space for that protein architecture. In a planar subspace of the parameter space of helical repeats we have identified points corresponding to both naturally occurring folds and potential folds not observed so far.

Beyond Mfold: Recent advances in RNA bioinformatics

Computational analysis of RNA secondary structure is a classical field of biosequence analysis, which has recently gained momentum due to the manyfold regulatory functions of RNA that have become apparent. We present five recent computational approaches that address the problems of synoptic folding space analysis, pseudoknot prediction, structure alignment, comparative structure prediction, and miRNA target prediction. All these programs are in current use and are available via the Bielefeld Bioinformatics Server at http://bibiserv.techfak.uni-bielefeld.de.

Complete probabilistic analysis of RNA shapes

BackgroundSoon after the first algorithms for RNA folding became available, it was recognised that the prediction of only one energetically optimal structure is insufficient to achieve reliable results. An in-depth analysis of the folding space as a whole appeared necessary to deduce the structural properties of a given RNA molecule reliably. Folding space analysis comprises various methods such as suboptimal folding, computation of base pair probabilities, sampling procedures and abstract shape analysis. Common to many approaches is the idea of partitioning the folding space into classes of structures, for which certain properties can be derived.ResultsIn this paper we extend the approach of abstract shape analysis. We show how to compute the accumulated probabilities of all structures that share the same shape. While this implies a complete (non-heuristic) analysis of the folding space, the computational effort depends only on the size of the shape space, which is much smaller. This approach has been integrated into the tool RNAshapes, and we apply it to various RNAs.ConclusionAnalyses of conformational switches show the existence of two shapes with probabilities approximately vs. , whereas the analysis of a microRNA precursor reveals one shape with a probability near to 1.0. Furthermore, it is shown that a shape can outperform an energetically more favourable one by achieving a higher probability. From these results, and the fact that we use a complete and exact analysis of the folding space, we conclude that this approach opens up new and promising routes for investigating and understanding RNA secondary structure.

Post-transcriptional nucleotide modification and alternative folding of RNA

Alternative foldings are an inherent property of RNA and a ubiquitous problem in scientific investigations. To a living organism, alternative foldings can be a blessing or a problem, and so nature has found both, ways to harness this property and ways to avoid the drawbacks. A simple and effective method employed by nature to avoid unwanted folding is the modulation of conformation space through post-transcriptional base modification. Modified nucleotides occur in almost all classes of natural RNAs in great chemical diversity. There are about 100 different base modifications known, which may perform a plethora of functions. The presumably most ancient and simple nucleotide modifications, such as methylations and uridine isomerization, are able to perform structural tasks on the most basic level, namely by blocking or reinforcing single base-pairs or even single hydrogen bonds in RNA. In this paper, functional, genomic and structural evidence on cases of folding space alteration by post-transcriptional modifications in native RNA are reviewed.

Diffusive searches in high-dimensional spaces and apparent ‘two-state’ behaviour in proteinfolding

We extend a simple model for protein folding as a high-dimensional diffusive search. Bysolving a steady-state diffusion equation on a hypersphere centred on an absorbing‘native state’ we find the general property that the kinetics of such a search willalways be nearly single exponential. This explains the common observation of suchsimple ‘two-state’ folding kinetics in models that contain considerable intermediatestructure. It also suggests that the experimental signature of single-exponentialfolding kinetics does not imply a simple two-state structure to the folding space.

Assessment of a Rigorous Transitive Profile Based Search Method to Detect Remotely Similar Proteins

Profile-based sequence search procedures are commonly employed to detect remote relationships between proteins. We provide an assessment of a Cascade PSI-BLAST protocol that rigorously employs intermediate sequences in detecting remote relationships between proteins. In this approach we detect using PSI-BLAST, which involves multiple rounds of iteration, an initial set of homologues for a protein in a ‘first generation’ search by querying a database. We propagate a ‘second generation’ search in the database, involving multiple runs of PSI-BLAST using each of the homologues identified in the previous generation as queries to recognize homologues not detected earlier. This non-directed search process can be viewed as an iteration of iterations that is continued to detect further homologues until no new hits are detectable. We present an assessment of the coverage of this ‘cascaded’ intermediate sequence search on diverse folds and find that searches for up to three generations detect most known homologues of a query. Our assessments show that this approach appears to perform better than the traditional use of PSI-BLAST by detecting 15% more relationships within a family and 35% more relationships within a superfamily. We show that such searches can be performed on generalized sequence databases and non-trivial relationships between proteins can be detected effectively. Such a propagation of searches maximizes the chances of detecting distant homologies by effectively scanning protein “fold space”.

Structural Genomics of Minimal Organisms and Protein Fold Space

The initial aim of the Berkeley Structural Genomics Center is to obtain a near-complete structural complement of two minimal organisms, closely related pathogens Mycoplasma genitalium and M. pneumoniae. The former has fewer than 500 genes and the latter fewer than 700 genes. To achieve this goal, the current protein targets have been selected starting with those predicted to be most tractable and likely to yield new structural and functional information. During the past 3 years, the semi-automated structural genomics pipeline has been set up from cloning, expression, purification, and ultimately to structural determination. The results from the pipeline substantially increased the coverage of the protein fold space of M. pneumoniae and M. genitalium. Furthermore, about 1/2 of the structures of 'unique' protein sequences revealed new and novel folds, and over 2/3 of the structures of previously annotated 'hypothetical proteins' inferred their molecular functions.

Visualization of Protein Fold Space via Nonmetric Multidimensional Scaling

In this paper, a 3D map of protein fold space was produced using Dali structure alignment and nonmetric multidimensional scaling. The fold space comprises four radial clusters, which correspond to the four classes of SCOP. The overall structure of the protein fold space is largely determined by three factors: secondary structure composition, topology of beta sheet, and domain size.

Fragnostic: walking through protein structure space

The Fragnostic () web tool implements a novel and useful view of protein structure space. We mined a non-redundant subset of the PDB for common fragments shared between proteins inhabiting different SCOP folds. Subsequently, we formulated an inter-fold similarity measure based on fragment sharing. Fold space is described as a graph whose nodes are folds between which the edges are drawn depending on the extent of fragment sharing. In this fashion, Fragnostic helps discover meaningful relationships between proteins belonging to different folds, based on sharing similar fragments in the proteins comprising those folds. Distant fold similarity information is supplemented by annotations taken from Gene Ontology, SCOP and CATH. Overall, Fragnostic is a tool which helps discover structural and functional relationships between proteins which are distantly related or seemingly unrelated.

Ribosome crystallography: catalysis and evolution of peptide-bond formation, nascent chain elongation and its co-translational folding

A ribosome is a ribozyme polymerizing amino acids, exploiting positional- and substrate-mediated chemical catalysis. We showed that peptide-bond formation is facilitated by the ribosomal architectural frame, provided by a sizable symmetry-related region in and around the peptidyl transferase centre, suggesting that the ribosomal active site was evolved by gene fusion. Mobility in tunnel components is exploited for elongation arrest as well as for trafficking nascent proteins into the folding space bordered by the bacterial chaperone, namely the trigger factor.

SH3-like Fold Proteins are Structurally Conserved and Functionally Divergent

The folding space for all the protein sequences is limited. Therefore it was observed that many proteins, whose sequences are not related, have similar fold characteristics. The fold databases like SCOP and CATH have classified various protein folds. However, in-depth analysis of the functional features of these folds was not done. We analyzed about twenty unique SH3-like folded proteins in their structural environment and functional characteristics. From our analysis it is apparent that the SH3-like folds could carry out various functions by modulation of loops and the functional region is restricted to one side of a particular sheet helped by two or three loops. The functions vary from oligonucleotide-binding to peptide-binding and other ligand binding. Although certain degree of sequence similarity was observed among the SH3-fold proteins, the similarity was restricted to the beta-strand regions of the proteins.

Trigger Factor in Complex with the Ribosome forms a Molecular Cradle for Nascent Proteins

During protein biosynthesis, nascent polypeptide chains that emerge from the ribosomal exit tunnel encounter ribosome-associated chaperones, which assist their folding to the native state1,2. Here we present a 2.7 A crystal structure of Escherichia coli trigger factor, the best-characterized chaperone of this type, together with the structure of its ribosome-binding domain in complex with the Haloarcula marismortui large ribosomal subunit. Trigger factor adopts a unique conformation resembling a crouching dragon with separated domains forming the amino-terminal ribosome-binding ‘tail’, the peptidyl-prolyl isomerase ‘head’, the carboxy-terminal ‘arms’ and connecting regions building up the ‘back’. From its attachment point on the ribosome, trigger factor projects the extended domains over the exit of the ribosomal tunnel, creating a protected folding space where nascent polypeptides may be shielded from proteases and aggregation. This study sheds new light on our understanding of co-translational protein folding, and suggests an unexpected mechanism of action for ribosome-associated chaperones.

NMR solution structure of Thermotoga maritima protein TM1509 reveals a Zn-metalloprotease-like tertiary structure

The 150-residue protein TM1509 is encoded in gene YF09_THEMA of Thermotoga maritima. TM1509 has so far no functional annotation and belongs to protein family UPF0054 (PFAM accession number: PF02130) which contains at least 146 members. The NMR structure of TM1509 reveals an alpha+beta fold comprising a four stranded beta-sheet with topology A( upward arrow), B( upward arrow), D( upward arrow), C( downward arrow) as well as five alpha-helices I-V. The structures of most members of family PF02130 can be reliably constructed using the TM1509 NMR structure, demonstrating high leverage for exploration of fold space. A multiple sequence alignment of TM1509 with homologues of family UPF0054 shows that three polypeptide segments, as well as a putative zinc-binding consensus motif HGXLHLXGYDH located at the C-terminal end of alpha-helix IV, are highly conserved. The spatial arrangement of the three His residues of this UPF0054 consensus motif is similar to the arrangement found for the His residues in the HEXXHXXGXXH zinc-binding consensus motif of matrix metallo-proteases (MMPs). Moreover, the other conserved polypeptide segments form a large cavity which encloses the putative Zn-binding pocket and might confer specificity during catalysis. However, TM1509 and the other members of the UPF0054 family do not have the crucial Glu residue in position 2 of the MMP consensus motif. Intriguingly, the TM1509 structure indicates that the Asp in the UPF0054 consensus motif (Asp 111 in TM1509) may overtake the catalytic role of the Glu. This suggests that protein family UPF0054 might contain members of a hitherto uncharacterized class of metalloproteases.

Graph-representation of oxidative folding pathways

BackgroundThe process of oxidative folding combines the formation of native disulfide bond with conformational folding resulting in the native three-dimensional fold. Oxidative folding pathways can be described in terms of disulfide intermediate species (DIS) which can also be isolated and characterized. Each DIS corresponds to a family of folding states (conformations) that the given DIS can adopt in three dimensions.ResultsThe oxidative folding space can be represented as a network of DIS states interconnected by disulfide interchange reactions that can either create/abolish or rearrange disulfide bridges. We propose a simple 3D representation wherein the states having the same number of disulfide bridges are placed on separate planes. In this representation, the shuffling transitions are within the planes, and the redox edges connect adjacent planes. In a number of experimentally studied cases (bovine pancreatic trypsin inhibitor, insulin-like growth factor and epidermal growth factor), the observed intermediates appear as part of contiguous oxidative folding pathways.ConclusionsSuch networks can be used to visualize folding pathways in terms of the experimentally observed intermediates. A simple visualization template written for the Tulip package can be obtained from V.A.

Optimal bundling of transmembrane helices using sparse distance constraints.

We present a two-step approach to modeling the transmembrane spanning helical bundles of integral membrane proteins using only sparse distance constraints, such as those derived from chemical cross-linking, dipolar EPR and FRET experiments. In Step 1, using an algorithm, we developed, the conformational space of membrane protein folds matching a set of distance constraints is explored to provide initial structures for local conformational searches. In Step 2, these structures refined against a custom penalty function that incorporates both measures derived from statistical analysis of solved membrane protein structures and distance constraints obtained from experiments. We begin by describing the statistical analysis of the solved membrane protein structures from which the theoretical portion of the penalty function was derived. We then describe the penalty function, and, using a set of six test cases, demonstrate that it is capable of distinguishing helical bundles that are close to the native bundle from those that are far from the native bundle. Finally, using a set of only 27 distance constraints extracted from the literature, we show that our method successfully recovers the structure of dark-adapted rhodopsin to within 3.2 A of the crystal structure.