Atomistic Insights into Photoprotein Formation: Computational Prediction of the Properties of Coelenterazine and Oxygen Binding in Obelin.

Bioluminescence in marine systems is dominated by the use of coelenterazine for light production. The bioluminescent reaction of coelenterazine is an enzyme catalyzed oxidative decarboxylation: coelenterazine reacts with molecular oxygen to form carbon dioxide, coelenteramide, and light. One such class is the Ca2+ -regulated photoproteins. These proteins bind coelenterazine and oxygen, and trap 2-hydroperoxycoelenterazine, an intermediate along the reaction pathway. The reaction is halted until Ca2+ binding triggers the completion of the reaction. There are currently no reported experimental, atomistic descriptions of this ternary Michaelis complex. This study utilized computational techniques to develop an atomistic model of the Michaelis complex. Extensive molecular dynamics simulations were carried out to study the interactions between four tautomeric/protonation states of coelenterazine and wide-type and mutant obelin. Only minor differences in binding modes were observed across all systems. Interestingly, no basic residues were identified in the vicinity of the N7-nitrogen of coelenterazine. This observation was surprising considering that deprotonation at this position is a key mechanistic step in the proposed bioluminescent reaction. This work suggests that coelenterazine binds either as the O10H tautomer, or in the deprotonated form. Implicit ligand sampling simulations were used to identify potential O2 binding and migration pathways within obelin. A key oxygen binding site was identified close to the coelenterazine imidazopyrazinone core. The O2 binding free energy was observed to be dependent on the protonation state of coelenterazine. Taken together, the description of the obelin-coelenterazine-O2 complexes established in this study provides the basis for future computational studies of the bioluminescent mechanism. © 2019 Wiley Periodicals, Inc.

pKa cycling of the general acid/base in glycoside hydrolase families 33 and 34.

Glycoside hydrolase families 33 and 34 catalyse the hydrolysis of terminal sialic acid residues from sialyl oligosaccharides and glycoconjugates with a net retention of the stereochemistry at the anomeric centre. It is generally believed that the conserved aspartic acid in the active site functions as a general acid to protonate the hydroxyl group of the departing aglycone during glycosylation, and then as a general base to facilitate the nucleophilic attack of the water molecule on the intermediate state during the deglycosylation reaction. The dual role of the general acid/base places specific demands upon its protonation state, and thus pKa values. However, it is not fully understood how this catalytic residue can achieve such pKa cycling during catalysis. We present both MM and combined QM/MM simulations to characterise the pKa values of the proposed catalytic general acid/base in the glycoside hydrolase families 33 and 34. Collectively, our study suggests that the binding of anionic substrates and the local solvation properties along with the neutralisation of the nearby glutamic acid upon glycosylation modulate the electrostatic environment around the general acid/base to achieve its proper protonation states.

Direct observation of photodissociation products from phenylperoxyl radicals isolated in the gas phase.

Gas phase peroxyl radicals are central to our chemical understanding of combustion and atmospheric processes and are typically characterized by strong absorption in the UV (λ(max) ≈ 240 nm). The analogous maximum absorption feature for arylperoxyl radicals is predicted to shift to the visible but has not previously been characterized nor have any photoproducts arising from this transition been identified. Here we describe the controlled synthesis and isolation in vacuo of an array of charge-substituted phenylperoxyl radicals at room temperature, including the 4-(N,N,N-trimethylammonium)methyl phenylperoxyl radical cation (4-Me3N([+])CH2-C6H4OO(•)), using linear ion-trap mass spectrometry. Photodissociation mass spectra obtained at wavelengths ranging from 310 to 500 nm reveal two major photoproduct channels corresponding to homolysis of aryl-OO and arylO-O bonds resulting in loss of O2 and O, respectively. Combining the photodissociation yields across this spectral window produces a broad (FWHM ≈ 60 nm) but clearly resolved feature centered at λ(max) = 403 nm (3.08 eV). The influence of the charge-tag identity and its proximity to the radical site are investigated and demonstrate no effect on the identity of the two dominant photoproduct channels. Electronic structure calculations have located the vertical B ← X transition of these substituted phenylperoxyl radicals within the experimental uncertainty and further predict the analogous transition for unsubstituted phenylperoxyl radical (C6H5OO(•)) to be 457 nm (2.71 eV), nearly 45 nm shorter than previous estimates and in good agreement with recent computational values.