Changes in hydrogen-bond strengths explain reduction potentials in 10 rubredoxin variants

Role of hydrogen-bond networks on the donor side of photosynthetic reaction centers from purple bacteria

For the last decades, significant progress has been made in studying the biological functions of H-bond networks in membrane proteins, proton transporters, receptors, and photosynthetic reaction centers. Increasing availability of the X-ray crystal and cryo-electron microscopy structures of photosynthetic complexes resolved with high atomic resolution provides a platform for their comparative analysis. It allows identifying structural factors that are ensuring the high quantum yield of the photochemical reactions and are responsible for the stability of the membrane complexes. The H-bond networks are known to be responsible for proton transport associated with electron transfer from the primary to the secondary quinone as well as in the processes of water oxidation in photosystem II. Participation of such networks in reactions proceeding on the periplasmic side of bacterial photosynthetic reaction centers is less studied. This review summarizes the current understanding of the role of H-bond networks on the donor side of photosynthetic reaction centers from purple bacteria. It is discussed that the networks may be involved in providing close association with mobile electron carriers, in light-induced proton transport, in regulation of the redox properties of bacteriochlorophyll cofactors, and in stabilization of the membrane protein structure at the interface of membrane and soluble phases.

Designed Rubredoxin miniature in a fully artificial electron chain triggered by visible light

Designing metal sites into de novo proteins has significantly improved, recently. However, identifying the minimal coordination spheres, able to encompass the necessary information for metal binding and activity, still represents a great challenge, today. Here, we test our understanding with a benchmark, nevertheless difficult, case. We assemble into a miniature 28-residue protein, the quintessential elements required to fold properly around a FeCys₄ redox center, and to function efficiently in electron-transfer. This study addresses a challenge in de novo protein design, as it reports the crystal structure of a designed tetra-thiolate metal-binding protein in sub-Å agreement with the intended design. This allows us to well correlate structure to spectroscopic and electrochemical properties. Given its high reduction potential compared to natural and designed FeCys₄-containing proteins, we exploit it as terminal electron acceptor of a fully artificial chain triggered by visible light.

Ensemble-function relationships to dissect mechanisms of enzyme catalysis

Decades of structure-function studies have established our current extensive understanding of enzymes. However, traditional structural models are snapshots of broader conformational ensembles of interchanging states. We demonstrate the need for conformational ensembles to understand function, using the enzyme ketosteroid isomerase (KSI) as an example. Comparison of prior KSI cryogenic x-ray structures suggested deleterious mutational effects from a misaligned oxyanion hole catalytic residue. However, ensemble information from room-temperature x-ray crystallography, combined with functional studies, excluded this model. Ensemble-function analyses can deconvolute effects from altering the probability of occupying a state (P-effects) and changing the reactivity of each state (k-effects); our ensemble-function analyses revealed functional effects arising from weakened oxyanion hole hydrogen bonding and substrate repositioning within the active site. Ensemble-function studies will have an integral role in understanding enzymes and in meeting the future goals of a predictive understanding of enzyme catalysis and engineering new enzymes.

Revisiting the concept of peptide bond planarity in an iron-sulfur protein by neutron structure analysis

The planarity of the peptide bond is important for the stability and structure formation of proteins. However, substantial distortion of peptide bonds has been reported in several high-resolution structures and computational analyses. To investigate the peptide bond planarity, including hydrogen atoms, we report a 1.2-Å resolution neutron structure of the oxidized form of high-potential iron-sulfur protein. This high-resolution neutron structure shows that the nucleus positions of the amide protons deviate from the peptide plane and shift toward the acceptors. The planarity of the H─N─C═O plane depends strongly on the pyramidalization of the nitrogen atom. Moreover, the orientation of the amide proton of Cys⁷⁵ is different in the reduced and oxidized states, possibly because of the electron storage capacity of the iron-sulfur cluster.

X-ray structure of the Rhodobacter sphaeroides reaction center with an M197 Phe→His substitution clarifies the properties of the mutant complex

The first steps of the global process of photosynthesis take place in specialized membrane pigment-protein complexes called photosynthetic reaction centers (RCs). The RC of the photosynthetic purple bacterium Rhodobacter sphaeroides, a relatively simple analog of the more complexly organized photosystem II in plants, algae and cyanobacteria, serves as a convenient model for studying pigment-protein interactions that affect photochemical processes. In bacterial RCs the bacteriochlorophyll (BChl) dimer P serves as the primary electron donor, and its redox potential is a critical factor in the efficient functioning of the RC. It has previously been shown that the replacement of Phe M197 by His strongly affects the oxidation potential of P (E _m P/P⁺), increasing its value by 125 mV, as well as increasing the thermal stability of RC and its stability in response to external pressure. The crystal structures of F(M197)H RC at high resolution obtained using various techniques presented in this report clarify the optical and electrochemical properties of the primary electron donor and the increased resistance of the mutant complex to denaturation. The electron-density maps are consistent with the donation of a hydrogen bond from the imidazole group of His M197 to the C2-acetyl carbonyl group of BChl P_B. The formation of this hydrogen bond leads to a considerable out-of-plane rotation of the acetyl carbonyl group and results in a 1.2 Å shift of the O atom of this group relative to the wild-type structure. Besides, the distance between BChl P_A and P_B in the area of pyrrole ring I was found to be increased by up to 0.17 Å. These structural changes are discussed in association with the spectral properties of BChl dimer P. The electron-density maps strongly suggest that the imidazole group of His M197 accepts another hydrogen bond from the nearest water molecule, which in turn appears to form two more hydrogen bonds to Asn M195 and Asp L155. As a result of the F(M197)H mutation, BChl P_B finds itself connected to the extensive hydrogen-bonding network that pre-existed in wild-type RC. Dissimilarities in the two hydrogen-bonding networks near the M197 and L168 sites may account for the different changes of the E _m P/P⁺ in F(M197)H and H(L168)F RCs. The involvement of His M197 in the hydrogen-bonding network also appears to be related to stabilization of the F(M197)H RC structure. Analysis of the experimental data presented here and of the data available in the literature points to the fact that the hydrogen-bonding networks in the vicinity of BChl dimer P may play an important role in fine-tuning the redox properties of the primary electron donor.

Electrochemical characterization of mutant forms of rubredoxin B from Mycobacterium tuberculosis

Electron transfer in metalloproteins is a driving force for many biological processes and widely distributed in nature. Rubredoxin B (RubB) from Mycobacterium tuberculosis is a first example among [1Fe-0S] proteins that support catalytic activity of terminal sterol-monooxygenases enabling its application in metabolic engineering. To explore the tolerance of RubB to the specific amino acid changes we evaluated the effect of surface mutations on its electrochemical properties. Based on the RubB fold we also designed the mutant with a putative additional site for protein-protein interactions to further evaluate electron transfer and electrochemical properties. The investigation of redox properties of mutant variants of RubB was done using screen-printed graphite electrodes (SPEs) modified with stable dispersion of multi-walled carbon nanotubes (MWCNTs). The redox potentials (midpoint potentials, E^0Ꞌ) of mutants did not significantly differ from the wild type protein and vary in the range of -264 to -231 mV vs. Ag/AgCl electrode. However, all mutations affect electron transfer rate between the protein and electrode. Notably, the modulation of the protein-protein interactions was observed for the insertion mutant suggesting the possibility of tailoring of rubredoxin for the selected redox-partner. Overall, RubB is tolerant to the significant modifications in its structure enabling rational engineering of novel redox proteins.

The role of thermodynamic features on the functional activity of electron bifurcating enzymes

Electron bifurcation is a biological mechanism to drive a thermodynamically unfavorable redox reaction through direct coupling with an exergonic reaction. This process allows microorganisms to generate high energy reducing equivalents in order to sustain life and is often found in anaerobic metabolism, where the energy economy of the cell is poor. Recent work has revealed details of the redox energy landscapes for a variety of electron bifurcating enzymes, greatly expanding the understanding of how energy is transformed by this unique mechanism. Here we highlight the plasticity of these emerging landscapes, what is known regarding their mechanistic underpinnings, and provide a context for interpreting their biochemical activity within the physiological framework. We conclude with an outlook for propelling the field toward an integrative understanding of the impact of electron bifurcation.

PRE-driven protein NMR structures: an alternative approach in highly paramagnetic systems

Metalloproteins play key roles across biology, and knowledge of their structure is essential to understand their physiological role. For those metalloproteins containing paramagnetic states, the enhanced relaxation caused by the unpaired electrons often makes signal detection unfeasible near the metal center, precluding adequate structural characterization right where it is more biochemically relevant. Here, we report a protein structure determination by NMR where two different sets of restraints, one containing Nuclear Overhauser Enhancements (NOEs) and another containing Paramagnetic Relaxation Enhancements (PREs), are used separately and eventually together. The protein PioC from Rhodopseudomonas palustris TIE-1 is a High Potential Iron-Sulfur Protein (HiPIP) where the [4Fe-4S] cluster is paramagnetic in both oxidation states at room temperature providing the source of PREs used as alternative distance restraints. Comparison of the family of structures obtained using NOEs only, PREs only, and the combination of both reveals that the pairwise root-mean-square deviation (RMSD) between them is similar and comparable with the precision within each family. This demonstrates that, under favorable conditions in terms of protein size and paramagnetic effects, PREs can efficiently complement and eventually replace NOEs for the structural characterization of small paramagnetic metalloproteins and de novo-designed metalloproteins by NMR. DATABASES: The 20 conformers with the lowest target function constituting the final family obtained using the full set of NMR restraints were deposited to the Protein Data Bank (PDB ID: 6XYV). The 20 conformers with the lowest target function obtained using NOEs only (PDB ID: 7A58) and PREs only (PDB ID: 7A4L) were also deposited to the Protein Data Bank. The chemical shift assignments were deposited to the BMRB (code 34487).

1H, 13C and 15N assignment of the paramagnetic high potential iron-sulfur protein (HiPIP) PioC from Rhodopseudomonas palustris TIE-1

High potential iron-sulfur proteins (HiPIPs) are a class of small proteins (50-100 aa residues), containing a 4Fe-4S iron-sulfur cluster. The 4Fe-4S cluster shuttles between the oxidation states [Fe₄S₄]^3+/2+, with a positive redox potential in the range (500-50 mV) throughout the different known HiPIPs. Both oxidation states are paramagnetic at room temperature. HiPIPs are electron transfer proteins, isolated from photosynthetic bacteria and usually provide electrons to the photosynthetic reaction-center. PioC, the HIPIP isolated from Rhodopseudomonas palustris TIE-1, is the smallest among all known HiPIPs. Despite their small dimensions, an extensive NMR assignment is only available for two of them, because paramagnetism prevents the straightforward assignment of all resonances. We report here the complete NMR assignment of ¹H, ¹³C and ¹⁵N signals for the reduced [Fe₄S₄]²⁺ state of the protein. A set of double and triple resonance experiments performed with standardized parameters/datasets provided the assignment of about 72% of the residues. The almost complete resonance assignment (99.5% of backbone and ca. 90% of side chain resonances) was achieved by combining the above information with those obtained using a second set of NMR experiments, in which acquisition and processing parameters, as well as pulse sequences design, were optimized to account for the peculiar features of this paramagnetic protein.

Paramagnetic NMR Spectroscopy Is a Tool to Address Reactivity, Structure, and Protein–Protein Interactions of Metalloproteins: The Case of Iron–Sulfur Proteins

The study of cellular machineries responsible for the iron–sulfur (Fe–S) cluster biogenesis has led to the identification of a large number of proteins, whose importance for life is documented by an increasing number of diseases linked to them. The labile nature of Fe–S clusters and the transient protein–protein interactions, occurring during the various steps of the maturation process, make their structural characterization in solution particularly difficult. Paramagnetic nuclear magnetic resonance (NMR) has been used for decades to characterize chemical composition, magnetic coupling, and the electronic structure of Fe–S clusters in proteins; it represents, therefore, a powerful tool to study the protein–protein interaction networks of proteins involving into iron–sulfur cluster biogenesis. The optimization of the various NMR experiments with respect to the hyperfine interaction will be summarized here in the form of a protocol; recently developed experiments for measuring longitudinal and transverse nuclear relaxation rates in highly paramagnetic systems will be also reviewed. Finally, we will address the use of extrinsic paramagnetic centers covalently bound to diamagnetic proteins, which contributed over the last twenty years to promote the applications of paramagnetic NMR well beyond the structural biology of metalloproteins.

The rupture mechanism of rubredoxin is more complex than previously thought

The surprisingly low rupture force and remarkable mechanical anisotropy of rubredoxin have been known for several years. Exploiting the first combination of steered molecular dynamics and the quantum chemical Judgement of Energy DIstribution (JEDI) analysis, the common belief that hydrogen bonds between neighboring amino acid backbones and the sulfur atoms of the central FeS₄ unit in rubredoxin determine the low mechanical resistance of the protein is invalidated. The distribution of strain energy in the central part of rubredoxin is elucidated in real-time with unprecedented detail, giving important insights into the mechanical unfolding pathway of rubredoxin. While structural anisotropy as well as the contribution of angle bendings in the FeS₄ unit have a significant influence on the mechanical properties of rubredoxin, these factors are insufficient to explain the experimentally observed low rupture force. Instead, the rupture mechanism of rubredoxin is far more complex than previously thought and requires more than just a hydrogen bond network.

Analysis of RNA Methylation by Phylogenetically Diverse Cfr Radical S-Adenosylmethionine Enzymes Reveals an Iron-Binding Accessory Domain in a Clostridial Enzyme

Cfr is a radical S-adenosylmethionine (SAM) RNA methylase linked to multidrug antibiotic resistance in bacterial pathogens. It catalyzes a chemically challenging C-C bond-forming reaction to methylate C8 of A2503 (Escherichia coli numbering) of 23S rRNA during ribosome assembly. The cfr gene has been identified as a mobile genetic element in diverse bacteria and in the genome of select Bacillales and Clostridiales species. Despite the importance of Cfr, few representatives have been purified and characterized in vitro. Here we show that Cfr homologues from Bacillus amyloliquefaciens, Enterococcus faecalis, Paenibacillus lautus, and Clostridioides difficile act as C8 adenine RNA methylases in biochemical assays. C. difficile Cfr contains an additional Cys-rich C-terminal domain that binds a mononuclear Fe²⁺ ion in a rubredoxin-type Cys₄ motif. The C-terminal domain can be truncated with minimal impact on C. difficile Cfr activity, but the rate of turnover is decreased upon disruption of the Fe²⁺-binding site by Zn²⁺ substitution or ligand mutation. These findings indicate an important purpose for the observed C-terminal iron in the native fusion protein. Bioinformatic analysis of the C. difficile Cfr Cys-rich domain shows that it is widespread (∼1400 homologues) as a stand-alone gene in pathogenic or commensal Bacilli and Clostridia, with >10% encoded adjacent to a predicted radical SAM RNA methylase. Although the domain is not essential for in vitro C. difficile Cfr activity, the genomic co-occurrence and high abundance in the human microbiome suggest a possible functional role for a specialized rubredoxin in certain radical SAM RNA methylases that are relevant to human health.

g-Tensor Directions in the Protein Structural Frame of Hyperthermophilic Archaeal Reduced Rieske-Type Ferredoxin Explored by 13C Pulsed Electron Paramagnetic Resonance

Interpretation of magnetic resonance data in the context of structural and chemical biology requires prior knowledge of the g-tensor directions for paramagnetic metallo-cofactors with respect to the protein structural frame. Access to this information is often limited by the strict requirement of suitable protein crystals for single-crystal electron paramagnetic resonance (EPR) measurements or the reliance on protons (with ambiguous locations in crystal structures) near the paramagnetic metal site. Here we develop a novel pulsed EPR approach with selective ¹³C_β-cysteine labeling of model [2Fe-2S] proteins to help bypass these problems. Analysis of the ¹³C_β-cysteine hyperfine tensors reproduces the g-tensor of the Pseudomonas putida ISC-like [2Fe-2S] ferredoxin (FdxB). Its application to the hyperthermophilic archaeal Rieske-type [2Fe-2S] ferredoxin (ARF) from Sulfolobus solfataricus, for which the single-crystal EPR approach was not feasible, supports the best-fit g _x-, g _z-, and g _y-tensor directions of the reduced cluster as nearly along Fe-Fe, S-S, and the cluster plane normal, respectively. These approximate principal directions of the reduced ARF g-tensor, explored by ¹³C pulsed EPR, are less skewed from the cluster molecular axes and are largely consistent with those previously determined by single-crystal EPR for the cytochrome bc₁-associated, reduced Rieske [2Fe-2S] center. This suggests the approximate g-tensor directions are conserved across the phylogenetically and functionally divergent Rieske-type [2Fe-2S] proteins.

Hydrogen Bonds: Simple after All?

Hydrogen bonds play integral roles in biological structure, function, and conformational dynamics and are fundamental to life as it has evolved on Earth. However, our understanding of these fundamental and ubiquitous interactions has seemed fractured and incomplete, and it has been difficult to extract generalities and principles about hydrogen bonds despite thousands of papers published on this topic, perhaps in part because of the expanse of this subject and the density of studies. Fortunately, recent hydrogen bond proposals, discussions, and debates have stimulated new tests and models and have led to a remarkably simple picture of the structure of hydrogen bonds. This knowledge also provides clarity concerning hydrogen bond energetics, limiting and simplifying the factors that need be considered. Herein we recount the advances that have led to this simpler view of hydrogen bond structure, dynamics, and energetics. A quantitative predictive model for hydrogen bond length can now be broadly and deeply applied to evaluate current proposals and to uncover structural features of proteins, their conformational restraints, and their correlated motions. In contrast, a quantitative energetic description of molecular recognition and catalysis by proteins remains an important ongoing challenge, although our improved understanding of hydrogen bonds may aid in testing predictions from current and future models. We close by codifying our current state of understanding into five "Rules for Hydrogen Bonding" that may provide a foundation for understanding and teaching about these vital interactions and for building toward a deeper understanding of hydrogen bond energetics.

Metabolic shift at the class level sheds light on adaptation of methanogens to oxidative environments

Methanogens have long been considered strictly anaerobic and oxygen-sensitive microorganisms, but their ability to survive oxygen stress has also been documented. Indeed, methanogens have been found in oxidative environments, and antioxidant genes have been detected in their genomes. How methanogens adapt to oxidative environments, however, remain poorly understood. Here, we systematically predicted and annotated antioxidant features from representative genomes across six well-established methanogen orders. Based on functional gene content involved in production of reactive oxygen species, Hierarchical Clustering analyses grouped methanogens into two distinct clusters, corresponding to the Class I and II methanogens, respectively. Comparative genomics suggested a systematic shift in metabolisms across the two classes, resulting in an enrichment of antioxidant features in the Class II. Moreover, meta-analysis of 16 S rRNA gene sequences obtained from EnvDB indicated that members of Class II were more frequently recovered from microaerophilic and even oxic environments than the Class I members. Phylogenomic analysis suggested that the Class I and II methanogens might have evolved before and around the Great Oxygenation Event, respectively. The enrichment of antioxidant features in the Class II methanogens may have played a key role in the adaption of this group to oxidative environments today and historically.

Biomolecular NMR: Past and future

The editors of this special volume suggested this topic, presumably because of the perspective lent by our combined >90-year association with biomolecular NMR. What follows is our personal experience with the evolution of the field, which we hope will illustrate the trajectory of change over the years. As for the future, one can confidently predict that it will involve unexpected advances. Our narrative is colored by our experience in using the NMR Facility for Biomedical Studies at Carnegie-Mellon University (Pittsburgh) and in developing similar facilities at Purdue (1977-1984) and the University of Wisconsin-Madison (1984-). We have enjoyed developing NMR technology and making it available to collaborators and users of these facilities. Our group's association with the Biological Magnetic Resonance data Bank (BMRB) and with the Worldwide Protein Data Bank (wwPDB) has also been rewarding. Of course, many groups contributed to the early growth and development of biomolecular NMR, and our brief personal account certainly omits many important milestones.

Experimental and DFT Investigations Reveal the Influence of the Outer Coordination Sphere on the Vibrational Spectra of Nickel-Substituted Rubredoxin, a Model Hydrogenase Enzyme

Nickel-substituted rubredoxin (NiRd) is a functional enzyme mimic of hydrogenase, highly active for electrocatalytic and solution-phase hydrogen generation. Spectroscopic methods can provide valuable insight into the catalytic mechanism, provided the appropriate technique is used. In this study, we have employed multiwavelength resonance Raman spectroscopy coupled with DFT calculations on an extended active-site model of NiRd to probe the electronic and geometric structures of the resting state of this system. Excellent agreement between experiment and theory is observed, allowing normal mode assignments to be made on the basis of frequency and intensity analyses. Both metal-ligand and ligand-centered vibrational modes are enhanced in the resonance Raman spectra. The latter provide information about the hydrogen bonding network and structural distortions due to perturbations in the secondary coordination sphere. To reproduce the resonance enhancement patterns seen for high-frequency vibrational modes, the secondary coordination sphere must be included in the computational model. The structure and reduction potential of the Ni^IIIRd state have also been investigated both experimentally and computationally. This work begins to establish a foundation for computational resonance Raman spectroscopy to serve in a predictive fashion for investigating catalytic intermediates of NiRd.

Synthesis and Speciation-Dependent Properties of a Multimetallic Model Complex of NiSOD That Exhibits Unique Hydrogen-Bonding

The complex Na₃[{Ni^II(nmp)}₃S₃BTA^alk)] (1) (nmp^2- = deprotonated form of N-(2-mercaptoethyl)picolinamide; H₃S₃BTA^alk = N¹,N³,N⁵-tris(2-mercaptoethyl)benzene-1,3,5-tricarboxamide, where H = dissociable protons), supported by the thiolate-benzenetricarboxamide scaffold (S₃BTA^alk), has been synthesized as a trimetallic model of nickel-containing superoxide dismutase (NiSOD). X-ray absorption spectroscopy (XAS) and ¹H NMR measurements on 1 indicate that the Ni^II centers are square-planar with N₂S₂ coordination, and Ni-N and Ni-S distances of 1.95 and 2.16 Å, respectively. Additional evidence from IR indicates the presence of H-bonds in 1 from the approximately -200 cm^-1 shift in ν_NH from free ligand. The presence of H-bonds allows for speciation that is temperature-, concentration-, and solvent-dependent. In unbuffered water and at low temperature, a dimeric complex (1^A; λ = 410 nm) that aggregates through intermolecular NH···O═C bonds of BTA units is observed. Dissolution of 1 in pH 7.4 buffer or in unbuffered water at temperatures above 50 °C results in monomeric complex (1^M; λ = 367 nm) linked through intramolecular NH···S bonds. DFT computations indicate a low energy barrier between 1^A and 1^M with nearly identical frontier MOs and Ni-ligand metrics. Notably, 1^A and 1^M exhibit remarkable stability in protic solvents such as MeOH and H₂O, in stark contrast to monometallic [Ni^II(nmp)(SR)]^- complexes. The reactivity of 1 with excess O₂, H₂O₂, and O₂^•- is species-dependent. IR and UV-vis reveal that 1^A in MeOH reacts with excess O₂ to yield an S-bound sulfinate, but does not react with O₂^•-. In contrast, 1^M is stable to O₂ in pH 7.4 buffer, but reacts with O₂^•- to yield a putative [Ni^II(nmp)(O₂)]^- complex from release of the BTA-thiolate based on EPR.

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.

Enhanced cancer cell killing of a Pt(IV) prodrug promoted by outer-sphere coordination with polyethyleneimines

The cisplatin prodrug c,c,t-[Pt(NH3)2Cl2(O2CCH2CH2CO2)2](2-) (1) forms outer-sphere coordination interactions with non-toxic low MW PEI, which results in enhanced cancer cell killing, also achieved using PEI-coated AuNPs.

Secondary Sphere Hydrogen Bonding in Monocopper Complexes of Potentially Dinucleating Bis(carboxamide) Ligands

Reaction of a macrocyclic ligand precursor comprising two bis(carboxamido)pyridine units (H4L4) connected by ethylene linkers with NMe4OH and CuX2 (X = Cl, OAc, or OTf) yielded monocopper complexes [NMe4][(H2L4)Cu(X)] (X = Cl (3), OAc (4), or OH (5)), in contrast to previous work using a related ligand with ortho-phenylene linkers wherein dicopper compounds were isolated. X-ray structures of the complexes revealed hydrogen bonding from the free carboxamide N-H groups in the doubly protonated form of the ligand (H2L4- ) to the monodentate fourth ligand coordinated to the Cu(II) ion. Similar secondary sphere hydrogen bonding interactions were identified in multinuclear compounds [NMe4]2[((H2L4)Cu)n(CO3)] (n = 2 or 3) that were isolated from exposure of 5 to air. Cyclic voltammetry revealed oxidations of 3 and 5 at potentials ~300 mV higher than analogous monocopper complexes of bis(arylcarboxamido)pyridine ligands which lack the intramolecular hydrogen bonds, consistent with removal of electron density from the metal center by the hydrogen bonding array. Another ligand variant (H4L5) with ortho-phenylene linkers and only one bis(carboxamido)pyridine moiety yielded monocopper complexes [NMe4][(H2L5)Cu(OAc)] • DMF (8) and [NMe4][(H2L5)CuCl)] • CH3CN (9), but the X-ray structures revealed a different hydrogen bonding arrangement to the solvate molecules. Nonetheless, a high redox potential for 9 was observed, consistent with intramolecular hydrogen bonding interactions in solution.

Determination of Hydrogen Bond Structure in Water versus Aprotic Environments To Test the Relationship Between Length and Stability

Hydrogen bonds profoundly influence the architecture and activity of biological macromolecules. Deep appreciation of hydrogen bond contributions to biomolecular function thus requires a detailed understanding of hydrogen bond structure and energetics and the relationship between these properties. Hydrogen bond formation energies (ΔGf) are enormously more favorable in aprotic solvents than in water, and two classes of contributing factors have been proposed to explain this energetic difference, focusing respectively on the isolated and hydrogen-bonded species: (I) water stabilizes the dissociated donor and acceptor groups much better than aprotic solvents, thereby reducing the driving force for hydrogen bond formation; and (II) water lengthens hydrogen bonds compared to aprotic environments, thereby decreasing the potential energy within the hydrogen bond. Each model has been proposed to provide a dominant contribution to ΔGf, but incisive tests that distinguish the importance of these contributions are lacking. Here we directly test the structural basis of model II. Neutron crystallography, NMR spectroscopy, and quantum mechanical calculations demonstrate that O-H···O hydrogen bonds in crystals, chloroform, acetone, and water have nearly identical lengths and very similar potential energy surfaces despite ΔGf differences >8 kcal/mol across these solvents. These results rule out a substantial contribution from solvent-dependent differences in hydrogen bond structure and potential energy after association (model II) and thus support the conclusion that differences in hydrogen bond ΔGf are predominantly determined by solvent interactions with the dissociated groups (model I). These findings advance our understanding of universal hydrogen-bonding interactions and have important implications for biology and engineering.

Hidden relationships between metalloproteins unveiled by structural comparison of their metal sites

Metalloproteins account for a substantial fraction of all proteins. They incorporate metal atoms, which are required for their structure and/or function. Here we describe a new computational protocol to systematically compare and classify metal-binding sites on the basis of their structural similarity. These sites are extracted from the MetalPDB database of minimal functional sites (MFSs) in metal-binding biological macromolecules. Structural similarity is measured by the scoring function of the available MetalS2 program. Hierarchical clustering was used to organize MFSs into clusters, for each of which a representative MFS was identified. The comparison of all representative MFSs provided a thorough structure-based classification of the sites analyzed. As examples, the application of the proposed computational protocol to all heme-binding proteins and zinc-binding proteins of known structure highlighted the existence of structural subtypes, validated known evolutionary links and shed new light on the occurrence of similar sites in systems at different evolutionary distances. The present approach thus makes available an innovative viewpoint on metalloproteins, where the functionally crucial metal sites effectively lead the discovery of structural and functional relationships in a largely protein-independent manner.

Importance of polarization effect in the study of metalloproteins: application of polarized protein specific charge scheme in predicting the reduction potential of azurin

Molecular dynamics (MD) simulation is commonly used in the study of protein dynamics, and in recent years, the extension of MD simulation to the study of metalloproteins is gaining much interest. Choice of force field is crucial in MD studies, and the inclusion of metal centers complicates the process of accurately describing the electrostatic environment that surrounds the redox centre. Herein, we would like to explore the importance of including electrostatic contribution from both protein and solvent in the study of metalloproteins. MD simulations with the implementation of thermodynamic integration will be conducted to model the reduction process of azurin from Pseudomonas aeruginosa. Three charge schemes will be used to derive the partial charges of azurin. These charge schemes differ in terms of the amount of immediate environment, respective to copper, considered during charge fitting, which ranges from the inclusion of copper and residues in the first coordination sphere during density functional theory charge fitting to the comprehensive inclusion of protein and solvent effect surrounding the metal centre using polarized protein-specific charge scheme. From the simulations conducted, the relative reduction potential of the mutated azurins respective to that of wild-type azurin (ΔEcal) were calculated and compared with experimental values. The ΔEcal approached experimental value with increasing consideration of environmental effect hence substantiating the importance of polarization effect in the study of metalloproteins. This study also attests the practicality of polarized protein-specific charge as a computational tool capable of incorporating both protein environment and solvent effect into MD simulations.

Investigating the function of [2Fe-2S] cluster N1a, the off-pathway cluster in complex I, by manipulating its reduction potential

NADH:quinone oxidoreductase (complex I) couples NADH oxidation and quinone reduction to proton translocation across an energy-transducing membrane. All complexes I contain a flavin to oxidize NADH, seven iron–sulfur clusters to transfer electrons from the flavin to quinone and an eighth cluster (N1a) on the opposite side of the flavin. The role of cluster N1a is unknown, but Escherichia coli complex I has an unusually high-potential cluster N1a and its reduced flavin produces H₂O₂, not superoxide, suggesting that cluster N1a may affect reactive oxygen species production. In the present study, we combine protein film voltammetry with mutagenesis in overproduced N1a-binding subunits to identify two residues that switch N1a between its high- (E. coli, valine and asparagine) and low- (Bos taurus and Yarrowia lipolytica, proline and methionine) potential forms. The mutations were incorporated into E. coli complex I: cluster N1a could no longer be reduced by NADH, but H₂O₂ and superoxide production were unaffected. The reverse mutations (that increase the potential by ~0.16 V) were incorporated into Y. lipolytica complex I, but N1a was still not reduced by NADH. We conclude that cluster N1a does not affect reactive oxygen species production by the complex I flavin; it is probably required for enzyme assembly or stability.

Two residues that determine the potential of cluster N1a in respiratory complex I were identified, and their effects on its flavin-site reactions were determined. Reduction of cluster N1a by NADH does not affect reactive oxygen species production by the flavin.

Metalloproteins diversified: the auracyanins are a family of cupredoxins that stretch the spectral and redox limits of blue copper proteins

The metal sites of electron transfer proteins are tuned for function. The type 1 copper site is one of the most utilized metal sites in electron transfer reactions. This site can be tuned by the protein environment from +80 mV to +680 mV in typical type 1 sites. Accompanying this huge variation in midpoint potentials are large changes in electronic structure, resulting in proteins that are blue, green, or even red. Here, we report a family of blue copper proteins, the auracyanins, from the filamentous anoxygenic phototroph Chloroflexus aurantiacus that display the entire known spectral and redox variations known in the type 1 copper site. C. aurantiacus encodes four auracyanins, labeled A-D. The midpoint potentials vary from +83 mV (auracyanin D) to +423 mV (auracyanin C). The electronic structures vary from classical blue copper UV-vis absorption spectra (auracyanin B) to highly perturbed spectra (auracyanins C and D). The spectrum of auracyanin C is temperature-dependent. The expansion and divergent nature of the auracyanins is a previously unseen phenomenon.

Parameters for molecular dynamics simulations of iron-sulfur proteins

Iron-sulfur proteins involved in electron transfer reactions have finely tuned redox potentials, which allow them to be highly efficient and specific. Factors such as metal center solvent exposure, interaction with charged residues, or hydrogen bonds between the ligand residues and amide backbone groups have all been pointed out to cause such specific redox potentials. Here, we derived parameters compatible with the AMBER force field for the metal centers of iron-sulfur proteins and applied them in the molecular dynamics simulations of three iron-sulfur proteins. We used density-functional theory (DFT) calculations and Seminario's method for the parameterization. Parameter validation was obtained by matching structures and normal frequencies at the quantum mechanics and molecular mechanics levels of theory. Having guaranteed a correct representation of the protein coordination spheres, the amide H-bonds and the water exposure to the ligands were analyzed. Our results for the pattern of interactions with the metal centers are consistent to those obtained by nuclear magnetic resonance spectroscopy (NMR) experiments and DFT calculations, allowing the application of molecular dynamics to the study of those proteins.

Measuring protein reduction potentials using 15N HSQC NMR spectroscopy

NMR spectroscopy was used to measure reduction potentials of four redox proteins by following multiple (15)N HSQC protein resonances across a titration series using mixtures of oxidised and reduced glutathione. Results for PDI a, PDI ab and DsbA agree with the literature and our result for ERp18 confirms this protein as an oxidoreductase of comparable or greater reducing strength than PDI a.

Dissection of hydrogen bond interaction network around an iron-sulfur cluster by site-specific isotope labeling of hyperthermophilic archaeal Rieske-type ferredoxin

The electronic structure and geometry of redox-active metal cofactors in proteins are tuned by the pattern of hydrogen bonding with the backbone peptide matrix. In this study we developed a method for selective amino acid labeling of a hyperthermophilic archaeal metalloprotein with engineered Escherichia coli auxotroph strains, and we applied this to resolve the hydrogen bond interactions with the reduced Rieske-type [2Fe-2S] cluster by two-dimensional pulsed electron spin resonance technique. Because deep electron spin-echo envelope modulation of two histidine (14)N(δ) ligands of the cluster decreased non-coordinating (15)N signal intensities via the cross-suppression effect, an inverse labeling strategy was employed in which (14)N amino acid-labeled archaeal Rieske-type ferredoxin samples were examined in an (15)N-protein background. This has directly identified Lys45 N(α) as providing the major pathway for the transfer of unpaired electron spin density from the reduced cluster by a "through-bond" mechanism. All other backbone peptide nitrogens interact more weakly with the reduced cluster. The extension of this approach will allow visualizing the three-dimensional landscape of preferred pathways for the transfer of unpaired spin density from a paramagnetic metal center onto the protein frame, and will discriminate specific interactions by a "through-bond" mechanism from interactions which are "through-space" in various metalloproteins.

Secondary interactions involving zinc-bound ligands: roles in structural stabilization and macromolecular interactions

A large number of proteins contain bound zinc ions. These zinc ions are frequently coordinated by a combination of histidine and cysteine residues. In addition to atoms that coordinate directly to the zinc ions, these side chains have groups that can donate or accept hydrogen bonds from other groups. These secondary interactions can help stabilize the zinc-binding sites, can contribute to protein folding and stability, and, on occasion, can participate in interactions with other macromolecules. Five examples of these secondary interactions are discussed: carbonic anhydrase (where secondary interactions involving histidine residues stabilize the zinc-binding site thermodynamically and kinetically), retroviral nucleocapsid proteins and TRAF proteins (where cysteinate sulfur to peptide NH hydrogen bonds contribute to the structural relationships between adjacent domains), and nucleic acid binding proteins, Zif268 and TIS11 where secondary interactions participate in protein-nucleic acid interactions.

Hyperfine-shifted 13C resonance assignments in an iron-sulfur protein with quantum chemical verification: aliphatic C-H···S 3-center-4-electron interactions

Although the majority of noncovalent interactions associated with hydrogen and heavy atoms in proteins and other biomolecules are classical hydrogen bonds between polar N−H or O−H moieties and O atoms or aromatic π electrons, high-resolution X-ray crystallographic models deposited in the Protein Data Bank show evidence for weaker C−H···O hydrogen bonds, including ones involving sp³-hybridized carbon atoms. Little evidence is available in proteins for the (even) weaker C−H···S interactions described in the crystallographic literature on small molecules. Here, we report experimental evidence and theoretical verification for the existence of nine aliphatic (sp³-hybridized) C−H···S 3-center−4-electron interactions in the protein Clostridium pasteurianum rubredoxin. Our evidence comes from the analysis of carbon-13 NMR chemical shifts assigned to atoms near the iron at the active site of this protein. We detected anomalous chemical shifts for these carbon-13 nuclei and explained their origin in terms of unpaired spin density from the iron atom being delocalized through interactions of the type: C−H···S−Fe, where S is the sulfur of one of the four cysteine side chains covalently bonded to the iron. These results suggest that polarized sulfur atoms in proteins can engage in multiple weak interactions with surrounding aliphatic groups. We analyze the strength and angular dependence of these interactions and conclude that they may contribute small, but significant, stabilization to the molecule.