Extended surfaces modulate and can catalyze hydrophobic effects

Interfaces are a most common motif in complex systems. To understand how the presence of interfaces affect hydrophobic phenomena, we use molecular simulations and theory to study hydration of solutes at interfaces. The solutes range in size from sub-nanometer to a few nanometers. The interfaces are self-assembled monolayers with a range of chemistries, from hydrophilic to hydrophobic. We show that the driving force for assembly in the vicinity of a hydrophobic surface is weaker than that in bulk water, and decreases with increasing temperature, in contrast to that in the bulk. We explain these distinct features in terms of an interplay between interfacial fluctuations and excluded volume effects---the physics encoded in Lum-Chandler-Weeks theory [J. Phys. Chem. B 103, 4570--4577 (1999)]. Our results suggest a catalytic role for hydrophobic interfaces in the unfolding of proteins, for example, in the interior of chaperonins and in amyloid formation.Comment: 22 pages, 5 figure

Model for the hydration of non-polar compounds and polymers

We introduce an exactly solvable statistical-mechanical model of the hydration of non-polar compounds, based on grouping water molecules in clusters where hydrogen bonds and isotropic interactions occur; interactions between clusters are neglected. Analytical results show that an effective strengthening of hydrogen bonds in the presence of the solute, together with a geometric reorganization of water molecules, are enough to yield hydrophobic behavior. We extend our model to describe a non-polar homopolymer in aqueous solution, obtaining a clear evidence of both ``cold'' and ``warm'' swelling transitions. This suggests that our model could be relevant to describe some features of protein folding.Comment: REVTeX, 6 pages, 3 figure

CAPS-DB: a structural classification of helix-capping motifs

The regions of the polypeptide chain immediately preceding or following an α-helix are known as Nt- and Ct cappings, respectively. Cappings play a central role stabilizing α-helices due to lack of intrahelical hydrogen bonds in the first and last turn. Sequence patterns of amino acid type preferences have been derived for cappings but the structural motifs associated to them are still unclassified. CAPS-DB is a database of clusters of structural patterns of different capping types. The clustering algorithm is based in the geometry and the (ϕ–ψ)-space conformation of these regions. CAPS-DB is a relational database that allows the user to search, browse, inspect and retrieve structural data associated to cappings. The contents of CAPS-DB might be of interest to a wide range of scientist covering different areas such as protein design and engineering, structural biology and bioinformatics. The database is accessible at: http://www.bioinsilico.org/CAPSDB

The correlation of RNase A enzymatic activity with the changes in the distance between Nepsilon2-His12 and N delta1-His119 upon addition of stabilizing and destabilizing salts.

The effect of stabilizing and destabilizing salts on the catalytic behavior of ribonuclease A (RNase A) was investigated at pH 7.5 and 25 degrees C, using spectrophotometric, viscometric and molecular dynamic methods. The changes in the distance between N(epsilon2) of His(12) and N(delta1) of His(119) at the catalytic center of RNase A upon the addition of sodium sulfate, sodium hydrogen sulfate and sodium thiocyanate were evaluated by molecular dynamic methods. The compactness and expansion in terms of Stokes radius of RNase A upon the addition of sulfate ions as kosmotropic salts, and thiocyanate ion as a chaotropic salt, were estimated by viscometric measurements. Enzyme activity was measured using cytidine 2', 3'-cyclic monophosphate as a substrate. The results from the measurements of distances between N(epsilon2) of His(12) and N(delta1) of His(119) and Stokes radius suggest (i) that the presence of sulfate ions decreases the distance between the catalytic His residues and increases the globular compactness, and (ii) that there is an expansion of the enzyme surface as well as elongation of the catalytic center in the presence of thiocyanate ion. These findings are in agreement with activity measurements

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

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

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