A Euclidean Perspective on the Unfolding of Azurin: Spatial Correlations

Structural stabilities of calcium proteins: Human intelectin-1 and frog lectin XEEL

We extend our study of the structural stability of helical and nonhelical regions in chain A of human intelectin-1 to include a second human intelectin (4WMY) and the frog protein "Xenopus embryonic epidermal lectin" (XEEL). These unique lectins have been shown to recognize carbohydrate residues found exclusively in microbes, thus they could potentially be developed into novel microbe detection and sequestration tools. We believe that by studying the structural stability of these proteins we can provide insights on their biological role and activities. Using a geometrical model introduced previously, we perform computational analyses of protein crystal structures that quantify the resiliency of the native state to steric perturbations. Based on these analyses, we conclude that differences in the resiliency of the human and frog proteins can be attributed primarily to differences in non-helical regions and to residues near Ca ions. Since these differences are particularly pronounced in the vicinity of the ligand binding site, they provide an explanation for the finding that human intelectin-1 has a higher affinity for a ligand than XEEL. We also present data on conserved and position-equivalent pairs of residues in 4WMY and XEEL. We identify residue pairs as well as regions in which the influence of neighboring residues is nearly uniform as the parent protein denatures. Since the structural signatures are conserved, this identification provides a basis for understanding why both proteins exhibit trimeric structures despite poor sequence conservation at the interface.

Relaxation of structural constraints during Amicyanin unfolding

We study the thermal unfolding of amicyanin by quantifying the resiliency of the native state to structural perturbations. Three signatures characterizing stages of unfolding are identified. The first signature, lateral extension of the polypeptide chain, is calculated directly from the reported crystallographic data. Two other signatures, the radial displacement of each residue from Cu(II) and the angular spread in the chain as the protein unfolds, are calculated using crystallographic data in concert with a geometrical model we introduced previously (J.J. Kozak, H. B. Gray, R. A. Garza-López, J. Inorg. Biochem. 155(2016) 44-55). Particular attention is paid to the resiliency of the two beta sheets in amicyanin. The resiliency of residues in the near neighborhood of the Cu center to destabilization provides information on the persistence of the entatic state. Similarly, examining the resiliency of residues intercalated between structured regions (beta sheets, the alpha helix) provides a basis for identifying a "hydrophobic core." A principal focus of our study is to compare results obtained using our geometrical model with the experimental results (C. La Rosa, D. Milardi, D. M. Grasso, M. P. Verbeet, G. W. Canters, L. Sportelli, R. Guzzi, Eur. Biophy. J.30(8),(2002) 559-570) on the denaturation of amicyanin, and we show that our results support a classical model proposed by these authors.

Structural Stability of Intelectin-1

We study the structural stability of helical and nonhelical regions in chain A of human intelectin-1. Using a geometrical model introduced previously, we carried out a computational analysis based on the recently reported crystal structure of this protein by Kiessling et al. to quantify the resiliency of the native state to steric perturbations. Response to these perturbations is characterized by calculating, relative to the native state, the lateral, radial, and angular displacements of n-residue segments of the polypeptide chain centered on each residue. By quantifying the stability of the protein through six stages of unfolding, we are able to identify regions in chain A of intelectin-1 that are markedly affected by structural perturbations versus those that are relatively unaffected, the latter suggesting that the native-state geometry of these regions is essentially conserved. Importantly, residues in the vicinity of calcium ions comprise a conserved region, suggesting that Ca ions play a role not only in the coordination of carbohydrate hydroxyl groups but also in preserving the integrity of the structure.

Cytochrome unfolding pathways from computational analysis of crystal structures

We have developed a model to study the role of geometrical factors in influencing the early stages of unfolding in three cytochromes: cyt c', cyt c-b562 and cyt c. Each stage in unfolding is quantified by the spatial extension λ̂i of n-residue segments, and by their angular extension 〈βn〉. Similarities and differences between and among the three cytochromes in the unfolding of helical and non-helical regions can be determined by analyzing the data for each signature separately. Definite conclusions can be drawn when spatial and angular changes are considered in tandem. To facilitate comparisons, we present graphical portraits of the three cytochromes at the same stage of unfolding, and in relation to their native state structures. We also display specific segments at different stages of unfolding to illustrate differences in stability of defined domains thereby allowing us to make specific predictions on the unfolding of corresponding internal and terminal helices in cyt c' and cyt c-b562. Our work accords with an earlier experimental report on the presence and persistence of a hydrophobic core in cyt c.

A Euclidean perspective on the unfolding of azurin: chain motion

We present a new approach to visualizing and quantifying the displacement of segments of Pseudomonas aeruginosa azurin in the early stages of denaturation. Our method is based on a geometrical method developed previously by the authors, and elaborated extensively for azurin. In this study, we quantify directional changes in three α-helical regions, two regions having β-strand residues, and three unstructured regions of azurin. Snapshots of these changes as the protein unfolds are displayed and described quantitatively by introducing a scaling diagnostic. In accord with molecular dynamics simulations, we show that the long α-helix in azurin (residues 54-67) is displaced from the polypeptide scaffolding and then pivots first in one direction, and then in the opposite direction as the protein continues to unfold. The two β-strand chains remain essentially intact and, except in the earliest stages, move in tandem. We show that unstructured regions 72-81 and 84-91, hinged by β-strand residues 82-83, pivot oppositely. The region comprising residues 72-91 (40 % hydrophobic and 16 % of the 128 total residues) forms an effectively stationary region that persists as the protein unfolds. This static behavior is a consequence of a dynamic balance between the competing motion of two segments, residues 72-81 and 84-91.