Novel strategy for protein exploration: high-throughput screening assisted with fuzzy neural network

To engineer proteins with desirable characteristics from a naturally occurring protein, high-throughput screening (HTS) combined with directed evolutional approach is the essential technology. However, most HTS techniques are simple positive screenings. The information obtained from the positive candidates is used only as results but rarely as clues for understanding the structural rules, which may explain the protein activity. In here, we have attempted to establish a novel strategy for exploring functional proteins associated with computational analysis. As a model case, we explored lipases with inverted enantioselectivity for a substrate p-nitrophenyl 3-phenylbutyrate from the wild-type lipase of Burkhorderia cepacia KWI-56, which is originally selective for (S)-configuration of the substrate. Data from our previous work on (R)-enantioselective lipase screening were applied to fuzzy neural network (FNN), bioinformatic algorithm, to extract guidelines for screening and engineering processes to be followed. FNN has an advantageous feature of extracting hidden rules that lie between sequences of variants and their enzyme activity to gain high prediction accuracy. Without any prior knowledge, FNN predicted a rule indicating that "size at position L167," among four positions (L17, F119, L167, and L266) in the substrate binding core region, is the most influential factor for obtaining lipase with inverted (R)-enantioselectivity. Based on the guidelines obtained, newly engineered novel variants, which were not found in the actual screening, were experimentally proven to gain high (R)-enantioselectivity by engineering the size at position L167. We also designed and assayed two novel variants, namely FIGV (L17F, F119I, L167G, and L266V) and FFGI (L17F, L167G, and L266I), which were compatible with the guideline obtained from FNN analysis, and confirmed that these designed lipases could acquire high inverted enantioselectivity. The results have shown that with the aid of bioinformatic analysis, high-throughput screening can expand its potential for exploring vast combinatorial sequence spaces of proteins.     

Inverting enantioselectivity of Burkholderia gladioli esterase EstB by directed and designed evolution

Esterase EstB from Burkholderia gladioli, showing moderate S-enantioselectivity (E(S)=6.1) in the hydrolytic kinetic resolution of methyl-beta-hydroxyisobutyrate, was subjected to directed evolution in order to reverse its enantioselectivity. After one round of ep-PCR, saturation mutagenesis and high-throughput screening, it was found that different mutations at position 152 (in the vicinity of the active site) increase, decrease and even reverse the natural enantioselectivity of this enzyme. The newly created R-enantioselectivity of the esterase mutein (E(Rapp)=1.5) has been further enhanced by a designed evolution strategy involving random mutations close to the active site. Based on the three-dimensional structure nineteen amino acid residues have been selected as mutation sites for saturation mutagenesis. Mutations at three sites (135, 253 and 351) were found to increase R-enantioselectivity. Successive rounds of saturation mutagenesis at these "hot spots" resulted in an increase in R-enantioselectivity from E(Rapp)=1.5 for the parent mutant to E(Rapp)=28.9 for the best variant which carried four amino acid substitutions. Our results prove designed evolution followed by high-throughput screening to be an efficient strategy for engineering enzyme enantioselectivity.     

Structure,mechanism, and enantioselectivity shifting of lipase LipK107 with a simple way

Because of the complex mechanisms of enzymatic reactions, no precise and simple method of understanding and controlling the chiral selectivity of enzymes has been developed. However, structure-based rational design is a powerful approach to engineering enzymes with desired catalytic activities. In this work, a simple, structure-based, large-scale in silico design and virtual screening strategy was developed and successfully applied to enzyme engineering. We first performed protein crystallization and X-ray diffraction to determine the structure of lipase LipK107, which is a novel family I.1 lipase displaying activity for both R and S isomers in chiral resolution reactions. The catalytic mechanism of family I.1, which includes LipK107, was ascertained first through comparisons of the sequences and structures of lipases from other families. The binding states of LipK107, including the energy and the conformation of complexes with the R and S enantiomers, have been evaluated by careful biocomputation to figure out the reason for the higher S selectivity. Based on this study, a simple strategy for manipulating the chiral selectivity by modulating a crucial distance in the enzyme–substrate complex and judging virtual mutations in silico is recommended. Then, a novel electrostatic interaction analysis protocol was used to design LipK107 mutants to validate our strategy. Both positive and negative mutations determined using this theoretical protocol have been implemented in wet experiments and were proved to produce the desired enantioselectivity, showing a 176% increase or 50% decrease in enantioselectivity as desired. Because of its accuracy and versatility, the strategy is promising for practical applications.     

Optical resolution of 2-(3-indolyl)propionic acid with Mucor javanicus and α-chymotrypsin

Resolution of 2-(3-indolyl)propionic acid was achieved via biocatalytic hydrolysis of its chloroethyl ester. Of the enzymes tested, Mucor javanicus lipase (R selectivity) and -chymotrypsin (S selectivity) had high reactivity and enantioselectivity (E value > 50). Neither enzyme showed enantioselectivity (E value = 1) for 2-phenylpropionic acid.     

Enantioselective inhibition of TNF-α release by thalidomide and thalidomide-analogues

The question whether the immunomodulating activity of rac-thalidomide resides in either the (−)-(S)- or the (+)-(R)-enantiomer was addressed by synthesis and separation of pure enantiomers of thalidomide-analogues which carry a methyl-group at the asymmetric carbon atom and are thus prevented from racemization. The effect of the pure enantiomers of the thalidomide-analogues and also of the enantiomers of thalidomide on relapse of TNF-α was tested in vitro by using stimulated peripheral mononuclear blood cells. Both enantiomers of thalidomide inhibited the release of TNF-α equally well at low concentrations (5 and 12.5 μg/ml) but at higher concentrations (25 and 50 μg/ml) there was a weak but statistically significant selectivity towards the (−)-(S)-enantiomer. In the case of the configuration-stable thalidomide-analogues there was a very pronounced and statistically significant enantioselectivity towards the (S)-form even at lower concentrations (≥5 μg/ml). The (S)-enantiomers of the thalidomide-analogues differed in their inhibitory potency from (−)-(S)-thalidomide suggesting that the introduction of the methyl-group increases the TNF-α-inhibitory activity while the reduction of one of the carbonyl-functions in the glutarimide-moiety to a methylene-group decreases activity. The effect of these small molecular alterations on activity and the enantioselectivity towards the (S)-enantiomers may indicate that thalidomide and its analogues directly interact with one or several cellular target-proteins. © 1996 Wiley-Liss, Inc.     

Novel route for the resolution of both enantiomers of dropropizine by using oxime esters and supported lipases of Pseudomonas cepacia

Resolution of (R)- and (S)-dropropizine which is an antitussive and central sedative therapeutic agent in high optical and chemical yields was achieved by lipases of Pseudomonas cepacia supported on ceramic particles (lipase PS-C) and on diatomite (lipase PS-D) with oxime esters in organic solvents. The influence of several factors (lipase source, structural variations in oxime esters, the amount of lipase and its recyclability) on the enantioselectivity have been investigated. Different properties were used to describe the solvents, namely the hydrophobicity (quantified by log P) and the dielectic constant (epsilon). This enzymatic acylation using oxime esters was significant as only (S)-dropropizine and (R)-dropropizine monoacetate was obtained. (R)-Dropropizine monoacetate was chemically hydrolyzed to obtain (R)-dropropizine. The highest enantioselectivity was observed when O-acetyl benzophenone oxime was used. This enzymatic resolution provides a versatile method for getting the pure enantiomers of dropropizine by effectively optimizing the various reaction parameters.     

Study on the relationship between structure and enantioselectivity of a hyperthermophilic esterase from archaeon Aeropyrum pernix K1

To enhance the enantioselectivity of a hyperthermophilic esterase from archaeon Aeropyrum pernix K1 (APE1547), a directed evolution approach is employed to generate mutant library from the native enzyme. A mutation (TBC26) is identified after one round of epPCR. The enantioselectivity of TBC26 is increased up to 2.6-fold compared to that of wild type enzyme. TBC26 contains five amino acid substitutions (R11G, L36P, V223A, I551L, A564T). The five mutation sites are spatially distant to the catalytic center. According to the published crystal structure of WT and considering the changes of secondary and tertiary structure, here we try to explain the change of enantioselectivity of the TBC26. The results suggest that the change of enantioselectivity of enzyme has a close relationship to the configuration of the enzyme.     

Partially purified Carica papaya lipase: a versatile biocatalyst for the hydrolytic resolution of (R,S)-2-arylpropionic thioesters in water-saturated organic solvents

With the hydrolytic resolution of (R,S)-naproxen 2,2,2-trifluoroethyl thioesters in water-saturated isooctane as a model system, improvements of the specific lipase activity and thermal stability were found when a crude Carica papaya lipase (CPL) was partially purified and employed as the biocatalyst. The partially purified Carica papaya lipase (PCPL) was furthermore explored as an effective enantioselective biocatalyst for the hydrolytic resolution of (R,S)-profen thioesters in water-saturated organic solvents. The kinetic analysis in water-saturated isooctane indicated that both acyl donor and acyl acceptor have profound influences on the lipase activity, E-value, and enantioselectivity. Inversion of the enantioselectivity from (S)- to (R)-thioester was found for (R,S)-fenoprofen and (R,S)-ketoprofen thioesters that contained a bulky substituent at the meta-position of 2-phenyl moiety of the acyl part. Kinetic constants for the acylation step were furthermore estimated for elucidating the kinetic data and postulating an active site model. The thermodynamic analysis indicated that the enantiomer discrimination was driven by the difference of activation enthalpy (DeltaDeltaH) and that of activation entropy (DeltaDeltaS), yet the latter was dominated for most of the reacting systems. The postulated active site model was supported from the variation of DeltaDeltaH and DeltaDeltaS with the acyl moiety, in which a good linear enthalpy-entropy compensation relationship was also illustrated. A comparison of the performances between Candida rugosa lipase (CRL) and PCPL indicated that PCPL was superior to CRL in terms of the better thermal stability, similar or better lipase activity for the fast-reacting substrate, time-course-stability, and lower enzyme cost.     

Stereospecific hydrolysis of 3-(4-methoxyphenyl)glycidic ester in supercritical carbon dioxide by immobilized lipase

In the hydrolysis of racemic 3-(4-methoxyphenyl)glycidic acid methyl ester by immobilized Mucor miehel lipase in supercritical CO₂ the initial hydrolysis rate of the (2S,3R)-form was faster than the rate of the (2R,3S) -form. The stereoisomeric excess of the (2R,3S)-form reached 87 % at 53 % total conversion level. The water content of the reaction mixture and the initial concentration of the 3-(4-methoxyphenyl) glycidic acid methyl ester had no effect on isomeric purity. The reaction rate in supercritical CO₂ was considerably faster than in toluene/water -mixture.     

Small molecule biaryl FSH receptor agonists. Part 1: Lead discovery via encoded combinatorial synthesis

High-throughput screening of two million compounds in 37 distinct encoded combinatorial libraries using FSH receptor transfected cells provided small molecule agonists such as 1 (EC(50)=3 microM) and 2 (EC(50)=3.9 microM), based on which a focused combinatorial library with a total of 31372 compounds was designed, synthesized, and screened to reveal 72 novel biaryl FSH receptor agonists such as 8a-c as well as a unique combinatorial SAR.     

Lipase-catalyzed kinetic resolutions of racemic β- and γ-thiolactones

Several racemic β- and γ-thiolactones were synthesized and kinetic resolutions of them were executed using lipases. While a lipase from Pseudomonas cepacia (PCL) showed the highest enantioselectivity for (S)-form (>99% ee_S at 53% conversion, E > 100) in the kinetic resolution of racemic -methyl-β-propiothiolactone (rac-MPTL), it showed no hydrolysis activity in the kinetic resolution of -benzyl--methyl-β-propiothiolactone (rac-BMPTL), suggesting that the changes in the size of alkyl group from rac-MPTL to rac-BMPTL leads to lower hydrolysis activity and enantioselectivity. In contrast, racemic γ-butyrothiolactones were hydrolyzed by several lipases with low enantioselectivity, whereas a lipase from Candida antarctica (CAL) showed moderate enantioselectivity for (S)-form (>99% ee_S at 76% conversion, E = 11) in the kinetic resolution of racemic -methyl-γ-butyrothiolactone (rac-MBTL). Computer-aided molecular modeling was also performed to investigate the enantioselectivites and activities of PCL toward β-propiothiolactones. The computer modeling results suggest that the alkyl side chains of β-propiothiolactones and γ-butyrothiolactones interact with amino acid residues around hydrophobic crevice, which affects the activity of PCL.     

Molecular engineering of Rhizopus oryzae lipase using a combinatorial protein library constructed on the yeast cell surface

We constructed a combinatorial yeast library through cell-surface display of the pro- and mature region of lipase from Rhizopus oryzae (ProROL) and obtained clones retaining lipase activity in fluorescent plate assay. The initial reaction rates of hydrolysis and methanolysis could be measured directly as whole-cell biocatalyst without complex treatments such as concentration, purification, and immobilization. The selected clones showed wide-ranging variation of reaction specificity. The K138R mutant showed a 1.3-fold shift of reaction specificity toward methanolysis compared to the wild type, while the V-95D, I53V, P-96S/F196Y, and Q128H/Q197L mutants showed shifts toward hydrolysis of 1.6–5.9-fold. Predictions of the mutants’ three-dimensional structure suggested that the hydrogen-bond distance between threonine 83 and aspartic acid 92 may influence reaction specificity, which shifted toward hydrolysis in mutants where this distance was shorter than in the wild type, but toward methanolysis where it was longer. The positions of amino acid residues (aa) 53, 138 and 196 in ProROL are considered the sites that influence hydrogen-bond distance and change reaction specificity. Construction of a surface-displayed combinatorial library in yeast cells is thus a powerful tool in accelerating the combinatorial approach to enzyme engineering and novel whole-cell biocatalyst development.     

An antiproliferative genetic screening identifies a peptide aptamer that targets calcineurin and up-regulates its activity

Peptide aptamers are combinatorial recognition molecules that consist of a constant scaffold protein displaying a doubly constrained variable peptide loop. They bind specifically target proteins and interfere with their function. We have built a peptide aptamer library in a lentiviral expression system to isolate aptamers that inhibit cell proliferation in vitro. Using one of the isolated aptamers (R5G42) as a bait protein, we have performed yeast two-hybrid screening of cDNA libraries and identified calcineurin A as a target protein candidate. R5G42 bound calcineurin A in vitro and stimulated its phosphatase activity. When expressed transiently in human cells, R5G42 induced the dephosphorylation of BAD. We have identified an antiproliferative peptide aptamer that binds calcineurin and stimulates its activity. The use of this ligand may help elucidate the still elusive structural mechanisms of activation and inhibition of calcineurin. Our work illustrates the power of phenotypic screening of combinatorial protein libraries to interrogate the proteome and chart molecular regulatory networks.     

易错PCR提高华根霉脂肪酶的热稳定性

为了提高华根霉Rhizopus chinensis CCTCC M201021脂肪酶的热稳定性,运用定向进化-易错PCR的方法,经两轮易错PCR引入突变,利用fast-blue RR顶层琼脂法对突变文库进行筛选,第一轮易错PCR后筛选到2株突变菌株,第二轮筛选到4株突变株。第二轮最佳突变株Ep2-4,其中3个氨基酸发生了突变：A129S、P168L和V329A。该突变酶ep2-4在60 ℃下半衰期相对原始酶r27RCL提高5.4倍,T50值提高7.8 ℃。酶学性质研究表明,突变酶ep2-4在热稳定性提高的基础上,仍保有良好的催化活性。蛋白质三维结构模拟显示,突变A129S可以和Gln133形成氢键,增加了酶表面的亲水性和极性;P168L可以与邻近的Leu164形成疏水键,导致突变酶的热稳定性提高。     

Enantioselective disposition of oral amlodipine in healthy volunteers

Plasma concentrations of (R)- and (S)-amlodipine were measured after single oral administrations to 18 healthy volunteers of 20 mg amlodipine racemate. The contribution of the pharmacologically active (S)-enantiomer to the concentrations of total amlodipine (sum of enantiomers) was significantly higher than that of the inactive (R)-enantiomer, with mean values of 47% R to 53% S for the C_max and 41% R to 59% S for the AUC (range between 24% R:76% S and 50% R:50% S). The oral clearance of the active (S)-form was subject to much less intersubject variation (25% CV) than that of the inactive (R)-form (52% CV). (R)-Amlodipine was more rapidly eliminated from plasma than (S)-amlodipine, with mean terminal half-lives of 34.9 h (R) and 49.6 h (S). The terminal half-lives of total amlodipine (mean 44.2 h) were strongly correlated with—and thus highly predictive for—the half-lives of the (S)-enantiomer. It is proposed that the observed enantioselectivity of oral amlodipine is due to differences in the systemic blood clearance of the enantiomers. © 1994 Wiley-Liss, Inc.     

Substitution of Val72 residue alters the enantioselectivity and activity of Penicillium expansum lipase

Error-prone PCR was used to create more active or enantioselective variants of Penicillium expansum lipase (PEL). A variant with a valine to glycine substitution at residue 72 in the lid structure exhibited higher activity and enantioselectivity than those of wild-type PEL. Site-directed saturation mutagenesis was used to explore the sequence-function relationship and the substitution of Val72 of P. expansum lipase changed both catalytic activity and enantioselectivity greatly. The variant V72A, displayed a highest enantioselectivity enhanced to about twofold for the resolution of (R, S)-naproxen (E value increased from 104 to 200.7 for wild-type PEL and V72A variant, respectively). In comparison to PEL, the variant V72A showed a remarkable increase in specific activity towards p-nitrophenyl palmitate (11- and 4-fold increase at 25 and 35?°C, respectively) whereas it had a decreased thermostability. The results suggest that the enantioselective variant V72A could be used for the production of pharmaceutical drugs such as enantiomerically pure (S)-naproxen and the residue Val 72 of P. expansum lipase plays a significant role in the enantioselectivity and activity of this enantioselective lipase.     

Sulfonamide derivatives as new potent and selective CB2 cannabinoid receptor agonists

A novel series of sulfonamide derivatives 3, the CB(2) receptor agonists, was synthesized and evaluated for activity against the human CB(2) receptor. We first identified sulfonamide 3a, which was obtained by random screening of our in-house chemical library as a moderately active (CB(2) IC(50)=340nM) CB(2) receptor agonist. We then attempted to test its analogues to identify compounds with a high affinity for the CB(2) receptor. One of these, compound 3f, exhibited high affinity for the human CB(2) receptor (IC(50)=16nM) and high selectivity for CB(2) over CB(1) (CB(1) IC(50)/CB(2)IC(50)=106), and behaved as a full CB(2) receptor agonist in the [(35)S]GTPgammaS binding assay (CB(2) EC(50)=7.2nM, E(max)=100%).     

Identification of a potent and selective σ₁ receptor agonist potentiating NGF-induced neurite outgrowth in PC12 cells

Herein we report the synthesis, drug-likeness evaluation, and in vitro studies of new sigma (σ) ligands based on arylalkenylaminic scaffold. For the most active olefin the corresponding arylalkylamine was studied. Novel arylalkenylamines generally possess high σ(1) receptor affinity (K(i) values <25 nM) and good σ(1)/σ(2) selectivity (K(i)σ(2) >100). Particularly, the piperidine derivative (E)-17 and its arylalkylamine analog (R,S)-33 were observed to be excellent σ(1) receptor ligands (K(i)=0.70 and 0.86 nM, respectively) and to display significantly high selectivity over σ(2), μ-, and κ-opioid receptors and phencyclidine (PCP) binding site of the N-methyl-d-aspartate (NMDA) receptors. Moreover in PC12 cells (R,S)-33 promoted the nerve growth factor (NGF)-induced neurite outgrowth and elongation. Co-administration of the selective σ(1) receptor antagonist BD-1063 totally counteracted this effect, confirming that σ(1) receptors are involved in the (R,S)-33 modulation of the NGF effect in PC12 cells and suggesting a σ(1) agonist profile. As a part of our work, a threedimensional σ(1) pharmacophore model was also developed employing GALAHAD methodology. Only active compounds were used for deriving this model. The model included two hydrophobes and a positive nitrogen as relevant features and it was able to discriminate between molecules with and without affinity toward σ(1) receptor subtype.