A network representation of protein structures: Implications for protein stability

This study views each protein structure as a network of noncovalent connections between amino acid side chains. Each amino acid in a protein structure is a node, and the strength of the noncovalent interactions between two amino acids is evaluated for edge determination. The protein structure graphs (PSGs) for 232 proteins have been constructed as a function of the cutoff of the amino acid interaction strength at a few carefully chosen values. Analysis of such PSGs constructed on the basis of edge weights has shown the following: 1), The PSGs exhibit a complex topological network behavior, which is dependent on the interaction cutoff chosen for PSG construction. 2), A transition is observed at a critical interaction cutoff, in all the proteins, as monitored by the size of the largest cluster (giant component) in the graph. Amazingly, this transition occurs within a narrow range of interaction cutoff for all the proteins, irrespective of the size or the fold topology. And 3), the amino acid preferences to be highly connected (hub frequency) have been evaluated as a function of the interaction cutoff. We observe that the aromatic residues along with arginine, histidine, and methionine act as strong hubs at high interaction cutoffs, whereas the hydrophobic leucine and isoleucine residues get added to these hubs at low interaction cutoffs, forming weak hubs. The hubs identified are found to play a role in bringing together different secondary structural elements in the tertiary structure of the proteins. They are also found to contribute to the additional stability of the thermophilic proteins when compared to their mesophilic counterparts and hence could be crucial for the folding and stability of the unique three-dimensional structure of proteins. Based on these results, we also predict a few residues in the thermophilic and mesophilic proteins that can be mutated to alter their thermal stability.

The Critical Role of N- and C-Terminal Contact in Protein Stability and Folding of a Family 10 Xylanase under Extreme Conditions

Background: Stabilization strategies adopted by proteins under extreme conditions are very complex and involve various kinds of interactions. Recent studies have shown that a large proportion of proteins have their N- and C-terminal elements in close contact and suggested they play a role in protein folding and stability. However, the biological significance of this contact remains elusive. Methodology: In the present study, we investigate the role of N- and C-terminal residue interaction using a family 10 xylanase (BSX) with a TIM-barrel structure that shows stability under high temperature,alkali pH, and protease and SDS treatment. Based on crystal structure,an aromatic cluster was identified that involves Phe4, Trp6 and Tyr343 holding the Nand C-terminus together; this is a unique and important feature of this protein that might be crucial for folding and stabilityunder poly-extreme conditions. Conclusion: A series of mutants was created to disrupt this aromatic cluster formation and study the loss of stability and function under given conditions. While the deletions of Phe4 resulted in loss of stability, removal of Trp6 and Tyr343 affected in vivo folding and activity. Alanine substitution with Phe4, Trp6 and Tyr343 drastically decreased stability under all parameters studied. Importantly,substitution of Phe4 with Trp increased stability in SDS treatment.Mass spectrometry results of limited proteolysis further demonstrated that the Arg344 residue is highly susceptible to trypsin digestion in sensitive mutants such as DF4, W6A and Y343A, suggesting again that disruption of the Phe4-Trp6-Tyr343 (F-W-Y) cluster destabilizes the N-and C-terminal interaction. Our results underscore the importance of N- and C-terminal contact through aromatic interactions in protein folding and stability under extreme conditions, and these results may be useful to improve the stability of other proteins under suboptimal conditions.

Enhanced J-protein interaction and compromised protein stability of mtHsp70 variants lead to mitochondrial dysfunction in Parkinsons disease

Parkinsons disease (PD) is the second most prevalent progressive neurological disorder commonly associated with impaired mitochondrial function in dopaminergic neurons. Although familial PD is multifactorial in nature, a recent genetic screen involving PD patients identified two mitochondrial Hsp70 variants (P509S and R126W) that are suggested in PD pathogenesis. However, molecular mechanisms underlying how mtHsp70 PD variants are centrally involved in PD progression is totally elusive. In this article, we provide mechanistic insights into the mitochondrial dysfunction associated with human mtHsp70 PD variants. Biochemically, the R126W variant showed severely compromised protein stability and was found highly susceptible to aggregation at physiological conditions. Strikingly, on the other hand, the P509S variant exhibits significantly enhanced interaction with J-protein cochaperones involved in folding and import machinery, thus altering the overall regulation of chaperone-mediated folding cycle and protein homeostasis. To assess the impact of mtHsp70 PD mutations at the cellular level, we developed yeast as a model system by making analogous mutations in Ssc1 ortholog. Interestingly, PD mutations in yeast (R103W and P486S) exhibit multiple in vivo phenotypes, which are associated with omitochondrial dysfunction', including compromised growth, impairment in protein translocation, reduced functional mitochondrial mass, mitochondrial DNA loss, respiratory incompetency and increased susceptibility to oxidative stress. In addition to that, R103W protein is prone to aggregate in vivo due to reduced stability, whereas P486S showed enhanced interaction with J-proteins, thus remarkably recapitulating the cellular defects that are observed in human PD variants. Taken together, our findings provide evidence in favor of direct involvement of mtHsp70 as a susceptibility factor in PD.

Effect of Crowding Agents, Signal Peptide, and Chaperone SecB on the Folding and Aggregation of E. coli Maltose Binding Protein

SecB, a soluble cytosolic chaperone component of the Secexport pathway, binds to newly synthesized precursor proteins and prevents their premature aggregation and folding and subsequently targets them to the translocation machinery on the membrane. PreMBP, the precursor form of maltose binding protein, has a 26-residue signal sequence attached to the N-terminus of MBP and is a physiological substrate of SecB. We examine the effect of macromolecular crowding and SecB on the stability and refolding of denatured preMBP and MBP. PreMBP was less stable than MBP (ΔTm =7( 0.5 K) in both crowded and uncrowded solutions. Crowding did not cause any substantial changes in the thermal stability ofMBP(ΔTm=1(0.4 K) or preMBP (ΔTm=0(0.6 K), as observed in spectroscopically monitored thermal unfolding experiments. However, both MBP and preMBP were prone to aggregation while refolding under crowded conditions. In contrast to MBP aggregates, which were amorphous, preMBP aggregates form amyloid fibrils.Under uncrowded conditions, a molar excess of SecB was able to completely prevent aggregation and promote disaggregation of preformed aggregates of MBP. When a complex of the denatured protein and SecB was preformed, SecB could completely prevent aggregation and promote folding of MBP and preMBP even in crowded solution. Thus, in addition to maintaining substrates in an unfolded, export-competent conformation, SecB also suppresses the aggregation of its substrates in the crowded intracellular environment. SecB is also able to promote passive disaggregation of macroscopic aggregates of MBP in the absence of an energy source such as ATP or additional cofactors. These experiments also demonstrate that signal peptide can reatly influence protein stability and aggregation propensity.

X-ray and molecular-dynamics studies on Mycobacterium leprae single-stranded DNA-binding protein and comparison with other eubacterial SSB structures

The crystal structures of two forms of Mycobacterium leprae single-stranded DNA-binding protein (SSB) have been determined at 2.05 and 2.8 A resolution. Comparison of these structures with the structures of other eubacterial SSBs indicates considerable variation in their quaternary association, although the DNA-binding domains in all of them exhibit the same OB-fold. This variation has no linear correlation with sequence variation, but could be related to variation in protein stability. Molecular-dynamics simulations have been carried out on tetrameric molecules derived from the two forms and the prototype Escherichia coli SSB and the individual subunits of both proteins. Together, the X-ray studies and molecular-dynamics simulations yield information on the relatively rigid and flexible regions of the molecule and on the effect of oligomerization on flexibility. The simulations provide insight into the changes in subunit structure on oligomerization. They also provide insight into the stability and time evolution of the hydrogen bonds/water bridges that connect the two pairs of monomers in the tetramer.

Effect of signal peptide on stability and folding of escherichia coli thioredoxin

The signal peptide plays a key role in targeting and membrane insertion of secretory and membrane proteins in both prokaryotes and eukaryotes. In E. coli, recombinant proteins can be targeted to the periplasmic space by fusing naturally occurring signal sequences to their N-terminus. The model protein thioredoxin was fused at its N-terminus with malE and pelB signal sequences. While WT and the pelB fusion are soluble when expressed, the malE fusion was targeted to inclusion bodies and was refolded in vitro to yield a monomeric product with identical secondary structure to WT thioredoxin. The purified recombinant proteins were studied with respect to their thermodynamic stability, aggregation propensity and activity, and compared with wild type thioredoxin, without a signal sequence. The presence of signal sequences leads to thermodynamic destabilization, reduces the activity and increases the aggregation propensity, with malE having much larger effects than pelB. These studies show that besides acting as address labels, signal sequences can modulate protein stability and aggregation in a sequence dependent manner.

PROTEIN DISULFIDE ANALYSIS AND DESIGN

The significant contribution of naturally occurring disulfide bonds to protein stability has encouraged development of methods to engineer non-native disulfides in proteins. These have yielded mixed results. We summarize applications of the program MODIP for disulfide engineering. The program predicts sites in proteins where disulfides can be stably introduced. The program has also been used as an aid in conformational analysis of naturally occurring disulfides in a-helices, antiparallel and parallel beta-strands. Disulfides in a-helices occur only at N-termini, where the first cysteine residue is the N-cap residue of the helix. The disulfide occurs as a CXXC motif and can possess redox activity. In antiparallel beta-strands, disulfides occur exclusively at non-hydrogen bonded (NHB) registered pairs of antiparallel beta-sheets with only 1 known natural example occurring at a hydrogen bonded (HB) registered pair. Conformational analysis suggests that disulfides between HB residue pairs are under torsional strain. A similar analysis to characterize disulfides in parallel beta-strands was carried out. We observed that only 9 instances of cross-strand disulfides exist in a non-redundant dataset. Stereochemical analysis shows that while tbe chi(square) angles are similar to those of other disulfides, the chi(1) and chi(2) angles show more variation and that one of tbe strands is generally an edge strand.

Impact of glycosylation on stability, structure and unfolding of soybean agglutinin (SBA): an insight from thermal perturbation molecular dynamics simulations

Glycosylation has been recognized as one of the most prevalent and complex post-translational modifications of proteins involving numerous enzymes and substrates. Its effect on the protein conformational transitions is not clearly understood yet. In this study, we have examined the effect of glycosylation on protein stability using molecular dynamics simulation of legume lectin soybean agglutinin (SBA). Its glycosylated moiety consists of high mannose type N-linked glycan (Man(9)GlcNAc(2)). To unveil the structural perturbations during thermal unfolding of these two forms, we have studied and compared them to the experimental results. From the perspective of dynamics, our simulations revealed that the nonglycosylated monomeric form is less stable than corresponding glycosylated form at normal and elevated temperatures. Moreover, at elevated temperature thermal destabilization is more prominent in solvent exposed loops, turns and ends of distinct beta sheets. SBA maintains it folded structure due to some important saltbridges, hydrogen bonds and hydrophobic interactions within the protein. The reducing terminal GlcNAc residues interact with the protein residues VAL161, PRO182 and SER225 via hydrophobic and via hydrogen bonding with ASN 9 and ASN 75. Our simulations also revealed that single glycosylation (ASN75) has no significant effect on corresponding cis peptide angle orientation. This atomistic description might have important implications for understanding the functionality and stability of Soybean agglutinin.

Unfolding of Plasmodium falciparum Triosephosphate Isomerase in Urea and Guanidinium Chloride: Evidence for a Novel Disulfide Exchange Reaction in a Covalently Cross-Linked Mutant†

The conformational stability of Plasmodium falciparum triosephosphate isomerase (TIMWT) enzyme has been investigated in urea and guanidinium chloride (GdmCl) solutions using circular dichroism, fluorescence, and size-exclusion chromatography. The dimeric enzyme is remarkably stable in urea solutions. It retains considerable secondary, tertiary, and quaternary structure even in 8 M urea. In contrast, the unfolding transition is complete by 2.4 M GdmCl. Although the secondary as well as the tertiary interactions melt before the perturbation of the quaternary structure, these studies imply that the dissociation of the dimer into monomers ultimately leads to the collapse of the structure, suggesting that the interfacial interactions play a major role in determining multimeric protein stability. The C-m(urea)/C-m(GdmCl) ratio (where C-m is the concentration of the denaturant required at the transition midpoint) is unusually high for triosephosphate isomerase as compared to other monomeric and dimeric proteins. A disulfide crosslinked mutant protein (Y74C) engineered to form two disulfide cross-links across the interface (13-74') and (13'-74) is dramatically destablized in urea. The unfolding transition is complete by 6 M urea and involves a novel mechanism of dimer dissociation through intramolecular thiol-disulfide exchange.

Thermodynamic and structural studies of cavity formation in proteins suggest that loss of packing interactions rather than the hydrophobic effect dominates the observed energetics

The hydrophobic effect is widely believed to be an important determinant of protein stability. However, it is difficult to obtain unambiguous experimental estimates of the contribution of the hydrophobic driving force to the overall free energy of folding. Thermodynamic and structural studies of large to small substitutions in proteins are the most direct method of measuring this contribution. We have substituted the buried residue Phe8 in RNase S with alanine, methionine, and norleucine, Binding thermodynamics and structures were characterized by titration calorimetry and crystallography, respectively. The crystal structures of the RNase S F8A, F8M, and F8Nle mutants indicate that the protein tolerates the changes without any main chain adjustments, The correlation of structural and thermodynamic parameters associated with large to small substitutions was analyzed for nine mutants of RNase S as well as 32 additional cavity-containing mutants of T4 lysozyme, human lysozyme, and barnase. Such substitutions were typically found to result in negligible changes in Delta C-p and positive values of both Delta Delta H degrees and aas of folding. Enthalpic effects were dominant, and the sign of Delta Delta S is the opposite of that expected from the hydrophobic effect. Values of Delta Delta G degrees and Delta Delta H degrees correlated better with changes in packing parameters such as residue depth or occluded surface than with the change in accessible surface area upon folding. These results suggest that the loss of packing interactions rather than the hydrophobic effect is a dominant contributor to the observed energetics for large to small substitutions. Hence, estimates of the magnitude of the hydrophobic driving force derived from earlier mutational studies are likely to be significantly in excess of the actual value.

Unfolding of Plasmodium falciparum triosephosphate isomerase in urea and guanidinium chloride solutions. Evidence for a novel disulfide exchange reaction in a covalently crosslinked mutant

The conformational stability of Plasmodium falciparum triosephosphate isomerase (TIMWT) enzyme has been investigated in urea and guanidinium chloride (GdmCl) solutions using circular dichroism, fluorescence, and size-exclusion chromatography. The dimeric enzyme is remarkably stable in urea solutions. It retains considerable secondary, tertiary, and quaternary structure even in 8 M urea. In contrast, the unfolding transition is complete by 2.4 M GdmCl. Although the secondary as well as the tertiary interactions melt before the perturbation of the quaternary structure, these studies imply that the dissociation of the dimer into monomers ultimately leads to the collapse of the structure, suggesting that the interfacial interactions play a major role in determining multimeric protein stability. The Cm(urea)/Cm(GdmCl) ratio (where Cm is the concentration of the denaturant required at the transition midpoint) is unusually high for triosephosphate isomerase as compared to other monomeric and dimeric proteins. A disulfide cross-linked mutant protein (Y74C) engineered to form two disulfide cross-links across the interface (13-74‘) and (13‘-74) is dramatically destablized in urea. The unfolding transition is complete by 6 M urea and involves a novel mechanism of dimer dissociation through intramolecular thiol−disulfide exchange.

Conformational analysis and design of cross-strand disulfides in antiparallel beta-sheets

Cross-strand disulfides bridge two cysteines in a registered pair of antiparallel beta-strands. A nonredundant data set comprising 5025 polypeptides containing 2311 disulfides was used to study cross-strand disulfides. Seventy-six cross-strand disulfides were found of which 75 and 1 occurred at non-hydrogen-bonded (NHB) and hydrogen-bonded (HB) registered pairs, respectively. Conformational analysis and modeling studies demonstrated that disulfide formation at HB pairs necessarily requires an extremely rare and positive chi(1) value for at least one of the cysteine residues. Disulfides at HB positions also have more unfavorable steric repulsion with the main chain. Thirteen pairs of disulfides were introduced in NHB and HB pairs in four model proteins: leucine binding protein (LBP), leucine, isoleucine, valine binding protein (LIVBP), maltose binding protein (MBP), and Top7. All mutants LIVBP T247C V331C showed disulfide formation either on purification, or on treatment with oxidants. Protein stability in both oxidized and reduced states of all mutants was measured. Relative to wild type, LBP and MBP mutants were destabilized with respect to chemical denaturation, although the sole exposed NHB LBP mutant showed an increase of 3.1 degrees C in T-m. All Top7 mutants were characterized for stability through guanidinium thiocyanate chemical denaturation. Both exposed and two of the three buried NHB mutants were appreciably stabilized. All four HB Top7 mutants were destabilized (Delta Delta G(0) = -3.3 to -6.7 kcal/mol). The data demonstrate that introduction of cross-strand disulfides at exposed NHB pairs is a robust method of improving protein stability. All four exposed Top7 disulfide mutants showed mild redox activity. Proteins 2011; 79: 244-260. (C) 2010 Wiley-Liss, Inc.

Thermodynamic and structural consequences of changing a sulfur atom to a methylene group in the M13Nle mutation in ribonuclease-S

Two fragments of pancreatic ribonuclease A, a truncated version of S-peptide (residues 1-15) and S-protein (residues 21-124), combine to give a catalytically active complex. We have substituted the wild-type residue at position 13, methionine (Met), with norleucine (Nle), where the only covalent change is the replacement of the sulfur atom with a methylene group. The thermodynamic parameters associated with the binding of this variant to S-protein, determined by titration calorimetry in the temperature range 10-40 degrees C, are reported and compared to values previously reported [Varadarajan, R., Connelly, P. R., Sturtevant, J. M., & Richards, F. M. (1992) Biochemistry 31, 1421-1426] for other position 13 analogs. The differences in the free energy and enthalpy of binding between the Met and Nle peptides are 0.6 and 7.9 kcal/mol at 25 degrees C, respectively. These differences are slightly larger than, but comparable to, the differences in the values for the Met/Ile and Met/Leu pairs. The structure of the mutant complex was determined to 1.85 Angstrom resolution and refined to an R-factor of 17.4% The structures of mutant and wild-type complexes are practically identical although the Nle side chain has a significantly higher average B-factor than the corresponding Met side chain. In contrast, the B-factors of the atoms of the cage of residues surrounding position 13 are all somewhat lower in the Nle variant than in the Met wild-type. Thus, the large differences in the binding enthalpy appear to reside entirely in the difference in chemical properties or dynamic behavior of the -S- and -CH2- groups and not in differences in the geometry of the side chains or the internal cavity surface. In addition, a novel method of obtaining protein stability data by means of isothermal titration calorimetry is introduced.

DEPTH: a web server to compute depth and predict small-molecule binding cavities in proteins

Depth measures the extent of atom/residue burial within a protein. It correlates with properties such as protein stability, hydrogen exchange rate, protein-protein interaction hot spots, post-translational modification sites and sequence variability. Our server, DEPTH, accurately computes depth and solvent-accessible surface area (SASA) values. We show that depth can be used to predict small molecule ligand binding cavities in proteins. Often, some of the residues lining a ligand binding cavity are both deep and solvent exposed. Using the depth-SASA pair values for a residue, its likelihood to form part of a small molecule binding cavity is estimated. The parameters of the method were calibrated over a training set of 900 high-resolution X-ray crystal structures of single-domain proteins bound to small molecules (molecular weight < 1.5 KDa). The prediction accuracy of DEPTH is comparable to that of other geometry-based prediction methods including LIGSITE, SURFNET and Pocket-Finder (all with Matthew's correlation coefficient of similar to 0.4) over a testing set of 225 single and multi-chain protein structures. Users have the option of tuning several parameters to detect cavities of different sizes, for example, geometrically flat binding sites. The input to the server is a protein 3D structure in PDB format. The users have the option of tuning the values of four parameters associated with the computation of residue depth and the prediction of binding cavities. The computed depths, SASA and binding cavity predictions are displayed in 2D plots and mapped onto 3D representations of the protein structure using Jmol. Links are provided to download the outputs. Our server is useful for all structural analysis based on residue depth and SASA, such as guiding site-directed mutagenesis experiments and small molecule docking exercises, in the context of protein functional annotation and drug discovery.

Packing in molten globules and native states

Close packing of hydrophobic residues in the protein interior is an important determinant of protein stability. Cavities introduced by large to small substitutions are known to destabilize proteins. Conversely, native states of proteins and protein fragments can be stabilized by filling in existing cavities. Molten globules (MGs) were initially used to describe a state of protein which has well-defined secondary structure but little or no tertiary packing. Subsequent studies have shown that MGs do have some degree of native-like topology and specific packing. Wet molten globules (WMGs) with hydrated cores and considerably decreased packing relative to the native state have been studied extensively. Recently there has been renewed interest in identification and characterization of dry molten globules (DMGs). These are slightly expanded forms of the native state which show increased conformational flexibility, native-like main-chain hydrogen bonding and dry interiors. The generality of occurrence of DMGs during protein unfolding and the extent and nature of packing in DMGs remain to be elucidated. Packing interactions in native proteins and MGs can be probed through mutations. Next generation sequencing technologies make it possible to determine relative populations of mutants in a large pool. When this is coupled to phenotypic screens or cell-surface display, it becomes possible to rapidly examine large panels of single-site or multi-site mutants. From such studies, residue specific contributions to protein stability and function can be estimated in a highly parallelized fashion. This complements conventional biophysical methods for characterization of packing in native states and molten globules.