Evolution of tunnels in α/ß-hydrolase fold proteins-What can we learn from studying epoxide hydrolases?

The evolutionary variability of a protein's residues is highly dependent on protein region and function. Solvent-exposed residues, excluding those at interaction interfaces, are more variable than buried residues whereas active site residues are considered to be conserved. The abovementioned rules apply also to α/ß-hydrolase fold proteins-one of the oldest and the biggest superfamily of enzymes with buried active sites equipped with tunnels linking the reaction site with the exterior. We selected soluble epoxide hydrolases as representative of this family to conduct the first systematic study on the evolution of tunnels. We hypothesised that tunnels are lined by mostly conserved residues, and are equipped with a number of specific variable residues that are able to respond to evolutionary pressure. The hypothesis was confirmed, and we suggested a general and detailed way of the tunnels' evolution analysis based on entropy values calculated for tunnels' residues. We also found three different cases of entropy distribution among tunnel-lining residues. These observations can be applied for protein reengineering mimicking the natural evolution process. We propose a 'perforation' mechanism for new tunnels design via the merging of internal cavities or protein surface perforation. Based on the literature data, such a strategy of new tunnel design could significantly improve the enzyme's performance and can be applied widely for enzymes with buried active sites.

Recent Advances in Studying Toll-like Receptors with the Use of Computational Methods.

Toll-like receptors (TLRs) are transmembrane proteins that recognize various molecular patterns and activate signaling that triggers the immune response. In this review, our goal is to summarize how, in recent years, various computational solutions have contributed to a better understanding of TLRs, regarding both their function and mechanism of action. We update the recent information about small-molecule modulators and expanded the topic toward next-generation vaccine design, as well as studies of the dynamic nature of TLRs. Also, we underline problems that remain unsolved.

Geometry-Based versus Small-Molecule Tracking Method for Tunnel Identification: Benefits and Pitfalls.

Different methods for tunnel identification, geometry-based and small-molecule tracking approaches, were compared to provide their benefits and pitfalls. Results obtained for both crystal structures and molecular dynamics (MD) simulations were analyzed to investigate if a more computationally demanding method would be beneficial. Careful examination of the results is essential for the low-diameter tunnel description, and assessment of the tunnel functionality based only on their geometrical parameters is challenging. We showed that the small-molecule tracking approach can provide a detailed description of the system; however, it can also be the most computationally demanding.

AQUA-DUCT 1.0: structural and functional analysis of macromolecules from an intramolecular voids perspective.

MOTIVATION: Tunnels, pores, channels, pockets and cavities contribute to proteins architecture and performance. However, analysis and characteristics of transportation pathways and internal binding cavities are performed separately. We aimed to provide universal tool for analysis of proteins integral interior with access to detailed information on the ligands transportation phenomena and binding preferences. RESULTS: AQUA-DUCT version 1.0 is a comprehensive method for macromolecules analysis from the intramolecular voids perspective using small ligands as molecular probes. This version gives insight into several properties of macromolecules and facilitates protein engineering and drug design by the combination of the tracking and local mapping approach to small ligands. AVAILABILITY AND IMPLEMENTATION: http://www.aquaduct.pl. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Computational Selectivity Assessment of Protease Inhibitors against SARS-CoV-2.

The pandemic of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) poses a serious global health threat. Since no specific therapeutics are available, researchers around the world screened compounds to inhibit various molecular targets of SARS-CoV-2 including its main protease (Mpro) essential for viral replication. Due to the high urgency of these discovery efforts, off-target binding, which is one of the major reasons for drug-induced toxicity and safety-related drug attrition, was neglected. Here, we used molecular docking, toxicity profiling, and multiple molecular dynamics (MD) protocols to assess the selectivity of 33 reported non-covalent inhibitors of SARS-CoV-2 Mpro against eight proteases and 16 anti-targets. The panel of proteases included SARS-CoV Mpro, cathepsin G, caspase-3, ubiquitin carboxy-terminal hydrolase L1 (UCHL1), thrombin, factor Xa, chymase, and prostasin. Several of the assessed compounds presented considerable off-target binding towards the panel of proteases, as well as the selected anti-targets. Our results further suggest a high risk of off-target binding to chymase and cathepsin G. Thus, in future discovery projects, experimental selectivity assessment should be directed toward these proteases. A systematic selectivity assessment of SARS-CoV-2 Mpro inhibitors, as we report it, was not previously conducted.

Simple Selection Procedure to Distinguish between Static and Flexible Loops.

Loops are the most variable and unorganized elements of the secondary structure of proteins. Their ability to shift their shape can play a role in the binding of small ligands, enzymatic catalysis, or protein-protein interactions. Due to the loop flexibility, the positions of their residues in solved structures show the largest B-factors, or in a worst-case scenario can be unknown. Based on the loops' movements' timeline, they can be divided into slow (static) and fast (flexible). Although most of the loops that are missing in experimental structures belong to the flexible loops group, the computational tools for loop reconstruction use a set of static loop conformations to predict the missing part of the structure and evaluate the model. We believe that these two loop types can adopt different conformations and that using scoring functions appropriate for static loops is not sufficient for flexible loops. We showed that common model evaluation methods, are insufficient in the case of flexible solvent-exposed loops. Instead, we recommend using the potential energy to evaluate such loop models. We provide a novel model selection method based on a set of geometrical parameters to distinguish between flexible and static loops without the use of molecular dynamics simulations. We have also pointed out the importance of water network and interactions with the solvent for the flexible loop modeling.

Structural and Evolutionary Analysis Indicate That the SARS-CoV-2 Mpro Is a Challenging Target for Small-Molecule Inhibitor Design.

The novel coronavirus whose outbreak took place in December 2019 continues to spread at a rapid rate worldwide. In the absence of an effective vaccine, inhibitor repurposing or de novo drug design may offer a longer-term strategy to combat this and future infections due to similar viruses. Here, we report on detailed classical and mixed-solvent molecular dynamics simulations of the main protease (Mpro) enriched by evolutionary and stability analysis of the protein. The results were compared with those for a highly similar severe acute respiratory syndrome (SARS) Mpro protein. In spite of a high level of sequence similarity, the active sites in both proteins showed major differences in both shape and size, indicating that repurposing SARS drugs for COVID-19 may be futile. Furthermore, analysis of the binding site's conformational changes during the simulation time indicated its flexibility and plasticity, which dashes hopes for rapid and reliable drug design. Conversely, structural stability of the protein with respect to flexible loop mutations indicated that the virus' mutability will pose a further challenge to the rational design of small-molecule inhibitors. However, few residues contribute significantly to the protein stability and thus can be considered as key anchoring residues for Mpro inhibitor design.

BALCONY: an R package for MSA and functional compartments of protein variability analysis.

BACKGROUND: Here, we present an R package for entropy/variability analysis that facilitates prompt and convenient data extraction, manipulation and visualization of protein features from multiple sequence alignments. BALCONY can work with residues dispersed across a protein sequence and map them on the corresponding alignment of homologous protein sequences. Additionally, it provides several entropy and variability scores that indicate the conservation of each residue. RESULTS: Our package allows the user to visualize evolutionary variability by locating the positions most likely to vary and to assess mutation candidates in protein engineering. CONCLUSION: In comparison to other R packages BALCONY allows conservation/variability analysis in context of protein structure with linkage of the appropriate metrics with physicochemical features of user choice. AVAILABILITY: CRAN project page: https://cran.r-project.org/package=BALCONY and our website: http://www.tunnelinggroup.pl/software/ for major platforms: Linux/Unix, Windows and Mac OS X.

AQUA-DUCT: a ligands tracking tool.

MOTIVATION: The identification and tracking of molecules which enter active site cavity requires screening the positions of thousands of single molecules along several thousand molecular dynamic steps. To fill the existing gap between tools searching for tunnels and pathways and advanced tools employed for accelerated water flux investigations, we have developed AQUA-DUCT. RESULTS: AQUA-DUCT is an easy-to-use tool that facilitates analysis of the behaviour of molecules that penetrate any selected region in a protein. It can be used for any type of molecules, e.g. water, oxygen, carbon dioxide, organic solvents, ions. AVAILABILITY AND IMPLEMENTATION: Linux, Windows, macOS, OpenBSD, http://www.aquaduct.pl . CONTACT: a.gora@tunnelinggroup.pl or info@aquaduct.pl. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

CAVER Analyst 1.0: graphic tool for interactive visualization and analysis of tunnels and channels in protein structures.

UNLABELLED: The transport of ligands, ions or solvent molecules into proteins with buried binding sites or through the membrane is enabled by protein tunnels and channels. CAVER Analyst is a software tool for calculation, analysis and real-time visualization of access tunnels and channels in static and dynamic protein structures. It provides an intuitive graphic user interface for setting up the calculation and interactive exploration of identified tunnels/channels and their characteristics. AVAILABILITY AND IMPLEMENTATION: CAVER Analyst is a multi-platform software written in JAVA. Binaries and documentation are freely available for non-commercial use at http://www.caver.cz.

Effects of γ-polyglutamic acid on grassland sandy soil properties and plant functional traits exposed to drought stress.

The current study provides field experimental data that support the use of γ-polyglutamic acid (γ-PGA) in drought stress and proposes its application in grassland management. We hypothesized that water treatment combined with PGA application to sandy soil would reduce drought stress in grasslands more effectively than watering alone. A randomized block design was used, with three replicate watering blocks (no watering, weekly watering, and monthly watering) and PGA treatments at four different concentrations (0%, 0.3%, 1%, and 2% PGA). The results showed that PGA acts as a biostimulant, alleviating the effects of stress in plants by: (1) increasing the availability of ions, especially K+, Zn2+, Mn2+, Fe2+/3+, Ca2+, and Mg2+, as well as N-NH4+, and N-NO3-, (2) elongating plant roots, (3) increasing the aboveground biomass, (4) improving the resprouting capacity of the dominant grass Nardus stricta, and (5) improving the regeneration of dicotyledons. In the case of meadows on sandy soils, the use of low PGA concentrations (0.3% or 1%) was the most beneficial for the availability of macro- and microelements and improving the functional traits of plants. Irrigation had a greater effect than using PGA only for the dicotyledon to monocotyledon ratio.

PBP-A, a cyanobacterial DD-peptidase with high specificity for amidated muropeptides, exhibits pH-dependent promiscuous activity harmful to Escherichia coli.

Penicillin binding proteins (PBPs) are involved in biosynthesis, remodeling and recycling of peptidoglycan (PG) in bacteria. PBP-A from Thermosynechococcus elongatus belongs to a cyanobacterial family of enzymes sharing close structural and phylogenetic proximity to class A ß-lactamases. With the long-term aim of converting PBP-A into a ß-lactamase by directed evolution, we simulated what may happen when an organism like Escherichia coli acquires such a new PBP and observed growth defect associated with the enzyme activity. To further explore the molecular origins of this harmful effect, we decided to characterize deeper the activity of PBP-A both in vitro and in vivo. We found that PBP-A is an enzyme endowed with DD-carboxypeptidase and DD-endopeptidase activities, featuring high specificity towards muropeptides amidated on the D-iso-glutamyl residue. We also show that a low promiscuous activity on non-amidated peptidoglycan deteriorates E. coli's envelope, which is much higher under acidic conditions where substrate discrimination is mitigated. Besides expanding our knowledge of the biochemical activity of PBP-A, this work also highlights that promiscuity may depend on environmental conditions and how it may hinder rather than promote enzyme evolution in nature or in the laboratory.

A single mutation in a tunnel to the active site changes the mechanism and kinetics of product release in haloalkane dehalogenase LinB.

Many enzymes have buried active sites. The properties of the tunnels connecting the active site with bulk solvent affect ligand binding and unbinding and also the catalytic properties. Here, we investigate ligand passage in the haloalkane dehalogenase enzyme LinB and the effect of replacing leucine by a bulky tryptophan at a tunnel-lining position. Transient kinetic experiments show that the mutation significantly slows down the rate of product release. Moreover, the mechanism of bromide ion release is changed from a one-step process in the wild type enzyme to a two-step process in the mutant. The rate constant of bromide ion release corresponds to the overall steady-state turnover rate constant, suggesting that product release became the rate-limiting step of catalysis in the mutant. We explain the experimental findings by investigating the molecular details of the process computationally. Analysis of trajectories from molecular dynamics simulations with a tunnel detection software reveals differences in the tunnels available for ligand egress. Corresponding differences are seen in simulations of product egress using a specialized enhanced sampling technique. The differences in the free energy barriers for egress of a bromide ion obtained using potential of mean force calculations are in good agreement with the differences in rates obtained from the transient kinetic experiments. Interactions of the bromide ion with the introduced tryptophan are shown to affect the free energy barrier for its passage. The study demonstrates how the mechanism of an enzymatic catalytic cycle and reaction kinetics can be engineered by modification of protein tunnels.

CAVER 3.0: a tool for the analysis of transport pathways in dynamic protein structures.

Tunnels and channels facilitate the transport of small molecules, ions and water solvent in a large variety of proteins. Characteristics of individual transport pathways, including their geometry, physico-chemical properties and dynamics are instrumental for understanding of structure-function relationships of these proteins, for the design of new inhibitors and construction of improved biocatalysts. CAVER is a software tool widely used for the identification and characterization of transport pathways in static macromolecular structures. Herein we present a new version of CAVER enabling automatic analysis of tunnels and channels in large ensembles of protein conformations. CAVER 3.0 implements new algorithms for the calculation and clustering of pathways. A trajectory from a molecular dynamics simulation serves as the typical input, while detailed characteristics and summary statistics of the time evolution of individual pathways are provided in the outputs. To illustrate the capabilities of CAVER 3.0, the tool was applied for the analysis of molecular dynamics simulation of the microbial enzyme haloalkane dehalogenase DhaA. CAVER 3.0 safely identified and reliably estimated the importance of all previously published DhaA tunnels, including the tunnels closed in DhaA crystal structures. Obtained results clearly demonstrate that analysis of molecular dynamics simulation is essential for the estimation of pathway characteristics and elucidation of the structural basis of the tunnel gating. CAVER 3.0 paves the way for the study of important biochemical phenomena in the area of molecular transport, molecular recognition and enzymatic catalysis. The software is freely available as a multiplatform command-line application at http://www.caver.cz.

Transient binding sites at the surface of haloalkane dehalogenase LinB as locations for fine-tuning enzymatic activity.

The haloalkane dehalogenase LinB is a well-known enzyme that contains buried active site and is used for many modelling studies. Using classical molecular dynamics simulations of enzymes and substrates, we searched for transient binding sites on the surface of the LinB protein by calculating maps of enzyme-ligand interactions that were then transformed into sparse matrices. All residues considered as functionally important for enzyme performance (e.g., tunnel entrances) were excluded from the analysis to concentrate rather on non-obvious surface residues. From a set of 130 surface residues, twenty-six were proposed as a promising improvement of enzyme performance. Eventually, based on rational selection and filtering out the potentially unstable mutants, a small library of ten mutants was proposed to validate the possibility of fine-tuning the LinB protein. Nearly half of the predicted mutant structures showed improved activity towards the selected substrates, which demonstrates that the proposed approach could be applied to identify non-obvious yet beneficial mutations for enzyme performance especially when obvious locations have already been explored.

Corrigendum to "Structure-function relationship between soluble epoxide hydrolases structure and their tunnel network" [Comput. Struct. Biotechnol. J. 20 (2022) 193-205].

[This corrects the article DOI: 10.1016/j.csbj.2021.10.042.].

Structure-function relationship between soluble epoxide hydrolases structure and their tunnel network.

Enzymes with buried active sites maintain their catalytic function via a single tunnel or tunnel network. In this study we analyzed the functionality of soluble epoxide hydrolases (sEHs) tunnel network, by comparing the overall enzyme structure with the tunnel's shape and size. sEHs were divided into three groups based on their structure and the tunnel usage. The obtained results were compared with known substrate preferences of the studied enzymes, as well as reported in our other work evolutionary analyses data. The tunnel network architecture corresponded well with the evolutionary lineage of the source organism and large differences between enzymes were observed from long fragments insertions. This strategy can be used during protein re-engineering process for large changes introduction, whereas tunnel modification can be applied for fine-tuning of enzyme.

Structural Analysis of the Effect of Asn107Ser Mutation on Alg13 Activity and Alg13-Alg14 Complex Formation and Expanding the Phenotypic Variability of ALG13-CDG.

Congenital Disorders of Glycosylation (CDG) are multisystemic metabolic disorders showing highly heterogeneous clinical presentation, molecular etiology, and laboratory results. Here, we present different transferrin isoform patterns (obtained by isoelectric focusing) from three female patients harboring the ALG13 c.320A>G mutation. Contrary to other known variants of type I CDGs, where transferrin isoelectric focusing revealed notably increased asialo- and disialotransferrin fractions, a normal glycosylation pattern was observed in the probands. To verify this data and give novel insight into this variant, we modeled the human Alg13 protein and analyzed the dynamics of the apo structure and the complex with the UDP-GlcNAc substrate. We also modeled the Alg13-Alg14 heterodimer and ran multiple simulations of the complex in the presence of the substrate. Finally, we proposed a plausible complex formation mechanism.