Molecular simulations reveal a new entry site in quercetin 2,3-dioxygenase. A pathway for dioxygen?

Molecular Dynamics and Free Energy Simulations of Phenylacetate and CO2 Release from AMDase and Its G74C/C188S Mutant: A Possible Rationale for the Reduced Activity of the Latter

Arylmalonate decarboxylase (AMDase) catalyzes the decarboxylation of α-aryl-α-methyl malonates to produce optically pure α-arylpropionates of industrial and medicinal importance. Herein, atomistic molecular dynamics simulations have been carried out to delineate the mechanism of the release of product molecules phenylacetate (PAC) and carbon dioxide (CO₂), from the wild-type (WT) and its G74C/C188S mutant enzymes. Both of the product molecules follow a crystallographically characterized solvent-accessible channel to come out of the protein interior. A higher free energy barrier for the release of PAC from G74C/C188S compared to that in the WT is consistent with the experimentally observed compromised efficiency of the mutant. The release of CO₂ precedes that of PAC; free energy barriers for CO₂ and PAC release in the WT enzyme are calculated to be ∼1-2 and ∼23 kcal/mol, respectively. Postdecarboxylation, CO₂ moves toward a hydrophobic pocket formed by Pro 14, Leu 38, Leu 40, Leu 77, and the side chain of Tyr 48 which serves as its temporary "reservoir". CO₂ releases following a channel mainly decorated by apolar residues, unlike in the case of oxalate decarboxylase where polar residues mediate its transport.

Kinetics and thermodynamics of gas diffusion in a NiFe hydrogenase

We have investigated O₂ and H₂ transport across a NiFe hydrogenase at the atomic scale by means of computational methods. The Wild Type protein has been compared with the V74Q mutant. Two distinct methodologies have been applied to study the gas access to the active site. Temperature locally enhanced sampling simulations have emphasized the importance of protein dynamics on gas diffusion. The O₂ diffusion free energy profiles, obtained by umbrella sampling, are in agreement with the known kinetic data and show that in the V74Q mutant, the inhibition process is lowered from both a kinetic and a thermodynamic point of view.

Finding molecular dioxygen tunnels in homoprotocatechuate 2,3-dioxygenase: implications for different reactivity of identical subunits

Extradiol dioxygenases facilitate microbial aerobic degradation of catechol and its derivatives by activating molecular dioxygen and incorporating both oxygen atoms into their substrates. Experimental and theoretical studies have focused on the mechanism of the reaction at the active site. However, whether the catalytic rate is limited by O(2) access to the active site has not yet been explored. Here, we choose a recently solved X-ray structure of homoprotocatechuate 2,3-dioxygenase as a typical example to determine potential pathways for O(2) migration from the solvent into the enzyme center. On the basis of the trajectories of two 10-ns molecular dynamics simulations, implicit ligand sampling was used to calculate the 3D free energy map for O(2) inside the protein. The energetically optimal routes for O(2) diffusion were identified for each subunit of the homotetrameric protein structure. The O(2) tunnels formed because of thermal fluctuations were also characterized by connecting elongated cavities inside the protein. By superimposing the favorable O(2) tunnels on to the free energy map, both energetically and geometrically preferred O(2) pathways were determined, as also were the amino acids that may be critical for O(2) passage along these paths. Our results demonstrate that identical subunits possess quite distinct O(2) tunnels. The order of O(2) affinity of these tunnels is generally consistent with the order of the catalytic rate of each subunit. As a consequence, the probability of finding the reaction product is highest in the subunit containing the highest O(2) affinity pathway.

A Rubisco mutant that confers growth under a normally "inhibitory" oxygen concentration

Ribulose-1,5-bisphosphate (RuBP) carboxylase/oxygenase (Rubisco) is a globally significant biocatalyst that facilitates the removal and sequestration of CO2 from the biosphere. Rubisco-catalyzed CO2 reduction thus provides virtually all of the organic carbon utilized by living organisms. Despite catalyzing the rate-limiting step of photosynthetic and chemoautotrophic CO2 assimilation, Rubisco is markedly inefficient as the competition between O2 and CO2 for the same substrate limits the ability of aerobic organisms to obtain maximum amounts of organic carbon for CO2-dependent growth. Random and site-directed mutagenesis procedures were coupled with genetic selection to identify an "oxygen-insensitive" mutant cyanobacterial (Synechococcus sp. strain PCC 6301) Rubisco that allowed for CO2-dependent growth of a host bacterium at an oxygen concentration that inhibited growth of the host containing wild-type Synechococcus Rubisco. The mutant substitution, A375V, was identified as an intragenic suppressor of D103V, a negative mutant enzyme incapable of supporting autotrophic growth. Ala-375 (Ala-378 of spinach Rubisco) is a conserved residue in all form I (plant-like) Rubiscos. Structure-function analyses indicate that the A375V substitution decreased the enzyme's oxygen sensitivity (and not CO2/O2 specificity), possibly by rearranging a network of interactions in a fairly conserved hydrophobic pocket near the active site. These studies point to the potential of engineering plants and other significant aerobic organisms to fix CO2 unfettered by the presence of O2.

Biochemical and molecular characterization of a quercetinase from Penicillium olsonii

Quercetinase (quercetin 2,3-dioxygenase, EC 1.13.11.24) is produced by various filamentous fungi when grown on rutin as the sole carbon and energy source. From a rutin based liquid culture of Penicillium olsonii, we purified a quercetinase with a specific activity of 175U mg(-1). The enzyme is a monomeric glycoprotein of approximately 55 kDa, containing 0.9+/-0.1 copper atoms per protein. Its substrate specificity is restricted to the flavonol family of flavonoids. It is completely inhibited by diethyldithiocarbamate at a concentration of 100 nM and 1H-2-benzyl-3-hydroxy-4-oxoquinolin is a competitive inhibitor with a K(I) of 4 microM. The cDNA poquer1 was cloned and sequenced. It encodes a 365 amino acids long enzyme with a strong sequence identity with the Aspergillus japonicus quercetinase (Q7SIC2). Like the enzyme from A. japonicus, only one of the two cupin domains of the Penicillium olsonii quercetinase is able to bind a metal atom.

Molecular simulations bring new insights into flavonoid/quercetinase interaction modes

Molecular dynamics simulations, using the AMBER force field, were performed to study Quercetin 2,3-Dioxygenase enzyme (Quercetinase or 2,3QD). We have analyzed the structural modifications of the active site and of the linker region between the native enzyme and the enzyme-substrate complex. New structural informations, such as an allosteric effect in the presence of the substrate, as well as description of the enzyme-substrate interactions and values of binding free energies were brought out. All these results confirm the idea that the linker encloses the substrate in the active site and also enlighten the recognition role of the substrate B-ring by the enzyme. Moreover, a specific interaction scheme has been proposed to explain the relative degradation rate of various flavonoid compounds under the oxygenolysis reaction catalyzed by the Quercetin 2,3-Dioxygenase enzyme.