Beyond the coupled distortion model: structural analysis of the single domain cupredoxin AcoP, a green mononuclear copper centre with original features.

Cupredoxins are widely occurring copper-binding proteins with a typical Greek-key beta barrel fold. They are generally described as electron carriers that rely on a T1 copper centre coordinated by four ligands provided by the folded polypeptide. The discovery of novel cupredoxins demonstrates the high diversity of this family, with variations in terms of copper-binding ligands, copper centre geometry, redox potential, as well as biological function. AcoP is a periplasmic cupredoxin belonging to the iron respiratory chain of the acidophilic bacterium Acidithiobacillus ferrooxidans. AcoP presents original features, including high resistance to acidic pH and a constrained green-type copper centre of high redox potential. To understand the unique properties of AcoP, we undertook structural and biophysical characterization of wild-type AcoP and of two Cu-ligand mutants (H166A and M171A). The crystallographic structures, including native reduced AcoP at 1.65 Å resolution, unveil a typical cupredoxin fold. The presence of extended loops, never observed in previously characterized cupredoxins, might account for the interaction of AcoP with physiological partners. The Cu-ligand distances, determined by both X-ray diffraction and EXAFS, show that the AcoP metal centre seems to present both T1 and T1.5 features, in turn suggesting that AcoP might not fit well to the coupled distortion model. The crystal structures of two AcoP mutants confirm that the active centre of AcoP is highly constrained. Comparative analysis with other cupredoxins of known structures, suggests that in AcoP the second coordination sphere might be an important determinant of active centre rigidity due to the presence of an extensive hydrogen bond network. Finally, we show that other cupredoxins do not perfectly follow the coupled distortion model as well, raising the suspicion that further alternative models to describe copper centre geometries need to be developed, while the importance of rack-induced contributions should not be underestimated.

Systematic elucidation of the second coordination sphere effect on the structure and properties of a blue copper protein, pseudoazurin.

The rational structural and computational studies of a blue copper protein, pseudoazurin (PAz), and its Met16X (X = Phe, Leu, Val, Ile) variants gave clear functional meanings of the noncovalent interaction (NCI) through the second coordination sphere. The high-resolution X-ray crystal structures of Met16X PAz demonstrated that the active site geometry is significantly affected by the substitution of Met16, which is located within the NCI distance from the His81 imidazole ring at the copper active site. The computational chemistry calculations based on the crystal structure analyses confirmed that the NCI of S-π/CH-π (wild-type), π-π (Met16Phe), double CH-π (Met16Leu), and single CH-π (Met16Val and Met16Ile). The estimated interaction energies for the NCI demonstrated that the fine-tuning of the protein stability and Cu site properties form the second coordination sphere of PAz.

Applicability of perturbed matrix method for charge transfer studies at bio/metallic interfaces: a case of azurin.

As the field of nanoelectronics based on biomolecules such as peptides and proteins rapidly grows, there is a need for robust computational methods able to reliably predict charge transfer properties at bio/metallic interfaces. Traditionally, hybrid quantum-mechanical/molecular-mechanical techniques are employed for systems where the electron hopping transfer mechanism is applicable to determine physical parameters controlling the thermodynamics and kinetics of charge transfer processes. However, these approaches are limited by a relatively high computational cost when extensive sampling of a configurational space is required, like in the case of soft biomatter. For these applications, semi-empirical approaches such as the perturbed matrix method (PMM) have been developed and successfully used to study charge-transfer processes in biomolecules. Here, we explore the performance of PMM on prototypical redox-active protein azurin in various environments, from solution to vacuum interfaces with gold surfaces and protein junction. We systematically benchmarked the robustness and convergence of the method with respect to the quantum-centre size, size of the Hamiltonian, number of samples, and level of theory. We show that PMM can adequately capture all the trends associated with the structural and electronic changes related to azurin oxidation at bio/metallic interfaces.

Structural basis of bacterial effector protein azurin targeting tumor suppressor p53 and inhibiting its ubiquitination.

Tumor suppressor p53 prevents tumorigenesis by promoting cell cycle arrest and apoptosis through transcriptional regulation. Dysfunction of p53 occurs frequently in human cancers. Thus, p53 becomes one of the most promising targets for anticancer treatment. A bacterial effector protein azurin triggers tumor suppression by stabilizing p53 and elevating its basal level. However, the structural and mechanistic basis of azurin-mediated tumor suppression remains elusive. Here we report the atomic details of azurin-mediated p53 stabilization by combining X-ray crystallography with nuclear magnetic resonance. Structural and mutagenic analysis reveals that the p28 region of azurin, which corresponds to a therapeutic peptide, significantly contributes to p53 binding. This binding stabilizes p53 by disrupting COP1-mediated p53 ubiquitination and degradation. Using the structure-based design, we obtain several affinity-enhancing mutants that enable amplifying the effect of azurin-induced apoptosis. Our findings highlight how the structure of the azurin-p53 complex can be leveraged to design azurin derivatives for cancer therapy.

Long-Range Charge Delocalization Mediates the Ultrafast Ligand-to-Metal Charge Transfer Dynamics at the Cu2+-Active Site in Azurin.

The blue color in metalloprotein azurin has traditionally been attributed to the intense cysteine-to-Cu2+ ligand-to-metal charge transfer transition centered at 628 nm. Although resonance Raman measurements of the Cu2+ active site have implied that the LMCT transition electronically couples to the protein scaffold well beyond its primary metal-ligand coordination shell, the structural extent of this electronic coupling and visualization of the protein-mediated charge transfer dynamics have remained elusive. Here, using femtosecond broadband transient absorption and impulsive Raman spectroscopy, we provide direct evidence for a rapid relaxation between two distinct charge transfer states, having different spatial delocalization, within ∼300 fs followed by recombination of charges in subpicosecond time scales. We invoke the formation of a protein-centered radical cation, possibly Trp48 or a Phe residue, within 100 fs substantiating the long-range electronic coupling for the first time beyond the traditional copper active site. The Raman spectra of the excited CT state show the presence of protein-centric vibrations along with the vibrational modes assigned to the copper active site. Our results demonstrate a large delocalization length scale of the initially populated CT state, thereby highlighting the possibility of exploiting azurin photochemistry for energy conversion techniques.

Thioester synthesis by a designed nickel enzyme models prebiotic energy conversion.

The formation of carbon-carbon bonds from prebiotic precursors such as carbon dioxide represents the foundation of all primordial life processes. In extant organisms, this reaction is carried out by the carbon monoxide dehydrogenase (CODH)/acetyl coenzyme A synthase (ACS) enzyme, which performs the cornerstone reaction in the ancient Wood-Ljungdahl metabolic pathway to synthesize the key biological metabolite, acetyl-CoA. Despite its significance, a fundamental understanding of this transformation is lacking, hampering efforts to harness analogous chemistry. To address these knowledge gaps, we have designed an artificial metalloenzyme within the azurin protein scaffold as a structural, functional, and mechanistic model of ACS. We demonstrate the intermediacy of the NiI species and requirement for ordered substrate binding in the bioorganometallic carbon-carbon bond-forming reaction from the one-carbon ACS substrates. The electronic and geometric structures of the nickel-acetyl intermediate have been characterized using time-resolved optical, electron paramagnetic resonance, and X-ray absorption spectroscopy in conjunction with quantum chemical calculations. Moreover, we demonstrate that the nickel-acetyl species is chemically competent for selective acyl transfer upon thiol addition to biosynthesize an activated thioester. Drawing an analogy to the native enzyme, a mechanism for thioester generation by this ACS model has been proposed. The fundamental insight into the enzymatic process provided by this rudimentary ACS model has implications for the evolution of primitive ACS-like proteins. Ultimately, these findings offer strategies for development of highly active catalysts for sustainable generation of liquid fuels from one-carbon substrates, with potential for broad applications across diverse fields ranging from energy storage to environmental remediation.

Tuning reduction potentials of type 1 copper center in azurin by replacing a histidine ligand with its isostructural analogues.

Type 1 copper proteins have a conserved ligand set of one cysteine and two histidines, with many proteins, such as azurin, also containing an axial methionine. While the cysteine and methionine in azurin have been replaced with their respective isostructural analogues of unnatural amino acids to reveal their roles in tuning electronic structures and functional properties, such as reduction potentials (E°'), the histidine ligands have not been probed in this way. We herein report the substitution of His117 in azurin with three unnatural isostructural analogues, 5-nitrohistidine(Ntr), thiazolylalanine(SHis) and 1-methylhistidine(MeH) by expressed protein ligation. While UV-vis absorption and electron paramagnetic resonance spectroscopies confirm that isostructural replacement results in minimal structural change in the Cu(II) state, the E°' of these variants increases with increasing pKa of the δ nitrogens of the imidazole. This counter-intuitive relationship between E°' of the protein and pKa of the sidechain group suggests additional factors may play a role in tuning E°'.

Biomolecular Electron Controller Composed of Nanobiohybrid with Electrically Released Complex for Spatiotemporal Control of Neuronal Differentiation.

In vitro spatiotemporal control of cell differentiation is a critical issue in several biomedical fields such as stem cell therapy and regenerative medicine, as it enables the generation of heterogeneous tissue structures similar to those of their native counterparts. However, the simultaneous control of both spatial and temporal cell differentiation poses important challenges, and therefore no previous studies have achieved this goal. Here, the authors develop a cell differentiation biomolecular electron controller ("Biomoletron") composed of recombinant proteins, DNA, Au nanoparticles, peptides, and an electrically released complex with retinoic acid (RA) to spatiotemporally control SH-SY5Y cell differentiation. RA is only released from the Biomoletron when the complex is electrically stimulated, thus demonstrating the temporal control of SH-SY5Y cell differentiation. Furthermore, by introducing a patterned Au substrate that allows controlling the area where the Biomoletron is immobilized, spatiotemporal differentiation of the SH-SY5Y cell is successfully achieved. Therefore, the proposed Biomoletron-mediated differentiation method provides a promising strategy for spatiotemporal cell differentiation control with applications in regenerative medicine and cell therapy.

Crystal structures of EfeB and EfeO in a bacterial siderophore-independent iron transport system.

EfeUOB is a siderophore-independent iron uptake mechanism in bacteria. EfeU, EfeO, and EfeB are a permease, an iron-binding or electron-transfer protein, and a peroxidase, respectively. A Gram-negative bacterium, Sphingomonas sp. strain A1, encodes EfeU, EfeO, EfeB together with alginate-binding protein Algp7, a truncated EfeO-like protein (EfeOII), in the genome. The typical EfeO (EfeOI) consists of N-terminal cupredoxin and C-terminal M75 peptidase domains. Here, we detail the structure and function of bacterial EfeB and EfeO. Crystal structures of strain A1 EfeB and Escherichia coli EfeOI were determined at 2.30 Å and 1.85 Å resolutions, respectively. A molecule of heme involved in oxidase activity was bound to the C-terminal Dyp peroxidase domain of EfeB. Two domains of EfeOI were connected by a short loop, and a zinc ion was bound to four residues, Glu156, Glu159, Asp173, and Glu255, in the C-terminal M75 peptidase domain. These residues formed tetrahedron geometry suitable for metal binding and are well conserved among various EfeO proteins including Algp7 (EfeOII), although the metal-binding site (HxxE) is proposed in the C-terminal M75 peptidase domain. This is the first report on structure of a typical EfeO with two domains, postulating a novel metal-binding motif "ExxE-//-D-//-E" in the EfeO C-terminal M75 peptidase domain.

Azurin-Derived Peptides: Comparison of Nickel- and Copper-Binding Properties.

Metalloproteins are an important class of proteins involved in metal uptake, transport, and electron-transfer reactions. Mimicking the active sites of these proteins through miniaturization is an active area of research with applications in biotechnology and medicine. Azurin is a 128-residue copper-binding cupredoxin protein involved in electron-transfer reactions. Previous studies have reported on the copper-binding-induced spectroscopic and structural properties of peptide loops (11 and 13 residues) from azurin. These azurin peptides exhibited novel stoichiometries. However, the underlying mechanism of fluorescence quenching upon copper binding remains to be understood, whether it is due to electron transfer, energy transfer, or both. Here, we report nickel-binding-associated spectroscopic and structural properties of the azurin peptides. They develop a ß-turn upon nickel binding as seen in circular dichroism and exhibit electronic transitions centered at 270 and 450 nm. Unlike copper, which exhibited 1:1 and 1:2 peptide:metal stoichiometries, nickel exhibited only a 1:1 stoichiometry. Tryptophan-containing peptides showed fluorescence quenching upon nickel binding, which is due to electron transfer. These results further suggest that the quenching in copper-bound peptides is also due to electron transfer, which could not be ascertained in previous studies. Overall, azurin peptides provide a platform for studying metal-induced structural and spectroscopic properties using transition-metal ions.

Copper-induced spectroscopic and structural changes in short peptides derived from azurin.

The active sites of metalloproteins may be mimicked by designing peptides that bind to their respective metal ions. Studying the binding of protein ligands to metal ions along with the associated structural changes is important in understanding metal uptake, transport and electron transfer functions of proteins. Copper-binding metalloprotein azurin is a 128-residue electron transfer protein with a redox-active copper cofactor. Here, we report the copper-binding associated spectroscopic and structural properties of peptide loops (11 and 13 residues) from the copper-binding site of azurin. These peptides develop a ß-turn upon copper-binding with a 1:1 Cu2+:peptide stoichiometry as seen in circular dichroism and exhibit electronic transitions centered at 340 nm and 540 nm. Further addition of copper develops a helical feature along with a shift in the absorption maxima to ~360 nm and ~580 nm at 2:1 Cu2+:peptide stoichiometry, indicating stoichiometric dependence of copper-binding geometry. Mass spectrometry indicates the copper-binding to cysteine, histidine and methionine in the peptide with 1:1 stoichiometry, and interestingly, dimerization through a disulfide linkage at 2:1 stoichiometry, as observed previously for denatured azurin. Fluorescence quenching studies on peptides with tryptophan further confirm the copper-binding induced changes in the two peptides are bi-phasic.

Tuning Structure and Dynamics of Blue Copper Azurin Junctions via Single Amino-Acid Mutations.

In the growing field of biomolecular electronics, blue-copper Azurin stands out as one of the most widely studied protein in single-molecule contacts. Interestingly, despite the paramount importance of the structure/dynamics of molecular contacts in their transport properties, these factors remain largely unexplored from the theoretical point of view in the context of single Azurin junctions. Here we address this issue using all-atom Molecular Dynamics (MD) of Pseudomonas Aeruginosa Azurin adsorbed to a Au(111) substrate. In particular, we focus on the structure and dynamics of the free/adsorbed protein and how these properties are altered upon single-point mutations. The results revealed that wild-type Azurin adsorbs on Au(111) along two well defined configurations: one tethered via cysteine groups and the other via the hydrophobic pocket surrounding the Cu 2 + . Surprisingly, our simulations revealed that single amino-acid mutations gave rise to a quenching of protein vibrations ultimately resulting in its overall stiffening. Given the role of amino-acid vibrations and reorientation in the dehydration process at the protein-water-substrate interface, we suggest that this might have an effect on the adsorption process of the mutant, giving rise to new adsorption configurations.

Solochrome cyanine: A histological stain for cobalt-chromium wear particles in metal-on-metal periprosthetic tissues.

Metal-on-metal (MoM) hip arthroplasties produce abundant implant-derived wear debris composed mainly of cobalt (Co) and chromium (Cr). Cobalt-chromium (Co-Cr) wear particles are difficult to identify histologically and need to be distinguished from other wear particle types and endogenous components (e.g., haemosiderin, fibrin) which may be present in MoM periprosthetic tissues. In this study we sought to determine whether histological stains that have an affinity for metals are useful in identifying Co-Cr wear debris in MoM periprosthetic tissues. Histological sections of periprosthetic tissue from 30 failed MoM hip arthroplasties were stained with haematoxylin-eosin (HE), Solochrome Cyanine (SC), Solochrome Azurine (SA) and Perls' Prussian Blue (PB). Sections of periprosthetic tissue from 10 cases of non-MoM arthroplasties using other implant biomaterials, including titanium, ceramic, polymethylmethacrylate (PMMA) and ultra-high molecular weight polyethylene (UHMWP) were similarly analysed. Sections of 10 cases of haemosiderin-containing knee tenosynovial giant cell tumour (TSGCT) were also stained with HE, SC, SA and PB. In MoM periprosthetic tissues, SC stained metal debris in phagocytic macrophages and in the superficial necrotic zone which exhibited little or no trichrome staining for fibrin. In non-MoM periprosthetic tissues, UHMWP, PMMA, ceramic and titanium particles were not stained by SC. Prussian Blue, but not SC or SA, stained haemosiderin deposits in MoM periprosthetic tissues and TSGT. Our findings show that SC staining (most likely Cr-associated) is useful in distinguishing Co-Cr wear particles from other metal/non-metal wear particles types in histological preparations of periprosthetic tissue and that SC reliably distinguishes haemosiderin from Co-Cr wear debris.

Raman Evidence of p53-DBD Disorder Decrease upon Interaction with the Anticancer Protein Azurin.

Raman spectroscopy, which is a suitable tool to elucidate the structural properties of intrinsically disordered proteins, was applied to investigate the changes in both the structure and the conformational heterogeneity of the DNA-binding domain (DBD) belonging to the intrinsically disordered protein p53 upon its binding to Azurin, an electron-transfer anticancer protein from Pseudomonas aeruginosa. The Raman spectra of the DBD and Azurin, isolated in solution or forming a complex, were analyzed by a combined analysis based on peak inspection, band convolution, and principal component analysis (PCA). In particular, our attention was focused on the Raman peaks of Tyrosine and Tryptophan residues, which are diagnostic markers of protein side chain environment, and on the Amide I band, of which the deconvolution allows us to extract information about α-helix, ß-sheet, and random coil contents. The results show an increase of the secondary structure content of DBD concomitantly with a decrease of its conformational heterogeneity upon its binding to Azurin. These findings suggest an Azurin-induced conformational change of DBD structure with possible implications for p53 functionality.

Synergistic effect of granzyme B-azurin fusion protein on breast cancer cells.

As one of the most prevalent malignancies, breast cancer still remains a significant risk for public health. Common therapeutic strategies include invasive surgery, chemotherapy and anti-herceptin antibodies. Adverse effects, drug resistance and low efficacy of current therapies necessitates the emergence of more effective platforms. Naturally released by the immune system, granzyme B activates multiple pro-apoptotic pathways by cleaving critical substrates. Bacterial cupredoxin, azurin, selectively targets cancer cells via a p53-dependent pathway. Fused by a linker, GrB-Azurin fusion protein was overexpressed in HEK293T cells, and purified by metal chromatography. SDS-PAGE, Western blotting and ELISA were performed to confirm successful expression, purification and analyze binding properties of the fusion protein. After treatment of various breast cancer cell lines with increasing concentrations of GrB-Azurin, quantitative real-time RT-PCR was used to measure relative expression of p21, Fas and DR5 pro-apoptotic genes. The results of DNA fragmentation and WST-1 cell viability assays indicated significant apoptosis induction in MDA-MB-231, MCF7 and SK-BR-3 cells, while insignificant cytotoxicity was detected on MCF 10A normal breast cells. Herein, we report the development of a novel biotherapeutic against breast cancer. Selective effectiveness of GrB-Azurin fusion protein on different breast cancer cells highlighted the potential of the designed construct as a candidate anti-cancer biodrug.

Clarifying the Copper Coordination Environment in a de Novo Designed Red Copper Protein.

Cupredoxins are copper-dependent electron-transfer proteins that can be categorized as blue, purple, green, and red depending on the spectroscopic properties of the Cu(II) bound forms. Interestingly, despite significantly different first coordination spheres and nuclearity, all cupredoxins share a common Greek Key ß-sheet fold. We have previously reported the design of a red copper protein within a completely distinct three-helical bundle protein, α3DChC2. (1) While this design demonstrated that a ß-barrel fold was not requisite to recapitulate the properties of a native cupredoxin center, the parent peptide α3D was not sufficiently stable to allow further study through additional mutations. Here we present the design of an elongated protein GRANDα3D (GRα3D) with Δ Gu = -11.4 kcal/mol compared to the original design's -5.1 kcal/mol. Diffraction quality crystals were grown of GRα3D (a first for an α3D peptide) and solved to a resolution of 1.34 Å. Examination of this structure suggested that Glu41 might interact with the Cu in our previously reported red copper protein. The previous bis(histidine)(cysteine) site (GRα3DChC2) was designed into this new scaffold and a series of variant constructs were made to explore this hypothesis. Mutation studies around Glu41 not only prove the proposed interaction, but also enabled tuning of the constructs' hyperfine coupling constant from 160 to 127 × 10-4 cm-1. X-ray absorption spectroscopy analysis is consistent with these hyperfine coupling differences being the result of variant 4p mixing related to coordination geometry changes. These studies not only prove that an Glu41-Cu interaction leads to the α3DChC2 construct's red copper protein like spectral properties, but also exemplify the exact control one can have in a de novo construct to tune the properties of an electron-transfer Cu site.

Azurin interaction with the lipid raft components ganglioside GM-1 and caveolin-1 increases membrane fluidity and sensitivity to anti-cancer drugs.

Membrane lipid rafts are highly ordered microdomains and essential components of plasma membranes. In this work, we demonstrate that azurin uptake by cancer cells is, in part, mediated by caveolin-1 and GM-1, lipid rafts' markers. This recognition is mediated by a surface exposed hydrophobic core displayed by azurin since the substitution of a phenylalanine residue in position 114 facing the hydrophobic cavity by alanine impacts such interactions, debilitating the uptake of azurin by cancer cells. Treating of cancer cells with azurin leads to a sequence of events: alters the lipid raft exposure at plasma membranes, causes a decrease in the plasma membrane order as examined by Laurdan two-photon imaging and leads to a decrease in the levels of caveolin-1. Caveolae, a subset of lipid rafts characterized by the presence of caveolin-1, are gaining increasing recognition as mediators in tumor progression and resistance to standard therapies. We show that azurin inhibits growth of cancer cells expressing caveolin-1, and this inhibition is only partially observed with mutant azurin. Finally, the simultaneous administration of azurin with anticancer therapeutic drugs (paclitaxel and doxorubicin) results in an enhancement in their activity, contrary to the mutated protein.

Electrochemical behavior of pyrite in sulfuric acid in presence of amino acids belonging to the amino acid sequence of rusticyanin.

The impact of different concentrations of three amino acids (cysteine, histidine and methionine) which are part of the amino acid sequence of rusticyanin on dissolution of pyrite is investigated by the application of electrochemical techniques. Cyclic voltammetric studies conducted in the anodic direction from corrosion potential have shown that in the vicinity of corrosion potential, histidine and methionine do not influence dissolution of pyrite independently on their concentrations. On the other hand, cysteine and solutions of these amino acids in the molar ratios Cys:His:Met/1:1:1 and Cys:His:Met/1:2:1 accelerate dissolution at concentrations 10-2 mol L-1 and 10-3 mol L-1. Potentiodynamic polarization measurements showed that methionine does not affect the anodic and cathodic dissolution at all concentrations, while histidine does not affect significantly on the anodic dissolution at all concentrations. Cysteine and solutions of three amino acids in the molar ratio Cys:His:Met/1:1:1 and Cys:His:Met/1:2:1 cause intensive cathodic inhibition and anodic activation at concentrations 10-2 mol L-1 and 10-3 mol L-1 respectively.

Direct Measurement of Charge Regulation in Metalloprotein Electron Transfer.

Determining whether a protein regulates its net electrostatic charge during electron transfer (ET) will deepen our mechanistic understanding of how polypeptides tune rates and free energies of ET (e.g., by affecting reorganization energy, and/or redox potential). Charge regulation during ET has never been measured for proteins because few tools exist to measure the net charge of a folded protein in solution at different oxidation states. Herein, we used a niche analytical tool (protein charge ladders analyzed with capillary electrophoresis) to determine that the net charges of myoglobin, cytochrome c, and azurin change by 0.62±0.06, 1.19±0.02, and 0.51±0.04 units upon single ET. Computational analysis predicts that these fluctuations in charge arise from changes in the pKa  values of multiple non-coordinating residues (predominantly histidine) that involve between 0.42-0.90 eV. These results suggest that ionizable residues can tune the reactivity of redox centers by regulating the net charge of the entire protein-cofactor-solvent complex.

Differences in the mechanical unfolding pathways of apo- and copper-bound azurins.

Metalloproteins carry out diverse biological functions including metal transport, electron transfer, and catalysis. At present, the influence of metal cofactors on metalloprotein stability is not well understood. Here, we report the mechanical stability and unfolding pathway of azurin, a cupredoxin family protein with ß-barrel topology and type I copper-binding centre. Single-molecule force spectroscopy (SMFS) experiments reveal 2-state and 3-state unfolding pathways for apo-azurin. The intermediate in the 3-state pathway occurs at an unfolding contour length of 7.5 nm from the native state. Steered molecular dynamics (SMD) simulations show that apo-azurin unfolds via a first transition state (TS) where ß2Β-ß8 and ß7-ß8 strand pairs rupture to form the intermediate, which subsequently unfolds by the collective rupture of remaining strands. SMFS experiments on holo-azurin exhibit an additional 4-state pathway besides the 2-state and 3-state pathways. The unfolding contour length leading to the first intermediate is 6.7 nm suggesting a sequestration of ~1 nm polypeptide chain length by the copper. SMD simulations reveal atomistic details of the copper sequestration and predict a combined ß4-ß7 pair and copper coordination sphere rupture to create the third TS in the 4-state pathway. Our systematic studies provide detailed mechanistic insights on modulation of protein mechanical properties by metal-cofactors.