Demonstrating the importance of porcine reproductive and respiratory syndrome virus papain-like protease 2 deubiquitinating activity in viral replication by structure-guided mutagenesis.

Deubiquitination of cellular substrates by viral proteases is a mechanism used to interfere with host cellular signaling processes, shared between members of the coronavirus- and arterivirus families. In the case of Arteriviruses, deubiquitinating and polyprotein processing activities are accomplished by the virus-encoded papain-like protease 2 (PLP2). Several studies have implicated the deubiquitinating activity of the porcine reproductive and respiratory syndrome virus (PRRSV) PLP2 in the downregulation of cellular interferon production, however to date, the only arterivirus PLP2 structure described is that of equine arteritis virus (EAV), a distantly related virus. Here we describe the first crystal structure of the PRRSV PLP2 domain both in the presence and absence of its ubiquitin substrate, which reveals unique structural differences in this viral domain compared to PLP2 from EAV. To probe the role of PRRSV PLP2 deubiquitinating activity in host immune evasion, we selectively removed this activity from the domain by mutagenesis and found that the viral domain could no longer downregulate cellular interferon production. Interestingly, unlike EAV, and also unlike the situation for MERS-CoV, we found that recombinant PRRSV carrying PLP2 DUB-specific mutations faces significant selective pressure to revert to wild-type virus in MARC-145 cells, suggesting that the PLP2 DUB activity, which in PRRSV is present as three different versions of viral protein nsp2 expressed during infection, is critically important for PRRSV replication.

Biochemical Insights into Imipenem Collateral Susceptibility Driven by ampC Mutations Conferring Ceftolozane/Tazobactam Resistance in Pseudomonas aeruginosa.

Several Pseudomonas aeruginosa AmpC mutants have emerged that exhibit enhanced activity against ceftazidime and ceftolozane, while also evading inhibition by avibactam. Interestingly, P. aeruginosa strains harboring these AmpC mutations fortuitously exhibit enhanced carbapenem susceptibility. This acquired susceptibility was investigated by comparing the degradation of imipenem by wild-type and cephalosporin-resistant AmpC. We show that cephalosporin-resistant AmpC enzymes lose their efficacy for hydrolyzing imipenem and suggest that this may be due to their increased flexibility and dynamics relative to the wild type.

Does Urea Preferentially Interact with Amide Moieties or Nonpolar Sidechains? A Question Answered Through a Judicious Selection of Model Systems.

The transfer model suggests that urea unfolds proteins mainly by increasing the solubility of the amide backbone, probably through urea-induced increase in hydrogen bonding. Other studies suggest that urea addition increases the magnitude of solvent-solute van der Waals interactions, which increases the solubility of nonpolar sidechains. More recent analyses hypothesize that urea has a similar effect in increasing the solubility of backbone and sidechain groups. In this work, we compare the effects of urea addition on the solvation of amides and alkyl groups. At first, we study the effects of urea addition upon solvent hydrogen bonding acidity and basicity through the perturbation in the fluorescence spectrum of probes 1-AN and 1-DMAN. Our results demonstrate that the solvent's hydrogen bonding properties are minimally affected by urea addition. Subsequently, we show that urea addition does not perturb the intra-molecular hydrogen bonding in salicylic acid significantly. Finally, we investigate how urea preferentially interacts with amide and alkyl groups moieties in water by comparing the effects of urea addition upon the solubility of acetaminophen and 4-tertbutylphenol. We show that urea affects amide and t-butyl solubility (lowers the transfer free energy of both amide (backbone) and alkyl (sidechain) groups) in a similar fashion. In other words, preferential interaction of urea with both moieties contributes to protein denaturation.

Molecular mechanism of quorum sensing inhibition in Streptococcus by the phage protein paratox.

Streptococcus pyogenes, or Group A Streptococcus, is a Gram-positive bacterium that can be both a human commensal and a pathogen. Central to this dichotomy are temperate bacteriophages that incorporate into the bacterial genome as prophages. These genetic elements encode both the phage proteins and the toxins harmful to the human host. One such conserved phage protein, paratox (Prx), is always found encoded adjacent to the toxin genes, and this linkage is preserved during all stages of the phage life cycle. Within S. pyogenes, Prx functions to inhibit the quorum-sensing receptor-signal pair ComRS, the master regulator of natural competence, or the ability to uptake endogenous DNA. However, the mechanism by which Prx directly binds and inhibits the receptor ComR is unknown. To understand how Prx inhibits ComR at the molecular level, we pursued an X-ray crystal structure of Prx bound to ComR. The structural data supported by solution X-ray scattering data demonstrate that Prx induces a conformational change in ComR to directly access its DNA-binding domain. Furthermore, electromobility shift assays and competition binding assays reveal that Prx effectively uncouples the interdomain conformational change required for activation of ComR via the signaling molecule XIP. Although to our knowledge the molecular mechanism of quorum-sensing inhibition by Prx is unique, it is analogous to the mechanism employed by the phage protein Aqs1 in Pseudomonas aeruginosa. Together, this demonstrates an example of convergent evolution between Gram-positive and Gram-negative phages to inhibit quorum-sensing and highlights the versatility of small phage proteins.

Hofmeister Effects of Group II Cations as Seen in the Unfolding of Ribonuclease A.

This work studies the effects of alkaline-earth cation addition on the unfolding free energy of a model protein, pancreatic Ribonuclease A (RNase A) by differential scanning calorimetry analysis. RNase A was chosen because: a) it does not specifically bind Mg2+ , Ca2+ and Sr2+ cations and b) maintains its structural integrity throughout a large pH range. We have measured and compared the effects of NaCl, MgCl2 , CaCl2 and SrCl2 addition on the melting point of RNase A. Our results show that even though the addition of group II cations to aqueous solvent reduces the solubility of nonpolar residues (and enhances the hydrophobic effect), their interactions with the amide moieties are strong enough to "salt-them-in" the solvent, thereby causing an overall protein stability reduction. We demonstrate that the amide-cation interactions are a major contributor to the observed "Hofmeister Effects" of group II cations in protein folding. Our analysis suggests that protein folding "Hofmeister Effects" of group II cations, are mostly the aggregate sum of how cation addition simultaneously salts-out hydrophobic moieties by increasing the cavitation free energy, while promoting the salting-in of amide moieties through contact pair formation.

Adding Insult to Injury: Mechanistic Basis for How AmpC Mutations Allow Pseudomonas aeruginosa To Accelerate Cephalosporin Hydrolysis and Evade Avibactam.

Pseudomonas aeruginosa is a leading cause of nosocomial infections worldwide and notorious for its broad-spectrum resistance to antibiotics. A key mechanism that provides extensive resistance to ß-lactam antibiotics is the inducible expression of AmpC ß-lactamase. Recently, a number of clinical isolates expressing mutated forms of AmpC have been found to be clinically resistant to the antipseudomonal ß-lactam-ß-lactamase inhibitor (BLI) combinations ceftolozane-tazobactam and ceftazidime-avibactam. Here, we compare the enzymatic activity of wild-type (WT) AmpC from PAO1 to those of four of these reported AmpC mutants, bearing mutations E247K (a change of E to K at position 247), G183D, T96I, and ΔG229-E247 (a deletion from position 229 to 247), to gain detailed insights into how these mutations allow the circumvention of these clinically vital antibiotic-inhibitor combinations. We found that these mutations exert a 2-fold effect on the catalytic cycle of AmpC. First, they reduce the stability of the enzyme, thereby increasing its flexibility. This appears to increase the rate of deacylation of the enzyme-bound ß-lactam, resulting in greater catalytic efficiencies toward ceftolozane and ceftazidime. Second, these mutations reduce the affinity of avibactam for AmpC by increasing the apparent activation barrier of the enzyme acylation step. This does not influence the catalytic turnover of ceftolozane and ceftazidime significantly, as deacylation is the rate-limiting step for the breakdown of these antibiotic substrates. It is remarkable that these mutations enhance the catalytic efficiency of AmpC toward ceftolozane and ceftazidime while simultaneously reducing susceptibility to inhibition by avibactam. Knowledge gained from the molecular analysis of these and other AmpC resistance mutants will, we believe, aid in the design of ß-lactams and BLIs with reduced susceptibility to mutational resistance.

Salt Effects on Hydrophobic Solvation: Is the Observed Salt Specificity the Result of Excluded Volume Effects or Water Mediated Ion-Hydrophobe Association?

The solubility of hydrophobic molecules in water is sensitive to salt addition in an ion-specific manner. Such "salting-out" and "salting-in" properties have been shown to be a major contributor to the measured ion-specific Hofmeister effects that are observed in many biophysical phenomena. Various theoretical models have suggested a number of disparate mechanisms for salting-out (salting-in) of hydrophobic moieties, the most popular of which include preferential interaction, water-mediated association, and electrostriction models. However, a complete molecular level description of this ion-specificity is not yet available. This work investigates the ion-specific nature of hydrophobic solvation by studying how sodium and chloride salts affect the thermodynamics of 1,2-hexanediol micellization. The results of this study are analyzed in terms of scaled-particle theory and we show that salt addition can affect hydrophobic solvation in two modalities: salt addition changes the cavitation free energy; salt addition also influences the solvent-solute interaction energy by changing the hydration of the hydrophobic solute. These two effects are salt specific in nature and we suggest that for small hydrophobic solutes these effects are the main cause of salt-specific Hofmeister effects on their solubility.

Deubiquitinase function of arterivirus papain-like protease 2 suppresses the innate immune response in infected host cells.

Protein ubiquitination regulates important innate immune responses. The discovery of viruses encoding deubiquitinating enzymes (DUBs) suggests they remove ubiquitin to evade ubiquitin-dependent antiviral responses; however, this has never been conclusively demonstrated in virus-infected cells. Arteriviruses are economically important positive-stranded RNA viruses that encode an ovarian tumor (OTU) domain DUB known as papain-like protease 2 (PLP2). This enzyme is essential for arterivirus replication by cleaving a site within the viral replicase polyproteins and also removes ubiquitin from cellular proteins. To dissect this dual specificity, which relies on a single catalytic site, we determined the crystal structure of equine arteritis virus PLP2 in complex with ubiquitin (1.45 Å). PLP2 binds ubiquitin using a zinc finger that is uniquely integrated into an exceptionally compact OTU-domain fold that represents a new subclass of zinc-dependent OTU DUBs. Notably, the ubiquitin-binding surface is distant from the catalytic site, which allowed us to mutate this surface to significantly reduce DUB activity without affecting polyprotein cleavage. Viruses harboring such mutations exhibited WT replication kinetics, confirming that PLP2-mediated polyprotein cleavage was intact, but the loss of DUB activity strikingly enhanced innate immune signaling. Compared with WT virus infection, IFN-ß mRNA levels in equine cells infected with PLP2 mutants were increased by nearly an order of magnitude. Our findings not only establish PLP2 DUB activity as a critical factor in arteriviral innate immune evasion, but the selective inactivation of DUB activity also opens unique possibilities for developing improved live attenuated vaccines against arteriviruses and other viruses encoding similar dual-specificity proteases.

Conformational itinerary of Pseudomonas aeruginosa 1,6-anhydro-N-acetylmuramic acid kinase during its catalytic cycle.

Anhydro-sugar kinases are unique from other sugar kinases in that they must cleave the 1,6-anhydro ring of their sugar substrate to phosphorylate it using ATP. Here we show that the peptidoglycan recycling enzyme 1,6-anhydro-N-acetylmuramic acid kinase (AnmK) from Pseudomonas aeruginosa undergoes large conformational changes during its catalytic cycle, with its two domains rotating apart by up to 32° around two hinge regions to expose an active site cleft into which the substrates 1,6-anhydroMurNAc and ATP can bind. X-ray structures of the open state bound to a nonhydrolyzable ATP analog (AMPPCP) and 1,6-anhydroMurNAc provide detailed insight into a ternary complex that forms preceding an operative Michaelis complex. Structural analysis of the hinge regions demonstrates a role for nucleotide binding and possible cross-talk between the bound ligands to modulate the opening and closing of AnmK. Although AnmK was found to exhibit similar binding affinities for ATP, ADP, and AMPPCP according to fluorescence spectroscopy, small angle x-ray scattering analyses revealed that AnmK adopts an open conformation in solution in the absence of ligand and that it remains in this open state after binding AMPPCP, as we had observed for our crystal structure of this complex. In contrast, the enzyme favored a closed conformation when bound to ADP in solution, consistent with a previous crystal structure of this complex. Together, our findings show that the open conformation of AnmK facilitates binding of both the sugar and nucleotide substrates and that large structural rearrangements must occur upon closure of the enzyme to correctly align the substrates and residues of the enzyme for catalysis.

Structural basis for the removal of ubiquitin and interferon-stimulated gene 15 by a viral ovarian tumor domain-containing protease.

The attachment of ubiquitin (Ub) and the Ub-like (Ubl) molecule interferon-stimulated gene 15 (ISG15) to cellular proteins mediates important innate antiviral responses. Ovarian tumor (OTU) domain proteases from nairoviruses and arteriviruses were recently found to remove these molecules from host proteins, which inhibits Ub and ISG15-dependent antiviral pathways. This contrasts with the Ub-specific activity of known eukaryotic OTU-domain proteases. Here we describe crystal structures of a viral OTU domain from the highly pathogenic Crimean-Congo haemorrhagic fever virus (CCHFV) bound to Ub and to ISG15 at 2.5-Å and 2.3-Å resolution, respectively. The complexes provide a unique structural example of ISG15 bound to another protein and reveal the molecular mechanism of an ISG15 cross-reactive deubiquitinase. To accommodate structural differences between Ub and ISG15, the viral protease binds the ß-grasp folds of Ub and C-terminal Ub-like domain of ISG15 in an orientation that is rotated nearly 75° with respect to that observed for Ub bound to a representative eukaryotic OTU domain from yeast. Distinct structural determinants necessary for binding either substrate were identified and allowed the reengineering of the viral OTU protease into enzymes with increased substrate specificity, either for Ub or for ISG15. Our findings now provide the basis to determine in vivo the relative contributions of deubiquitination and deISGylation to viral immune evasion tactics, and a structural template of a promiscuous deubiquitinase from a haemorrhagic fever virus that can be targeted for inhibition using small-molecule-based strategies.

The Streptococcus phage protein paratox is an intrinsically disordered protein.

The bacteriophage protein paratox (Prx) blocks quorum sensing in its streptococcal host by directly binding the signal receptor and transcription factor ComR. This reduces the ability of Streptococcus to uptake environmental DNA and protects phage DNA from damage by recombination. Past work characterizing the Prx:ComR molecular interaction revealed that paratox adopts a well-ordered globular fold when bound to ComR. However, solution-state biophysical measurements suggested that Prx may be conformationally dynamic. To address this discrepancy, we investigated the stability and dynamic properties of Prx in solution using circular dichroism, nuclear magnetic resonance, and several fluorescence-based protein folding assays. Our work shows that under dilute buffer conditions Prx is intrinsically disordered. We also show that the addition of kosmotropic salts or protein stabilizing osmolytes induces Prx folding. However, the solute stabilized fold is different from the conformation Prx adopts when it is bound to ComR. Furthermore, we have characterized Prx folding thermodynamics and folding kinetics through steady-state fluorescence and stopped flow kinetic measurements. Our results show that Prx is a highly dynamic protein in dilute solution, folding and refolding within the 10 ms timescale. Overall, our results demonstrate that the streptococcal phage protein Prx is an intrinsically disordered protein in a two-state equilibrium with a solute-stabilized folded form. Furthermore, the solute-stabilized fold is likely the predominant form of Prx in a solute-crowded bacterial cell. Finally, our work suggests that Prx binds and inhibits ComR, and thus quorum sensing in Streptococcus, by a combination of conformational selection and induced-fit binding mechanisms.

The effect of urea on aqueous hydrophobic contact-pair interactions.

Urea is perhaps the most common denaturant used for studying proteins. However the mechanism of denaturation is still not well understood. Recent theoretical work suggests that van der Waals interactions between urea and non-polar amino acid residues are a major contributor to the protein denaturation process. However, there are few experimental data measuring the effect of urea on hydrophobic interactions. In this work we have determined how the addition of urea to the aqueous solvent affects the contact-pair formed between alkyl and phenyl groups in model compounds: the data indicate that for solutes having a radius smaller than 2.87 Å cavity formation energetics dominate, therefore the addition of urea promotes the formation of hydrophobic contact pairs; while for larger solutes, van der Waals interactions have the largest magnitude, causing urea to disrupt the formation of contact pairs. The influence of urea on hydrophobic interactions is shown to be continuous in the 1-8 M concentration range and is well-correlated with the predictions of scaled particle theory. This demonstrates that the effect of urea on hydrophobic contact pairs can be explained by the changes observed in the solvent packing density, without having to invoke changes in the hydrogen bonding network of water.

Increased phosphorylation of HexM improves lysosomal uptake and potential for managing GM2 gangliosidoses.

Tay-Sachs and Sandhoff diseases are genetic disorders resulting from mutations in HEXA or HEXB, which code for the α- and ß-subunits of the heterodimer ß-hexosaminidase A (HexA), respectively. Loss of HexA activity results in the accumulation of GM2 ganglioside (GM2) in neuronal lysosomes, culminating in neurodegeneration and death, often by age 4. Previously, we combined critical features of the α- and ß-subunits of HexA into a single subunit to create a homodimeric enzyme known as HexM. HexM is twice as active as HexA and degrades GM2 in vivo, making it a candidate for enzyme replacement therapy (ERT). Here we show HexM production is scalable to meet ERT requirements and we describe an approach that enhances its cellular uptake via co-expression with an engineered GlcNAc-1-phosphotransferase that highly phosphorylates lysosomal proteins. Further, we developed a HexA overexpression system and functionally compared the recombinant enzyme to HexM, revealing the kinetic differences between the enzymes. This study further advances HexM as an ERT candidate and provides a convenient system to produce HexA for comparative studies.

Antiseptic quaternary ammonium compound tolerance by gram-negative bacteria can be rapidly detected using an impermeant fluorescent dye-based assay.

Biocides such as quaternary ammonium compounds (QACs) are potentially important contributors towards bacterial antimicrobial resistance development, however, their contributions are unclear due to a lack of internationally recognized biocide testing standards. Methods to detect QAC tolerance are limited to laborious traditional antimicrobial susceptibility testing (AST) methods. Here, we developed a rapid fluorescent dye-based membrane impermeant assay (RFDMIA) to discriminate QAC susceptibility among Gram-negative Enterobacterales and Pseudomonadales species. RFDMIA uses a membrane impermeant fluorescent dye, propidium iodide, in a 30-min 96-well fluorescent microplate-based assay where cell suspensions are exposed to increasing QAC concentrations. Our results demonstrate that RFDMIA can discriminate between QAC-susceptible and QAC-adapted Escherichia coli tolerant phenotypes and predict benzalkonium and cetrimide tolerance in all species tested except for intrinsically fluorescent Pseudomonas aeruginosa. RFDMIA identified a close association to minimum inhibitory concentration values determined by broth microdilution AST and increasing fluorescent dye emission values. RFDMIA emission values and scanning electron microscopy results also suggest that CET-adapted E. coli isolates have a CET dependence, where cells require sub-inhibitory CET concentrations to maintain bacilliform cell integrity. Overall, this study generates a new, rapid, sensitive fluorescent assay capable of detecting QAC-susceptible Gram-negative bacteria phenotypes and cell membrane perturbations.

Do soft anions promote protein denaturation through binding interactions? A case study using ribonuclease A.

It has long been known that large soft anions like bromide, iodide and thiocyanate are protein denaturing agents, but their mechanism of action is still unclear. In this work we have investigated the protein denaturing properties of these anions using Ribonuclease A (RNase A) as a model protein system. Salt-induced perturbations to the protein folding free energy were determined using differential scanning calorimetry and the results demonstrate that the addition of sodium iodide and sodium thiocyanate significantly decreases the melting temperature of the protein. In order to account for this reduction in protein stability, we show that the introduction of salts that contain soft anions to the aqueous solvent perturbs the protein unfolding free energy through three mechanisms: (a) screening Coulomb interactions that exist between charged protein residues, (b) Hofmeister effects, and (c) specific anion binding to CH and CH2 moieties in the protein polypeptide backbone. Using the micellization of 1,2-hexanediol as a ruler for hydrophobicity, we have devised a practical methodology that separates the Coulomb and Hofmeister contributions of salts to the protein unfolding free energy. This allowing us to isolate the contribution of soft anion binding interactions to the unfolding process. The analysis shows that binding contributions have the largest magnitude, confirming that it is the binding of soft anions to the polypeptide backbone that is the main promoter of protein unfolding.

Rapid mixing of colliding picoliter liquid droplets delivered through-space from piezoelectric-actuated pipettes characterized by time-resolved fluorescence monitoring.

Rapid mixing of aqueous solutions is a crucial first step to study the kinetics of fast biochemical reactions with high temporal resolution. Remarkable progress toward this goal has been made through the development of advanced stopped-flow mixing techniques resulting in reduced dead times, and thereby extending reaction monitoring capabilities to numerous biochemical systems. Concurrently, piezoelectric actuators for through-space liquid droplet sample delivery have also been applied in several experimental systems, providing discrete picoliter sample volume delivery and precision sample deposition onto a surface, free of confinement within microfluidic devices, tubing, or other physical constraints. Here, we characterize the inertial mixing kinetics of two aqueous droplets (130 pl) produced by piezoelectric-actuated pipettes, following droplet collision in free space and deposition on a surface in a proof of principle experiment. A time-resolved fluorescence system was developed to monitor the mixing and fluorescence quenching of 5-carboxytetramethylrhodamine (5-Tamra) and N-Bromosuccinimide, which we show to occur in less than 10 ms. In this respect, this methodology is unique in that it offers millisecond mixing capabilities for very small quantities of discrete sample volumes. Furthermore, the use of discrete droplets for sample delivery and mixing in free space provides potential advantages, including the elimination of the requirement for a physical construction as with microfluidic systems, and thereby makes possible and extends the experimental capabilities of many systems.

Proteinaceous Nano container Encapsulate Polycyclic Aromatic Hydrocarbons.

Polycyclic aromatic hydrocarbons (PAHs) are toxic, mutagenic and among the most damaging chemical compounds with regard to living organisms. Because of their persistence and wide distribution removal from the environment is an important challenge. Here we report a new Nano container matrix based on the deep sea archaea-derived RHCC-Nanotube (RHCC-NT), which rapidly and preferentially binds low molecular weight PAHs. Under controlled-laboratory conditions and using fluorescence spectroscopy in combination with X-ray crystallography and MD simulations, we quantified the real-time binding of low molecular weight PAHs (2-4 rings) to our substrate. Binding coefficients ranged from 5.4 ± 1.6 (fluorene) to 32 ± 7.0 µM (acenaphthylene) and a binding capacity of 85 pmoles PAH per mg RHCC-NT, or 2.12 µmoles in a standard 25 mg sampler. The uptake rate of pyrene was calculated to be 1.59 nmol/hr∙mol RHCC-NT (at 10  C). Our results clearly show that RHCC-NT is uniquely suited as a monitoring matrix for low molecular weight PAHs.

Effects of the protein denaturant guanidinium chloride on aqueous hydrophobic contact-pair interactions.

Guanidinium chloride (GdmCl) is one of the most common protein denaturants. Although GdmCl is well known in the field of protein folding, the mechanism by which it denatures proteins is not well understood. In fact, there are few studies looking at its effects on hydrophobic interactions. In this work the effect of GdmCl on hydrophobic interactions has been studied by observing how the denaturant influences model systems of phenyl and alkyl hydrophobic contact pairs. Contact pair formation is monitored through the use of fluorescence spectroscopy, i.e., measuring the intrinsic phenol fluorescence being quenched by carboxylate ions. Hydrophobic interactions are isolated from other interactions through a previously developed methodology. The results show that GdmCl does not significantly affect hydrophobic interactions between small moieties such as methyl groups and phenol; while on the other hand, the interaction of larger hydrophobes such as hexyl and heptyl groups with phenol is significantly destabilized.

Melittin as model system for probing interactions between proteins and cyclodextrins.

Cylcodextrin sugars are cyclic sugars that have a hydrophilic exterior and a hydrophobic center. This enables cyclodextrins to solubilize hydrophobic molecules in aqueous media. Cyclodextrins may inhibit aggregation by intercalating surface aromatic residues and competing with interprotein aromatic clusters (pi-pi interactions). In order to investigate this concept, the interaction of hydroxypropyl-beta-cyclodextrin (HPBCD) with melittin is studied with steady-state and time-resolved fluorescence, fluorescence polarization, circular dichroism, and IR spectroscopy. HPBCD inhibits the aggregation of melittin. This inhibition and the spectroscopic results are consistent with the lone aromatic tryptophan of the peptide being intercalated within HPBCD.

Accessibility of oxygen with respect to the heme pocket in horseradish peroxidase.

Oxygen and other molecules of similar size take part in a variety of protein reactions. Therefore, it is critical to understand how these small molecules penetrate the protein matrix. The protein system studied in this case is horseradish peroxidase (HRP). We have converted the native HRP into a phosphorescent analog by replacing the heme prosthetic group by Pd-mesoporphyrin. Oxygen readily quenches the phosphorescence of Pd porphyrins, and this can be used to characterize oxygen diffusion through the protein matrix. Our measurements indicate that solvent viscosity and pH modulate the accessibility of the heme pocket relative to small molecules. The binding of the substrate benzohydroxamic acid (BHA) to the protein drastically impedes oxygen access to the heme pocket. These results indicate that, first, the penetration of small molecules through the protein matrix is a function of protein dynamics, and second, there are specific pathways for the diffusion of these molecules. The effect of substrate and pH on protein dynamics has been investigated with the use of molecular dynamics calculations. We demonstrate that the model of a "fluctuating entry point," as suggested by Zwanzig (J Chem Phys 1992;97:3587-3589), properly describes the diffusion of oxygen through the protein matrix.