ChannelsDB 2.0: a comprehensive database of protein tunnels and pores in AlphaFold era.

ChannelsDB 2.0 is an updated database providing structural information about the position, geometry and physicochemical properties of protein channels-tunnels and pores-within deposited biomacromolecular structures from PDB and AlphaFoldDB databases. The newly deposited information originated from several sources. Firstly, we included data calculated using a popular CAVER tool to complement the data obtained using original MOLE tool for detection and analysis of protein tunnels and pores. Secondly, we added tunnels starting from cofactors within the AlphaFill database to enlarge the scope of the database to protein models based on Uniprot. This has enlarged available channel annotations ∼4.6 times as of 1 September 2023. The database stores information about geometrical features, e.g. length and radius, and physico-chemical properties based on channel-lining amino acids. The stored data are interlinked with the available UniProt mutation annotation data. ChannelsDB 2.0 provides an excellent resource for deep analysis of the role of biomacromolecular tunnels and pores. The database is available free of charge: https://channelsdb2.biodata.ceitec.cz.

AggreProt: a web server for predicting and engineering aggregation prone regions in proteins.

Recombinant proteins play pivotal roles in numerous applications including industrial biocatalysts or therapeutics. Despite the recent progress in computational protein structure prediction, protein solubility and reduced aggregation propensity remain challenging attributes to design. Identification of aggregation-prone regions is essential for understanding misfolding diseases or designing efficient protein-based technologies, and as such has a great socio-economic impact. Here, we introduce AggreProt, a user-friendly webserver that automatically exploits an ensemble of deep neural networks to predict aggregation-prone regions (APRs) in protein sequences. Trained on experimentally evaluated hexapeptides, AggreProt compares to or outperforms state-of-the-art algorithms on two independent benchmark datasets. The server provides per-residue aggregation profiles along with information on solvent accessibility and transmembrane propensity within an intuitive interface with interactive sequence and structure viewers for comprehensive analysis. We demonstrate AggreProt efficacy in predicting differential aggregation behaviours in proteins on several use cases, which emphasize its potential for guiding protein engineering strategies towards decreased aggregation propensity and improved solubility. The webserver is freely available and accessible at https://loschmidt.chemi.muni.cz/aggreprot/.

FireProt 2.0: web-based platform for the fully automated design of thermostable proteins.

Thermostable proteins find their use in numerous biomedical and biotechnological applications. However, the computational design of stable proteins often results in single-point mutations with a limited effect on protein stability. However, the construction of stable multiple-point mutants can prove difficult due to the possibility of antagonistic effects between individual mutations. FireProt protocol enables the automated computational design of highly stable multiple-point mutants. FireProt 2.0 builds on top of the previously published FireProt web, retaining the original functionality and expanding it with several new stabilization strategies. FireProt 2.0 integrates the AlphaFold database and the homology modeling for structure prediction, enabling calculations starting from a sequence. Multiple-point designs are constructed using the Bron-Kerbosch algorithm minimizing the antagonistic effect between the individual mutations. Users can newly limit the FireProt calculation to a set of user-defined mutations, run a saturation mutagenesis of the whole protein or select rigidifying mutations based on B-factors. Evolution-based back-to-consensus strategy is complemented by ancestral sequence reconstruction. FireProt 2.0 is significantly faster and a reworked graphical user interface broadens the tool's availability even to users with older hardware. FireProt 2.0 is freely available at http://loschmidt.chemi.muni.cz/fireprotweb.

PredictONCO: a web tool supporting decision-making in precision oncology by extending the bioinformatics predictions with advanced computing and machine learning.

PredictONCO 1.0 is a unique web server that analyzes effects of mutations on proteins frequently altered in various cancer types. The server can assess the impact of mutations on the protein sequential and structural properties and apply a virtual screening to identify potential inhibitors that could be used as a highly individualized therapeutic approach, possibly based on the drug repurposing. PredictONCO integrates predictive algorithms and state-of-the-art computational tools combined with information from established databases. The user interface was carefully designed for the target specialists in precision oncology, molecular pathology, clinical genetics and clinical sciences. The tool summarizes the effect of the mutation on protein stability and function and currently covers 44 common oncological targets. The binding affinities of Food and Drug Administration/ European Medicines Agency -approved drugs with the wild-type and mutant proteins are calculated to facilitate treatment decisions. The reliability of predictions was confirmed against 108 clinically validated mutations. The server provides a fast and compact output, ideal for the often time-sensitive decision-making process in oncology. Three use cases of missense mutations, (i) K22A in cyclin-dependent kinase 4 identified in melanoma, (ii) E1197K mutation in anaplastic lymphoma kinase 4 identified in lung carcinoma and (iii) V765A mutation in epidermal growth factor receptor in a patient with congenital mismatch repair deficiency highlight how the tool can increase levels of confidence regarding the pathogenicity of the variants and identify the most effective inhibitors. The server is available at https://loschmidt.chemi.muni.cz/predictonco.

pyCaverDock: Python implementation of the popular tool for analysis of ligand transport with advanced caching and batch calculation support.

SUMMARY: Access pathways in enzymes are crucial for the passage of substrates and products of catalysed reactions. The process can be studied by computational means with variable degrees of precision. Our in-house approximative method CaverDock provides a fast and easy way to set up and run ligand binding and unbinding calculations through protein tunnels and channels. Here we introduce pyCaverDock, a Python3 API designed to improve user experience with the tool and further facilitate the ligand transport analyses. The API enables users to simplify the steps needed to use CaverDock, from automatizing setup processes to designing screening pipelines. AVAILABILITY AND IMPLEMENTATION: pyCaverDock API is implemented in Python 3 and is freely available with detailed documentation and practical examples at https://loschmidt.chemi.muni.cz/caverdock/.

LoopGrafter: a web tool for transplanting dynamical loops for protein engineering.

The transplantation of loops between structurally related proteins is a compelling method to improve the activity, specificity and stability of enzymes. However, despite the interest of loop regions in protein engineering, the available methods of loop-based rational protein design are scarce. One particular difficulty related to loop engineering is the unique dynamism that enables them to exert allosteric control over the catalytic function of enzymes. Thus, when engaging in a transplantation effort, such dynamics in the context of protein structure need consideration. A second practical challenge is identifying successful excision points for the transplantation or grafting. Here, we present LoopGrafter (https://loschmidt.chemi.muni.cz/loopgrafter/), a web server that specifically guides in the loop grafting process between structurally related proteins. The server provides a step-by-step interactive procedure in which the user can successively identify loops in the two input proteins, calculate their geometries, assess their similarities and dynamics, and select a number of loops to be transplanted. All possible different chimeric proteins derived from any existing recombination point are calculated, and 3D models for each of them are constructed and energetically evaluated. The obtained results can be interactively visualized in a user-friendly graphical interface and downloaded for detailed structural analyses.

Virtual screening of potential anticancer drugs based on microbial products.

The development of microbial products for cancer treatment has been in the spotlight in recent years. In order to accelerate the lengthy and expensive drug development process, in silico screening tools are systematically employed, especially during the initial discovery phase. Moreover, considering the steadily increasing number of molecules approved by authorities for commercial use, there is a demand for faster methods to repurpose such drugs. Here we present a review on virtual screening web tools, such as publicly available databases of molecular targets and libraries of ligands, with the aim to facilitate the discovery of potential anticancer drugs based on microbial products. We provide an entry-level step-by-step description of the workflow for virtual screening of microbial metabolites with known protein targets, as well as two practical examples using freely available web tools. The first case presents a virtual screening study of drugs developed from microbial products using Caver Web, a web tool that performs docking along a tunnel. The second case comprises a comparative analysis between a wild type isocitrate dehydrogenase 1 and a mutant that results in cancer, using the recently developed web tool PredictSNPOnco. In summary, this review provides the basic and essential background information necessary for virtual screening experiments, which may accelerate the discovery of novel anticancer drugs.

FireProtASR: A Web Server for Fully Automated Ancestral Sequence Reconstruction.

There is a great interest in increasing proteins' stability to widen their usability in numerous biomedical and biotechnological applications. However, native proteins cannot usually withstand the harsh industrial environment, since they are evolved to function under mild conditions. Ancestral sequence reconstruction is a well-established method for deducing the evolutionary history of genes. Besides its applicability to discover the most probable evolutionary ancestors of the modern proteins, ancestral sequence reconstruction has proven to be a useful approach for the design of highly stable proteins. Recently, several computational tools were developed, which make the ancestral reconstruction algorithms accessible to the community, while leaving the most crucial steps of the preparation of the input data on users' side. FireProtASR aims to overcome this obstacle by constructing a fully automated workflow, allowing even the unexperienced users to obtain ancestral sequences based on a sequence query as the only input. FireProtASR is complemented with an interactive, easy-to-use web interface and is freely available at https://loschmidt.chemi.muni.cz/fireprotasr/.

FireProtDB: database of manually curated protein stability data.

The majority of naturally occurring proteins have evolved to function under mild conditions inside the living organisms. One of the critical obstacles for the use of proteins in biotechnological applications is their insufficient stability at elevated temperatures or in the presence of salts. Since experimental screening for stabilizing mutations is typically laborious and expensive, in silico predictors are often used for narrowing down the mutational landscape. The recent advances in machine learning and artificial intelligence further facilitate the development of such computational tools. However, the accuracy of these predictors strongly depends on the quality and amount of data used for training and testing, which have often been reported as the current bottleneck of the approach. To address this problem, we present a novel database of experimental thermostability data for single-point mutants FireProtDB. The database combines the published datasets, data extracted manually from the recent literature, and the data collected in our laboratory. Its user interface is designed to facilitate both types of the expected use: (i) the interactive explorations of individual entries on the level of a protein or mutation and (ii) the construction of highly customized and machine learning-friendly datasets using advanced searching and filtering. The database is freely available at https://loschmidt.chemi.muni.cz/fireprotdb.

Hidden Potential of Highly Efficient and Widely Accessible Thrombolytic Staphylokinase.

Stroke burden is substantially increasing but current therapeutic drugs are still far from ideal. Here we highlight the vast potential of staphylokinase as an efficient, fibrin-selective, inexpensive, and evolvable thrombolytic agent. The emphasis is escalated by new recent findings. Staphylokinase nonimmunogenic variant was proven noninferior to alteplase in a clinical trial, with decreased risk of intracranial hemorrhage and the advantage of single bolus administration. Furthermore, our detailed kinetic analysis revealed a new staphylokinase limiting bottleneck whose elimination might provide up to 1000-fold higher activity than the clinically approved alteplase. This knowledge of limitations unlocks new possibilities for improvements that are now achievable by the community of protein engineers who have the required expertise and are ready to transform staphylokinase into a powerful molecule. Collectively, the noninferiority and safety of nonimmunogenic staphylokinase together with the newly identified effectivity limitation make staphylokinase a perfect candidate for further exploration, modification, and advancement to make it the next-generation widely accessible thrombolytic drug effectively treating stroke all around the world, including middle- and low-income countries.

SoluProt: prediction of soluble protein expression in Escherichia coli.

MOTIVATION: Poor protein solubility hinders the production of many therapeutic and industrially useful proteins. Experimental efforts to increase solubility are plagued by low success rates and often reduce biological activity. Computational prediction of protein expressibility and solubility in Escherichia coli using only sequence information could reduce the cost of experimental studies by enabling prioritization of highly soluble proteins. RESULTS: A new tool for sequence-based prediction of soluble protein expression in E.coli, SoluProt, was created using the gradient boosting machine technique with the TargetTrack database as a training set. When evaluated against a balanced independent test set derived from the NESG database, SoluProt's accuracy of 58.5% and AUC of 0.62 exceeded those of a suite of alternative solubility prediction tools. There is also evidence that it could significantly increase the success rate of experimental protein studies. SoluProt is freely available as a standalone program and a user-friendly webserver at https://loschmidt.chemi.muni.cz/soluprot/. AVAILABILITY AND IMPLEMENTATION: https://loschmidt.chemi.muni.cz/soluprot/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Follistatin-like 1 and its paralogs in heart development and cardiovascular disease.

Cardiovascular diseases (CVDs) are a group of disorders affecting the heart and blood vessels and a leading cause of death worldwide. Thus, there is a need to identify new cardiokines that may protect the heart from damage as reported in GBD 2017 Causes of Death Collaborators (2018) (The Lancet 392:1736-1788). Follistatin-like 1 (FSTL1) is a cardiokine that is highly expressed in the heart and released to the serum after cardiac injury where it is associated with CVD and predicts poor outcome. The action of FSTL1 likely depends not only on the tissue source but also post-translation modifications that are target tissue- and cell-specific. Animal studies examining the effect of FSTL1 in various models of heart disease have exploded over the past 15 years and primarily report a protective effect spanning from inhibiting inflammation via transforming growth factor, preventing remodeling and fibrosis to promoting angiogenesis and hypertrophy. A better understanding of FSTL1 and its homologs is needed to determine whether this protein could be a useful novel biomarker to predict poor outcome and death and whether it has therapeutic potential. The aim of this review is to provide a comprehensive description of the literature for this family of proteins in order to better understand their role in normal physiology and CVD.

PepSeA: Peptide Sequence Alignment and Visualization Tools to Enable Lead Optimization.

Therapeutic peptides offer potential advantages over small molecules in terms of selectivity, affinity, and their ability to target "undruggable" proteins that are associated with a wide range of pathologies. Despite their importance, current molecular design capabilities that inform medicinal chemistry decisions on peptide programs are limited. More specifically, there are unmet needs for structure-activity relationship (SAR) analysis and visualization of linear, cyclic, and cross-linked peptides containing non-natural motifs, which are widely used in drug discovery. To bridge this gap, we developed PepSeA (Peptide Sequence Alignment and Visualization), an open-source, freely available package of sequence-based tools (https://github.com/Merck/PepSeA). PepSeA enables multiple sequence alignment of non-natural amino acids and enhanced visualization with the hierarchical editing language for macromolecules (HELM). Via stepwise SAR analysis of a ChEMBL peptide data set, we demonstrate the utility of PepSeA to accelerate decision making in lead optimization campaigns in pharmaceutical setting. PepSeA represents an initial attempt to expand cheminformatics capabilities for therapeutic peptides and to enable rapid and more efficient design-make-test cycles.

EnzymeMiner: automated mining of soluble enzymes with diverse structures, catalytic properties and stabilities.

Millions of protein sequences are being discovered at an incredible pace, representing an inexhaustible source of biocatalysts. Despite genomic databases growing exponentially, classical biochemical characterization techniques are time-demanding, cost-ineffective and low-throughput. Therefore, computational methods are being developed to explore the unmapped sequence space efficiently. Selection of putative enzymes for biochemical characterization based on rational and robust analysis of all available sequences remains an unsolved problem. To address this challenge, we have developed EnzymeMiner-a web server for automated screening and annotation of diverse family members that enables selection of hits for wet-lab experiments. EnzymeMiner prioritizes sequences that are more likely to preserve the catalytic activity and are heterologously expressible in a soluble form in Escherichia coli. The solubility prediction employs the in-house SoluProt predictor developed using machine learning. EnzymeMiner reduces the time devoted to data gathering, multi-step analysis, sequence prioritization and selection from days to hours. The successful use case for the haloalkane dehalogenase family is described in a comprehensive tutorial available on the EnzymeMiner web page. EnzymeMiner is a universal tool applicable to any enzyme family that provides an interactive and easy-to-use web interface freely available at https://loschmidt.chemi.muni.cz/enzymeminer/.

Caver Web 1.0: identification of tunnels and channels in proteins and analysis of ligand transport.

Caver Web 1.0 is a web server for comprehensive analysis of protein tunnels and channels, and study of the ligands' transport through these transport pathways. Caver Web is the first interactive tool allowing both the analyses within a single graphical user interface. The server is built on top of the abundantly used tunnel detection tool Caver 3.02 and CaverDock 1.0 enabling the study of the ligand transport. The program is easy-to-use as the only required inputs are a protein structure for a tunnel identification and a list of ligands for the transport analysis. The automated guidance procedures assist the users to set up the calculation in a way to obtain biologically relevant results. The identified tunnels, their properties, energy profiles and trajectories for ligands' passages can be calculated and visualized. The tool is very fast (2-20 min per job) and is applicable even for virtual screening purposes. Its simple setup and comprehensive graphical user interface make the tool accessible for a broad scientific community. The server is freely available at https://loschmidt.chemi.muni.cz/caverweb.

CaverDock: a molecular docking-based tool to analyse ligand transport through protein tunnels and channels.

MOTIVATION: Protein tunnels and channels are key transport pathways that allow ligands to pass between proteins' external and internal environments. These functionally important structural features warrant detailed attention. It is difficult to study the ligand binding and unbinding processes experimentally, while molecular dynamics simulations can be time-consuming and computationally demanding. RESULTS: CaverDock is a new software tool for analysing the ligand passage through the biomolecules. The method uses the optimized docking algorithm of AutoDock Vina for ligand placement docking and implements a parallel heuristic algorithm to search the space of possible trajectories. The duration of the simulations takes from minutes to a few hours. Here we describe the implementation of the method and demonstrate CaverDock's usability by: (i) comparison of the results with other available tools, (ii) determination of the robustness with large ensembles of ligands and (iii) the analysis and comparison of the ligand trajectories in engineered tunnels. Thorough testing confirms that CaverDock is applicable for the fast analysis of ligand binding and unbinding in fundamental enzymology and protein engineering. AVAILABILITY AND IMPLEMENTATION: User guide and binaries for Ubuntu are freely available for non-commercial use at https://loschmidt.chemi.muni.cz/caverdock/. The web implementation is available at https://loschmidt.chemi.muni.cz/caverweb/. The source code is available upon request. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A Haloalkane Dehalogenase from Saccharomonospora viridis Strain DSM 43017, a Compost Bacterium with Unusual Catalytic Residues, Unique (S)-Enantiopreference, and High Thermostability.

Haloalkane dehalogenases can cleave a carbon-halogen bond in a broad range of halogenated aliphatic compounds. However, a highly conserved catalytic pentad composed of a nucleophile, a catalytic base, a catalytic acid, and two halide-stabilizing residues is required for their catalytic activity. Only a few family members, e.g., DsaA, DmxA, or DmrB, remain catalytically active while employing a single halide-stabilizing residue. Here, we describe a novel haloalkane dehalogenase, DsvA, from a mildly thermophilic bacterium, Saccharomonospora viridis strain DSM 43017, possessing one canonical halide-stabilizing tryptophan (W125). At the position of the second halide-stabilizing residue, DsvA contains the phenylalanine F165, which cannot stabilize the halogen anion released during the enzymatic reaction by a hydrogen bond. Based on the sequence and structural alignments, we identified a putative second halide-stabilizing tryptophan (W162) located on the same α-helix as F165, but on the opposite side of the active site. The potential involvement of this residue in DsvA catalysis was investigated by the construction and biochemical characterization of the three variants, DsvA01 (F165W), DsvA02 (W162F), and DsvA03 (W162F and F165W). Interestingly, DsvA exhibits a preference for the (S)- over the (R)-enantiomers of ß-bromoalkanes, which has not been reported before for any characterized haloalkane dehalogenase. Moreover, DsvA shows remarkable operational stability at elevated temperatures. The present study illustrates that protein sequences possessing an unconventional composition of catalytic residues represent a valuable source of novel biocatalysts.IMPORTANCE The present study describes a novel haloalkane dehalogenase, DsvA, originating from a mildly thermophilic bacterium, Saccharomonospora viridis strain DSM 43017. We report its high thermostability, remarkable operational stability at high temperatures, and an (S)-enantiopreference, which makes this enzyme an attractive biocatalyst for practical applications. Sequence analysis revealed that DsvA possesses an unusual composition of halide-stabilizing tryptophan residues in its active site. We constructed and biochemically characterized two single point mutants and one double point mutant and identified the noncanonical halide-stabilizing residue. Our study underlines the importance of searching for noncanonical catalytic residues in protein sequences.

HotSpot Wizard 3.0: web server for automated design of mutations and smart libraries based on sequence input information.

HotSpot Wizard is a web server used for the automated identification of hotspots in semi-rational protein design to give improved protein stability, catalytic activity, substrate specificity and enantioselectivity. Since there are three orders of magnitude fewer protein structures than sequences in bioinformatic databases, the major limitation to the usability of previous versions was the requirement for the protein structure to be a compulsory input for the calculation. HotSpot Wizard 3.0 now accepts the protein sequence as input data. The protein structure for the query sequence is obtained either from eight repositories of homology models or is modeled using Modeller and I-Tasser. The quality of the models is then evaluated using three quality assessment tools-WHAT_CHECK, PROCHECK and MolProbity. During follow-up analyses, the system automatically warns the users whenever they attempt to redesign poorly predicted parts of their homology models. The second main limitation of HotSpot Wizard's predictions is that it identifies suitable positions for mutagenesis, but does not provide any reliable advice on particular substitutions. A new module for the estimation of thermodynamic stabilities using the Rosetta and FoldX suites has been introduced which prevents destabilizing mutations among pre-selected variants entering experimental testing. HotSpot Wizard is freely available at http://loschmidt.chemi.muni.cz/hotspotwizard.

CalFitter: a web server for analysis of protein thermal denaturation data.

Despite significant advances in the understanding of protein structure-function relationships, revealing protein folding pathways still poses a challenge due to a limited number of relevant experimental tools. Widely-used experimental techniques, such as calorimetry or spectroscopy, critically depend on a proper data analysis. Currently, there are only separate data analysis tools available for each type of experiment with a limited model selection. To address this problem, we have developed the CalFitter web server to be a unified platform for comprehensive data fitting and analysis of protein thermal denaturation data. The server allows simultaneous global data fitting using any combination of input data types and offers 12 protein unfolding pathway models for selection, including irreversible transitions often missing from other tools. The data fitting produces optimal parameter values, their confidence intervals, and statistical information to define unfolding pathways. The server provides an interactive and easy-to-use interface that allows users to directly analyse input datasets and simulate modelled output based on the model parameters. CalFitter web server is available free at https://loschmidt.chemi.muni.cz/calfitter/.

Conformational changes allow processing of bulky substrates by a haloalkane dehalogenase with a small and buried active site.

Haloalkane dehalogenases catalyze the hydrolysis of halogen-carbon bonds in organic halogenated compounds and as such are of great utility as biocatalysts. The crystal structures of the haloalkane dehalogenase DhlA from the bacterium from Xanthobacter autotrophicus GJ10, specifically adapted for the conversion of the small 1,2-dichloroethane (DCE) molecule, display the smallest catalytic site (110 Å3) within this enzyme family. However, during a substrate-specificity screening, we noted that DhlA can catalyze the conversion of far bulkier substrates, such as the 4-(bromomethyl)-6,7-dimethoxy-coumarin (220 Å3). This large substrate cannot bind to DhlA without conformational alterations. These conformational changes have been previously inferred from kinetic analysis, but their structural basis has not been understood. Using molecular dynamic simulations, we demonstrate here the intrinsic flexibility of part of the cap domain that allows DhlA to accommodate bulky substrates. The simulations displayed two routes for transport of substrates to the active site, one of which requires the conformational change and is likely the route for bulky substrates. These results provide insights into the structure-dynamics function relationships in enzymes with deeply buried active sites. Moreover, understanding the structural basis for the molecular adaptation of DhlA to 1,2-dichloroethane introduced into the biosphere during the industrial revolution provides a valuable lesson in enzyme design by nature.