Greek Key Topology Research Articles

How oxidized EGCG remodels α-synuclein fibrils into non-toxic aggregates: insights from computational simulations.

The misfolding and aggregation of the presynaptic protein α-synuclein (α-syn) is a pathological hallmark of Parkinson's disease (PD). Targeting α-syn has emerged as a promising therapeutic strategy for PD. Emerging in vitro evidence supports a dual action of epigallocatechin-3-gallate (EGCG) against amyloid neurotoxicity. EGCG can halt the formation of toxic aggregates by redirecting the amyloid fibril aggregation pathway toward non-toxic aggregates and remodeling the existing toxic fibrils into non-toxic aggregates. Moreover, EGCG oxidation can enhance fibril's remodeling by forming Schiff bases, leading to crosslinking of the fibril. However, this covalent modification is not required for amyloid remodeling, and establishing non-specific hydrophobic interactions with sidechains seems to be the main driver of amyloid remodeling by EGCG. Thioflavin (ThT) is a gold standard probe to detect amyloid fibrils in vitro, and oxidized EGCG competes with ThT for amyloid fibrils' binding sites. In this work, we performed docking and molecular dynamics (MD) simulations to gain insights into the intermolecular interactions of oxidized EGCG and ThT with a mature α-syn fibril. We find that oxidized EGCG moves within lysine-rich sites within the hydrophobic core of the α-syn fibril, forming aromatic and hydrogen-bonding (H-bond) interactions with different residues during the whole MD simulation time. In contrast, ThT, which does not remodel amyloid fibrils, was docked to the same sites but only via aromatic interactions. Our findings suggest that non-covalent interactions play a role in oxidized EGCG binding into the hydrophobic core, including H-bond and aromatic interactions with some residues in the amyloid remodeling processes. These interactions would ultimately lead to a disturbance of structural features as determinants for stabilizing this fibril into a compact and pathogenic Greek key topology.

Recent High-Resolution Structures of Amyloids Involved in Neurodegenerative Diseases.

Amyloids are highly ordered aggregates composed of proteins or peptides. They are involved in several pathologies, including hallmark neurodegenerative disorders such as Alzheimer’s (AD) and Parkinson’s (PD). Individuals affected by these diseases accumulate in their brains amyloids inclusions composed of misfolded forms of a peptide (Aβ) and a protein (Tau) in AD and α-synuclein protein (α-Sn) in PD. Tau and α-Sn aggregates are also present in other neurodegenerative diseases. The insoluble nature and heterogeneity of amyloids have hampered their study at the molecular level. However, the use of solid state NMR and Cryogenic-electron microscopy along with fine-tuned modulation of the aggregation in vitro and improved isolation methods of brain-derived amyloids has allowed the elucidation of these elusive conformations at high resolution. In this work, we review the latest progress on the recent amyloid structures reported for Aβ, Tau, and α-Sn. The two-fold symmetry emerges as a convergent feature in the tridimensional arrangement of the protofilaments in the fibrillary structure of these pathological amyloids, with many of them exhibiting a Greek-key topology as part of their overall architecture. These specific features can serve as novel guides to seek potential molecular targets in drug design efforts.

Conformational distortion in a fibril-forming oligomer arrests alpha-Synuclein fibrillation and minimizes its toxic effects

The fibrillation pathway of alpha-Synuclein, the causative protein of Parkinson’s disease, encompasses transient, heterogeneous oligomeric forms whose structural understanding and link to toxicity are not yet understood. We report that the addition of the physiologically-available small molecule heme at a sub-stoichiometric ratio to either monomeric or aggregated α-Syn, targets a His50 residue critical for fibril-formation and stabilizes the structurally-heterogeneous populations of aggregates into a minimally-toxic oligomeric state. Cryo-EM 3D reconstruction revealed a ‘mace’-shaped structure of this monodisperse population of oligomers, which is comparable to a solid-state NMR Greek key-like motif (where the core residues are arranged in parallel in-register sheets with a Greek key topology at the C terminus) that forms the fundamental unit/kernel of protofilaments. Further structural analyses suggest that heme binding induces a distortion in the Greek key-like architecture of the mace oligomers, which impairs their further appending into protofilaments and fibrils. Additionally, our study reports a novel mechanism of prevention as well as reclamation of amyloid fibril formation by blocking an inter-protofilament His50 residue using a small molecule.

Network Connectivity, Centrality and Fragmentation in the Greek-Key Protein Topology.

Understanding and computationally predicting the protein folding process remains one of the most challenging scientific problems and has uniquely garnered the interdisciplinary efforts of researchers from both the biological, chemical, physical and computational disciplines. Previous studies have demonstrated the importance of long-range interactions in guiding the native structure. However, predicting how the native long-range interaction network forms to generate a specific topology from among all other conformations remains unresolved. The present research study conducts an exploratory study to identify amino acids and long-range interactions that have the potential to play a key role in building and maintaining the protein topology. Towards this end, the application of network science is utilized and developed to analyze the structures of a group of proteins that share a common Greek-key topology but differ in sequence, secondary structure and function. We investigate the idea that the residues with high betweeness centrality score are potentially significant in maintaining the protein network and in governing the Greek-key topology. This hypothesis is tested by two different computational methods: through a fragmentation test and by the analysis of diameter impacts. In summary, we find a subset of selected residues in similar geographical positions in all model proteins, which demonstrates the role of these specific residues and regions in governing the Greek-key topology from a network perspective.

Solid-state NMR structure of a pathogenic fibril of full-length human α-synuclein.

Misfolded α-synuclein amyloid fibrils are the principal components of Lewy bodies and neurites, hallmarks of Parkinson’s disease (PD). Here we present a high-resolution structure of an α-synuclein fibril, in a form that induces robust pathology in primary neuronal culture, determined by solid-state NMR spectroscopy and validated by electron microscopy and X-ray fiber diffraction. Over 200 unique long-range distance restraints define a consensus structure with common amyloid features including parallel in-register β-sheets and hydrophobic core residues, but also substantial complexity, arising from diverse structural features: an intermolecular salt bridge, a glutamine ladder, close backbone interactions involving small residues, and several steric zippers stabilizing a novel, orthogonal Greek-key topology. These characteristics contribute to the robust propagation of this fibril form, as evidenced by structural similarity of early-onset PD mutants. The structure provides a framework for understanding the interactions of α-synuclein with other proteins and small molecules to diagnose and treat PD.

Biophysical Analysis of the Transition of an all α-Helical Greek-Key Protein into Amyloid Fibrils Composed of β-Sheet Structure

Amyloidosis resulting from the deposition of aggregated protein has been linked to many debilitating degenerative diseases which include most notably Alzheimer's and Parkinson's. The tendency for a protein to alternatively form highly ordered amyloid fibrils is dependent on many biological factors. Mutations, temperature, concentration, translational motion and pH play a pivotal role in inducing fibril aggregate assembly in vitro. The key feature appears to be the need to destabilize the native state structure as a required first step. In this paper we report on the detailed conversion of the death domain of the human Fas-associated death domain, an all α-helical protein with a Greek-key topology, into an all β-sheet amyloid fibril, using a comprehensive range of spectroscopic techniques that provide insight into this process. This transition from α-helical to β-sheet seems to require destabilization but not complete loss of the secondary structure to explore alternative conformations. This is a fascinating transition that supports the hypothesis that all proteins have the innate ability to form a fibril-like structure. Thus, the primary structure can encode two alternative three-dimensional structures: the native, functional state and the β-amyloid state. The Fas-associated death domain does not appear to naturally form amyloid fibrils in vivo. Our results clearly indicate that proteins evolved to avoid amyloid fibril formation because we find that the conditions required for formation in our model system are very specific and far from physiological.

Folding of an all-helical Greek-key protein monitored by quenched-flow hydrogen–deuterium exchange and NMR spectroscopy

To advance our understanding of the protein folding process, we use stopped-flow far-ultraviolet (far-UV) circular dichroism and quenched-flow hydrogen-deuterium exchange coupled with nuclear magnetic resonance (NMR) spectroscopy to monitor the formation of hydrogen-bonded secondary structure in the C-terminal domain of the Fas-associated death domain (Fadd-DD). The death domain superfamily fold consists of six α-helices arranged in a Greek-key topology, which is shared by the all-β-sheet immunoglobulin and mixed α/β-plait superfamilies. Fadd-DD is selected as our model death domain protein system because the structure of this protein has been solved by NMR spectroscopy, and both thermodynamic and kinetic analysis indicate it to be a stable, monomeric protein with a rapidly formed hydrophobic core. Stopped-flow far-UV circular dichroism spectroscopy revealed that the folding process was monophasic and the rate is 23.4s(-1). Twenty-two amide hydrogens in the backbone of the helices and two in the backbone of the loops were monitored, and the folding of all six helices was determined to be monophasic with rate constants between 19 and 22s(-1). These results indicate that the formation of secondary structure is largely cooperative and concomitant with the hydrophobic collapse. This study also provides unprecedented insight into the formation of secondary structure within the highly populated Greek-key fold more generally. Additional insights are gained by calculating the exchange rates of 23 residues from equilibrium hydrogen-deuterium exchange experiments. The majority of protected amide protons are found on helices2, 4, and 5, which make up core structural elements of the Greek-key topology.

Genomic Evidence for Elements of a Programmed Cell Death Pathway In Plasmodium: Exploiting Programmed Parasite Death for Malaria Control?

AbstractAbstract 4226Malaria, caused by species of the Plasmodium parasite, remains a major global health burden. To combat the disease, new areas of investigation are constantly developed with a view to understanding the biology of the organism and identifying potential drug targets. A new line of enquiry has emerged concerning the occurrence of parasite programmed cell death (PCD); however, there are conflicting experimental reports regarding the presence and phenotype of PCD in Plasmodium species. Furthermore, very little genomic evidence exists for a PCD pathway in this genus with only caspase-like domains being detected. One reason for the limited availability of genomic data may be the unusual characteristics of the P. falciparum genome, such as the extreme AT bias. In this study we examined the complete genome sequences of four Plasmodium species (P. falciparum, P. vivax, P. yoelii and P. knowlesi) for evidence of a p53-dependant PCD pathway. This pathway is well characterized in animals, and elements of the molecular machinery have been identified in protozoa. In addition, p53-like responses leading to apoptosis (a common PCD phenotype) have been demonstrated in phylogenetically diverse eukaryotes, including animals, protozoa, green algae and plants, suggesting that the origin of p53-initiated PCD is evolutionary ancient leading us to believe that a similar genomic toolkit may exist in Plasmodium. A detailed analysis of the four Plasmodium genomes was performed using an array of computational approaches, which included standard homology methods, phylogenetics, structural models and a novel evolutionary rate-based alignment algorithm FIRE (Functional Inference using the Rates of Evolution), which was developed to identify homologous and analogous genes in organisms with unusual genomes, such as P. falciparum, and hence low sequence similarity. Homology methods uncovered key elements of a classical PCD pathway, including ATM, MDM2, CR6 and three peptidase C14 (catalytic caspase) domains in each of the four Plasmodium genomes. Phylogenetic analysis demonstrated that the peptidase C14 (caspase) domains are evolutionary ancient and cluster with other PCD-related caspases suggesting that they are involved in a PCD pathway as opposed to other cellular functions that may use similar domains. In addition, highly sensitive hidden Markov models retrieved 15 sequences with low similarity to the eukaryotic p53 DNA-binding domain (DBD). Further analysis revealed that two P. falciparum sequences (accession numbers PFE1120w and PF11_0091 in the Plasmodium database www.PlasmoDB.org version 6.5) may be p53 DBD-like sequences. Both have similar evolutionary rates across codons, as well as anti-parallel beta sheet structures with Greek key topology, which are a characteristic feature of known p53 DBDs. Whether either of these proteins has functional significance for the PCD pathway in P. falciparum requires empirical verification. However, these data and the existence of p53-like activity in plants and algae, suggest that this pathway has an ancient origin and provide the first evidence for a p53-dependant PCD pathway in malaria parasites. The presence of this genomic toolkit raises questions regarding the potential role of this pathway in the P. falciparum lifecycle. The phenomenon of programmed parasite death may have adaptive significance and proteins in the PCD pathway may potentially be sufficiently different from the human orthologs, given the low level of sequence similarity. These findings could thus open a new line of investigation for novel malaria drug targets. Disclosures:No relevant conflicts of interest to declare.

Analysis of conservation in the Fas-associated death domain protein and the importance of conserved tryptophans in structure, stability and folding

Computational and experimental studies focusing on the role of conserved amino acids for folding and stability is an active and promising area of research. To date however, only a small fraction of the potential superfamilies have been investigated. To further expand our understanding we present the results of a bioinformatics analysis of the death domain superfamily. The fold consists of a six helical bundle with a Greek-key topology. Our sequence and structural studies have identified a group of conserved hydrophobic residues and corresponding long-range interactions, which we propose are important in the formation and stabilization of the hydrophobic core and native topology. Six conserved hydrophilic residues were additionally identified and may play a functional role during apoptosis. We also report the establishment of an experimental system that will facilitate studies to test the role of the conserved residues in folding and stability. Equilibrium unfolding and refolding studies of a model superfamily member, Fas-associated death domain protein indicate that the process is cooperative, two-state and reversible. Stopped-flow fluorescence studies reveal that the folding is rapid and biphasic with the majority of the hydrophobic core forming in the first phase. Site-directed mutagenesis studies indicate that conserved tryptophans 112 and 148 are important to structure, native state stability and folding. These results present the earliest conservation analysis and biophysical characterization of the Fas-associated death domain.

Crystal structure of human IPS-1/MAVS/VISA/Cardif caspase activation recruitment domain.

BackgroundIPS-1/MAVS/VISA/Cardif is an adaptor protein that plays a crucial role in the induction of interferons in response to viral infection. In the initial stage of the intracellular antiviral response two RNA helicases, retinoic acid inducible gene-I (RIG-I) and melanoma differentiation-association gene 5 (MDA5), are independently able to bind viral RNA in the cytoplasm. The 62 kDa protein IPS-1/MAVS/VISA/Cardif contains an N-terminal caspase activation and recruitment (CARD) domain that associates with the CARD regions of RIG-I and MDA5, ultimately leading to the induction of type I interferons. As a first step towards understanding the molecular basis of this important adaptor protein we have undertaken structural studies of the IPS-1 MAVS/VISA/Cardif CARD region.ResultsThe crystal structure of human IPS-1/MAVS/VISA/Cardif CARD has been determined to 2.1Å resolution. The protein was expressed and crystallized as a maltose-binding protein (MBP) fusion protein. The MBP and IPS-1 components each form a distinct domain within the structure. IPS-1/MAVS/VISA/Cardif CARD adopts a characteristic six-helix bundle with a Greek-key topology and, in common with a number of other known CARD structures, contains two major polar surfaces on opposite sides of the molecule. One face has a surface-exposed, disordered tryptophan residue that may explain the poor solubility of untagged expression constructs.ConclusionThe IPS-1/MAVS/VISA/Cardif CARD domain adopts the classic CARD fold with an asymmetric surface charge distribution that is typical of CARD domains involved in homotypic protein-protein interactions. The location of the two polar areas on IPS-1/MAVS/VISA/Cardif CARD suggest possible types of associations that this domain makes with the two CARD domains of MDA5 or RIG-I. The N-terminal CARD domains of RIG-I and MDA5 share greatest sequence similarity with IPS-1/MAVS/VISA/Cardif CARD and this has allowed modelling of their structures. These models show a very different charge profile for the equivalent surfaces compared to IPS-1/MAVS/VISA/Cardif CARD.

Three-dimensional Structure of the EphB2 Receptor in Complex with an Antagonistic Peptide Reveals a Novel Mode of Inhibition

The Eph family of receptor tyrosine kinases has been implicated in tumorigenesis as well as pathological forms of angiogenesis. Understanding how to modulate the interaction of Eph receptors with their ephrin ligands is therefore of critical interest for the development of therapeutics to treat cancer. Previous work identified a set of 12-mer peptides that displayed moderate binding affinity but high selectivity for the EphB2 receptor. The SNEW antagonistic peptide inhibited the interaction of EphB2 with ephrinB2, with an IC50 of approximately 15 microm. To gain a better molecular understanding of how to inhibit Eph/ephrin binding, we determined the crystal structure of the EphB2 receptor in complex with the SNEW peptide to 2.3-A resolution. The peptide binds in the hydrophobic ligand-binding cleft of the EphB2 receptor, thus competing with the ephrin ligand for receptor binding. However, the binding interactions of the SNEW peptide are markedly different from those described for the TNYL-RAW peptide, which binds to the ligand-binding cleft of EphB4, indicating a novel mode of antagonism. Nevertheless, we identified a conserved structural motif present in all known receptor/ligand interfaces, which may serve as a scaffold for the development of therapeutic leads. The EphB2-SNEW complex crystallized as a homodimer, and the residues involved in the dimerization interface are similar to those implicated in mediating tetramerization of EphB2-ephrinB2 complexes. The structure of EphB2 in complex with the SNEW peptide reveals novel binding determinants that could serve as starting points in the development of compounds that modulate Eph receptor/ephrin interactions and biological activities.

Plasticity Within the Obligatory Folding Nucleus of an Immunoglobulin-like Domain

A number of β-sandwich immunoglobulin-like domains have been shown to fold using a set of structurally equivalent residues that form a folding nucleus deep within the core of the protein. Formation of this nucleus is sufficient to establish the complex Greek key topology of the native state. These nucleating residues are highly conserved within the immunoglobulin superfamily, but are less well conserved in the fibronectin type III (fnIII) superfamily, where the requirement is simply to have four interacting hydrophobic residues. However, there are rare examples where this nucleation pattern is absent. In this study, we have investigated the folding of a novel member of the fnIII superfamily whose nucleus appears to lack one of the four buried hydrophobic residues. We show that the folding mechanism is unaltered, but the folding nucleus has moved within the hydrophobic core.

Rapid Folding and Unfolding of Apaf-1 CARD

Caspase recruitment domains (CARDs) are members of the death domain superfamily and contain six antiparallel helices in an α-helical Greek key topology. We have examined the equilibrium and kinetic folding of the CARD of Apaf-1 (apoptotic protease activating factor 1), which consists of 97 amino acid residues, at pH 6 and pH 8. The results showed that an apparent two state equilibrium mechanism is not adequate to describe the folding of Apaf-1 CARD at either pH, suggesting the presence of intermediates in equilibrium unfolding. Interestingly, the results showed that the secondary structure is less stable than the tertiary structure, based on the transition mid-points for unfolding. Single mixing and sequential mixing stopped-flow studies showed that Apaf-1 CARD folds and unfolds rapidly and suggest a folding mechanism that contains parallel channels with two unfolded conformations folding to the native conformation. Kinetic simulations show that a slow folding phase is described by a third conformation in the unfolded ensemble that interconverts with one or both unfolded species. Overall, the native ensemble is formed rapidly upon refolding. This is in contrast to other CARDs in which folding appears to be dominated by formation of kinetic traps.

Structural and Biophysical Characterization of the EphB4·EphrinB2 Protein-Protein Interaction and Receptor Specificity

Increasing evidence implicates the interaction of the EphB4 receptor with its preferred ligand, ephrinB2, in pathological forms of angiogenesis and in tumorigenesis. To identify the molecular determinants of the unique specificity of EphB4 for ephrinB2, we determined the crystal structure of the ligand binding domain of EphB4 in complex with the extracellular domain of ephrinB2. This structural analysis suggested that one amino acid, Leu-95, plays a particularly important role in defining the structural features that confer the ligand selectivity of EphB4. Indeed, all other Eph receptors, which promiscuously bind many ephrins, have a conserved arginine at the position corresponding to Leu-95 of EphB4. We have also found that amino acid changes in the EphB4 ligand binding cavity, designed based on comparison with the crystal structure of the more promiscuous EphB2 receptor, yield EphB4 variants with altered binding affinity for ephrinB2 and an antagonistic peptide. Isothermal titration calorimetry experiments with an EphB4 Leu-95 to arginine mutant confirmed the importance of this amino acid in conferring high affinity binding to both ephrinB2 and the antagonistic peptide ligand. Isothermal titration calorimetry measurements also revealed an interesting thermodynamic discrepancy between ephrinB2 binding, which is an entropically driven process, and peptide binding, which is an enthalpically driven process. These results provide critical information on the EphB4*ephrinB2 protein interfaces and their mode of interaction, which will facilitate development of small molecule compounds inhibiting the EphB4*ephrinB2 interaction as novel cancer therapeutics.

Elucidation of conserved long-range interaction networks in proteins and their significance in determining protein topology

Investigations into the nature of sequence and structural conservation underlying protein folds have recently yielded profound insights into the mechanism of protein folding and stability. Combining this avenue of research with principles being pioneered in the field of network science holds the promise to further extend the boundaries of our knowledge. In this report we propose that critical determinants of a protein's native topology are encoded by a conserved network of interactions between amino acids from geographically important positions. This hypothesis is based on the novel elucidation of a conserved network of long-range interactions within a set of proteins that share a Greek-key topology and similar chain length, but differ in secondary structure composition, function and sequence. Exploratory macromolecular simulations using the conserved networks as constraints were successful in generating the gross native-like topology from a random linear coil for each of the model proteins. The results indicate that the conserved network contains governing features and supports the idea that the geographical location of these residue interactions is a pivotal feature underlying their conservation. The partially folded model proteins also display a clear scale-free distribution of long-range interactions. To further test the hypothesis, the network parameter betweeness-centrality was calculated for the protein structure networks of our model proteins and highlights two structural elements as particularly vital to the structural stability of the network topology.

Solution structure of ribosomal protein S28E from Methanobacterium thermoautotrophicum.

The ribosomal protein S28E from the archaeon Methanobacterium thermoautotrophicum is a component of the 30S ribosomal subunit. Sequence homologs of S28E are found only in archaea and eukaryotes. Here we report the three-dimensional solution structure of S28E by NMR spectroscopy. S28E contains a globular region and a long C-terminal tail protruding from the core. The globular region consists of four antiparallel beta-strands that are arranged in a Greek-key topology. Unique features of S28E include an extended loop L2-3 that folds back onto the protein and a 12-residue charged C-terminal tail with no regular secondary structure and greater flexibility relative to the rest of the protein. The structural and surface resemblance to OB-fold family of proteins and the presence of highly conserved basic residues suggest that S28E may bind to RNA. A broad positively charged surface extending over one side of the beta-barrel and into the flexible C terminus may present a putative binding site for RNA.

Solution NMR structure of the 30S ribosomal protein S28E from Pyrococcus horikoshii.

We report NMR assignments and solution structure of the 71-residue 30S ribosomal protein S28E from the archaean Pyrococcus horikoshii, target JR19 of the Northeast Structural Genomics Consortium. The structure, determined rapidly with the aid of automated backbone resonance assignment (AutoAssign) and automated structure determination (AutoStructure) software, is characterized by a four-stranded beta-sheet with a classic Greek-key topology and an oligonucleotide/oligosaccharide beta-barrel (OB) fold. The electrostatic surface of S28E exhibits positive and negative patches on opposite sides, the former constituting a putative binding site for RNA. The 13 C-terminal residues of the protein contain a consensus sequence motif constituting the signature of the S28E protein family. Surprisingly, this C-terminal segment is unstructured in solution.

Approaching the speed limit for Greek Key β-barrel formation: transition-state movement tunes folding rate of zinc-substituted azurin

Azurin is a blue-copper protein with a β-barrel structure of Greek Key topology. In vitro, copper can be substituted with zinc without change in protein structure. We here analyze the kinetic folding behavior of zinc-substituted Pseudomonas aeruginosa azurin. Our findings can be summarized in three key conclusions: First, zinc remains strongly bound to the polypeptide upon unfolding, suggesting that the cofactor may bind to the protein before polypeptide folding in vivo. Second, the semi-logarithmic plot of folding and unfolding rates for zinc-substituted azurin as a function of denaturant concentration exhibits curvature due to a changing transition-state structure. Third, the extrapolated folding speed in water for zinc-substituted azurin is similar to that of other proteins with the same topology, implying that there is a speed limit that can be modulated by stability-driven transition-state movement for formation of β-barrel structures with Greek Key topology.

A simple lattice model that exhibits a protein-like cooperative all-or-none folding transition.

In a recent paper (D. Gront et al., Journal of Chemical Physics, Vol. 115, pp. 1569, 2001) we applied a simple combination of the Replica Exchange Monte Carlo and the Histogram methods in the computational studies of a simplified protein lattice model containing hydrophobic and polar units and sequence-dependent local stiffness. A well-defined, relatively complex Greek-key topology, ground (native) conformations was found; however, the cooperativity of the folding transition was very low. Here we describe a modified minimal model of the same Greek-key motif for which the folding transition is very cooperative and has all the features of the "all-or-none" transition typical of real globular proteins. It is demonstrated that the all-or-none transition arises from the interplay between local stiffness and properly defined tertiary interactions. The tertiary interactions are directional, mimicking the packing preferences seen in proteins. The model properties are compared with other minimal protein-like models, and we argue that the model presented here captures essential physics of protein folding (structurally well-defined protein-like native conformation and cooperative all-or-none folding transition).

Effects of metal ligation and oxygen on the reversibility of the thermal denaturation of Pseudomonas aeruginosa azurin.

Thermodynamic equilibrium transition models in DSC are only applicable to reversible processes, but reversibility of the thermal transitions of proteins is comparatively rare because of intermolecular aggregation of denatured proteins and the degradation that occurs at high temperatures. The cupredoxin azurin from Pseudomonas aeruginosa has previously been found to exhibit irreversible thermal denaturation, both as holo- and apoprotein [Engeseth, H. R., and McMillin, D. R. (1986) Biochemistry 25, 2448-2455]. In this study, however, we demonstrate that this beta-barrel protein of Greek key topology in fact unfolds reversibly in anaerobic solutions when nonreducible metal ions are ligated to the protein. We show that it is the metal-coordinating cysteine residue (C112) that becomes exclusively oxidized in a transition metal catalyzed oxidation reaction with dissolved O(2) at high temperatures. Both Cu(I)- and Zn(II)-coordinating wild-type azurin therefore unfold reversibly in anaerobic solutions, as well as the Zn(II)-coordinating disulfide-deficient C3A/C26A mutant. Correspondingly, apoazurin mutants C112A and C112S unfold reversibly, even in aerobic solutions, and exhibit nearly perfect two-state transitions. Unfolding of Cu(II)-coordinating azurin is, on the other hand, always irreversible due to autoxidation of the thiolate resulting in Cu(I) and a thiyl radical prone to oxidation.