Thermodynamic and Evolutionary Coupling between the Native and Amyloid State of Globular Proteins.
Ontology highlight
ABSTRACT:
SUBMITTER: Langenberg T
PROVIDER: S-EPMC7175379 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:The presence of macromolecules in cells geometrically restricts the available space for poplypeptide chains. To study the effects of macromolecular crowding on folding thermodynamics and kinetics, we used an off-lattice model of the all-beta-sheet WW domain in the presence of large spherical particles whose interaction with the polypeptide chain is purely repulsive. At all volume fractions, phi(c), of the crowding agents the stability of the native state is enhanced. Remarkably, the refolding rates, which are larger than the value at phi(c) = 0, increase nonmonotonically as phi(c) increases, reaching a maximum at phi(c)=phi(c)(*). At high values of phi(c), the depletion-induced intramolecular attraction produces compact structures with considerable structure in the denatured state. Changes in native state stability and folding kinetics at phi(c) can be quantitatively mapped onto confinement in a volume-fraction-dependent spherical pore with radius R(s) approximately (4pi/3phi(c))(1/3) R(c) (R(c) is the radius of the crowding particles) as long as phi(c)< or =phi(c)(*). We show that the extent of native state stabilization at finite phi(c) is comparable with that in a spherical pore. In both situations, rate enhancement is due to destabilization of the denatured states with respect to phi(c) = 0.
| S-EPMC555696 | biostudies-literature
Project description:We show that a framework derived from the common character of globular proteins can be used to understand the design of protein sequences, the behavior of intrinsically unstructured proteins, and the formation of amyloid fibrils in a unified manner. Our studies provide compelling support for the idea that protein native-state structures, the structures adopted by intrinsically unstructured proteins on binding as well as those of amyloid aggregates, all reside in a physical state of matter in which the free energy landscape is sculpted not by the specific sequence of amino acids, but rather by considerations of geometry and symmetry. We elucidate the key role played by sequence design in selecting the structure of choice from the predetermined menu of putative native-state structures.
| S-EPMC1458988 | biostudies-literature
Project description:The energetic contributions of hydrogen bonding to protein folding are still unclear, despite more than 70 years of study. This is due partly to the difficulty of extracting thermodynamic information about specific interactions from protein mutagenesis data and partly to the context dependence of hydrogen bond strengths. Herein, we test the hypothesis that hydrogen bond strengths depend on the polarity of their microenvironment, with stronger hydrogen bonds forming in nonpolar surroundings. Double-mutant cycle analysis using a combination of amide-to-ester backbone mutagenesis and traditional side chain mutagenesis revealed that hydrogen bonds can be stronger by up to 1.2 kcal mol(-1) when they are sequestered in hydrophobic surroundings than when they are solvent exposed. Such large coupling energies between hydrogen bond strengths and local polarity suggest that the context dependence of hydrogen bond strengths must be accounted for in any comprehensive account of the forces responsible for protein folding.
| S-EPMC2754385 | biostudies-literature
Project description:Single-molecule manipulation techniques reveal that stretching unravels individually folded domains in the muscle protein titin and the extracellular matrix protein tenascin. These elastic proteins contain tandem repeats of folded domains with beta-sandwich architecture. Herein, we propose by stretching two model sequences (S1 and S2) with four-stranded beta-barrel topology that unfolding forces and pathways in folded domains can be predicted by using only the structure of the native state. Thermal refolding of S1 and S2 in the absence of force proceeds in an all-or-none fashion. In contrast, phase diagrams in the force-temperature (f,T) plane and steered Langevin dynamics studies of these sequences, which differ in the native registry of the strands, show that S1 unfolds in an allor-none fashion, whereas unfolding of S2 occurs via an obligatory intermediate. Force-induced unfolding is determined by the native topology. After proving that the simulation results for S1 and S2 can be calculated by using native topology alone, we predict the order of unfolding events in Ig domain (Ig27) and two fibronectin III type domains ((9)FnIII and (10)FnIII). The calculated unfolding pathways for these proteins, the location of the transition states, and the pulling speed dependence of the unfolding forces reflect the differences in the way the strands are arranged in the native states. We also predict the mechanisms of force-induced unfolding of the coiled-coil spectrin (a three-helix bundle protein) for all 20 structures deposited in the Protein Data Bank. Our approach suggests a natural way to measure the phase diagram in the (f,C) plane, where C is the concentration of denaturants.
| S-EPMC16532 | biostudies-literature
Project description:The issue of changing the structure of globular proteins into an amyloid form is in the focus of researchers' attention. Numerous experimental studies are carried out, and mathematical models to define the essence of amyloid transformation are sought. The present work focuses on the issue of the hydrophobic core structure in amyloids. The form of ordering the hydrophobic core in globular proteins is described by a 3D Gaussian distribution analog to the distribution of hydrophobicity in a spherical micelle. Amyloid fibril is a ribbon-like micelle made up of numerous individual chains, each representing a flat structure. The distribution of hydrophobicity within a single chain included in the fibril describes the 2D Gaussian distribution. Such a description expresses the location of polar residues on a circle with a center with a high level of hydrophobicity. The presence of this type of order in the amyloid forms available in Preotin Data Bank (PDB) (both in proto- and superfibrils) is demonstrated in the present work. In this system, it can be assumed that the amyloid transformation is a chain transition from 3D Gauss ordering to 2D Gauss ordering. This means changing the globular structure to a ribbon-like structure. This observation can provide a simple mathematical model for simulating the amyloid transformation of proteins.
| S-EPMC7370054 | biostudies-literature
Project description:The presence of beta-sheets in the core of amyloid fibrils raised questions as to whether or not beta-sheet-containing proteins, such as transthyretin, are predisposed to form such fibrils. However, we show here that the molecular structure of amyloid fibrils differs more generally from the beta-sheets in native proteins. This difference is evident from the amide I region of the infrared spectrum and relates to the distribution of the phi/psi dihedral angles within the Ramachandran plot, the average number of strands per sheet, and possibly, the beta-sheet twist. These data imply that amyloid fibril formation from native beta-sheet proteins can involve a substantial structural reorganization.
| S-EPMC2287307 | biostudies-literature
Project description:Prediction of the native tertiary structure of a globular protein from the primary sequence will require a potential energy model that can discriminate all nonnative structures from the native structure(s). A successful model must distinguish not only alternate structures that are very nonnative but also alternate structures that are compact and near-native. We describe here a method, based on molecular dynamics simulation, that allows generation of hundreds of compact alternate structures that are arbitrarily close to the native structure. In this way, a significant amount of conformational space in the neighborhood of the native structure can be sampled and these alternate structures can be used as a stringent test of protein folding models. We have used two sets of these alternate structures generated for six crystallographically characterized small globular proteins (1200 alternate structures in all) to test eight empirical energy models for their ability to discriminate alternate from native structures. Seven of the models fail to correctly identify at least some of the alternate structures as nonnative. An atomic solvation model is presented that succeeds in discriminating all 1200 alternate structures from native.
| S-EPMC42689 | biostudies-other
Project description:DisProt is the primary repository of Intrinsically Disordered Proteins (IDPs). This database is manually curated and the annotations there have strong experimental support. Currently, DisProt contains a relatively small number of proteins highlighting the importance of transferring annotations regarding verified disorder state and corresponding functions to homologous proteins in other species. In such a way, providing them with highly valuable information to better understand their biological roles. While the principles and practicalities of homology transfer are well-established for globular proteins, these are largely lacking for disordered proteins. We used DisProt to evaluate the transferability of the annotation terms to orthologous proteins. For each protein, we looked for their orthologs, with the assumption that they will have a similar function. Then, for each protein and their orthologs, we made multiple sequence alignments (MSAs). Disordered sequences are fast evolving and can be hard to align, therefore, we implemented alignment quality control steps ensuring robust alignments before mapping the annotations. We have designed a pipeline to obtain good-quality MSAs and to transfer annotations from any protein to their orthologs. Applying the pipeline to DisProt proteins, from the 1731 entries with 5623 annotations, we can reach 97,555 orthologs and transfer a total of 301,190 terms by homology. We also provide a web server for consulting the results of DisProt proteins and execute the pipeline for any other protein. The server Homology Transfer IDP (HoTIDP) is accessible at http://hotidp.leloir.org.ar.
| S-EPMC10278204 | biostudies-literature
Project description:Ionization-coupled conformational phenomena are ubiquitous in biology. However, quantitative characterization of the underlying thermodynamic cycle comprised of protonation and conformational equilibria has remained an elusive goal. Here we use theory and continuous constant pH molecular dynamics (CpHMD) simulations to provide a thermodynamic description for the coupling of proton titration and conformational exchange between two distinct states of a protein. CpHMD simulations with a hybrid-solvent scheme and the pH-based replica-exchange (REX) protocol are applied to obtain the equilibrium constants and atomic details of the ionization-coupled conformational exchange between open and closed states of an engineered mutant of staphylococcal nuclease. Although the coupling of protonation and conformational equilibria is not exact in the simulation, the results are encouraging. They demonstrate that REX-CpHMD simulations can be used to study thermodynamics of ionization-coupled conformational processes--which has not possible using present experimental techniques or traditional simulations based on fixed protonation states.
| S-EPMC3318114 | biostudies-literature
Project description:An amino acid sequence, in the context of the solvent environment, contains all of the thermodynamic information necessary to encode a three-dimensional protein structure. To investigate the relationship between an amino acid sequence and its corresponding protein fold, a database of thermodynamic stability information was assembled that spanned 2951 residues from 44 nonhomologous proteins. This information was obtained using the COREX algorithm, which computes an ensemble-based description of the native state of a protein. It was observed that amino acid types partitioned unequally into high, medium, and low thermodynamic stability environments. Furthermore, these distributions were reproducible and were significantly different than those expected from random partitioning. To assess the structural importance of the distributions, simple fold-recognition experiments were performed based on a 3D-1D scoring matrix containing only COREX residue stability information. This procedure was able to recover amino acid sequences corresponding to correct target structures more effectively than scoring matrices derived from randomized data. High-scoring sequences were often aligned correctly with their corresponding target profiles, suggesting that calculated thermodynamic stability profiles have the potential to encode sequence information. As a control, identical fold-recognition experiments were performed on the same database of proteins using DSSP secondary structure information in the scoring matrix, instead of COREX residue stability information. The comparable performance of both approaches suggested that COREX residue stability information and secondary structure information could be of equivalent utility in more sophisticated fold-recognition techniques. The results of this work are a consequence of the idea that amino acid sequences fold not into single, rigidly stable structures but rather into thermodynamic ensembles best represented by a time-averaged structure.
| S-EPMC2374190 | biostudies-literature