FeNO structure in distal pocket mutants of myoglobin based on resonance Raman spectroscopy

FeNO vibrational frequencies were investigated for a series of myoglobin mutants using isotope-edited resonance Raman spectra of (15/14)NO adducts, which reveal the FeNO and NO stretching modes. The latter give rise to doublet bands, as a result of Fermi resonances with coincident porphyrin vibrations; these doublets were analyzed by curve-fitting to obtain the nuNO frequencies. Variations in nuNO among the mutants correlate with the reported nuCO variations for the CO adducts of the same mutants. The correlation has a slope near unity, indicating equal sensitivity of the NO and CO bonds to polar influences in the heme pocket. A few mutants deviate from the correlation, indicating that distal interactions differ for the NO and CO adducts, probably because of the differing distal residue geometries. In contrast to the strong and consistent nuFeC/nuCO correlation found for the CO adducts, nuFeN correlates only weakly with nuNO, and the slope of the correlation depends on which residue is being mutated. This variability is suggested to arise from steric interactions, which change the FeNO angle and therefore alter the Fe-NO and N-O bond orders. This effect is modeled with Density Functional Theory (DFT) and is rationalized on the basis of a valence isomer bonding model. The FeNO unit, which is naturally bent, is a more sensitive reporter of steric interactions than the FeCO unit, which is naturally linear. An important additional factor is the strength of the bond to the proximal ligand, which modulates the valence isomer equilibrium. The FeNO unit is bent more strongly in MbNO than in protein-free heme-NO complexes because of a combination of a strengthened proximal bond and distal interactions.     

DFT analysis of axial and equatorial effects on heme-CO vibrational modes: applications to CooA and H-NOX heme sensor proteins

Determinants of the Fe-CO and C-O stretching frequencies in (imidazole)heme-CO adducts have been investigated via density functional theory (DFT) analysis, in connection with puzzling characteristics of the heme sensor protein CooA and of the H-NOX (Heme-Nitric Oxide and/or OXygen binding) family of proteins, including soluble guanylate cyclase (sGC). The computations show that two mechanisms of Fe-histidine bond weakening have opposite effects on the nuFeC/nuCO pattern. Mechanical tension is expected to raise nuFeC with little change in nuCO whereas the weakening of H-bond donation from the imidazole ligand has the opposite effect. Data on CooA indicate imidazole H-bond weakening associated with heme displacement, as part of the activation mechanism. The computations also reveal that protein-induced distortion of the porphyrin ring, a prominent structural feature of the H-NOX protein TtTar4H (Thermoanaerobacter tengcongensis Tar4 protein heme domain), has surprisingly little effect on nuFeC or nuCO. However, another structural feature, strong H-bonding to the propionates, is suggested to account for the weakened back bonding that is evident in sGC. TtTar4H-CO itself has an elevated nuFeC, which is successfully modeled as a compression effect, resulting from steric crowding in the distal pocket. nuFeC/nuCO data, in conjunction with modeling, can provide valuable insight into mechanisms for heme-protein modulation.     

High pressure fourier transform infrared (FT-IR) spectroscopic studies on inducible nitric oxide (NO) synthase active site: a comparison to cytochrome p450CAM.

High pressure Fourier transform infrared (FT-IR) spectroscopy is performed for the first time to analyse the active site of inducible nitric oxide synthase (iNOSox) using the carbon monoxide (CO) heme iron ligand stretch mode (nuCO) as spectroscopic probe. A membrane-driven sapphire anvil high-pressure cell is used. Three major conformational substates exist in substrate-free iNOSox which are characterized by nuCO at approximately 1936, 1945 and 1952 cm(-1). High pressure favors the 1936 cm(-1) substate with a volume difference to the 1945 substate of approximately -21 cm3/mol. The pressure induced cytochrome P420 formation with a reaction volume of approximately -80 cm3/mol is observed. Arginine binding produces a very low nuCO at approximately 1905 cm(-1) caused by the H-bond from the substrate to CO. nuCO for the substates in the substrate-free and arginine-bound proteins shift linearly with pressure which is qualitatively similar to the observation on cytochrome P450cam. The slightly smaller positive slope of the shift in substrate-free iNOSox compared to substrate-free P450cam is interpreted as a slightly lesser compressible heme pocket. In contrast, the significant slower negative slope for arginine-bound iNOSox compared to camphor-bound P450cam results from the different kind of interactions to the CO ligand (electrostatic interaction in P450cam, H-bond in iNOSox).     

Thr-E11 regulates O2 affinity in Cerebratulus lacteus mini-hemoglobin

The mini-hemoglobin from Cerebratulus lacteus (CerHb) belongs to a class of globins containing the polar Tyr-B10/Gln-E7 amino acid pair that normally causes low rates of O2 dissociation and ultra-high O2 affinity, which suggest O2 sensing or NO scavenging functions. CerHb, however, has high rates of O2 dissociation (kO2 = 200-600 s(-1)) and moderate O2 affinity (KO2) approximately 1 microm(-1)) as a result of a third polar amino acid in its active site, Thr-E11. When Thr-E11 is replaced by Val, kO2 decreases 1000-fold and KO2 increases 130-fold at pH 7.0, 20 degrees C. The mutation also shifts the stretching frequencies of both heme-bound and photodissociated CO, indicating marked changes of the electrostatic field at the active site. The crystal structure of Thr-E11 --> Val CerHbO2 at 1.70 A resolution is almost identical to that of the wild-type protein (root mean square deviation of 0.12 A). The dramatic functional and spectral effects of the Thr-E11 --> Val mutation are due exclusively to changes in the hydrogen bonding network in the active site. Replacing Thr-E11 with Val "frees" the Tyr-B10 hydroxyl group to rotate toward and donate a strong hydrogen bond to the heme-bound ligand, causing a selective increase in O2 affinity, a decrease of the rate coefficient for O2 dissociation, a 40 cm(-1) decrease in nuCO of heme-bound CO, and an increase in ligand migration toward more remote intermediate sites.     

Novel heme ligand displacement by CO in the soluble hemophore HasA and its proximal ligand mutants: implications for heme uptake and release

HasASM, a hemophore secreted by the Gram-negative bacteria Serratia marcescens, extracts heme from host hemoproteins and shuttles it to HasRSM, a specific hemophore outer membrane receptor. Heme iron in HasASM is in a six-coordinate ferric state. It is linked to the protein by the heretofore uncommon axial ligand set, His32 and Tyr75. A third residue of the heme pocket, His83, plays a crucial role in heme ligation through hydrogen bonding to Tyr75. The vibrational frequencies of coordinated carbon monoxide constitute a sensitive probe of trans ligand field, FeCO structure, and electrostatic landscape of the distal heme pockets of heme proteins. In this study, carbonyl complexes of wild-type (WT) HasASM and its heme pocket mutants His32Ala, Tyr75Ala, and His83Ala were characterized by resonance Raman spectroscopy. The CO complexes of WT HasASM, HasASM(His32Ala), and HasASM(His83Ala) exhibit similar spectral features and fall above the line that correlates nuFe-CO and nuC-O for proteins having a proximal imidazole ligand. This suggests that the proximal ligand field in these CO adducts is weaker than that for heme-CO proteins bearing a histidine axial ligand. In contrast, the CO complex of HasASM(Tyr75Ala) has resonance Raman signatures consistent with ImH-Fe-CO ligation. These results reveal that in WT HasASM, the axial ImH side chain of His32 is displaced by CO. This is in contrast to other heme proteins known to have the His/Tyr axial ligand set, wherein the phenolic side chain of the Tyr ligand dissociates upon CO addition. The displacement of His32 and its stabilization in an unbound state is postulated to be relevant to heme uptake and/or release.     

A photolysis-triggered heme ligand switch in H93G myoglobin

Resonance Raman spectroscopy and step-scan Fourier transform infrared (FTIR) spectroscopy have been used to identify the ligation state of ferrous heme iron for the H93G proximal cavity mutant of myoglobin in the absence of exogenous ligand on the proximal side. Preparation of the H93G mutant of myoglobin has been previously reported for a variety of axial ligands to the heme iron (e.g., substituted pyridines and imidazoles) [DePillis, G., Decatur, S. M., Barrick, D., and Boxer, S. G. (1994) J. Am. Chem. Soc. 116, 6981-6982]. The present study examines the ligation states of heme in preparations of the H93G myoglobin with no exogenous ligand. In the deoxy form of H93G, resonance Raman spectroscopic evidence shows water to be the axial (fifth) ligand to the deoxy heme iron. Analysis of the infrared C-O and Raman Fe-C stretching frequencies for the CO adduct indicates that it is six-coordinate with a histidine trans ligand. Following photolysis of CO, a time-dependent change in ligation is evident in both step-scan FTIR and saturation resonance Raman spectra, leading to the conclusion that a conformationally driven ligand switch exists in the H93G protein. In the absence of exogenous nitrogenous ligands, the CO trans effect stabilizes endogenous histidine ligation, while conformational strain favors the dissociation of histidine following photolysis of CO. The replacement of histidine by water in the five-coordinate complex is estimated to occur in < 5 micros. The results demonstrate that the H93G myoglobin cavity mutant has potential utility as a model system for studying the conformational energetics of ligand switching in heme proteins such as those observed in nitrite reductase, guanylyl cyclase, and possibly cytochrome c oxidase.     

Contact order dependent protein folding rates: kinetic consequences of a cooperative interplay between favorable nonlocal interactions and local conformational preferences

Physical mechanisms underlying the empirical correlation between relative contact order (CO) and folding rate among naturally occurring small single-domain proteins are investigated by evaluating postulated interaction schemes for a set of three-dimensional 27mer lattice protein models with 97 different CO values. Many-body interactions are constructed such that contact energies become more favorable when short chain segments sequentially adjacent to the contacting residues adopt native-like conformations. At a given interaction strength, this scheme leads to folding rates that are logarithmically well correlated with CO (correlation coefficient r = 0.914) and span more than 2.5 orders of magnitude, whereas folding rates of the corresponding Gō models with additive contact energies have much less logarithmic correlation with CO and span only approximately one order of magnitude. The present protein chain models also exhibit calorimetric cooperativity and linear chevron plots similar to that observed experimentally for proteins with apparent simple two-state folding/unfolding kinetics. Thus, our findings suggest that CO-dependent folding rates of real proteins may arise partly from a significant positive coupling between nonlocal contact favorabilities and local conformational preferences.     

Cytochrome c peroxidase mutant active site structures probed by resonance Raman and infrared signatures of the CO adducts

Vibrational frequencies associated with FeC and CO stretching and FeCO bending modes have been determined via resonance Raman (RR) and infrared (IR) spectroscopy for cytochrome c peroxidase (CCP) mutants prepared by site-directed mutagenesis. These include the bacterial "wild type", CCP(MI), and mutations involving groups on the proximal (Asp-235----Asn; Trp-191---Phe) and distal (Trp-51----Phe; Arg-48----Leu and Lys) side of the heme. The data were analyzed with the aid of a recently established correlation between nu FeC and nu CO, which can be used to distinguish between back-bonding and axial ligand donor effects. At high pH all adducts showed essentially the same vibrational pattern (form I') with nu FeC approximately 505 cm-1, nu CO approximately 1948 cm-1, and delta FeCO (weak RR band) approximately 576 cm-1. These frequencies are very similar to those shown by the myoglobin CO adduct and imply a "normal" H-bond of the proximal histidine. At pH 7 (pH 6 for Asn-235 and Leu-48), different forms are seen for different proteins: form I (nu FeC approximately 500 cm-1, nu CO = 1922-1941 cm-1, and delta FeCO approximately 580 cm-1, very weak) in the case of CCP(MI) and Phe-191, as well as bakers' yeast CCP, or form II (nu FeC approximately 530 cm-1, nu CO = 1922-1933 cm-1, and delta FeCO = 585 cm-1, moderately strong) for Asn-235 and Phe-51.(ABSTRACT TRUNCATED AT 250 WORDS)     

Computational design of protein-small molecule interfaces

The computational design of proteins that bind small molecule ligands is one of the unsolved challenges in protein engineering. It is complicated by the relatively small size of the ligand which limits the number of intermolecular interactions. Furthermore, near-perfect geometries between interacting partners are required to achieve high binding affinities. For apolar, rigid small molecules the interactions are dominated by short-range van der Waals forces. As the number of polar groups in the ligand increases, hydrogen bonds, salt bridges, cation–π, and π–π interactions gain importance. These partial covalent interactions are longer ranged, and additionally, their strength depends on the environment (e.g. solvent exposure). To assess the current state of protein-small molecule interface design, we benchmark the popular computer algorithm Rosetta on a diverse set of 43 protein–ligand complexes. On average, we achieve sequence recoveries in the binding site of 59% when the ligand is allowed limited reorientation, and 48% when the ligand is allowed full reorientation. When simulating the redesign of a protein binding site, sequence recovery among residues that contribute most to binding was 52% when slight ligand reorientation was allowed, and 27% when full ligand reorientation was allowed. As expected, sequence recovery correlates with ligand displacement.     

DFT study on the reactivity of iron porphyrins tuned by ring substitution

The effect of beta-substituents (-NO2, -Br, -OCH3) in the reactivity of Fe(II) and Fe(III) porphyrins is studied by means of density functional theory (DFT) calculations. The binding of nitric oxide, carbon monoxide and dioxygen (NO, CO, O2) was explored due to the relevance of their interactions in the chemistry of heme proteins and in biomimetic catalysis. The binding capability (BC) of the porphyrins was found to be strongly modulated both by the donor and attractor substituents used in the work. Unexpectedly, we found that the BC of Fe(II) porphyrins is mainly decreased for the diatomic ligands, when both donor or withdrawing substituents were considered. This effect was particularly significant when the ligand was oxygen. The correlation of Fe-X and X-O (X=N, C, O) bond distances is explained in terms of backdonation effects.     

Systematic and quantitative analysis of protein-protein recognition between nonribosomal peptide synthetases investigated in the tyrocidine biosynthetic template

We present a systematic and quantitative study of the protein-protein recognition between the three tyrocidine synthetases TycA, TycB, and TycC investigated with two artificial in trans assay systems, which had been previously developed: the "DKP assay system" for the interaction of TycA with TycB and the "L/D-Phe-L-Asn assay system" for the interaction of TycB with TycC. TycA-A(Phe)TE and TycB(3)-A(Phe)TE, which are used as donor enzymes, both provide D-Phe-S-Ppant, so that no substrate specificities interfered with the quantification of protein-protein recognition. We tested all donor/acceptor enzyme combinations between the two artificial assay systems for product formation activities as well as two hybrid enzymes, where the E-domains of TycA and TycB(3) had been exchanged against each other. Furthermore, four donor/acceptor protein fusions were constructed on gene level, resulting in dimodular proteins. We were able to show that the E-domains mediate protein-protein recognition in trans. Product formation of the different donor assayed with the two acceptor enzymes TycB(1)-CA(Pro)T and TycC(1)-CA(Asn)T/Te in trans was only obtained if the donor enzyme harbored the cognate E-domain. Interestingly, all in cis fusions (dimodular proteins) were active, giving strong evidence that unnatural protein-protein interactions can be "forced" by fusion of the distinct enzymes. Finally, we were able to detect product formation in the "DKP system" with engineered hybrid proteins where the A-domain of TycA had been exchanged against the isoleucine-activating A-domain of BacA(1) and the valine-activating A-domain of TycC(4), respectively. All of these findings are of high relevance for future NRPS engineering approaches.     

Structural dynamics of myoglobin: ligand migration and binding in valine 68 mutants

We have combined Fourier transform infrared/temperature derivative (FTIR-TDS) spectroscopy at cryogenic temperatures and flash photolysis at ambient temperature to examine the effects of polar and bulky amino acid replacements of the highly conserved distal valine 68 in sperm whale myoglobin. In FTIR-TDS experiments, the CO ligand can serve as an internal voltmeter that monitors the local electrostatic field not only at the active site but also at intermediate ligand docking sites. Mutations of residue 68 alter size, shape, and electric field of the distal pocket, especially in the vicinity of the primary docking site (state B). As a consequence, the infrared bands associated with the ligand at site B are shifted. The effect is most pronounced in mutants with large aromatic side chains. Polar side chains (threonine or serine) have only little effect on the peak frequencies. Ligands that migrate toward more remote sites C and D give rise to IR bands with altered frequencies. TDS experiments separate the photoproducts according to their recombination temperatures. The rates and extent of ligand migration among internal cavities at cryogenic temperatures can be used to interpret geminate and bimolecular O2 and CO recombination at room temperature. The kinetics of geminate recombination can be explained by steric arguments alone, whereas both the polarity and size of the position 68 side chain play major roles in regulating bimolecular ligand binding from the solvent.     

H-bonding networks of the distal residues and water molecules in the active site of Thermobifida fusca hemoglobin

The ferric form of truncated hemoglobin II from Thermobifida fusca (Tf-trHb) and its triple mutant WG8F-YB10F-YCD1F at neutral and alkaline pH, and in the presence of CN⁻ have been characterized by resonance Raman spectroscopy, electron paramagnetic resonance spectroscopy, and molecular dynamics simulations. Tf-trHb contains three polar residues in the distal site, namely TrpG8, TyrCD1 and TyrB10. Whereas TrpG8 can act as a potential hydrogen-bond donor, the tyrosines can act as donors or acceptors. Ligand binding in heme-containing proteins is determined by a number of factors, including the nature and conformation of the distal residues and their capability to stabilize the heme-bound ligand via hydrogen-bonding and electrostatic interactions. Since both the RR Fe–OH⁻ and Fe–CN⁻ frequencies are very sensitive to the distal environment, detailed information on structural variations has been obtained. The hydroxyl ligand binds only the WT protein giving rise to two different conformers. In form 1 the anion is stabilized by H-bonds with TrpG8, TyrCD1 and a water molecule, in turn H-bonded to TyrB10. In form 2, H-bonding with TyrCD1 is mediated by a water molecule. Unlike the OH⁻ ligand, CN⁻ binds both WT and the triple mutant giving rise to two forms with similar spectroscopic characteristics. The overall results clearly indicate that H-bonding interactions both with distal residues and water molecules are important structural determinants in the active site of Tf-trHb. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.     

A survey of amino acid side-chain interactions in 21 proteins

Based on the atomic co-ordinate data for 21 representative proteins, the frequencies of long-range interactions between side-chain groups and 15 different types of side-chain atoms have been determined. The observed frequencies are compared to the results expected for random association in order to define a scale of relative affinities. Thirty-five residue-atom pairs exhibit frequencies of interaction that differ by at least 50% from the expected values. The amino acids tend to fall into three classes: non-polar, neutral and polar amino acids. The data are regrouped in a different way to determine the average affinity of each amino acid side-chain group for all other types of side-chain groups. Fourteen side-chain pairs have at least 50% fewer interactions than expected, while 21 side-chain pairs have at least 50% more interactions than expected. Unusual patterns of association are discussed and compared with current ideas about the organization of protein structure.     

Multinuclear (13C, 17O, 57Fe) NMR studies of carbonmonoxy heme proteins and synthetic model compounds

13C, 17O and 57Fe NMR spectra of several carbonmonoxy hemoprotein models with varying polar and steric effects of the distal organic superstructure, constraints of the proximal side, and porphyrin ruffling are reported. Both heme models and heme proteins obey a similar excellent linear delta(13C) versus nu(C-O) relationship which is primarily due to modulation of pi-back-bonding from the Fe d(pi) to CO pi* orbital by the distal pocket polar interactions. The lack of correlation between delta(13C) and delta(17O) suggests that the two probes do not reflect a similar type of electronic and structural perturbation. delta(17O) is not primarily influenced by the local distal field interactions and does not correlate with any single structural property of the Fe-C-O unit; however, atropisomerism and deformation of the porphyrin geometry appear to play a significant role. 57Fe shieldings vary by nearly 900 ppm among various hemes and an excellent correlation was found between delta(57Fe) and the absolute crystallographic average displacement of the meso carbon atoms, /Cm/, relative to the porphyrin core mean plane. The excellent correlation between iron-57 shieldings and the average shieldings of the meso carbons of the porphyrin skeleton of TPP derivatives suggests that the two probes reflect a similar type of electronic and structural perturbation which is primarily porphyrin ruffling.     

A comparative resonance Raman analysis of heme-binding PAS domains: heme iron coordination structures of the BjFixL, AxPDEA1, EcDos, and MtDos proteins

The heme-PAS is a specialized domain with which a broad class of signal-transducing heme proteins detect physiological heme ligands. Such domains exhibit a wide range of ligand binding parameters, yet they are all expected to feature an alpha-beta heme binding fold and a predominantly hydrophobic heme distal pocket without a distal histidine. We have compared, for the first time, the resonance Raman spectra of several heme-PASs: the heme-binding domains of Bradyrhizobium japonicum FixL, Escherichia coli Dos, Acetobacter xylinum PDEA1, and Methanobacterium thermoautotrophicum Dos. In all cases, the nu(Fe)-(CO) and nu(C-O) values of the carbonmonoxy forms were consistent with coordination of the heme iron to histidine on the proximal side and binding of the CO without electrostatic interaction with the heme distal pocket. EcDos was unusual in having predominantly hexacoordinate heme iron in the deoxy and met forms. Despite an evident lack of CO interaction with the EcDos heme pocket, relatively low Fe-O(2) (562 cm(-1)) and N-O (1576 cm(-1)) stretching frequencies indicated that strong polar interactions with that heme distal pocket are possible for highly bent ligands such as O(2) or NO. None of the newly studied NO adducts exhibited evidence of the Fe-His rupture and pentacoordination previously noted for Sinorhizobium meliloti FixL. A low Fe-His stretching frequency, formerly interpreted as a strained Fe-His bond, and the slow association of O(2) with S. meliloti FixL failed to correlate with the newly studied proteins having low association rate or low equilibrium association constants for binding of O(2). We conclude that although heme-PASs share some features, they represent distinct signal transduction mechanisms.     

Unique cyanide nitrogen-15 nuclear magnetic resonance chemical shift values for cyano-peroxidase complexes. Relevance to the heme active-site structure and mechanism of peroxide activation

Cyanide ion has been utilized to probe the heme environment of the ferric states of horseradish peroxidase, lactoperoxidase and chloroperoxidase. The 15N-NMR signal for cyanide bound to these enzymes is located in the downfield region from 578 to 412 ppm (with respect to the nitrate ion reference). The corresponding signal for met-forms of hemoglobin, myoglobin and cytochrome c is much further downfield in the 1047-847 ppm region. The signal position for peroxidases is quite invariant with pH in the physiological ranges. The upfield bias for peroxidase chemical shifts must reflect unique trans iron(III) ligand types and/or proximal-group hydrogen bonding or steric effects. Model compound studies reveal a significant upfield cyanide 15N shift with addition of agents capable of hydrogen-bonding to the coordinated cyanide ion. An even more striking upfield shift of 277 ppm is associated with deprotonation of a trans imidazole residue. The distinctive chemical shifts observed for the cyano ligand in peroxidases support the hypothesis that a distal hydrogen-bonding network and perhaps a polar, basic trans ligand are essential for O-O bond activation by peroxidases.