Computing the relative stabilities and the per-residue components in protein conformational changes.
Ontology highlight
ABSTRACT:
SUBMITTER: Roy A
PROVIDER: S-EPMC3905753 | biostudies-literature | 2014 Jan
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:A new method termed "Relative Principal Components Analysis" (RPCA) is introduced that extracts optimal relevant principal components to describe the change between two data samples representing two macroscopic states. The method is widely applicable in data-driven science. Calculating the components is based on a physical framework that introduces the objective function (the Kullback-Leibler divergence) appropriate for quantifying the change of the macroscopic state affected by the changes in the microscopic features. To demonstrate the applicability of RPCA, we analyze the thermodynamically relevant conformational changes of the protein HIV-1 protease upon binding to different drug molecules. In this case, the RPCA method provides a sound thermodynamic foundation for analyzing the binding process and thus characterizing both the collective and the locally relevant conformational changes. Moreover, the relevant collective conformational changes can be reconstructed from the informative latent variables to exhibit both the enhanced and the restricted conformational fluctuations upon ligand association.
| S-EPMC6728065 | biostudies-literature
Project description:New amino acid sequences of proteins are being learned at a rapid rate, thanks to modern genomics. The native structures and functions of those proteins can often be inferred using bioinformatics methods. We show here that it is also possible to infer the stabilities and thermal folding properties of proteins, given only simple genomics information: the chain length and the numbers of charged side chains. In particular, our model predicts DeltaH(T), DeltaS(T), DeltaC(p), and DeltaF(T)--the folding enthalpy, entropy, heat capacity, and free energy--as functions of temperature T; the denaturant m values in guanidine and urea; the pH-temperature-salt phase diagrams, and the energy of confinement F(s) of the protein inside a cavity of radius s. All combinations of these phase equilibria can also then be computed from that information. As one illustration, we compute the pH and salt conditions that would denature a protein inside a small confined cavity. Because the model is analytical, it is computationally efficient enough that it could be used to automatically annotate whole proteomes with protein stability information.
| S-EPMC2705543 | biostudies-literature
Project description:Intrinsically disordered proteins (IDPs) adopt heterogeneous ensembles of conformations under physiological conditions. Understanding the relationship between amino acid sequence and conformational ensembles of IDPs can help clarify the role of disorder in physiological function. Recent studies revealed that polar IDPs favor collapsed ensembles in water despite the absence of hydrophobic groups--a result that holds for polypeptide backbones as well. By studying highly charged polypeptides, a different archetype of IDPs, we assess how charge content modulates the intrinsic preference of polypeptide backbones for collapsed structures. We characterized conformational ensembles for a set of protamines in aqueous milieus using molecular simulations and fluorescence measurements. Protamines are arginine-rich IDPs involved in the condensation of chromatin during spermatogenesis. Simulations based on the ABSINTH implicit solvation model predict the existence of a globule-to-coil transition, with net charge per residue serving as the discriminating order parameter. The transition is supported by quantitative agreement between simulation and experiment. Local conformational preferences partially explain the observed trends of polymeric properties. Our results lead to the proposal of a schematic protein phase diagram that should enable prediction of polymeric attributes for IDP conformational ensembles using easily calculated physicochemical properties of amino acid sequences. Although sequence composition allows the prediction of polymeric properties, interresidue contact preferences of protamines with similar polymeric attributes suggest that certain details of conformational ensembles depend on the sequence. This provides a plausible mechanism for specificity in the functions of IDPs.
| S-EPMC2889596 | biostudies-literature
Project description:Protein interactions are often accompanied by significant changes in conformation. We have analyzed the relationships between protein structures and the conformational changes they undergo upon binding. Based upon this, we introduce a simple measure, the relative solvent accessible surface area, which can be used to predict the magnitude of binding-induced conformational changes from the structures of either monomeric proteins or bound subunits. Applying this to a large set of protein complexes suggests that large conformational changes upon binding are common. In addition, we observe considerable enrichment of intrinsically disordered sequences in proteins predicted to undergo large conformational changes. Finally, we demonstrate that the relative solvent accessible surface area of monomeric proteins can be used as a simple proxy for protein flexibility. This reveals a powerful connection between the flexibility of unbound proteins and their binding-induced conformational changes, consistent with the conformational selection model of molecular recognition.
| S-EPMC3145976 | biostudies-literature
Project description:The prediction of changes in protein stability and structure resulting from single amino acid substitutions is both a fundamental test of macromolecular modeling methodology and an important current problem as high throughput sequencing reveals sequence polymorphisms at an increasing rate. In principle, given the structure of a wild-type protein and a point mutation whose effects are to be predicted, an accurate method should recapitulate both the structural changes and the change in the folding-free energy. Here, we explore the performance of protocols which sample an increasing diversity of conformations. We find that surprisingly similar performances in predicting changes in stability are achieved using protocols that involve very different amounts of conformational sampling, provided that the resolution of the force field is matched to the resolution of the sampling method. Methods involving backbone sampling can in some cases closely recapitulate the structural changes accompanying mutations but not surprisingly tend to do more harm than good in cases where structural changes are negligible. Analysis of the outliers in the stability change calculations suggests areas needing particular improvement; these include the balance between desolvation and the formation of favorable buried polar interactions, and unfolded state modeling.
| S-EPMC3760476 | biostudies-literature
Project description:We have previously studied the unfolding equilibrium of bacterioopsin in a single phase solvent, using Förster mechanism fluorescence resonance energy transfer (FRET) as a probe, from tryptophan donors to a dansyl acceptor. We observed an apparent unfolding transition in bacterioopsin perturbed by increasing ethanol concentrations [Nannepaga et al. (2004) Biochemistry 43, 50-59]. We have further investigated this transition and find that the unfolding is pH-dependent. We have now measured the apparent pK of acid-induced unfolding of bacterioopsin in 90% ethanol. When the acceptor is on helix B (Lys 41), the apparent pK for unfolding is 4.75; on the EF connecting loop (Cys 163), 5.15; and on helix G (Cys 222), 5.75. Five-helix proteolytic fragments are less stable. The apparent unfolding pKs are 5.46 for residues 72-248 (Cys 163) and 7.36 for residues 1-166 (Lys 41). When interpreted in terms of a simple equilibrium model for unfolding, the apparent pKs give relative free energies of unfolding in the range of -0.54 to -3.5 kcal/mol. The results suggest that the C-terminal helix of bacterioopsin is less stably folded than the N-terminal helices. We analyzed the pairwise helix-helix interaction surfaces of bacteriorhodopsin and three other seven-transmembrane-helix proteins on the basis of crystal structures. The results show that the interaction surfaces are smoother and the helix axis separations are closer in the amino-terminal two-thirds of the proteins compared with the carboxyl-terminal one-third. However, the F helix is important in stabilizing the folded structure, as shown by the instability of the 1-166 fragment. Considering the high-resolution crystal structure of bacteriorhodopsin, there are no obvious helix-helix interactions involving protein side chains which would be destabilized by protonation at the estimated pH of the unfolding transitions. However, a number of helix-bridging water molecules could become protonated, thereby weakening the helix-helix interactions.
| S-EPMC2519877 | biostudies-literature
Project description:The analysis of evolutionary amino acid correlations has recently attracted a surge of renewed interest, also due to their successful use in de novo protein native structure prediction. However, many aspects of protein function, such as substrate binding and product release in enzymatic activity, can be fully understood only in terms of an equilibrium ensemble of alternative structures, rather than a single static structure. In this paper we combine coevolutionary data and molecular dynamics simulations to study protein conformational heterogeneity. To that end, we adapt the Boltzmann-learning algorithm to the analysis of homologous protein sequences and develop a coarse-grained protein model specifically tailored to convert the resulting contact predictions to a protein structural ensemble. By means of exhaustive sampling simulations, we analyze the set of conformations that are consistent with the observed residue correlations for a set of representative protein domains, showing that (i) the most representative structure is consistent with the experimental fold and (ii) the various regions of the sequence display different stability, related to multiple biologically relevant conformations and to the cooperativity of the coevolving pairs. Moreover, we show that the proposed protocol is able to reproduce the essential features of a protein folding mechanism as well as to account for regions involved in conformational transitions through the correct sampling of the involved conformers.
| S-EPMC4640757 | biostudies-literature
Project description:Alexander disease (AxD) is an often fatal astrogliopathy caused by dominant gain-of-function missense mutations in the glial fibrillary acidic protein (GFAP) gene. The mechanism by which the mutations produce the AxD phenotype is not known. However, the observation that features of AxD are displayed by mice that express elevated levels of GFAP from a human WT GFAP transgene has contributed to the notion that the mutations produce AxD by increasing accumulation of total GFAP above some toxic threshold rather than the mutant GFAP being inherently toxic. A possible mechanism for accumulation of GFAP in AxD patients is that the mutated GFAP variants are more stable than the WT, an attribution abetted by observations that GFAP complexes containing GFAP variants are more resistant to solvent extraction. Here we tested this hypothesis by determining the relative levels of WT and mutant GFAP in three individuals with AxD, each of whom carried a common but different GFAP mutation (R79C, R239H, or R416W). Mass spectrometry analysis identified a peptide specific to the mutant or WT GFAP in each patient, and we quantified this peptide by comparing its signal to that of an added [15N]GFAP standard. In all three individuals, the level of mutant GFAP was less than that of the WT. This finding suggests that AxD onset is due to an intrinsic toxicity of the mutant GFAP instead of it acting indirectly by being more stable than WT GFAP and thereby increasing the total GFAP level.
| S-EPMC6816090 | biostudies-literature
Project description:In this paper, a novel method to compute side chain conformational variations for a protein molecule tunnel (or channel) is proposed. From the conformational variations, we compute the flexibly deformed shapes of the initial tunnel, and present a way to compute the maximum size of the ligand that can pass through the deformed tunnel. By using the two types of graphs corresponding to amino acids and their side chain rotamers, the suggested algorithm classifies amino acids and rotamers which possibly have collisions. Based on the divide and conquer technique, local side chain conformations are computed first, and then a global conformation is generated by combining them. With the exception of certain cases, experimental results show that the algorithm finds up to 327,680 valid side chain conformations from 12⁸~1233 conformation candidates within three seconds.
| S-EPMC6222877 | biostudies-literature
Project description:Conformational changes upon protein-protein association are the key element of the binding mechanism. The study presents a systematic large-scale analysis of such conformational changes in the side chains. The results indicate that short and long side chains have different propensities for the conformational changes. Long side chains with three or more dihedral angles are often subject to large conformational transition. Shorter residues with one or two dihedral angles typically undergo local conformational changes not leading to a conformational transition. A relationship between the local readjustments and the equilibrium fluctuations of a side chain around its unbound conformation is suggested. Most of the side chains undergo larger changes in the dihedral angle most distant from the backbone. The frequencies of the core-to-surface interface transitions of six nonpolar residues and Tyr are larger than the frequencies of the opposite surface-to-core transitions. The binding increases both polar and nonpolar interface areas. However, the increase of the nonpolar area is larger for all considered classes of protein complexes, suggesting that the protein association perturbs the unbound interfaces to increase the hydrophobic contribution to the binding free energy. To test modeling approaches to side-chain flexibility in protein docking, conformational changes in the X-ray set were compared with those in the docking decoy sets. The results lead to a better understanding of the conformational changes in proteins and suggest directions for efficient conformational sampling in docking protocols.
| S-EPMC3073032 | biostudies-literature