Preparative Biocatalytic Synthesis of α-Ketomethylselenobutyrate-A Putative Agent for Cancer Therapy.

Biomedical studies of the role of organic selenium compounds indicate that the amino acid derivative of L-selenomethionine, α-ketomethylselenobutyrate (KMSB), can be considered a potential anticancer therapeutic agent. It was noted that, in addition to a direct effect on redox signaling molecules, α-ketoacid metabolites of organoselenium compounds are able to change the status of histone acetylation and suppress the activity of histone deacetylases in cancer cells. However, the wide use of KMSB in biomedical research is hindered not only by its commercial unavailability, but also by the fact that there is no detailed information in the literature on possible methods for the synthesis of this compound. This paper describes in detail the procedure for obtaining a high-purity KMSB preparation (purity ≥ 99.3%) with a yield of the target product of more than 67%. L-amino acid oxidase obtained from C. adamanteus was used as a catalyst for the conversion of L-selenomethionine to KMSB. If necessary, this method can be used as a basis both for scaling up the synthesis of KMSB and for developing cost-effective biocatalytic technologies for obtaining other highly purified drugs.

Natural Guanine Derivatives Exert PARP-Inhibitory and Cytoprotective Effects in a Model of Cardiomyocyte Damage under Oxidative Stress.

Inhibitors of human poly(ADP-ribose) polymerase (PARP) are considered as promising agents for treatment of cardiovascular, neurological, and other diseases accompanied by inflammation and oxidative stress. Previously, the ability of natural compounds 7-methylguanine (7mGua) and 8-hydroxy-7-methylguanine (8h7mGua) to suppress activity of the recombinant PARP protein was demonstrated. In the present work, we have investigated the possibility of PARP-inhibitory and cytoprotective action of 7mGua and 8h7mGua against the rat cardiomyoblast cultures (undifferentiated and differentiated H9c2). It was found that 7mGua and 8h7mGua rapidly penetrate into the cells and effectively suppress the H2O2-stimulated PARP activation (IC50 = 270 and 55 µM, respectively). The pronounced cytoprotective effects of 7mGua and 8h7mGua were shown in a cellular model of oxidative stress, and effectiveness of 8h7mGua exceeded the classic PARP inhibitor 3-aminobenzamide. The obtained data indicate promise for the development of PARP inhibitors based on guanine derivatives and their testing using the models of ischemia-reperfusion tissue damage.

Inhibitory Effects of 7-Methylguanine and Its Metabolite 8-Hydroxy-7-Methylguanine on Human Poly(ADP-Ribose) Polymerase 1.

Previously, we have found that a nucleic acid metabolite, 7-methylguanine (7mGua), produced in the body can have an inhibitory effect on the poly(ADP-ribose) polymerase 1 (PARP1) enzyme, an important pharmacological target in anticancer therapy. In this work, using an original method of analysis of PARP1 activity based on monitoring fluorescence anisotropy, we studied inhibitory properties of 7mGua and its metabolite, 8-hydroxy-7-methylguanine (8h7mGua). Both compounds inhibited PARP1 enzymatic activity in a dose-dependent manner, however, 8h7mGua was shown to be a stronger inhibitor. The IC50 values for 8h7mGua at different concentrations of the NAD+ substrate were found to be 4 times lower, on average, than those for 7mGua. The more efficient binding of 8h7mGua in the PARP1 active site is explained by the presence of an additional hydrogen bond with the Glu988 catalytic residue. Experimental and computational studies did not reveal the effect of 7mGua and 8h7mGua on the activity of other DNA repair enzymes, indicating selectivity of their inhibitory action.

Modeling the Structure of Human tRNA-Guanine Transglycosylase in Complex with 7-Methylguanine and Revealing the Factors that Determine the Enzyme Interaction with Inhibitors.

tRNA-guanine transglycosylase, an enzyme catalyzing replacement of guanine with queuine in human tRNA and participating in the translation mechanism, is involved in the development of cancer. However, information on the small-molecule inhibitors that can suppress activity of this enzyme is very limited. Molecular dynamics simulations were used to determine the amino acid residues that provide efficient binding of inhibitors in the active site of tRNA-guanine transglycosylase. It was demonstrated using 7-methylguanine molecule as a probe that the ability of the inhibitor to adopt a charged state in the environment of hydrogen bond acceptors Asp105 and Asp159 plays a key role in complex formation. Formation of the hydrogen bonds and hydrophobic contacts with Gln202, Gly229, Phe109, and Met259 residues are also important. It has been predicted that introduction of the substituents would have a different effect on the ability to inhibit tRNA-guanine transglycosylase, as well as the DNA repair protein poly(ADP-ribose) polymerase 1, which can contribute to the development of more efficient and selective compounds.

Guide tree optimization with genetic algorithm to improve multiple protein 3D-structure alignment.

MOTIVATION: With the increasing availability of 3D-data, the focus of comparative bioinformatic analysis is shifting from protein sequence alignments toward more content-rich 3D-alignments. This raises the need for new ways to improve the accuracy of 3D-superimposition. RESULTS: We proposed guide tree optimization with genetic algorithm (GA) as a universal tool to improve the alignment quality of multiple protein 3D-structures systematically. As a proof of concept, we implemented the suggested GA-based approach in popular Matt and Caretta multiple protein 3D-structure alignment (M3DSA) algorithms, leading to a statistically significant improvement of the TM-score quality indicator by up to 220-1523% on 'SABmark Superfamilies' (in 49-77% of cases) and 'SABmark Twilight' (in 59-80% of cases) datasets. The observed improvement in collections of distant homologies highlights the potentials of GA to optimize 3D-alignments of diverse protein superfamilies as one plausible tool to study the structure-function relationship. AVAILABILITY AND IMPLEMENTATION: The source codes of patched gaCaretta and gaMatt programs are available open-access at https://github.com/n-canter/gamaps. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

vsFilt: A Tool to Improve Virtual Screening by Structural Filtration of Docking Poses.

The ability of ligands to form crucial interactions with a protein target, characteristic for the substrate and/or inhibitors, could be considered a structural criterion for identifying potent binders among docked compounds. Structural filtration of predicted poses improves the performance of virtual screening and helps in recovering specifically bound ligands. Here, we present vsFilt-a highly automated and easy-to-use Web server for postdocking structural filtration. The new tool can detect various types of interactions that are known to be involved in the molecular recognition, including hydrogen and halogen bonds, ionic interactions, hydrophobic contacts, π-stacking, and cation-π interactions. A case study for poly(ADP-ribose) polymerase 1 ligands illustrates the utility of the software. The Web server is freely available at https://biokinet.belozersky.msu.ru/vsfilt.

parMATT: parallel multiple alignment of protein 3D-structures with translations and twists for distributed-memory systems.

MOTIVATION: Accurate structural alignment of proteins is crucial at studying structure-function relationship in evolutionarily distant homologues. Various software tools were proposed to align multiple protein 3D-structures utilizing one CPU and thus are of limited productivity at large-scale analysis of protein families/superfamilies. RESULTS: The parMATT is a hybrid MPI/pthreads/OpenMP parallel re-implementation of the MATT algorithm to align multiple protein 3D-structures by allowing translations and twists. The parMATT can be faster than MATT on a single multi-core CPU, and provides a much greater speedup when executed on distributed-memory systems, i.e. computing clusters and supercomputers hosting memory-independent computing nodes. The most computationally demanding steps of the MATT algorithm-the initial construction of pairwise alignments between all input structures and further iterative progression of the multiple alignment-were parallelized using MPI and pthreads, and the concluding refinement step was optimized by introducing the OpenMP support. The parMATT can significantly accelerate the time-consuming process of building a multiple structural alignment from a large set of 3D-records of homologous proteins. AVAILABILITY AND IMPLEMENTATION: The source code is available at https://biokinet.belozersky.msu.ru/parMATT. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Mustguseal: a server for multiple structure-guided sequence alignment of protein families.

Motivation: Comparative analysis of homologous proteins in a functionally diverse superfamily is a valuable tool at studying structure-function relationship, but represents a methodological challenge. Results: The Mustguseal web-server can automatically build large structure-guided sequence alignments of functionally diverse protein families that include thousands of proteins basing on all available information about their structures and sequences in public databases. Superimposition of protein structures is implemented to compare evolutionarily distant relatives, whereas alignment of sequences is used to compare close homologues. The final alignment can be downloaded for a local use or operated on-line with the built-in interactive tools and further submitted to the integrated sister web-servers of Mustguseal to analyze conserved, subfamily-specific and co-evolving residues at studying a protein function and regulation, designing improved enzyme variants for practical applications and selective ligands to modulate functional properties of proteins. Availability and implementation: Freely available on the web at https://biokinet.belozersky.msu.ru/mustguseal. Contact: vytas@belozersky.msu.ru. Supplementary information: Supplementary data are available at Bioinformatics online.

Molecular modeling of different substrate-binding modes and their role in penicillin acylase catalysis.

Molecular modeling was addressed to understand different substrate-binding modes and clarify the role of two positively charged residues of the penicillin G acylase active site - ßR263 and αR145 - in binding of negatively charged substrates. Although the electrostatic contribution to productive substrate binding was dominated by ßR263 rather than αR145, it was found that productive binding was not the only possible mode of substrate placement in the active site. Two extra binding modes - nonproductive and preproductive - were located by means of molecular docking and dynamics with binding affinities comparable with the productive one. A unique feature of nonproductive and preproductive complexes was that the substrate's acyl group did not penetrate the hydrophobic pocket, but occupied a patch on the protein interface spanning from ßR263 to αR145. Nonproductive and preproductive complexes competed with each other and productive binding mode, giving rise to increased apparent substrate binding. Preproductive complex revealed an ability to switch to a productive one during molecular dynamics simulations, and conformational plasticity of the penicillin G acylase active site was shown to be crucial for that. Nonproductive binding observed at molecular modeling corresponded well with experimentally observed substrate inhibition in penicillin acylase catalysis. By combining estimated free energies of substrate binding in each mode, and accounting for two possible conformations of the penicillin G acylase active site (closed and open) quantitative agreement with experimentally measured K(M) values was achieved. Calculated near-attack conformation frequencies from corresponding molecular dynamics simulations were in a quantitative correlation with experimental k(cat) values and demonstrated adequate application of molecular modeling methods.

Molecular modeling of formate dehydrogenase: the formation of the Michaelis complex.

The formation of the reactive enzyme-substrate complex of formate dehydrogenase has been investigated by molecular dynamics techniques accounting for different conformational states of the enzyme. Simulations revealed that the transport of substrate to the active site through the substrate channel proceeds in the open conformation of enzyme due to the crucial role of the Arg284 residue acting as a vehicle. However, formate binding in the active site of the open conformation leads to the formation of a nonproductive enzyme-substrate complex. The productive Michaelis complex is formed only in the closed enzyme conformation after the substrate and coenzyme have bound, when required rigidity of the binding site and reactive formate orientation due to interactions with Arg284, Asn146, Ile122, and His332 residues is attained. Then, the high occupancy (up to 75%) of the reactive substrate-coenzyme conformation is reached, which was demonstrated by hybrid quantum mechanics/molecular mechanics simulations using various semiempirical Hamiltonians.

Guidelines for reporting of biocatalytic reactions.

Enzymes and whole cells are being increasingly applied in research and industry, but the adoption of biocatalysis relies strongly on useful scientific literature. Unfortunately, too many published papers lack essential information needed to reproduce and understand the results. Here, members of the scientific committee of the European Federation of Biotechnology Section on Applied Biocatalysis (ESAB) provide practical guidelines for reporting experiments. The document embraces the recommendations of the STRENDA initiative (Standards for Reporting Enzymology Data) in the context of pure enzymology and provides further guidelines and explanations on topics of crucial relevance for biocatalysis. In particular, guidelines are given on issues such as the selectivity, specificity, productivity and stability of biocatalysts, as well as on methodological problems related to reactions in multiphase systems. We believe that adoption and use of these guidelines could greatly increase the value and impact of published work in biocatalysis, and hence promote the further growth of applications.

Thermodynamic and kinetic stability of penicillin acylase from Escherichia coli.

Thermal denaturation of penicillin acylase (PA) from Escherichia coli has been studied by high-sensitivity differential scanning calorimetry as a function of heating rate, pH and urea concentration. It is shown to be irreversible and kinetically controlled. Upon decrease in the heating rate from 2 to 0.1 K min(-1) the denaturation temperature of PA at pH 6.0 decreases by about 6 degrees C, while the denaturation enthalpy does not change notably giving an average value of 31.6+/-2.1 J g(-1). The denaturation temperature of PA reaches a maximum value of 64.5 degrees C at pH 6.0 and decreases by about of 15 degrees C at pH 3.0 and 9.5. The pH induced changes in the denaturation enthalpy follow changes in the denaturation temperature. Increasing the urea concentration causes a decrease in both denaturation temperature and enthalpy of PA, where denaturation temperature obeys a linear relation. The heat capacity increment of PA is not sensitive to the heating rate, nor to pH, and neither to urea. Its average value is of 0.58+/-0.02 J g(-1) K(-1). The denaturation transition of PA is approximated by the Lumry-Eyring model. The first stage of the process is assumed to be a reversible unfolding of the alpha-subunit. It activates the second stage involving dissociation of two subunits and subsequent denaturation of the beta-subunit. This stage is irreversible and kinetically controlled. Using this model the temperature, enthalpy and free energy of unfolding of the alpha-subunit, and a rate constant of the irreversible stage are determined as a function of pH and urea concentration. Structural features of the folded and unfolded conformation of the alpha-subunit as well as of the transition state of the PA denaturation in aqueous and urea solutions are discussed.

Saturation mutagenesis reveals the importance of residues alphaR145 and alphaF146 of penicillin acylase in the synthesis of beta-lactam antibiotics.

Penicillin acylase (PA) from Escherichia coli can catalyze the coupling of an acyl group to penicillin- and cephalosporin-derived beta-lactam nuclei, a conversion that can be used for the industrial synthesis of beta-lactam antibiotics. The modest synthetic properties of the wild-type enzyme make it desirable to engineer improved mutants. Analysis of the crystal structure of PA has shown that residues alphaR145 and alphaF146 undergo extensive repositioning upon binding of large ligands to the active site, suggesting that these residues may be good targets for mutagenesis aimed at improving the catalytic performance of PA. Therefore, site-saturation mutagenesis was performed on both positions and a complete set of all 38 variants was subjected to rapid HPLC screening for improved ampicillin synthesis. Not less than 33 mutants showed improved synthesis, indicating the importance of the mutated residues in PA-catalyzed acyl transfer kinetics. In several mutants at low substrate concentrations, the maximum level of ampicillin production was increased up to 1.5-fold, and the ratio of the synthetic rate over the hydrolytic rate was increased 5-15-fold. Moreover, due to increased tendency of the acyl-enzyme intermediate to react with beta-lactam nucleophile instead of water, mutants alphaR145G, alphaR145S and alphaR145L demonstrated an enhanced synthetic yield over wild-type PA at high substrate concentrations. This was accompanied by an increased conversion of 6-APA to ampicillin as well as a decreased undesirable hydrolysis of the acyl donor. Therefore, these mutants are interesting candidates for the enzymatic production of semi-synthetic beta-lactam antibiotics.

Efficient enantiomeric analysis of primary amines and amino alcohols by high-performance liquid chromatography with precolumn derivatization using novel chiral SH-reagent N-(R)-mandelyl-(S)-cysteine.

Novel N-acylated-(S)-cysteine derivative-N-(R)-mandelyl-(S)-cysteine (R-NMC), containing additional chiral center, aromatic and polar alpha-substituents in contrast to the traditionally used enantiomerically pure thiols, has been demonstrated to be an efficient SH-reagent for enantiomeric HPLC analysis of primary nonfunctionalized amines and amino alcohols after precolumn derivatization with o-phthalaldehyde. The R-NMC-derived isoindoles as well as adducts formed using traditional SH-reagents had a characteristic absorption maximum at 340 nm with a molar absorbance 6000 M(-1) cm(-1), were stable during the HPLC-analysis and highly fluorescent allowing to detect 1 fmol of amino compound. Using diastereomeric R-NMC all tested amino alcohols were resolved effectively as well as nonfunctionalized amines, some of which were not resolved by a direct method on a chiral phase. Applying traditional enantiomeric N-acetyl-(S)-cysteine (NAC) only some isoindoles formed by aliphatic amino alcohols have been separated satisfactorily. The enhanced selectivity for R-NMC-derived isomers has been achieved, obviously, due to the involvement of the substituents at an extra chiral center into additional intramolecular interactions.

Penicillin acylase-catalyzed synthesis of beta-lactam antibiotics in highly condensed aqueous systems: beneficial impact of kinetic substrate supersaturation.

Advantages of performing penicillin acylase-catalyzed synthesis of new penicillins and cephalosporins by enzymatic acyl transfer to the beta-lactam antibiotic nuclei in the supersaturated solutions of substrates have been demonstrated. It has been shown that the effective nucleophile reactivity of 6-aminopenicillanic (6-APA) and 7-aminodesacetoxycephalosporanic (7-ADCA) acids in their supersaturated solutions continue to grow proportionally to the nucleophile concentration. As a result, synthesis/hydrolysis ratio in the enzymatic synthesis can be significantly (up to three times) increased due to the nucleophile supersaturation. In the antibiotic nuclei conversion to the target antibiotic the remarkable improvement (up to 14%) has been gained. Methods of obtaining relatively stable supersaturated solutions of 6-APA, 7-ADCA, and D-p-hydroxyphenylglycine amide (D-HPGA) have been developed and syntheses of ampicillin, amoxicillin, and cephalexin starting from the supersaturated homogeneous solutions of substrates were performed. Higher synthetic efficiency and increased productivity of these reactions compared to the heterogeneous "aqueous solution-precipitate" systems were observed. The suggested approach seems to be an effective solution for the aqueous synthesis of the most widely requested beta-lactam antibiotics (i.e., amoxicillin, cephalexin, cephadroxil, cephaclor, etc.).

Quantitative characterization of the nucleophile reactivity in penicillin acylase-catalyzed acyl transfer reactions.

Nucleophile reactivity of two most known nuclei of penicillins and cephalosporins, 6-aminopenicillanic (6-APA) and 7-aminodesacetoxycephalosporanic (7-ADCA) acids, was quantitatively characterized. In penicillin acylase (PA)-catalyzed acyl transfer reactions the relative reactivity of the added nucleophile compared to the water (i.e. nucleophile reactivity) is defined by two complex kinetic parameters beta(0) and gamma, and depends on the nucleophile concentration. In turn, parameters beta(0) and gamma were shown to be dependent on the structure of both reactants involved: nucleophile and acyl donor. Analysis of the kinetic scheme revealed that nucleophile reactivity is one of a few key parameters controlling efficiency of PA-catalyzed acyl transfer to the added nucleophile in an aqueous medium. Computation of the maximum nucleophile conversion to the product using determined nucleophile reactivity parameters in the synthesis of three different antibiotics, ampicillin, amoxicillin and cephalexin, showed good correlation with the results of corresponding synthetic experiments. Suggested approach can be extended to the quantitative description and optimization of PA-catalyzed acyl transfer reactions in a wide range of experimental conditions.

Penicillin acylase-catalyzed ampicillin synthesis using a pH gradient: a new approach to optimization.

The penicillin acylase-catalyzed synthesis of ampicillin by acyl transfer from D-(-)-phenylglycine amide (D-PGA) to 6-aminopenicillanic acid (6-APA) becomes more effective when a judiciously chosen pH gradient is applied in the course of the process. This reaction concept is based on two experimental observations: 1) The ratio of the initial synthesis and hydrolysis rates (V(S)/V(H)) is pH-dependent and exhibits a maximum at pH 6.5-7.0 for a saturated solution of 6-APA; 2) at a fixed 6-APA concentration below saturation, V(S)/V(H) increases with decreasing pH. Optimum synthetic efficiency could, therefore, be achieved by starting with a concentrated 6-APA solution at pH 7 and gradually decreasing the pH to 6.3 in the course of 6-APA consumption. A conversion of 96% of 6-APA and 71% of D-PGA into ampicillin was accomplished in an optimized procedure, which significantly exceeds the efficiency of enzymatic synthesis performed at a constant pH of either 7.0 or 6.3.

Active site titration as a tool for the evaluation of immobilization procedures of penicillin acylase.

Native and immobilized preparations of penicillin acylase from Escherichia coli and Alcaligenes faecalis were studied using an active site titration technique. Knowledge of the number of active sites allowed the calculation of the average turnover rate of the enzyme in the various preparations and allowed us to quantify the contribution of irreversible inactivation of the enzyme to the loss of catalytic activity during the immobilization procedure. In most cases a loss of active sites as well as a decrease of catalytic activity per active site (turnover rate) was observed upon immobilization. Immobilization techniques affected the enzymes differently. The effect of increased loading of penicillin acylase on the average turnover rate was determined by active site titration to assess diffusion limitations in the carrier.