Solution structures of proteins from NMR data and modeling: alternative folds for neutrophil peptide 5

The structure of neutrophil peptide 5 in solution has recently been reported (Pardi et al., 1988). The structure determination was accomplished by using a distance geometry algorithm and 107 interproton distance constraints obtained from 2D NMR data. In each of the eight independent solutions to the distance geometry equations, the overall fold of the polypeptide backbone was identical and the root mean square (rms) deviation between backbone atoms of the superimposed structures was small (approximately 2.4 A). In this paper we report additional NP-5 structures obtained by using a new structure generation algorithm: a Monte Carlo search in torsion angle space. These structures have a large rms backbone deviation from the distance geometry structures (approximately 5.0 A). The backbone topologies differ in significant respects from the distance geometry structures and from each other. Structures are found that are pseudo mirror images of part or all of the fold corresponding to that first obtained with the distance geometry procedure. For small proteins, the problem of distinguishing the correct structure among pseudo mirror images is likely to be greater than previously recognized. When a set of test distance constraints constructed from a novel Monte Carlo structure is used as input in the distance geometry algorithm, the fold of the resulting structure does not correspond to that of the target. The results also demonstrate that the previously accepted criteria (the magnitude of the rms deviation between multiple solutions of the distance geometry equations) for defining the accuracy and precision of a peptide structure generated from NMR data are inadequate. An energetic analysis of structures corresponding to the different folding topologies has been carried out. The molecular mechanics energies obtained by minimization and molecular dynamics refinement provide sufficient information to eliminate certain alternative structures. On the basis of a careful comparison of the different trial structures with the experimental data, it is concluded that the NP-5 peptide fold which was originally reported is most consistent with the data. An alternative fold corresponding to structures with low energies and small total distance violations is ruled out because for this fold predicted NOEs are not observed experimentally.     

An efficient high-throughput resonance assignment procedure for structural genomics and protein folding research by NMR.

Sequence specific resonance assignment is the primary requirement for all investigations of proteins by NMR methods. In the present postgenomic era where structural genomics and protein folding have occupied the center stage of NMR research, there is a high demand on the speed of resonance assignment, whereas the presently available methods based either on NOESY or on some triple-resonance experiments are rather slow. They also have limited success with unfolded proteins because of the lack of NOEs, and poor dispersion of amide and carbon chemical shifts. This paper describes an efficient approach to rapid resonance assignment that is suitable for both folded and unfolded proteins, making use of the triple-resonance experiments described recently [HNN and HN(C)N]. It has three underlying principles. First, the experiments exploit the (15)N chemical shift dispersions which are generally very good for both folded and unfolded proteins, along two of the three dimensions; second, they directly display sequential amide and (15)N correlations along the polypeptide chain, and third, the sign patterns of the diagonal and the sequential peaks originating from any residue are dependent on the nature of the adjacent residues, especially the glycines and the prolines. These lead to so-called "triplet fixed points" which serve as starting points and/or check points during the course of sequential walks, and explicit side chains assignment becomes less crucial for unambiguous backbone assignment. These features significantly enhance the speed of data analysis, reduce the amount of experimentation required, and thus result in a substantially faster and unambiguous assignment. Following the amide and (15)N assignments, the other proton and carbon assignments can be obtained in a straightforward manner, from the well-established three-dimensional triple-resonance experiments. We have successfully tested the new approach with different proteins in the molecular mass range of 10-22 kDa, and for illustration, we present here the backbone results on the HIV-1 protease-tethered dimer (molecular mass approximately 22 kDa), both in the folded and in the unfolded forms, the two ends of the folding funnel. We believe that the new assignment approach will be of great value for both structural genomics and protein folding research by NMR.     

Time-resolved NMR studies of RNA folding

The application of real-time NMR experiments to the study of RNA folding, as reviewed in this article, is relatively new. For many RNA folding events, current investigations suggest that the time scales are in the second to minute regime. In addition, the initial investigations suggest that different folding rates are observed for one structural transition may be due to the hierarchical folding units of RNA. Many of the experiments developed in the field of NMR of protein folding cannot directly be transferred to RNA: hydrogen exchange experiments outside the spectrometer cannot be applied since the intrinsic exchange rates are too fast in RNA, relaxation dispersion experiments on the other require faster structural transitions than those observed in RNA. On the other hand, information derived from time-resolved NMR experiments, namely the acquisition of native chemical shifts, can be readily interpreted in light of formation of a single long-range hydrogen bonding interaction. Together with mutational data that can readily be obtained for RNA and new ligation technologies that enhance site resolution even further, time-resolved NMR may become a powerful tool to decipher RNA folding. Such understanding will be of importance to understand the functions of coding and non-coding RNAs in cells.     

Highly compact folding of chromatin induced by cellular cation concentrations. Evidence from atomic force microscopy studies in aqueous solution

We have performed a very extensive investigation of chromatin folding in different buffers over a wide range of ionic conditions similar to those found in eukaryotic cells. Our results show that in the presence of physiological concentrations of monovalent cations and/or low concentrations of divalent cations, small chicken erythrocyte chromatin fragments and chromatin from HeLa cells observed by transmission electron microscopy (TEM) show a compact folding, forming circular bodies of approximately 35 nm in diameter that were found previously in our laboratory in studies performed under very limited conditions. Since TEM images are obtained with dehydrated samples, we have performed atomic force microscopy (AFM) experiments to analyze chromatin structure in the presence of solutions containing different cation concentrations. The highly compact circular structures (in which individual nucleosomes are not visible as separated units) produced by small chromatin fragments in interphase ionic conditions observed by AFM are equivalent to the structures observed by TEM with chromatin samples prepared under the same ionic conditions. We have also carried out experiments of sedimentation and trypsin digestion of chromatin fragments; the results obtained confirm our AFM observations. Our results suggest that the compaction of bulk interphase chromatin in solution at room temperature is considerably higher than that generally considered in current literature. The dense chromatin folding observed in this study is consistent with the requirement of compact chromatin structures as starting elements for the building of metaphase chromosomes, but poses a difficult physical problem for gene expression during interphase.     

The "first in-last out" hypothesis on protein folding revisited

We calculated profiles for mean residue depth, contact order, and number of contacts in the native structure of a series of proteins for which folding has been studied extensively, the chymotrypsin inhibitor 2, the SH3 module from the src tyrosine kinase, the small ribonuclease barnase, the bacterial immunity protein Im7, and apomyoglobin. We compared these profiles with experimental data from equilibrium or pulse labeling hydrogen-deuterium exchange obtained from NMR and phi values obtained from the protein engineering approach. We find a good qualitative agreement between the hierarchy of formation of topological elements during the folding process and the ranking of secondary structure elements in terms of residue depth. Residues that are most deeply buried in the core of the native protein usually belong to stretches of secondary structure elements that are formed early in the folding pathway. Residue depth can thus provide a useful and simple tool for the design of folding experiments.     

Assortative mixing in Protein Contact Networks and protein folding kinetics

MOTIVATION: Starting from linear chains of amino acids, the spontaneous folding of proteins into their elaborate 3D structures is one of the remarkable examples of biological self-organization. We investigated native state structures of 30 single-domain, two-state proteins, from complex networks perspective, to understand the role of topological parameters in proteins' folding kinetics, at two length scales--as 'Protein Contact Networks (PCNs)' and their corresponding 'Long-range Interaction Networks (LINs)' constructed by ignoring the short-range interactions. RESULTS: Our results show that, both PCNs and LINs exhibit the exceptional topological property of 'assortative mixing' that is absent in all other biological and technological networks studied so far. We show that the degree distribution of these contact networks is partly responsible for the observed assortativity. The coefficient of assortativity also shows a positive correlation with the rate of protein folding at both short- and long-contact scale, whereas, the clustering coefficients of only the LINs exhibit a negative correlation. The results indicate that the general topological parameters of these naturally evolved protein networks can effectively represent the structural and functional properties required for fast information transfer among the residues facilitating biochemical/kinetic functions, such as, allostery, stability and the rate of folding. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.     

COCO: a simple tool to enrich the representation of conformational variability in NMR structures

NMR structures are typically deposited in databases such as the PDB in the form of an ensemble of structures. Generally, each of the models in such an ensemble satisfies the experimental data and is equally valid. No unique solution can be calculated because the experimental NMR data is insufficient, in part because it reflects the conformational variability and dynamical behavior of the molecule in solution. Even for relatively rigid molecules, the limited number of structures that are typically deposited cannot completely encompass the structural diversity allowed by the observed NMR data, but they can be chosen to try and maximize its representation. We describe here the adaptation and application of techniques more commonly used to examine large ensembles from molecular dynamics simulations, to the analysis of NMR ensembles. The approach, which is based on principal component analysis, we call COCO ("Complementary Coordinates"). The COCO approach analyses the distribution of an NMR ensemble in conformational space, and generates a new ensemble that fills "gaps" in the distribution. The method is very rapid, and analysis of a 25-member ensemble and generation of a new 25 member ensemble typically takes 1-2 min on a conventional workstation. Applied to the 545 structures in the RECOORD database, we find that COCO generates new ensembles that are as structurally diverse-both from each other and from the original ensemble-as are the structures within the original ensemble. The COCO approach does not explicitly take into account the NMR restraint data, yet in tests on selected structures from the RECOORD database, the COCO ensembles are frequently good matches to this data, and certainly are structures that can be rapidly refined against the restraints to yield high-quality, novel solutions. COCO should therefore be a useful aid in NMR structure refinement and in other situations where a richer representation of conformational variability is desired-for example in docking studies. COCO is freely accessible via the website www.ccpb.ac.uk/COCO.     

Folding dynamics and conformational heterogeneity of human telomeric G-quadruplex structures in Na+ solutions by single molecule FRET microscopy

G-quadruplex structures can occur throughout the genome, including at telomeres. They are involved in cellular regulation and are potential drug targets. Human telomeric G-quadruplex structures can fold into a number of different conformations and show large conformational diversity. To elucidate the different G-quadruplex conformations and their dynamics, we investigated telomeric G-quadruplex folding using single molecule FRET microscopy in conditions where it was previously believed to yield low structural heterogeneity. We observed four FRET states in Na⁺ buffers: an unfolded state and three G-quadruplex related states that can interconvert between each other. Several of these states were almost equally populated at low to medium salt concentrations. These observations appear surprising as previous studies reported primarily one G-quadruplex conformation in Na⁺ buffers. Our results permit, through the analysis of the dynamics of the different observed states, the identification of a more stable G-quadruplex conformation and two transient G-quadruplex states. Importantly these results offer a unique view into G-quadruplex topological heterogeneity and conformational dynamics.     

Requirements on paramagnetic relaxation enhancement data for membrane protein structure determination by NMR

Nuclear magnetic resonance (NMR) structure calculations of the α-helical integral membrane proteins DsbB, GlpG, and halorhodopsin show that distance restraints from paramagnetic relaxation enhancement (PRE) can provide sufficient structural information to determine their structure with an accuracy of about 1.5?? in the absence of other long-range conformational restraints. Our systematic study with simulated NMR data shows that about one spin label per transmembrane helix is necessary for obtaining enough PRE distance restraints to exclude wrong topologies, such as pseudo mirror images, if only limited other NMR restraints are available. Consequently, an experimentally realistic amount of PRE data enables α-helical membrane protein structure determinations that would not be feasible with the very limited amount of conventional NOESY data normally available for these systems. These findings are in line with our recent first de novo NMR structure determination of a heptahelical integral membrane protein, proteorhodopsin, that relied extensively on PRE data.     

The metaphase scaffold is helically folded: sister chromatids have predominantly opposite helical handedness

We have studied the three-dimensional folding of the scaffolding in histone H1-depleted chromosomes by immunofluorescence with an antibody specific for topoisomerase II. Two different types of decondensed chromosomes are observed. The majority of the chromosomes are expanded, and the central fluorescence signal is surrounded by a large halo of chromatin. A much smaller number of chromosomes are more compact in length; they contain a smaller halo of chromatin and their scaffolds are not extended but folded into a genuine, quite regular helical coil. This conclusion is based on a three-dimensional structural analysis by optical sectioning. The number of helical coils is related to chromosome length. Surprisingly, sister chromatids have predominantly opposite helical handedness; that is, they are related by mirror symmetry.     

OB(oligonucleotide/oligosaccharide binding)-fold: common structural and functional solution for non-homologous sequences.

A novel folding motif has been observed in four different proteins which bind oligonucleotides or oligosaccharides: staphylococcal nuclease, anticodon binding domain of asp-tRNA synthetase and B-subunits of heat-labile enterotoxin and verotoxin-1. The common fold of the four proteins, which we call the OB-fold, has a five-stranded beta-sheet coiled to form a closed beta-barrel. This barrel is capped by an alpha-helix located between the third and fourth strands. The barrel-helix frameworks can be superimposed with r.m.s. deviations of 1.4-2.2 A, but no similarities can be observed in the corresponding alignment of the four sequences. The nucleotide or sugar binding sites, known for three of the four proteins, are located in nearly the same position in each protein: on the side surface of the beta-barrel, where three loops come together. Here we describe the determinants of the OB-fold, based on an analysis of all four structures. These proposed determinants explain how very different sequences adopt the OB-fold. They also suggest a reinterpretation of the controversial structure of gene 5 ssDNA binding protein, which exhibits some topological and functional similarities with the OB-fold proteins.     

An experimental investigation of conformational fluctuations in proteins G and L

The B1 domains of streptococcal proteins G and L are structurally similar, but they have different sequences and they fold differently. We have measured their NMR spectra at variable temperature using a range of concentrations of denaturant. Many residues have curved amide proton temperature dependence, indicating that they significantly populate alternative, locally unfolded conformations. The results, therefore, provide a view of the locations of low-lying, locally unfolded conformations. They indicate approximately 4-6 local minima for each protein, all within ca. 2.5 kcal/mol of the native state, implying a locally rough energy landscape. Comparison with folding data for these proteins shows that folding involves most molecules traversing a similar path, once a transition state containing a beta hairpin has been formed, thereby defining a well-populated pathway down the folding funnel. The hairpin that directs the folding pathway differs for the two proteins and remains the most stable part of the folded protein.     

Folding simulations and computer redesign of protein A three-helix bundle motifs

In solution, the B domain of protein A from Staphylococcus aureus (B domain) possesses a three-helix bundle structure. This simple motif has been previously reproduced by Kolinski and Skolnick (Proteins 18: 353–366, 1994) using a reduced representation lattice model of proteins with a statistical interaction scheme. In this paper, an improved version of the potential has been used, and the robustness of this result has been tested by folding from the random state a set of three-helix bundle proteins that are highly homologous to the B domain of protein A. Furthermore, an attempt to redesign the B domain native structure to its topological mirror image fold has been made by multiple mutations of the hydrophobic core and the turn region between helices I and II. A sieve method for scanning a large set of mutations to search for this desired property has been proposed. It has been shown that mutations of native B domain hydrophobic core do not introduce significant changes in the protein motif. Mutations in the turn region were also very conservative; nevertheless, a few mutants acquired the desired topological mirror image motif. A set of all atom models of the most probable mutant was reconstructed from the reduced models and refined using a molecular dynamics algorithm in the presence of water. The packing of all atom structures obtained corroborates the lattice model results. We conclude that the change in the handedness of the turn induced by the mutations, augmented by the repacking of hydrophobic core and the additional burial of the second helix N-cap side chain, are responsible for the predicted preferential adoption of the mirror image structure. © 1996 Wiley-Liss, Inc.     

Modeling disordered regions in proteins using Rosetta

Protein structure prediction methods such as Rosetta search for the lowest energy conformation of the polypeptide chain. However, the experimentally observed native state is at a minimum of the free energy, rather than the energy. The neglect of the missing configurational entropy contribution to the free energy can be partially justified by the assumption that the entropies of alternative folded states, while very much less than unfolded states, are not too different from one another, and hence can be to a first approximation neglected when searching for the lowest free energy state. The shortcomings of current structure prediction methods may be due in part to the breakdown of this assumption. Particularly problematic are proteins with significant disordered regions which do not populate single low energy conformations even in the native state. We describe two approaches within the Rosetta structure modeling methodology for treating such regions. The first does not require advance knowledge of the regions likely to be disordered; instead these are identified by minimizing a simple free energy function used previously to model protein folding landscapes and transition states. In this model, residues can be either completely ordered or completely disordered; they are considered disordered if the gain in entropy outweighs the loss of favorable energetic interactions with the rest of the protein chain. The second approach requires identification in advance of the disordered regions either from sequence alone using for example the DISOPRED server or from experimental data such as NMR chemical shifts. During Rosetta structure prediction calculations the disordered regions make only unfavorable repulsive contributions to the total energy. We find that the second approach has greater practical utility and illustrate this with examples from de novo structure prediction, NMR structure calculation, and comparative modeling.     

Specific non-native hydrophobic interactions in a hidden folding intermediate: implications for protein folding

Structures of intermediates and transition states in protein folding are usually characterized by amide hydrogen exchange and protein engineering methods and interpreted on the basis of the assumption that they have native-like conformations. We were able to stabilize and determine the high-resolution structure of a partially unfolded intermediate that exists after the rate-limiting step of a four-helix bundle protein, Rd-apocyt b(562), by multidimensional NMR methods. The intermediate has partial native-like secondary structure and backbone topology, consistent with our earlier native state hydrogen exchange results. However, non-native hydrophobic interactions exist throughout the structure. These and other results in the literature suggest that non-native hydrophobic interactions may occur generally in partially folded states. This can alter the interpretation of mutational protein engineering results in terms of native-like side chain interactions. In addition, since the intermediate exists after the rate-limiting step and Rd-apocyt b(562) folds very rapidly (k(f) approximately 10(4) s(-1)), these results suggest that non-native hydrophobic interactions, in the absence of topological misfolding, are repaired too rapidly to slow folding and cause the accumulation of folding intermediates. More generally, these results illustrate an approach for determining the high-resolution structure of folding intermediates.     

Amphipols for Each Season

Amphipols (APols) are short amphipathic polymers that can substitute for detergents at the transmembrane surface of membrane proteins (MPs) and, thereby, keep them soluble in detergent free aqueous solutions. APol-trapped MPs are, as a rule, more stable biochemically than their detergent-solubilized counterparts. APols have proven useful to produce MPs, most noticeably by assisting their folding from the denatured state obtained after solubilizing MP inclusion bodies in either SDS or urea. They facilitate the handling in aqueous solution of fragile MPs for the purpose of proteomics, structural and functional studies, and therapeutics. Because APols can be chemically labeled or functionalized, and they form very stable complexes with MPs, they can also be used to functionalize those indirectly, which opens onto many novel applications. Following a brief recall of the properties of APols and MP/APol complexes, an update is provided of recent progress in these various fields.     

Protein folding studied by real-time NMR spectroscopy

Real-time NMR spectroscopy developed to a generally applicable method to follow protein folding reactions. It combines the access to high resolution data with kinetic experiments allowing very detailed insights into the development of the protein structure during different steps of folding. The present review concentrates mainly on the progress of real-time NMR during the last 5 years. Starting from simple 1D experiments, mainly changes of the chemical shifts and line widths of the resonances have been used to analyze the different states populated during the folding reactions. Today, we have a broad spectrum of 1D, 2D, and even 3D NMR methods focusing on different characteristics of the folding polypeptide chains. More than 20 proteins have been investigated so far by these time-resolved experiments and the main results and conclusions are discussed in this report. Real-time NMR provides comprehensive contributions for joining experiment and theory within the 'new view' of protein folding.     

Predicting knot and catenane type of products of site-specific recombination on twist knot substrates

Site-specific recombination on supercoiled circular DNA molecules can yield a variety of knots and catenanes. Twist knots are some of the most common conformations of these products, and they can act as substrates for further rounds of site-specific recombination. They are also one of the simplest families of knots and catenanes. Yet, our systematic understanding of their implication in DNA and important cellular processes such as site-specific recombination is very limited. Here, we present a topological model of site-specific recombination characterizing all possible products of this reaction on twist knot substrates, extending the previous work of Buck and Flapan. We illustrate how to use our model to examine previously uncharacterized experimental data. We also show how our model can help determine the sequence of products in multiple rounds of processive recombination and distinguish between products of processive and distributive recombinations.This model studies generic site-specific recombination on arbitrary twist knot substrates, a subject for which there is limited global understanding. We also provide a systematic method of applying our model to a variety of different recombination systems.     

Fast and faster: a designed variant of the B-domain of protein A folds in 3 microsec

We have introduced the mutation glycine 29 to alanine, designed to increase the rate of protein folding, into the B-domain of protein A (BdpA). From NMR lineshape analysis, we find the G29A mutation increases the folding rate constant by threefold; the folding time is 3 microsec. Although wild-type BdpA folds extremely fast, simple-point mutations can still speed up the folding; thus, the folding rate is not evolutionarily maximized. The short folding time of G29A BdpA (the shortest time yet reported) makes it an attractive candidate for an all-atom molecular dynamics simulation that could potentially show a complete folding reaction starting from an extended chain. We also constructed a fluorescent variant of BdpA by mutating phenylalanine 13 to tryptophan, allowing fluorescence-based time-resolved temperature-jump measurements. Temperature jumps and NMR complement each other, and give a very complete picture of the folding kinetics.