Coevolution of function and the folding landscape: correlation with density of native contacts.
Ontology highlight
ABSTRACT:
SUBMITTER: Hills RD
PROVIDER: S-EPMC2567933 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:The recent availability of long equilibrium simulations of protein folding in atomistic detail for more than 10 proteins allows us to identify the key interactions driving folding. We find that the collective fraction of native amino acid contacts, Q, captures remarkably well the transition states for all the proteins with a folding free energy barrier. Going beyond this global picture, we devise two different measures to quantify the importance of individual interresidue contacts in the folding mechanism: (i) the log-ratio of lifetimes of contacts during folding transition paths and in the unfolded state and (ii) a Bayesian measure of how predictive the formation of each contact is for being on a transition path. Both of these measures indicate that native, or near-native, contacts are important for determining mechanism, as might be expected. More remarkably, however, we found that for almost all the proteins, with the designed protein α3D being a notable exception, nonnative contacts play no significant part in determining folding mechanisms.
| S-EPMC3816414 | biostudies-literature
Project description:The similarity in the three-dimensional structures of homologous proteins imposes strong constraints on their sequence variability. It has long been suggested that the resulting correlations among amino acid compositions at different sequence positions can be exploited to infer spatial contacts within the tertiary protein structure. Crucial to this inference is the ability to disentangle direct and indirect correlations, as accomplished by the recently introduced direct-coupling analysis (DCA). Here we develop a computationally efficient implementation of DCA, which allows us to evaluate the accuracy of contact prediction by DCA for a large number of protein domains, based purely on sequence information. DCA is shown to yield a large number of correctly predicted contacts, recapitulating the global structure of the contact map for the majority of the protein domains examined. Furthermore, our analysis captures clear signals beyond intradomain residue contacts, arising, e.g., from alternative protein conformations, ligand-mediated residue couplings, and interdomain interactions in protein oligomers. Our findings suggest that contacts predicted by DCA can be used as a reliable guide to facilitate computational predictions of alternative protein conformations, protein complex formation, and even the de novo prediction of protein domain structures, contingent on the existence of a large number of homologous sequences which are being rapidly made available due to advances in genome sequencing.
| S-EPMC3241805 | biostudies-literature
Project description:The NusG protein family is structurally and functionally conserved in all domains of life. Its members directly bind RNA polymerases and regulate transcription processivity and termination. RfaH, a divergent sub-family in its evolutionary history, is known for displaying distinct features than those in NusG proteins, which allows them to regulate the expression of virulence factors in enterobacteria in a DNA sequence-dependent manner. A striking feature is its structural interconversion between an active fold, which is the canonical NusG three-dimensional structure, and an autoinhibited fold, which is distinctively novel. How this novel fold is encoded within RfaH sequence to encode a metamorphic protein remains elusive. In this work, we used publicly available genomic RfaH protein sequences to construct a complete multiple sequence alignment, which was further augmented with metagenomic sequences and curated by predicting their secondary structure propensities using JPred. Coevolving pairs of residues were calculated from these sequences using plmDCA and GREMLIN, which allowed us to detect the enrichment of key metamorphic contacts after sequence filtering. Finally, we combined our coevolutionary predictions with molecular dynamics to demonstrate that these interactions are sufficient to predict the structures of both native folds, where coevolutionary-derived non-native contacts may play a key role in achieving the compact RfaH novel fold. All in all, emergent coevolutionary signals found within RfaH sequences encode the autoinhibited and active folds of this protein, shedding light on the key interactions responsible for the action of this metamorphic protein.
| S-EPMC9123645 | biostudies-literature
Project description:According to the 'thermodynamic hypothesis', the sequence of a biological macromolecule defines its folded, active (or 'native') structure as a global energy minimum in the folding landscape. However, the enormous complexity of folding landscapes of large macromolecules raises the question of whether there is in fact a unique global minimum corresponding to a unique native conformation or whether there are deep local minima corresponding to alternative active conformations. The folding of many proteins is well described by two-state models, leading to highly simplified representations of protein folding landscapes with a single native conformation. Nevertheless, accumulating experimental evidence suggests a more complex topology of folding landscapes with multiple active conformations that can take seconds or longer to interconvert. Here we demonstrate, using single-molecule experiments, that an RNA enzyme folds into multiple distinct native states that interconvert on a timescale much longer than that of catalysis. These data demonstrate that severe ruggedness of RNA folding landscapes extends into conformational space occupied by native conformations.
| S-EPMC2818749 | biostudies-literature
Project description:RfaH is a virulence factor from Escherichia coli whose C-terminal domain (CTD) undergoes a dramatic α-to-β conformational transformation. The CTD in its α-helical fold is stabilized by interactions with the N-terminal domain (NTD), masking an RNA polymerase binding site until a specific recruitment site is encountered. Domain dissociation is triggered upon binding to DNA, allowing the NTD to interact with RNA polymerase to facilitate transcription while the CTD refolds into the β-barrel conformation that interacts with the ribosome to activate translation. However, structural details of this transformation process in the context of the full protein remain to be elucidated. Here, we explore the mechanism of the α-to-β conformational transition of RfaH in the full-length protein using a dual-basin structure-based model. Our simulations capture several features described experimentally, such as the requirement of disruption of interdomain contacts to trigger the α-to-β transformation, confirms the roles of previously indicated residues E48 and R138, and suggests a new important role for F130, in the stability of the interdomain interaction. These native basins are connected through an intermediate state that builds up upon binding to the NTD and shares features from both folds, in agreement with previous in silico studies of the isolated CTD. We also examine the effect of RNA polymerase binding on the stabilization of the β fold. Our study shows that native-biased models are appropriate for interrogating the detailed mechanisms of structural rearrangements during the dramatic transformation process of RfaH.
| S-EPMC4521827 | biostudies-literature
Project description:Experiments have shown that the folding rate constants of two dozen structurally unrelated, small, single-domain proteins can be expressed in terms of one quantity (the contact order) that depends exclusively on the topology of the folded state. Such dependence is unique in chemical kinetics. Here we investigate its physical origin and derive the approximate formula ln(k) = ln(N) + a + bN, were N is the number of contacts in the folded state, and a and b are constants whose physical meaning is understood. This formula fits well the experimentally determined folding rate constants of the 24 proteins, with single values for a and b.
| S-EPMC122558 | biostudies-literature
Project description:The bacterial colicin-immunity proteins Im7 and Im9 fold by different mechanisms. Experimentally, at pH 7.0 and 10°C, Im7 folds in a three-state manner via an intermediate but Im9 folding is two-state-like. Accordingly, Im7 exhibits a chevron rollover, whereas the chevron arm for Im9 folding is linear. Here we address the biophysical basis of their different behaviors by using native-centric models with and without additional transferrable, sequence-dependent energies. The Im7 chevron rollover is not captured by either a pure native-centric model or a model augmented by nonnative hydrophobic interactions with a uniform strength irrespective of residue type. By contrast, a more realistic nonnative interaction scheme that accounts for the difference in hydrophobicity among residues leads simultaneously to a chevron rollover for Im7 and an essentially linear folding chevron arm for Im9. Hydrophobic residues identified by published experiments to be involved in nonnative interactions during Im7 folding are found to participate in the strongest nonnative contacts in this model. Thus our observations support the experimental perspective that the Im7 folding intermediate is largely underpinned by nonnative interactions involving large hydrophobics. Our simulation suggests further that nonnative effects in Im7 are facilitated by a lower local native contact density relative to that of Im9. In a one-dimensional diffusion picture of Im7 folding with a coordinate- and stability-dependent diffusion coefficient, a significant chevron rollover is consistent with a diffusion coefficient that depends strongly on native stability at the conformational position of the folding intermediate.
| S-EPMC4446218 | biostudies-literature
Project description:Escherichia coli RfaH activates gene expression by tethering the elongating RNA polymerase to the ribosome. This bridging action requires a complete refolding of the RfaH C-terminal domain (CTD) from an α-helical hairpin, which binds to the N-terminal domain (NTD) in the free protein, to a β-barrel, which interacts with the ribosomal protein S10 following RfaH recruitment to its target operons. The CTD forms a β-barrel when expressed alone or proteolytically separated from the NTD, indicating that the α-helical state is trapped by the NTD, perhaps co-translationally. Alternatively, the interdomain contacts may be sufficient to drive the formation of the α-helical form. Here, we use functional and NMR analyses to show that the denatured RfaH refolds into the native state and that RfaH in which the order of the domains is reversed is fully functional in vitro and in vivo. Our results indicate that all information necessary to determine its fold is encoded within RfaH itself, whereas accessory factors or sequential folding of NTD and CTD during translation are dispensable. These findings suggest that universally conserved RfaH homologs may change folds to accommodate diverse interaction partners and that context-dependent protein refolding may be widespread in nature.
| S-EPMC3905879 | biostudies-literature
Project description:Are all protein interactions fully optimized? Do suboptimal interactions compromise specificity? What is the functional impact of frustration? Why does evolution not optimize some contacts? Proteins and their complexes are best described as ensembles of states populating an energy landscape. These ensembles vary in breadth from narrow ensembles clustered around a single average X-ray structure to broader ensembles encompassing a few different functional "taxonomic" states on to near continua of rapidly interconverting conformations, which are called "fuzzy" or even "intrinsically disordered". Here we aim to provide a comprehensive framework for confronting the structural and dynamical continuum of protein assemblies by combining the concepts of energetic frustration and interaction fuzziness. The diversity of the protein structural ensemble arises from the frustrated conflicts between the interactions that create the energy landscape. When frustration is minimal after folding, it results in a narrow ensemble, but residual frustrated interactions result in fuzzy ensembles, and this fuzziness allows a versatile repertoire of biological interactions. Here we discuss how fuzziness and frustration play off each other as proteins fold and assemble, viewing their significance from energetic, functional, and evolutionary perspectives.We demonstrate, in particular, that the common physical origin of both concepts is related to the ruggedness of the energy landscapes, intramolecular in the case of frustration and intermolecular in the case of fuzziness. Within this framework, we show that alternative sets of suboptimal contacts may encode specificity without achieving a single structural optimum. Thus, we demonstrate that structured complexes may not be optimized, and energetic frustration is realized via different sets of contacts leading to multiplicity of specific complexes. Furthermore, we propose that these suboptimal, frustrated, or fuzzy interactions are under evolutionary selection and expand the biological repertoire by providing a multiplicity of biological activities. In accord, we show that non-native interactions in folding or interaction landscapes can cooperate to generate diverse functional states, which are essential to facilitate adaptation to different cellular conditions. Thus, we propose that not fully optimized structures may actually be beneficial for biological activities of proteins via an alternative set of suboptimal interactions. The importance of such variability has not been recognized across different areas of biology.This account provides a modern view on folding, function, and assembly across the protein universe. The physical framework presented here is applicable to the structure and dynamics continuum of proteins and opens up new perspectives for drug design involving not fully structured, highly dynamic protein assemblies.
| S-EPMC8023570 | biostudies-literature
Project description:Protein folding is described conceptually in terms of diffusion over a configurational free-energy landscape, typically reduced to a one-dimensional profile along a reaction coordinate. In principle, kinetic properties can be predicted directly from the landscape profile using Kramers theory for diffusive barrier crossing, including the folding rates and the transition time for crossing the barrier. Landscape theory has been widely applied to interpret the time scales for protein conformational dynamics, but protein folding rates and transition times have not been calculated directly from experimentally measured free-energy profiles. We characterized the energy landscape for native folding of the prion protein using force spectroscopy, measuring the change in extension of a single protein molecule at high resolution as it unfolded/refolded under tension. Key parameters describing the landscape profile were first recovered from the distributions of unfolding and refolding forces, allowing the diffusion constant for barrier crossing and the transition path time across the barrier to be calculated. The full landscape profile was then reconstructed from force-extension curves, revealing a double-well potential with an extended, partially unfolded transition state. The barrier height and position were consistent with the previous results. Finally, Kramers theory was used to predict the folding rates from the landscape profile, recovering the values observed experimentally both under tension and at zero force in ensemble experiments. These results demonstrate how advances in single-molecule theory and experiment are harnessing the power of landscape formalisms to describe quantitatively the mechanics of folding.
| S-EPMC3437844 | biostudies-other