Beerens K, Mazurenko S, Kunka A, Marques SM, Hansen N, Musil M, Chaloupkova R, Waterman J, Brezovsky J, Bednar D, Prokop Z, Damborsky J, 2018: Evolutionary Analysis is a Powerful Complement to Energy Calculations for Protein Stabilization. ACS Catalysis 8: 9420-9428. full text
Stability is one of the most important characteristics of proteins employed as biocatalysts, biotherapeutics, and biomaterials, and the role of computational approaches in modifying protein stability is rapidly expanding. We have recently identified stabilizing mutations in haloalkane dehalogenase DhaA using phylogenetic analysis but were not able to reproduce the effects of these mutations using force-field calculations. Here we tested four different hypotheses to explain the molecular basis of stabilization using structural, biochemical, biophysical, and computational analyses. We demonstrate that stabilization of DhaA by the mutations identified using the phylogenetic analysis is driven by both entropy and enthalpy contributions, in contrast to primarily enthalpy-driven stabilization by mutations designed by the force-field calculations. Comprehensive bioinformatics analysis revealed that more than half (53%) of 1 099 evolution-based stabilizing mutations would be evaluated as destabilizing by force-field calculations. Thermodynamic integration considers both folded and unfolded states and can describe the entropic component of stabilization, yet it is not suitable for predictive purposes due to its high computational demands. Altogether, our results strongly suggest that energetic calculations should be complemented by a phylogenetic analysis in protein-stabilization endeavors.
Jurcik A, Bednar D, Byska J, Marques SM, Furmanova K, Daniel L, Kokkonen P, Brezovsky J, Strnad O, Stourac J, Pavelka A, Manak M, Damborsky J, Kozlikova B, 2018: CAVER Analyst 2.0: Analysis and Visualization of Channels and Tunnels in Protein Structures and Molecular Dynamics Trajectories. Bioinformatics (just accepted): doi:10.1093/bioinformatics/bty386. full text
Studying the transport paths of ligands, solvents, or ions in transmembrane proteins and proteins with buried binding sites is fundamental to the understanding of their biological function. A detailed analysis of the structural features influencing the transport paths is also important for engineering proteins for biomedical and biotechnological applications.
CAVER Analyst 2.0 is a software tool for quantitative analysis and real-time visualization of tunnels and channels in static and dynamic structures. This version provides the users with many new functions, including advanced techniques for intuitive visual inspection of the spatiotemporal behavior of tunnels and channels. Novel integrated algorithms allow an efficient analysis and data reduction in large protein structures and molecular dynamic simulations.
CAVER Analyst 2.0 is a multi-platform standalone Java-based application. Binaries and documentation are freely available at www.caver.cz.
Kaushik S, Marques SM, Khirsariya P, Paruch K, Libichova L, Brezovsky J, Prokop Z, Chaloupkova R, Damborsky J, 2018: Impact of the Access Tunnel Engineering on Catalysis is Strictly Ligand-Specific.FEBS Journal (just accepted): doi: 10.1111/febs.14418. full text
The traditional way of rationally engineering enzymes to change their biocatalytic properties utilizes the modifications of their active sites. Another emerging approach is the engineering of structural features involved in the exchange of ligands between buried active sites and the surrounding solvent. However, surprisingly little is known about the effects of mutations that alter the access tunnels on the enzymes’ catalytic properties, and how these tunnels should be redesigned to allow fast passage of cognate substrates and products. Thus, we have systematically studied the effects of single-point mutations in a tunnel-lining residue of a haloalkane dehalogenase on the binding kinetics and catalytic conversion of both linear and branched haloalkanes. The hotspot residue Y176 was identified using computer simulations and randomized through saturation mutagenesis, and the resulting variants were screened for shifts in binding rates. Strikingly, opposite effects of the substituted residues on the catalytic efficiency towards linear and branched substrates were observed, which was found to be due to substrate-specific requirements in the critical steps of the respective catalytic cycles. We conclude that not only the catalytic sites but also the access pathways must be tailored specifically for each individual ligand, which is a new paradigm in protein engineering and de novo protein design. A rational approach is proposed here to address more effectively the task of designing ligand-specific tunnels using computational tools.
Vanacek P, Sebestova E, Babkova P, Bidmanova S, Daniel L, Dvorak P, Stepankova V, Chaloupkova R, Brezovsky J, Prokop Z, Damborsky J, 2018: Exploration of Enzyme Diversity by Integrating Bioinformatics with Expression Analysis and Biochemical Characterization.ACS Catalysis 8: 2402–2412. full text
Millions of protein sequences are being discovered at an incredible pace, representing an inexhaustible source of biocatalysts. Here, we describe an integrated system for automated in silico screening and systematic characterization of diverse family members. The workflow consists of: (i) identification and computational characterization of relevant genes by sequence/structural bioinformatics, (ii) expression analysis and activity screening of selected proteins, and (iii) complete biochemical/biophysical characterization, was validated against the haloalkane dehalogenase family. The sequence-based search identified 658 potential dehalogenases. The subsequent structural bioinformatics prioritized and selected 20 candidates for exploration of protein functional diversity. Out of these twenty, the expression analysis and the robotic screening of enzymatic activity provided 8 soluble proteins with dehalogenase activity. The enzymes discovered originated from genetically unrelated Bacteria and Eukaryota, and, for the first time, from Archaea. Overall, the integrated system provided biocatalysts with broad catalytic diversity showing unique substrate specificity profiles, covering a wide range of optimal operational temperature from 20 to 70 °C and an unusually broad pH range from 5.7 to 10. We obtained the most catalytically proficient native haloalkane dehalogenase enzyme to date (kcat/K0.5 = 96.8 mM-1.s-1), the most thermostable enzyme with melting temperature 71 °C, three different cold-adapted enzymes showing dehalogenase activity at near-to-zero temperatures and a biocatalyst degrading the warfare chemical sulfur mustard. The established strategy can be adapted to other enzyme families for exploration of their biocatalytic diversity in a large sequence space continuously growing due to the use of next-generation sequencing technologies.
Dvorak P, Bednar D, Vanacek P, Balek L, Eiselleova L, Stepankova V, Sebestova E, Kunova Bosakova M, Konecna Z, Mazurenko S, Kunka A, Vanova T, Zoufalova K, Chaloupkova R, Brezovsky J, Krejci P, Prokop Z, Dvorak P, Damborsky J, 2018: Computer-Assisted Engineering of Hyperstable Fibroblast Growth Factor 2.Biotechnology and Bioengineering (just accepted): doi: 10.1002/bit.26531. full text
Fibroblast growth factors (FGFs) serve numerous regulatory functions in complex organisms, and their corresponding therapeutic potential is of growing interest to academics and industrial researchers alike. However, applications of these proteins are limited due to their low stability. Here we tackle this problem using a generalizable computer-assisted protein engineering strategy to create a unique modified FGF2 with nine mutations displaying unprecedented stability and uncompromised biological function. The data from the characterization of stabilized FGF2 showed a remarkable prediction potential of in silico methods and provided insight into the unfolding mechanism of the protein. The molecule holds a considerable promise for stem cell research and medical or pharmaceutical applications.