On the Pathways of Biologically Relevant Diatomic Gases through Proteins. Dioxygen and Heme Oxygenase from the Perspective of Molecular Dynamics

This work deals with dioxygen (O₂) binding sites and pathways through inducible human heme oxygenase (HO‐1). The experimentally known distal binding site 1, and sites 2–3 above it, could be reproduced by means of non‐deterministic random‐acceleration molecular‐dynamics (RAMD) simulations. In addition, RAMD revealed the proximal binding site 5, a deeply‐seated binding site 4, which lies behind heme, as well as a few gates communicating with the external medium. In getting from site 1 to the main gate, which lies on the protein front opposed to site 4, O₂ follows chiefly the shortest direct pathway. Less frequently, O₂ visits intermediate sites 2, 4, or 5 along longer pathways. A similarity between HO‐1, myoglobin, and cytoglobin in using, for diatomic gas delivery, the direct shortest pathway from the heme center to the surrounding medium, is emphasized. Otherwise, comparing other proteins and diatomic gases, each system reveals its peculiarities as to sites, gates, and pathways. Thus, relating these properties to the physiological functions of the proteins remains in general a challenge for future studies.     

On the Quest of Dioxygen by Monomeric Sarcosine Oxidase. A Molecular Dynamics Investigation

It is reported here on random acceleration molecular dynamics (RAMD) simulations with the 2GF3 bacterial monomeric sarcosine oxidase (MSOX), O₂, and furoic acid in place of sarcosine, solvated by TIP3 H₂O in a periodic box. An external tiny force, acting randomly on O₂, accelerated its relocation, from the center of activation between residue K265 and the si face of the flavin ring of the flavin adenine dinucleotide cofactor, to the surrounding solvent. Only three of the four O₂ gates previously described for this system along a composite method technique were identified, while two more major O₂ gates were found. The RAMD simulations also revealed that the same gate can be reached by O₂ along different pathways, often involving traps for O₂. Both the residence time of O₂ in the traps, and the total trajectory time for O₂ getting to the solvent, could be evaluated. The new quick pathways discovered here suggest that O₂ exploits all nearby interstices created by the thermal fluctuations of the protein, not having necessarily to look for the permanent large channel used for uptake of the FADH cofactor. To this regard, MSOX resembles closely KijD3 N‐oxygenase. These observations solicit experimental substantiation, in a long term aim at discovering whether gates and pathways for the small gaseous ligands inside the proteins are under Darwinian functional evolution or merely stochastic control operates.     

Binding Pockets and Permeation Channels for Dioxygen through Cofactorless 3‐Hydroxy‐2‐methylquinolin‐4‐one 2,4‐Dioxygenase in Association with Its Natural Substrate, 3‐Hydroxy‐2‐methylquinolin‐4(1H)‐one. A Perspective from Molecular Dynamics Simulations

This work describes an investigation of pathways and binging pockets (BPs) for dioxygen (O₂) through the cofactorless oxygenase 3‐hydroxy‐2‐methylquinolin‐4‐one 2,4‐dioxygenase in complex with its natural substrate, 3‐hydroxy‐2‐methylquinolin‐4(1H)‐one, in aqueous solution. The investigation tool was random‐acceleration molecular dynamics (RAMD), whereby a tiny, randomly oriented external force is applied to O₂ in order to accelerate its movements. In doing that, care was taken that the external force only continues, if O₂ moves along a direction for a given period of time, otherwise the force changed direction randomly. Gates for expulsion of O₂ from the protein, which can also be taken as gates for O₂ uptake, were found throughout almost the whole external surface of the protein, alongside a variety of BPs for O₂. The most exploited gates and BPs were not found to correspond to the single gate and BP proposed previously from the examination of the static model from X‐ray diffraction analysis of this system. Therefore, experimental investigations of this system that go beyond the static model are urgently needed.     

New Vistas on the Recruiting of Ferrous Iron and Dioxygen by Ferritins: A Case Study of the Escherichia coli 24‐mer Ferritin by All‐Atom Molecular Dynamics in Aqueous Medium

It is shown here that Fe²⁺ and O₂ ligands are displaced from the ferroxidase center of the C1 four‐helix bundle of E. coli 24‐mer ferritin under molecular dynamics (MD) aided by a randomly oriented external force applied to the ligand. Under these conditions, ligand egress toward the external aqueous medium occurs preferentially from the same four‐helix bundle, in the case of O₂, or other bundle, in the case of Fe². Viewing ligand egress from the protein as the microscopic reverse of ligand influx into the protein under unbiased MD, these findings challenge current views that preferential gates for recruitment of Fe²⁺ are 3‐fold channels with human ferritin, or the short path from the ferroxidase center to H93 with bacterial ferritins.     

Unveiling the Pathways of Dioxygen Through the C2 Component of the Environmentally Relevant Monooxygenase p‐Hydroxyphenylacetate Hydroxylase from Acinetobacter baumannii: A Molecular Dynamics Investigation

In this work, models of the homotetrameric C2 component of the monooxygenase p‐hydroxyphenylacetate hydroxylase from Acinetobacter baumannii, in complex with dioxygen (O₂) and, or not, the substrate p‐hydroxyphenylacetate (HPA) were built. Both models proved to be amenable to random‐acceleration molecular dynamics (RAMD) simulations, whereby a tiny randomly oriented external force, acting on O₂ at the active site in front of flavin mononucleotide (FMNH^?), accelerated displacement of O₂ toward the bulk solvent. This allowed us to carry out a sufficiently large number of RAMD simulations to be of statistical significance. The two systems behaved very similarly under RAMD, except for O₂ leaving the active site more easily in the absence of HPA, but then finding similar obstacles in getting to the gate as when the active site was sheltered by HPA. This challenges previous conclusions that HPA can only reach the active center after that the C4aOOH derivative of FMNH^? is formed, requiring uptake of O₂ at the active site before HPA. According to these RAMD simulations, O₂ could well get to FMNH^? also in the presence of the substrate at the active site.     

Structural similarity and functional diversity in diiron-oxo proteins

 Diiron-oxo proteins currently represent one of the most rapidly developing areas of bioinorganic chemistry. All of these proteins contain a four-helix bundle protein fold surrounding a (μ-carboxylato)diiron core, and most, if not all, of the diiron(II) sites appear to react with O₂ as part of their functional processes. Despite these common characteristics, an emerging functional diversity is one of the most striking aspects of this class of proteins. X-ray crystal structures of diiron(II) sites are now available for four of these proteins: hemerythrin (Hr), the hydroxylase protein of methane monooxygenase (MMOH), the R2 protein of Escherichia coli ribonucleotide reductase (RNR-R2), and a plant acyl-carrier protein Δ⁹-desaturase. The structure of the diiron(II) site in Hr, the sole O₂ carrier in the group, is clearly distinct from the other three, whose function is oxygen activation. The Hr diiron site is more histidine rich, and the oxygen-activating diiron sites contain a pair of (D/E)X_30–37EX₂H ligand sequence motifs, which is clearly not found in Hr. The Hr diiron site apparently permits only terminal O₂ coordination to a single iron, whereas the oxygen-activating diiron(II) centers present open or labile coordination sites on both irons of the center, and show a much greater coordinative flexibility upon oxidation to the diiron(III) state. Intermediates at the formal Fe^IIIFe^III and Fe^IVFe^IV oxidation levels for MMOH and formal Fe^IIIFe^IV oxidation level for RNR-R2 have been identified during reactions of the diiron(II) sites with O₂. An [Fe₂(μ-O)₂]^4+, 3+ "diamond core" structure has been proposed for the latter two oxidation levels. The intermediate at the Fe^IIIFe^IV oxidation level in RNR-R2 is kinetically competent to generate a stable, functionally essential tyrosyl radical. The Fe^IVFe^IV oxidation level is presumed to effect hydroxylation of hydrocarbons in MMOH, but the mechanism of this hydroxylation, particularly the involvement of discrete radicals, is currently controversial. The biological function of diiron sites in three members of this class, rubrerythrin, ferritin and bacterioferritin, remains enigmatic. Received: 31 July 1996 / Accepted: 4 October 1996     

On Dioxygen Permeation of MaL Laccase from the Thermophilic Fungus Melanocarpus albomyces: An all‐Atom Molecular Dynamics Investigation

In this article, biased molecular dynamics (MD) simulations of O₂ egress from the active center of MaL laccase toward the bulk solvent were described. Parameterization of the set of four Cu(II) ions, in the framework of CHARMM‐36 FF, was carried out on a recent dummy‐atom model that takes into account the Jahn–Teller effect. By carrying out a number of statistically relevant MD simulations, under a tiny randomly oriented external force applied to the molecule O₂, three preferred gates for O₂ egress from the enzyme and a few intermediate binding pockets (BPs) for O₂ were visible; all the gates and pockets were located on two of the three domains. This wide distribution of preferred gates notwithstanding the molecule O₂ was seen to follow specific pathways, exploiting consistently the interstices created by the enzyme thermal fluctuations. These are features that can be imagined to have evolved to make MaL laccase extremely efficient as catalysts in various reactions that require O₂.     

On Dioxygen Permeation through a Dehydrogenase‐Pyrroloquinoline Quinone Complex. A Molecular‐Dynamics Investigation

In this work, an all atom model of the quinoprotein dehydrogenase PqqC in complex with the PQQ (=4,5‐dihydro‐4,5‐dioxo‐1H‐pyrrolo[2,3‐f]quinoline‐2,7,9‐tricarboxylic acid) cofactor and dioxygen (O₂), solvated with TIP3 water in periodic boxes, was subjected to random‐acceleration molecular dynamics (RAMD). It was found that O₂ leaves the active binding pocket, in front of PQQ, to get to the solvent, as easily as with a variety of other O₂‐activating enzymes, O₂ carriers, and gas‐sensing proteins. The shortest pathway, orthogonal to the center of the mean plane of PQQ, was largely preferred by O₂ over pathways slightly deviating from this line. These observations challenge the interpretation of an impermeable active binding pocket of PqqC‐PQQ, as drawn from both X‐ray diffraction data of the crystal at low temperature and physiological experimentation.     

On the Permeation by Dioxygen of the Cofactor‐Independent Unusual Oxygenase RhCC,in Complex with Substrate 4‐Hydroxyphenylenolpyruvate. A Molecular Dynamics Investigation

This work deals with a trimeric bacterial protein, RhCC, which, although belonging to the tautomerase superfamily, shows oxygenase activity. A model of the complex from RhCC and substrate 4‐hydroxyphenylenolpyruvate (4HPP), fitting the observation of extra electron densities from X‐ray diffraction of the crystal, could be built by autodocking. When subjected to molecular dynamics (MD) aided by an external random force applied to a O₂ molecule placed above 4HPP, this model evolved with O₂ egressing toward the bulk solvent from two nearly opposite gates. These were located between the nearly parallel helices 75 – 91 and 15 – 33 of either chain C (gate SE) or chain B (gate FL). Alternatively, with four O₂ molecules in the bulk solvent, unbiased MD led to O₂ entering the protein from gate SE and getting to 4HPP, while forming a stabilizing salt bridge between the 4HPP carboxylate and P1.C ⁺NH₂, thus providing scientific ground for a refined model of the complex.     

On the Permeation by Dioxygen of Urate Oxidase from Aspergillus flavus in Complex with Xanthine Anion: Dioxygen Pathways and a Portrait of the Enzyme Cavities from Molecular Dynamics Simulations in Water Solution

This work describes molecular dynamics (MD) simulations in aqueous media for the complex of the homotetrameric urate oxidase (UOX) from Aspergillus flavus with xanthine anion ( 5 ) in the presence of dioxygen (O₂). After 196.6 ns of trajectory from unrestrained MD, a O₂ molecule was observed leaving the bulk solvent to penetrate the enzyme between two subunits, A/C. From here, the same O₂ molecule was observed migrating, across subunit C, to the hydrophobic cavity that shares residue V227 with the active site. The latter was finally attained, after 378.3 ns of trajectory, with O₂ at a bonding distance from 5 . The reverse same O₂ pathway, from 5 to the bulk solvent, was observed as preferred pathway under random acceleration MD (RAMD), where an external, randomly oriented force was acting on O₂. Both MD and RAMD simulations revealed several cavities populated by O₂ during its migration from the bulk solvent to the active site or backwards. Paying attention to the last hydrophobic cavity that apparently serves as O₂ reservoir for the active site, it was noticed that its volume undergoes ample fluctuations during the MD simulation, as expected from the thermal motion of a flexible protein, independently from the particular subunit and no matter whether the cavity is filled or not by O₂.     

Effects of elevated ozone concentration on methane emission from a rice paddy in Yangtze River Delta,China

Few investigations have been made on the impact of elevated ozone (O₃) concentration on methane (CH₄) emission from rice paddies. Using open‐top chambers in situ with different O₃ treatments, CH₄ emissions were measured in a rice paddy in Yangtze River Delta, China in 2007 and 2008. There were four treatments applied: charcoal‐filtered air (CF), nonfiltered air (NF), and charcoal‐filtered air with different O₃ additions (O₃‐1 and O₃‐2). The mean O₃ concentrations during the O₃ fumigation were 19.7, 22.6, 69.6 and 118.6 ppb in 2007 and 7.0, 17.4, 82.2 and 138.3 ppb in 2008 for treatments CF, NF, O₃‐1 and O₃‐2, respectively. The rice yields, as compared with CF, were reduced by 32.8% and 37.1%, 58.3% and 52.1% in treatments O₃‐1 and O₃‐2 in 2007 and 2008, respectively. The diurnal patterns of CH₄ emission varied temporally with treatments and there was inconsistence in diurnal variations in CH₄ emissions from the paddy field. The daily mean CH₄ emissions were significantly lower in treatments O₃‐1 and O₃‐2 than those in treatments CF and NF. Compared with CF treatment, CH₄ emissions from the paddy field were decreased to 46.5% and 38.3%, 50.6% and 46.8% under treatments O₃‐1 and O₃‐2 in the whole growing seasons of 2007 and 2008, respectively. The seasonal mean CH₄ emissions were negatively related with AOT40 (accumulative O₃ concentration above 40 ppb; P < 0.01 in both years), but positively related to the relative rice yield (reference to CF; P < 0.01 in 2007 and P < 0.001 in 2008), aboveground biomass (P < 0.01 in both years) and underground biomass (P < 0.01 in 2007 and P < 0.05 in 2008). The decreased CH₄ emission from the rice paddy due to an increased O₃ exposure might partially mitigate the global warming potential induced by soil carbon loss under elevated O₃ concentrations.     

Gates and Binding Pockets for Nitric Oxide with Cytochrome c′, According to Molecular Dynamics

Random‐acceleration molecular‐dynamics (RAMD) simulations with models of homodimeric 6‐ligated distal‐NO and 5‐ligated proximal‐NO cytochrome c′ complexes, in TIP3 H₂O, showed two distinct, non‐intercommunicating worlds. In the framework of a long cavity formed by four protein helices with heme at one extremity, NO was observed to follow different pathways with the two complexes to reach the solvent. With the 6‐ligated complex, NO was observed to progress by exploiting protein internal channels created by thermal fluctuations, and be temporarily trapped into binding pockets before reaching the preferred gate at the heme end of the cavity. In contrast, with the 5‐ligated complex, NO was observed to surface the solvent‐exposed helix 7, up to a gate at the other extremity of the protein, only occasionally finding an earlier, direct way out toward the solvent. That only bulk NO gets involved in forming the 5‐ligated proximal‐NO complex is in agreement with previous experimental observations, while the occurrence of binding pockets suggests that also reservoir NO might play a role with the distal‐NO complex.     

Effects of elevated ozone concentration on CH4 and N2O emission from paddy soil under fully open‐air field conditions

We investigated the effects of elevated ozone concentration (E‐O₃) on CH₄ and N₂O emission from paddies with two rice cultivars: an inbred Indica cultivar Yangdao 6 (YD6) and a hybrid one II‐you 084 (IIY084), under fully open‐air field conditions in China. A mean 26.7% enhancement of ozone concentration above the ambient level (A‐O₃) significantly reduced CH₄ emission at tillering and flowering stages leading to a reduction of seasonal integral CH₄ emission by 29.6% on average across the two cultivars. The reduced CH₄ emission is associated with O₃‐induced reduction in the whole‐plant biomass (?13.2%), root biomass (?34.7%), and maximum tiller number (?10.3%), all of which curbed the carbon supply for belowground CH₄ production and its release from submerged soil to atmosphere. Although no significant difference was detected between the cultivars in the CH₄ emission response to E‐O₃, a larger decrease in CH₄ emission with IIY084 (?33.2%) than that with YD6 (?7.0%) was observed at tillering stage, which may be due to the larger reduction in tiller number in IIY084 by E‐O₃. Additionally, E‐O₃ reduced seasonal mean NO_x flux by 5.7% and 11.8% with IIY084 and YD6, respectively, but the effects were not significant statistically. We found that the relative response of CH₄ emission to E‐O₃ was not significantly different from those reported in open‐top chamber experiments. This study has thus confirmed that increasing ozone concentration would mitigate the global warming potential of CH₄ and suggested consideration of the feedback mechanism between ozone and its precursor emission into the projection of future ozone effects on terrestrial ecosystem.     

Oxygen effects on methane production and oxidation in humid tropical forest soils

We investigated the effects of oxygen (O₂) concentration on methane (CH₄) production and oxidation in two humid tropical forests that differ in long‐term, time‐averaged soil O₂ concentrations. We identified sources and sinks of CH₄ through the analysis of soil gas concentrations, surface emissions, and carbon isotope measurements. Isotope mass balance models were used to calculate the fraction of CH₄ oxidized in situ. Complementary laboratory experiments were conducted to determine the effects of O₂ concentration on gross and net rates of methanogenesis. Field and laboratory experiments indicated that high levels of CH₄ production occurred in soils that contained between 9±1.1% and 19±0.2% O₂. For example, we observed CH₄ concentrations in excess of 3% in soils with 9±1.1% O₂. CH₄ emissions from the lower O₂ sites were high (22–101 nmol CH₄ m^?2 s^?1), and were equal in magnitude to CH₄ emissions from natural wetlands. During peak periods of CH₄ efflux, carbon dioxide (CO₂) emissions became enriched in ¹³C because of high methanogenic activity. Gross CH₄ production was probably greater than flux measurements indicated, as isotope mass balance calculations suggested that 48–78% of the CH₄ produced was oxidized prior to atmospheric egress. O₂ availability influenced CH₄ oxidation more strongly than methanogenesis. Gross CH₄ production was relatively insensitive to O₂ concentrations in laboratory experiments. In contrast, methanotrophic bacteria oxidized a greater fraction of total CH₄ production with increasing O₂ concentration, shifting the δ¹³C composition of CH₄ to values that were more positive. Isotopic measurements suggested that CO₂ was an important source of carbon for methanogenesis in humid forests. The δ¹³C value of methanogenesis was between ?84‰ and ?98‰, which is well within the range of CH₄ produced from CO₂ reduction, and considerably more depleted in ¹³C than CH₄ formed from acetate.     

Molecular‐Dynamics Simulation of Dioxygen Egress from 12/15‐Lipoxygenase?Arachidonic Acid Complex

Extensive random‐acceleration molecular‐dynamics (RAMD) simulations of the egress of dioxygen (O₂) from a model of rabbit 12/15‐lipoxygenase? arachidonic acid complex disclosed several exit portals in addition to those previously described from implicit ligand sampling calculations and limited MD simulations.     

From Dioxygen Storing to Dioxygen Sensing with Neuroglobins: An Insight from Molecular Mechanics

This work deals with two neuroglobins from phylogenetically distant organisms. Deriving from the acoelomorph Symsagittifera roscoffensis, SrNgb is functionally pentacoordinated, and is assumed to function as a reserve of dioxygen (O₂). Obtained from mice, mNgb is functionally hexacoordinated, and presumably triggers signals from sensing O₂. Here, it is investigated how these two globins are permeated by diatomic gases, SrNgb by O₂ and mNgb by CO. With protein atomic coordinates available from high‐resolution X‐ray diffraction analysis, O₂ and CO pathways were traced from molecular‐dynamics simulations in H₂O solution, which makes no difference between the two gases, accelerated by applying an external randomly‐oriented minimal force to the center of mass of the diatomic gas molecule. This allowed us to explore a statistically significant large number of trajectories. It emerged that CO leaves mNgb from preferentially peripheral gates located on the side of the heme propionate chains, whereas O₂ leaves SrNgb from the opposite side. This shows no analogy with either the functionally pentacoordinated, O₂‐transporting, myoglobin (Mgb), or the hexacoordinated, O₂‐sensing, cytoglobin, despite the same three‐over‐three typical α‐helical globin folding. The sole analogy that could be observed was a preference for the shortest diatomic gas pathways with both SrNgb and Mgb. It is tempting to speculate that this fulfills the need of being quick in delivering O₂ to depleted organs.     

X-ray absorption spectroscopic studies of the methane monooxygenase hydroxylase component from Methylosinus trichosporium OB3b

The conversion from methane to methanol is catalyzed by methane monooxygenase (MMO) in methanotrophic bacteria. Earlier work on the crystal structures of the MMO hydroxylase component (MMOH) from Methylococcus capsulatus (Bath) at 4??°C and –160??°C has revealed two different core arrangements for the diiron active site. To ascertain the generality of these results, we have now carried out the first structural characterization on MMOH from Methylosinus trichosporium OB3b. Our X-ray absorption spectroscopic (XAS) analysis suggests the presence of two Fe-Fe distances of about 3?Å and 3.4?Å, which are proposed to reflect two populations of MMOH molecules with either a bis(μ-hydroxo)(μ-carboxylato)- or a (μ-hydroxo)(μ-carboxylato)diiron(III) core structure, respectively. The observation of these two different core structures, together with the crystallographic results of the MMOH from Methylococcus capsulatus (Bath), suggests the presence of an equilibrium that may reflect a core flexibility that is required to accommodate the various intermediates in the catalytic cycle of the enzyme. XAS studies on the binding of component B (MMOB) to the hydroxylase component show that MMOB does not perturb either this equilibrium or the gross structure of the oxidized diiron site in MMOH.     

Crystal structure of the secreted protein HP1454 from the human pathogen Helicobacter pylori

HP1454 is a protein of 303 amino acids found in the extracellular milieu of Helicobacter pylori. The protein structure, crystallized in the orthorhombic C222₁ space group with one molecule per asymmetric unit, has been determined using the single‐wavelength anomalous dispersion method. HP1454 exhibits an elongated bent shape, composed of three distinct domains. Each domain possesses a fold already present in other structures: Domain I contains a three‐strand antiparallel β‐barrel flanked by a long α‐helix, Domain II is an anti‐parallel three‐helix bundle, and Domain III a β‐sheet flanked by two α‐helices. The overall assembly of the protein does not bear any similarity with known structures. Proteins 2014; 82:2868–2873. © 2014 Wiley Periodicals, Inc.     

Ligand Entry and Exit Pathways in the β2-Adrenergic Receptor

The recently determined crystal structure of the human β₂-adrenergic (β₂AR) G-protein-coupled receptor provides an excellent structural basis for exploring β₂AR-ligand binding and dissociation process. Based on this crystal structure, we simulated ligand exit from the β₂AR receptor by applying the random acceleration molecular dynamics (RAMD) simulation method. The simulation results showed that the extracellular opening on the receptor surface was the most frequently observed egress point (referred to as pathway A), and a few other pathways through interhelical clefts were also observed with significantly lower frequencies. In the egress trajectories along pathway A, the D192-K305 salt bridge between the extracellular loop 2 (ECL2) and the apex of the transmembrane helix 7 (TM7) was exclusively broken. The spatial occupancy maps of the ligand computed from the 100 RAMD simulation trajectories indicated that the receptor-ligand interactions that restrained the ligand in the binding pocket were the major resistance encountered by the ligand during exit and no second barrier was notable. We next performed RAMD simulations by using a putative ligand-free conformation of the receptor as input structure. This conformation was obtained in a standard molecular dynamics simulation in the absence of the ligand and it differed from the ligand-bound conformation in a hydrophobic patch bridging ECL2 and TM7 due to the rotation of F193 of ECL2. Results from the RAMD simulations with this putative ligand-free conformation suggest that the cleft formed by the hydrophobic bridge, TM2, TM3, and TM7 on the extracellular surface likely serves as a more specific ligand-entry site and the ECL2-TM7 hydrophobic junction can be partially interrupted upon the entry of ligand that pushes F193 to rotate, resulting in a conformation as observed in the ligand-bound crystal structure. These results may help in the design of β₂AR-targeting drugs with improved efficacy, as well as in understanding the receptor subtype selectivity of ligand binding in the β family of the adrenergic receptors that share almost identical ligand-binding pockets, but show notable amino acid sequence divergence in the putative ligand-entry site, including ECL2 and the extracellular end of TM7.     

On the Dynamical Behavior of the Cysteine Dioxygenase‐l‐Cysteine Complex in the Presence of Free Dioxygen and l‐Cysteine

In this work, viable models of cysteine dioxygenase (CDO) and its complex with l ‐cysteine dianion were built for the first time, under strict adherence to the crystal structure from X‐ray diffraction studies, for all atom molecular dynamics (MD). Based on the CHARMM36 FF, the active site, featuring an octahedral dummy Fe(II) model, allowed us observing water exchange, which would have escaped attention with the more popular bonded models. Free dioxygen (O₂) and l ‐cysteine, added at the active site, could be observed being expelled toward the solvating medium under Random Accelerated Molecular Dynamics (RAMD) along major and minor pathways. Correspondingly, free dioxygen (O₂), added to the solvating medium, could be observed to follow the same above pathways in getting to the active site under unbiased MD. For the bulky l ‐cysteine, 600 ns of trajectory were insufficient for protein penetration, and the molecule was stuck at the protein borders. These models pave the way to free energy studies of ligand associations, devised to better clarify how this cardinal enzyme behaves in human metabolism.