Equilibrium O2 distribution in the Zn2+-protoporphyrin IX deoxymyoglobin mimic: application to oxygen migration pathway analysis

Proton spin relaxation induced by the triplet ground state of O(2) in the zinc-containing diamagnetic analogue of sperm whale deoxymyoglobin has been measured as a function of oxygen concentration. As no covalent binding of oxygen to the metal occurs in the zinc species, the relaxation effects of O(2) on the protein (1)H resonances arise exclusively via much weaker noncovalent interactions. The relaxation effects at the amide proton sites are found to be highly localized and are derived almost exclusively from O(2) binding at the four previously identified xenon binding sites. Relative binding constants of 1.0, 0.08, 0.07, and 0.23 were determined for the Xe 1, Xe 2, Xe 3, and Xe 4 sites, respectively. In combination with earlier measurements of the kinetics of the heme binding of oxygen, these equilibria measurements enable a more detailed analysis of models characterizing O(2) entry and egress. A correlation is established between the fraction of O(2) which enters the Fe(2+)-binding site via rotation of the distal histidine side chain (so-called "histidine gate") and the experimentally observable O(2) (or CO) lifetime in the Xe 1 site. A physiological role for these secondary oxygen binding sites is proposed in enhancing the efficiency of the O(2) association reaction by rendering more favorable its competition with water binding in the distal heme pocket.     

A comparison of predicted and experimental vibrational spectra in some small fluorocarbons

The vibrational assignments of difluoromethane, cis and trans 1,2-difluoroethylene, trifluoroethylene, 2-fluoropropene, 1,1-difluoroethylene, and tetrafluoroethylene have been reassessed in the light of recent theoretical work. Inconsistent experimental wave number values and assignments have been clarified, and some changes to the experimental assignments are proposed. The assignments compare favorably with recent scaled density functional theory calculations using the hybrid three-parameter B3-PW91 density functional.     

Near-field diffraction in a two-dimensional V-groove and its role in SERS

Optical field distribution in micro-nano geometries of miniaturized optical devices is often significantly different from that in identical but macroscopic geometries. Plasmon effects and near-field diffraction can enhance the local field intensity, leading to enhanced cross section for light absorption and scattering, which can be utilized in substrate-enhanced spectroscopies for the detection of trace amounts of adsorbed chemicals. A specific problem is an ingenious but only empirically described way to enhance signal intensity in Raman spectroscopy by the use of a substrate patterned with gold coated micron size pyramidal pits. While Raman enhancement on nanostructured substrates is generally attributed to surface plasmons, here the micron size, and thus the sub-wavelength to near-wavelength dimensions suggest that resonant enhancement emanating from optical near-field diffraction might also play a role. To answer this question, light diffraction in a projection of the pyramidal pit: a V-groove, was modelled with a modified Neerhoff-Mur formalism suitable to calculate electromagnetic field distribution in sub-wavelength structures. Under the boundary conditions a perfect conductor screen was assumed, which excludes plasmon effects. The calculations show that interference in the cavity causes a modest resonant increase in local intensity and that near-field diffraction strongly influences the field distribution, which is explained with the electrodynamic edge effect. The magnitude of the resonant electric field on its own cannot account for the experimentally observed Raman enhancement. However, a resonant enhancement of a similar magnitude is expected for the emitted Stokes frequencies. In this case the geometry implements the conditions for the classical electromagnetic Raman enhancement, ~E(4), in a good agreement with experimental results.     

Comparing the electronic properties of the low-spin cyano-ferric [Fe(N4)(Cys)] active sites of superoxide reductase and p450cam using ENDOR spectroscopy and DFT calculations

Superoxide reductase (SOR) and P450 enzymes contain similar [Fe(N)4(SCys)] active sites and, although they catalyze very different reactions, are proposed to involve analogous low-spin (hydro)peroxo-Fe(III) intermediates in their respective mechanisms that can be modeled by cyanide binding. The equatorial FeN4 ligation by four histidine ligands in CN-SOR and the heme in CN-P450cam is directly compared by 14N ENDOR, while the axial Fe-CN and Fe-S bonding is probed by 13C ENDOR of the cyanide ligand and 1Hbeta ENDOR measurements to determine the spin density delocalization onto the cysteine sulfur. There are small, but notable, differences in the bonding between Fe(III) and its ligands in the two enzymes. The ENDOR measurements are complemented by DFT computations that support the semiempirical equation used to compute spin densities on metal-coordinated cysteinyl and shed light on bonding changes as the Fe-C-N linkage bends. They further indicate that H bonds to the cysteinyl thiolate sulfur ligand reduce the spin density on the sulfur in both active sites to a degree that exceeds the difference induced by the alternative sets of "in-plane" nitrogen ligands.     

Maximising rovibrational assignments in the ν1 band of NSCl by spectral analysis by subtraction of simulated intensities (SASSI)

The ν₁ band of thiazyl chloride (NSCl) has been measured by high resolution FTIR spectroscopy, and studied using the spectral analysis by subtraction of simulated intensities (SASSI)-technique. This involves assignment and fitting of a component, creating a global simulation of the entire component to be subtracted from the experimental spectrum, and subsequent analysis of the resultant spectrum. Through this iterative procedure it has been possible to assign rovibrational transitions for seven components, with populations in the lower vibrational state as low as 2% of the total. The components are the fundamental and , , and hotbands of ¹⁴N³²S³⁵Cl, the fundamental and hotband of ¹⁴N³²S³⁷Cl, and the fundamental of ¹⁴N³⁴S³⁵Cl. Rotational and centrifugal distortion constants beyond the quartic level have been obtained by fitting to Watson’s A-reduced Hamiltonian in upper and lower states.     

Sequential Injection/Mid-Infrared Spectroscopic Analysis of an Acetone‐Butanol‐Ethanol Fermentation: Analyte Cross‐Correlation Effects

Mid‐infrared spectroscopy together with sequential injection analysis (SIA) and partial least squares (PLS) regression analysis was used to monitor acetone‐butanol‐ethanol (ABE) fermentations under different fermentation conditions. Five analytes were simultaneously predicted (acetone, acetate, butyrate, n‐butanol, and glucose). In order to compare the overall model prediction ability, a relative average of the root mean square error of prediction (RMSEP) across all five analytes was employed. To form a PLS model devoid of any cross‐correlations between analytes, a synthetic calibration data set was created by the SIA system. As a test of their robustness, PLS models from synthetic samples and those from real fermentation samples were compared and used to predict samples from the opposite data set and from independent “acid‐crash” fermentations. The PLS model developed from the synthetic samples proved to be far more robust and accurate and used fewer factors than PLS models from the real fermentations, which were found to contain analyte cross‐correlations. The use of synthetic data enabled more accurate selection of factors and showed the importance of investigating spectral regression coefficients plots to aid and confirm appropriate factor selection. In addition, an alternative method of factor selection was proposed, using a “similarity measure” between the regression coefficient plots of factors for certain analytes and their standard spectra. Predictions using this method of factor selection over the common “minimum from an error vs. factor” plot proved to be more accurate and used far fewer factors.     

A programmable "build-couple" approach to the synthesis of heterofunctionalized polyvalent molecules

A maximally divergent "build-couple" synthesis of heterofunctionalized polyvalent molecules is described. This strategic approach enables the synthesis of highly diverse polyvalent structures from a pre-programmed combinatorial set of modules.     

Tuning through-space interactions via the secondary coordination sphere of an artificial metalloenzyme leads to enhanced Rh(iii)-catalysis

We report computationally-guided protein engineering of monomeric streptavidin Rh(iii) artificial metalloenzyme to enhance catalysis of the enantioselective coupling of acrylamide hydroxamate esters and styrenes. Increased TON correlates with calculated distances between the Rh(iii) metal and surrounding residues, underscoring an artificial metalloenzyme''s propensity for additional control in metal-catalyzed transformations by through-space interactions.

We report computationally-guided protein engineering of monomeric streptavidin Rh(iii) artificial metalloenzyme to enhance catalysis of the enantioselective coupling of acrylamide hydroxamate esters and styrenes.

Artificial metalloenzymes (ArMs) can be made by anchoring a non-natural (metal) cofactor into a protein scaffold, with the goal of imbuing new-to-nature reactivity.¹ One of the most common ArM platforms is the biotin-tetrameric(strept)avidin (biotin-tSav) system pioneered by Whitesides and Ward.^2,3 These ArMs utilize high-affinity (up to K_D ∼10⁻¹⁴ M) interactions between tSav and biotin–metal conjugates. tSav-based ArMs have appeared in an increasing number of transition-metal catalyzed transformations.^4–6 In collaboration with the Ward group, we have previously described a tetrameric streptavidin (tSav) system containing a biotinylated Rh(iii) cofactor for the asymmetric synthesis of dihydroisoquinolones using benzhydroxamate esters and acrylate partners.⁷ Monomeric streptavidin (mSav), a streptavidin/rhizavidin hybrid designed to resist tetramerization, retains its high affinity for biotin (K_D ∼10⁻⁹ M).^8,9 We recently described the use of mSav as a new ArM,¹⁰ whose simpler topology encourages protein engineering via a site-directed mutagenesis approach.Traditional manipulation of a metal''s reactivity has been accomplished by modification of the electronic and steric properties of the bound ligands (Fig. 1a).^11,12 For example, we have documented and parsed the impact of Cp electronics and sterics on a number of Rh(iii) catalyzed transformations, by structural changes to the ligand in the primary coordination sphere of Rh.¹³ On the other hand, ArMs have traditionally been used as modifiers of a metal''s steric environment largely focusing on inducing asymmetry in the bond-forming events. Less broadly appreciated is the fact that any mutations in residues proximal to the active site may also impact the metal''s electronic properties via changes to the secondary coordination sphere (Fig. 1b), with the prospect of delivering more active catalysts for a given transformation.Open in a separate windowFig. 1Methods to modify the (a) primary and (b) secondary coordination sphere of a Rh(iii) catalyst.Previously, we described a mSav·Rh(iii) catalyst and demonstrated its use in the direct enantioselective coupling of acrylamide hydroxamate esters and styrenes.¹⁰ The reaction allows rapid access to piperidines – the most common N-heterocycle found in FDA-approved pharmaceuticals.¹⁴ One of the most interesting aspects of this reaction was our observation of a 7-fold increase in turnover number (TON) by embedding the cofactor into mSav''s active site.¹⁵ It has been a long-standing goal of ArMs to not only enable new-to-nature reactivity, but also for them to achieve the stellar kinetics of a native metalloenzyme. As these systems lack the evolutionary privilege of a natural metalloenzyme, extensive mutation of the protein scaffold may be required to find the optimal environment of the metal cofactor.Predicting the effects of specific mutations can prove very challenging, as any alterations to the protein conformation and charge distribution can impact reactivity regardless of the mutation''s distance from the active site.^16–19 In order to design a better mutant, we embarked on a collaborative experimental and computational study to define the role of the protein scaffold and how single point mutations affect reactivity. We identified two key residues that play a pivotal role in mSav·Rh(iii) ArM''s secondary coordination sphere, and have used this insight to design a more active mutant.For the purposes of this study, we focused on the mSav·Rh(iii) ArM-catalyzed coupling of methacrylamide with 4-methoxystyrene as our model reaction (Fig. 2a). Using a small model of the catalyst, the lowest energy pathway of this reaction''s proposed mechanism was generated (Fig. S9^†). The calculations were performed in Turbomole^20–32 with the M06 density functional.³³ Geometries were optimized with the def2-SVP basis set, and final electronic energies were calculated with the def2-TZVP basis set.³⁴ The conductor-like screening model (COSMO)³⁵ was used as implicit solvent with a dielectric of 80 to simulate water. These calculations predicted similar barriers for the N–H activation, the C–H activation, and the migratory insertion (differences less than 3 kcal mol⁻¹). Isotope-exchange experiments revealed that the C–H activation step is reversible, implicating the migratory insertion step as turnover-limiting.¹⁰Open in a separate windowFig. 2(a) Model transformation. (b) Snapshot of the transition state for alkene insertion illustrating key nearby residues Y112 (red), E124 (blue), and S119 (purple). (c) Computed barrier to alkene insertion in the presence and absence of phenol and acetate (shown in blue).The Cp* moiety of the Cp*^biotinRhX₂ cofactor is non-covalently localized in the active site likely due to a π–π stacking interaction with Y112 (Fig. 2b). This assignment is supported by the observation that mutant Y112A leads to lower yield and enantioselectivity.¹⁰ We hypothesized that we could further manipulate both the sterics and electronics of the Cp* moiety by either directly mutating Y112 or indirectly by mutating other residues that affect the Y112-Cp* interaction.To generate a model of mSav''s protein scaffold and active-site we used QM/DMD³⁶ – a hybrid quantum mechanics/molecular mechanics method that simulates proteins piecewise. Discrete molecular dynamics (DMD) equilibrates the entire system except for the metal and part of the substrate.³⁶ After a trajectory of ∼0.5 ns, quantum mechanics (QM) is used to optimize the metal region plus sidechains and residues immediately surrounding it. This process is repeated, providing efficient sampling of the entire protein scaffold while treating the metal environment quantum-mechanically. For this study, the migratory insertion transition state was modeled in WT by freezing the coordinates of the rhodium atom and the two carbon atoms forming a bond. For each system, five replicate simulations were run for ∼20 ns each.Residues E124 and S119 both hydrogen bond to Y112 and are in close proximity to the RhCp* catalytic site (Fig. 2b).³⁷ To estimate the electronic effects of these three residues on the reaction, an acetate ion, methanol molecule, and 4-methylphenol (p-cresol) molecule were added to a small catalyst model without constraints but initially positioned to mimic the sidechains of these residues (Fig. 2c). The migratory insertion energy barrier decreases by 2 kcal mol⁻¹ with incorporation of the three residues. However, this energy barrier decreases by an additional 3 kcal mol⁻¹ upon the deletion of the methanol molecule representing S119. Not only does this imply that these amino acid sidechains may be the primary reason for the increased activity of the protein-installed catalyst, but also suggest that a longer Y112–S119 distance is favorable, so long as no water can insert in this region and replace S119 in its H-bond with Y112. We hypothesize that the carboxylate group of E124 acts as a hydrogen bond acceptor, donating electron density to the Y112 phenol ring, which in turn donates electron density to the catalyst via π–π charge transfer. This could enhance the electron donation of the metal and decrease the energy barrier to the migratory insertion step. On the other hand, S119 acts as a hydrogen bond donor which would remove electron density from Y112 and subsequently the Rh(iii) moiety.Unfortunately, mutation of Y112 (Y112F and Y112W) results in negligible protein yields. We thus identified three flanking residues (T111, E113, H87) that may be expected to have a significant impact on Y112''s position, and one distal (T32) residue, chosen as distal mutations sometimes have significant impact (Fig. 3). Through this subset of mutants, we attempted to increase TON and establish a correlation between the Y112–Rh distance and Y112–S119 distance of the mutants and their reactivities.Open in a separate windowFig. 3Structure of mSav from two different views highlighting some of the mutated residues including their TON and enantioselectivity.We used QM/DMD to simulate a representative set of these mutants spanning a wide range of TONs measured in the experiment. The Y112–Rh and Y112–S119 distances were measured every ∼0.5 ps for every simulation. The results can be represented by a 3-dimensional plot with Y112–Rh distance on the X axis, Y112–S119 distance on the Y axis, and probability density on the Z axis (Fig. 4). We find the best correlation between TON and probability density in the conformational region where the Y112–Rh distance is the shortest and the Y112–S119 interaction is not energetically relevant.³⁸Open in a separate windowFig. 4Three-dimensional probability distributions from select mutants by simultaneous sampling of Rh–Y112 and S119–Y112 distances. Probabilities for the outlined regions are also shown.To clarify this correlation, we calculated the probability of having a Y112–S119 distance between 3.5–6 Å and a Y112–Rh distance less than 5.65 Å. This Y112–S119 distance corresponds to negligible hydrogen bonding.³⁹ Additionally, we constrained the small model catalyst shown in Fig. 5b (ref. 40) and calculated the corresponding energy barriers at different Y112–Rh distances (Fig. 5a). Since rate increases exponentially as the barrier decreases,³⁵ differences in probabilities in the region where the Y112–Rh distance is between 5.4–5.65 Å have the greatest impact on the relative TONs of our model methacrylamide styrene coupling. We conclude that mutants with increasing probability in this region provide increasing TON.Open in a separate windowFig. 5a) Theoretical dependence of migratory insertion barrier on Rh-phenol distance. (b) Small-model catalyst.Theoretically, a shorter Y112–Rh distance relative to WT would result in increased reactivity. Residue G49 is located under the Rh(iii) moiety (Fig. 6). We hypothesize that by mutating the glycine into an alanine, steric congestion would force the biotinylated Rh(iii) cofactor to shift upwards closer to the electron donating phenol side chain of residue Y112. Analyzing the critical portions of the Y112–Rh and Y112–S119 distances in tandem reveals that G49A has the highest probability density in this region (Fig. 4). Indeed, experimentally, this mutant gives 97 TON and 91% ee (Fig. 6). The combination of a short Y112–Rh distance and long Y112–S119 distance leads to an increase in reactivity. This is an approximate 3-fold improvement in the TON relative to WT. The G49A mutant serves as an experimental proof of concept that a computational analysis of an ArMs secondary coordination sphere can lead to the design of a more efficient ArM.Open in a separate windowFig. 6Snapshot of the transition state for alkene insertion highlighting the position of G49 (purple) relative to Rh. Y112 is shown in red and E124 is shown in blue.In summary, we have identified three key residues that contribute to accelerating the rate of a Rh(iii)-catalyzed reaction by electronic communication to the metal via the secondary coordination sphere. E124 hydrogen bonds to Y112 transferring electron density via π–π charge transfer, an effect that is attenuated by hydrogen bonding from S119. Optimal interaction of these residues can be described computationally by finding mutants that have multiple conformations bearing short Y112–Rh distances coupled with negligible bonding between Y112 and S119. This hypothesis was experimentally verified by a mutant that enforces a closer Y112–Rh distance leading to improved TON. This result demonstrates the use of a hypothesis-based site-directed mutagenesis of the secondary sphere residues, to optimize the metal''s electronic environment within the protein scaffold and enhance an ArM''s activity.