The Cluster of Hydrophobic Residues Controls the Entrance to the Active Site of Choline Oxidase

The in-silico study of the structural changes in the Arthrobacter globiformis choline oxidase induced by high temperature

BACKGROUND

Choline oxidase, a flavoprotein, is an enzyme that catalyzes the reaction which converts choline into glycine betaine. Choline oxidase started its journey way back in 1933. However, the impact of the high temperature on its structure has not been explored despite the long history and availability of its crystal structure. Both choline oxidase and its product, glycine betaine, have enormous applications spanning across multiple industries. Understanding how the 3D structure of the enzyme will change with the temperature change can open new ways to make it more stable and useful for industry.

PROCESS

This research paper presents the in-silico study and analysis of the structural changes of A. globiformis choline oxidase at temperatures from 25 °C to 60 °C. A step-wise process is depicted in Fig. 1.

RESULTS

Multiple sequence alignment (MSA) of 11 choline oxidase sequences from different bacteria vs Arthrobacter globiformis choline oxidase showed that active site residues are highly conserved. The available crystal structure of A. globiformis choline oxidase with cofactor Flavin Adenine Dinucleotide (FAD) in the dimeric state (PDB ID: 4MJW)1 was considered for molecular dynamics simulations. A simulated annealing option was used to gradually increase the temperature of the system from 25 °C to 60 °C. Analysis of the conserved residues, as well as residues involved in Flavin Adenine Dinucleotide (FAD) binding, substrate binding, substate gating, and dimer formationwas done. At high temperatures, the formation of the inter-chain salt bridge between Arg50 and Glu63 was a significant observation near the active site of choline oxidase.

CONCLUSION

Molecular dynamics studies suggest that an increase in temperature has a significant impact on the extended Flavin Adenine Dinucleotide (FAD) binding region. These changes interfere with the entry of substrate to the active site of the enzyme and make the enzyme inactive.

Enzyme engineering of choline oxidase for improving stability

Functioning as a flavoprotein, choline oxidase facilitates the transformation of choline into glycine betaine. Notably, choline oxidase and its resultant product, glycine betaine, find extensive applications across various industries and fields of study. However, its high sensitivity and tendency to lose functional activity at high temperatures reduces its industrial usage. MD simulation and mutation studies have revealed the role of certain residues responsible for the enzyme's thermal instability. This study focuses on inducing thermal stability to choline oxidase of A. globiformis through computational approaches at a maximum temperature of 60 °C. MD simulation analysis showed that Trp 331, Val 464 and Ser 101 contribute to structural instability, leading to the instability at 60 °C. Mutation of these residues with phenylalanine residues and simulation of the mutated enzyme at 60 °C exhibited thermostability and insignificant residual fluctuation. The re-docking and MM/GBSA analyses further validated the mutated enzyme's binding affinity and catalytic activity.Communicated by Ramaswamy H. Sarma.

Simplifying the complexity: Single enzyme (choline oxidase) inhibition-based biosensor with dual-readout method for organophosphorus pesticide detection

Organophosphorus pesticides (OPs) are widely used in agricultural production, but their residues could cause pollution to the environment and living organisms. In this paper, a simple dual-readout method for OPs detection was proposed based on ChOx single enzyme inhibition. Firstly, ChOx can catalyze the production of H₂O₂ from choline chloride (Ch-Cl). Bifunctional iron-doped carbon dots (Fe-CDs) with good peroxidase-like activity and superior fluorescence properties can catalyze the oxidation of 3,3',5,5'-tetramethylbenzidine (TMB) to blue oxidized TMB (oxTMB) by H₂O₂ formed, and oxTMB could quench the fluorescence of Fe-CDs. In light of the fact that OPs exhibited activity in inhibiting ChOx, less H₂O₂ and the decreasing oxTMB led to a result that the fluorescence of the system recovered and the solution became lighter in blue color. Moreover, the process of ChOx inhibition by OPs was analyzed by molecular docking technique and it was found that OPs interact with key amino acid residues catalyzed by ChOx (Asn510, His466, Ser101, His351, Phe357, Trp331, Glu312). Finally, a dual-mode (colorimetry and fluorescence) sensor was created for the detection of OPs with the detection limit of 6 ng/L, and was successfully used in the quantitative determination of OPs in actual samples with satisfactory results.

Reaction mechanisms and applications of aryl-alcohol oxidase

Aryl-alcohol oxidases (AAO) constitute a family of FAD-containing enzymes, included in the glucose-methanol-choline oxidase/dehydrogenase superfamily of proteins. They are commonly found in fungi, where their eco-physiological role is to produce hydrogen peroxide that activates ligninolytic peroxidases in white-rot (lignin-degrading) basidiomycetes or to trigger the Fenton reactions in brown-rot (carbohydrate-degrading) basidiomycetes. These enzymes catalyze the oxidation of a plethora of aromatic, and some aliphatic, polyunsaturated alcohols bearing conjugated primary hydroxyl group. Besides, the enzymes show activity on the hydrated forms of the corresponding aldehydes. Some AAO features, such as the broad range of substrates that it can oxidize (with the only need of molecular oxygen as co-substrate) and its stereoselective mechanism, confer good properties to these enzymes as industrial biocatalysts. In fact, AAO can be used for different biotechnological applications, such as flavor synthesis, secondary alcohol deracemization and oxidation of furfurals for the production of furandicarboxylic acid as a chemical building block. Also, AAO can participate in processes of interest in the wood biorefinery and textile industries as an auxiliary enzyme providing hydrogen peroxide to ligninolytic or dye-decolorizing peroxidases. Both rational design and directed molecular evolution have been employed to engineer AAO for some of the above biotechnological applications.

Choline oxidases

Choline oxidase catalyzes the four-electron, two-step, flavin-mediated oxidation of choline to glycine betaine. The enzyme is important both for medical and biotechnological reasons, because glycine betaine is one among a limited number of compatible solutes used by cells to counteract osmotic pressure. From a fundamental standpoint, choline oxidase has emerged as one of the paradigm enzymes for the oxidation of alcohols catalyzed by flavoproteins. Mechanistic, structural, and computational studies have elucidated the mechanism of action of the enzyme from Arthrobacter globiformis at the molecular level. Both choline and oxygen access to the active site cavity are gated and tightly controlled. Amino acid residues involved in substrate binding, and their contribution, have been identified. The mechanism of choline oxidation, with a hydride transfer reaction, an asynchronous transition state, the formation and stabilization of an alkoxide transient species, and a quantum mechanical mode of reaction, has been elucidated. The importance of nonpolar side chains for oxygen localization and of the positive charge harbored on the substrate for activation of oxygen for reaction with the reduced flavin have been recognized. Interesting phenomena, like the formation of a metastable photoinduced flavin-protein adduct, the reversible formation of a bicovalent flavoprotein, and the trapping of the enzyme in inactive conformations, have been described. This review summarizes the current status of our understanding on the structure-function-dynamics of choline oxidase.

Multiple implications of an active site phenylalanine in the catalysis of aryl-alcohol oxidase

Aryl-alcohol oxidase (AAO) has demonstrated to be an enzyme with a bright future ahead due to its biotechnological potential in deracemisation of chiral compounds, production of bioplastic precursors and other reactions of interest. Expanding our understanding on the AAO reaction mechanisms, through the investigation of its structure-function relationships, is crucial for its exploitation as an industrial biocatalyst. In this regard, previous computational studies suggested an active role for AAO Phe397 at the active-site entrance. This residue is located in a loop that partially covers the access to the cofactor forming a bottleneck together with two other aromatic residues. Kinetic and affinity spectroscopic studies, complemented with computational simulations using the recently developed adaptive-PELE technology, reveal that the Phe397 residue is important for product release and to help the substrates attain a catalytically relevant position within the active-site cavity. Moreover, removal of aromaticity at the 397 position impairs the oxygen-reduction activity of the enzyme. Experimental and computational findings agree very well in the timing of product release from AAO, and the simulations help to understand the experimental results. This highlights the potential of adaptive-PELE to provide answers to the questions raised by the empirical results in the study of enzyme mechanisms.

An integrated perspective on RNA aptamer ligand-recognition models: clearing muddy waters

Riboswitches are short RNA motifs that sensitively and selectively bind cognate ligands to modulate gene expression. Like protein receptor-ligand pairs, their binding dynamics are traditionally categorized as following one of two paradigmatic mechanisms: conformational selection and induced fit. In conformational selection, ligand binding stabilizes a particular state already present in the receptor's dynamic ensemble. In induced fit, ligand-receptor interactions enable the system to overcome the energetic barrier into a previously inaccessible state. In this article, we question whether a polarized division of RNA binding mechanisms truly meets the conceptual needs of the field. We will review the history behind this classification of RNA-ligand interactions, and the way induced fit in particular has been rehabilitated by single-molecule studies of RNA aptamers. We will highlight several recent results from single-molecule experimental studies of riboswitches that reveal gaps or even contradictions between common definitions of the two terms, and we will conclude by proposing a more robust framework that considers the range of RNA behaviors unveiled in recent years as a reality to be described, rather than an increasingly unwieldy set of exceptions to the traditional models.

Substitutions of S101 decrease proton and hydride transfers in the oxidation of betaine aldehyde by choline oxidase

Choline oxidase oxidizes choline to glycine betaine, with two flavin-mediated reactions to convert the alcohol substrate to the carbon acid product. Proton abstraction from choline or hydrated betaine aldehyde in the wild-type enzyme occurs in the mixing time of the stopped-flow spectrophotometer, thereby precluding a mechanistic investigation. Mutagenesis of S101 rendered the proton transfer reaction amenable to study. Here, we have investigated the aldehyde oxidation reaction catalyzed by the mutant enzymes using steady-state and rapid kinetics with betaine aldehyde. Stopped-flow traces for the reductive half-reaction of the S101T/V/C variants were biphasic, corresponding to the reactions of proton abstraction and hydride transfer. In contrast, the S101A enzyme yielded monophasic traces like wild-type choline oxidase. The rate constants for proton transfer in the S101T/C/V variants decreased logarithmically with increasing hydrophobicity of residue 101, indicating a behavior different from that seen previously with choline for which no correlation was determined. The rate constants for hydride transfer also showed a logarithmic decrease with increasing hydrophobicity at position 101, which was similar to previous results with choline as a substrate for the enzyme. Thus, the hydrophilic character of S101 is necessary not only for efficient hydride transfer but also for the proton abstraction reaction.

Evidence for proton tunneling and a transient covalent flavin-substrate adduct in choline oxidase S101A

The effect of temperature on the reaction of alcohol oxidation catalyzed by choline oxidase was investigated with the S101A variant of choline oxidase. Anaerobic enzyme reduction in a stopped-flow spectrophotometer was biphasic using either choline or 1,2-[²H₄]-choline as a substrate. The limiting rate constants k_lim1 and k_lim2 at saturating substrate were well separated (k_lim1/k_lim2>9), and were >15-fold slower than for wild-type choline oxidase. Solvent deuterium kinetic isotope effects (KIEs) ~4 established that k_lim1 probes the proton transfer from the substrate hydroxyl to a catalytic base. Primary substrate deuterium KIEs ≥7 demonstrated that k_lim2 reports on hydride transfer from the choline alkoxide to the flavin. Between 15°C and 39°C the k_lim1 and k_lim2 values increased with increasing temperature, allowing for the analyses of H⁺ and H^- transfers using Eyring and Arrhenius formalisms. Temperature-independent KIE on the k_lim1 value (^H2Ok_lim1/^D2Ok_lim1) suggests that proton transfer occurs within a highly reorganized tunneling-ready-state with a narrow distribution of donor-acceptor distances. Eyring analysis of the k_lim2 value gave lines with the slope_(choline)>slope_(D-choline), suggesting kinetic complexity. Spectral evidence for the transient occurrence of a covalent flavin-substrate adduct during the first phase of the anaerobic reaction of S101A CHO with choline is presented, supporting the notion that an important role of amino acid residues in the active site of flavin-dependent enzymes is to eliminate alternative reactions of the versatile enzyme-bound flavin for the reaction that needs to be catalyzed.

Enzyme Tunnels and Gates As Relevant Targets in Drug Design

Many enzymes contain tunnels and gates that are essential to their function. Gates reversibly switch between open and closed conformations and thereby control the traffic of small molecules-substrates, products, ions, and solvent molecules-into and out of the enzyme's structure via molecular tunnels. Many transient tunnels and gates undoubtedly remain to be identified, and their functional roles and utility as potential drug targets have received comparatively little attention. Here, we describe a set of general concepts relating to the structural properties, function, and classification of these interesting structural features. In addition, we highlight the potential of enzyme tunnels and gates as targets for the binding of small molecules. The different types of binding that are possible and the potential pharmacological benefits of such targeting are discussed. Twelve examples of ligands bound to the tunnels and/or gates of clinically relevant enzymes are used to illustrate the different binding modes and to explain some new strategies for drug design. Such strategies could potentially help to overcome some of the problems facing medicinal chemists and lead to the discovery of more effective drugs.

Protein Ensembles: How Does Nature Harness Thermodynamic Fluctuations for Life? The Diverse Functional Roles of Conformational Ensembles in the Cell

All soluble proteins populate conformational ensembles that together constitute the native state. Their fluctuations in water are intrinsic thermodynamic phenomena, and the distributions of the states on the energy landscape are determined by statistical thermodynamics; however, they are optimized to perform their biological functions. In this review we briefly describe advances in free energy landscape studies of protein conformational ensembles. Experimental (nuclear magnetic resonance, small-angle X-ray scattering, single-molecule spectroscopy, and cryo-electron microscopy) and computational (replica-exchange molecular dynamics, metadynamics, and Markov state models) approaches have made great progress in recent years. These address the challenging characterization of the highly flexible and heterogeneous protein ensembles. We focus on structural aspects of protein conformational distributions, from collective motions of single- and multi-domain proteins, intrinsically disordered proteins, to multiprotein complexes. Importantly, we highlight recent studies that illustrate functional adjustment of protein conformational ensembles in the crowded cellular environment. We center on the role of the ensemble in recognition of small- and macro-molecules (protein and RNA/DNA) and emphasize emerging concepts of protein dynamics in enzyme catalysis. Overall, protein ensembles link fundamental physicochemical principles and protein behavior and the cellular network and its regulation.

Role of F357 as an Oxygen Gate in the Oxidative Half-Reaction of Choline Oxidase

Choline oxidase from Arthrobacter globiformis catalyzes the oxidation of choline to glycine betaine by using oxygen as an electron acceptor. A partially rate limiting isomerization of the reduced wild-type enzyme during the reaction with oxygen was previously detected using solvent viscosity effects. In this study, we hypothesized that the side chains of M62 and F357, located at the entrance to the active site of choline oxidase, may be related to the slow isomerization detected. We engineered a double-variant enzyme M62A/F357A. The kinetic characterization of the double-variant enzyme showed a lack of the isomerization detected in wild-type choline oxidase, and a lack of saturation with an oxygen concentration as high as 1 mM, while most other kinetic parameters were similar to those of wild-type choline oxidase. The kinetic characterization of the single-variant enzymes established that only the side chain of F357 plays a role in the isomerization of choline oxidase in the oxidative half-reaction. Molecular dynamics studies suggest that the slow isomerization related to F357 is possibly due to the participation of the phenyl ring in a newly proposed gating mechanism for a narrow tunnel, assumed to regulate the access of oxygen to the reduced cofactor.

Structure of choline oxidase in complex with the reaction product glycine betaine

Choline oxidase fromArthrobacter globiformis, which is involved in the biosynthesis of glycine betaine from choline, has been extensively characterized in its mechanistic and structural properties. Despite the knowledge gained on the enzyme, the details of substrate access to the active site are not fully understood. The `loop-and-lid' mechanism described for the glucose–methanol–choline enzyme superfamily has not been confirmed for choline oxidase. Instead, a hydrophobic cluster on the solvent-accessible surface of the enzyme has been proposed by molecular dynamics to control substrate access to the active site. Here, the crystal structure of the enzyme was solved in complex with glycine betaine at pH 6.0 at 1.95 Å resolution, allowing a structural description of the ligand–enzyme interactions in the active site. This structure is the first of choline oxidase in complex with a physiologically relevant ligand. The protein structures with and without ligand are virtually identical, with the exception of a loop at the dimer interface, which assumes two distinct conformations. The different conformations of loop 250–255 define different accessibilities of the proposed active-site entrance delimited by the hydrophobic cluster on the other subunit of the dimer, suggesting a role in regulating substrate access to the active site.

Human choline dehydrogenase: medical promises and biochemical challenges

Human choline dehydrogenase (CHD) is located in the inner membrane of mitochondria primarily in liver and kidney and catalyzes the oxidation of choline to glycine betaine. Its physiological role is to regulate the concentrations of choline and glycine betaine in the blood and cells. Choline is important for regulation of gene expression, the biosynthesis of lipoproteins and membrane phospholipids and for the biosynthesis of the neurotransmitter acetylcholine; glycine betaine plays important roles as a primary intracellular osmoprotectant and as methyl donor for the biosynthesis of methionine from homocysteine, a required step for the synthesis of the ubiquitous methyl donor S-adenosyl methionine. Recently, CHD has generated considerable medical attention due to its association with various human pathologies, including male infertility, homocysteinuria, breast cancer and metabolic syndrome. Despite the renewed interest, the biochemical characterization of the enzyme has lagged behind due to difficulties in the obtainment of purified, active and stable enzyme. This review article summarizes the medical relevance and the physiological roles of human CHD, highlights the biochemical knowledge on the enzyme, and provides an analysis based on the comparison of the protein sequence with that of bacterial choline oxidase, for which structural and biochemical information is available.

Modulating O2 reactivity in a fungal flavoenzyme: involvement of aryl-alcohol oxidase Phe-501 contiguous to catalytic histidine

Aryl-alcohol oxidase (AAO) is a flavoenzyme responsible for activation of O(2) to H(2)O(2) in fungal degradation of lignin. The AAO crystal structure shows a buried active site connected to the solvent by a hydrophobic funnel-shaped channel, with Phe-501 and two other aromatic residues forming a narrow bottleneck that prevents the direct access of alcohol substrates. However, ligand diffusion simulations show O(2) access to the active site following this channel. Site-directed mutagenesis of Phe-501 yielded a F501A variant with strongly reduced O(2) reactivity. However, a variant with increased reactivity, as shown by kinetic constants and steady-state oxidation degree, was obtained by substitution of Phe-501 with tryptophan. The high oxygen catalytic efficiency of F501W, ∼2-fold that of native AAO and ∼120-fold that of F501A, seems related to a higher O(2) availability because the turnover number was slightly decreased with respect to the native enzyme. Free diffusion simulations of O(2) inside the active-site cavity of AAO (and several in silico Phe-501 variants) yielded >60% O(2) population at 3-4 Å from flavin C4a in F501W compared with 44% in AAO and only 14% in F501A. Paradoxically, the O(2) reactivity of AAO decreased when the access channel was enlarged and increased when it was constricted by introducing a tryptophan residue. This is because the side chain of Phe-501, contiguous to the catalytic histidine (His-502 in AAO), helps to position O(2) at an adequate distance from flavin C4a (and His-502 Nε). Phe-501 substitution with a bulkier tryptophan residue resulted in an increase in the O(2) reactivity of this flavoenzyme.

Substrate diffusion and oxidation in GMC oxidoreductases: an experimental and computational study on fungal aryl-alcohol oxidase

AAO (aryl-alcohol oxidase) provides H2O2 in fungal degradation of lignin, a process of high biotechnological interest. The crystal structure of AAO does not show open access to the active site, where different aromatic alcohols are oxidized. In the present study we investigated substrate diffusion and oxidation in AAO compared with the structurally related CHO (choline oxidase). Cavity finder and ligand diffusion simulations indicate the substrate-entrance channel, requiring side-chain displacements and involving a stacking interaction with Tyr92. Mixed QM (quantum mechanics)/MM (molecular mechanics) studies combined with site-directed mutagenesis showed two active-site catalytic histidine residues, whose substitution strongly decreased both catalytic and transient-state reduction constants for p-anisyl alcohol in the H502A (over 1800-fold) and H546A (over 35-fold) variants. Combination of QM/MM energy profiles, protonation predictors, molecular dynamics, mutagenesis and pH profiles provide a robust answer regarding the nature of the catalytic base. The histidine residue in front of the FAD ring, AAO His502 (and CHO His466), acts as a base. For the two substrates assayed, it was shown that proton transfer preceded hydride transfer, although both processes are highly coupled. No stable intermediate was observed in the energy profiles, in contrast with that observed for CHO. QM/MM, together with solvent KIE (kinetic isotope effect) results, suggest a non-synchronous concerted mechanism for alcohol oxidation by AAO.

Enzyme dynamics point to stepwise conformational selection in catalysis

Recent data increasingly reveal that conformational dynamics are indispensable to enzyme function throughout the catalytic cycle, in substrate recruiting, chemical transformation, and product release. Conformational transitions may involve conformational selection and induced fit, which can be viewed as a special case in the catalytic network. NMR, X-ray crystallography, single-molecule FRET, and simulations clearly demonstrate that the free enzyme dynamics already encompass all the conformations necessary for substrate binding, preorganization, transition-state stabilization, and product release. Conformational selection and substate population shift at each step of the catalytic turnover can accommodate enzyme specificity and efficiency. Within such a framework, entropy can have a larger role in conformational dynamics than in direct energy transfer in dynamically promoted catalysis.