Challenging Density Functional Theory Calculations with Hemes and Porphyrins

In this paper we review recent advances in computational chemistry and specifically focus on the chemical description of heme proteins and synthetic porphyrins that act as both mimics of natural processes and technological uses. These are challenging biochemical systems involved in electron transfer as well as biocatalysis processes. In recent years computational tools have improved considerably and now can reproduce experimental spectroscopic and reactivity studies within a reasonable error margin (several kcal·mol(-1)). This paper gives recent examples from our groups, where we investigated heme and synthetic metal-porphyrin systems. The four case studies highlight how computational modelling can correctly reproduce experimental product distributions, predicted reactivity trends and guide interpretation of electronic structures of complex systems. The case studies focus on the calculations of a variety of spectroscopic features of porphyrins and show how computational modelling gives important insight that explains the experimental spectra and can lead to the design of porphyrins with tuned properties.

Destructive interactions of dirhodium(ii) tetraacetate with β metallothionein rh1a

The four bridging acetates of dirhodium(ii) tetraacetate are removed in a stepwise manner by human-metallothionein 1a β-fragment, a reaction captured by electrospray ionization mass spectrometry.

A Heme-responsive Regulator Controls Synthesis of Staphyloferrin B in Staphylococcus aureus

Staphylococcus aureus possesses a multitude of mechanisms by which it can obtain iron during growth under iron starvation conditions. It expresses an effective heme acquisition system (the iron-regulated surface determinant system), it produces two carboxylate-type siderophores staphyloferrin A and staphyloferrin B (SB), and it expresses transporters for many other siderophores that it does not synthesize. The ferric uptake regulator protein regulates expression of genes encoding all of these systems. Mechanisms of fine-tuning expression of iron-regulated genes, beyond simple iron regulation via ferric uptake regulator, have not been uncovered in this organism. Here, we identify the ninth gene of the sbn operon, sbnI, as encoding a ParB/Spo0J-like protein that is required for expression of genes in the sbn operon from sbnD onward. Expression of sbnD-I is drastically decreased in an sbnI mutant, and the mutant does not synthesize detectable SB during early phases of growth. Thus, SB-mediated iron acquisition is impaired in an sbnI mutant strain. We show that the protein forms dimers and tetramers in solution and binds to DNA within the sbnC coding region. Moreover, we show that SbnI binds heme and that heme-bound SbnI does not bind DNA. Finally, we show that providing exogenous heme to S. aureus growing in an iron-free medium results in delayed synthesis of SB. This is the first study in S. aureus that identifies a DNA-binding regulatory protein that senses heme to control gene expression for siderophore synthesis.

Kinetics of Zinc and Cadmium Exchanges between Metallothionein and Carbonic Anhydrase

The flexible coordination stoichiometry of a relatively high number of metal ions is a property unique to the metallothionein (MT) family of proteins. Mammalian MTs, for example, accommodate up to seven divalent metal ions in tetrahedral coordination geometries, using its complement of 20 cysteine ligands. The lability of the metals from these metalloclusters has been used to support the proposal of MTs acting as metal chaperones, by donating to other metal-binding proteins. The metal exchange kinetics between human MT1A and carbonic anhydrase (CA) were examined using time-dependent electrospray ionization mass spectrometry (ESI-MS). The time dependence of three different reaction conditions were studied: (i) zinc donation from partially metalated zinc-MT to apoCA; (ii) metal exchange between zinc saturated MTs and cadmium saturated CA (Cd-CA); and (iii) metal exchange between partially metalated zinc-MTs and Cd-CA. The results show that zinc donation from Zn-MTs to apo-zinc-dependent enzymes is dependent on the metal loading of the Znn-MT (where n = 1-7) and that this is a direct consequence of the increasing metal affinity for smaller values of n. Partially metalated MTs are also shown to extract cadmium from Cd-CA with significantly faster rates than metal saturated MTs and that even under zinc limiting conditions, mammalian Cd-CA would not coexist with MT. On the basis of these and previously published results, we suggest that protein-protein interactions between MT and CA facilitate metal transfers through favorable electrostatic interactions and hypothesize that the metal could be transferred between the MT and the enzyme active site using nearby metal-binding functionalities along the transfer pathway.

The zinc balance: competitive zinc metalation of carbonic anhydrase and metallothionein 1A

The small, cysteine-rich metallothionein family of proteins is currently considered to play a critical role in the provision of metals to metalloenzymes. However, there is limited information available on the mechanisms of these fundamentally important interactions. We report on the competitive zinc metalation of apocarbonic anhydrase in the presence of apometallothionein 1A using electrospray-ionization mass spectrometry. These experiments revealed the relative affinities of zinc to all species in solution. The carbonic anhydrase is shown to compete efficiently only against Zn5-7MT. The calculated equilibrium zinc binding constants of each of the 7 zinc metallothionein 1A species ranged from a high of (log(KF)) 12.5 to a low of 11.8. The 8 equilibrium constants connecting the 10 active species in competition for the zinc were modeled by fitting the KF values of the 8 competitive bimolecular reactions to the ESI-mass spectral data. These modeled K values are shown to be experimentally connected to the metalation efficiency of the carbonic anhydrase. The series of 7 metallothionein binding affinities for zinc highlight the buffering role of zinc metallothioneins that permit simultaneously zinc storage and zinc sensing. Finally, the significance of the multiple zinc binding affinities of zinc metallothionein is discussed in relation to zinc homeostasis.

Pentacene-fused diporphyrins

In this work, we report the synthesis, spectroscopic characterization, and theoretical analysis of a linearly conjugated pentacene-fused porphyrin dimer and cross-conjugated quinone-fused dinaphtho[2,3]porphyrins. These multichromophoric systems display non-typical UV-visible absorptions of either porphyrins or pentacenes/quinones. UV-visible, emission and magnetic circular dichroism (MCD) spectroscopy suggest strong electronic interactions among the multichromophores in the system. DFT calculations revealed the delocalization of the HOMOs and LUMOs spanning the entire dimer and linker assembly. The pentacene-fused porphyrin dimer is significantly more stable than both the corresponding pentacene and the heptacene derivatives. The availability of these huge π-extended and electronically highly interactive multichromophoric systems promises unprecedented electronic and photophysical properties.

Rational design of a zinc phthalocyanine binding protein

Phthalocyanines have long been used as primary donor molecules in synthetic light-powered devices due to their superior properties when compared to natural light activated molecules such as chlorophylls. Their use in biological contexts, however, has been severely restricted due to their high degree of self-association, and its attendant photoquenching, in aqueous environments. To this end we report the rational redesign of a de novo four helix bundle di-heme binding protein into a heme and Zinc(II) phthalocyanine (ZnPc) dyad in which the ZnPc is electronically and photonically isolated. The redesign required transformation of the homodimeric protein into a single chain four helix bundle and the addition of a negatively charge sulfonate ion to the ZnPc macrocycle. To explore the role of topology on ZnPc binding two constructs were made and the resulting differences in affinity can be explained by steric interference of the newly added connecting loop. Singular binding of ZnPc was verified by absorption, fluorescence, and magnetic circular dichroism spectroscopy. The engineering guidelines determined here, which enable the simple insertion of a monomeric ZnPc binding site into an artificial helical bundle, are a robust starting point for the creation of functional photoactive nanodevices.

Challenging conventional wisdom: single domain metallothioneins

Metallation studies of human metallothioneins support the role of single metal-binding-domains as commonplace with the typical two-domain-cluster structure as exceptional.

Single-domain metallothioneins: evidence of the onset of clustered metal binding domains in Zn-rhMT 1a

Mammalian metallothioneins bind up to seven Zn(2+) ions in two distinct domains: an N-terminal β-domain that binds three Zn(2+) ions and a C-terminal α-domain that binds four Zn(2+) ions. Domain specificity has been invoked in the metalation mechanism with cluster formation and bridging of the 20 Cys residues taking place prior to saturation with seven Zn(2+) ions. We report a novel experiment that examines Zn(2+) metalation by exploiting the expected decrease in K(F) at the onset of clustering using electrospray ionization mass spectrometry (ESI-MS). During the titration with Zn(2+), the ESI-MS data show that several metalated species coexist until the fully saturated proteins are formed. The relative Zn binding affinities of the seven total sites in the α- and β-fragments were determined through direct competition for added Zn(2+). The K(F) values for each Zn(2+) are expected to decrease as a function of the remaining available sites and the onset of clustering. Analysis shows that Zn(2+) binds to β-rhMT with a greater affinity than α-rhMT. The incremental distribution of Zn(2+) between the competing fragments and apo-βα-rhMT (essentially three and four sites competing with seven sites) identifies the exact point at which clustering begins in the full protein. Analysis of the speciation data shows that Zn(5)-MT forms before clustering begins. This means that all 20 Cys residues of apo-βα-rhMT are bound terminally to Zn(2+) as [Zn(Cys)(4)](2-) units before clustering begins; there is no domain preference in this first metalation stage. Preferential binding of Zn(2+) to β- and α-rhMT at the point where βα-rhMT must form clusters is caused by a significant decrease in the affinity of βα-rhMT for further Zn(2+). The single-domain Zn(5)-rhMT, in which there are no exposed cysteine sulfurs, is a key component of the metalation pathway because the lower affinities of the two clustered Zn(2+) ions allow donation to apoenzymes.

Cysteine accessibility during As3+ metalation of the α- and β-domains of recombinant human MT1a

Metallothionein is a ubiquitous metal binding protein that plays an important role in metal ion homeostasis and redox chemistry within cells. Mammalian metallothioneins bind a wide variety of metals including the metalloid As3+ in two domains (β and α) connected by a short linker sequence. Three As3+ bind in each domain for a total of 6 As3+ per protein. In recombinant human metallothionein (rh-MT1a) each As3+ binds three cysteine residues to form As3Cys9(CysSH)2-α-rhMT1a in the 11 Cys α-domain and As3Cys9-β-rhMT1a in the 9 Cys β-domain. This means that there should be 2 free cysteines in the α-domain but no free cysteines in the β-domain. By using benzoquinone, the number and relative accessibility of the free cysteinyl thiols during the metalation reactions were determined. The electrospray ionization mass spectrometry (ESI-MS) data confirmed that each As3+ binds using exactly 3 cysteine thiols and showed that there was a significant difference in the reactivity of the free cysteines during the metalation reaction. After a reaction with two molar equivalents of As3+ to form As2Cys6(CysSH)3-αβ-rhMT1a, the remaining 3 Cys in the 9 Cys β-domain were far less reactive than those in the α-domain. Molecular dynamics calculations for the metalation reactions with As3+ measured by ESI-MS allowed an interpretation of the mass spectral data in terms of the relative location of the cysteine thiols that were not involved in As3+ coordination. Together, these data provide insight into the selection of a specific cysteinyl thiol by the incoming metals during the stepwise metalation of metallothioneins.

Soluble diamagnetic model for malaria pigment: coordination chemistry of gallium(III)protoporphyrin-IX

The facile axial ligand exchange properties of gallium(III) protoporphyrin IX in methanol solution were utilized to explore self-association interactions by NMR techniques. Structural changes were observed, as well as competitive behavior with the ligands acetate and fluoride, which differed from that seen with the synthetic analogue gallium(III) octaethylporphyrin which lacks acid groups in its side-chains and has less solution heterogeneity as indicated by absorption and MCD spectroscopies. The propionic acid side chains of protoporphyrin IX are implicated in all such interactions of PPIX, and both dynamic metal-propionic interactions and the formation of propionate-bridged dimers are observed. Fluoride coordination provides an unusual example of slow ligand exchange, and this allows for the identification of a fluoride bridged dimer in solution. An improved synthesis of the chloride and hydroxide complexes of gallium(III) protoporphyrin IX is reported. An insoluble gallium analogue of hematin anhydride is described. In general, the interactions between solvent and the metal are found to confer very high solubility, making [Ga(PPIX)](+) a useful model for ferric heme species.

Structural properties of metal-free apometallothioneins

The metalated forms of metallothionein are well studied (particularly Zn-MT, Cu-MT and Cd-MT), but almost nothing is known about the chemical and structural properties of apometallothioneins despite their importance in initial metalation and subsequent demetalation. Electrospray ionization mass spectrometry was used to provide a detailed view of the structural properties of the metal-free protein. Mass spectra of Zn(7)-MT and apo-MT at pH 7 exhibit the same charge state distribution, indicating that apo-MT is tightly folded like the metallated protein, whereas apo-MT at pH 3 exhibits a charge state spectrum associated with unfolding or denaturation. Benzoquinone was used to modify the cysteines in the β-MT (9 Bq), and α-MT (11 Bq) fragments, and the full βα-MT (20 Bq) protein. ESI-MS showed that the overall volume and, therefore, the extent of folding for the modified proteins is similar to that of Zn-MT. Molecular modeling using MM3-MD methods provided the volume of each modified protein. The volumes of the partially modified proteins follow the same trend as the charge states, showing that ESI-MS is an excellent method with which to follow small changes in protein folding as a function of applied chemical stress. The data suggest that the structure of apo-βα-MT is more organized than previously considered.

Insight into blocking heme transfer by exploiting molecular interactions in the core Isd heme transporters IsdA-NEAT, IsdC-NEAT, and IsdE of Staphylococcus aureus

The pathogenic bacterium Staphylococcus aureus has adopted specialized mechanisms for scavenging iron from its host. The nine cell wall and membrane-associated iron regulated surface determinant (Isd) proteins (IsdH, IsdB, IsdA, IsdC, IsdDEF, IsdG and IsdI) allow Staphylococcus aureus to scavenge iron from the heme in hemoglobin and haptoglobin-hemoglobin. Of these, it is IsdE that chaperones the heme to the ATP binding cassette-type transmembrane transporter (IsdF). IsdH, IsdB, IsdA and IsdC contain at least one heme binding Near Transporter (NEAT) domain. Previous studies have shown that ferric heme is transferred unidirectionally in the sequence IsdA-NEAT (Tyr - proximal amino acid) → IsdC-NEAT (Tyr) → IsdE (His). IsdA-NEAT does not transfer heme directly to IsdE. In this paper we investigated PPIX transfer through the core cell wall proteins of the Isd system (IsdA-NEAT, IsdC-NEAT and IsdE) with FePPIX-dimethylester, and the metal substituted CoPPIX and MnPPIX using ESI-MS, UV-visible absorption and MCD spectroscopy. IsdA binds each of the rings but the subsequent transfer properties to IsdC-N or IsdE are not the same as found with heme. FePPIX-DME transfers from IsdA-N to IsdC-N but neither protein transfers the ring to IsdE. IsdA-N does not transfer CoPPIX to IsdC-N or IsdE. IsdA-N does transfer MnPPIX to both IsdC-N and IsdE. Significantly, it is possible that since CoPPIX and FePPIX-DME bind to IsdA-N, the lack of transfer to IsdC-N and subsequently to IsdE for CoPPIX could prove to be used as a potential disruption agent to the S. aureus heme transfer system and may identify a possible anti-microbial.

Application of magnetic circular dichroism spectroscopy to porphyrins, phthalocyanines and hemes

This review of the application of the MCD spectroscopic technique to porphyrinoids aims to illustrate the features and unique information that can be obtained and used in the assignment of the optical absorption spectrum and also the electronic structures of porphyrins and phthalocyanines. By illustrating each of the major and readily available spectral features with simple compounds we hope this review will act to guide new users in formulating the questions that room temperature MCD measurements can answer. The review is in three parts: (i) a brief introduction to the instrumental and theoretical background with an emphasis on the parameters that can be readily extracted without computational or theoretical analyses. (ii) Illustration of each MCD spectral feature observed in room temperature spectra. The MCD spectra Zn -porphyrins and phthalocyanines, and chlorophyll a are described to provide models for measurements of new compounds. (iii) The use of "finger printing" to determine the oxidation and spin states of the central iron in heme proteins, and to identify the proximal (5th) and distal (6th) iron binding ligands from the MCD spectral envelope observed between 250 nm and 1000 nm. The value of this approach is immeasurable because through room temperature measurement of the MCD spectra from dilute solutions under conditions that allow the heme protein to be stabilized it is possible to identify accurately the proximal and distal ligands. Identifying the heme binding environment for heme proteins newly purified and for heme enzymes during their reaction cycle allows function to be considered long before X-ray structural analysis is available.

Spectroscopic and theoretical studies of Ga(III)protoporphyrin-IX and its reactions with myoglobin

Ga(III)protoporphyrin-IX (Ga-PP) has been proposed as a model for the key interporphyrin interactions in malaria pigment. Unlike the paramagnetic parent iron heme derivatives, Ga-PP is readily soluble in methanol (MeOH). We report optical, mass spectroscopic, and theoretical results for Ga-PP as well as its reactions with myoglobin. UV-visible absorption and MCD spectroscopy show that Ga-PP exhibits a typical spectrum for a main group metal: a Q-band at 539 nm and a B band at 406 nm when dissolved in MeOH. We also report optical data for Zn(II)protoporphyrin IX (Zn-PP) dissolved in MeOH, which exhibits a Q-band at 545 nm and a B band at 415 nm. ESI mass spectral data for Ga-PP dissolved in MeOH show the presence of predominantly monomers, with smaller fractions of dimers [(Ga-PP)(2)] and trimers. UV-visible and MCD absorption spectroscopy and ESI mass spectral data demonstrate the successful insertion of monomeric Ga-PP into apo-Mb. Ga-PP-Mb exhibits a B band at 417 nm and Q bands at 545 and 584 nm, which are all red-shifted from the free Ga-PP values. The calculated electronic structures and frontier molecular orbitals of Ga-PP, (Ga-PP)(2) and Zn-PP fit the previously reported trends in band energies and oscillator strengths as a function of molecular orbital energies. These new data can be applied to explain the experimentally observed optical spectroscopy. The observed Q-band energies are accounted for by calculated (HOMO-LUMO) gap of the frontier MOs, while the split in the two top occupied MOs accounts for the magnitude of the Q-band oscillator strength as well as the experimentally observed Q to B band energy separation. Although Ga-PP shares more spectroscopic properties with Zn-PP than it does with Fe(III)PPIX, the trivalent oxidation state allows this molecule to be used as a model for ferric hemes in heme proteins.