Non-linearity of the collagen triple helix in solution and implications for collagen function

Collagen adopts a characteristic supercoiled triple helical conformation which requires a repeating (Xaa-Yaa-Gly)n sequence. Despite the abundance of collagen, a combined experimental and atomistic modelling approach has not so far quantitated the degree of flexibility seen experimentally in the solution structures of collagen triple helices. To address this question, we report an experimental study on the flexibility of varying lengths of collagen triple helical peptides, composed of six, eight, ten and twelve repeats of the most stable Pro-Hyp-Gly (POG) units. In addition, one unblocked peptide, (POG)10unblocked, was compared with the blocked (POG)10 as a control for the significance of end effects. Complementary analytical ultracentrifugation and synchrotron small angle X-ray scattering data showed that the conformations of the longer triple helical peptides were not well explained by a linear structure derived from crystallography. To interpret these data, molecular dynamics simulations were used to generate 50 000 physically realistic collagen structures for each of the helices. These structures were fitted against their respective scattering data to reveal the best fitting structures from this large ensemble of possible helix structures. This curve fitting confirmed a small degree of non-linearity to exist in these best fit triple helices, with the degree of bending approximated as 4–17° from linearity. Our results open the way for further studies of other collagen triple helices with different sequences and stabilities in order to clarify the role of molecular rigidity and flexibility in collagen extracellular and immune function and disease

A Folding Pathway-Dependent Score to Recognize Membrane Proteins

While various approaches exist to study protein localization, it is still a challenge to predict where proteins localize. Here, we consider a mechanistic viewpoint for membrane localization. Taking into account the steps for the folding pathway of α-helical membrane proteins and relating biophysical parameters to each of these steps, we create a score capable of predicting the propensity for membrane localization and call it FP3mem. This score is driven from the principal component analysis (PCA) of the biophysical parameters related to membrane localization. FP3mem allows us to rationalize the colocalization of a number of channel proteins with the Cav1.2 channel by their fewer propensities for membrane localization

Structure and Behavior of Human α-Thrombin upon Ligand Recognition: Thermodynamic and Molecular Dynamics Studies

Thrombin is a serine proteinase that plays a fundamental role in coagulation. In this study, we address the effects of ligand site recognition by alpha-thrombin on conformation and energetics in solution. Active site occupation induces large changes in secondary structure content in thrombin as shown by circular dichroism. Thrombin-D-Phe-Pro-Arg-chloromethyl ketone (PPACK) exhibits enhanced equilibrium and kinetic stability compared to free thrombin, whose difference is rooted in the unfolding step. Small-angle X-ray scattering (SAXS) measurements in solution reveal an overall similarity in the molecular envelope of thrombin and thrombin-PPACK, which differs from the crystal structure of thrombin. Molecular dynamics simulations performed with thrombin lead to different conformations than the one observed in the crystal structure. These data shed light on the diversity of thrombin conformers not previously observed in crystal structures with distinguished catalytic and conformational behaviors, which might have direct implications on novel strategies to design direct thrombin inhibitors

Universality of Thermodynamic Constants Governing Biological Growth Rates

Background: Mathematical models exist that quantify the effect of temperature on poikilotherm growth rate. One family of such models assumes a single rate-limiting ‘master reaction ’ using terms describing the temperature-dependent denaturation of the reaction’s enzyme. We consider whether such a model can describe growth in each domain of life. Methodology/Principal Findings: A new model based on this assumption and using a hierarchical Bayesian approach fits simultaneously 95 data sets for temperature-related growth rates of diverse microorganisms from all three domains of life, Bacteria, Archaea and Eukarya. Remarkably, the model produces credible estimates of fundamental thermodynamic parameters describing protein thermal stability predicted over 20 years ago. Conclusions/Significance: The analysis lends support to the concept of universal thermodynamic limits to microbial growth rate dictated by protein thermal stability that in turn govern biological rates. This suggests that the thermal stability of proteins is a unifying property in the evolution and adaptation of life on earth. The fundamental nature of this conclusion has importance for many fields of study including microbiology, protein chemistry, thermal biology, and ecological theory including, for example, the influence of the vast microbial biomass and activity in the biosphere that is poorly described in current climate models

New Binding Mode to TNF-Alpha Revealed by Ubiquitin-Based Artificial Binding Protein

A variety of approaches have been employed to generate binding proteins from non-antibody scaffolds. Utilizing a beta-sheet of the human ubiquitin for paratope creation we obtained binding proteins against tumor necrosis factor (TNF)-alpha. The bioactive form of this validated pharmacological target protein is a non-covalently linked homo-trimer. This structural feature leads to the observation of a certain heterogeneity concerning the binding mode of TNF-alpha binding molecules, for instance in terms of monomer/trimer specificity. We analyzed a ubiquitin-based TNF-alpha binder, selected by ribosome display, with a particular focus on its mode of interaction. Using enzyme-linked immunosorbent assays, specific binding to TNF-alpha with nanomolar affinity was observed. In isothermal titration calorimetry we obtained comparable results regarding the affinity and detected an exothermic reaction with one ubiquitin-derived binding molecule binding one TNF-alpha trimer. Using NMR spectroscopy and other analytical methods the 1∶3 stoichiometry could be confirmed. Detailed binding analysis showed that the interaction is affected by the detergent Tween-20. Previously, this phenomenon was reported only for one other type of alternative scaffold-derived binding proteins – designed ankyrin repeat proteins – without further investigation. As demonstrated by size exclusion chromatography and NMR spectroscopy, the presence of the detergent increases the association rate significantly. Since the special architecture of TNF-alpha is known to be modulated by detergents, the access to the recognized epitope is indicated to be restricted by conformational transitions within the target protein. Our results suggest that the ubiquitin-derived binding protein targets a new epitope on TNF-alpha, which differs from the epitopes recognized by TNF-alpha neutralizing antibodies

Rational Design of Protein Stability: Effect of (2S,4R)-4-Fluoroproline on the Stability and Folding Pathway of Ubiquitin

BACKGROUND: Many strategies have been employed to increase the conformational stability of proteins. The use of 4-substituted proline analogs capable to induce pre-organization in target proteins is an attractive tool to deliver an additional conformational stability without perturbing the overall protein structure. Both, peptides and proteins containing 4-fluorinated proline derivatives can be stabilized by forcing the pyrrolidine ring in its favored puckering conformation. The fluorinated pyrrolidine rings of proline can preferably stabilize either a C(γ)-exo or a C(γ)-endo ring pucker in dependence of proline chirality (4R/4S) in a complex protein structure. To examine whether this rational strategy can be generally used for protein stabilization, we have chosen human ubiquitin as a model protein which contains three proline residues displaying C(γ)-exo puckering. METHODOLOGY/PRINCIPAL FINDINGS: While (2S,4R)-4-fluoroproline ((4R)-FPro) containing ubiquitinin can be expressed in related auxotrophic Escherichia coli strain, all attempts to incorporate (2S,4S)-4-fluoroproline ((4S)-FPro) failed. Our results indicate that (4R)-FPro is favoring the C(γ)-exo conformation present in the wild type structure and stabilizes the protein structure due to a pre-organization effect. This was confirmed by thermal and guanidinium chloride-induced denaturation profile analyses, where we observed an increase in stability of -4.71 kJ·mol(-1) in the case of (4R)-FPro containing ubiquitin ((4R)-FPro-ub) compared to wild type ubiquitin (wt-ub). Expectedly, activity assays revealed that (4R)-FPro-ub retained the full biological activity compared to wt-ub. CONCLUSIONS/SIGNIFICANCE: The results fully confirm the general applicability of incorporating fluoroproline derivatives for improving protein stability. In general, a rational design strategy that enforces the natural occurring proline puckering conformation can be used to stabilize the desired target protein

A Didactic Model of Macromolecular Crowding Effects on Protein Folding

A didactic model is presented to illustrate how the effect of macromolecular crowding on protein folding and association is modeled using current analytical theory and discrete molecular dynamics. While analytical treatments of crowding may consider the effect as a potential of average force acting to compress a polypeptide chain into a compact state, the use of simulations enables the presence of crowding reagents to be treated explicitly. Using an analytically solvable toy model for protein folding, an approximate statistical thermodynamic method is directly compared to simulation in order to gauge the effectiveness of current analytical crowding descriptions. Both methodologies are in quantitative agreement under most conditions, indication that both current theory and simulation methods are capable of recapitulating aspects of protein folding even by utilizing a simplistic protein model

Hydrogen bond networks determine emergent mechanical and thermodynamic properties across a protein family

<p>Abstract</p> <p>Background</p> <p>Gram-negative bacteria use periplasmic-binding proteins (bPBP) to transport nutrients through the periplasm. Despite immense diversity within the recognized substrates, all members of the family share a common fold that includes two domains that are separated by a conserved hinge. The hinge allows the protein to cycle between open (apo) and closed (ligated) conformations. Conformational changes within the proteins depend on a complex interplay of mechanical and thermodynamic response, which is manifested as an increase in thermal stability and decrease of flexibility upon ligand binding.</p> <p>Results</p> <p>We use a distance constraint model (DCM) to quantify the give and take between thermodynamic stability and mechanical flexibility across the bPBP family. Quantitative stability/flexibility relationships (QSFR) are readily evaluated because the DCM links mechanical and thermodynamic properties. We have previously demonstrated that QSFR is moderately conserved across a mesophilic/thermophilic RNase H pair, whereas the observed variance indicated that different enthalpy-entropy mechanisms allow similar mechanical response at their respective melting temperatures. Our predictions of heat capacity and free energy show marked diversity across the bPBP family. While backbone flexibility metrics are mostly conserved, cooperativity correlation (long-range couplings) also demonstrate considerable amount of variation. Upon ligand removal, heat capacity, melting point, and mechanical rigidity are, as expected, lowered. Nevertheless, significant differences are found in molecular cooperativity correlations that can be explained by the detailed nature of the hydrogen bond network.</p> <p>Conclusion</p> <p>Non-trivial mechanical and thermodynamic variation across the family is explained by differences within the underlying H-bond networks. The mechanism is simple; variation within the H-bond networks result in altered mechanical linkage properties that directly affect intrinsic flexibility. Moreover, varying numbers of H-bonds and their strengths control the likelihood for energetic fluctuations as H-bonds break and reform, thus directly affecting thermodynamic properties. Consequently, these results demonstrate how unexpected large differences, especially within cooperativity correlation, emerge from subtle differences within the underlying H-bond network. This inference is consistent with well-known results that show allosteric response within a family generally varies significantly. Identifying the hydrogen bond network as a critical determining factor for these large variances may lead to new methods that can predict such effects.</p

Stabilizing Salt-Bridge Enhances Protein Thermostability by Reducing the Heat Capacity Change of Unfolding

Most thermophilic proteins tend to have more salt bridges, and achieve higher thermostability by up-shifting and broadening their protein stability curves. While the stabilizing effect of salt-bridge has been extensively studied, experimental data on how salt-bridge influences protein stability curves are scarce. Here, we used double mutant cycles to determine the temperature-dependency of the pair-wise interaction energy and the contribution of salt-bridges to ΔCp in a thermophilic ribosomal protein L30e. Our results showed that the pair-wise interaction energies for the salt-bridges E6/R92 and E62/K46 were stabilizing and insensitive to temperature changes from 298 to 348 K. On the other hand, the pair-wise interaction energies between the control long-range ion-pair of E90/R92 were negligible. The ΔCp of all single and double mutants were determined by Gibbs-Helmholtz and Kirchhoff analyses. We showed that the two stabilizing salt-bridges contributed to a reduction of ΔCp by 0.8–1.0 kJ mol−1 K−1. Taken together, our results suggest that the extra salt-bridges found in thermophilic proteins enhance the thermostability of proteins by reducing ΔCp, leading to the up-shifting and broadening of the protein stability curves

Polar or Apolar—The Role of Polarity for Urea-Induced Protein Denaturation

Urea-induced protein denaturation is widely used to study protein folding and stability; however, the molecular mechanism and driving forces of this process are not yet fully understood. In particular, it is unclear whether either hydrophobic or polar interactions between urea molecules and residues at the protein surface drive denaturation. To address this question, here, many molecular dynamics simulations totalling ca. 7 µs of the CI2 protein in aqueous solution served to perform a computational thought experiment, in which we varied the polarity of urea. For apolar driving forces, hypopolar urea should show increased denaturation power; for polar driving forces, hyperpolar urea should be the stronger denaturant. Indeed, protein unfolding was observed in all simulations with decreased urea polarity. Hyperpolar urea, in contrast, turned out to stabilize the native state. Moreover, the differential interaction preferences between urea and the 20 amino acids turned out to be enhanced for hypopolar urea and suppressed (or even inverted) for hyperpolar urea. These results strongly suggest that apolar urea–protein interactions, and not polar interactions, are the dominant driving force for denaturation. Further, the observed interactions provide a detailed picture of the underlying molecular driving forces. Our simulations finally allowed characterization of CI2 unfolding pathways. Unfolding proceeds sequentially with alternating loss of secondary or tertiary structure. After the transition state, unfolding pathways show large structural heterogeneity