Cytochrome unfolding pathways from computational analysis of crystal structures

Unfolding cytochromes c-b562 and Rd apo b562

We have analyzed the early stages of unfolding of cytochromes c-b₅₆₂ (PDB ID: 2BC5) and Rd apo b₅₆₂ (PDB ID: 1YYJ). Our geometrical approach proceeds from an analysis of the crystal structure reported for each protein. We quantify, residue-by-residue and region-by-region, the spatial and angular changes in the structure as the protein denatures, and quantify differences that result from the seven residues that differ in the two proteins. Using two independent analyses, one based on spatial metrics and the second on angular metrics, we establish the order of unfolding of the five helices in cyt c-b₅₆₂ and the four helices in the apo protein. For the two helices nearest the N-terminal end of both proteins, the ones in the apo protein unfold first. For the two helices nearest the C-terminal end, the interior helix of the apo protein unfolds first, whereas the terminal helix of the holo protein unfolds first. Excluded-volume effects (repulsive interactions) are minimized in turning regions; the overall range in Δ values is Δ = 36.3 Å³ for cyt c-b₅₆₂ and Δ = 36.6 Å³ for the apo protein, whereas the span for all 20 amino acids is Δ = 167.7 Å³. As our work indicates that the interior helix of cytochrome c-b₅₆₂ is the first to fold, we suggest that this helix protects the heme from misligation, consistent with ultrafast folding over a minimally frustrated funneled landscape.

Funneled angle landscapes for helical proteins

We use crystallographic data for four helical iron proteins (cytochrome c-b₅₆₂, cytochrome c', sperm whale myoglobin, human cytoglobin) to calculate radial and angular signatures as each unfolds from the native state stepwise though four unfolded states. From these data we construct an angle phase diagram to display the evolution of each protein from its native state; and, in turn, the phase diagram is used to construct a funneled angle landscape for comparison with the topography of its folding energy landscape. We quantify the departure of individual helical and turning regions from the areal, angular profile of corresponding regions of the native state. This procedure allows us to identify the similarities and differences among individual helical and turning regions in the early stages of unfolding of the four helical heme proteins.

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.