Protein design: reengineering cellular retinoic acid binding protein II into a rhodopsin protein mimic

Excited-State Dynamics of a CRABPII-Based Microbial Rhodopsin Mimic

Microbial rhodopsin, a pivotal photoreceptor protein, has garnered widespread application in diverse fields such as optogenetics, biotechnology, biodevices, etc. However, current microbial rhodopsins are all transmembrane proteins, which both complicates the investigation on the photoreaction mechanism and limits their further applications. Therefore, a specific mimic for microbial rhodopsin can not only provide a better model for understanding the mechanism but also can extend the applications. The human protein CRABPII turns out to be a good template for design mimics on rhodopsin due to the convenience in synthesis and the stability after mutations. Recently, Geiger et al. designed a new CRABPII-based mimic M1-L121E on microbial rhodopsin with the 13-cis, syn (13C) isomerization after irradiation. However, it still remains a question as to how similar it is compared with the natural microbial rhodopsin, in particular, in the aspect of the photoreaction dynamics. In this article, we investigate the excited-state dynamics of this mimic by measuring its transient absorption spectra. Our results reveal that there are two components in the solution of mimic M1-L121E at pH 8, known as protonated Schiff base (PSB) and unprotonated Schiff base (USB) states. In both states, the photoreaction process from 13-cis, syn(13C) to all-trans,anti (AT) is faster than that from the inverse direction. In addition, the photoreaction process in the PSB state is faster than that in the USB state. We compared the isomerization time of the PSB state to that of microbial rhodopsin. Our findings indicate that M1-L121E exhibits behaviors similar to those of microbial rhodopsins in the general pattern of PSB isomerization, where the isomerization from 13C to AT is much faster than its inverse direction. However, our results also reveal significant differences in the excited-state dynamics of the mimic relative to the native microbial rhodopsin, including the slower PSB isomerization rates as well as the unusual USB photoreaction dynamics at pH = 8. By elucidating the distinctive characteristics of mimics M1-L121E, this study enhances our understanding of microbial rhodopsin mimics and their potential applications.

Control of Protonated Schiff Base Excited State Decay within Visual Protein Mimics: A Unified Model for Retinal Chromophores

Artificial biomimetic chromophore-protein complexes inspired by natural visual pigments can feature color tunability across the full visible spectrum. However, control of excited state dynamics of the retinal chromophore, which is of paramount importance for technological applications, is lacking due to its complex and subtle photophysics/photochemistry. Here, ultrafast transient absorption spectroscopy and quantum mechanics/molecular mechanics simulations are combined for the study of highly tunable rhodopsin mimics, as compared to retinal chromophores in solution. Conical intersections and transient fluorescent intermediates are identified with atomistic resolution, providing unambiguous assignment of their ultrafast excited state absorption features. The results point out that the electrostatic environment of the chromophore, modified by protein point mutations, affects its excited state properties allowing control of its photophysics with same power of chemical modifications of the chromophore. The complex nature of such fine control is a fundamental knowledge for the design of bio-mimetic opto-electronic and photonic devices.

NeoR, a near-infrared absorbing rhodopsin

The Rhizoclosmatium globosum genome encodes three rhodopsin-guanylyl cyclases (RGCs), which are predicted to facilitate visual orientation of the fungal zoospores. Here, we show that RGC1 and RGC2 function as light-activated cyclases only upon heterodimerization with RGC3 (NeoR). RGC1/2 utilize conventional green or blue-light-sensitive rhodopsins (λmax = 550 and 480 nm, respectively), with short-lived signaling states, responsible for light-activation of the enzyme. The bistable NeoR is photoswitchable between a near-infrared-sensitive (NIR, λmax = 690 nm) highly fluorescent state (QF = 0.2) and a UV-sensitive non-fluorescent state, thereby modulating the activity by NIR pre-illumination. No other rhodopsin has been reported so far to be functional as a heterooligomer, or as having such a long wavelength absorption or high fluorescence yield. Site-specific mutagenesis and hybrid quantum mechanics/molecular mechanics simulations support the idea that the unusual photochemical properties result from the rigidity of the retinal chromophore and a unique counterion triad composed of two glutamic and one aspartic acids. These findings substantially expand our understanding of the natural potential and limitations of spectral tuning in rhodopsin photoreceptors.

Engineering of a Red Fluorogenic Protein/Merocyanine Complex for Live-Cell Imaging

A reengineered human cellular retinol binding protein II (hCRBPII), a 15-kDa protein belonging to the intracellular lipid binding protein (iLBP) family, generates a highly fluorescent red pigment through the covalent linkage of a merocyanine aldehyde to an active site lysine residue. The complex exhibits "turn-on" fluorescence, due to a weakly fluorescent aldehyde that "lights up" with subsequent formation of a strongly fluorescent merocyanine dye within the binding pocket of the protein. Cellular penetration of merocyanine is rapid, and fluorophore maturation is nearly instantaneous. The hCRBPII/merocyanine complex displays high quantum yield, low cytotoxicity, specificity in labeling organelles, and compatibility in both cancer cell lines and yeast cells. The hCRBPII/merocyanine tag is brighter than most common red fluorescent proteins.

Selective Colorimetric "Turn-On" Probe for Efficient Engineering of Iminium Biocatalysis

The efficient engineering of iminium biocatalysis has drawn considerable attention, with many applications in pharmaceutical synthesis. Here, we report a tailor-made iminium-activated colorimetric "turn-on" probe, specifically designed as a prescreening tool to facilitate engineering of iminium biocatalysis. Upon complexation of the probe with the catalytic Pro-1 residue of the model enzyme 4-oxalocrotonate tautomerase (4-OT), a brightly colored merocyanine-dye-type structure is formed. 4-OT mutants that formed this brightly colored species upon incubation with the probe proved to have a substantial activity for the iminium-based Michael-type addition of nitromethane to cinnamaldehyde, whereas mutants that showed no staining by the probe exhibited no or very low-level "Michaelase" activity. This system was exploited in a solid-phase prescreening assay termed as activated iminium colony staining (AICS) to enrich libraries for active mutants. AICS prescreening reduced the screening effort up to 20-fold. After two rounds of directed evolution, two artificial Michaelases were identified with up to 39-fold improvement in the activity for the addition of nitromethane to cinnamaldehyde, yielding the target γ-nitroaldehyde product with excellent isolated yield (up to 95%) and enantiopurity (up to >99% ee). The colorimetric activation of the turn-on probe could be extended to the class I aldolase 2-deoxy-d-ribose 5-phosphate aldolase, implicating a broader application of AICS in engineering iminium biocatalysis.

Structural Basis of Design and Engineering for Advanced Plant Optogenetics

In optogenetics, light-sensitive proteins are specifically expressed in target cells and light is used to precisely control the activity of these proteins at high spatiotemporal resolution. Optogenetics initially used naturally occurring photoreceptors to control neural circuits, but has expanded to include carefully designed and engineered photoreceptors. Several optogenetic constructs are based on plant photoreceptors, but their application to plant systems has been limited. Here, we present perspectives on the development of plant optogenetics, considering different levels of design complexity. We discuss how general principles of light-driven signal transduction can be coupled with approaches for engineering protein folding to develop novel optogenetic tools. Finally, we explore how the use of computation, networks, circular permutation, and directed evolution could enrich optogenetics.

Engineering the hCRBPII Domain-Swapped Dimer into a New Class of Protein Switches

Protein conformational switches or allosteric proteins play a key role in the regulation of many essential biological pathways. Nonetheless, the implementation of protein conformational switches in protein design applications has proven challenging, with only a few known examples that are not derivatives of naturally occurring allosteric systems. We have discovered that the domain-swapped (DS) dimer of hCRBPII undergoes a large and robust conformational change upon retinal binding, making it a potentially powerful template for the design of protein conformational switches. Atomic resolution structures of the apo- and holo-forms illuminate a simple, mechanical movement involving sterically driven torsion angle flipping of two residues that drive the motion. We further demonstrate that the conformational "readout" can be altered by addition of cross-domain disulfide bonds, also visualized at atomic resolution. Finally, as a proof of principle, we have created an allosteric metal binding site in the DS dimer, where ligand binding results in a reversible 5-fold loss of metal binding affinity. The high resolution structure of the metal-bound variant illustrates a well-formed metal binding site at the interface of the two domains of the DS dimer and confirms the design strategy for allosteric regulation.

a-ARM: Automatic Rhodopsin Modeling with Chromophore Cavity Generation, Ionization State Selection, and External Counterion Placement

The Automatic Rhodopsin Modeling (ARM) protocol has recently been proposed as a tool for the fast and parallel generation of basic hybrid quantum mechanics/molecular mechanics (QM/MM) models of wild type and mutant rhodopsins. However, in its present version, input preparation requires a few hours long user's manipulation of the template protein structure, which also impairs the reproducibility of the generated models. This limitation, which makes model building semiautomatic rather than fully automatic, comprises four tasks: definition of the retinal chromophore cavity, assignment of protonation states of the ionizable residues, neutralization of the protein with external counterions, and finally congruous generation of single or multiple mutations. In this work, we show that the automation of the original ARM protocol can be extended to a level suitable for performing the above tasks without user's manipulation and with an input preparation time of minutes. The new protocol, called a-ARM, delivers fully reproducible (i.e., user independent) rhodopsin QM/MM models as well as an improved model quality. More specifically, we show that the trend in vertical excitation energies observed for a set of 25 wild type and 14 mutant rhodopsins is predicted by the new protocol better than when using the original. Such an agreement is reflected by an estimated (relative to the probed set) trend deviation of 0.7 ± 0.5 kcal mol-1 (0.03 ± 0.02 eV) and mean absolute error of 1.0 kcal mol-1 (0.04 eV).

Bile Acid Binding Protein Functionalization Leads to a Fully Synthetic Rhodopsin Mimic

Rhodopsins are photoreceptive proteins using light to drive a plethora of biological functions such as vision, proton and ion pumping, cation and anion channeling, and gene and enzyme regulation. Here we combine organic synthesis, NMR structural studies, and photochemical characterization to show that it is possible to prepare a fully synthetic mimic of rhodopsin photoreceptors. More specifically, we conjugate a bile acid binding protein with a synthetic mimic of the rhodopsin protonated Schiff base chromophore to achieve a covalent complex featuring an unnatural protein host, photoswitch, and photoswitch-protein linkage with a reverse orientation. We show that, in spite of its molecular-level diversity, light irradiation of the prepared mimic fuels a photochromic cycle driven by sequential photochemical and thermal Z/E isomerizations reminiscent of the photocycles of microbial rhodopsins.

A Genetically Encoded Ratiometric pH Probe: Wavelength Regulation-Inspired Design of pH Indicators

Mutants of human cellular retinol-binding protein II (hCRBPII) were engineered to bind a julolidine retinal analogue for the purpose of developing a ratiometric pH sensor. The design relied on the electrostatic influence of a titratable amino acid side chain, which affects the absorption and, thus, the emission of the protein/fluorophore complex. The ratio of emissions obtained at two excitation wavelengths that correspond to the absorption of the two forms of the protein/fluorophore complex, leads to a concentration-independent measure of pH.

Free-Energy-Based Protein Design: Re-Engineering Cellular Retinoic Acid Binding Protein II Assisted by the Moveable-Type Approach

How to fine-tune the binding free energy of a small-molecule to a receptor site by altering the amino acid residue composition is a key question in protein engineering. Indeed, the ultimate solution to this problem, to chemical accuracy (±1 kcal/mol), will result in profound and wide-ranging applications in protein design. Numerous tools have been developed to address this question using knowledge-based models to more computationally intensive molecular dynamics simulations-based free energy calculations, but while some success has been achieved there remains room for improvement in terms of overall accuracy and in the speed of the methodology. Here we report a fast, knowledge-based movable-type (MT)-based approach to estimate the absolute and relative free energy of binding as influenced by mutations in a small-molecule binding site in a protein. We retrospectively validate our approach using mutagenesis data for retinoic acid binding to the Cellular Retinoic Acid Binding Protein II (CRABPII) system and then make prospective predictions that are borne out experimentally. The overall performance of our approach is supported by its success in identifying mutants that show high or even sub-nano-molar binding affinities of retinoic acid to the CRABPII system.

Predictive methods for computational metalloenzyme redesign - a test case with carboxypeptidase A

Computational metalloenzyme design is a multi-scale problem. It requires treating the metal coordination quantum mechanically, extensive sampling of the protein backbone, and additionally accounting for the polarization of the active site by both the metal cation and the surrounding protein (a phenomenon called electrostatic preorganization). We bring together a combination of theoretical methods that jointly offer these desired qualities: QM/DMD for mixed quantum-classical dynamic sampling, quantum theory of atoms in molecules (QTAIM) for the assessment of electrostatic preorganization, and Density Functional Theory (DFT) for mechanistic studies. Within this suite of principally different methods, there are both complementarity of capabilities and cross-validation. Using these methods, predictions can be made regarding the relative activities of related enzymes, as we show on the native Zn²⁺-dependent carboxypeptidase A (CPA), and its mutant proteins, which are hypothesized to hydrolyze modified substrates. For the native CPA, we replicated the catalytic mechanism and the rate in close agreement with the experiment, giving validity to the QM/DMD predicted structure, the DFT mechanism, and the QTAIM assessment of catalytic activity. For most sequences of the modified substrate and tried CPA mutants, substantially worsened activity is predicted. However, for the substrate mutant that contains Asp instead of Phe at the C-terminus, one CPA mutant exhibits a reasonable activity, as predicted across the theoretical methods. CPA is a well-studied system, and here it serves as a testing ground for the offered methods.

Domain-Swapped Dimers of Intracellular Lipid-Binding Proteins: Evidence for Ordered Folding Intermediates

Human Cellular Retinol Binding Protein II (hCRBPII), a member of the intracellular lipid-binding protein family, is a monomeric protein responsible for the intracellular transport of retinol and retinal. Herein we report that hCRBPII forms an extensive domain-swapped dimer during bacterial expression. The domain-swapped region encompasses almost half of the protein. The dimer represents a novel structural architecture with the mouths of the two binding cavities facing each other, producing a new binding cavity that spans the length of the protein complex. Although wild-type hCRBPII forms the dimer, the propensity for dimerization can be substantially increased via mutation at Tyr60. The monomeric form of the wild-type protein represents the thermodynamically more stable species, making the domain-swapped dimer a kinetically trapped entity. Hypothetically, the wild-type protein has evolved to minimize dimerization of the folding intermediate through a critical hydrogen bond (Tyr60-Glu72) that disfavors the dimeric form.

A Photoisomerizing Rhodopsin Mimic Observed at Atomic Resolution

The members of the rhodopsin family of proteins are involved in many essential light-dependent processes in biology. Specific photoisomerization of the protein-bound retinylidene PSB at a specified wavelength range of light is at the heart of all of these systems. Nonetheless, it has been difficult to reproduce in an engineered system. We have developed rhodopsin mimics, using intracellular lipid binding protein family members as scaffolds, to study fundamental aspects of protein/chromophore interactions. Herein we describe a system that specifically isomerizes the retinylidene protonated Schiff base both thermally and photochemically. This isomerization has been characterized at atomic resolution by quantitatively interconverting the isomers in the crystal both thermally and photochemically. This event is accompanied by a large pKa change of the imine similar to the pKa changes observed in bacteriorhodopsin and visual opsins during isomerization.

Light-Activated Reversible Imine Isomerization: Towards a Photochromic Protein Switch

Mutants of cellular retinoic acid-binding protein II (CRABPII), engineered to bind all-trans-retinal as an iminium species, demonstrate photochromism upon irradiation with light at different wavelengths. UV light irradiation populates the cis-imine geometry, which has a high pKa , leading to protonation of the imine and subsequent "turn-on" of color. Yellow light irradiation yields the trans-imine isomer, which has a depressed pKa , leading to loss of color because the imine is not protonated. The protein-bound retinylidene chromophore undergoes photoinduced reversible interconversion between the colored and uncolored species, with excellent fatigue resistance.

"Turn-on" protein fluorescence: in situ formation of cyanine dyes

Protein reengineering of cellular retinoic acid binding protein II (CRABPII) has yielded a genetically addressable system, capable of binding a profluorophoric chromophore that results in fluorescent protein/chromophore complexes. These complexes exhibit far-red emission, with high quantum efficiencies and brightness and also exhibit excellent pH stability spanning the range of 2-11. In the course of this study, it became evident that single mutations of L121E and R59W were most effective in improving the fluorescent characteristics of CRABPII mutants as well as the kinetics of complex formation. The readily crystallizable nature of these proteins was invaluable to provide clues for the observed spectroscopic behavior that results from single mutation of key residues.

Structures of holo wild-type human cellular retinol-binding protein II (hCRBPII) bound to retinol and retinal

Cellular retinol-binding proteins (CRBPs) I and II, which are members of the intracellular lipid-binding protein (iLBP) family, are retinoid chaperones that are responsible for the intracellular transport and delivery of both retinol and retinal. Although structures of retinol-bound CRBPI and CRBPII are known, no structure of a retinal-bound CRBP has been reported. In addition, the retinol-bound human CRBPII (hCRBPII) structure shows partial occupancy of a noncanonical conformation of retinol in the binding pocket. Here, the structure of retinal-bound hCRBPII and the structure of retinol-bound hCRBPII with retinol fully occupying the binding pocket are reported. It is further shown that the retinoid derivative seen in both the zebrafish CRBP and the hCRBPII structures is likely to be the product of flux-dependent and wavelength-dependent X-ray damage during data collection. The structures of retinoid-bound CRBPs are compared and contrasted, and rationales for the differences in binding affinities for retinal and retinol are provided.

Adherence to Bürgi-Dunitz stereochemical principles requires significant structural rearrangements in Schiff-base formation: insights from transaldolase complexes

The Bürgi-Dunitz angle (αBD) describes the trajectory of approach of a nucleophile to an electrophile. The adoption of a stereoelectronically favorable αBD can necessitate significant reactive-group repositioning over the course of bond formation. In the context of enzyme catalysis, interactions with the protein constrain substrate rotation, which could necessitate structural transformations during bond formation. To probe this theoretical framework vis-à-vis biocatalysis, Schiff-base formation was analysed in Francisella tularensis transaldolase (TAL). Crystal structures of wild-type and Lys→Met mutant TAL in covalent and noncovalent complexes with fructose 6-phosphate and sedoheptulose 7-phosphate clarify the mechanism of catalysis and reveal that substrate keto moieties undergo significant conformational changes during Schiff-base formation. Structural changes compelled by the trajectory considerations discussed here bear relevance to bond formation in a variety of constrained enzymic/engineered systems and can inform the design of covalent therapeutics.

Rational design of a colorimetric pH sensor from a soluble retinoic acid chaperone

Reengineering of cellular retinoic acid binding protein II (CRABPII) to be capable of binding retinal as a protonated Schiff base is described. Through rational alterations of the binding pocket, electrostatic perturbations of the embedded retinylidene chromophore that favor delocalization of the iminium charge lead to exquisite control in the regulation of chromophoric absorption properties, spanning the visible spectrum (474-640 nm). The pKa of the retinylidene protonated Schiff base was modulated from 2.4 to 8.1, giving rise to a set of proteins of varying colors and pH sensitivities. These proteins were used to demonstrate a concentration-independent, ratiometric pH sensor.

The photochemical determinants of color vision: revealing how opsins tune their chromophore's absorption wavelength

The evolution of a variety of important chromophore-dependent biological processes, including microbial light sensing and mammalian color vision, relies on protein modifications that alter the spectral characteristics of a bound chromophore. Three different color opsins share the same chromophore, but have three distinct absorptions that together cover the entire visible spectrum, giving rise to trichromatic vision. The influence of opsins on the absorbance of the chromophore has been studied through methods such as model compounds, opsin mutagenesis, and computational modeling. The recent development of rhodopsin mimic that uses small soluble proteins to recapitulate the binding and wavelength tuning of the native opsins provides a new platform for studying protein-regulated spectral tuning. The ability to achieve far-red shifted absorption in the rhodopsin mimic system was attributed to a combination of the lack of a counteranion proximal to the iminium, and a uniformly neutral electrostatic environment surrounding the chromophore.

An empirical test of convergent evolution in rhodopsins

Rhodopsins are photochemically reactive membrane proteins that covalently bind retinal chromophores. Type I rhodopsins are found in both prokaryotes and eukaryotic microbes, whereas type II rhodopsins function as photoactivated G-protein coupled receptors (GPCRs) in animal vision. Both rhodopsin families share the seven transmembrane α-helix GPCR fold and a Schiff base linkage from a conserved lysine to retinal in helix G. Nevertheless, rhodopsins are widely cited as a striking example of evolutionary convergence, largely because the two families lack detectable sequence similarity and differ in many structural and mechanistic details. Convergence entails that the shared rhodopsin fold is so especially suited to photosensitive function that proteins from separate origins were selected for this architecture twice. Here we show, however, that the rhodopsin fold is not required for photosensitive activity. We engineered functional bacteriorhodopsin variants with novel folds, including radical noncircular permutations of the α-helices, circular permutations of an eight-helix construct, and retinal linkages relocated to other helices. These results contradict a key prediction of convergence and thereby provide an experimental attack on one of the most intractable problems in molecular evolution: how to establish structural homology for proteins devoid of discernible sequence similarity.

Relocating the active-site lysine in rhodopsin and implications for evolution of retinylidene proteins

Type I and type II rhodopsins share several structural features including a G protein-coupled receptor fold and a highly conserved active-site Lys residue in the seventh transmembrane segment of the protein. However, the two families lack significant sequence similarity that would indicate common ancestry. Consequently, the rhodopsin fold and conserved Lys are widely thought to have arisen from functional constraints during convergent evolution. To test for the existence of such a constraint, we asked whether it were possible to relocate the highly conserved Lys296 in the visual pigment bovine rhodopsin. We show here that the Lys can be moved to three other locations in the protein while maintaining the ability to form a pigment with 11-cis-retinal and activate the G protein transducin in a light-dependent manner. These results contradict the convergent hypothesis and support the homology of type I and type II rhodopsins by divergent evolution from a common ancestral protein.

Tuning the electronic absorption of protein-embedded all-trans-retinal

Protein-chromophore interactions are a central component of a wide variety of critical biological processes such as color vision and photosynthesis. To understand the fundamental elements that contribute to spectral tuning of a chromophore inside the protein cavity, we redesigned human cellular retinol binding protein II (hCRBPII) to fully encapsulate all-trans-retinal and form a covalent bond as a protonated Schiff base. This system, using rational mutagenesis designed to alter the electrostatic environment within the binding pocket of the host protein, enabled regulation of the absorption maximum of the pigment in the range of 425 to 644 nanometers. With only nine point mutations, the hCRBPII mutants induced a systematic shift in the absorption profile of all-trans-retinal of more than 200 nanometers across the visible spectrum.

Förster resonance energy transfer competitive displacement assay for human soluble epoxide hydrolase

The soluble epoxide hydrolase (sEH), responsible for the hydrolysis of various fatty acid epoxides to their corresponding 1,2-diols, is becoming an attractive pharmaceutical target. These fatty acid epoxides, particularly epoxyeicosatrienoic acids (EETs), play an important role in human homeostatic and inflammation processes. Therefore, inhibition of human sEH, which stabilizes EETs in vivo, brings several beneficial effects to human health. Although there are several catalytic assays available to determine the potency of sEH inhibitors, measuring the in vitro inhibition constant (K(i)) for these inhibitors using catalytic assay is laborious. In addition, k(off), which has been recently suggested to correlate better with the in vivo potency of inhibitors, has never been measured for sEH inhibitors. To better measure the potency of sEH inhibitors, a reporting ligand, 1-(adamantan-1-yl)-3-(1-(2-(7-hydroxy-2-oxo-2H-chromen-4-yl)acetyl) piperidin-4-yl)urea (ACPU), was designed and synthesized. With ACPU, we have developed a Förster resonance energy transfer (FRET)-based competitive displacement assay using intrinsic tryptophan fluorescence from sEH. In addition, the resulting assay allows us to measure the K(i) values of very potent compounds to the picomolar level and to obtain relative k(off) values of the inhibitors. This assay provides additional data to evaluate the potency of sEH inhibitors.

Systematic interaction analysis of human lipocalin-type prostaglandin D synthase with small lipophilic ligands

L-PGDS [lipocalin-type PG (prostaglandin) D synthase] is a multi-functional protein, acting as a PGD2-producing enzyme and a lipid-transporter. In the present study, we focus on the function of L-PGDS as an extracellular transporter for small lipophilic molecules. We characterize the binding mechanism of human L-PGDS for the molecules, especially binding affinity stoichiometry and driving force, using tryptophan fluorescence quenching, ICD (induced circular dichroism) and ITC (isothermal titration calorimetry). The tryptophan fluorescence quenching measurements revealed that haem metabolites such as haemin, biliverdin and bilirubin bind to L-PGDS with significantly higher affinities than the other small lipophilic ligands examined, showing dissociation constant (Kd) values from 17.0 to 20.9 nM. We focused particularly on the extra-specificities of haem metabolites and L-PGDS. The ITC and ICD data revealed that two molecules of the haem metabolites bind to L-PGDS with high and low affinities, showing Kd values from 2.8 to 18.1 nM and from 0.209 to 1.63 μM respectively. The thermodynamic parameters for the interactions revealed that the contributions of enthalpy and entropy change were considerably different for each haem metabolite even when the Gibbs energy change was the same. Thus we believe that the binding energy of haem metabolites to L-PGDS is optimized by balancing enthalpy and entropy change.

Selectivity of labeled bromoethylamine for protein alkylation

Alkylation of cysteine residues has been used extensively for characterization of proteins and their mode of action in biological systems, research endeavors that are at the core of proteomics. Treatment with a simple alkylating agent such as [2-(13)C] bromoethylamine would result in labeled thialysine at the ε-position. This chemical modification of proteins would allow investigations via both (13)C NMR spectroscopy and mass spectrometry. However [2-(13)C] labeled bromoethylamine is not available commercially. We investigated its synthesis at acid pH with the goal of obtaining singly labeled bromoethylamine and understanding the mechanistic details of the reaction. Based on our experimental and theoretical results, bromination of [2-(13)C] labeled ethanolamine in acidic conditions takes place via exclusive attack of the nucleophile (HBr) at the hydroxyl bearing C. Moreover, hydrogen bonding guides the nucleophilic attack, resulting in no label scrambling of the bromoethylamine product. Protein alkylation at cysteine residue with the synthesized Br(13)CH(2)CH(2)NH(2)-HBr is successful. Ab initio calculations in which CH(3)SH serves as a model for the cysteine residue suggest that in gas phase intermolecular attack by the sulfur bearing nucleophile is favored over the intramolecular substitution by the amino group by 15.4 kJ mol(-1). Solution modeling shows that the trend is preserved at basic pH, which is the experimental one, but is reversed at neutral pH.

QM/MM refinement and analysis of protein bound retinoic acid

Retinoic acid (RA) is a vitamin A derivative, which modifies the appearance of fine wrinkles and roughness of facial skin and treats acne and activates gene transcription by binding to heterodimers of the retinoic acid receptor (RAR) and the retinoic X receptor (RXR). There are series of protein bound RA complexes available in the protein databank (PDB), which provides a broad range of information about the different bioactive conformations of RA. To gain further insights into the observed bioactive RA conformations we applied quantum mechanic (QM)/molecular mechanic (MM) approaches to re-refine the available RA protein-ligand complexes. MP2 complete basis set (CBS) extrapolations single energy calculations are also carried out for both the experimental conformations and QM optimized geometries of RA in the gas as well as solution phase. The results demonstrate that the re-refined structures show better geometries for RA than seen in the originally deposited PDB structures through the use of QMs for the ligand in the X-ray refinement procedure. QM/MM re-refined conformations also reduced the computed strain energies found in the deposited crystal conformations for RA. Finally, the dependence of ligand strain on resolution is analyzed. It is shown that ligand strain is not converged in our calculations and is likely an artifact of the typical resolutions employed to study protein-ligand complexes.

p-Coumaric acid decarboxylase from Lactobacillus plantarum: structural insights into the active site and decarboxylation catalytic mechanism

p-Coumaric acid decarboxylases (PDCs) catalyze the nonoxidative decarboxylation of hydroxycinnamic acids to generate the corresponding vinyl derivatives. Despite the biotechnological relevance of PDCs in food industry, their catalytic mechanism remains largely unknown. Here, we report insights into the structural basis of catalysis for the homodimeric PDC from Lactobacillus plantarum (LpPDC). The global fold of LpPDC is based on a flattened beta-barrel surrounding an internal cavity. Crystallographic and functional analyses of single-point mutants of residues located within this cavity have permitted identifying a potential substrate-binding pocket and also to provide structural evidences for rearrangements of surface loops so that they can modulate the accessibility to the active site. Finally, combination of the structural and functional data with in silico results enables us to propose a two-step catalytic mechanism for decarboxylation of p-coumaric acid by PDCs where Glu71 is involved in proton transfer, and Tyr18 and Tyr20 are involved in the proper substrate orientation and in the release of the CO(2) product.

Dissection of the critical binding determinants of cellular retinoic acid binding protein II by mutagenesis and fluorescence binding assay

The binding of retinoic acid to mutants of Cellular Retinoic Acid Binding Protein II (CRABPII) was evaluated to better understand the importance of the direct protein/ligand interactions. The important role of Arg111 for the correct structure and function of the protein was verified and other residues that directly affect retinoic acid binding have been identified. Furthermore, retinoic acid binding to CRABPII mutants that lack all previously identified interacting amino acids was rescued by providing a carboxylic acid dimer partner in the form of a Glu residue.

Elucidating the exact role of engineered CRABPII residues for the formation of a retinal protonated Schiff base

Cellular Retinoic Acid Binding Protein II (CRABPII) has been reengineered to specifically bind and react with all-trans-retinal to form a protonated Schiff base. Each step of this process has been dissected and four residues (Lys132, Tyr134, Arg111, and Glu121) within the CRABPII binding site have been identified as crucial for imine formation and/or protonation. The precise role of each residue has been examined through site directed mutagenesis and crystallographic studies. The crystal structure of the R132K:L121E-CRABPII (PDB-3I17) double mutant suggests a direct interaction between engineered Glu121 and the native Arg111, which is critical for both Schiff base formation and protonation.

Catalytic mechanism and performance of computationally designed enzymes for Kemp elimination

A series of enzymes for Kemp elimination of 5-nitrobenzisoxazole has been recently designed and tested. In conjunction with the design process, extensive computational analyses were carried out to evaluate the potential performance of four of the designs, as presented here. The enzyme-catalyzed reactions were modeled using mixed quantum and molecular mechanics (QM/MM) calculations in the context of Monte Carlo (MC) statistical mechanics simulations. Free-energy perturbation (FEP) calculations were used to characterize the free-energy surfaces for the catalyzed reactions as well as for reference processes in water. The simulations yielded detailed information about the catalytic mechanisms, activation barriers, and structural evolution of the active sites over the course of the reactions. The catalytic mechanism for the designed enzymes KE07, KE10(V131N), and KE15 was found to be concerted with proton transfer, generally more advanced in the transition state than breaking of the isoxazolyl N-O bond. On the basis of the free-energy results, all three enzymes were anticipated to be active. Ideas for further improvement of the enzyme designs also emerged. On the technical side, the synergy of parallel QM/MM and experimental efforts in the design of artificial enzymes is well illustrated.

Origin of the activity drop with the E50D variant of catalytic antibody 34E4 for Kemp elimination

In enzymes, multiple structural effects cooperatively lead to the high catalytic activity, while individually these effects can be small. The design of artificial enzymes requires the understanding and ability to manipulate such subtle effects. The 34E4 catalytic antibody, catalyzing Kemp elimination of 5-nitrobenzisoxazole, and its Glu50Asp (E50D) variant are the subject of the present investigation. This removal of only a methylene group yields an approximately 30-fold reduction in the rate for the catalyzed Kemp elimination. Here, the aim is to understand this difference in the catalytic performance. The mechanism of Kemp elimination catalyzed by 34E4 and the E50D mutant is elucidated using QM/MM Monte Carlo simulations and free energy perturbation theory. In both proteins, the reaction is shown to follow a single-step, concerted mechanism. In the mutant, the activation barrier rises by 2.4 kcal/mol, which corresponds to a 62-fold rate deceleration, which is in good agreement with the experimental data. The positions and functionality of the residues in the active site are monitored throughout the reaction. It is concluded that the looser contact with the base, shorter base-Asn58 contact, less favorable pi-stacking with Trp91 in the transition state of the reaction, and different solvation pattern all contribute to the reduction of the reaction rate in the E50D variant of 34E4.

Structural analysis of site-directed mutants of cellular retinoic acid-binding protein II addresses the relationship between structural integrity and ligand binding

The structural integrity of cellular retinoic acid-binding protein II (CRABPII) has been investigated using the crystal structures of CRABPII mutants. The overall fold was well maintained by these CRABPII mutants, each of which carried multiple different mutations. A water-mediated network is found to be present across the large binding cavity, extending from Arg111 deep inside the cavity to the alpha2 helix at its entrance. This chain of interactions acts as a ;pillar' that maintains the integrity of the protein. The disruption of the water network upon loss of Arg111 leads to decreased structural integrity of the protein. A water-mediated network can be re-established by introducing the hydrophilic Glu121 inside the cavity, which results in a rigid protein with the alpha2 helix adopting an altered conformation compared with wild-type CRABPII.