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.

Probing the role of the residues in the active site of the transaminase from Thermobaculum terrenum.

Creating biocatalysts for (R)-selective amination effectively is highly desirable in organic synthesis. Despite noticeable progress in the engineering of (R)-amine activity in pyridoxal-5'-phosphate-dependent transaminases of fold type IV, the specialization of the activity is still an intuitive task, as there is poor understanding of sequence-structure-function relationships. In this study, we analyzed this relationship in transaminase from Thermobaculum terrenum, distinguished by expanded substrate specificity and activity in reactions with L-amino acids and (R)-(+)-1-phenylethylamine using α-ketoglutarate and pyruvate as amino acceptors. We performed site-directed mutagenesis to create a panel of the enzyme variants, which differ in the active site residues from the parent enzyme to a putative transaminase specific to (R)-primary amines. The variants were examined in the overall transamination reactions and half-reaction with (R)-(+)-1-phenylethylamine. A structural analysis of the most prominent variants revealed a spatial reorganization in the active sites, which caused changes in activity. Although the specialization to (R)-amine transaminase was not implemented, we succeeded in understanding the role of the particular active site residues in expanding substrate specificity of the enzyme. We showed that the specificity for (R)-(+)-1-phenylethylamine in transaminase from T. terrenum arises without sacrificing the specificity for L-amino acids and α-ketoglutarate and in consensus with it.

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.

Biochemical and structural insights into PLP fold type IV transaminase from Thermobaculum terrenum.

The high catalytic efficiency of enzymes under reaction conditions is one of the main goals in biocatalysis. Despite the dramatic progress in the development of more efficient biocatalysts by protein design, the search for natural enzymes with useful properties remains a promising strategy. The pyridoxal 5'-phosphate (PLP)-dependent transaminases represent a group of industrially important enzymes due to their ability to stereoselectively transfer amino groups between diverse substrates; however, the complex mechanism of substrate recognition and conversion makes the design of transaminases a challenging task. Here we report a detailed structural and kinetic study of thermostable transaminase from the bacterium Thermobaculum terrenum (TaTT) using the methods of enzyme kinetics, X-ray crystallography and molecular modeling. TaTT can convert L-branched-chain and L-aromatic amino acids as well as (R)-(+)-1-phenylethylamine at a high rate and with high enantioselectivity. The structures of TaTT in complex with the cofactor pyridoxal 5'-phosphate covalently bound to enzyme and in complex with its reduced form, pyridoxamine 5'-phosphate, were determined at resolutions of 2.19 Å and 1.5 Å, and deposited in the Protein Data Bank as entries 6GKR and 6Q8E, respectively. TaTT is a fold type IV PLP-dependent enzyme. In terms of structural similarity, the enzyme is close to known branched-chain amino acid aminotransferases, but differences in characteristic sequence motifs in the active site were observed in TaTT compared to canonical branched-chain amino acid aminotransferases, which can explain the improved binding of aromatic amino acids and (R)-(+)-1-phenylethylamine. This study has shown for the first time that high substrate specificity towards both various l-amino acids and (R)-primary amines can be implemented within one pyridoxal 5'-phosphate-dependent active site of fold type IV. These results complement our knowledge of the catalytic diversity of transaminases and indicate the need for further biochemical and bioinformatic studies to understand the sequence-structure-function relationship in these enzymes.

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.

Experimental and computational studies on the unusual substrate specificity of branched-chain amino acid aminotransferase from Thermoproteus uzoniensis.

PLP-Dependent fold-type IV branched-chain amino acid aminotransferases (BCATs) from archaea have so far been poorly characterized. A new BCAT from the hyperthermophilic archaeon Thermoproteus uzoniensis (TUZN1299) has been studied. TUZN1299 was found to be highly active toward branched-chain amino acids (BCAAs), positively charged amino acids, l-methionine, l-threonine, l-homoserine, l-glutamine, as well as toward 2-oxobutyrate and keto analogs of BCAAs, whereas l-glutamate and α-ketoglutarate were not converted in the overall reaction. According to stopped-flow experiments, the enzyme showed the highest specificity to BCAAs and their keto analogs. In order to explain the molecular mechanism of the unusual specificity of TUZN1299, bioinformatic analysis was implemented to identify the subfamily-specific positions in the aminotransferase class IV superfamily of enzymes. The role of the selected residues in binding of various ligands in the active site was further studied using molecular modeling. The results indicate that Glu188 forms a novel binding site for positively charged and polar side-chains of amino acids. Lack of accommodation for α-ketoglutarate and l-glutamate is due to the unique orientation and chemical properties of residues 102-106 in the loop forming the A-pocket. The likely functional roles of TUZN1299 in cellular metabolism - in the synthesis and degradation of BCAAs - are discussed.