Physicochemical evaluation of protein folds predicted by threading

Protein structure prediction remains an unsolved problem. Since prediction of the native structure seems very difficult, one usually tries to predict the correct fold of a protein. Here the "fold" is defined by the approximate backbone structure of the protein. However, physicochemical factors that determine the correct fold are not well understood. It has recently been reported that molecular mechanics energy functions combined with effective solvent terms can discriminate the native structures from misfolded ones. Using such a physicochemical energy function, we studied the factors necessary for discrimination of correct and incorrect folds. We first selected correct and incorrect folds by a conventional threading method. Then, all-atom models of those folds were constructed by simply minimizing the atomic overlaps. The constructed correct model representing the native fold has almost the same backbone structure as the native structure but differs in side-chain packing. Finally, the energy values of the constructed models were compared with that of the experimentally determined native structure. The correct model as well as the native structure showed lower energy than misfolded models. However, a large energy gap was found between the native structure and the correct model. By decomposing the energy values into their components, it was found that solvent effects such as the hydrophobic interaction or solvent shielding and the Born energy stabilized the correct model rather than the native structure. The large energetic stabilization of the native structure was attained by specific side-chain packing. The stabilization by solvent effects is small compared to that by side-chain packing. Therefore, it is suggested that in order to confidently predict the correct fold of a protein, it is also necessary to predict correct side-chain packing.     

Protein conformational transitions coupled to binding in molecular recognition of unstructured proteins: deciphering the effect of intermolecular interactions on computational structure prediction of the p27Kip1 protein bound to the cyclin A-cyclin-dependent kinase 2 complex

The relationship between folding mechanism coupled to binding and structure prediction of the tertiary complexes is studied for the p27(Kip) (1) protein which has an intrinsically disordered unbound form and undergoes a functional folding transition during complex formation with the phosphorylated cyclin A-cyclin-dependent kinase 2 (Cdk2) binary complex. Hierarchy of p27(Kip1) structural loss determined in our earlier studies from temperature-induced Monte Carlo simulations and subsequent characterization of the transition state ensemble (TSE) for the folding reaction have shown that simultaneous ordering of the p27(Kip1) native intermolecular interface for the beta-hairpin and beta-strand secondary structure elements is critical for nucleating a rapid kinetic transition to the native tertiary complex. In the present study, we investigate the effect of forming specific intermolecular interactions on structure prediction of the p27(Kip1) tertiary complex. By constraining different secondary structure elements of p27(Kip1) in their native bound conformations and conducting multiple simulated annealing simulations, we analyze differences in the success rate of predicting the native structure of p27(Kip1) in the tertiary complex. In accordance with the nucleation-condensation mechanism, we have found that further stabilization of the native intermolecular interface for the beta-hairpin and beta-strand elements of p27(Kip1), that become ordered in the TSE, but are hardly populated in the unbound state, results in a consistent acquisition of the native bound structure. Conversely, the excessive stablization of the local secondary structure elements, which are rarely detected in the TSE, has a detrimental effect on convergence to the native bound structure.     

Rescore protein–protein docked ensembles with an interface contact statistics

The recently developed statistical measure for the type of residue–residue contact at protein complex interfaces, based on a parameter‐free definition of contact, has been used to define a contact score that is correlated with the likelihood of correctness of a proposed complex structure. Comparing the proposed contact scores on the native structure and on a set of model structures the proposed measure was shown to generally favor the native structure but in itself was not able to reliably score the native structure to be the best. Adjusting the scores of redocking experiments with the contact score showed that the adjusted score was able to move up the ranking of the native‐like structure among the proposed complexes when the native‐like was not ranked the best by the respective program. Tests on docking of unbound proteins compared the contact scores of the complexes with the contact score of the crystal structure again showing the tendency of the contact score to favor native‐like conformations. The possibility of using the contact score to improve the determination of biological dimers in a crystal structure was also explored. Proteins 2017; 85:235–241. © 2016 Wiley Periodicals, Inc.     

Non‐native species modify the isotopic structure of freshwater fish communities across the globe

Multiple anthropogenic pressures including the widespread introductions of non‐native species threaten biodiversity and ecosystem functioning notably by modifying the trophic structure of communities. Here, we provided a global evaluation of the impacts of non‐native species on the isotopic structure (δ¹³C and δ¹⁵N) of freshwater fish communities. We gathered the stable isotope values (n = 4030) of fish species in 496 fish communities in lentic (lakes, backwaters, reservoirs) and lotic (running waters such as streams, rivers) ecosystems throughout the world and quantified the isotopic structure of communities. Overall, we found that communities containing non‐native species had a different isotopic structure than communities without non‐native species. However, these differences varied between ecosystem types and the trophic positions of non‐native species. In lotic ecosystems, communities containing non‐native species had a larger total isotopic niche than communities without non‐native species. This was primarily driven by the addition of non‐native predators at the top of the food chain that increased δ¹⁵N range without modifying the isotopic niche size of native species. In lentic ecosystems, non‐native primary consumers increased δ¹⁵N range and this was likely driven by an increase of resource availability for species at higher trophic levels, increasing food chain length. The introduction of non‐native secondary consumers at the centre of the isotopic niche of recipient communities decreased the core isotopic niche size, the δ¹³C range of recipient communities and the total isotopic niche of coexisting native species. These results suggested a modified contribution of the basal resources consumed (e.g. multi‐chain omnivory) and an increase level of competition with native species. Our results notably imply that, by affecting the isotopic structure of freshwater fish communities at a global scale, non‐native species represent an important source of perturbations that should be accounted for when investigating macro‐ecological patterns of community structure and biotic interactions.     

How well can we predict native contacts in proteins based on decoy structures and their energies?

One strategy for ab initio protein structure prediction is to generate a large number of possible structures (decoys) and select the most fitting ones based on a scoring or free energy function. The conformational space of a protein is huge, and chances are rare that any heuristically generated structure will directly fall in the neighborhood of the native structure. It is desirable that, instead of being thrown away, the unfitting decoy structures can provide insights into native structures so prediction can be made progressively. First, we demonstrate that a recently parameterized physics-based effective free energy function based on the GROMOS96 force field and a generalized Born/surface area solvent model is, as several other physics-based and knowledge-based models, capable of distinguishing native structures from decoy structures for a number of widely used decoy databases. Second, we observe a substantial increase in correlations of the effective free energies with the degree of similarity between the decoys and the native structure, if the similarity is measured by the content of native inter-residue contacts in a decoy structure rather than its root-mean-square deviation from the native structure. Finally, we investigate the possibility of predicting native contacts based on the frequency of occurrence of contacts in decoy structures. For most proteins contained in the decoy databases, a meaningful amount of native contacts can be predicted based on plain frequencies of occurrence at a relatively high level of accuracy. Relative to using plain frequencies, overwhelming improvements in sensitivity of the predictions are observed for the 4_state_reduced decoy sets by applying energy-dependent weighting of decoy structures in determining the frequency. There, approximately 80% native contacts can be predicted at an accuracy of approximately 80% using energy-weighted frequencies. The sensitivity of the plain frequency approach is much lower (20% to 40%). Such improvements are, however, not observed for the other decoy databases. The rationalization and implications of the results are discussed.     

Flexibility and bioactivity of insulin: an NMR investigation of the solution structure and folding of an unusually flexible human insulin mutant with increased biological activity

The structure and folding of a novel human insulin mutant, [Thr(B27) --> Pro, Pro(B28) --> Thr]insulin (PT insulin), in aqueous solution and in mixtures of water and 2,2,2-trifluoroethanol (TFE) have been studied by NMR spectroscopy. It was found that PT insulin has a highly flexible structure in pure water and is present in at least two different conformations, although with an overall tertiary structure similar to that of native insulin. Furthermore, the native helical structures are poorly defined. Surprisingly, the mutant has a biological activity about 50% higher than native insulin. In contrast, in TFE/water solution the mutant reveals a propensity of forming a well-defined structure at the secondary structure level, similar to monomeric native insulin. Thus, as shown by a detailed determination of the structure from 208 distance restraints and 52 torsion angle restraints by distance geometry, simulated annealing, and restrained energy minimization, the native insulin helices (A2-A7, A13-A19, and B10-B19) as well as the beta-turn (B20-B23) are formed in 35% TFE. However, the amount of tertiary structure is decreased significantly in TFE/water solution. The obtained results suggest that only an overall tertiary fold, as observed for PT insulin in pure water, is necessary for expressing the biological activity of insulin, as long as the molecule is flexible and retains the propensity to form the secondary structure required for its receptor binding. In contrast, a compact secondary structure, as found for native insulin in solution, is unnecessary for the biological activity. A model for the receptor binding of insulin is suggested that relates the increased bioactivity to the enhanced flexibility of the mutant.     

Solvent accessibility, protein surfaces, and protein folding.

Studies of the native structures of proteins, together with measurements of the thermodynamic properties of the transition between unfolded and native states, have defined the major components of the forces that stabilize native protein structures. However, the nature of the intermediates in the folding process remains largely hypothetical. It is a fairly widespread and not implausible assumption that the intermediates in the folding of a monomeric protein contain the same kinds of secondary and tertiary structures that appear in the native conformation, and that, although unstable, their lifetimes are prolonged by forces similar to those that stabilize the native structure. We wished to examine what happens if, during the folding of a monomeric protein, regions of secondary structure come together to form an intermediate of reduced instability. We applied calculations of accessible surface area (a measure of hydrophobic stabilization) and parameterized nonbonded energy calculations (measuring the strengths of van der Waals forces) to identify the kinds of stabilizing interactions that might be available to such an intermediate. First, we analyzed the total buried surface area of two types of proteins into contributions from formation of secondary structure alone, interaction of pairs of secondary-structural elements, the formation of the structure alone, interaction of pairs of secondary-structural elements, the formation of the complete secondary structure without the turns, and the complete native structure. The formation of secondary structure alone, without tertiary-structural interactions, buries roughly half the surface that the complete structure does. We then analyzed in more detail the approach of two alpha-helices to form a complex, as an illustrative example of the nature of the interaction between compact structural units which remain fairly rigid during their interaction. Many features of the results are not limited to the interaction of alpha-helices. (The results therefore neither confirm nor refute the hypothesis that alpha-helices are intermediates in the folding proteins). We find that the first forces to be felt upon approach arise from solvent conditions on the relative position and orientation of the two helices as does the close packing which optimizes the van der Waals interactions at shorter distances apart. Therefore there appears to be a range of distances in which hydrophobic interactions could create a nonspecific complex between two helices in which the side chains might have sufficient time to seek the proper interdigitation observed in the native structure, where the two helices are in intimate contact. Indeed, we find that only in the final stages of approach is the native geometry the most stable; in the region in which solvent-exclusion forces predominate, the conformation with helix axes parallel is more stable than the native conformation, in the cases we examined...     

The paradoxical behavior of a highly structured misfolded intermediate in RNA folding

Like many structured RNAs, the Tetrahymena group I ribozyme is prone to misfolding. Here we probe a long-lived misfolded species, referred to as M, and uncover paradoxical aspects of its structure and folding. Previous work indicated that a non-native local secondary structure, termed alt P3, led to formation of M during folding in vitro. Surprisingly, hydroxyl radical footprinting, fluorescence measurements with site-specifically incorporated 2-aminopurine, and functional assays indicate that the native P3, not alt P3, is present in the M state. The paradoxical behavior of alt P3 presumably arises because alt P3 biases folding toward M, but, after commitment to this folding pathway and before formation of M, alt P3 is replaced by P3. Further, structural and functional probes demonstrate that the misfolded ribozyme contains extensive native structure, with only local differences between the two states, and the misfolded structure even possesses partial catalytic activity. Despite the similarity of these structures, re-folding of M to the native state is very slow and is strongly accelerated by urea, Na+, and increased temperature and strongly impeded by Mg2+ and the presence of native peripheral contacts. The paradoxical observations of extensive native structure within the misfolded species but slow conversion of this species to the native state are readily reconciled by a model in which the misfolded state is a topological isomer of the native state, and computational results support the feasibility of this model. We speculate that the complex topology of RNA secondary structures and the inherent rigidity of RNA helices render kinetic traps due to topological isomers considerably more common for RNA than for proteins.     

Denaturation and reassembly of chaperonin GroEL studied by solution X-ray scattering

We measured the denaturation and reassembly of Escherichia coli chaperonin GroEL using small-angle solution X-ray scattering, which is a powerful technique for studying the overall structure and assembly of a protein in solution. The results of the urea-induced unfolding transition show that GroEL partially dissociates in the presence of more than 2 M urea, cooperatively unfolds at around 3 M urea, and is in a monomeric random coil-like unfolded structure at more than 3.2 M urea. Attempted refolding of the unfolded GroEL monomer by a simple dilution procedure is not successful, leading to formation of aggregates. However, the presence of ammonium sulfate and MgADP allows the fully unfolded GroEL to refold into a structure with the same hydrodynamic dimension, within experimental error, as that of the native GroEL. Moreover, the X-ray scattering profiles of the GroEL thus refolded and the native GroEL are coincident with each other, showing that the refolded GroEL has the same structure and the molecular mass as the native GroEL. These results demonstrate that the fully unfolded GroEL monomer can refold and reassemble into the native tetradecameric structure in the presence of ammonium sulfate and MgADP without ATP hydrolysis and preexisting chaperones. Therefore, GroEL can, in principle, fold and assemble into the native structure according to the intrinsic characteristic of its polypeptide chain, although preexisting GroEL would be important when the GroEL folding takes place under in vivo conditions, in order to avoid misfolding and aggregation.     

What is the shape of the distribution of protein conformations at equilibrium?

According to the thermodynamic hypothesis, the native state of proteins is that in which the free energy of the system is at its lowest, so that at normal temperature and pressure, proteins evolve to that state. We selected four proteins representative of each of the four classes, and for each protein make four simulations, one starting from the native structure and the other three starting from the structure obtained by threading the sequence of one protein onto the native backbone fold of the other three proteins. Because of their large conformational distances with respect to the native structure, the three alternative initial structures cannot be considered as local minima within the native ensemble of the corresponding protein. As expected, the initial native states are preserved in the .5?μs simulations performed here and validate the simulations. On the other hand, when the initial state is not native, an analysis of the trajectories does not reveal any evolution towards the native state, during that time. These results indicate that the distribution of protein conformations is multipeak shaped, so that apart from the peak corresponding to the native state, there are other peaks associated with average structures that are very different from the native and that can last as long as the native state.     

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.     

Synthesis, structure and imaging of oligodeoxyribonucleotides with tellurium-nucleobase derivatization

We report here the first synthesis of 5-phenyl-telluride-thymidine derivatives and the Te-phosphoramidite. We also report here the synthesis, structure and STM current-imaging studies of DNA oligonucleotides containing the nucleobases (thymine) derivatized with 5-phenyl-telluride functionality (5-Te). Our results show that the 5-Te-DNA is stable, and that the Te-DNA duplex has the thermo-stability similar to the corresponding native duplex. The crystal structure indicates that the 5-Te-DNA duplex structure is virtually identical to the native one, and that the Te-modified T and native A interact similarly to the native T and A pair. Furthermore, while the corresponding native showed weak signals, the DNA duplex modified with electron-rich tellurium functionality showed strong topographic and current peaks by STM imaging, suggesting a potential strategy to directly image DNA without structural perturbation.     

New pathways in folding of the Tetrahymena group I RNA enzyme.

Previous studies have shown that the earliest detectable step in folding of the Tetrahymena ribozyme is tertiary structure formation of the peripheral element P5abc. This, along with other results, has suggested that P5abc may serve as a scaffold upon which additional tertiary structure is built. Herein we use the onset of oligonucleotide cleavage activity as a readout for native state formation and investigate the effect of P5abc on the rate of folding to the native structure. Despite the early folding of P5abc, its removal to give the E delta P5abc variant decreases the rate of attainment of an active structure less than fivefold (20-100 mM Mg2+, 15-50 degrees C). Furthermore, P5abc added in trans is able to bind the folded E delta P5abc ribozyme and promote oligonucleotide cleavage at least tenfold more rapidly than folding of the wild-type ribozyme, indicating that E delta P5abc does not have to first unfold before productively binding P5abc to form the true native state. This suggests that a state with the overall tertiary structure formed but with P5abc unfolded represents a viable on-pathway intermediate for the wild-type ribozyme. These results provide strong evidence for the existence of two pathways to the native state: in one pathway P5abc forms tertiary structure first, and in another it forms late. The pathway in which P5abc forms first is favored because P5abc can fold quickly and because its tertiary structure is stable in the absence of additional structured elements, not because P5abc formation is required for subsequent folding steps. In the course of these experiments, we also found that most of the ribozyme population does not reach the native state directly under standard conditions in vitro, but instead forms an inactive structure that is stable for hours. Finally, the fraction that does fold to the native state folds with a single rate constant of 1 min-1, suggesting that there are no significantly populated "fast-track" pathways that reach the native state directly by avoiding slow folding steps.     

D2N: Distance to the native

Root-mean-square-deviation (RMSD), of computationally-derived protein structures from experimentally determined structures, is a critical index to assessing protein-structure-prediction-algorithms (PSPAs). The development of PSPAs to obtain 0 Å RMSD from native structures is considered central to computational biology. However, till date it has been quite challenging to measure how far a predicted protein structure is from its native — in the absence of a known experimental/native structure. In this work, we report the development of a metric “D2N” (distance to the native) — that predicts the “RMSD” of any structure without actually knowing the native structure. By combining physico-chemical properties and known universalities in spatial organization of soluble proteins to develop D2N, we demonstrate the ability to predict the distance of a proposed structure to within ± 1.5 ? error with a remarkable average accuracy of 93.6% for structures below 5 ? from the native. We believe that this work opens up a completely new avenue towards assigning reliable structures to whole proteomes even in the absence of experimentally determined native structures. The D2N tool is freely available at http://www.scfbio-iitd.res.in/software/d2n.jsp.     

Capturing native/native like structures with a physico-chemical metric (pcSM) in protein folding

Specification of the three dimensional structure of a protein from its amino acid sequence, also called a “Grand Challenge” problem, has eluded a solution for over six decades. A modestly successful strategy has evolved over the last couple of decades based on development of scoring functions (e.g. mimicking free energy) that can capture native or native-like structures from an ensemble of decoys generated as plausible candidates for the native structure. A scoring function must be fast enough in discriminating the native from unfolded/misfolded structures, and requires validation on a large data set(s) to generate sufficient confidence in the score. Here we develop a scoring function called pcSM that detects true native structure in the top 5 with 93% accuracy from an ensemble of candidate structures. If we eliminate the native from ensemble of decoys then pcSM is able to capture near native structure (RMSD < = 5 ?) in top 10 with 86% accuracy. The parameters considered in pcSM are a C-alpha Euclidean metric, secondary structural propensity, surface areas and an intramolecular energy function. pcSM has been tested on 415 systems consisting 142,698 decoys (public and CASP—largest reported hitherto in literature). The average rank for the native is 2.38, a significant improvement over that existing in literature. In-silico protein structure prediction requires robust scoring technique(s). Therefore, pcSM is easily amenable to integration into a successful protein structure prediction strategy. The tool is freely available at http://www.scfbio-iitd.res.in/software/pcsm.jsp.     

Reconstitution of chromatin higher-order structure from histone H5 and depleted chromatin

Reconstitution of the 30 nm filament of chromatin from pure histone H5 and chromatin depleted of H1 and H5 has been studied using small-angle neutron-scattering. We find that depleted, or stripped, chromatin is saturated by H5 at the same stoichiometry as that of linker histone in native chromatin. The structure and condensation behavior of fully reconstituted chromatin is indistinguishable from that of native chromatin. Both native and reconstituted chromatin condense continuously as a function of salt concentration, to reach a limiting structure that has a mass per unit length of 6.4 nucleosomes per 11 nm. Stripped chromatin at all ionic strengths appears to be a 10 nm filament, or a random coil of nucleosomes. In contrast, both native and reconstituted chromatin have a quite different structure, showing that H5 imposes a spatial correlation between neighboring nucleosomes even at low ionic strength. Our data also suggest that five to seven contiguous nucleosomes must have H5 bound in order to be able to form a higher-order structure.     

Characterization of a partly folded protein by NMR methods: studies on the molten globule state of guinea pig alpha-lactalbumin

NMR spectroscopy has been used to investigate the structure of a partially folded state of a protein, the molten globule or A-state of alpha-lactalbumin. The 1H NMR spectrum of this species differs substantially from those of both the native and fully unfolded states, reflecting the intermediate level of order. The resolution in the spectrum is limited by the widespread overlap and substantial line widths of many of the resonances. Methods have therefore been developed that exploit the well-resolved spectrum of the native protein to probe indirectly the A-state. A number of resonances of the A-state have been found to be substantially shifted from their positions in the spectrum of the unfolded state and have been identified through magnetization transfer with the native state, under conditions where the two states are interconverting. The most strongly perturbed residues in the A-state were found to be among those that form a hydrophobic core to the native structure. A number of amides were found to be highly protected from solvent exchange in the A-state. These have been identified through pH-jump experiments, which label them in the spectrum of the native protein. They were found to occur mainly in segments that are helical in the native structure. These results enable a model of the A-state to be proposed in which significant conformational freedom exists but where specific elements of native-like structure are preserved.     

Conformational transition between native and reactive center cleaved forms of alpha 1-antitrypsin by Fourier transform infrared spectroscopy and small-angle neutron scattering

alpha 1-Antitrypsin (alpha 1-AT) is the best-characterized member of the serpin superfamily of plasma proteins. Protease inhibitor members of this family undergo a characteristic reactive-center cleavage during expression of their inhibitory activity. The physical basis of this transition in alpha 1-AT from the stressed native conformation to the more stable reactive center cleaved (split) form was studied by Fourier transform infrared (FT-IR) spectroscopy and neutron scattering. The FT-IR spectra show that, while split alpha 1-AT has three intense well-resolved components associated with the presence of antiparallel beta-sheet and alpha-helix conformations, the amide I band of native alpha 1-AT has only one intense component, associated with the presence of beta-sheet structure. 1H-2H exchange within the polypeptide backbone, studied by FT-IR and NMR spectroscopy, shows that the native form undergoes greater exchange than the split form. Under the same conditions, neutron scattering shows no differences in the radius of gyration RG of the native and the split forms. In contrast, in high concentrations of phosphate approaching those used for crystallization, the native form (unlike the split form) undergoes dimerization. These data indicate that the conformational transition largely involves localized secondary and tertiary structure rearrangements. We propose that the energetically stressed native alpha 1-AT structure is the consequence of a significantly reduced number of hydrogen bonds in secondary structure components and that reactive-site cleavage between Met358 and Ser359 is the key for the development of the fully hydrogen bonded more stable serpin structure.     

Enhanced specificity against misfolding in a thermostable mutant of the Tetrahymena ribozyme

Structured RNAs encode native conformations that are more stable than the vast ensembles of alternative conformations, but how this specificity is evolved is incompletely understood. Here we show that a variant of the Tetrahymena group I intron ribozyme that was generated previously by in vitro selection for enhanced thermostability also displays modestly enhanced specificity against a stable misfolded structure that is globally similar to the native state, despite the absence of selective pressure to increase the energy gap between these structures. The enhanced specificity for native folding arises from mutations in two nucleotides that are close together in space in the native structure, and additional experiments show that these two mutations do not affect the stability of the misfolded conformation relative to the largely unstructured transition state ensemble for interconversion between the native and misfolded conformers. Thus, they selectively stabilize the native state, presumably by strengthening a local tertiary contact network that cannot form in the misfolded conformation. The stabilization is larger in the presence of the peripheral element P5abc, suggesting that cooperative tertiary structure formation plays a key role in the enhanced stability. The increased specificity in the absence of explicit selection suggests that the large energy gap in the wild-type RNA may have arisen analogously, a consequence of selective pressure for stability of the functional structure. More generally, the structural rigidity and intricate networks of contacts in structured RNAs may allow them to evolve substantial structural specificity without explicit negative selection, even against closely related alternative structures.     

The denatured state of Engrailed Homeodomain under denaturing and native conditions

Protein folding starts from the elusive form of the denatured state that is present under conditions that favour the native state. We have studied the denatured state of Engrailed Homeodomain (En-HD) under mildly and strongly denaturing conditions at the level of individual residues by NMR and more globally by conventional spectroscopy and solution X-ray scattering. We have compared these states with a destabilized mutant, L16A, which is predominantly denatured under conditions where the wild-type is native. This engineered denatured state, which could be directly studied under native conditions, was in genuine equilibrium with the native state, which could be observably populated by changing the conditions or introducing a stabilizing mutation. The denatured state had extensive native secondary structure and was significantly compact and globular. But, the side-chains and backbone were highly mobile. Non-cooperative melting of the residual structure on the denatured state of En-HD was observed, both at the residue and the molecular level, with increasingly denaturing conditions. The absence of a co-operative transition could result from the denatured state ensemble progressing through a series of intermediates or from a more general slide (second-order transition) from the compact form under native conditions to the more extended at highly denaturing conditions. In either case, the starting point for folding under native conditions is highly structured and already poised to adopt the native structure.