Stability of protein structure and hydrophobic interaction

Temperature stability of proteins: Analysis of irreversible denaturation using isothermal calorimetry

The structural stability of proteins has been traditionally studied under conditions in which the folding/unfolding reaction is reversible, since thermodynamic parameters can only be determined under these conditions. Achieving reversibility conditions in temperature stability experiments has often required performing the experiments at acidic pH or other nonphysiological solvent conditions. With the rapid development of protein drugs, the fastest growing segment in the pharmaceutical industry, the need to evaluate protein stability under formulation conditions has acquired renewed urgency. Under formulation conditions and the required high protein concentration (∼100 mg/mL), protein denaturation is irreversible and frequently coupled to aggregation and precipitation. In this article, we examine the thermal denaturation of hen egg white lysozyme (HEWL) under irreversible conditions and concentrations up to 100 mg/mL using several techniques, especially isothermal calorimetry which has been used to measure the enthalpy and kinetics of the unfolding and aggregation/precipitation at 12°C below the transition temperature measured by DSC. At those temperatures the rate of irreversible protein denaturation and aggregation of HEWL is measured to be on the order of 1 day^-1 . Isothermal calorimetry appears a suitable technique to identify buffer formulation conditions that maximize the long term stability of protein drugs.

Role of water in the formation of macromolecular structures

This review shows that water in biological systems is not just a passive liquid solvent but also a partner in the formation of the structure of proteins, nucleic acids and their complexes, thereby contributing to the stability and flexibility required for their proper function. Reciprocally, biological macromolecules affect the state of the water contacting them, so that it is only partly in the normal liquid state, being somewhat ordered when bound to macromolecules. While the compaction of globular proteins results from the reluctance of their hydrophobic groups to interact with water, the collagen superhelix is maintained by water forming a hydroxyproline-controlled frame around this coiled-coil macromolecule. As for DNA, its stability and rigidity are linked to water fixed by AT pairs in the minor groove: this leads to the enthalpic contribution of AT pairs exceeding that of GC pairs, but this is overbalanced by their greater entropy contribution, with the result that AT pairs melt at lower temperatures than GCs. Loss of this water drives transcription factor binding to the minor groove.

Multifactorial level of extremostability of proteins: can they be exploited for protein engineering?

Research on extremostable proteins has seen immense growth in the past decade owing to their industrial importance. Basic research of attributes related to extreme-stability requires further exploration. Modern mechanistic approaches to engineer such proteins in vitro will have more impact in industrial biotechnology economy. Developing a priori knowledge about the mechanism behind extreme-stability will nurture better understanding of pathways leading to protein molecular evolution and folding. This review is a vivid compilation about all classes of extremostable proteins and the attributes that lead to myriad of adaptations divulged after an extensive study of 6495 articles belonging to extremostable proteins. Along with detailing on the rationale behind extreme-stability of proteins, emphasis has been put on modern approaches that have been utilized to render proteins extremostable by protein engineering. It was understood that each protein shows different approaches to extreme-stability governed by minute differences in their biophysical properties and the milieu in which they exist. Any general rule has not yet been drawn regarding adaptive mechanisms in extreme environments. This review was further instrumental to understand the drawback of the available 14 stabilizing mutation prediction algorithms. Thus, this review lays the foundation to further explore the biophysical pleiotropy of extreme-stable proteins to deduce a global prediction model for predicting the effect of mutations on protein stability.

Magnetic nanoparticles as double-edged swords: concentration-dependent ordering or disordering effects on lysozyme

The study of the interaction of nanoparticles (NPs) with proteins is of great importance due to its relevance in several fields including nano-biosafety, nano-bioscience, nano-biomedicine, and nano-biotechnology.

Effects of ligand binding on the stability of aldo-keto reductases: Implications for stabilizer or destabilizer chaperones

Ligands such as enzyme inhibitors stabilize the native conformation of a protein upon binding to the native state, but some compounds destabilize the native conformation upon binding to the non-native state. The former ligands are termed "stabilizer chaperones" and the latter ones "destabilizer chaperones." Because the stabilization effects are essential for the medical chaperone (MC) hypothesis, here we have formulated a thermodynamic system consisting of a ligand and a protein in its native- and non-native state. Using the differential scanning fluorimetry and the circular dichroism varying the urea concentration and temperature, we found that when the coenzyme NADP+ was absent, inhibitors such as isolithocholic acid stabilized the aldo-keto reductase AKR1A1 upon binding, which showed actually the three-state folding, but destabilized AKR1B10. In contrast, in the presence of NADP+ , they destabilized AKR1A1 and stabilized AKR1B10. To explain these phenomena, we decomposed the free energy of stabilization (ΔΔG) into its enthalpy (ΔΔH) and entropy (ΔΔS) components. Then we found that in a relatively unstable protein showing the three-state folding, native conformation was stabilized by the negative ΔΔH in association with the negative ΔΔS, suggesting that the stabilizer chaperon decreases the conformational fluctuation of the target protein or increase its hydration. However, in other cases, ΔΔG was essentially determined by the delicate balance between ΔΔH and ΔΔS. The proposed thermodynamic formalism is applicable to the system including multiple ligands with allosteric interactions. These findings would promote the development of screening strategies for MCs to regulate the target conformations.

How cooperative are protein folding and unfolding transitions?

A thermodynamically and kinetically simple picture of protein folding envisages only two states, native (N) and unfolded (U), separated by a single activation free energy barrier, and interconverting by cooperative two-state transitions. The folding/unfolding transitions of many proteins occur, however, in multiple discrete steps associated with the formation of intermediates, which is indicative of reduced cooperativity. Furthermore, much advancement in experimental and computational approaches has demonstrated entirely non-cooperative (gradual) transitions via a continuum of states and a multitude of small energetic barriers between the N and U states of some proteins. These findings have been instrumental towards providing a structural rationale for cooperative versus noncooperative transitions, based on the coupling between interaction networks in proteins. The cooperativity inherent in a folding/unfolding reaction appears to be context dependent, and can be tuned via experimental conditions which change the stabilities of N and U. The evolution of cooperativity in protein folding transitions is linked closely to the evolution of function as well as the aggregation propensity of the protein. A large activation energy barrier in a fully cooperative transition can provide the kinetic control required to prevent the accumulation of partially unfolded forms, which may promote aggregation. Nevertheless, increasing evidence for barrier-less "downhill" folding, as well as for continuous "uphill" unfolding transitions, indicate that gradual non-cooperative processes may be ubiquitous features on the free energy landscape of protein folding.

Equation of State for Phospholipid Self-Assembly

Phospholipid self-assembly is the basis of biomembrane stability. The entropy of transfer from water to self-assembled micelles of lysophosphatidylcholines and diacyl phosphatidylcholines with different chain lengths converges to a common value at a temperature of 44°C. The corresponding enthalpies of transfer converge at ∼-18°C. An equation of state for the free energy of self-assembly formulated from this thermodynamic data depends on the heat capacity of transfer as the sole parameter needed to specify a particular lipid. For lipids lacking calorimetric data, measurement of the critical micelle concentration at a single temperature suffices to define an effective heat capacity according to the model. Agreement with the experimental temperature dependence of the critical micelle concentration is then good. The predictive powers should extend also to amphiphile partitioning and the kinetics of lipid-monomer transfer.

Thermal protein unfolding by differential scanning calorimetry and circular dichroism spectroscopy Two-state model versus sequential unfolding

Thermally-induced protein unfolding is commonly described with the two-state model. This model assumes only two types of protein molecules in solution, the native (N) and the denatured, unfolded (U) protein. In reality, protein unfolding is a multistep process, even if intermediate states are only sparsely populated. As an alternative approach we explore the Zimm-Bragg theory, originally developed for the α-helix-to-random coil transition of synthetic polypeptides. The theory includes intermediate structures with concentrations determined by the cooperativity of the unfolding reaction. We illustrate the differences between the two-state model and the Zimm-Bragg theory with measurements of apolipoprotein A-1 and lysozyme by differential scanning calorimetry (DSC) and CD spectroscopy. Nine further protein examples are taken from the literature. The Zimm-Bragg theory provides a perfect fit of the calorimetric unfolding transitions for all proteins investigated. In contrast, the transition curves and enthalpies predicted by the two-state model differ considerably from the experimental results. Apolipoprotein A-1 is ~50% α-helical at ambient temperature and its unfolding follows the classical α-helix-to-random coil equilibrium. The unfolding of proteins with little α-helix content, such as lysozyme, can also be analyzed with the Zimm-Bragg theory by introducing the concept of 'folded' and 'unfolded' peptide units assuming an average unfolding enthalpy per peptide unit. DSC is the method of choice to measure the unfolding enthalpy, , but CD spectroscopy in combination with the two-state model is often used to deduce the unfolding enthalpy. This can lead to erroneous result. Not only are different enthalpies required to describe the CD and DSC transition curves but these values deviate distinctly from the experimental result. In contrast, the Zimm-Bragg theory predicts the DSC and CD unfolding transitions with the same set of parameters.

Analysis of self-assembly of S-layer protein slp-B53 from Lysinibacillus sphaericus

The formation of stable and functional surface layers (S-layers) via self-assembly of surface-layer proteins on the cell surface is a dynamic and complex process. S-layers facilitate a number of important biological functions, e.g., providing protection and mediating selective exchange of molecules and thereby functioning as molecular sieves. Furthermore, S-layers selectively bind several metal ions including uranium, palladium, gold, and europium, some of them with high affinity. Most current research on surface layers focuses on investigating crystalline arrays of protein subunits in Archaea and bacteria. In this work, several complementary analytical techniques and methods have been applied to examine structure-function relationships and dynamics for assembly of S-layer protein slp-B53 from Lysinibacillus sphaericus: (1) The secondary structure of the S-layer protein was analyzed by circular dichroism spectroscopy; (2) Small-angle X-ray scattering was applied to gain insights into the three-dimensional structure in solution; (3) The interaction with bivalent cations was followed by differential scanning calorimetry; (4) The dynamics and time-dependent assembly of S-layers were followed by applying dynamic light scattering; (5) The two-dimensional structure of the paracrystalline S-layer lattice was examined by atomic force microscopy. The data obtained provide essential structural insights into the mechanism of S-layer self-assembly, particularly with respect to binding of bivalent cations, i.e., Mg²⁺ and Ca²⁺. Furthermore, the results obtained highlight potential applications of S-layers in the fields of micromaterials and nanobiotechnology by providing engineered or individual symmetric thin protein layers, e.g., for protective, antimicrobial, or otherwise functionalized surfaces.

Role of denatured-state properties in chaperonin action probed by single-molecule spectroscopy

The bacterial chaperonin GroEL/GroES assists folding of a broad spectrum of denatured and misfolded proteins. Here, we explore the limits of this remarkable promiscuity by mapping two denatured proteins with very different conformational properties, rhodanese and cyclophilin A, during binding and encapsulation by GroEL/GroES with single-molecule spectroscopy, microfluidic mixing, and ensemble kinetics. We find that both proteins bind to GroEL with high affinity in a reaction involving substantial conformational adaptation. However, whereas the compact denatured state of rhodanese is encapsulated efficiently upon addition of GroES and ATP, the more expanded and unstructured denatured cyclophilin A is not encapsulated but is expelled into solution. The origin of this surprising disparity is the weaker interactions of cyclophilin A with a transiently formed GroEL-GroES complex, which may serve as a crucial checkpoint for substrate discrimination.

The kinetic stability of cytochrome C oxidase: effect of bound phospholipid and dimerization

Thermally induced transitions of the 13-subunit integral membrane protein bovine cytochrome c oxidase (CcO) have been studied by differential scanning calorimetry (DSC) and circular dichroism (CD). Thermal denaturation of dodecyl maltoside solubilized CcO proceeds in two consecutive, irreversible, kinetically driven steps with the apparent transition temperatures at ∼ 51°C and ∼ 61°C (5μM CcO at scan rate of 1.5 K/min). The thermal denaturation data were analyzed according to the Lyubarev and Kurganov model of two consecutive irreversible steps. However, because of the limitation of the model to describe the complex mechanism of the thermal denaturation of CcO, the obtained results were utilized only for comparison purposes of kinetic stabilities of CcO under specific protein concentration (5μM) and scan rate (1.5 K/min). This enabled us to show that both the amphiphilic environment and the self-association state of CcO affect its kinetic stability. Kinetic stabilities of both steps are significantly decreased when all of the phospholipids are removed from CcO by phospholipase A2 (the half-life decreases at 37°C). Conversely, dimerization of CcO induced by sodium cholate significantly increases its kinetic stability of only the first step (the half-life increases at 37°C). Protein concentration-dependent nonspecific oligomerization also indicate mild stabilization of CcO. Both, reversed-phase high-performance liquid chromatography (HPLC) and SDS-PAGE subunit analysis reveal that the first step of thermal denaturation involves dissociation of subunits III, VIa, VIb, and VIIa, whereas the second step is less well defined and most likely involves global unfold and aggregation of the remaining subunits. Electron transport activity of CcO decreases in a sigmoidal manner during the first transition and this dependence is very well described by kinetic parameters for the first step of the thermal transition. Therefore, dissociation of subunit III and/or VIIa is responsible for temperature-induced inactivation of CcO because VIa and VIb can be removed from CcO without affecting the enzyme activity. These results demonstrate an important role of tightly bound phospholipids and oligomeric state (particularly the dimeric form) of CcO for kinetic stability of the protein.

The Biokinetic Spectrum for Temperature

We identify and describe the distribution of temperature-dependent specific growth rates for life on Earth, which we term the biokinetic spectrum for temperature. The spectrum has the potential to provide for more robust modeling in thermal ecology since any conclusions derived from it will be based on observed data rather than using theoretical assumptions. It may also provide constraints for systems biology model predictions and provide insights in physiology. The spectrum has a Δ-shape with a sharp peak at around 42°C. At higher temperatures up to 60°C there was a gap of attenuated growth rates. We found another peak at 67°C and a steady decline in maximum rates thereafter. By using Bayesian quantile regression to summarise and explore the data we were able to conclude that the gap represented an actual biological transition between mesophiles and thermophiles that we term the Mesophile-Thermophile Gap (MTG). We have not identified any organism that grows above the maximum rate of the spectrum. We used a thermodynamic model to recover the Δ-shape, suggesting that the growth rate limits arise from a trade-off between activity and stability of proteins. The spectrum provides underpinning principles that will find utility in models concerned with the thermal responses of biological processes.

Protein stability, conformational change and binding mechanism of human serum albumin upon binding of embelin and its role in disease control

Here, we present the inclusive binding mode of phytochemical embelin, an anticancer drug with human serum albumin (HSA) established under physiological condition. Also, to understand the pharmacological role of embelin molecule, here, we have studied the anti-cancer activity of embelin on human cervical cancer cell line (HeLa cell line), which revealed that embelin showed dose dependent inhibition in the growth of cancer cells and also induces 26.3% of apoptosis at an IC50 value of 29μM. Further, embelin was titrated with HSA and the fluorescence emission quenching of HSA due to the formation of the HSA-embelin complex was observed. The binding constant of this complex is 5.9±.01×10(4)M(-1) and the number of bound embelin molecules is approximately 1.0. Consequently, molecular displacement and computational docking experiments show that the embelin is binding to subdomain IB to HSA. Further evidence from microTOF-Q mass spectrometry showed an increase in mass from 66,563Da to 66,857Da observed for free HSA and HSA+embelin complex, signifying that there is robust binding of embelin with HSA. In addition, the variations of HSA secondary structural elements in presence of embelin were confirmed by circular dichroism which indicates partial unfolding of protein. Furthermore, the transmission electron micrographs established that complex formation leads to aggregation of HSA plus embelin. Molecular dynamics simulations revealed that the stability of the HSA-embelin complexes and results suggests that at around 3500ps the complex reaches equilibration state which clearly contributes to the understanding of the stability of the HSA-embelin complexes.

Hydration Properties and Solvent Effects for All-Atom Solutes in Polarizable Coarse-Grained Water

Due to the importance of water in chemical and biological systems, a coarse-grained representation of the solvent can greatly simplify the description of the system while retaining key thermodynamic properties of the medium. A multiscale solvation model that couples all-atom solutes and polarizable Martini coarse-grained water (AAX/CGS) is developed to reproduce free energies of hydration of organic solutes. Using Monte Carlo/free energy perturbation (MC/FEP) calculations, results from multiscale and all-atom simulations are compared. Improved accuracy is obtained with the AAX/CGS approach for hydrophobic and sulfur- or halogen-containing solutes, but larger deviations are found for polar solute molecules where hydrogen bonding is featured. Furthermore, solvent effects on conformational and tautomeric equilibria of AA solutes were investigated using AA, CG, and GB/SA solvent models. It is found that the CG solvent model can reproduce well the medium effects from experiment and AA simulations; however, the GB/SA solvent model fails in some cases. A 7-30-fold reduction in computational cost is found for the present AAX/CGS multiscale simulations compared to the AA alternative.

Molecular mechanism of thermosensory function of human heat shock transcription factor Hsf1

The heat shock response is a universal homeostatic cell autonomous reaction of organisms to cope with adverse environmental conditions. In mammalian cells, this response is mediated by the heat shock transcription factor Hsf1, which is monomeric in unstressed cells and upon activation trimerizes, and binds to promoters of heat shock genes. To understand the basic principle of Hsf1 activation we analyzed temperature-induced alterations in the conformational dynamics of Hsf1 by hydrogen exchange mass spectrometry. We found a temperature-dependent unfolding of Hsf1 in the regulatory region happening concomitant to tighter packing in the trimerization region. The transition to the active DNA binding-competent state occurred highly cooperative and was concentration dependent. Surprisingly, Hsp90, known to inhibit Hsf1 activation, lowered the midpoint temperature of trimerization and reduced cooperativity of the process thus widening the response window. Based on our data we propose a kinetic model of Hsf1 trimerization.

Cold Denaturation Unveiled: Molecular Mechanism of the Asymmetric Unfolding of Yeast Frataxin

What is the mechanism that determines the denaturation of proteins at low temperatures, which is, by now, recognized as a fundamental property of all proteins? We present experimental evidence that clarifies the role of specific interactions that favor the entrance of water into the hydrophobic core, a mechanism originally proposed by Privalov but never proved experimentally. By using a combination of molecular dynamics simulation, molecular biology, and biophysics, we identified a cluster of negatively charged residues that represents a preferential gate for the entrance of water molecules into the core. Even single-residue mutations in this cluster, from acidic to neutral residues, affect cold denaturation much more than heat denaturation, suppressing cold denaturation at temperatures above zero degrees. The molecular mechanism of the cold denaturation of yeast frataxin is intrinsically different from that of heat denaturation.

Thermodynamics of protein denaturation at temperatures over 100 °C: CutA1 mutant proteins substituted with hydrophobic and charged residues

Although the thermodynamics of protein denaturation at temperatures over 100 °C is essential for the rational design of highly stable proteins, it is not understood well because of the associated technical difficulties. We designed certain hydrophobic mutant proteins of CutA1 from Escherichia coli, which have denaturation temperatures (Td) ranging from 101 to 113 °C and show a reversible heat denaturation. Using a hydrophobic mutant as a template, we successfully designed a hyperthermostable mutant protein (Td = 137 °C) by substituting six residues with charged ones. Thermodynamic analyses of these mutant proteins indicated that the hydrophobic mutants were stabilized by the accumulation of denaturation enthalpy (ΔH) with no entropic gain from hydrophobic solvation around 100 °C, and that the stabilization due to salt bridges resulted from both the increase in ΔH from ion-ion interactions and the entropic effect of the electrostatic solvation over 113 °C. This is the first experimental evidence that has successfully overcome the typical technical difficulties.

"Hydration Shells" of CH2 Groups of ω-Amino Acids as Studied by Deuteron NMR Relaxation

Hydration phenomena play a very important role in various processes, in particular in biological systems. Water molecules in aqueous solutions of organic compounds can be distributed among the following substructures: (i) hydration shells of hydrophilic functional groups of molecules, (ii) water in the environment of nonpolar moieties, and (iii) bulk water. Up to now, the values of hydration parameters suggested for the description of various solutions of organic compounds were not thoroughly analyzed in the aspect of the consideration of the total molecular composition. The temperature and concentration dependences of relaxation rates of water deuterons were studied in a wide range of concentration and temperature in aqueous (D2O) solutions of a set of ω-amino acids. Assuming the coordination number of the CH2 group equal to 7, which was determined from quantum-chemical calculations, it was found that the rotational correlation times of water molecules near the methylene group is 1.5-2 times greater than one for pure water. The average rotational mobility of water molecules in the hydration shells of hydrophilic groups of ω-amino acids is a bit slower than that in pure solvent at temperatures higher that 60 °C, but at lower temperatures, it is 0.8-1.0 of values of correlation times for bulk water. The technique suggested provides the basis for the characterization of different hydrophobic and hydrophilic species in the convenient terms of the rotational correlation times for the nearest water molecules.

The role of intra-domain disulfide bonds in heat-induced irreversible denaturation of camelid single domain VHH antibodies

Camelid-derived single domain VHH antibodies are highly heat resistant, and the mechanism of heat-induced VHH denaturation predominantly relies on the chemical modification of amino acids. Although chemical modification of disulfide bonds has been recognized as a cause for heat-induced denaturation of many proteins, there have been no mutagenesis studies, in which the number of disulfide bonds was controlled. In this article, we examined a series of mutants of two different VHHs with single, double or no disulfide bonds, and scrutinized the effects of these disulfide bond modifications on VHH denaturation. With the exception of one mutant, the heat resistance of VHHs decreased when the number of disulfide bonds increased. The effect of disulfide bonds on heat denaturation was more striking if the VHH had a second disulfide bond, suggesting that the contribution of disulfide shuffling is significant in proteins with multiple disulfide bonds. Furthermore, our results directly indicate that removal of a disulfide bond can indeed increase the heat resistance of a protein, irrespective of the negative impact on equilibrium thermodynamic stability.

Proline substitutions in a Mip-like peptidyl-prolyl cis-trans isomerase severely affect its structure, stability, shape and activity

FKBP22, an Escherichia coli-specific peptidyl-prolyl cis-trans isomerase, shows substantial homology with the Mip-like virulence factors. Mip-like proteins are homodimeric and possess a V-shaped conformation. Their N-terminal domains form dimers, whereas their C-terminal domains bind protein/peptide substrates and distinct inhibitors such as rapamycin and FK506. Interestingly, the two domains of the Mip-like proteins are separated by a lengthy, protease-susceptible α-helix. To delineate the structural requirement of this domain-connecting region in Mip-like proteins, we have investigated a recombinant FKBP22 (rFKBP22) and its three point mutants I65P, V72P and A82P using different probes. Each mutant harbors a Pro substitution mutation at a distinct location in the hinge region. We report that the three mutants are not only different from each other but also different from rFKBP22 in structure and activity. Unlike rFKBP22, the three mutants were unfolded by a non-two state mechanism in the presence of urea. In addition, the stabilities of the mutants, particularly I65P and V72P, differed considerably from that of rFKBP22. Conversely, the rapamycin binding affinity of no mutant was different from that of rFKBP22. Of the mutants, I65P showed the highest levels of structural/functional loss and dissociated partly in solution. Our computational study indicated a severe collapse of the V-shape in I65P due to the anomalous movement of its C-terminal domains. The α-helical nature of the domain-connecting region is, therefore, critical for the Mip-like proteins.

Structural Rigidity and Protein Thermostability in Variants of Lipase A from Bacillus subtilis

Understanding the origin of thermostability is of fundamental importance in protein biochemistry. Opposing views on increased or decreased structural rigidity of the folded state have been put forward in this context. They have been related to differences in the temporal resolution of experiments and computations that probe atomic mobility. Here, we find a significant (p = 0.004) and fair (R² = 0.46) correlation between the structural rigidity of a well-characterized set of 16 mutants of lipase A from Bacillus subtilis (BsLipA) and their thermodynamic thermostability. We apply the rigidity theory-based Constraint Network Analysis (CNA) approach, analyzing directly and in a time-independent manner the statics of the BsLipA mutants. We carefully validate the CNA results on macroscopic and microscopic experimental observables and probe for their sensitivity with respect to input structures. Furthermore, we introduce a robust, local stability measure for predicting thermodynamic thermostability. Our results complement work that showed for pairs of homologous proteins that raising the structural stability is the most common way to obtain a higher thermostability. Furthermore, they demonstrate that related series of mutants with only a small number of mutations can be successfully analyzed by CNA, which suggests that CNA can be applied prospectively in rational protein design aimed at higher thermodynamic thermostability.

Regulation of human 3-beta-hydroxysteroid dehydrogenase type-2 (3βHSD2) by molecular chaperones and the mitochondrial environment affects steroidogenesis

Human 3-β-hydroxysteroid dehydrogenase/isomerase types 1 and 2 (3βHSD1 and 3βHSD2, respectively) are expressed in a tissue-specific pattern by different genes. Site-directed mutagenesis studies have confirmed the function of the catalytic amino acids (Tyr154, Lys 158, Ser124 in both isoenzymes), substrate/inhibitor isoform-specific residues (His156 and Arg195 in 3βHSD1) and cofactor binding residues (Asp36 provides NAD(+) specificity in both isoenzymes). However, detailed analysis of isoform-specific organelle localization and characterization is difficult due to the 93% amino acid identity between the two isoforms. With recent advances in the knowledge of mitochondrial architecture and localization of the various translocases, our laboratory has studied the mechanisms regulating mitochondrial 3βHSD2 localization. The mitochondrial N-terminal leader sequence of 3βHSD2 directs its entry into the mitochondria where it is localized to the intermembrane space (IMS). Unlike other mitochondrial proteins, the N-terminal signal sequence of 3βHSD2 is not cleaved upon mitochondrial import. 3βHSD2 interacts with the mitochondrial translocase, Tim50, to regulate progesterone and androstenedione formation. Our studies suggest that its activity at the IMS is facilitated in a partially unfolded "molten globule" conformation by the proton pump between the matrix and IMS. The unfolded protein is refolded by the mitochondrial chaperones. The protons at the IMS are absorbed by the lipid vesicles, to maintain the proton pump and recycle 3βHSD2. As a result, one molecule of 3βHSD2 may participate in multiple catalytic reactions. In summary, the steroidogenic cell recycles 3βHSD2 to catalyze the reactions needed to produce androstenedione, progesterone and 17α-hydroxyprogesterone on demand in coordination with the mitochondrial translocase, Tim50. This article is part of a Special Issue entitled 'Steroid/Sterol signaling'.

The Thermodynamic Mechanism of Peptide-MHC Class II Complex Formation Is a Determinant of Susceptibility to HLA-DM

Peptides bind MHC class II molecules through a thermodynamically nonadditive process consequent to the flexibility of the reactants. Currently, how the specific outcome of this binding process affects the ensuing epitope selection needs resolution. Calorimetric assessment of binding thermodynamics for hemagglutinin 306-319 peptide variants to the human MHC class II HLA-DR1 (DR1) and a mutant DR1 reveals that peptide/DR1 complexes can be formed with different enthalpic and entropic contributions. Complexes formed with a smaller entropic penalty feature circular dichroism spectra consistent with a non-compact form, and molecular dynamics simulation shows a more flexible structure. The opposite binding mode, compact and less flexible, is associated with greater entropic penalty. These structural variations are associated with rearrangements of residues known to be involved in HLA-DR (DM) binding, affinity of DM for the complex, and complex susceptibility to DM-mediated peptide exchange. Thus, the thermodynamic mechanism of peptide binding to DR1 correlates with the structural rigidity of the complex, and DM mediates peptide exchange by "sensing" flexible complexes in which the aforementioned residues are rearranged at a higher frequency than in more rigid ones.

Carnitine in bacterial physiology and metabolism

Carnitine is a quaternary amine compound found at high concentration in animal tissues, particularly muscle, and is most well studied for its contribution to fatty acid transport into mitochondria. In bacteria, carnitine is an important osmoprotectant, and can also enhance thermotolerance, cryotolerance and barotolerance. Carnitine can be transported into the cell or acquired from metabolic precursors, where it can serve directly as a compatible solute for stress protection or be metabolized through one of a few distinct pathways as a nutrient source. In this review, we summarize what is known about carnitine physiology and metabolism in bacteria. In particular, recent advances in the aerobic and anaerobic metabolic pathways as well as the use of carnitine as an electron acceptor have addressed some long-standing questions in the field.

The Effect of Surface Charge Saturation on Heat-induced Aggregation of Firefly Luciferase

We present here the effect of firefly luciferase surface charge saturation and the presence of some additives on its thermal-induced aggregation. Three mutants of firefly luciferase prepared by introduction of surface Arg residues named as 2R, 3R and 5R have two, three and five additional arginine residues substituted at their surface compared to native luciferase; respectively. Turbidimetric study of heat-induced aggregation indicates that all three mutants were reproducibly aggregated at higher rates relative to wild type in spite of their higher thermostability. Among them, 2R had most evaluated propensity to heat-induced aggregation. Therefore, the hydrophilization followed by appearing of more substituted arginine residues with positive charge on the firefly luciferase surface was not reduced its thermal aggregation. Nevertheless, at the same condition in the presence of charged amino acids, e.g. Arg, Lys and Glu, as well as a hydrophobic amino acid, e.g. Val, the heat-induced aggregation of wild type and mutants of firefly luciferases was markedly decelerated than those in the absence of additives. On the basis of obtained results it seems, relinquishment of variety in charge of amino acid side chains, they via local interactions with proteins cause to decrease rate and extent of their thermal aggregation.

Biophysics of protein evolution and evolutionary protein biophysics

The study of molecular evolution at the level of protein-coding genes often entails comparing large datasets of sequences to infer their evolutionary relationships. Despite the importance of a protein's structure and conformational dynamics to its function and thus its fitness, common phylogenetic methods embody minimal biophysical knowledge of proteins. To underscore the biophysical constraints on natural selection, we survey effects of protein mutations, highlighting the physical basis for marginal stability of natural globular proteins and how requirement for kinetic stability and avoidance of misfolding and misinteractions might have affected protein evolution. The biophysical underpinnings of these effects have been addressed by models with an explicit coarse-grained spatial representation of the polypeptide chain. Sequence-structure mappings based on such models are powerful conceptual tools that rationalize mutational robustness, evolvability, epistasis, promiscuous function performed by 'hidden' conformational states, resolution of adaptive conflicts and conformational switches in the evolution from one protein fold to another. Recently, protein biophysics has been applied to derive more accurate evolutionary accounts of sequence data. Methods have also been developed to exploit sequence-based evolutionary information to predict biophysical behaviours of proteins. The success of these approaches demonstrates a deep synergy between the fields of protein biophysics and protein evolution.

Universal distribution of mutational effects on protein stability, uncoupling of protein robustness from sequence evolution and distinct evolutionary modes of prokaryotic and eukaryotic proteins

Robustness to destabilizing effects of mutations is thought of as a key factor of protein evolution. The connections between two measures of robustness, the relative core size and the computationally estimated effect of mutations on protein stability (ΔΔG), protein abundance and the selection pressure on protein-coding genes (dN/dS) were analyzed for the organisms with a large number of available protein structures including four eukaryotes, two bacteria and one archaeon. The distribution of the effects of mutations in the core on protein stability is universal and indistinguishable in eukaryotes and bacteria, centered at slightly destabilizing amino acid replacements, and with a heavy tail of more strongly destabilizing replacements. The distribution of mutational effects in the hyperthermophilic archaeon Thermococcus gammatolerans is significantly shifted toward strongly destabilizing replacements which is indicative of stronger constraints that are imposed on proteins in hyperthermophiles. The median effect of mutations is strongly, positively correlated with the relative core size, in evidence of the congruence between the two measures of protein robustness. However, both measures show only limited correlations to the expression level and selection pressure on protein-coding genes. Thus, the degree of robustness reflected in the universal distribution of mutational effects appears to be a fundamental, ancient feature of globular protein folds whereas the observed variations are largely neutral and uncoupled from short term protein evolution. A weak anticorrelation between protein core size and selection pressure is observed only for surface residues in prokaryotes but a stronger anticorrelation is observed for all residues in eukaryotic proteins. This substantial difference between proteins of prokaryotes and eukaryotes is likely to stem from the demonstrable higher compactness of prokaryotic proteins.

Influence of Hydrophobic Face Amino Acids on the Hydrogelation of β-Hairpin Peptide Amphiphiles

Hydrophobic residues provide much of the thermodynamic driving force for the folding, self-assembly, and consequent hydrogelation of amphiphilic β-hairpin peptides. We investigate how the identity of hydrophobic side chains displayed from the hydrophobic face of these amphiphilic peptides influences their behavior to expound on the design criteria important to gel formation. Six peptides were designed that globally incorporate valine, aminobutyric acid, norvaline, norleucine, phenylalanine, or isoleucine on the hydrophobic face of the hairpin to study how systematic changes in hydrophobic content, β-sheet propensity, and aromaticity affect gelation. Circular dichroism (CD) spectroscopy indicates that hydrophobic content, rather than β-sheet propensity, dictates the temperature- and pH-dependent folding and assembly behavior of these peptides. Transmission electron microscopy (TEM) and small-angle neutron scattering (SANS) show that the local morphology of the fibrils formed via self-assembly is little affected by amino acid type. However, residue type does influence the propensity of peptide fibrils to undergo higher order assembly events. Oscillatory rheology shows that the mechanical rigidity of the peptide gels is highly influenced by residue type, but there is no apparent correlation between rigidity and residue hydrophobicity nor β-sheet propensity. Lastly, the large planar aromatic side chain of phenylalanine supports hairpin folding and assembly, affording a gel characterized by a rate of formation and storage modulus similar to the parent valine-containing peptide.

Protein denaturation at a single-molecule level: the effect of nonpolar environments and its implications on the unfolding mechanism by proteases

Most proteins are typically folded into predetermined three-dimensional structures in the aqueous cellular environment. However, proteins can be exposed to a nonpolar environment under certain conditions, such as inside the central cavity of chaperones and unfoldases during protein degradation. It remains unclear how folded proteins behave when moved from an aqueous solvent to a nonpolar one. Here, we employed single-molecule atomic force microscopy and molecular dynamics (MD) simulations to investigate the structural and mechanical variations of a polyprotein, I278, during the change from a polar to a nonpolar environment. We found that the polyprotein was unfolded into an unstructured polypeptide spontaneously when pulled into nonpolar solvents. This finding was corroborated by MD simulations where I27 was dragged from water into a nonpolar solvent, revealing details of the unfolding process at the water/nonpolar solvent interface. These results highlight the importance of water in maintaining folding stability, and provide insights into the response of folded proteins to local hydrophobic environments.