The extended environment of mononuclear metal centers in protein structures

The energetics and evolution of oxidoreductases in deep time

The core metabolic reactions of life drive electrons through a class of redox protein enzymes, the oxidoreductases. The energetics of electron flow is determined by the redox potentials of organic and inorganic cofactors as tuned by the protein environment. Understanding how protein structure affects oxidation-reduction energetics is crucial for studying metabolism, creating bioelectronic systems, and tracing the history of biological energy utilization on Earth. We constructed ProtReDox (https://protein-redox-potential.web.app), a manually curated database of experimentally determined redox potentials. With over 500 measurements, we can begin to identify how proteins modulate oxidation-reduction energetics across the tree of life. By mapping redox potentials onto networks of oxidoreductase fold evolution, we can infer the evolution of electron transfer energetics over deep time. ProtReDox is designed to include user-contributed submissions with the intention of making it a valuable resource for researchers in this field.

Important Structural Features of Thiolate-Rich Four-Helix Bundles for Cu(I) Uptake and Removal

A family of bacterial copper storage proteins (the Csps) possess thiolate-lined four-helix bundles whose cores can be filled with Cu(I) ions. The majority of Csps are cytosolic (Csp3s), and in vitro studies carried out to date indicate that the Csp3s from Methylosinus trichosporium OB3b (MtCsp3), Bacillus subtilis (BsCsp3), and Streptomyces lividans (SlCsp3) are alike. Bioinformatics have highlighted homologues with potentially different Cu(I)-binding properties from these characterized "classical" Csp3s. Determination herein of the crystal structure of the protein (RkCsp3) from the methanotroph Methylocystis sp. strain Rockwell with Cu(I) bound identifies this as the first studied example of a new subgroup of Csp3s. The most significant structural difference from classical Csp3s is the presence of only two Cu(I) sites at the mouth of the bundle via which Cu(I) ions enter and leave. This is due to the absence of three Cys residues and a His-containing motif, which allow classical Csp3s to bind five to six Cu(I) ions in this region. Regardless, RkCsp3 exhibits rapid Cu(I) binding and the fastest measured Cu(I) removal rate for a Csp3 when using high-affinity ligands as surrogate partners. New experiments on classical Csp3s demonstrate that their His-containing motif is not essential for fast Cu(I) uptake and removal. Other structural features that could be important for these functionally relevant in vitro properties are discussed.

Classifying metal-binding sites with neural networks

To advance our ability to predict impacts of the protein scaffold on catalysis, robust classification schemes to define features of proteins that will influence reactivity are needed. One of these features is a protein's metal-binding ability, as metals are critical to catalytic conversion by metalloenzymes. As a step toward realizing this goal, we used convolutional neural networks (CNNs) to enable the classification of a metal cofactor binding pocket within a protein scaffold. CNNs enable images to be classified based on multiple levels of detail in the image, from edges and corners to entire objects, and can provide rapid classification. First, six CNN models were fine-tuned to classify the 20 standard amino acids to choose a performant model for amino acid classification. This model was then trained in two parallel efforts: to classify a 2D image of the environment within a given radius of the central metal binding site, either an Fe ion or a [2Fe-2S] cofactor, with the metal visible (effort 1) or the metal hidden (effort 2). We further used two sub-classifications of the [2Fe-2S] cofactor: (1) a standard [2Fe-2S] cofactor and (2) a Rieske [2Fe-2S] cofactor. The accuracy for the model correctly identifying all three defined features was >95%, despite our perception of the increased challenge of the metalloenzyme identification. This demonstrates that machine learning methodology to classify and distinguish similar metal-binding sites, even in the absence of a visible cofactor, is indeed possible and offers an additional tool for metal-binding site identification in proteins.

Catalysis and Electron Transfer in De Novo Designed Metalloproteins

One of the hallmark advances in our understanding of metalloprotein function is showcased in our ability to design new, non-native, catalytically active protein scaffolds. This review highlights progress and milestone achievements in the field of de novo metalloprotein design focused on reports from the past decade with special emphasis on de novo designs couched within common subfields of bioinorganic study: heme binding proteins, monometal- and dimetal-containing catalytic sites, and metal-containing electron transfer sites. Within each subfield, we highlight several of what we have identified as significant and important contributions to either our understanding of that subfield or de novo metalloprotein design as a discipline. These reports are placed in context both historically and scientifically. General suggestions for future directions that we feel will be important to advance our understanding or accelerate discovery are discussed.

Second Sphere Effects on Oxygen Reduction and Peroxide Activation by Mononuclear Iron Porphyrins and Related Systems

Activation and reduction of O₂ and H₂O₂ by synthetic and biosynthetic iron porphyrin models have proved to be a versatile platform for evaluating second-sphere effects deemed important in naturally occurring heme active sites. Advances in synthetic techniques have made it possible to install different functional groups around the porphyrin ligand, recreating artificial analogues of the proximal and distal sites encountered in the heme proteins. Using judicious choices of these substituents, several of the elegant second-sphere effects that are proposed to be important in the reactivity of key heme proteins have been evaluated under controlled environments, adding fundamental insight into the roles played by these weak interactions in nature. This review presents a detailed description of these efforts and how these have not only demystified these second-sphere effects but also how the knowledge obtained resulted in functional mimics of these heme enzymes.

Proton Relay in Iron Porphyrins for Hydrogen Evolution Reaction

The efficiency of the hydrogen evolution reaction (HER) can be facilitated by the presence of proton-transfer groups in the vicinity of the catalyst. A systematic investigation of the nature of the proton-transfer groups present and their interplay with bulk proton sources is warranted. The HERs electrocatalyzed by a series of iron porphyrins that vary in the nature and number of pendant amine groups are investigated using proton sources whose pK_a values vary from ∼9 to 15 in acetonitrile. Electrochemical data indicate that a simple iron porphyrin (FeTPP) can catalyze the HER at this Fe^I state where the rate-determining step is the intermolecular protonation of a Fe^III-H^- species produced upon protonation of the iron(I) porphyrin and does not need to be reduced to its formal Fe⁰ state. A linear free-energy correlation of the observed rate with pK_a of the acid source used suggests that the rate of the HER becomes almost independent of pK_a of the external acid used in the presence of the protonated distal residues. Protonation to the Fe^III-H^- species during the HER changes from intermolecular in FeTPP to intramolecular in FeTPP derivatives with pendant basic groups. However, the inclusion of too many pendant groups leads to a decrease in HER activity because the higher proton binding affinity of these residues slows proton transfer for the HER. These results enrich the existing understanding of how second-sphere proton-transfer residues alter both the kinetics and thermodynamics of transition-metal-catalyzed HER.

How the Local Environment of Functional Sites Regulates Protein Function

Proteins form complex biological machineries whose functions in the cell are highly regulated at both the cellular and molecular levels. Cellular regulation of protein functions involves differential gene expressions, post-translation modifications, and signaling cascades. Molecular regulation, on the other hand, involves tuning an optimal local protein environment for the functional site. Precisely how a protein achieves such an optimal environment around a given functional site is not well understood. Herein, by surveying the literature, we first summarize the various reported strategies used by certain proteins to ensure their correct functioning. We then formulate three key physicochemical factors for regulating a protein's functional site, namely, (i) its immediate interactions, (ii) its solvent accessibility, and (iii) its conformational flexibility. We illustrate how these factors are applied to regulate the functions of free/metal-bound Cys and Zn sites in proteins.

Nature-inspired electrocatalysts and devices for energy conversion

A NICE approach for the design of nature-inspired electrocatalysts and electrochemical devices for energy conversion.

PCuAC domains from methane-oxidizing bacteria use a histidine brace to bind copper

Copper is critically important for methanotrophic bacteria because their primary metabolic enzyme, particulate methane monooxygenase (pMMO), is copper-dependent. In addition to pMMO, many other copper proteins are encoded in the genomes of methanotrophs, including proteins that contain periplasmic copper-Achaperone (PCu_AC) domains. Using bioinformatics analyses, we identified three distinct classes of PCu_AC domain-containing proteins in methanotrophs, termed PmoF1, PmoF2, and PmoF3. PCu_AC domains from other types of bacteria bind a single Cu(I) ion via an HX_nMX_21/22HXM motif, which is also present in PmoF3, but PmoF1 and PmoF2 lack this motif entirely. Instead, the PCu_AC domains of PmoF1 and PmoF2 bind only Cu(II), and PmoF1 binds additional Cu(II) ions in a His-rich extension to its PCu_AC domain. Crystal structures of the PmoF1 and PmoF2 PCu_AC domains reveal that Cu(II) is coordinated by an N-terminal histidine brace HX₁₀H motif. This binding site is distinct from those of previously characterized PCu_AC domains but resembles copper centers in CopC proteins and lytic polysaccharide monooxygenase (LPMO) enzymes. Bioinformatics analysis of the entire PCu_AC family reveals previously unappreciated diversity, including sequences that contain both the HX_nMX_21/22HXM and HX₁₀H motifs, and sequences that lack either set of copper-binding ligands. These findings provide the first characterization of an additional class of copper proteins from methanotrophs, further expand the PCu_AC family, and afford new insight into the biological significance of histidine brace-mediated copper coordination.

Noncoded Amino Acids in de Novo Metalloprotein Design: Controlling Coordination Number and Catalysis

The relationship between structure and function has long been one of the major points of investigation in Biophysics. Understanding how much, or how little, of a protein's often complicated structure is necessary for its function can lead to directed therapeutic strategies and would allow one to design proteins for specific desired functions. Studying protein function by de novo design builds the functionality from the ground up in a completely unrelated and noncoded protein scaffold. Our lab has used this strategy to study heavy and transition metal binding within the TRI family of three stranded coiled coil (3SCC) constructs to understand coordination geometry and metalloenzyme catalytic control within a protein environment. These peptides contain hydrophobic layers within the interior of the 3SCC, which one can mutate to metal binding residues to create a minimal metal binding site, while solid phase synthesis allows our lab to easily incorporate a number of noncoded amino acids including d enantiomers of binding or secondary coordination sphere amino acids, penicillamine, or methylated versions of histidine. Our studies of Cd(II) binding to Cys₃ environments have determined, largely through the use of ¹¹³Cd NMR and ^111mCd PAC, that the coordination environment around a heavy metal can be controlled by incorporating noncoded amino acids in either the primary or secondary coordination spheres. We found mutating the metal binding amino acids to l-Pen can enforce trigonal Cd(II)S₃ geometry exclusively compared to the mixed coordination determined for l-Cys coordination. The same result can be achieved with secondary sphere mutations as well by incorporating d-Leu above a Cys₃. We hypothesize this latter effect is due to the increased steric packing above the metal binding site that occurs when the l-Leu oriented toward the N-terminus of the scaffold is mutated to d-Leu and oriented toward the C-terminus. Mutating the layer below Cys₃ to d-Leu instead formed a mixed 4- and 5-coordinate Cd(II)S₃(H₂O) and Cd(II)S₃(H₂O)₂ construct as steric bulk was decreased below the metal binding site. We have also applied noncoded amino acids to metalloenzyme systems by incorporating His residues that are methylated at the δ- or ε-nitrogen to enforce Cu(I) ligation to the opposite open nitrogen of His and found a 2 orders of magnitude increased catalytic efficiency for nitrite reductase activity with ε-nitrogen coordination compared to δ-nitrogen. These results exemplify the ability to tune coordination environment and catalytic efficiency within a de novo scaffold as well as the utility of noncoded amino acids to increase the chemist's toolbox. By furthering our understanding of metalloprotein design one could envision, through our use of amino acids not normally available to nature, that protein design laboratories will soon be capable of outperforming the native systems previously used as their benchmark of successful design. The ability to design proteins at this level would have far reaching and exciting benefits within various fields including medical and industrial applications.

Methylated Histidines Alter Tautomeric Preferences that Influence the Rates of Cu Nitrite Reductase Catalysis in Designed Peptides

Copper proteins have the capacity to serve as both redox active catalysts and purely electron transfer centers. A longstanding question in this field is how the function of histidine ligated Cu centers are modulated by δ vs ε-nitrogen ligation of the imidazole. Evaluating the impact of these coordination modes on structure and function by comparative analysis of deposited crystal structures is confounded by factors such as differing protein folds and disparate secondary coordination spheres that make direct comparison of these isomers difficult. Here, we present a series of de novo designed proteins using the noncanonical amino acids 1-methyl-histidine and 3-methyl-histidine to create Cu nitrite reductases where δ- or ε-nitrogen ligation is enforced by the opposite nitrogen's methylation as a means of directly comparing these two ligation states in the same protein fold. We find that ε-nitrogen ligation allows for a better nitrite reduction catalyst, displaying 2 orders of magnitude higher activity than the δ-nitrogen ligated construct. Methylation of the δ nitrogen, combined with a secondary sphere mutation we have previously published, has produced a new record for efficiency within a homogeneous aqueous system, improving by 1 order of magnitude the previously published most efficient construct. Furthermore, we have measured Michaelis-Menten kinetics on these highly active constructs, revealing that the remaining barriers to matching the catalytic efficiency ( kcat/ KM) of native Cu nitrite reductase involve both substrate binding ( KM) and catalysis ( kcat).

Benchmarking a computational design method for the incorporation of metal ion-binding sites at symmetric protein interfaces

The design of novel metal-ion binding sites along symmetric axes in protein oligomers could provide new avenues for metalloenzyme design, construction of protein-based nanomaterials and novel ion transport systems. Here, we describe a computational design method, symmetric protein recursive ion-cofactor sampling (SyPRIS), for locating constellations of backbone positions within oligomeric protein structures that are capable of supporting desired symmetrically coordinated metal ion(s) chelated by sidechains (chelant model). Using SyPRIS on a curated benchmark set of protein structures with symmetric metal binding sites, we found high recovery of native metal coordinating rotamers: in 65 of the 67 (97.0%) cases, native rotamers featured in the best scoring model while in the remaining cases native rotamers were found within the top three scoring models. In a second test, chelant models were crossmatched against protein structures with identical cyclic symmetry. In addition to recovering all native placements, 10.4% (8939/86013) of the non-native placements, had acceptable geometric compatibility scores. Discrimination between native and non-native metal site placements was further enhanced upon constrained energy minimization using the Rosetta energy function. Upon sequence design of the surrounding first-shell residues, we found further stabilization of native placements and a small but significant (1.7%) number of non-native placement-based sites with favorable Rosetta energies, indicating their designability in existing protein interfaces. The generality of the SyPRIS approach allows design of novel symmetric metal sites including with non-natural amino acid sidechains, and should enable the predictive incorporation of a variety of metal-containing cofactors at symmetric protein interfaces.

Simplest Monodentate Imidazole Stabilization of the oxy-Tyrosinase Cu2 O2 Core: Phenolate Hydroxylation through a Cu(III) Intermediate

Tyrosinases are ubiquitous binuclear copper enzymes that oxygenate to Cu(II) 2 O2 cores bonded by three histidine Nτ-imidazoles per Cu center. Synthetic monodentate imidazole-bonded Cu(II) 2 O2 species self-assemble in a near quantitative manner at -125 °C, but Nπ-ligation has been required. Herein, we disclose the syntheses and reactivity of three Nτ-imidazole bonded Cu(II) 2 O2 species at solution temperatures of -145 °C, which was achieved using a eutectic mixture of THF and 2-MeTHF. The addition of anionic phenolates affords a Cu(III) 2 O2 species, where the bonded phenolates hydroxylate to catecholates in high yields. Similar Cu(III) 2 O2 intermediates are not observed using Nπ-bonded Cu(II) 2 O2 species, hinting that Nτ-imidazole ligation, conserved in all characterized Ty, has functional advantage beyond active-site flexibility. Substrate accessibility to the oxygenated Cu2 O2 core and stabilization of a high oxidation state of the copper centers are suggested from these minimalistic models.

Minimal Functional Sites in Metalloproteins and Their Usage in Structural Bioinformatics

Metal ions play a functional role in numerous biochemical processes and cellular pathways. Indeed, about 40% of all enzymes of known 3D structure require a metal ion to be able to perform catalysis. The interactions of the metals with the macromolecular framework determine their chemical properties and reactivity. The relevant interactions involve both the coordination sphere of the metal ion and the more distant interactions of the so-called second sphere, i.e., the non-bonded interactions between the macromolecule and the residues coordinating the metal (metal ligands). The metal ligands and the residues in their close spatial proximity define what we call a minimal functional site (MFS). MFSs can be automatically extracted from the 3D structures of metal-binding biological macromolecules deposited in the Protein Data Bank (PDB). They are 3D templates that describe the local environment around a metal ion or metal cofactor and do not depend on the overall macromolecular structure. MFSs provide a different view on metal-binding proteins and nucleic acids, completely focused on the metal. Here we present different protocols and tools based upon the concept of MFS to obtain deeper insight into the structural and functional properties of metal-binding macromolecules. We also show that structure conservation of MFSs in metalloproteins relates to local sequence similarity more strongly than to overall protein similarity.

Statistical investigation of surface bound ions and further development of BION server to include pH and salt dependence

Ions are engaged in multiple biological processes in cells. By binding to the macromolecules or being mobile in the solvent, they maintain the integrity of the structure of macromolecules; participate in their enzymatic activity; or screen electrostatic interactions. While experimental methods are not always able to assign the exact location of ions, computational methods are in demand. Although the majority of computational methods are successful in predicting the position of ions buried inside macromolecules, they are less effective in deciphering positions of surface bound ions. Here, we propose the new BION algorithm (http://compbio.clemson.edu/bion_server_ph/) that predicts the location of the surface bound ions. It is more efficient and accurate compared to the previous version since it uses more advanced clustering algorithm in combination with pairing rules. In addition, the BION webserver allows specifying the pH and the salt concentration in predicting ions positions.

Structural principles for computational and de novo design of 4Fe-4S metalloproteins

Iron-sulfur centers in metalloproteins can access multiple oxidation states over a broad range of potentials, allowing them to participate in a variety of electron transfer reactions and serving as catalysts for high-energy redox processes. The nitrogenase FeMoCO cluster converts di-nitrogen to ammonia in an eight-electron transfer step. The 2(Fe4S4) containing bacterial ferredoxin is an evolutionarily ancient metalloprotein fold and is thought to be a primordial progenitor of extant oxidoreductases. Controlling chemical transformations mediated by iron-sulfur centers such as nitrogen fixation, hydrogen production as well as electron transfer reactions involved in photosynthesis are of tremendous importance for sustainable chemistry and energy production initiatives. As such, there is significant interest in the design of iron-sulfur proteins as minimal models to gain fundamental understanding of complex natural systems and as lead-molecules for industrial and energy applications. Herein, we discuss salient structural characteristics of natural iron-sulfur proteins and how they guide principles for design. Model structures of past designs are analyzed in the context of these principles and potential directions for enhanced designs are presented, and new areas of iron-sulfur protein design are proposed. This article is part of a Special issue entitled Biodesign for Bioenergetics--the design and engineering of electronic transfer cofactors, protein networks, edited by Ronald L. Koder and J.L Ross Anderson.

Methionine

This review focuses on the steps unique to methionine biosynthesis, namely the conversion of homoserine to methionine. The past decade has provided a wealth of information concerning the details of methionine metabolism and the review focuses on providing a comprehensive overview of the field, emphasizing more recent findings. Details of methionine biosynthesis are addressed along with key cellular aspects, including regulation, uptake, utilization, AdoMet, the methyl cycle, and growing evidence that inhibition of methionine biosynthesis occurs under stressful cellular conditions. The first unique step in methionine biosynthesis is catalyzed by the metA gene product, homoserine transsuccinylase (HTS, or homoserine O-succinyltransferase). Recent experiments suggest that transcription of these genes is indeed regulated by MetJ, although the repressor-binding sites have not yet been verified. Methionine also serves as the precursor of S-adenosylmethionine, which is an essential molecule employed in numerous biological processes. S-adenosylhomocysteine is produced as a consequence of the numerous AdoMet-dependent methyl transfer reactions that occur within the cell. In E. coli and Salmonella, this molecule is recycled in two discrete steps to complete the methyl cycle. Cultures challenged by oxidative stress appear to experience a growth limitation that depends on methionine levels. E. coli that are deficient for the manganese and iron superoxide dismutases (the sodA and sodB gene products, respectively) require the addition of methionine or cysteine for aerobic growth. Modulation of methionine levels in response to stressful conditions further increases the complexity of its regulation.

Insights into the posttranslational assembly of the Mo-, S- and Cu-containing cluster in the active site of CO dehydrogenase of Oligotropha carboxidovorans

Oligotropha carboxidovorans is characterized by the aerobic chemolithoautotrophic utilization of CO. CO oxidation by CO dehydrogenase proceeds at a unique bimetallic [CuSMoO2] cluster which matures posttranslationally while integrated into the completely folded apoenzyme. Kanamycin insertional mutants in coxE, coxF and coxG were characterized with respect to growth, expression of CO dehydrogenase, and the type of metal center present. These data along with sequence information were taken to delineate a model of metal cluster assembly. Biosynthesis starts with the MgATP-dependent, reductive sulfuration of [Mo(VI)O3] to [Mo(V)O2SH] which entails the AAA+-ATPase chaperone CoxD. Then Mo(V) is reoxidized and Cu(1+)-ion is integrated. Copper is supplied by the soluble CoxF protein which forms a complex with the membrane-bound von Willebrand protein CoxE through RGD-integrin interactions and enables the reduction of CoxF-bound Cu(2+), employing electrons from respiration. Copper appears as Cu(2+)-phytate, is mobilized through the phytase activity of CoxF and then transferred to the CoxF putative copper-binding site. The coxG gene does not participate in the maturation of the bimetallic cluster. Mutants in coxG retained the ability to utilize CO, although at a lower growth rate. They contained a regular CO dehydrogenase with a functional catalytic site. The presence of a pleckstrin homology (PH) domain on CoxG and the observed growth rates suggest a role of the PH domain in recruiting CO dehydrogenase to the cytoplasmic membrane enabling electron transfer from the enzyme to the respiratory chain. CoxD, CoxE and CoxF combine motifs of a DEAD-box RNA helicase which would explain their mutual translation.

Ion binding to biological macromolecules

Biological macromolecules carry out their functions in water and in the presence of ions. The ions can bind to the macromolecules either specifically or non-specifically, or can simply to be a part of the water phase providing physiological gradient across various membranes. This review outlines the differences between specific and non-specific ion binding in terms of the function and stability of the corresponding macromolecules. Furthermore, the experimental techniques to identify ion positions and computational methods to predict ion binding are reviewed and their advantages compared. It is indicated that specifically bound ions are relatively easier to be revealed while non-specifically associated ions are difficult to predict. In addition, the binding and the residential time of non-specifically bound ions are very much sensitive to the environmental factors in the cells, specifically to the local pH and ion concentration. Since these characteristics differ among the cellular compartments, the non-specific ion binding must be investigated with respect to the sub-cellular localization of the corresponding macromolecule.

Copper-peptide complex structure and reactivity when found in conserved His-X(aa)-His sequences

Oxygen-activating copper proteins may possess His-X(aa)-His chelating sequences at their active sites and additionally exhibit imidiazole group δN vs εN tautomeric preferences. As shown here, such variations strongly affect copper ion's coordination geometry, redox behavior, and oxidative reactivity. Copper(I) complexes bound to either δ-HGH or ε-HGH tripeptides were synthesized and characterized. Structural investigations using X-ray absorption spectroscopy, density functional theory calculations, and solution conductivity measurements reveal that δ-HGH forms the Cu(I) dimer complex [{Cu(I)(δ-HGH)}2](2+) (1) while ε-HGH binds Cu(I) to give the monomeric complex [Cu(I)(ε-HGH)](+) (2). Only 2 exhibits any reactivity, forming a strong CO adduct, [Cu(I)(ε-HGH)(CO)](+), with properties closely matching those of the copper monooxygenase PHM. Also, 2 is reactive toward O2 or H2O2, giving a new type of O2-adduct or Cu(II)-OOH complex, respectively.

Beyond the active site: the impact of the outer coordination sphere on electrocatalysts for hydrogen production and oxidation

Redox active metalloenzymes play a major role in energy transformation reactions in biological systems. Examples include formate dehydrogenases, nitrogenases, CO dehydrogenase, and hydrogenases. Many of these reactions are also of interest to humans as potential energy storage or utilization reactions for photoelectrochemical, electrolytic, and fuel cell applications. These metalloenzymes consist of redox active metal centers where substrates are activated and undergo transformation to products accompanied by electron and proton transfer to or from the substrate. These active sites are typically buried deep within a protein matrix of the enzyme with channels for proton transport, electron transport, and substrate/product transport between the active site and the surface of the protein. In addition, there are amino acid residues that lie in close proximity to the active site that are thought to play important roles in regulating and enhancing enzyme activity. Directly studying the outer coordination sphere of enzymes can be challenging due to their complexity, and the use of modified molecular catalysts may allow us to provide some insight. There are two fundamentally different approaches to understand these important interactions. The "bottom-up" approach involves building an amino acid or peptide containing outer coordination sphere around a functional molecular catalyst, and the "top-down" approach involves attaching molecular catalyst to a structured protein. Both of these approaches have been undertaken for hydrogenase mimics and are the emphasis of this Account. Our focus has been to utilize amino acid or peptide based scaffolds on an active functional enzyme mimic for H2 oxidation and production, [Ni(P(R)2N(R('))2)2](2+). This "bottom-up" approach has allowed us to evaluate individual functional group and structural contributions to electrocatalysts for H2 oxidation and production. For instance, using amine, ether, and carboxylic acid functionalities in the outer coordination sphere enhances proton movement and results in lower catalytic overpotentials for H2 oxidation, while achieving water solubility in some cases. Amino acids with acidic and basic side chains concentrate substrate around catalysts for H2 production, resulting in up to 5-fold enhancements in rate. The addition of a structured peptide in an H2 production catalyst limited the structural freedom of the amino acids nearest the active site, while enhancing the overall rate. Enhanced stability to oxygen or extreme conditions such as strongly acidic or basic conditions has also resulted from an amino acid based outer coordination sphere. From the "top-down" approach, others have achieved water solubility and photocatalytic activity by associating this core complex with photosystem-I. Collectively, by use of this well understood core, the role of individual and combined features of the outer coordination sphere are starting to be understood at a mechanistic level. Common mechanisms have yet to be defined to predictably control these processes, but our growing knowledge in this area is essential for the eventual mimicry of enzymes by efficient molecular catalysts for practical use.

Latex clearing protein-an oxygenase cleaving poly(cis-1,4-isoprene) rubber at the cis double bonds

Gordonia polyisoprenivorans strain VH2, a potent rubber-degrading actinomycete, harbors two latex clearing proteins (Lcps), which are known to be essential for the microbial degradation of rubber. However, biochemical information on the exact role of this protein in the degradation of polyisoprene was lacking. In this study, the gene encoding Lcp1VH2 was heterologously expressed in strains of Escherichia coli, the corresponding protein was purified, and its role in rubber degradation was examined by measurement of oxygen consumption as well as by chromatographic and spectroscopic methods. It turned out that active Lcp1VH2 is a monomer and is responsible for the oxidative cleavage of poly(cis-1,4-isoprene) in synthetic as well as in natural rubber by the addition of oxygen (O2) to the cis double bonds. The resulting oligomers possess repetitive isoprene units with aldehyde (CHO-CH2-) and ketone (-CH2-CO-CH3) functional groups at the termini. Two fractions with average isoprene contents of 18 and 10, respectively, were isolated, thus indicating an endocleavage mechanism. The activity of Lcp1VH2 was determined by applying a polarographic assay. Alkenes, acyclic terpenes, or other rubber-like polymers, such as poly(cis-1,4-butadiene) or poly(trans-1,4-isoprene), are not oxidatively cleaved by Lcp1VH2. The pH and temperature optima of the enzyme are at pH 7 and 30°C, respectively. Furthermore, it was demonstrated that active Lcp1VH2 is a Cu(II)-containing oxygenase that exhibits a conserved domain of unknown function which cannot be detected in any other hitherto-characterized enzyme. The results presented here indicate that this domain might represent a new protein family of oxygenases.

Analyses of cobalt-ligand and potassium-ligand bond lengths in metalloproteins: trends and patterns

Cobalt and potassium are biologically important metal elements that are present in a large array of proteins. Cobalt is mostly found in vivo associated with a corrin ring, which represents the core of the vitamin B12 molecule. Potassium is the most abundant metal in the cytosol, and it plays a crucial role in maintaining membrane potential as well as correct protein function. Here, we report a thorough analysis of the geometric properties of cobalt and potassium coordination spheres that was performed with high resolution on a representative set of structures from the Protein Data Bank and complemented by quantum mechanical calculations realized at the DFT level of theory (B3LYP/ SDD) on mononuclear model systems. The results allowed us to draw interesting conclusions on the structural characteristics of both Co and K centers, and to evaluate the importance of effects such as their association energies and intrinsic thermodynamic stabilities. Overall, the results obtained provide useful data for enhancing the atomic models normally applied in theoretical and computational studies of Co or K proteins performed at the quantum mechanical level, and for developing molecular mechanical parameters for treating Co or K coordination spheres in molecular mechanics or molecular dynamics studies.

Putative cobalt- and nickel-binding proteins and motifs in Streptococcus pneumoniae

Cobalt and nickel play important roles in various biological processes. The present work focuses on the enrichment and identification of Co- and Ni-binding motifs and proteins in Gram-positive bacteria. Immobilized metal affinity column (IMAC) was used to partially enrich putative metal-binding proteins and peptides from Streptococcus pneumoniae, and then LTQ-Orbitrap mass spectrometry (MS) was applied to identify and characterize the metal-binding motifs and proteins. In total, 208 and 223 proteins were isolated by Co- and Ni-IMAC columns respectively, in which 129 proteins were present in both preparations. Based on the gene ontology (GO) analysis, the putative metal-binding proteins were found to be mainly involved in protein metabolism, gene expression regulation and carbohydrate metabolism. These putative metal-binding proteins form a highly connected network, indicating that they may synergistically work together to achieve specific biological functions. Putative Co- and Ni-binding motifs were identified with H(X)nH, M(X)nH and H(X)nM derived from the identified 51 Co-binding peptides and 66 Ni-binding peptides. Statistics of frequency of amino acids in the metal-binding motifs showed that cobalt and nickel prefer to bind histidine and methionine, but not cysteine. These results obtained by a systematic metalloproteomic approach provide important clues for the further investigation of metal homeostasis and metal-related virulence of bacteria.

Enzyme design from the bottom up: an active nickel electrocatalyst with a structured peptide outer coordination sphere

Catalytic, peptide-containing metal complexes with a well-defined peptide structure have the potential to enhance molecular catalysts through an enzyme-like outer coordination sphere. Here, we report the synthesis and characterization of an active, peptide-based metal complex built upon the well-characterized hydrogen production catalyst [Ni(P(Ph)2N(Ph))2](2+) (P(Ph)2N(Ph)=1,3,6-triphenyl-1-aza-3,6-diphosphacycloheptane). The incorporated peptide maintains its β-hairpin structure when appended to the metal core, and the electrocatalytic activity of the peptide-based metal complex (≈100,000 s(-1)) is enhanced compared to the parent complex ([Ni(P(Ph)2N(APPA))2](2+); ≈50,500 s(-1)). The combination of an active molecular catalyst with a structured peptide provides a scaffold that permits the incorporation of features of an enzyme-like outer-coordination sphere necessary to create molecular electrocatalysts with enhanced functionality.

Role of a remote leucine residue in the catalytic function of polyol dehydrogenase

This study examined the role of remote residues on the structure and function of zinc-dependent polyol dehydrogenases.

Minimal proton channel enables H2 oxidation and production with a water-soluble nickel-based catalyst

Hydrogenase enzymes use first-row transition metals to interconvert H2 with protons and electrons, reactions that are important for the storage and recovery of energy from intermittent sources such as solar, hydroelectric, and wind. Here we present Ni(P(Cy)2N(Gly)2)2, a water-soluble molecular electrocatalyst with the amino acid glycine built into the diphosphine ligand framework. Proton transfer between the outer coordination sphere carboxylates and the second coordination sphere pendant amines is rapid, as observed by cyclic voltammetry and FTIR spectroscopy, indicating that the carboxylate groups may participate in proton transfer during catalysis. This complex oxidizes H2 (1-33 s(-1)) at low overpotentials (150-365 mV) over a range of pH values (0.1-9.0) and produces H2 under identical solution conditions (>2400 s(-1) at pH 0.5). Enzymes employ proton channels for the controlled movement of protons over long distances-the results presented here demonstrate the effects of a simple two-component proton channel in a synthetic molecular electrocatalyst.

MetalS2: a tool for the structural alignment of minimal functional sites in metal-binding proteins and nucleic acids

We developed a new software tool, MetalS(2), for the structural alignment of Minimal Functional Sites (MFSs) in metal-binding biological macromolecules. MFSs are 3D templates that describe the local environment around the metal(s) independently of the larger context of the macromolecular structure. Such local environment has a determinant role in tuning the chemical reactivity of the metal, ultimately contributing to the functional properties of the whole system. On our example data sets, MetalS(2) unveiled structural similarities that other programs for protein structure comparison do not consistently point out and overall identified a larger number of structurally similar MFSs. MetalS(2) supports the comparison of MFSs harboring different metals and/or with different nuclearity and is available both as a stand-alone program and a Web tool ( http://metalweb.cerm.unifi.it/tools/metals2/).

A modular, energy-based approach to the development of nickel containing molecular electrocatalysts for hydrogen production and oxidation

This review discusses the development of molecular electrocatalysts for H2 production and oxidation based on nickel. A modular approach is used in which the structure of the catalyst is divided into first, second, and outer coordination spheres. The first coordination sphere consists of the ligands bound directly to the metal center, and this coordination sphere can be used to control such factors as the presence or absence of vacant coordination sites, redox potentials, hydride donor abilities and other important thermodynamic parameters. The second coordination sphere includes functional groups such as pendent acids or bases that can interact with bound substrates such as H2 molecules and hydride ligands, but that do not form strong bonds with the metal center. These functional groups can play diverse roles such as assisting the heterolytic cleavage of H2, controlling intra- and intermolecular proton transfer reactions, and providing a physical pathway for coupling proton and electron transfer reactions. By controlling both the hydride donor ability of the catalysts using the first coordination sphere and the proton donor abilities of the functional groups in the second coordination sphere, catalysts can be designed that are biased toward H2 production, oxidation, or bidirectional (catalyzing both H2 oxidation and production). The outer coordination sphere is defined as that portion of the catalytic system that is beyond the second coordination sphere. This coordination sphere can assist in the delivery of protons and electrons to and from the catalytically active site, thereby adding another important avenue for controlling catalytic activity. Many features of these simple catalytic systems are good models for enzymes, and these simple systems provide insights into enzyme function and reactivity that may be difficult to probe in enzymes. This article is part of a Special Issue entitled: Metals in Bioenergetics and Biomimetics Systems.

The role of a dipeptide outer-coordination sphere on H2-production catalysts: influence on catalytic rates and electron transfer

The outer-coordination sphere of enzymes acts to fine-tune the active site reactivity and control catalytic rates, suggesting that incorporation of analogous structural elements into molecular catalysts may be necessary to achieve rates comparable to those observed in enzyme systems at low overpotentials. In this work, we evaluate the effect of an amino acid and dipeptide outer-coordination sphere on [Ni(P(Ph)(2)N(Ph-R)(2))(2)](2+) hydrogen production catalysts. A series of 12 new complexes containing non-natural amino acids or dipeptides was prepared to test the effects of positioning, size, polarity and aromaticity on catalytic activity. The non-natural amino acid was either 3-(meta- or para-aminophenyl)propionic acid terminated as an acid, an ester or an amide. Dipeptides consisted of one of the non-natural amino acids coupled to one of four amino acid esters: alanine, serine, phenylalanine or tyrosine. All of the catalysts are active for hydrogen production, with rates averaging ∼1000 s(-1), 40 % faster than the unmodified catalyst. Structure and polarity of the aliphatic or aromatic side chains of the C-terminal peptide do not strongly influence rates. However, the presence of an amide bond increases rates, suggesting a role for the amide in assisting catalysis. Overpotentials were lower with substituents at the N-phenyl meta position. This is consistent with slower electron transfer in the less compact, para-substituted complexes, as shown in digital simulations of catalyst cyclic voltammograms and computational modeling of the complexes. Combining the current results with insights from previous results, we propose a mechanism for the role of the amino acid and dipeptide based outer-coordination sphere in molecular hydrogen production catalysts.

A new vector for identification of prokaryotes and their variable-size genomes

A large number of prokaryotes have been produced, so how to provide a means to describe and distinguish them accurately is becoming a key issue of prokaryotic taxonomy. We proposed an efficient algorithm to filter out most genome fragments that are horizontally transferred, and extracted a new genome vector (GV). To highlight the power of GV, we applied it to identify prokaryotes and their variable-size genome fragments. The result indicated that the new vector as species tags can accurately identify genome fragments as short as 3,000 bp at species level.

Predicting nonspecific ion binding using DelPhi

Ions are an important component of the cell and affect the corresponding biological macromolecules either via direct binding or as a screening ion cloud. Although some ion binding is highly specific and frequently associated with the function of the macromolecule, other ions bind to the protein surface nonspecifically, presumably because the electrostatic attraction is strong enough to immobilize them. Here, we test such a scenario and demonstrate that experimentally identified surface-bound ions are located at a potential that facilitates binding, which indicates that the major driving force is the electrostatics. Without taking into consideration geometrical factors and structural fluctuations, we show that ions tend to be bound onto the protein surface at positions with strong potential but with polarity opposite to that of the ion. This observation is used to develop a method that uses a DelPhi-calculated potential map in conjunction with an in-house-developed clustering algorithm to predict nonspecific ion-binding sites. Although this approach distinguishes only the polarity of the ions, and not their chemical nature, it can predict nonspecific binding of positively or negatively charged ions with acceptable accuracy. One can use the predictions in the Poisson-Boltzmann approach by placing explicit ions in the predicted positions, which in turn will reduce the magnitude of the local potential and extend the limits of the Poisson-Boltzmann equation. In addition, one can use this approach to place the desired number of ions before conducting molecular-dynamics simulations to neutralize the net charge of the protein, because it was shown to perform better than standard screened Coulomb canned routines, or to predict ion-binding sites in proteins. This latter is especially true for proteins that are involved in ion transport, because such ions are loosely bound and very difficult to detect experimentally.

Minimal functional sites allow a classification of zinc sites in proteins

Zinc is indispensable to all forms of life as it is an essential component of many different proteins involved in a wide range of biological processes. Not differently from other metals, zinc in proteins can play different roles that depend on the features of the metal-binding site. In this work, we describe zinc sites in proteins with known structure by means of three-dimensional templates that can be automatically extracted from PDB files and consist of the protein structure around the metal, including the zinc ligands and the residues in close spatial proximity to the ligands. This definition is devised to intrinsically capture the features of the local protein environment that can affect metal function, and corresponds to what we call a minimal functional site (MFS). We used MFSs to classify all zinc sites whose structures are available in the PDB and combined this classification with functional annotation as available in the literature. We classified 77% of zinc sites into ten clusters, each grouping zinc sites with structures that are highly similar, and an additional 16% into seven pseudo-clusters, each grouping zinc sites with structures that are only broadly similar. Sites where zinc plays a structural role are predominant in eight clusters and in two pseudo-clusters, while sites where zinc plays a catalytic role are predominant in two clusters and in five pseudo-clusters. We also analyzed the amino acid composition of the coordination sphere of zinc as a function of its role in the protein, highlighting trends and exceptions. In a period when the number of known zinc proteins is expected to grow further with the increasing awareness of the cellular mechanisms of zinc homeostasis, this classification represents a valuable basis for structure-function studies of zinc proteins, with broad applications in biochemistry, molecular pharmacology and de novo protein design.

A comparative molecular dynamics study of methylation state specificity of JMJD2A

Histone modifications have great importance in epigenetic regulation. JMJD2A is a histone demethylase which is selective for di- and trimethyl forms of residues Lys9 and Lys36 of Histone 3 tail (H3K9 and H3K36). We present a molecular dynamics simulations of mono-, di- and trimethylated histone tails in complex with JMJD2A catalytic domain to gain insight into how JMJD2A discriminates between the methylation states of H3K9. The methyl groups are located at specific distances and orientations with respect to Fe(II) in methylammonium binding pocket. For the trimethyllysine the mechanism which provides the effectual orientation of methyl groups is the symmetry, whereas for the dimethyllysine case the determining factors are the interactions between methyllysine head and its environment and subsequently the restriction on angular motion. The occurrence frequency of methyl groups in a certain proximity of Fe(II) comes out as the explanation of the enzyme activity difference on di- and tri-methylated peptides. Energy analysis suggests that recognition is mostly driven by van der Waals and followed by Coulombic interactions in the enzyme-substrate interface. The number (mono, di or tri) and orientations of methyl groups and water molecules significantly affect the extent of van der Waals interaction strengths. Hydrogen bonding analysis suggests that the interaction between JMJD2A and its substrates mainly comes from main chain-side chain interactions. Binding free energy analysis points out Arg8 as an important residue forming an intra-substrate hydrogen bond with tri and dimethylated Lys9 of the H3 chain. Our study provides new insights into how JMJD2A discriminates between its substrates from both a structural and dynamical point of view.

Putative copper- and zinc-binding motifs in Streptococcus pneumoniae identified by immobilized metal affinity chromatography and mass spectrometry

The aim of metalloproteomics is to identify and characterize putative metal-binding proteins and metal-binding motifs. In this study, we performed a systematical metalloproteomic analysis on Streptococcus pneumoniae through the combined use of efficient immobilized metal affinity chromatography enrichment and high-accuracy linear ion trap-Orbitrap MS to identify metal-binding proteins and metal-binding peptides. In total, 232 and 166 putative metal-binding proteins were respectively isolated by Cu- and Zn-immobilized metal affinity chromatography columns, in which 133 proteins were present in both preparations. The putative metalloproteins are mainly involved in protein, nucleotide and carbon metabolisms, oxidation and cell cycle regulation. Based on the sequence of the putative Cu- and Zn-binding peptides, putative Cu-binding motifs were identified: H(X)mH (m=0-11), C(X)(2) C, C(X)nH (n=2-4, 6, 9), H(X)iM (i=0-10) and M(X)tM (t=8 or 12), while putative Zn-binding motifs were identified as follows: H(X)mH (m=1-12), H(X)iM (i=0-12), M(X)tM (t=0, 3 and 4), C(X)nH (n=1, 2, 7, 10 and 11). Equilibrium dialysis and inductively coupled plasma-MS experiments confirmed that the artificially synthesized peptides harboring differential identified metal-binding motifs interacted directly with the metal ions. The metalloproteomic study presented here suggests that the comparably large size and diverse functions of the S. pneumoniae metalloproteome may play important roles in various biological processes and thus contribute to the bacterial pathologies.