In Silico-Directed Evolution Using CADEE.

Recent years have seen an explosion of interest in both sequence- and structure-based approaches toward in silico-directed evolution. We recently developed a novel computational toolkit, CADEE, which facilitates the computer-aided directed evolution of enzymes. Our initial work (Amrein et al., IUCrJ 4:50-64, 2017) presented a pedagogical example of the application of CADEE to triosephosphate isomerase, to illustrate the CADEE workflow. In this contribution, we describe this workflow in detail, including code input/output snippets, in order to allow users to set up and execute CADEE simulations on any system of interest.

CADEE: Computer-Aided Directed Evolution of Enzymes.

The tremendous interest in enzymes as biocatalysts has led to extensive work in enzyme engineering, as well as associated methodology development. Here, a new framework for computer-aided directed evolution of enzymes (CADEE) is presented which allows a drastic reduction in the time necessary to prepare and analyze in silico semi-automated directed evolution of enzymes. A pedagogical example of the application of CADEE to a real biological system is also presented in order to illustrate the CADEE workflow.

Expanding the Catalytic Triad in Epoxide Hydrolases and Related Enzymes.

Potato epoxide hydrolase 1 exhibits rich enantio- and regioselectivity in the hydrolysis of a broad range of substrates. The enzyme can be engineered to increase the yield of optically pure products as a result of changes in both enantio- and regioselectivity. It is thus highly attractive in biocatalysis, particularly for the generation of enantiopure fine chemicals and pharmaceuticals. The present work aims to establish the principles underlying the activity and selectivity of the enzyme through a combined computational, structural, and kinetic study using the substrate trans-stilbene oxide as a model system. Extensive empirical valence bond simulations have been performed on the wild-type enzyme together with several experimentally characterized mutants. We are able to computationally reproduce the differences between the activities of different stereoisomers of the substrate and the effects of mutations of several active-site residues. In addition, our results indicate the involvement of a previously neglected residue, H104, which is electrostatically linked to the general base H300. We find that this residue, which is highly conserved in epoxide hydrolases and related hydrolytic enzymes, needs to be in its protonated form in order to provide charge balance in an otherwise negatively charged active site. Our data show that unless the active-site charge balance is correctly treated in simulations, it is not possible to generate a physically meaningful model for the enzyme that can accurately reproduce activity and selectivity trends. We also expand our understanding of other catalytic residues, demonstrating in particular the role of a noncanonical residue, E35, as a "backup base" in the absence of H300. Our results provide a detailed view of the main factors driving catalysis and regioselectivity in this enzyme and identify targets for subsequent enzyme design efforts.

Recent advances in QM/MM free energy calculations using reference potentials.

BACKGROUND: Recent years have seen enormous progress in the development of methods for modeling (bio)molecular systems. This has allowed for the simulation of ever larger and more complex systems. However, as such complexity increases, the requirements needed for these models to be accurate and physically meaningful become more and more difficult to fulfill. The use of simplified models to describe complex biological systems has long been shown to be an effective way to overcome some of the limitations associated with this computational cost in a rational way. SCOPE OF REVIEW: Hybrid QM/MM approaches have rapidly become one of the most popular computational tools for studying chemical reactivity in biomolecular systems. However, the high cost involved in performing high-level QM calculations has limited the applicability of these approaches when calculating free energies of chemical processes. In this review, we present some of the advances in using reference potentials and mean field approximations to accelerate high-level QM/MM calculations. We present illustrative applications of these approaches and discuss challenges and future perspectives for the field. MAJOR CONCLUSIONS: The use of physically-based simplifications has shown to effectively reduce the cost of high-level QM/MM calculations. In particular, lower-level reference potentials enable one to reduce the cost of expensive free energy calculations, thus expanding the scope of problems that can be addressed. GENERAL SIGNIFICANCE: As was already demonstrated 40 years ago, the usage of simplified models still allows one to obtain cutting edge results with substantially reduced computational cost. This article is part of a Special Issue entitled Recent developments of molecular dynamics.

Force field independent metal parameters using a nonbonded dummy model.

The cationic dummy atom approach provides a powerful nonbonded description for a range of alkaline-earth and transition-metal centers, capturing both structural and electrostatic effects. In this work we refine existing literature parameters for octahedrally coordinated Mn(2+), Zn(2+), Mg(2+), and Ca(2+), as well as providing new parameters for Ni(2+), Co(2+), and Fe(2+). In all the cases, we are able to reproduce both M(2+)-O distances and experimental solvation free energies, which has not been achieved to date for transition metals using any other model. The parameters have also been tested using two different water models and show consistent performance. Therefore, our parameters are easily transferable to any force field that describes nonbonded interactions using Coulomb and Lennard-Jones potentials. Finally, we demonstrate the stability of our parameters in both the human and Escherichia coli variants of the enzyme glyoxalase I as showcase systems, as both enzymes are active with a range of transition metals. The parameters presented in this work provide a valuable resource for the molecular simulation community, as they extend the range of metal ions that can be studied using classical approaches, while also providing a starting point for subsequent parametrization of new metal centers.

Modeling catalytic promiscuity in the alkaline phosphatase superfamily.

In recent years, it has become increasingly clear that promiscuity plays a key role in the evolution of new enzyme function. This finding has helped to elucidate fundamental aspects of molecular evolution. While there has been extensive experimental work on enzyme promiscuity, computational modeling of the chemical details of such promiscuity has traditionally fallen behind the advances in experimental studies, not least due to the nearly prohibitive computational cost involved in examining multiple substrates with multiple potential mechanisms and binding modes in atomic detail with a reasonable degree of accuracy. However, recent advances in both computational methodologies and power have allowed us to reach a stage in the field where we can start to overcome this problem, and molecular simulations can now provide accurate and efficient descriptions of complex biological systems with substantially less computational cost. This has led to significant advances in our understanding of enzyme function and evolution in a broader sense. Here, we will discuss currently available computational approaches that can allow us to probe the underlying molecular basis for enzyme specificity and selectivity, discussing the inherent strengths and weaknesses of each approach. As a case study, we will discuss recent computational work on different members of the alkaline phosphatase superfamily (AP) using a range of different approaches, showing the complementary insights they have provided. We have selected this particular superfamily, as it poses a number of significant challenges for theory, ranging from the complexity of the actual reaction mechanisms involved to the reliable modeling of the catalytic metal centers, as well as the very large system sizes. We will demonstrate that, through current advances in methodologies, computational tools can provide significant insight into the molecular basis for catalytic promiscuity, and, therefore, in turn, the mechanisms of protein functional evolution.

PEGylation of lysine residues reduces the pro-migratory activity of IGF-I.

BACKGROUND: The insulin-like growth factor (IGF) system is composed of ligands and receptors which regulate cell proliferation, survival, differentiation and migration. Some of these functions involve regulation by the extracellular milieu, including binding proteins and other extracellular matrix proteins. However, the functions and exact nature of these interactions remain incomplete. METHODS: IGF-I variants PEGylated at lysines K27, K65 and K68, were assessed for binding to IGFBPs using BIAcore, and for phosphorylation of the IGF-IR. Furthermore, functional consequences of PEGylation were investigated using cell viability and migration assays. In addition, downstream signaling pathways were analyzed using phospho-AKT and phospho-ERK1/2 assays. RESULTS: IGF-I PEGylated at lysines 27 (PEG-K27), 65 (PEG-K65) or 68 (PEG-K68) was employed. Receptor phosphorylation was similarly reduced 2-fold with PEG-K65 and PEG-K68 in 3T3 fibroblasts and MCF-7 breast cancer cells, whereas PEG-K27 showed a more than 10- and 3-fold lower activation for 3T3 and MCF-7 cells, respectively. In addition, all PEG-IGF-I variants had a 10-fold reduced association rate to IGF binding proteins (IGFBPs). Functionally, all PEG variants lost their ability to induce cell migration in the presence of IGFBP-3/vitronectin (VN) complexes, whereas cell viability was fully preserved. Analysis of downstream signaling revealed that AKT was preferentially affected upon treatment with PEG-IGF-I variants whereas MAPK signaling was unaffected by PEGylation. CONCLUSION: PEGylation of IGF-I has an impact on cell migration but not on cell viability. GENERAL SIGNIFICANCE: PEG-IGF-I may differentially modulate IGF-I mediated functions that are dependent on receptor interaction as well as key extracellular proteins such as VN and IGFBPs.

Identification of two-histidines one-carboxylate binding motifs in proteins amenable to facial coordination to metals.

Among natural metalloenzymes, the facial two-histidines one-carboxylate binding motif (FTM) is a widely represented first coordination sphere motif present in the active site of a variety of metalloenzymes. A PDB search revealed a total of 1685 structures bearing such FTMs bound to a metal. Sixty statistically representative FTMs were selected and used as template for the identification of structurally characterized proteins bearing these three amino acids in a propitious environment for binding to a transition metal. This geometrical superposition search, carried out using the STAMPS software, returned 2320 hits. While most consisted of either apo-FTMs or bore strong sequence homology to known FTMs, seven such structures lying within a cavity were identified as novel and viable scaffolds for the creation of artificial metalloenzymes bearing an FTM.