The response of Ω-loop D dynamics to truncation of trimethyllysine 72 of yeast iso-1-cytochrome c depends on the nature of loop deformation

Applications of the Newly Developed Force-Field Parameters Uncover a Dynamic Nature of Ω-Loop C in the Lys-Ligated Alkaline Form of Cytochrome c

Lys-ligated cytochromes make up an emerging family of heme proteins. Density functional theory calculations on the amine/imidazole-ligated c-type ferric heme were employed to develop force-field parameters for molecular dynamics (MD) simulations of structural and dynamic features of these proteins. The new force-field parameters were applied to the alkaline form of yeast iso-1 cytochrome c to rationalize discrepancies resulting from distinct experimental conditions in prior structural studies and to provide insights into the mechanisms of the alkaline transition. Our simulations have revealed the dynamic nature of Ω-loop C in the Lys-ligated protein and its unfolding in the Lys-ligated conformer having this loop in the same position as in the native Met-ligated protein. The proximity of Tyr67 or Tyr74 to the Lys ligand of ferric heme iron suggests a possible mechanism of the backward alkaline transition where a proton donor Tyr assists in Lys dissociation. The developed force-field parameters will be useful in structural and dynamic characterization of other native or engineered Lys-ligated heme proteins.

Effect of proline content and histidine ligation on the dynamics of Ω-loop D and the peroxidase activity of iso-1-cytochrome c

To study how proline residues affect the dynamics of Ω-loop D (residues 70 to 85) of cytochrome c, we prepared G83P and G83A variants of yeast iso-1-cytochrome c (iso-1-Cytc) in the presence and absence of a K73H mutation. Ω-loop D is important in controlling both the electron transfer function of Cytc and the peroxidase activity of Cytc used in apoptosis because it provides the Met80 heme ligand. The G83P and G83A mutations have no effect on the global stability of iso-1-Cytc in presence or absence of the K73H mutation. However, both mutations destabilize the His73-mediated alkaline conformer relative to the native state. pH jump stopped-flow experiments show that the dynamics of the His73-mediated alkaline transition are significantly enhanced by the G83P mutation. Gated electron transfer studies show that the enhanced dynamics result from an increased rate of return to the native state, whereas the rate of loss of Met80 ligation is unchanged by the G83P mutation. Thus, the G83P substitution does not stiffen the conformation of the native state. Because bis-His heme ligation occurs when Cytc binds to cardiolipin-containing membranes, we studied the effect of His73 ligation on the peroxidase activity of Cytc, which acts as an early signal in apoptosis by causing oxygenation of cardiolipin. We find that the His73 alkaline conformer suppresses the peroxidase activity of Cytc. Thus, the bis-His ligated state of Cytc formed upon binding to cardiolipin is a negative effector for the peroxidase activity of Cytc early in apoptosis.

Insights into structure and activity relationship of clinically mutated PER1 and PER2 class A β-lactamase enzymes

PER1 and PER2 are among the class A β-lactamase enzymes, which have evolved clinically to form antibiotic resistance and have significantly expanded their spectrum of activity. Hence, there is a need to study the clinical mutation responsible for such β-lactamase mediated antibiotic resistance. Alterations in catalytic centre and Ω-loop structure could be the cause of antibiotic resistance in these β-lactamase enzymes. Structural and functional alterations are caused due to mutations on or near the catalytic centre, which results in active site plasticity and are responsible for its expanded spectrum of activity in these class A β-lactamase enzymes. Multiple sequence alignment, structure, kinetic, molecular docking, MMGBSA and molecular dynamic simulation comparisons were done on 38 clinically mutated and wild class A β-lactamase enzymes. This work shows that PER1 and PER2 enzymes contains most unique mutations and have altered Ω-loop structure, which could be responsible for altering the structure-activity relationship and extending the spectrum of activity of these enzymes. Alterations in molecular docking, MMGBSA, kinetic values reveals the modification in the binding and activity of these clinically mutated enzymes with antibiotics. Further, the cause of these alterations can be revealed by active site interactions and H-bonding pattern of these enzymes with antibiotics. Met69Gln, Glu104Thr, Tyr105Trp, Met129His, Pro167Ala, Glu168Gln, Asn170His, Ile173Asp and Asp176Gln mutations were uniquely found in PER1 and PER2 enzymes. These mutations occurs at catalytic important residues and results in altered interactions with β-lactam antibiotics. Hence, these mutations could be responsible for altering the structure-activity of PER1 and PER2 enzymes.

Influence of cysteine-directed mutations at the Ω-loops on peroxidase activity of human cytochrome c

Cytochrome c (Cytc) is a multifunctional protein associated with electron shuttling in the inner membrane of mitochondria and also involving in the apoptotic pathway. It has been identified that mutations located in the flexible central 40-57 Ω-loop including the naturally occurring G41S, Y48H, and A51V mutants, which are found in patients with thrombocytopenia 4, a platelet disorder, alter the structural properties of human Cytc (hCytc) that associated to enhanced peroxidase activity. In this work we compared the cysteine-directed mutants of hCytc located in three different parts of Ω-loops, i.e., T28C and G34C (proximal Ω-loop), and A50C (central Ω-loop), with respect to the wild-type (WT) hCytc. The mutants and WT hCytc were structurally characterized by circular dichroism, heating and chemical denaturations, and fluorescence spectroscopy. The flexibility at the cysteine mutated sites was directly determined by site-directed spin-labeling Electron Spin Resonance. Alkaline transitions were determined by pH titration and the alkaline conformers were related to peroxidase activity of all hCytc proteins. Structural and dynamic characterizations were rationally correlated to the modulation of peroxidase activity in these mutants in comparison to the WT hCytc. We found that the cysteine mutations at residues T28 and G34, both located in the same region of Ω-loop, developed different conformations and dynamical properties that lead to different effects on the rates of peroxidase activity (G34C was ~2.6 folds higher), whereas the rate of G34C was closer to that of A50C mutant. The results implied that the flexibility and local structures of the proximal Ω-loop could also play an important role in modulating the peroxidase activity which can be associated to apoptosis.

The K79G Mutation Reshapes the Heme Crevice and Alters Redox Properties of Cytochrome c

The two roles of cytochrome c (cyt c), in oxidative phosphorylation and apoptosis, critically depend on redox properties of its heme iron center. The K79G mutant has served as a parent protein for a series of mutants of yeast iso-1 cyt c. The mutation preserves the Met80 coordination to the heme iron, as found in WT* (K72A/C102S), and many spectroscopic properties of K79G and WT* are indistinguishable. The K79G mutation does not alter the global stability, fold, rate of Met80 dissociation, or thermodynamics of the alkaline transition (p K_a) of the protein. However, the reduction potential of the heme iron decreases; further, the p K_H of the trigger group and the rate of the Met-to-Lys ligand exchange associated with the alkaline transition decrease, suggesting changes in the environment of the heme. The rates of electron self-exchange and bimolecular electron transfer (ET) with positively charged inorganic complexes increase, as does the intrinsic peroxidase activity. Analysis of the reaction rates suggests that there is increased accessibility of the heme edge in K79G and supports the importance of the Lys79 site for bimolecular ET reactions of cyt c, including those with some of its native redox partners. Structural modeling rationalizes the observed effects to arise from changes in the volume of the heme pocket and solvent accessibility of the heme group. Kinetic and structural analyses of WT* characterize the properties of the heme crevice of this commonly employed reference variant. This study highlights the important role of Lys79 for defining functional redox properties of cyt c.

Probing Denatured State Conformational Bias in a Three-Helix Bundle, UBA(2), Using a Cytochrome c Fusion Protein

Previous work with the four-helix-bundle protein cytochrome c' from Rhodopseudomonas palustris using histidine-heme loop formation methods revealed fold-specific deviations from random coil behavior in its denatured state ensemble. To examine the generality of this finding, we extend this work to a three-helix-bundle polypeptide, the second ubiquitin-associated domain, UBA(2), of the human DNA excision repair protein. We use yeast iso-1-cytochrome c as a scaffold, fusing the UBA(2) domain at the N-terminus of iso-1-cytochrome c. We have engineered histidine into highly solvent accessible positions of UBA(2), creating six single histidine variants. Guanidine hydrochloride denaturation studies show that the UBA(2)-cytochrome c fusion protein unfolds in a three-state process with iso-1-cytochrome c unfolding first. Furthermore, engineered histidine residues in UBA(2) strongly destabilize the iso-1-cytochrome c domain. Equilibrium and kinetic histidine-heme loop formation measurements in the denatured state at 4 and 6 M guanidine hydrochloride show that loop stability decreases as the size of the histidine-heme loop increases, in accord with the Jacobson-Stockmayer equation. However, we observe that the His27-heme loop is both more stable than expected from the Jacobson-Stockmayer relationship and breaks more slowly than expected. These results show that the sequence near His27, which is in the reverse turn between helices 2 and 3 of UBA(2), is prone to persistent interactions in the denatured state. Therefore, consistent with our results for cytochrome c', this reverse turn sequence may help to establish the topology of this fold by biasing the conformational distribution of the denatured state.

Humanlike substitutions to Ω-loop D of yeast iso-1-cytochrome c only modestly affect dynamics and peroxidase activity

Structural studies of yeast iso-1-cytochrome c (L.J. McClelland, T.-C. Mou, M.E. Jeakins-Cooley, S.R. Sprang, B.E. Bowler, Proc. Natl. Acad. Sci. U.S.A. 111 (2014) 6648-6653) show that modest movement of Ω-loop D (residues 70-85, average RMSD versus the native structure: 0.81 Å) permits loss of Met80-heme ligation creating an available coordination site to catalyze the peroxidase activity mediated by cytochrome c early in apoptosis. However, Ala81 and Gly83 move significantly (RMSDs of 2.18 and 1.26 Å, respectively). Ala81 and Gly83 evolve to Ile and Val, respectively, in human cytochrome c and peroxidase activity decreases 25-fold relative to the yeast protein at pH 7. To test the hypothesis that these residues evolved to restrict the peroxidase activity of cytochrome c, A81I and G83V variants of yeast iso-1-cytochrome c were prepared. For both variants, the apparent pK_a of the alkaline transition increases by 0.2 to 0.3 relative to the wild type (WT) protein and the rate of opening the heme crevice is slowed. The cooperativity of acid unfolding is decreased for the G83V variant. At pH 7 and 8, the catalytic rate constant, k_cat, for the peroxidase activity of both variants decreases relative to WT, consistent with the effects on alkaline isomerization. Below pH 7, the loss in the cooperativity of acid unfolding causes k_cat for peroxidase activity to increase for the G83V variant relative to WT. Neither variant decreases k_cat to the level of the human protein, indicating that other residues also contribute to the low peroxidase activity of human cytochrome c.

Effect of an Imposed Contact on Secondary Structure in the Denatured State of Yeast Iso-1-cytochrome c

There is considerable evidence that long-range interactions stabilize residual protein structure under denaturing conditions. However, evaluation of the effect of a specific contact on structure in the denatured state has been difficult. Iso-1-cytochrome c variants with a Lys54 → His mutation form a particularly stable His-heme loop in the denatured state, suggestive of loop-induced residual structure. We have used multidimensional nuclear magnetic resonance methods to assign ¹H and ¹⁵N backbone amide and ¹³C backbone and side chain chemical shifts in the denatured state of iso-1-cytochrome c carrying the Lys54 → His mutation in 3 and 6 M guanidine hydrochloride and at both pH 6.4, where the His54-heme loop is formed, and pH 3.6, where the His54-heme loop is broken. Using the secondary structure propensity score, with the 6 M guanidine hydrochloride chemical shift data as a random coil reference state for data collected in 3 M guanidine hydrochloride, we found residual helical structure in the denatured state for the 60s helix and the C-terminal helix, but not in the N-terminal helix in the presence or absence of the His54-heme loop. Non-native helical structure is observed in two regions that form Ω-loops in the native state. There is more residual helical structure in the C-terminal helix at pH 6.4 when the loop is formed. Loop formation also appears to stabilize helical structure near His54, consistent with induction of helical structure observed when His-heme bonds form in heme-peptide model systems. The results are discussed in the context of the folding mechanism of cytochrome c.