Computational design of protein-ligand interfaces: potential in therapeutic development

Computational Redesign of an ω-Transaminase from Pseudomonas jessenii for Asymmetric Synthesis of Enantiopure Bulky Amines

ω-Transaminases (ω-TA) are attractive biocatalysts for the production of chiral amines from prochiral ketones via asymmetric synthesis. However, the substrate scope of ω-TAs is usually limited due to steric hindrance at the active site pockets. We explored a protein engineering strategy using computational design to expand the substrate scope of an (S)-selective ω-TA from Pseudomonas jessenii (PjTA-R6) toward the production of bulky amines. PjTA-R6 is attractive for use in applied biocatalysis due to its thermostability, tolerance to organic solvents, and acceptance of high concentrations of isopropylamine as amino donor. PjTA-R6 showed no detectable activity for the synthesis of six bicyclic or bulky amines targeted in this study. Six small libraries composed of 7–18 variants each were separately designed via computational methods and tested in the laboratory for ketone to amine conversion. In each library, the vast majority of the variants displayed the desired activity, and of the 40 different designs, 38 produced the target amine in good yield with >99% enantiomeric excess. This shows that the substrate scope and enantioselectivity of PjTA mutants could be predicted in silico with high accuracy. The single mutant W58G showed the best performance in the synthesis of five structurally similar bulky amines containing the indan and tetralin moieties. The best variant for the other bulky amine, 1-phenylbutylamine, was the triple mutant W58M + F86L + R417L, indicating that Trp58 is a key residue in the large binding pocket for PjTA-R6 redesign. Crystal structures of the two best variants confirmed the computationally predicted structures. The results show that computational design can be an efficient approach to rapidly expand the substrate scope of ω-TAs to produce enantiopure bulky amines.

Approaching protein design with multisite λ dynamics: Accurate and scalable mutational folding free energies in T4 lysozyme

The estimation of changes in free energy upon mutation is central to the problem of protein design. Modern protein design methods have had remarkable success over a wide range of design targets, but are reaching their limits in ligand binding and enzyme design due to insufficient accuracy in mutational free energies. Alchemical free energy calculations have the potential to supplement modern design methods through more accurate molecular dynamics based prediction of free energy changes, but suffer from high computational cost. Multisite λ dynamics (MSλD) is a particularly efficient and scalable free energy method with potential to explore combinatorially large sequence spaces inaccessible with other free energy methods. This work aims to quantify the accuracy of MSλD and demonstrate its scalability. We apply MSλD to the classic problem of calculating folding free energies in T4 lysozyme, a system with a wealth of experimental measurements. Single site mutants considering 32 mutations show remarkable agreement with experiment with a Pearson correlation of 0.914 and mean unsigned error of 1.19 kcal/mol. Multisite mutants in systems with up to five concurrent mutations spanning 240 different sequences show comparable agreement with experiment. These results demonstrate the promise of MSλD in exploring large sequence spaces for protein design.

An Evaluation of the Potential of NMR Spectroscopy and Computational Modelling Methods to Inform Biopharmaceutical Formulations

Protein-based therapeutics are considered to be one of the most important classes of pharmaceuticals on the market. The growing need to prolong stability of high protein concentrations in liquid form has proven to be challenging. Therefore, significant effort is being made to design formulations which can enable the storage of these highly concentrated protein therapies for up to 2 years. Currently, the excipient selection approach involves empirical high-throughput screening, but does not reveal details on aggregation mechanisms or the molecular-level effects of the formulations under storage conditions. Computational modelling approaches have the potential to elucidate such mechanisms, and rapidly screen in silico prior to experimental testing. Nuclear Magnetic Resonance (NMR) spectroscopy can also provide complementary insights into excipient–protein interactions. This review will highlight the underpinning principles of molecular modelling and NMR spectroscopy. It will also discuss the advancements in the applications of computational and NMR approaches in investigating excipient–protein interactions.

Insights from molecular dynamics simulations for computational protein design

A grand challenge in the field of structural biology is to design and engineer proteins that exhibit targeted functions. Although much success on this front has been achieved, design success rates remain low, an ever-present reminder of our limited understanding of the relationship between amino acid sequences and the structures they adopt. In addition to experimental techniques and rational design strategies, computational methods have been employed to aid in the design and engineering of proteins. Molecular dynamics (MD) is one such method that simulates the motions of proteins according to classical dynamics. Here, we review how insights into protein dynamics derived from MD simulations have influenced the design of proteins. One of the greatest strengths of MD is its capacity to reveal information beyond what is available in the static structures deposited in the Protein Data Bank. In this regard simulations can be used to directly guide protein design by providing atomistic details of the dynamic molecular interactions contributing to protein stability and function. MD simulations can also be used as a virtual screening tool to rank, select, identify, and assess potential designs. MD is uniquely poised to inform protein design efforts where the application requires realistic models of protein dynamics and atomic level descriptions of the relationship between dynamics and function. Here, we review cases where MD simulations was used to modulate protein stability and protein function by providing information regarding the conformation(s), conformational transitions, interactions, and dynamics that govern stability and function. In addition, we discuss cases where conformations from protein folding/unfolding simulations have been exploited for protein design, yielding novel outcomes that could not be obtained from static structures.

Designer protein delivery: From natural to engineered affinity-controlled release systems

Exploiting binding affinities between molecules is an established practice in many fields, including biochemical separations, diagnostics, and drug development; however, using these affinities to control biomolecule release is a more recent strategy. Affinity-controlled release takes advantage of the reversible nature of noncovalent interactions between a therapeutic protein and a binding partner to slow the diffusive release of the protein from a vehicle. This process, in contrast to degradation-controlled sustained-release formulations such as poly(lactic-co-glycolic acid) microspheres, is controlled through the strength of the binding interaction, the binding kinetics, and the concentration of binding partners. In the context of affinity-controlled release--and specifically the discovery or design of binding partners--we review advances in in vitro selection and directed evolution of proteins, peptides, and oligonucleotides (aptamers), aided by computational design.

Rosetta and the Design of Ligand Binding Sites

Proteins that bind small molecules (ligands) can be used as biosensors, signal modulators, and sequestering agents. When naturally occurring proteins for a particular target ligand are not available, artificial proteins can be computationally designed. We present a protocol based on RosettaLigand to redesign an existing protein pocket to bind a target ligand. Starting with a protein structure and the structure of the ligand, Rosetta can optimize both the placement of the ligand in the pocket and the identity and conformation of the surrounding sidechains, yielding proteins that bind the target compound.

Binding pocket optimization by computational protein design

Engineering specific interactions between proteins and small molecules is extremely useful for biological studies, as these interactions are essential for molecular recognition. Furthermore, many biotechnological applications are made possible by such an engineering approach, ranging from biosensors to the design of custom enzyme catalysts. Here, we present a novel method for the computational design of protein-small ligand binding named PocketOptimizer. The program can be used to modify protein binding pocket residues to improve or establish binding of a small molecule. It is a modular pipeline based on a number of customizable molecular modeling tools to predict mutations that alter the affinity of a target protein to its ligand. At its heart it uses a receptor-ligand scoring function to estimate the binding free energy between protein and ligand. We compiled a benchmark set that we used to systematically assess the performance of our method. It consists of proteins for which mutational variants with different binding affinities for their ligands and experimentally determined structures exist. Within this test set PocketOptimizer correctly predicts the mutant with the higher affinity in about 69% of the cases. A detailed analysis of the results reveals that the strengths of PocketOptimizer lie in the correct introduction of stabilizing hydrogen bonds to the ligand, as well as in the improved geometric complemetarity between ligand and binding pocket. Apart from the novel method for binding pocket design we also introduce a much needed benchmark data set for the comparison of affinities of mutant binding pockets, and that we use to asses programs for in silico design of ligand binding.

Two birds with one stone? Possible dual-targeting H1N1 inhibitors from traditional Chinese medicine

The H1N1 influenza pandemic of 2009 has claimed over 18,000 lives. During this pandemic, development of drug resistance further complicated efforts to control and treat the widespread illness. This research utilizes traditional Chinese medicine Database@Taiwan (TCM Database@Taiwan) to screen for compounds that simultaneously target H1 and N1 to overcome current difficulties with virus mutations. The top three candidates were de novo derivatives of xylopine and rosmaricine. Bioactivity of the de novo derivatives against N1 were validated by multiple machine learning prediction models. Ability of the de novo compounds to maintain CoMFA/CoMSIA contour and form key interactions implied bioactivity within H1 as well. Addition of a pyridinium fragment was critical to form stable interactions in H1 and N1 as supported by molecular dynamics (MD) simulation. Results from MD, hydrophobic interactions, and torsion angles are consistent and support the findings of docking. Multiple anchors and lack of binding to residues prone to mutation suggest that the TCM de novo derivatives may be resistant to drug resistance and are advantageous over conventional H1N1 treatments such as oseltamivir. These results suggest that the TCM de novo derivatives may be suitable candidates of dual-targeting drugs for influenza.

Determination of binding kinetics using whispering gallery mode microcavities

Silica optical microcavity sensors show great promise in the kinetic evaluation of binding pairs, fundamental in understanding biomolecular interactions. Here, we develop and demonstrate a novel platform, based on bioconjugated silica microsphere resonators, to study the binding kinetics of the biotin-streptavidin system. We characterize the optical performance, verify the covalent attachment of biotin to the surface, and perform streptavidin detection experiments. We perform preliminary kinetic analysis of the detection data which shows the potential of whispering gallery mode resonators in the determination of the dissociation constant of the binding pair, which is in good agreement with previously published values.