Multiple wavelength anomalous diffraction (MAD) crystal structure of rusticyanin: a highly oxidizing cupredoxin with extreme acid stability

Kinetic, electrochemical and spectral characterization of bacterial and archaeal rusticyanins; unexpected stability issues and consequences for applications in biotechnology

Motivated by the ambition to establish an enzyme-driven bioleaching pathway for copper extraction, properties of the Type-1 copper protein rusticyanin from Acidithiobacillus ferrooxidans (AfR) were compared with those from an ancestral form of this enzyme (N0) and an archaeal enzyme identified in Ferroplasma acidiphilum (FaR). While both N0 and FaR show redox potentials similar to that of AfR their electron transport rates were significantly slower. The lack of a correlation between the redox potentials and electron transfer rates indicates that AfR and its associated electron transfer chain evolved to specifically facilitate the efficient conversion of the energy of iron oxidation to ATP formation. In F. acidiphilum this pathway is not as efficient unless it is up-regulated by an as of yet unknown mechanism. In addition, while the electrochemical properties of AfR were consistent with previous data, previously unreported behavior was found leading to a form that is associated with a partially unfolded form of the protein. The cyclic voltammetry (CV) response of AfR immobilized onto an electrode showed limited stability, which may be connected to the presence of the partially unfolded state of this protein. Insights gained in this study may thus inform the engineering of optimized rusticyanin variants for bioleaching processes as well as enzyme-catalyzed solubilization of copper-containing ores such as chalcopyrite.

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.

AI-driven pan-proteome analyses reveal insights into the biohydrometallurgical properties of Acidithiobacillia

Microorganism-mediated biohydrometallurgy, a sustainable approach for metal recovery from ores, relies on the metabolic activity of acidophilic bacteria. Acidithiobacillia with sulfur/iron-oxidizing capacities are extensively studied and applied in biohydrometallurgy-related processes. However, only 14 distinct proteins from Acidithiobacillia have experimentally determined structures currently available. This significantly hampers in-depth investigations of Acidithiobacillia's structure-based biological mechanisms pertaining to its relevant biohydrometallurgical processes. To address this issue, we employed a state-of-the-art artificial intelligence (AI)-driven approach, with a median model confidence of 0.80, to perform high-quality full-chain structure predictions on the pan-proteome (10,458 proteins) of the type strain Acidithiobacillia. Additionally, we conducted various case studies on de novo protein structural prediction, including sulfate transporter and iron oxidase, to demonstrate how accurate structure predictions and gene co-occurrence networks can contribute to the development of mechanistic insights and hypotheses regarding sulfur and iron utilization proteins. Furthermore, for the unannotated proteins that constitute 35.8% of the Acidithiobacillia proteome, we employed the deep-learning algorithm DeepFRI to make structure-based functional predictions. As a result, we successfully obtained gene ontology (GO) terms for 93.6% of these previously unknown proteins. This study has a significant impact on improving protein structure and function predictions, as well as developing state-of-the-art techniques for high-throughput analysis of large proteomic data.

Engineering Polyhistidine Tags on Surface Proteins of Acidithiobacillus ferrooxidans: Impact of Localization on the Binding and Recovery of Divalent Metal Cations

Metal processing using microorganisms has many advantages including the potential for reduced environmental impacts as compared to conventional technologies.Acidithiobacillus ferrooxidansis an iron- and sulfur-oxidizing chemolithoautotroph that is known to participate in metal bioleaching, and its metabolic capabilities have been exploited for industrial-scale copper and gold biomining. In addition to bioleaching, microorganisms could also be engineered for selective metal binding, enabling new opportunities for metal bioseparation and recovery. Here, we explored the ability of polyhistidine (polyHis) tags appended to two recombinantly expressed endogenous proteins to enhance the metal binding capacity of A. ferrooxidans. The genetically engineered cells achieved enhanced cobalt and copper binding capacities, and the Langmuir isotherm captures their interaction behavior with these divalent metals. Additionally, the cellular localization of the recombinant proteins correlated with kinetic modeling of the binding interactions, where the outer membrane-associated polyHis-tagged licanantase peptide bound the metals faster than the periplasmically expressed polyHis-tagged rusticyanin protein. The selectivity of the polyHis sequences for cobalt over copper from mixed metal solutions suggests potential utility in practical applications, and further engineering could be used to create metal-selective bioleaching microorganisms.

Symmetry Breaking in Solution-Phase [Cu(tsc)2(H2O)2]2+: Emergent Asymmetry in Cu-S Distances and in Covalence

The structure of aqueous Cu(II)-bis-thiosemicarbazide, [Cu(tsc)₂]²⁺, is reported following EXAFS and MXAN analyses of the copper K-edge X-ray absorption (XAS) spectrum. The rising K-edge feature at 8987.1 eV is higher energy than those of crystalline models, implying unique electronic and structural solution states. EXAFS analysis (k = 2-13 Å^-1; 2 × Cu-N = 2.02 ± 0.01 Å; 2 × Cu-S = 2.27 ± 0.01 Å; Cu-O_ax = 2.41 ± 0.04 Å) could not resolve 5- versus 6-coordinate models. However, MXAN fits converged to an asymmetric broken symmetry 6-coordinate model with cis-disposed TSC ligands (Cu-O_ax = 2.07 and 2.54 Å; Cu-N = 1.94 Å, 1.98 Å; Cu-S = 2.20 Å, 2.41 Å). Transition dipole integral evaluation of the sulfur K-edge XAS 1s → 3p valence transition feature at 2470.7 eV yielded a Cu-S covalence of 0.66 e^-, indicating Cu^1.34+. The high Cu-S covalence and short Cu-S bond in aqueous [Cu(tsc)₂(H₂O)₂]²⁺ again contradict the need for a protein rack to explain the unique structure of the blue copper active site. MXAN models of dissolved Cu(II) complex ions have invariably featured broken centrosymmetry. The potential energy ground state for dissolved Cu(II) evidently includes the extended solvation field, providing a target for improved physical theory. A revised solvation model for aqueous Cu(II), |[Cu(H₂O)₅]·14H₂O|²⁺, is presented.

Multiscale Quantum Refinement Approaches for Metalloproteins

Biomolecules with metal ion(s) (e.g., metalloproteins) play many important biological roles. However, accurate structural determination of metalloproteins, particularly those containing transition metal ion(s), is challenging due to their complicated electronic structure, complex bonding of metal ions, and high number of conformations in biomolecules. Quantum refinement, which was proposed to combine crystallographic data with computational chemistry methods by several groups, can improve the local structures of some proteins. In this study, a quantum refinement method combining several multiscale computational schemes with experimental (X-ray diffraction) information was developed for metalloproteins. Various quantum refinement approaches using different ONIOM (our own N-layered integrated molecular orbital and molecular mechanics) combinations of quantum mechanics (QM), semiempirical (SE), and molecular mechanics (MM) methods were conducted to assess the performance and reliability on the refined local structure in two metalloproteins. The structures for two (Cu- or Zn-containing) metalloproteins were refined by combining two-layer ONIOM2(QM1/QM2) and ONIOM2(QM/MM) and three-layer ONIOM3(QM1/QM2/MM) schemes with experimental data. The accuracy of the quantum-refined metal binding sites was also examined and compared in these multiscale quantum refinement calculations. ONIOM3(QM/SE/MM) schemes were found to give good results with lower computational costs and were proposed to be a good choice for the multiscale computational scheme for quantum refinement calculations of metal binding site(s) in metalloproteins with high efficiency. Additionally, a two-center ONIOM approach was employed to speed up the quantum refinement calculations for the Zn metalloprotein with two remote active sites/ligands. Moreover, a recent quantum-embedding wavefunction-in-density functional theory (WF-in-DFT) method was also adopted as the high-level method in unprecedented ONIOM2(CCSD-in-B3LYP/MM) and ONIOM3(CCSD-in-B3LYP/SE/MM) calculations, which can be regarded as novel pseudo-three- and pseudo-four-layer ONIOM methods, respectively, to refine the key Zn binding site at the coupled-cluster singles and doubles (CCSD) level. These refined results indicate that multiscale quantum refinement schemes can be used to improve the structural accuracy obtained for local metal binding site(s) in metalloproteins with high efficiency.

Computational structure prediction provides a plausible mechanism for electron transfer by the outer membrane protein Cyc2 from Acidithiobacillus ferrooxidans

Cyc2 is the key protein in the outer membrane of Acidithiobacillus ferrooxidans that mediates electron transfer between extracellular inorganic iron and the intracellular central metabolism. This cytochrome c is specific for iron and interacts with periplasmic proteins to complete a reversible electron transport chain. A structure of Cyc2 has not yet been characterized experimentally. Here we describe a structural model of Cyc2, and associated proteins, to highlight a plausible mechanism for the ferrous iron electron transfer chain. A comparative modeling protocol specific for trans membrane beta barrel (TMBB) proteins in acidophilic conditions (pH ~ 2) was applied to the primary sequence of Cyc2. The proposed structure has three main regimes: Extracellular loops exposed to low-pH conditions, a TMBB, and an N-terminal cytochrome-like region within the periplasmic space. The Cyc2 model was further refined by identifying likely iron and heme docking sites. This represents the first computational model of Cyc2 that accounts for the membrane microenvironment and the acidity in the extracellular matrix. This approach can be used to model other TMBBs which can be critical for chemolithotrophic microbial growth.

CuA-based chimeric T1 copper sites allow for independent modulation of reorganization energy and reduction potential

Attaining rational modulation of thermodynamic and kinetic redox parameters of metalloproteins is a key milestone towards the (re)design of proteins with new or improved redox functions. Here we report that implantation of ligand loops from natural T1 proteins into the scaffold of a CuA protein leads to a series of distorted T1-like sites that allow for independent modulation of reduction potentials (E°') and electron transfer reorganization energies (λ). On the one hand E°' values could be fine-tuned over 120 mV without affecting λ. On the other, λ values could be modulated by more than a factor of two while affecting E°' only by a few millivolts. These results are in sharp contrast to previous studies that used T1 cupredoxin folds, thus highlighting the importance of the protein scaffold in determining such parameters.

Effect of the Met148Leu mutation on the structure and dynamics of the rusticyanin protein from Acidithiobacillus sp. FJ2

The rusticyanin protein, a blue monomeric copper protein type-1, is one of the main components in the iron-electron transfer chain of the Acidithiobacillus ferrooxidans, and is the product of the rus gene expression. Herein, first the bacterial DNA of Acidithiobacillus sp. FJ2 was extracted. Then, the rus gene sequence and the sequence amino acid rusticyanin protein were determined. The Met148Leu mutation increased the oxidase activity of the rusticyanin protein, thereby enhancing the efficiency of the bioleaching process by bacteria Acidithiobacillus ferroxidans. Met148Leu mutation was created in the rusticyanin protein, then molecular dynamics (MD) simulations and structural analysis were performed. The MD analysis of the wild-type and mutant protein demonstrated a slight instability in the mutant protein and significant instability in the active site of the mutant protein. The usefulness of this study is the genetic manipulation of the native Acidithiobacillus sp. FJ2 bacterium, which can boost the bioleaching efficiency of the bacterium to some extent, and investigating its effects on the structure of a mutant protein using computational methods.

In situ absorbance measurements: a new means to study respiratory electron transfer in chemolithotrophic microorganisms

Absorbance measurements on intact chemolithotrophic microorganisms that respire aerobically on soluble iron are described that used a novel integrating cavity absorption meter to eliminate the effects of light scattering on the experimental results. Steady state kinetic measurements on ferric iron production by intact cells revealed that the Michaelis Menten equation described the initial rates of product formation for at least 8 different chemolithotrophic microorganisms in 6 phyla distributed equally among the archaea and the Gram negative and Gram positive eubacteria. Cell-monitored turnover measurements during aerobic respiration on soluble iron by the same 12 intact microorganisms revealed six different patterns of iron-dependent absorbance changes, suggesting that there may be at least six different sets of prosthetic groups and biomolecules that can accomplish aerobic respiration on soluble iron. Detailed kinetic studies revealed that the 3-component iron respiratory chain of Acidithiobacillus ferrooxidans functioned as an ensemble with a single macroscopic rate constant when the iron-reduced proteins were oxidized in the presence of excess molecular oxygen. The principal member of this 3-component system was a cupredoxin called rusticyanin that was present in the periplasm of At. ferrooxidans at an approximate concentration of 350 mg/mL, an observation that provides new insights into the crowded environments in the periplasms of Gram negative eubacteria that conduct electrons across their periplasm. The ability to conduct direct spectrophotometric measurements under noninvasive physiological conditions represents a new and powerful approach to examine the rates and extents of biological events in situ without disrupting the complexity of the live cellular environment.

Protein stability governed by its structural plasticity is inferred by physicochemical factors and salt bridges

Several organisms, specifically microorganisms survive in a wide range of harsh environments including extreme temperature, pH, and salt concentration. We analyzed systematically a large number of protein sequences with their structures to understand their stability and to discriminate extremophilic proteins from their non-extremophilic orthologs. Our results highlighted that the strategy for the packing of the protein core was influenced by the environmental stresses through substitutive structural events through better ionic interaction. Statistical analysis showed that a significant difference in number and composition of amino acid exist among them. The negative correlation of pairwise sequence alignments and structural alignments indicated that most of the extremophile and non-extremophile proteins didn’t contain any association for maintaining their functional stability. A significant numbers of salt bridges were noticed on the surface of the extremostable proteins. The Ramachandran plot data represented more occurrences of amino acids being present in helix and sheet regions of extremostable proteins. We also found that a significant number of small nonpolar amino acids and moderate number of charged amino acids like Arginine and Aspartic acid represented more nonplanar Omega angles in their peptide bond. Thus, extreme conditions may predispose amino acid composition including geometric variability for molecular adaptation of extremostable proteins against atmospheric variations and associated changes under natural selection pressure. The variation of amino acid composition and structural diversifications in proteins play a major role in evolutionary adaptation to mitigate climate change.

Iron and sulfur oxidation pathways of Acidithiobacillus ferrooxidans

Acidithiobacillus ferrooxidans is a gram-negative, autotrophic and rod-shaped bacterium. It can gain energy through the oxidation of Fe(II) and reduced inorganic sulfur compounds for bacterial growth when oxygen is sufficient. It can be used for bio-leaching and bio-oxidation and contributes to the geobiochemical circulation of metal elements and nutrients in acid mine drainage environments. The iron and sulfur oxidation pathways of A. ferrooxidans play key roles in bacterial growth and survival under extreme circumstances. Here, the electrons transported through the thermodynamically favourable pathway for the reduction to H₂O (downhill pathway) and against the redox potential gradient reduce to NAD(P)(H) (uphill pathway) during the oxidation of Fe(II) were reviewed, mainly including the electron transport carrier, relevant operon and regulation of its expression. Similar to the electron transfer pathway, the sulfur oxidation pathway of A. ferrooxidans, related genes and operons, sulfur oxidation mechanism and sulfur oxidase system are systematically discussed.

Comparative Genomic Analysis Reveals the Distribution, Organization, and Evolution of Metal Resistance Genes in the Genus Acidithiobacillus

Members of the genus Acidithiobacillus, which can adapt to extremely high concentrations of heavy metals, are universally found at acid mine drainage (AMD) sites. Here, we performed a comparative genomic analysis of 37 strains within the genus Acidithiobacillus to answer the untouched questions as to the mechanisms and the evolutionary history of metal resistance genes in Acidithiobacillus spp. The results showed that the evolutionary history of metal resistance genes in Acidithiobacillus spp. involved a combination of gene gains and losses, horizontal gene transfer (HGT), and gene duplication. Phylogenetic analyses revealed that metal resistance genes in Acidithiobacillus spp. were acquired by early HGT events from species that shared habitats with Acidithiobacillus spp., such as Acidihalobacter, Thiobacillus, Acidiferrobacter, and Thiomonas species. Multicopper oxidase genes involved in copper detoxification were lost in iron-oxidizing Acidithiobacillus ferridurans, Acidithiobacillus ferrivorans, and Acidithiobacillus ferrooxidans and were replaced by rusticyanin genes during evolution. In addition, widespread purifying selection and the predicted high expression levels emphasized the indispensable roles of metal resistance genes in the ability of Acidithiobacillus spp. to adapt to harsh environments. Altogether, the results suggested that Acidithiobacillus spp. recruited and consolidated additional novel functionalities during the adaption to challenging environments via HGT, gene duplication, and purifying selection. This study sheds light on the distribution, organization, functionality, and complex evolutionary history of metal resistance genes in Acidithiobacillus spp.IMPORTANCE Horizontal gene transfer (HGT), natural selection, and gene duplication are three main engines that drive the adaptive evolution of microbial genomes. Previous studies indicated that HGT was a main adaptive mechanism in acidophiles to cope with heavy-metal-rich environments. However, evidences of HGT in Acidithiobacillus species in response to challenging metal-rich environments and the mechanisms addressing how metal resistance genes originated and evolved in Acidithiobacillus are still lacking. The findings of this study revealed a fascinating phenomenon of putative cross-phylum HGT, suggesting that Acidithiobacillus spp. recruited and consolidated additional novel functionalities during the adaption to challenging environments via HGT, gene duplication, and purifying selection. Altogether, the insights gained in this study have improved our understanding of the metal resistance strategies of Acidithiobacillus spp.

Geometrical Description of Protein Structural Motifs

We present a geometrical method that can identify secondary structural motifs in proteins via angular correlations. The method uses crystal structure coordinates to calculate angular and radial signatures of each residue relative to an external reference point as the number of nearest-neighbor residues increases. We apply our approach to the blue copper protein amicyanin using the copper cofactor as the external reference point. We define a signature termed Δβ which describes the change in angular correlation as the span of nearest neighbor residues increases. We find that three turn regions of amicyanin harbor residues with Δβ near zero, while residues in other secondary structures have Δβ greater than zero: for β-strands, Δβ changes gradually residue by residue along the strand. Extension of our analysis to other blue copper proteins demonstrated that the noted structural trends are general. Importantly, our geometrical description of the folded protein accounts for all forces holding the structure together. Through this analysis, we identified some of the turns in amicyanin as symmetrical anchor points.

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.

Prediction of Reduction Potentials of Copper Proteins with Continuum Electrostatics and Density Functional Theory

Blue copper proteins, such as azurin, show dramatic changes in Cu²⁺/Cu⁺ reduction potential upon mutation over the full physiological range. Hence, they have important functions in electron transfer and oxidation chemistry and have applications in industrial biotechnology. The details of what determines these reduction potential changes upon mutation are still unclear. Moreover, it has been difficult to model and predict the reduction potential of azurin mutants and currently no unique procedure or workflow pattern exists. Furthermore, high‐level computational methods can be accurate but are too time consuming for practical use. In this work, a novel approach for calculating reduction potentials of azurin mutants is shown, based on a combination of continuum electrostatics, density functional theory and empirical hydrophobicity factors. Our method accurately reproduces experimental reduction potential changes of 30 mutants with respect to wildtype within experimental error and highlights the factors contributing to the reduction potential change. Finally, reduction potentials are predicted for a series of 124 new mutants that have not yet been investigated experimentally. Several mutants are identified that are located well over 10 Å from the copper center that change the reduction potential by more than 85 mV. The work shows that secondary coordination sphere mutations mostly lead to long‐range electrostatic changes and hence can be modeled accurately with continuum electrostatics.

Multiple Osmotic Stress Responses in Acidihalobacter prosperus Result in Tolerance to Chloride Ions

Extremely acidophilic microorganisms (pH optima for growth of ≤3) are utilized for the extraction of metals from sulfide minerals in the industrial biotechnology of “biomining.” A long term goal for biomining has been development of microbial consortia able to withstand increased chloride concentrations for use in regions where freshwater is scarce. However, when challenged by elevated salt, acidophiles experience both osmotic stress and an acidification of the cytoplasm due to a collapse of the inside positive membrane potential, leading to an influx of protons. In this study, we tested the ability of the halotolerant acidophile Acidihalobacter prosperus to grow and catalyze sulfide mineral dissolution in elevated concentrations of salt and identified chloride tolerance mechanisms in Ac. prosperus as well as the chloride susceptible species, Acidithiobacillus ferrooxidans. Ac. prosperus had optimum iron oxidation at 20 g L⁻¹ NaCl while At. ferrooxidans iron oxidation was inhibited in the presence of 6 g L⁻¹ NaCl. The tolerance to chloride in Ac. prosperus was consistent with electron microscopy, determination of cell viability, and bioleaching capability. The Ac. prosperus proteomic response to elevated chloride concentrations included the production of osmotic stress regulators that potentially induced production of the compatible solute, ectoine uptake protein, and increased iron oxidation resulting in heightened electron flow to drive proton export by the F₀F₁ ATPase. In contrast, At. ferrooxidans responded to low levels of Cl⁻ with a generalized stress response, decreased iron oxidation, and an increase in central carbon metabolism. One potential adaptation to high chloride in the Ac. prosperus Rus protein involved in ferrous iron oxidation was an increase in the negativity of the surface potential of Rus Form I (and Form II) that could help explain how it can be active under elevated chloride concentrations. These data have been used to create a model of chloride tolerance in the salt tolerant and susceptible species Ac. prosperus and At. ferrooxidans, respectively.

The Rise of Radicals in Bioinorganic Chemistry

Prior to 1950, the consensus was that biological transformations occurred in two-electron steps, thereby avoiding the generation of free radicals. Dramatic advances in spectroscopy, biochemistry, and molecular biology have led to the realization that protein-based radicals participate in a vast array of vital biological mechanisms. Redox processes involving high-potential intermediates formed in reactions with O₂ are particularly susceptible to radical formation. Clusters of tyrosine (Tyr) and tryptophan (Trp) residues have been found in many O₂-reactive enzymes, raising the possibility that they play an antioxidant protective role. In blue copper proteins with plastocyanin-like domains, Tyr/Trp clusters are uncommon in the low-potential single-domain electron-transfer proteins and in the two-domain copper nitrite reductases. The two-domain muticopper oxidases, however, exhibit clusters of Tyr and Trp residues near the trinuclear copper active site where O₂ is reduced. These clusters may play a protective role to ensure that reactive oxygen species are not liberated during O₂ reduction.

The Multicenter Aerobic Iron Respiratory Chain of Acidithiobacillus ferrooxidans Functions as an Ensemble with a Single Macroscopic Rate Constant

Electron transfer reactions among three prominent colored proteins in intact cells of Acidithiobacillus ferrooxidans were monitored using an integrating cavity absorption meter that permitted the acquisition of accurate absorbance data in suspensions of cells that scattered light. The concentrations of proteins in the periplasmic space were estimated to be 350 and 25 mg/ml for rusticyanin and cytochrome c, respectively; cytochrome a was present as one molecule for every 91 nm(2) in the cytoplasmic membrane. All three proteins were rapidly reduced to the same relative extent when suspensions of live bacteria were mixed with different concentrations of ferrous ions at pH 1.5. The subsequent molecular oxygen-dependent oxidation of the multicenter respiratory chain occurred with a single macroscopic rate constant, regardless of the proteins' in vitro redox potentials or their putative positions in the aerobic iron respiratory chain. The crowded electron transport proteins in the periplasm of the organism constituted an electron conductive medium where the network of protein interactions functioned in a concerted fashion as a single ensemble with a standard reduction potential of 650 mV. The appearance of product ferric ions was correlated with the reduction levels of the periplasmic electron transfer proteins; the limiting first-order catalytic rate constant for aerobic respiration on iron was 7,400 s(-1). The ability to conduct direct spectrophotometric studies under noninvasive physiological conditions represents a new and powerful approach to examine the extent and rates of biological events in situ without disrupting the complexity of the live cellular environment.

Bioleaching of a low-grade nickel-copper sulfide by mixture of four thermophiles

This study investigated thermophilic bioleaching of a low grade nickel-copper sulfide using mixture of four acidophilic thermophiles. Effects of 0.2g/L l-cysteine on the bioleaching process were further evaluated. It aimed at offering new alternatives for enhancing metal recoveries from nickel-copper sulfide. Results showed a recovery of 80.4% nickel and 68.2% copper in 16-day bioleaching without l-cysteine; while 83.7% nickel and 81.4% copper were recovered in the presence of l-cysteine. Moreover, nickel recovery was always higher than copper recovery. l-Cysteine was found contributing to lower pH value, faster microbial growth, higher Oxidation-Reduction Potential (ORP), higher zeta potential and absorbing on the sulfide surfaces through amino, carboxyl and sulfhydryl groups. X-ray Diffraction (XRD) patterns of leached residues showed generation of S, jarosite and ammoniojarosite. Denaturing Gradient Gel Electrophoresis (DGGE) results revealed that l-cysteine could have variant impacts on different microorganisms and changed the microbial community composition dramatically during nickel-copper sulfide bioleaching.

Highly active antibacterial ferrocenoylated or ruthenocenoylated Arg-Trp peptides can be discovered by anl-to-dsubstitution scan

By taking a systematic approach several short organometallic AMP conjugates were discovered which have very low hemolytic activity but virtually the same antimicrobial activity against MRSA as that of vancomycin.

Mineral respiration under extreme acidic conditions: from a supramolecular organization to a molecular adaptation in Acidithiobacillus ferrooxidans

Acidithiobacillus ferrooxidans is an acidophilic chemolithoautotrophic Gram-negative bacterium that can derive energy from the oxidation of ferrous iron at pH 2 using oxygen as electron acceptor. The study of this bacterium has economic and fundamental biological interest because of its use in the industrial extraction of copper and uranium from ores. For this reason, its respiratory chain has been analysed in detail in recent years. Studies have shown the presence of a functional supercomplex that spans the outer and the inner membranes and allows a direct electron transfer from the extracellular Fe2+ ions to the inner membrane cytochrome c oxidase. Iron induces the expression of two operons encoding proteins implicated in this complex as well as in the regeneration of the reducing power. Most of these are metalloproteins that have been characterized biochemically, structurally and biophysically. For some of them, the molecular basis of their adaptation to the periplasmic acidic environment has been described. Modifications in the metal surroundings have been highlighted for cytochrome c and rusticyanin, whereas, for the cytochrome c oxidase, an additional partner that maintains its stability and activity has been demonstrated recently.

Inner- and outer-sphere metal coordination in blue copper proteins

Blue copper proteins (BCPs) comprise classic cases of Nature's profound control over the electronic structures and chemical reactivity of transition metal ions. Early studies of BCPs focused on their inner coordination spheres, that is, residues that directly coordinate Cu. Equally important are the electronic and geometric perturbations to these ligands provided by the outer coordination sphere. In this tribute to Hans Freeman, we review investigations that have advanced the understanding of how inner-sphere and outer-sphere coordination affects biological Cu properties.

In silico classification of proteins from acidic and neutral cytoplasms

Protein acidostability is a common problem in biopharmaceutical and other industries. However, it remains a great challenge to engineer proteins for enhanced acidostability because our knowledge of protein acidostabilization is still very limited. In this paper, we present a comparative study of proteins from bacteria with acidic (AP) and neutral cytoplasms (NP) using an integrated statistical and machine learning approach. We construct a set of 393 non-redundant AP-NP ortholog pairs and calculate a total of 889 sequence based features for these proteins. The pairwise alignments of these ortholog pairs are used to build a residue substitution propensity matrix between APs and NPs. We use Gini importance provided by the Random Forest algorithm to rank the relative importance of these features. A scoring function using the 10 most significant features is developed and optimized using a hill climbing algorithm. The accuracy of the score function is 86.01% in predicting AP-NP ortholog pairs and is 76.65% in predicting non-ortholog AP-NP pairs, suggesting that there are significant differences between APs and NPs which can be used to predict relative acidostability of proteins. The overall trends uncovered in the study can be used as general guidelines for designing acidostable proteins. To best of our knowledge, this work represents the first systematic comparative study of the acidostable proteins and their non-acidostable orthologs.

Enhanced stereoselectivity of a Cu(II) complex chiral auxiliary in the synthesis of Fmoc-L-γ-carboxyglutamic acid

L-γ-Carboxyglutamic acid (Gla) is an uncommon amino acid that binds avidly to mineral surfaces and metal ions. Herein, we report the synthesis of N-α-Fmoc-L-γ-carboxyglutamic acid γ,γ'-tert-butyl ester (Fmoc-Gla(O(t)Bu)(2)-OH), a suitably protected analogue for Fmoc-based solid-phase peptide synthesis. The residue was synthesized using a novel chiral Cu(II) complex, whose structure-based design was inspired by the blue copper protein rusticyanin. The five-coordinate complex is formed by Shiff base formation between glycine and the novel ligand (S)-2-(N-(2-methylthio)benzylprolyl)aminobenzophenone in the presence of copper. Michael addition of di-tert-butyl methylenemalonate to the α-carbon of the glycine portion of the complex occurs in a diastereoselective fashion. The resulting (S,S)-complex diastereomer can be easily purified by chromatography. Metal complex decomposition followed by Fmoc protection affords the enantiomerically pure amino acid. With the use of this novel chiral complex, the asymmetric synthesis of Fmoc-Gla(O(t)Bu)(2)-OH was completed in nine steps from thiosalicylic acid in 14.5% overall yield.

An unconventional copper protein required for cytochrome c oxidase respiratory function under extreme acidic conditions

Very little is known about the processes used by acidophile organisms to preserve stability and function of respiratory pathways. Here, we reveal a potential strategy of these organisms for protecting and keeping functional key enzymes under extreme conditions. Using Acidithiobacillus ferrooxidans, we have identified a protein belonging to a new cupredoxin subfamily, AcoP, for "acidophile CcO partner," which is required for the cytochrome c oxidase (CcO) function. We show that it is a multifunctional copper protein with at least two roles as follows: (i) as a chaperone-like protein involved in the protection of the Cu(A) center of the CcO complex and (ii) as a linker between the periplasmic cytochrome c and the inner membrane cytochrome c oxidase. It could represent an interesting model for investigating the multifunctionality of proteins known to be crucial in pathways of energy metabolism.

The 1.4 A resolution structure of Paracoccus pantotrophus pseudoazurin

Pseudoazurins are small type 1 copper proteins that are involved in the flow of electrons between various electron donors and acceptors in the bacterial periplasm, mostly under denitrifying conditions. The previously determined structure of Paracoccus pantotrophus pseudoazurin in the oxidized form was improved to a nominal resolution of 1.4 A, with R and R(free) values of 0.188 and 0.206, respectively. This high-resolution structure makes it possible to analyze the interactions between the monomers and the solvent structure in detail. Analysis of the high-resolution structure revealed the structural regions that are responsible for monomer-monomer recognition during dimer formation and for protein-protein interaction and that are important for partner recognition. The pseudoazurin structure was compared with other structures of various type 1 copper proteins and these were grouped into families according to similarities in their secondary structure; this may be useful in the annotation of copper proteins in newly sequenced genomes and in the identification of novel copper proteins.

Rationally tuning the reduction potential of a single cupredoxin beyond the natural range

Redox processes are at the heart of numerous functions in chemistry and biology, from long-range electron transfer in photosynthesis and respiration to catalysis in industrial and fuel cell research. These functions are accomplished in nature by only a limited number of redox-active agents. A long-standing issue in these fields is how redox potentials are fine-tuned over a broad range with little change to the redox-active site or electron-transfer properties. Resolving this issue will not only advance our fundamental understanding of the roles of long-range, non-covalent interactions in redox processes, but also allow for design of redox-active proteins having tailor-made redox potentials for applications such as artificial photosynthetic centres or fuel cell catalysts for energy conversion. Here we show that two important secondary coordination sphere interactions, hydrophobicity and hydrogen-bonding, are capable of tuning the reduction potential of the cupredoxin azurin over a 700 mV range, surpassing the highest and lowest reduction potentials reported for any mononuclear cupredoxin, without perturbing the metal binding site beyond what is typical for the cupredoxin family of proteins. We also demonstrate that the effects of individual structural features are additive and that redox potential tuning of azurin is now predictable across the full range of cupredoxin potentials.

Low temperature 65Cu NMR spectroscopy of the Cu+ site in azurin

(65)Cu central-transition NMR spectroscopy of the blue copper protein azurin in the reduced Cu(I) state, conducted at 18.8 T and 10 K, gave a strongly second order quadrupole perturbed spectrum, which yielded a (65)Cu quadrupole coupling constant of +/-71.2 +/- 1 MHz, corresponding to an electric field gradient of +/-1.49 atomic units at the copper site, and an asymmetry parameter of approximately 0.2. Quantum chemical calculations employing second order Møller-Plesset perturbation theory and large basis sets successfully reproduced these experimental results. Sensitivity and relaxation times were quite favorable, suggesting that NMR may be a useful probe of the electronic state of copper sites in proteins.

NMR hyperfine shifts in blue copper proteins: a quantum chemical investigation

We present the results of the first quantum chemical investigations of 1H NMR hyperfine shifts in the blue copper proteins (BCPs): amicyanin, azurin, pseudoazurin, plastocyanin, stellacyanin, and rusticyanin. We find that very large structural models that incorporate extensive hydrogen bond networks, as well as geometry optimization, are required to reproduce the experimental NMR hyperfine shift results, the best theory vs experiment predictions having R2 = 0.94, a slope = 1.01, and a SD = 40.5 ppm (or approximately 4.7% of the overall approximately 860 ppm shift range). We also find interesting correlations between the hyperfine shifts and the bond and ring critical point properties computed using atoms-in-molecules theory, in addition to finding that hyperfine shifts can be well-predicted by using an empirical model, based on the geometry-optimized structures, which in the future should be of use in structure refinement.

The role of ligand-containing loops at copper sites in proteins

Many approaches are being used to engineer metalloproteins, with most of these informed by, and aiming to further elucidate, the basic structural requirements for biological metal centers. Cupredoxins are type 1 (T1) copper-containing electron transfer (ET) proteins with a -barrel fold that is thought to constrain metal site structure. The T1 copper ion is bound by ligands mainly originating from a single active site loop whose length and structure varies. This Highlight article will focus on protein engineering studies which have investigated the role of the metal-binding loop for active site integrity and functionality. Scaffold differences are present within the cupredoxin family and their influence has also been assessed. Given the widespread occurrence of -barrel domains in nature, and the array of metal sites in proteins composed of loop regions, the studies described on this model system have implications for a variety of metalloproteins.

Modeling Study of Rusticyanin-Cytochrome C4Complex: An Insight to Possible H-Bond Mediated Recognition and Electron—Transfer Process

Rusticyanin (RCy) mediated transfer of electron to Cytochrome C(4) (Cytc(4)) from the extracellular Fe(+2) ion is primarily involved in the Thiobacillus ferrooxidans induced bio-leaching of pyrite ore and also in the metabolism of this acidophilic bacteria. The modeling studies have revealed the two possible mode of RCy-Cytc(4) complexation involving nearly the same stabilization energy approximately -15 x 10(3) kJ/mol, one through N-terminal Asp 15 and another -C terminal Glu 121 of Cytc(4) with the Cu-bonded His 143 of RCy. The Asp 15:His 143 associated complex (DH) of Cytc(4)-RCy was stabilized by the intermolecular H-bonds of the carboxyl oxygen atoms O(delta1) and O(delta2) of Asp 15 with the Nepsilon-atom of His 143 and O(b) atoms of Ala 8 and Asp 5 (of Cytc(4)) with the Thr 146 and Phe 51 (of RCy). But the other Glu 121:His 143 associated complex (EH) of Cytc(4)-RCy was stabilized by the H-bonding interaction of the oxygen atoms O(epsilon1) and O(epsilon2) of Glu 121 with the Nepsilon and Ogamma atoms of His 143 and Thr 146 of RCy. The six water molecules were present in the binding region of the two proteins in the energy minimized autosolvated DH and EH-complexes. The MD studies also revealed the presence of six interacting water molecules at the binding region between the two proteins in both the complexes. Several residues Gly 82 and 84, His 143 (RCy) were participated through the water mediated (W 389, W 430, W 413, W 431, W 373, and W 478) interaction with the Asp 15, Ile 82, and 62, Tyr 63 (Cytc(4)) in DH complex, whereas in EH complex the Phe 51, Asn 80, Tyr 146 (RCy) residues were observed to interact with Asn 108, Met 120, Glu 121 (of Cytc(4)) through the water molecules W 507, W 445, W 401, W 446, and W 440. The direct water mediated (W 478) interaction of His 143 (RCy) to Asp 15 (of Cytc(4)) was observed only in the DH complex but not in EH. These direct and water mediated H-bonding between the two respective proteins and the binding free energy with higher interacting buried surface area of the DH complex compare to other EH complex have indicated an alternative possibility of the electron transfer route through the interaction of His 143 of RCy and the N-terminal Asp 15 of Cytc(4).

Correlation of rhombic distortion of the type 1 copper site of M98Q amicyanin with increased electron transfer reorganization energy

Mutation of the axial Met ligand of the type 1 copper site of amicyanin to Ala or Gln yielded M98A amicyanin, which exhibits typical axial type 1 ligation geometry but with a water molecule providing the axial ligand, and M98Q amicyanin, which exhibits significant rhombic distortion of the type 1 site (Carrell, C. J., Ma, J. K., Antholine, W. E., Hosler, J. P., Mathews, F. S., and Davidson, V. L. (2007) Biochemistry 46, 1900-1912). Despite the change of the axial ligand, the M98Q and M98A mutations had little effect on the redox potential of copper. The true electron transfer (ET) reactions from O-quinol methylamine dehydrogenase to oxidized native and mutant amicyanins revealed that the M98A mutation had little effect on kET, but the M98Q mutation reduced kET 45-fold. Thermodynamic analysis of the latter showed that the decrease in kET was due to an increase of 0.4 eV in the reorganization energy (lambda) associated with the ET reaction to M98Q amicyanin. No change in the experimentally determined electronic coupling or ET distance was observed, confirming that the mutation had not altered the rate-determining step for ET and that this was still a true ET reaction. The basis for the increased lambda is not the nature of the atom that provides the axial ligand because each uses an oxygen from Gln in M98Q amicyanin and from water in M98A amicyanin. Comparisons of the distance of the axial copper ligand from the equatorial plane that is formed by the other three copper ligands in isomorphous crystals of native and mutant amicyanins at atomic resolution indicate an increase in distance from 0.20 A in the native to 0.42 A in M98Q amicyanin and a slight decrease in distance for M98A amicyanin. This correlates with the rhombic distortion caused by the M98Q mutation that is clearly evident in the EPR and visible absorption spectra of the protein and suggests that the extent of rhombicity of the type 1 copper site influences the magnitude of lambda.

Computational studies of Cu(II)/Met and Cu(I)/Met binding motifs relevant for the chemistry of Alzheimer's disease

A systematic study of the binding motifs of Cu(II) and Cu(I) to a methionine model peptide, namely, N-formylmethioninamide 1, has been carried out by quantum chemical computations. Geometries of the coordination modes obtained at the B3LYP/6-31G(d) level of theory are discussed in the context of copper coordination by the peptide backbone and the S atom of a methionine residue in peptides with special emphasis on Met35 of the amyloid-beta peptide (Abeta) of Alzheimer's disease. The relative binding free energies in the gas phase, DeltaG(g), are calculated at the B3LYP/6-311+G(2df,2p)//B3LYP/6-31G(d) level of theory, and the solvation affects are included by means of the COSMO model to obtain the relative binding energies in solution, DeltaG(aq). A free energy of binding, DeltaG(aq) = -19.4 kJ mol(-1), relative to aqueous Cu(II) and the free peptide is found for the most stable Cu(II)/Met complex, 12. The most stable Cu(I)/Met complex, 23, is bound by -15.6 kJ mol(-1) relative to the separated species. The reduction potential relative to the standard hydrogen electrode is estimated to be E degrees (12/23) = 0.41 V. On the basis of these results, the participation of Met35 as a low affinity binding site of Cu(II) in Abeta, and its role in the redox chemistry underlying Alzheimer's disease is discussed.

Active site comparison of CoII blue and green nitrite reductases

Copper-containing nitrite reductases (NiRs) possess type 1 (T1) and type 2 (T2) copper sites and can be either green or blue in color owing to differences at their T1 centers. The active sites of a green and a blue NiR were studied by utilizing their T1CuI/T2CoII and T1CoII/T2CoII-substituted forms. The UV/Vis spectra of these derivatives highlight the similarity of the T2 centers in these enzymes and that T1 site differences are also present in the CoII forms. The paramagnetic NMR spectra of T1CuI/T2CoII enzymes allow hyperfine shifted resonances from the three T2 His ligands to be assigned: these exhibit remarkably similar positions in the spectra of both NiRs, emphasizing the homology of the T2 centers. The addition of nitrite results in subtle alterations in the paramagnetic NMR spectra of the T1CuI/T2CoII forms at pH<7, which indicate a geometry change upon the binding of substrate. Shifted resonances from all of the T1 site ligands have been assigned and the CoII--N(His) interactions are alike, whereas the CbetaH proton resonances of the Cys ligand exhibit subtle chemical shift differences in the blue and green NiRs. The strength of the axial CoII--S(Met) interaction is similar in the two NiRs studied, but the altered conformation of the side chain of this ligand results in a dramatically different chemical shift pattern for the CgammaH protons. This indicates an alteration in the bonding of the axial ligand in these derivatives, which could be influential in the CuII proteins.

Engineering copper sites in proteins: loops confer native structures and properties to chimeric cupredoxins

The ligand-containing loops of two copper-binding electron-transfer proteins (cupredoxins) have been swapped. In the azurin (AZ) variant in which the plastocyanin (PC) sequence is introduced (AZPC), the loop adopts a conformation identical to that in PC. The reduction potential of AZPC is raised as compared to AZ and matches that of PC. In the previously published AZAMI variant (AMI = amicyanin), the shorter introduced loop adopts the same conformation as in AMI, and the reduction potential is lowered to equal that of AMI (Yanagisawa, S.; Dennison, C. J. Am. Chem. Soc. 2004, 126, 15711-15719. Li, C.; et al. Proc. Natl. Acad. Sci. U.S.A. 2006, 103, 7258-7263). Thus, the loop structure plays an important role in tuning the reduction potential of a type 1 copper site with contributions from protein dipoles in this region probably the most important feature. The structure of the loop also seems to be a major factor in controlling dissociation and protonation of the C-terminal His ligand, which can act as a switch to regulate electron-transfer reactivity. The PCAZ variant (PC with the AZ loop) possesses an active site, which is different from those of both PC and AZ, and it is assumed that the introduced loop does not adopt a structure as in AZ. This contributes to the observed instability of PCAZ and highlights that loop-scaffold interactions are important for stabilizing the active site of a cupredoxin.

The sulfhydryl group of Cys138 of rusticyanin from Acidithiobacillus ferrooxidans is crucial for copper binding

Rusticyanin is a small blue copper protein isolated from Acidithiobacillus ferrooxidans with extreme acid stability and redox potential. The protein is thought to be a principal component in the iron respiratory electron transport chain in this microorganism, but its exact role in electron transfer remains controversial. The gene of rusticyanin was cloned then overexpressed in Escherichia coli, the soluble protein was purified by one-step affinity chromatography to apparent homogeneity. It was reported that Cys138, His85 and His143 were important residues for copper binding, but the significance of Cys138 was not verified so far. We constructed the mutant expression plasmids of these three residues using site-directed mutagenesis. Mutant proteins were expressed in E. coli and purified with a nickel metal affinity column. The EPR and atomic absorption spectroscopy results confirmed that Cys138 was crucial for copper binding. Removal of the sulfhydryl group of Cys138 resulted in copper loss. Mutations of His85 and His143 showed little effect on copper binding.

Structural and mechanistic exploration of acid resistance: kinetic stability facilitates evolution of extremophilic behavior

Kinetically stable proteins are unique in that their stability is determined solely by kinetic barriers rather than by thermodynamic equilibria. To better understand how kinetic stability promotes protein survival under extreme environmental conditions, we analyzed the unfolding behavior and determined the structure of Nocardiopsis alba Protease A (NAPase), an acid-resistant, kinetically stable protease, and compared these results with a neutrophilic homolog, alpha-lytic protease (alphaLP). Although NAPase and alphaLP have the same number of acid-titratable residues, kinetic studies revealed that the height of the unfolding free energy barrier for NAPase is less sensitive to acid than that of alphaLP, thereby accounting for NAPase's improved tolerance of low pH. A comparison of the alphaLP and NAPase structures identified multiple salt-bridges in the domain interface of alphaLP that were relocated to outer regions of NAPase, suggesting a novel mechanism of acid stability in which acid-sensitive electrostatic interactions are rearranged to similarly affect the energetics of both the native state and the unfolding transition state. An acid-stable variant of alphaLP in which a single interdomain salt-bridge is replaced with a corresponding intradomain NAPase salt-bridge shows a dramatic >15-fold increase in acid resistance, providing further evidence for this hypothesis. These observations also led to a general model of the unfolding transition state structure for alphaLP protease family members in which the two domains separate from each other while remaining relatively intact themselves. These results illustrate the remarkable utility of kinetic stability as an evolutionary tool for developing longevity over a broad range of harsh conditions.

Structural analysis of the HiPIP from the acidophilic bacteria: Acidithiobacillus ferrooxidans

Hip is a high-potential iron-sulfur protein (HiPIP) isolated from the acidophilic bacterium, Acidithiobacillus ferrooxidans. In the present work, a structural model of Hip suggests that the role of proline residues is essential to stabilize the protein folding at very low pH. The presence of an unusual disulfide bridge in Hip is demonstrated using mass spectrometry and nuclear magnetic resonance. This disulfide bridge is necessary to anchor the N-terminal extremity of the protein, but is not involved in the acid stability of Hip. The structural parameters correlated with the pH dependence of Hip redox potential are also analysed on the basis of this model. Given that the same structural features can enhance acidic stability and lead to elevated redox potentials, modulation of the redox potentials of electron carriers may be necessary to achieve electron transfer at very low pH.

Atomic Resolution Crystal Structures, EXAFS, and Quantum Chemical Studies of Rusticyanin and Its Two Mutants Provide Insight into Its Unusual Properties,

Rusticyanin from the extremophile Thiobacillus ferrooxidans is a blue copper protein with unusually high redox potential and acid stability. We present the crystal structures of native rusticyanin and of its Cu site mutant His143Met at 1.27 and 1.10 A, respectively. The very high resolution of these structures allows a direct comparison with EXAFS data and with quantum chemical models of the oxidized and reduced forms of the proteins, based upon both isolated and embedded clusters and density functional theory (DFT) methods. We further predict the structure of the Cu(II) form of the His143Met mutant which has been experimentally inaccessible due to its very high redox potential. We also present metrical EXAFS data and quantum chemical calculations for the oxidized and reduced states of the Met148Gln mutant, this protein having the lowest redox potential of all currently characterized mutants of rusticyanin. These data offer new insights into the structural factors which affect the redox potential in this important class of proteins. Calculations successfully predict the structure and the order of redox potentials for the three proteins. The calculated redox potential of H143M ( approximately 400 mV greater than native rusticyanin) is consistent with the failure of readily available chemical oxidants to restore a Cu(II) species of this mutant. The structural and energetic effects of mutating the equatorial cysteine to serine, yet to be studied experimentally, are predicted to be considerable by our calculations.

An NMR view of the unfolding process of rusticyanin: Structural elements that maintain the architecture of a beta-barrel metalloprotein

The unfolding process of the blue copper protein rusticyanin (Rc) as well as its dynamic and D(2)O/H(2)O exchange properties in an incipient unfolded state have been studied by heteronuclear NMR spectroscopy. Titrations of apo, Cu(I), and Cu(II)Rc with guanidinium chloride (GdmCl) show that the copper ion stabilizes the folded species and remains bound in the completely unfolded state. The oxidized state of the copper ion is more efficient than the reduced form in this respect. The long loop of Rc (where the first ligand of the copper ion is located) is one of the most mobile domains of the protein. This region has no defined secondary structure elements and is prone to exchange its amide protons. In contrast, the last loop (including a short alpha-helix) and the last beta-strand (where the other three ligands of the metal ion are located) form the most rigid domain of the protein. The results taken as a whole suggest that the first ligand detaches from the metal ion when the protein unfolds, while the other three ligands remain bound to it. The implications of these findings for the biological folding process of Rc are also discussed.