Mehdi D. Davari

Photophysics of the LOV-Based Fluorescent Protein Variant iLOV-Q489K Determined by Simulation and Experiment

Light, oxygen, voltage (LOV) based fluorescent proteins (FPs) represent a promising alternative to fluorescent reporters of the green fluorescent protein family. For certain applications like multicolor imaging or the design of FRET-based biosensors, the generation of spectrally shifted LOV-based FPs would be required. In a recent theoretical study ( Khrenova J. Phys. Chem. B 2015 , 119 ( 16 ), pp 5176 - 5183 ), the photophysical properties of a variant of the LOV-based fluorescent protein iLOV were predicted using quantum mechanics/molecular mechanics (QM/MM) approaches. The variant contained a lysine residue at the position of a highly conserved glutamine residue (Q489K), which directly interacts with the O4 and N5 atom of the flavin mononucleotide (FMN) chromophore. On the basis of QM/MM calculations, iLOV-Q489K was suggested to possess substantially red-shifted absorption and fluorescence-emission maxima with respect to parental iLOV. Here, we describe the experimental characterization of this variant, which, surprisingly contrary to the theoretical prediction, shows blue-shifted absorption and fluorescence-emission maxima. Using molecular dynamics (MD) simulations and QM/MM calculations, the molecular basis for the contradictory theoretical and experimental results is presented. Essentially, our computational analysis suggests that, in the Q489K variant, two possible side-chain conformers exist: (i) a least populated conformer K489in forming a hydrogen bond with the O4 atom of FMN chromophore and (ii) a most populated conformer K489out with the side-chain amino group flipped away from the FMN chromophore forming a new hydrogen bond with the backbone oxygen of G487. QM/MM calculated spectra of the K489out conformer are blue-shifted compared to the calculated spectra of parental iLOV, which is in accordance with experimental data. This suggests that the change in the conformation of K489 from K498in to K489out accounts for the change in the direction of the spectral shift from red to blue, thus reconciling theory and experiment.

Investigation of Structural Determinants for the Substrate Specificity in the Zinc-Dependent Alcohol Dehydrogenase CPCR2 from Candida parapsilosis.

Zinc‐dependent alcohol dehydrogenases (ADHs) are a class of enzymes applied in different biocatalytic processes ranging from lab to industrial scale. However, one drawback is the limited substrate range, necessitating a whole array of different ADHs for the relevant substrate classes. In this study, we investigated structural determinants of the substrate spectrum in the zinc‐dependent ADH carbonyl reductase 2 from Candida parapsilosis (CPCR2), combining methods of mutational analysis with in silico substrate docking. Assigned active site residues were genetically randomized, and the resulting mutant libraries were screened with a selection of challenging carbonyl substrates. Three variants (C57A, W116K, and L119M) with improved activities toward different substrates were detected at neighboring positions in the active site. Thus, all possible combinations of the mutations were generated and characterized for their substrate specificity, yielding several improved variants. The most interesting were a C57A variant, with a 27‐fold increase in specific activity for 4′‐acetamidoacetophenone, and the double mutant CPCR2 B16‐(C57A, L119M), with a 45‐fold improvement in the kcat⋅KM−1 value. The obtained variants were further investigated by in silico docking experiments. The results indicate that the mentioned residues are structural determinants of the substrate specificity of CPCR2, being major players in the definition of the active site. Comparison of these results with closely related enzymes suggests that these might even be transferred to other ADHs.

Activity prediction of substrates in NADH-dependent carbonyl reductase by docking requires catalytic constraints and charge parameterization of catalytic zinc environment

Molecular docking of substrates is more challenging compared to inhibitors as the reaction mechanism has to be considered. This becomes more pronounced for zinc-dependent enzymes since the coordination state of the catalytic zinc ion is of greater importance. In order to develop a predictive substrate docking protocol, we have performed molecular docking studies of diketone substrates using the catalytic state of carbonyl reductase 2 from Candida parapsilosis (CPCR2). Different docking protocols using two docking methods (AutoDock Vina and AutoDock4.2) with two different sets of atomic charges (AM1-BCC and HF-RESP) for catalytic zinc environment and substrates as well as two sets of vdW parameters for zinc ion were examined. We have selected the catalytic binding pose of each substrate by applying mechanism based distance criteria. To compare the performance of the docking protocols, the correlation plots for the binding energies of these catalytic poses were obtained against experimental Vmax values of the 11 diketone substrates for CPCR2. The best correlation of 0.73 was achieved with AutoDock4.2 while treating catalytic zinc ion in optimized non-bonded (NBopt) state with +1.01 charge on the zinc ion, compared to 0.36 in non-bonded (+2.00 charge on the zinc ion) state. These results indicate the importance of catalytic constraints and charge parameterization of catalytic zinc environment for the prediction of substrate activity in zinc-dependent enzymes by molecular docking. The developed predictive docking protocol described here is in principle generally applicable for the efficient in silico substrate spectra characterization of zinc-dependent ADH.Graphical Abstract

What’s My Substrate? Computational Function Assignment of Candida parapsilosis ADH5 by Genome Database Search, Virtual Screening, and QM/MM Calculations

Zinc-dependent medium chain reductase from Candida parapsilosis can be used in the reduction of carbonyl compounds to pharmacologically important chiral secondary alcohols. To date, the nomenclature of cpADH5 is differing (CPCR2/RCR/SADH) in the literature, and its natural substrate is not known. In this study, we utilized a substrate docking based virtual screening method combined with KEGG, MetaCyc pathway, and Candida genome databases search for the discovery of natural substrates of cpADH5. The virtual screening of 7834 carbonyl compounds from the ZINC database provided 94 aldehydes or methyl/ethyl ketones as putative carbonyl substrates. Out of which, 52 carbonyl substrates of cpADH5 with catalytically active docking pose were identified by employing mechanism based substrate docking protocol. Comparison of the virtual screening results with KEGG, MetaCyc database search, and Candida genome pathway analysis suggest that cpADH5 might be involved in the Ehrlich pathway (reduction of fusel aldehydes in leucine, isoleucine, and valine degradation). Our QM/MM calculations and experimental activity measurements affirmed that butyraldehyde substrates are the potential natural substrates of cpADH5, suggesting a carbonyl reductase role for this enzyme in butyraldehyde reduction in aliphatic amino acid degradation pathways. Phylogenetic tree analysis of known ADHs from Candida albicans shows that cpADH5 is close to caADH5. We therefore propose, according to the experimental substrate identification and sequence similarity, the common name butyraldehyde dehydrogenase cpADH5 for Candida parapsilosis CPCR2/RCR/SADH.

Disulfide-induced self-assembled targets: A novel strategy for the label free colorimetric detection of DNAs/RNAs via unmodified gold nanoparticles

A modified non-cross-linking gold-nanoparticles (Au-NPs) aggregation strategy has been developed for the label free colorimetric detection of DNAs/RNAs based on self-assembling target species in the presence of thiolated probes. Two complementary thiol- modified probes, each of which specifically binds at one half of the target introduced SH groups at both ends of dsDNA. Continuous disulfide bond formation at 3′ and 5′ terminals of targets leads to the self-assembly of dsDNAs into the sulfur- rich and flexible products with different lengths. These products have a high affinity for the surface of Au-NPs and efficiently protect the surface from salt induced aggregation. To evaluate the assay efficacy, a small part of the citrus tristeza virus (CTV) genome was targeted, leading to a detection limit of about 5 × 10−9 mol.L−1 over a linear ranged from 20 × 10−9 to 10 × 10−7 mol.L−1. This approach also exhibits good reproducibility and recovery levels in the presence of plant total RNA or human plasma total circulating RNA extracts. Self-assembled targets can be then sensitively distinguished from non-assembled or mismatched targets after gel electrophoresis. The disulfide reaction method and integrating self-assembled DNAs/RNAs targets with bare AuNPs as a sensitive indicator provide us a powerful and simple visual detection tool for a wide range of applications.

A loop engineering strategy improves laccase lcc2 activity in ionic liquid and aqueous solution

Laccases, especially high redox potential laccases, play an important role in lignin degradation. The fungal laccase lcc2 from Trametes versicolor has a high redox potential and EMIM- and BMIM-based ionic liquids show excellent solubilization of wood-derived biomass. Concentrations of EMIM and BMIM to efficiently dissolve lignin impede laccase activity. Protein engineering to improve the activity and resistance of laccases in ionic liquids offers a promising opportunity for lignin valorization for the sustainable production of fuels and bulk high-value chemicals. In this work, we have performed computational assisted protein engineering of lcc2 to increase its performance in the presence of ionic liquid and aqueous solution. We showed that the loop L1 (amino acid residues 284–320) is highly important for improving lcc2 activity in EMIM EtSO4 and aqueous solutions. Lcc2 activity was improved based on a KnowVolution campaign through site saturation of seven amino acid positions identified by computational modeling. Simultaneous site saturation of four amino acid positions by OmniChange yielded variants OM1 (A285P/A310R/A312E/A318G) and OM3 (A310D/A312P/A318R) with a 3.9-fold (535.8 ± 36.9 U mg−1) and 1.6-fold (216.8 ± 5.3 U mg−1) increased specific activity in aqueous solution (lcc2 WT, 138.9 ± 6.5 U mg−1), respectively. High conservation of the loop L1 in fungal laccases suggests that computational assisted loop engineering might be used as a general strategy to improve their activity in ionic liquids and aqueous solutions.

Directed Evolution of Hyaluronic Acid Synthase from Pasteurella multocida towards High-Molecular-Weight Hyaluronic Acid

Hyaluronic acid (HA), with diverse cosmetic and medical applications, is the natural glycosaminoglycan product of HA synthases. Although process and/or metabolic engineering are used for industrial HA production, the potential of protein engineering has barely been realised. Herein, knowledge‐gaining directed evolution (KnowVolution) was employed to generate an HA synthase variant from Pasteurella multocida (pmHAS) with improved chain‐length specificity and a twofold increase in mass‐based turnover number. Seven improved pmHAS variants out of 1392 generated by error‐prone PCR were identified; eight prospective positions were saturated and the most beneficial amino acid substitutions were recombined. After one round of KnowVolution, the longest HA polymer (<4.7 MDa), through an engineered pmHAS variant in a cell‐free system, was synthesised. Computational studies showed that substitutions from the best variant (T40L, V59M and T104A) are distant from the glycosyltransferase sites and increase the flexibility of the N‐terminal region of pmHAS. Taken together, these findings suggest that the N terminus may be involved in HA synthesis and demonstrate the potential of protein engineering towards improved HA synthase activity.

A Theoretical Model of the Protochlorophyllide Oxidoreductase from a Hierarchy of Protocols

The enzyme protochlorophyllide oxidoreductase (LPOR) catalyzes the light-driven reduction of protochlorophyllide (Pchlide), a crucial step in chlorophyll biosynthesis. Molecular understanding of the photocatalytic mechanism of LPOR is essential for harnessing light energy to mediate enzymatic reactions. The absence of X-ray crystal structure has promoted the development of LPOR homology models that lack a catalytically competent active site and could not explain the variously reported spectroscopic evidence, including time-resolved optical spectroscopy data. We have refined previous structural models to account for the catalytic active site and the characteristic experimental spectral features of Pchlide binding, including the 26 cm-1 red shift of the C13(1) carbonyl stretch vibration in the mid-infrared (IR) and the 12 nm red shift of the Q x electronic band. A hierarchy of theoretical methods, including homology modeling, molecular dynamics simulations, hybrid quantum mechanics [(TD-)DFT]/molecular mechanics [AMBER] calculations, and computational vibrational and electronic spectroscopies, have been combined in an iterative protocol to reproduce experimental evidence and to predict ultrafast transient IR spectroscopic fingerprints associated with the catalytic process. The successful application to the LPOR enzyme indicates that the presented hierarchical protocol provides a general workflow to protein structure refinement.

KnowVolution Campaign of an Aryl Sulfotransferase Increases Activity toward Cellobiose

Sulfated polysaccharides such as cellulose can mimic the functionalities of pathophysiologically important glycosaminoglycans. Enzymatic sulfation offers a green chemistry route to selective (mono)sulfation of oligosaccharides (e.g., cellobiose as a building block of cellulose) in aqueous solution, at ambient temperature, and high chemoselectivity. Here, we report the first KnowVolution campaign for the aryl sulfotransferase B (ASTB) from Desulfitobacterium hafniense to advance ASTB toward a synthetically attractive biocatalyst. The generated final recombination variant (ASTB-M5) carries two amino acid substitutions (Leu446Pro and Val579Lys) leading to an up to 7.6-fold increase in specific activity (6.15 U mg-1 ) that was obtained with one round of KnowVolution. Mass spectrometry analysis confirmed a monosulfated product of cellobiose and structure elucidation by NMR confirmed the sulfation at the positions C-3 or C-4 of GlcNAc-linker-tBoc as opposed to the preferred C-6 by chemical means. Computational analysis suggested an important role of Leu446Pro in substrate-binding and recognized Val579Lys as a distal substitution.

A Comparative Reengineering Study of cpADH5 through Iterative and Simultaneous Multisite Saturation Mutagenesis

Positions identified in directed evolution campaigns or by (semi)rational design can be recombined iteratively or simultaneously. Iterative recombination has yielded many success stories and is beneficially used if screening capabilities are limited (four iterative SSMs generate 20×4=80 different enzyme variants). Simultaneous site saturation mutagenesis offers significantly higher diversity (204=160 000 variants) and enables greater improvements to be found, especially if the selected positions are in close proximity to each other (cooperative effects). Here we report a first comprehensive comparison of iterative and simultaneous saturation of four residues in Candida parapsilosis alcohol dehydrogenase 5 (cpADH5) with methyl 3‐hydroxyhexanoate as substrate. Screening of 7200 clones from 33 site saturation mutagenesis libraries (exploring 17 recombination paths) yielded the cpADH5 W286A variant, with a 82‐fold improved initial activity toward methyl 3‐hydroxyhexanoate. Screening 3500 clones from a single OmniChange library with four positions (C57, W116, L119, and W286; 1.8 % of the generated sequence space) yielded the cpADH5 C57V/W286S variant, with a 108‐fold improvement in initial activity toward methyl 3‐hydroxyhexanoate. A 1.8 % coverage of the sequence space of the simultaneous multisite saturation library was, in comparison to the investigated 17 recombination paths, sufficient to identify a cpADH5 variant with improved activity.