Dissecting the roles of local packing density and longer‐range effects in protein sequence evolution

What are the structural determinants of protein sequence evolution? A number of site‐specific structural characteristics have been proposed, most of which are broadly related to either the density of contacts or the solvent accessibility of individual residues. Most importantly, there has been disagreement in the literature over the relative importance of solvent accessibility and local packing density for explaining site‐specific sequence variability in proteins. We show that this discussion has been confounded by the definition of local packing density. The most commonly used measures of local packing, such as contact number and the weighted contact number, represent the combined effects of local packing density and longer‐range effects. As an alternative, we propose a truly local measure of packing density around a single residue, based on the Voronoi cell volume. We show that the Voronoi cell volume, when calculated relative to the geometric center of amino‐acid side chains, behaves nearly identically to the relative solvent accessibility, and each individually can explain, on average, approximately 34% of the site‐specific variation in evolutionary rate in a data set of 209 enzymes. An additional 10% of variation can be explained by nonlocal effects that are captured in the weighted contact number. Consequently, evolutionary variation at a site is determined by the combined effects of the immediate amino‐acid neighbors of that site and effects mediated by more distant amino acids. We conclude that instead of contrasting solvent accessibility and local packing density, future research should emphasize on the relative importance of immediate contacts and longer‐range effects on evolutionary variation. Proteins 2016; 84:841–854. © 2016 Wiley Periodicals, Inc.     

Temperature-sensitive mutations of bacteriophage T4 lysozyme occur at sites with low mobility and low solvent accessibility in the folded protein

Twenty-five different temperature-sensitive point mutations at 20 sites in the lysozyme gene of bacteriophage T4 have been identified. All of the mutations alter amino acid side chains that have lower than average crystallographic thermal factors and reduced solvent accessibility in the folded protein. This suggests that the amino acids with well-defined conformations can form specific intramolecular interactions that make relatively large contributions to the thermal stability of the protein. Residues with high mobility or high solvent accessibility are much less susceptible to destabilizing substitutions, suggesting that, in general, such amino acids contribute less to protein stability. The pattern of the sites of ts substitutions observed in the folded conformation of T4 lysozyme suggests that severe destabilizing mutations that primarily affect the free energy of the unfolded state are rare. These results indicate that proteins can be stabilized by adding new interactions to regions that are rigid or buried in the folded conformation.     

Prediction of local structural stabilities of proteins from their amino acid sequences

Hydrogen exchange experiments provide detailed information about the local stability and the solvent accessibility of different regions of the structures of folded proteins, protein complexes, and amyloid fibrils. We introduce an approach to predict protection factors from hydrogen exchange in proteins based on the knowledge of their amino acid sequences without the inclusion of any additional structural information. These results suggest that the propensity of different regions of the structures of globular proteins to undergo local unfolding events can be predicted from their amino acid sequences with an accuracy of 80% or better.     

Structural adaptation of enzymes to low temperatures

A systematic comparative analysis of 21 psychrophilic enzymes belonging to different structural families from prokaryotic and eukaryotic organisms is reported. The sequences of these enzymes were multiply aligned to 427 homologous proteins from mesophiles and thermophiles. The net flux of amino acid exchanges from meso/thermophilic to psychrophilic enzymes was measured. To assign the observed preferred exchanges to different structural environments, such as secondary structure, solvent accessibility and subunit interfaces, homology modeling was utilized to predict the secondary structure and accessibility of amino acid residues for the psychrophilic enzymes for which no experimental three-dimensional structure is available. Our results show a clear tendency for the charged residues Arg and Glu to be replaced at exposed sites on alpha-helices by Lys and Ala, respectively, in the direction from 'hot' to 'cold' enzymes. Val is replaced by Ala at buried regions in alpha-helices. Compositional analysis of psychrophilic enzymes shows a significant increase in Ala and Asn and a decrease in Arg at exposed sites. Buried sites in beta-strands tend to be depleted of VAL: Possible implications of the observed structural variations for protein stability and engineering are discussed.     

A comparative analysis of the spatial structure of nonspecific porins from <Emphasis Type="Italic">Yersinia ruckeri</Emphasis> using optical spectroscopy and molecular modeling

The spatial organization of outer-membrane porins is studied by optical spectroscopy and molecular modeling. It was found that the OmpF and OmpC porins from Yеrsiniа ruckeri are β-structured membrane proteins typical of the pore-forming proteins of other Gram-negative bacteria. The spatial structures of monomers and trimers of the OmpC and OmpF porins from Y. ruckeri are simulated using methods of structural bioinformatics. It was found that the structural stability of the more thermostable OmpF trimer is sustained by a greater number of hydrogen bonds and hydrophobic interactions. The main differences of the spatial structures of the test porins are observed in the structure of their outer loops. There are three tryptophan residues in the molecules of the OmpC and OmpF porins of Y. ruckeri. It is demonstrated by moleculardynamics methods that after thermal denaturation the solvent accessibility of the Trp212 residue in OmpF porin increased by two times, while the solvent accessibility of a Trp184 residue in OmpC porin was not increased. It is hypothesized that the red-shifted tryptophan fluorescence spectrum of OmpF porin during thermal denaturation is due to the behavior of the Trp212 residue.     

Increasing the thermostability of staphylococcal nuclease: implications for the origin of protein thermostability

Seven hyper-stable multiple mutants have been constructed in staphylococcal nuclease by various combinations of eight different stabilizing single mutants. The stabilities of these multiple mutants determined by guanidine hydrochloride denaturation were 3.4 to 5.6 kcal/mol higher than that of the wild-type. Their thermal denaturation midpoint temperatures were 12.6 to 22.9 deg. C higher than that of the wild-type. These are among the greatest increases in protein stability and thermal denaturation midpoint temperature relative to the wild-type yet attained. There has been great interest in understanding how proteins found in thermophilic organisms are stabilized. One frequently cited theory is that the packing of hydrophobic side-chains is improved in the cores of proteins isolated from thermophiles when compared to proteins from mesophiles. The crystal structures of four single and five multiple stabilizing mutants of staphylococcal nuclease were solved to high resolution. No large overall structural change was found, with most changes localized around the sites of mutation. Rearrangements were observed in the packing of side-chains in the major hydrophobic core, although none of the mutations was in the core. It is surprising that detailed structural analysis showed that packing had improved, with the volume of the mutant protein's hydrophobic cores decreasing as protein stability increased. Further, the number of van der Waals interactions in the entire protein showed an experimentally significant increase correlated with increasing stability. These results indicate that optimization of packing follows as a natural consequence of increased protein thermostability and that good packing is not necessarily the proximate cause of high stability. Another popular theory is that thermostable proteins have more electrostatic and hydrogen bonding interactions and these are responsible for the high stabilities. The mutants here show that increased numbers of electrostatic and hydrogen bonding interactions are not obligatory for large increases in protein stability.     

Important inter-residue contacts for enhancing the thermal stability of thermophilic proteins

Proteins from thermophilic organisms exhibit high thermal stability, but have structures that are very similar to their mesophilic homologues. In order to gain insight into the basis of thermostability, we have analyzed the medium- and long-range contacts in mesophilic and thermophilic proteins of 16 different families. We found that the thermophiles prefer to have contacts between residues with hydrogen-bond-forming capability. Apart from hydrophobic contacts, more contacts are observed between polar and non-polar residues in thermophiles than mesophiles. Residue-wise analysis showed that Tyr has good contacts with several other residues, and Cys has considerably higher long-range contacts in thermophiles compared with mesophiles. Furthermore, the residues occurring in the range of 31-34 residues apart in the sequence contribute significant long-range contacts to the stability of thermophilic proteins.     

Fragmentation of monoclonal antibodies

Fragmentation is a degradation pathway ubiquitously observed in proteins despite the remarkable stability of peptide bond; proteins differ only by how much and where cleavage occurs. The goal of this review is to summarize reports regarding the non-enzymatic fragmentation of the peptide backbone of monoclonal antibodies (mAbs). The sites in the polypeptide chain susceptible to fragmentation are determined by a multitude of factors. Insights are provided on the intimate chemical mechanisms that can make some bonds prone to cleavage due to the presence of specific side-chains. In addition to primary structure, the secondary, tertiary and quaternary structures have a significant impact in modulating the distribution of cleavage sites by altering local flexibility, accessibility to solvent or bringing in close proximity side chains that are remote in sequence. This review focuses on cleavage sites observed in the constant regions of mAbs, with special emphasis on hinge fragmentation. The mechanisms responsible for backbone cleavage are strongly dependent on pH and can be catalyzed by metals or radicals. The distribution of cleavage sites are different under acidic compared to basic conditions, with fragmentation rates exhibiting a minimum in the pH range 5 to 6; therefore, the overall fragmentation pattern observed for a mAb is a complex result of structural and solvent conditions. A critical review of the techniques used to monitor fragmentation is also presented; usually a compromise has to be made between a highly sensitive method with good fragment separation and the capability to identify the cleavage site. The effect of fragmentation on the function of a mAb must be evaluated on a case-by-case basis depending on whether cleavage sites are observed in the variable or constant regions, and on the mechanism of action of the molecule.     

Protein rigidity and thermophilic adaptation

We probe the hypothesis of corresponding states, according to which homologues from mesophilic and thermophilic organisms are in corresponding states of similar rigidity and flexibility at their respective optimal temperatures. For this, the local distribution of flexible and rigid regions in 19 pairs of homologous proteins from meso- and thermophilic organisms is analyzed and related to activity characteristics of the enzymes by constraint network analysis (CNA). Two pairs of enzymes are considered in more detail: 3-isopropylmalate dehydrogenase and thermolysin-like protease. By comparing microscopic stability features of homologues with the help of stability maps, introduced for the first time, we show that adaptive mutations in enzymes from thermophilic organisms maintain the balance between overall rigidity, important for thermostability, and local flexibility, important for activity, at the appropriate working temperature. Thermophilic adaptation in general leads to an increase of structural rigidity but conserves the distribution of functionally important flexible regions between homologues. This finding provides direct evidence for the hypothesis of corresponding states. CNA thereby implicitly captures and unifies many different mechanisms that contribute to increased thermostability and to activity at high temperatures. This allows to qualitatively relate changes in the flexibility of active site regions, induced either by a temperature change or by the introduction of mutations, to experimentally observed losses of the enzyme function. As for applications, the results demonstrate that exploiting the principle of corresponding states not only allows for successful thermostability optimization but also for guiding experiments in order to improve enzyme activity in protein engineering.     

Quantitative Mapping of Protein Structure by Hydroxyl Radical Footprinting-Mediated Structural Mass Spectrometry: A Protection Factor?Analysis

Measurements from hydroxyl radical footprinting (HRF) provide rich information about the solvent accessibility of amino acid side chains of a protein. Traditional HRF data analyses focus on comparing the difference in the modification/footprinting rate of a specific site to infer structural changes across two protein states, e.g., between a free and ligand-bound state. However, the rate information itself is not fully used for the purpose of comparing different protein sites within a protein on an absolute scale. To provide such a cross-site comparison, we present a new, to our knowledge, data analysis algorithm to convert the measured footprinting rate constant to a protection factor (PF) by taking into account the known intrinsic reactivity of amino acid side chain. To examine the extent to which PFs can be used for structural interpretation, this PF analysis is applied to three model systems where radiolytic footprinting data are reported in the literature. By visualizing structures colored with the PF values for individual peptides, a rational view of the structural features of various protein sites regarding their solvent accessibility is revealed, where high-PF regions are buried and low-PF regions are more exposed to the solvent. Furthermore, a detailed analysis correlating solvent accessibility and local structural contacts for gelsolin shows a statistically significant agreement between PF values and various structure measures, demonstrating that the PFs derived from this PF analysis readily explain fundamental HRF rate measurements. We also tested this PF analysis on alternative, chemical-based HRF data, showing improved correlations of structural properties of a model protein barstar compared to examining HRF rate data alone. Together, this PF analysis not only permits a novel, to our knowledge, approach of mapping protein structures by using footprinting data, but also elevates the use of HRF measurements from a qualitative, cross-state comparison to a quantitative, cross-site assessment of protein structures in the context of individual conformational states of interest.     

Dynamic alpha-helices: conformations that do not conform

Structural transitions are important for the stability and function of proteins, but these phenomena are poorly understood. An extensive analysis of Protein Data Bank entries reveals 103 regions in proteins with a tendency to transform from helical to nonhelical conformation and vice versa. We find that these dynamic helices, unlike other helices, are depleted in hydrophobic residues. Furthermore, the dynamic helices have higher surface accessibility and conformational mobility (P-value = 3.35e-07) than the rigid helices. Contact analyses show that these transitions result from protein-ligand, protein-nucleic acid, and crystal-contacts. The immediate structural environment differs quantitatively (P-value = 0.003) as well as qualitatively in the two alternate conformations. Often, dynamic helix experiences more contacts in its helical conformation than in the nonhelical counterpart (P-value = 0.001). There is differential preference for the type of short contacts observed in two conformational states. We also demonstrate that the regions in protein that can undergo such large conformational transitions can be predicted with a reasonable accuracy using logistic regression model of supervised learning. Our findings have implications in understanding the molecular basis of structural transitions that are coupled with binding and are important for the function and stability of the protein. Based on our observations, we propose that several functionally relevant regions on the protein surface can switch over their conformation from coil to helix and vice-versa, to regulate the recognition and binding of their partner and hence these may work as "molecular switches" in the proteins to regulate certain biological process. Our results supports the idea that protein structure-function paradigm should transform from static to a highly dynamic one.     

Structural differences between mesophilic, moderately thermophilic and extremely thermophilic protein subunits: results of a comprehensive survey

BACKGROUND: Proteins from thermophilic organisms usually show high intrinsic thermal stability but have structures that are very similar to their mesophilic homologues. From prevous studies it is difficult to draw general conclusions about the structural features underlying the increased thermal stability of thermophilic proteins. RESULTS: In order to reveal the general evolutionary strategy for changing the heat stability of proteins, a non-redundant data set was compiled comprising all high-quality structures of thermophilic proteins and their mesophilic homologues from the Protein Data Bank. The selection (quality) criteria were met by 64 mesophilic and 29 thermophilic protein subunits, representing 25 protein families. From the atomic coordinates, 13 structural parameters were calculated, compared and evaluated using statistical methods. This study is distinguished from earlier ones by the strict quality control of the structures used and the size of the data set. CONCLUSIONS: Different protein families adapt to higher temperatures by different sets of structural devices. Regarding the structural parameters, the only generally observed rule is an increase in the number of ion pairs with increasing growth temperature. Other parameters show just a trend, whereas the number of hydrogen bonds and the polarity of buried surfaces exhibit no clear-cut tendency to change with growth temperature. Proteins from extreme thermophiles are stabilized in different ways to moderately thermophilic ones. The preferences of these two groups are different with regards to the number of ion pairs, the number of cavities, the polarity of exposed surface and the secondary structural composition.     

Role of solvation barriers in protein kinetic stability

The stability of several protein systems of interest has been shown to have a kinetic basis. Besides the obvious biotechnological implications, the general interest of understanding protein kinetic stability is emphasized by the fact that some emerging molecular approaches to the inhibition of amyloidogenesis focus on the increase of the kinetic stability of protein native states. Lipases are among the most important industrial enzymes. Here, we have studied the thermal denaturation of the wild-type form, four single-mutant variants and two highly stable, multiple-mutant variants of lipase from Thermomyces lanuginosa. In all cases, thermal denaturation was irreversible, kinetically controlled and conformed to the two-state irreversible model. This result supports that the novel molecular-dynamics-focused, directed-evolution approach involved in the preparation of the highly stable variants is successful likely because it addresses kinetic stability and, in particular, because heated molecular dynamics simulations possibly identify regions of disrupted native interactions in the transition state for irreversible denaturation. Furthermore, we find very large mutation effects on activation enthalpy and entropy, which were not accompanied by similarly large changes in kinetic urea m-value. From this we are led to conclude that these mutation effects are associated to some structural feature of the transition state for the irreversible denaturation process that is not linked to large changes in solvent accessibility. Recent computational studies have suggested the existence of solvation/desolvation barriers in at least some protein folding/unfolding processes. We thus propose that a solvation barrier (arising from the asynchrony between breaking of internal contacts and water penetration) may contribute to the kinetic stability of lipase from T. lanuginosa (and, possibly, to the kinetic stability of other proteins as well).     

Characterization of Protein–Protein Interfaces

We analyze the characteristics of protein–protein interfaces using the largest datasets available from the Protein Data Bank (PDB). We start with a comparison of interfaces with protein cores and non-interface surfaces. The results show that interfaces differ from protein cores and non-interface surfaces in residue composition, sequence entropy, and secondary structure. Since interfaces, protein cores, and non-interface surfaces have different solvent accessibilities, it is important to investigate whether the observed differences are due to the differences in solvent accessibility or differences in functionality. We separate out the effect of solvent accessibility by comparing interfaces with a set of residues having the same solvent accessibility as the interfaces. This strategy reveals residue distribution propensities that are not observable by comparing interfaces with protein cores and non-interface surfaces. Our conclusions are that there are larger numbers of hydrophobic residues, particularly aromatic residues, in interfaces, and the interactions apparently favored in interfaces include the opposite charge pairs and hydrophobic pairs. Surprisingly, Pro-Trp pairs are over represented in interfaces, presumably because of favorable geometries. The analysis is repeated using three datasets having different constraints on sequence similarity and structure quality. Consistent results are obtained across these datasets. We have also investigated separately the characteristics of heteromeric interfaces and homomeric interfaces.     

Inter-residue interactions in protein folding and stability

During the process of protein folding, the amino acid residues along the polypeptide chain interact with each other in a cooperative manner to form the stable native structure. The knowledge about inter-residue interactions in protein structures is very helpful to understand the mechanism of protein folding and stability. In this review, we introduce the classification of inter-residue interactions into short, medium and long range based on a simple geometric approach. The features of these interactions in different structural classes of globular and membrane proteins, and in various folds have been delineated. The development of contact potentials and the application of inter-residue contacts for predicting the structural class and secondary structures of globular proteins, solvent accessibility, fold recognition and ab initio tertiary structure prediction have been evaluated. Further, the relationship between inter-residue contacts and protein-folding rates has been highlighted. Moreover, the importance of inter-residue interactions in protein-folding kinetics and for understanding the stability of proteins has been discussed. In essence, the information gained from the studies on inter-residue interactions provides valuable insights for understanding protein folding and de novo protein design.     

Replacement of staphylococcal nuclease hydrophobic core residues with those from thermophilic homologues indicates packing is improved in some thermostable proteins

The importance of tight hydrophobic core packing in stabilizing proteins found in thermophilic organisms has been vigorously disputed. Here, portions of the cores found in three thermophilic homologues were transplanted into the core of staphylococcal nuclease, a protein of modest stability. Packing of the core was evaluated by comparing interaction energy of the three mutants to the comprehensive mutant library built up previously at these same sites in staphylococcal nuclease. It was found that the interaction energy of one thermophilic sequence is extraordinarily favorable and the interaction energies of other two transplanted thermophilic sequences are good, comparable to the interaction energies of mutant cores based on cores found in mesophilic homologues. As expected when transferring just a portion of the core sequence, the mutant proteins were destabilized overall relative to wild-type staphylococcal nuclease. The overall conclusion is that improvement of packing interactions is a mechanism to confer stability employed in some proteins from thermophiles, but not all.     

Fragmentation of monoclonal antibodies

Fragmentation is a degradation pathway ubiquitously observed in proteins despite the remarkable stability of peptide bond; proteins differ only by how much and where cleavage occurs. The goal of this review is to summarize reports regarding the non-enzymatic fragmentation of the peptide backbone of monoclonal antibodies (mAbs). The sites in the polypeptide chain susceptible to fragmentation are determined by a multitude of factors. Insights are provided on the intimate chemical mechanisms that can make some bonds prone to cleavage due to the presence of specific side-chains. In addition to primary structure, the secondary, tertiary and quaternary structures have a significant impact in modulating the distribution of cleavage sites by altering local flexibility, accessibility to solvent or bringing in close proximity side chains that are remote in sequence. This review focuses on cleavage sites observed in the constant regions of mAbs, with special emphasis on hinge fragmentation. The mechanisms responsible for backbone cleavage are strongly dependent on pH and can be catalyzed by metals or radicals. The distribution of cleavage sites are different under acidic compared to basic conditions, with fragmentation rates exhibiting a minimum in the pH range 5–6; therefore, the overall fragmentation pattern observed for a mAb is a complex result of structural and solvent conditions. A critical review of the techniques used to monitor fragmentation is also presented; usually a compromise has to be made between a highly sensitive method with good fragment separation and the capability to identify the cleavage site. The effect of fragmentation on the function of a mAb must be evaluated on a case-by-case basis depending on whether cleavage sites are observed in the variable or constant regions, and on the mechanism of action of the molecule.Key words: fragmentation, cleavage, clipping, hinge region, peptide bond hydrolysis, IgG1, IgG2     

Latency and substrate binding globally reduce solvent accessibility of plasminogen activator inhibitor type 1 (PAI-1). An adaptation of PAI-1 conformer crystal structures by hydrogen-deuterium exchange

Plasminogen activator inhibitor type 1 (PAI-1) plays key regulatory roles in fibrinolysis, cell migration, and tissue remodeling. A regulatory protein without known catalytic activity, PAI-1 modulates plasminogen activators through protein-protein interactions. Although global conformational alterations that occur in PAI-1 determine its regulatory activity, comprehensive assessments of concurrent dynamic, structural, and functional alterations of this critical regulatory protein have not yet been clearly defined. X-ray crystallographic studies have described four distinct PAI-1 conformational states: active, latent, reactive center loop peptide-annealed (RCL-PA), and cleaved mutant. In this study, backbone amide hydrogen-deuterium exchange detected by mass spectrometry was used to characterize dynamic and structural alterations of human PAI-1 (hPAI-1) in relation to its function. hPAI-1 conformers were defined by surface mapping the solvent-accessible sites for strategic secondary structural components of the protein. We observed a global protection from solvent for a majority of peptides in the latent conformer relative to the active conformer. Significant differences were observed in the RCL, helix A, helix D, and sheet 1C, and these regions were markedly more dynamic or solvent-exposed in the active conformation. The RCL-PA form adopts an intermediate conformational state between the active and the latent conformers. Our results demonstrate that the most dynamic regions of PAI-1 (the RCL, helices D and A, and sheet 5A) are flexible in the transition toward latency. They also show that the dynamic surface structures of the active, latent, and peptide-annealed conformers of PAI-1 are underestimated by theoretical solvent accessibility calculations derived from crystallographic data.     

Allosteric changes in solvent accessibility observed in thrombin upon active site occupation

The solvent accessibility of thrombin in its substrate-free and substrate-bound forms has been compared by amide hydrogen/deuterium (H/(2)H) exchange. The optimized inhibitor peptide dPhe-Pro-Arg chloromethyl ketone (PPACK) was used to simulate the substrate-bound form of thrombin. These studies were motivated by the lack of observed changes in the active site of thrombin in the crystal structure of the thrombin-thrombomodulin complex. This result appeared to contradict amide exchange studies on the thrombin-thrombomodulin complex that suggested subtle changes occur in the active site loops upon thrombomodulin binding. Our results show that two active site loops, residues 214-222 and residues 126-132, undergo decreases in solvent accessibility due to steric contacts with PPACK substrate. However, we also observe two regions outside the active site undergoing solvent protection upon substrate binding. The first region corresponds to anion binding exosite 1, and the second is a beta-strand-containing loop which runs through the core of the molecule and contains Trp141 which makes critical contacts with anion binding exosite 1. These results indicate two pathways of allosteric change that connect the active site to the distal anion binding exosite 1.     

Investigations on the role of π–π interactions and π–π networks in eNOS and nNOS proteins

π–π Interactions play an important role in the stability of protein structures. In the present study, we have analyzed the influence of π–π interactions in eNOS and nNOS proteins. The contribution of these π–π interacting residues in sequential separation, secondary structure involvement, solvent accessibility and stabilization centers has been evaluated. π–π interactions stabilize the core regions within eNOS and nNOS proteins. π–π interacting residues are evolutionary conserved. There is a significant number of π–π interactions in spite of the lesser natural occurrences of π-residues in eNOS and nNOS proteins. In addition to π–π interactions, π residues also form π–π networks in both eNOS and nNOS proteins which might play an important role in the structural stability of these protein structures.