Conformational flexibility decreased due to Y67F and F82H mutations in cytochrome c: Molecular dynamics simulation studies

Heme-iron ligand (M80-Fe) in cytochrome c is destabilizing: combined in vitro and in silico approaches to monitor changes in structure, stability and dynamics of the protein on mutation

Structure, stability and dynamics properties of horse cytochrome c (cyt c) and its genetically engineered M80G mutant have been investigated. The nature of the Met80 axial ligation to heme iron is believed to be the major determinant of the oxidation-reduction reactions inside and outside the cell of a particular cytochrome. This ligation has played an important role in the studies of protein structure, stability and protein folding/unfolding. To understand this ligation better, Met80 of horse cyt c has been mutated to Gly that is unable to bind to the heme iron. We have examined the effect of the M80G mutation on the structure and stability of the WT (wild type) protein by using absorbance spectroscopy, far-UV, near-UV and Soret circular dichroism, fluorescence spectroscopy and differential scanning calorimetry. We have observed that mutation caused a partial loss of secondary and tertiary structure with slightly increased overall stability of the protein. We have also measured the dynamic behavior of WT cyt c and its M80G mutant in the oxidized form (Fe³⁺) using the essential dynamics (ED) method. A 400 ns MD simulations were run for WT cyt c and its mutant M80G in water using GROMOS96 force field. MD results revealed that the stability and flexibility increased in mutant M80G (Fe…S (Met80) bond removed). Essential dynamics analysis revealed that the first five eigenvectors were mainly involved in overall motions of WT cyt c and its M80G mutant but the amplitude of concerted motions decreased in M80G mutant relative to WT cyt c.Communicated by Ramaswamy H. Sarma.

Met80 and Tyr67 affect the chemical unfolding of yeast cytochromec: comparing the solutionvs.immobilized state

The motional regime affects the unfolding propensity and axial heme coordination of the Met80Ala and Met80Ala/Tyr67Ala variants of yeast iso-1 cytochromec.

Temperature-dependent conformational dynamics of cytochrome c: Implications in apoptosis

Heat, electric shock, and burn injuries induce apoptosis by releasing cytochrome c (cyt-c) from mitochondria and by subsequently activating the death protease, caspases-3. During apoptosis, cyt-c undergoes changes in the secondary structure that have been suggested to increase its peroxidase activity. Information about these structural changes will provide better understanding of the apoptotic mechanism. Hence, temperature-dependent conformational dynamics of cyt-c has been investigated through molecular dynamics (MD) simulations to explain the structural changes and to correlate them with its apoptotic behavior. We observe that, at lower temperatures (223, 248, and 300K), the secondary structure of cyt-c, remains stable, while at higher temperatures (323, 373, 423, and 473K), the secondary structural regions change significantly. Further, our MD results indicate that these structural changes are mainly localized on α-helices, turns, β-sheets, and important loops that were involved in the stabilization of the heme conformation. This conformational transition between specific regions of secondary structure of cyt-c directly affects the electron tunneling properties of the proteins as observed experimentally. We quantify and compare these changes and explain that the temperature plays a vital role in assuring the structural stability of cyt-c and thus its functions. Our findings from this MD study reproduce experimental results at high temperatures and provide evidence for the alteration of the heme through the disruption of the H-bonding interactions between specific regions of cyt-c, thereby enhancing its peroxidase activity which plays a crucial role in the apoptotic process.

Investigation of structural dynamics of Thrombocytopenia Cargeeg mutants of human apoptotic cytochrome c: A molecular dynamics simulation approach

Naturally occurring mutations to cytochrome c (cyt-c) have been identified recently in patients with mild autosomal dominant thrombocytopenia (low platelet levels), which yield cyt-c mutants with enhanced apoptotic activity. However, the molecular mechanism underlying this low platelet production and enhanced apoptosis remain unclear. Therefore, an attempt is made herein for the first time to investigate the effects of mutations of glycine 41 by serine (G41S) and tyrosine 48 by histidine (Y48H) on the conformational and dynamic changes of apoptotic (Fe³⁺) cyt-c using all atom molecular dynamics (MD) simulations in explicit water solvent. Our 30ns MD simulations demonstrate considerable structural differences in G41S and Y48H compared to wild type (WT) cyt-c, such as increasing distances between the critical electron transfer residues results in open conformation at the heme active site, large fluctuations in β-turns and α-helices. Additionally, although the β-sheets remain mostly unaffected in all the three cyt-c simulations, the α-helices undergo conformational switch to β-turns in both the mutant simulations. Importantly, this conformational switch of α-helix to β-turn around heme active site should attributes to the loss of intraprotein H-bonds in the mutant simulations especially between NE2 (His26) and O (Pro44) in agreement with the experimental report. Further, essential dynamics analysis reveals that overall motions of WT cyt-c is mainly involved only in the first eigenvector, but in G41S and Y48H the overall motions are mainly in three and two eigenvectors respectively. Overall, the detailed atomistic level information provide a unifying description for the molecular mechanism of structural destabilization, disregulation of platelet formation and enhanced peroxidase activity of the mutant cyt-c's in the pathology of intrinsic apoptosis.

Topography of tyrosine residues and their involvement in peroxidation of polyunsaturated cardiolipin in cytochrome c/cardiolipin peroxidase complexes

Formation of cytochrome c (cyt c)/cardiolipin (CL) peroxidase complex selective toward peroxidation of polyunsaturated CLs is a pre-requisite for mitochondrial membrane permeabilization. Tyrosine residues - via the generation of tyrosyl radicals (Tyr) - are likely reactive intermediates of the peroxidase cycle leading to CL peroxidation. We used mutants of horse heart cyt c in which each of the four Tyr residues was substituted for Phe and assessed their contribution to the peroxidase catalysis. Tyr67Phe mutation was associated with a partial loss of the oxygenase function of the cyt c/CL complex and the lowest concentration of H(2)O(2)-induced Tyr radicals in electron paramagnetic resonance (EPR) spectra. Our MS experiments directly demonstrated decreased production of CL-hydroperoxides (CL-OOH) by Tyr67Phe mutant. Similarly, oxidation of a phenolic substrate, Amplex Red, was affected to a greater extent in Tyr67Phe than in three other mutants. Tyr67Phe mutant exerted high resistance to H(2)O(2)-induced oligomerization. Measurements of Tyr fluorescence, hetero-nuclear magnetic resonance (NMR) and computer simulations position Tyr67 in close proximity to the porphyrin ring heme iron and one of the two axial heme-iron ligand residues, Met80. Thus, the highly conserved Tyr67 is a likely electron-donor (radical acceptor) in the oxygenase half-reaction of the cyt c/CL peroxidase complex.

Binding conformation prediction between human acetylcholinesterase and cytochrome c using molecular modeling methods

The acetylcholinesterase (AChE) is important to terminate acetylcholine-mediated neurotransmission at cholinergic synapses. The pivotal role of AChE in apoptosome formation through the interactions with cytochrome c (Cyt c) was demonstrated in recent study. In order to investigate the proper binding conformation between the human AChE (hAChE) and human Cyt c (hCyt c), macro-molecular docking simulation was performed using DOT 2.0 program. The hCyt c was bound to peripheral anionic site (PAS) on hAChE and binding mode of the docked conformation was very similar to the reported crystal structure of the AChE and fasciculin-II (Fas-II) complex. Two 10ns molecular dynamics (MD) simulations were carried out to refine the binding mode of docked structure and to observe the differences of the binding conformations between the absent (Apo) and presence (Holo) of heme group. The key hydrogen bonding residues between hAChE and hCyt c proteins were found in Apo and Holo systems, as well as each Tyr341 and Trp286 residue of hAChE was participated in cation-pi (π) interactions with Lys79 of hCyt c in Apo and Holo systems, respectively. From the present study, although the final structures of the Apo and Holo systems have similar binding pattern, several differences were investigated in flexibilities, interface interactions, and interface accessible surface areas. Based on these results, we were able to predict the reasonable binding conformation which is indispensable for apoptosome formation.