Design of protein-ligand binding based on the molecular-mechanics energy model

Ribose-Binding Protein Mutants With Improved Interaction Towards the Non-natural Ligand 1,3-Cyclohexanediol

Bioreporters consist of genetically modified living organisms that respond to the presence of target chemical compounds by production of an easily measurable signal. The central element in a bioreporter is a sensory protein or aptamer, which, upon ligand binding, modifies expression of the reporter signal protein. A variety of naturally occurring or modified versions of sensory elements has been exploited, but it has proven to be challenging to generate elements that recognize non-natural ligands. Bacterial periplasmic binding proteins have been proposed as a general scaffold to design receptor proteins for non-natural ligands, but despite various efforts, with only limited success. Here, we show how combinations of randomized mutagenesis and reporter screening improved the performance of a set of mutants in the ribose binding protein (RbsB) of Escherichia coli, which had been designed based on computational simulations to bind the non-natural ligand 1,3-cyclohexanediol (13CHD). Randomized mutant libraries were constructed that used the initially designed mutants as scaffolds, which were cloned in an appropriate E. coli bioreporter system and screened for improved induction of the GFPmut2 reporter fluorescence in presence of 1,3-cyclohexanediol. Multiple rounds of library screening, sorting, renewed mutagenesis and screening resulted in 4.5-fold improvement of the response to 1,3-cyclohexanediol and a lower detection limit of 0.25 mM. All observed mutations except one were located outside the direct ligand-binding pocket, suggesting they were compensatory and helping protein folding or functional behavior other than interaction with the ligand. Our results thus demonstrate that combinations of ligand-binding-pocket redesign and randomized mutagenesis can indeed lead to the selection and recovery of periplasmic-binding protein mutants with non-natural compound recognition. However, current lack of understanding of the intermolecular movement and ligand-binding in periplasmic binding proteins such as RbsB are limiting the rational production of further and better sensory mutants.

Toward the Accuracy and Speed of Protein Side-Chain Packing: A Systematic Study on Rotamer Libraries

Protein rotamers refer to the conformational isomers taken by the side-chains of amino acids to accommodate specific structural folding environments. Since accurate modeling of atomic interactions is difficult, rotamer information collected from experimentally solved protein structures is often used to guide side-chain packing in protein folding and sequence design studies. Many rotamer libraries have been built in the literature but there is little quantitative guidance on which libraries should be chosen for different structural modeling studies. Here, we performed a comparative study of six widely used rotamer libraries and systematically examined their suitability for protein folding and sequence design in four aspects: (1) side-chain match accuracy, (2) side-chain conformation prediction, (3) de novo protein sequence design, and (4) computational time cost. We demonstrated that, compared to the backbone-dependent rotamer libraries (BBDRLs), the backbone-independent rotamer libraries (BBIRLs) generated conformations that more closely matched the native conformations due to the larger number of rotamers in the local rotamer search spaces. However, more practically, using an optimized physical energy function incorporated into a simulated annealing Monte Carlo searching scheme, we showed that utilization of the BBDRLs could result in higher accuracies in side-chain prediction and higher sequence recapitulation rates in protein design experiments. Detailed data analyses showed that the major advantage of BBDRLs lies in the energy term derived from the rotamer probabilities that are associated with the individual backbone torsion angle subspaces. This term is important for distinguishing between amino acid identities as well as the rotamer conformations of an amino acid. Meanwhile, the backbone torsion angle subspace-specific rotamer search drastically speeds up the searching time, despite the significantly larger number of total rotamers in the BBDRLs. These results should provide important guidance for the development and selection of rotamer libraries for practical protein design and structure prediction studies.

Computational redesign of the Escherichia coli ribose-binding protein ligand binding pocket for 1,3-cyclohexanediol and cyclohexanol

Bacterial periplasmic-binding proteins have been acclaimed as general biosensing platform, but their range of natural ligands is too limited for optimal development of chemical compound detection. Computational redesign of the ligand-binding pocket of periplasmic-binding proteins may yield variants with new properties, but, despite earlier claims, genuine changes of specificity to non-natural ligands have so far not been achieved. In order to better understand the reasons of such limited success, we revisited here the Escherichia coli RbsB ribose-binding protein, aiming to achieve perceptible transition from ribose to structurally related chemical ligands 1,3-cyclohexanediol and cyclohexanol. Combinations of mutations were computationally predicted for nine residues in the RbsB binding pocket, then synthesized and tested in an E. coli reporter chassis. Two million variants were screened in a microcolony-in-bead fluorescence-assisted sorting procedure, which yielded six mutants no longer responsive to ribose but with 1.2-1.5 times induction in presence of 1 mM 1,3-cyclohexanediol, one of which responded to cyclohexanol as well. Isothermal microcalorimetry confirmed 1,3-cyclohexanediol binding, although only two mutant proteins were sufficiently stable upon purification. Circular dichroism spectroscopy indicated discernable structural differences between these two mutant proteins and wild-type RbsB. This and further quantification of periplasmic-space abundance suggested most mutants to be prone to misfolding and/or with defects in translocation compared to wild-type. Our results thus affirm that computational design and library screening can yield RbsB mutants with recognition of non-natural but structurally similar ligands. The inherent arisal of protein instability or misfolding concomitant with designed altered ligand-binding pockets should be overcome by new experimental strategies or by improved future protein design algorithms.

Rosetta:MSF: a modular framework for multi-state computational protein design

Computational protein design (CPD) is a powerful technique to engineer existing proteins or to design novel ones that display desired properties. Rosetta is a software suite including algorithms for computational modeling and analysis of protein structures and offers many elaborate protocols created to solve highly specific tasks of protein engineering. Most of Rosetta’s protocols optimize sequences based on a single conformation (i. e. design state). However, challenging CPD objectives like multi-specificity design or the concurrent consideration of positive and negative design goals demand the simultaneous assessment of multiple states. This is why we have developed the multi-state framework MSF that facilitates the implementation of Rosetta’s single-state protocols in a multi-state environment and made available two frequently used protocols. Utilizing MSF, we demonstrated for one of these protocols that multi-state design yields a 15% higher performance than single-state design on a ligand-binding benchmark consisting of structural conformations. With this protocol, we designed de novo nine retro-aldolases on a conformational ensemble deduced from a (βα)₈-barrel protein. All variants displayed measurable catalytic activity, testifying to a high success rate for this concept of multi-state enzyme design.

Protein engineering, i. e. the targeted modification or design of proteins has tremendous potential for medical and industrial applications. One generally applicable strategy for protein engineering is rational protein design: based on detailed knowledge of structure and function, computer programs like Rosetta propose the sequence of a protein possessing the desired properties. So far, most computer protocols have used rigid structures for design, which is a simplification because a protein’s structure is more accurately specified by a conformational ensemble. We have now implemented a framework for computational protein design that allows certain design protocols of Rosetta to make use of multiple design states like structural ensembles. An in silico assessment simulating ligand-binding design showed that this new approach generates more reliably native-like sequences than a single-state approach. As a proof-of-concept, we introduced de novo retro-aldolase activity into a scaffold protein and characterized nine variants experimentally, all of which were catalytically active.

Small Molecule-Induced Domain Swapping as a Mechanism for Controlling Protein Function and Assembly

Domain swapping is the process by which identical proteins exchange reciprocal segments to generate dimers. Here we introduce induced domain swapping (INDOS) as a mechanism for regulating protein function. INDOS employs a modular design consisting of the fusion of two proteins: a recognition protein that binds a triggering molecule, and a target protein that undergoes a domain swap in response to binding of the triggering ligand. The recognition protein (FK506 binding protein) is inserted into functionally-inactivated point mutants of two target proteins (staphylococcal nuclease and ribose binding protein). Binding of FK506 to the FKBP domain causes the target domain to first unfold, then refold via domain swap. The inactivating mutations become ‘swapped out’ in the dimer, increasing nuclease and ribose binding activities by 100-fold and 15-fold, respectively, restoring them to near wild-type values. INDOS is intended to convert an arbitrary protein into a functional switch, and is the first example of rational design in which a small molecule is used to trigger protein domain swapping and subsequent activation of biological function.

Engineered Domain Swapping as an On/Off Switch for Protein Function

Domain swapping occurs when identical proteins exchange segments in reciprocal fashion. Natural swapping mechanisms remain poorly understood, and engineered swapping has the potential for creating self-assembling biomaterials that encode for emergent functions. We demonstrate that induced swapping can be used to regulate the function of a target protein. Swapping is triggered by inserting a "lever" protein (ubiquitin) into one of four loops of the ribose binding protein (RBP) target. The lever splits the target, forcing RBP to refold in trans to generate swapped oligomers. Identical RBP-ubiquitin fusions form homo-swapped complexes with the ubiquitin domain acting as the hinge. Surprisingly, some pairs of non-identical fusions swap more efficiently with each other than they do with themselves. Nuclear magnetic resonance experiments reveal that the hinge of these hetero-swapped complexes maps to a region of RBP distant from both ubiquitins. This design is expected to be applicable to other proteins to convert them into functional switches.

Protein:Ligand binding free energies: A stringent test for computational protein design

A computational protein design method is extended to allow Monte Carlo simulations where two ligands are titrated into a protein binding pocket, yielding binding free energy differences. These provide a stringent test of the physical model, including the energy surface and sidechain rotamer definition. As a test, we consider tyrosyl-tRNA synthetase (TyrRS), which has been extensively redesigned experimentally. We consider its specificity for its substrate l-tyrosine (l-Tyr), compared to the analogs d-Tyr, p-acetyl-, and p-azido-phenylalanine (ac-Phe, az-Phe). We simulate l- and d-Tyr binding to TyrRS and six mutants, and compare the structures and binding free energies to a more rigorous "MD/GBSA" procedure: molecular dynamics with explicit solvent for structures and a Generalized Born + Surface Area model for binding free energies. Next, we consider l-Tyr, ac- and az-Phe binding to six other TyrRS variants. The titration results are sensitive to the precise rotamer definition, which involves a short energy minimization for each sidechain pair to help relax bad contacts induced by the discrete rotamer set. However, when designed mutant structures are rescored with a standard GBSA energy model, results agree well with the more rigorous MD/GBSA. As a third test, we redesign three amino acid positions in the substrate coordination sphere, with either l-Tyr or d-Tyr as the ligand. For two, we obtain good agreement with experiment, recovering the wildtype residue when l-Tyr is the ligand and a d-Tyr specific mutant when d-Tyr is the ligand. For the third, we recover His with either ligand, instead of wildtype Gln.

Protein engineering: a new frontier for biological therapeutics

Protein engineering holds the potential to transform the metabolic drug landscape through the development of smart, stimulusresponsive drug systems. Protein therapeutics are a rapidly expanding segment of Food and Drug Administration approved drugs that will improve clinical outcomes over the long run. Engineering of protein therapeutics is still in its infancy, but recent general advances in protein engineering capabilities are being leveraged to yield improved control over both pharmacokinetics and pharmacodynamics. Stimulus- responsive protein therapeutics are drugs which have been designed to be metabolized under targeted conditions. Protein engineering is being utilized to develop tailored smart therapeutics with biochemical logic. This review focuses on applications of targeted drug neutralization, stimulus-responsive engineered protein prodrugs, and emerging multicomponent smart drug systems (e.g., antibody-drug conjugates, responsive engineered zymogens, prospective biochemical logic smart drug systems, drug buffers, and network medicine applications).

Computational design of enzyme-ligand binding using a combined energy function and deterministic sequence optimization algorithm

Enzyme amino-acid sequences at ligand-binding interfaces are evolutionarily optimized for reactions, and the natural conformation of an enzyme-ligand complex must have a low free energy relative to alternative conformations in native-like or non-native sequences. Based on this assumption, a combined energy function was developed for enzyme design and then evaluated by recapitulating native enzyme sequences at ligand-binding interfaces for 10 enzyme-ligand complexes. In this energy function, the electrostatic interaction between polar or charged atoms at buried interfaces is described by an explicitly orientation-dependent hydrogen-bonding potential and a pairwise-decomposable generalized Born model based on the general side chain in the protein design framework. The energy function is augmented with a pairwise surface-area based hydrophobic contribution for nonpolar atom burial. Using this function, on average, 78% of the amino acids at ligand-binding sites were predicted correctly in the minimum-energy sequences, whereas 84% were predicted correctly in the most-similar sequences, which were selected from the top 20 sequences for each enzyme-ligand complex. Hydrogen bonds at the enzyme-ligand binding interfaces in the 10 complexes were usually recovered with the correct geometries. The binding energies calculated using the combined energy function helped to discriminate the active sequences from a pool of alternative sequences that were generated by repeatedly solving a series of mixed-integer linear programming problems for sequence selection with increasing integer cuts.

Improved energy bound accuracy enhances the efficiency of continuous protein design

Flexibility and dynamics are important for protein function and a protein's ability to accommodate amino acid substitutions. However, when computational protein design algorithms search over protein structures, the allowed flexibility is often reduced to a relatively small set of discrete side-chain and backbone conformations. While simplifications in scoring functions and protein flexibility are currently necessary to computationally search the vast protein sequence and conformational space, a rigid representation of a protein causes the search to become brittle and miss low-energy structures. Continuous rotamers more closely represent the allowed movement of a side chain within its torsional well and have been successfully incorporated into the protein design framework to design biomedically relevant protein systems. The use of continuous rotamers in protein design enables algorithms to search a larger conformational space than previously possible, but adds additional complexity to the design search. To design large, complex systems with continuous rotamers, new algorithms are needed to increase the efficiency of the search. We present two methods, PartCR and HOT, that greatly increase the speed and efficiency of protein design with continuous rotamers. These methods specifically target the large errors in energetic terms that are used to bound pairwise energies during the design search. By tightening the energy bounds, additional pruning of the conformation space can be achieved, and the number of conformations that must be enumerated to find the global minimum energy conformation is greatly reduced.

Escherichia [corrected] coli ribose binding protein based bioreporters revisited

Bioreporter bacteria, i.e., strains engineered to respond to chemical exposure by production of reporter proteins, have attracted wide interest because of their potential to offer cheap and simple alternative analytics for specified compounds or conditions. Bioreporter construction has mostly exploited the natural variation of sensory proteins, but it has been proposed that computational design of new substrate binding properties could lead to completely novel detection specificities at very low affinities. Here we reconstruct a bioreporter system based on the native Escherichia coli ribose binding protein RbsB and one of its computationally designed variants, reported to be capable of binding 2,4,6-trinitrotoluene (TNT). Our results show in vivo reporter induction at 50 nM ribose, and a 125 nM affinity constant for in vitro ribose binding to RbsB. In contrast, the purified published TNT-binding variant did not bind TNT nor did TNT cause induction of the E. coli reporter system.

Systematic optimization model and algorithm for binding sequence selection in computational enzyme design

A systematic optimization model for binding sequence selection in computational enzyme design was developed based on the transition state theory of enzyme catalysis and graph-theoretical modeling. The saddle point on the free energy surface of the reaction system was represented by catalytic geometrical constraints, and the binding energy between the active site and transition state was minimized to reduce the activation energy barrier. The resulting hyperscale combinatorial optimization problem was tackled using a novel heuristic global optimization algorithm, which was inspired and tested by the protein core sequence selection problem. The sequence recapitulation tests on native active sites for two enzyme catalyzed hydrolytic reactions were applied to evaluate the predictive power of the design methodology. The results of the calculation show that most of the native binding sites can be successfully identified if the catalytic geometrical constraints and the structural motifs of the substrate are taken into account. Reliably predicting active site sequences may have significant implications for the creation of novel enzymes that are capable of catalyzing targeted chemical reactions.

Computational enzyme design

Recent developments in computational chemistry and biology have come together in the "inside-out" approach to enzyme engineering. Proteins have been designed to catalyze reactions not previously accelerated in nature. Some of these proteins fold and act as catalysts, but the success rate is still low. The achievements and limitations of the current technology are highlighted and contrasted to other protein engineering techniques. On its own, computational "inside-out" design can lead to the production of catalytically active and selective proteins, but their kinetic performances fall short of natural enzymes. When combined with directed evolution, molecular dynamics simulations, and crowd-sourced structure-prediction approaches, however, computational designs can be significantly improved in terms of binding, turnover, and thermal stability.

Stepwise conversion of a binding protein to a fluorescent switch: application to Thermoanaerobacter tengcongensis ribose binding protein

Alternate frame folding (AFF) is a protein engineering methodology the purpose of which is to convert an ordinary binding protein into a molecular switch. The AFF modification entails duplicating an amino- or carboxy-terminal segment of the protein and appending it to the opposite end of the molecule. This duplication allows the protein to interconvert, in a ligand-dependent fashion, between two mutually exclusive native folds: the wild-type structure and a circularly permuted form. The fold shift can be detected by placement of extrinsic fluorophores at sites sensitive to the engineered conformational change. Here, we apply the AFF mechanism to create several ribose-sensing proteins derived from Thermoanaerobacter tengcongensis ribose binding protein. Our purpose is to systematically explore the parameters of the AFF design. These considerations include the site of circular permutation, the length and location of the duplicated segment, thermodynamic and kinetic optimization of the switching mechanism, and placement of extrinsic fluorophores. Three of the four AFF variants created here undergo the expected conformational shift and exhibit a ribose-dependent fluorescence change. The fourth construct fails to switch folds upon addition of ribose, likely because the circularly permuted form folds much more slowly than the nonpermuted form. This disparity apparently introduces a kinetic barrier that partitions the refolding molecules to the nonpermuted structure. The results of this study serve as a guideline for applying the AFF modification to other proteins of biomedical, diagnostic, and industrial interest.

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.

Multistate approaches in computational protein design

Computational protein design (CPD) is a useful tool for protein engineers. It has been successfully applied towards the creation of proteins with increased thermostability, improved binding affinity, novel enzymatic activity, and altered ligand specificity. Traditionally, CPD calculations search and rank sequences using a single fixed protein backbone template in an approach referred to as single-state design (SSD). While SSD has enjoyed considerable success, certain design objectives require the explicit consideration of multiple conformational and/or chemical states. Cases where a "multistate" approach may be advantageous over the SSD approach include designing conformational changes into proteins, using native ensembles to mimic backbone flexibility, and designing ligand or oligomeric association specificities. These design objectives can be efficiently tackled using multistate design (MSD), an emerging methodology in CPD that considers any number of protein conformational or chemical states as inputs instead of a single protein backbone template, as in SSD. In this review article, recent examples of the successful design of a desired property into proteins using MSD are described. These studies employing MSD are divided into two categories--those that utilized multiple conformational states, and those that utilized multiple chemical states. In addition, the scoring of competing states during negative design is discussed as a current challenge for MSD.

QM/MM refinement and analysis of protein bound retinoic acid

Retinoic acid (RA) is a vitamin A derivative, which modifies the appearance of fine wrinkles and roughness of facial skin and treats acne and activates gene transcription by binding to heterodimers of the retinoic acid receptor (RAR) and the retinoic X receptor (RXR). There are series of protein bound RA complexes available in the protein databank (PDB), which provides a broad range of information about the different bioactive conformations of RA. To gain further insights into the observed bioactive RA conformations we applied quantum mechanic (QM)/molecular mechanic (MM) approaches to re-refine the available RA protein-ligand complexes. MP2 complete basis set (CBS) extrapolations single energy calculations are also carried out for both the experimental conformations and QM optimized geometries of RA in the gas as well as solution phase. The results demonstrate that the re-refined structures show better geometries for RA than seen in the originally deposited PDB structures through the use of QMs for the ligand in the X-ray refinement procedure. QM/MM re-refined conformations also reduced the computed strain energies found in the deposited crystal conformations for RA. Finally, the dependence of ligand strain on resolution is analyzed. It is shown that ligand strain is not converged in our calculations and is likely an artifact of the typical resolutions employed to study protein-ligand complexes.

Generation of recombinant antibodies and means for increasing their affinity

Highly specific interaction with foreign molecules is a unique feature of antibodies. Since 1975, when Keller and Milstein proposed the method of hybridoma technology and prepared mouse monoclonal antibodies, many antibodies specific to various antigens have been obtained. Recent development of methods for preparation of recombinant DNA libraries and in silico bioinformatics approaches for protein structure analysis makes possible antibody preparation using gene engineering approaches. The development of gene engineering methods allowed creating recombinant antibodies and improving characteristics of existing antibodies; this significantly extends the applicability of antibodies. By modifying biochemical and immunochemical properties of antibodies by changing their amino acid sequences it is possible to create antibodies with properties optimal for certain tasks. For example, application of recombinant technologies resulted in antibody preparation of high affinity significantly exceeding the initial affinity of natural antibodies. In this review we summarize information about the structure, modes of preparation, and application of recombinant antibodies and their fragments and also consider the main approaches used to increase antibody affinity.

Computational protein design with a generalized Born solvent model: application to Asparaginyl-tRNA synthetase

Computational Protein Design (CPD) is a promising method for high throughput protein and ligand mutagenesis. Recently, we developed a CPD method that used a polar-hydrogen energy function for protein interactions and a Coulomb/Accessible Surface Area (CASA) model for solvent effects. We applied this method to engineer aspartyl-adenylate (AspAMP) specificity into Asparaginyl-tRNA synthetase (AsnRS), whose substrate is asparaginyl-adenylate (AsnAMP). Here, we implement a more accurate function, with an all-atom energy for protein interactions and a residue-pairwise generalized Born model for solvent effects. As a first test, we compute aminoacid affinities for several point mutants of Aspartyl-tRNA synthetase (AspRS) and Tyrosyl-tRNA synthetase and stability changes for three helical peptides and compare with experiment. As a second test, we readdress the problem of AsnRS aminoacid engineering. We compare three design criteria, which optimize the folding free-energy, the absolute AspAMP affinity, and the relative (AspAMP-AsnAMP) affinity. The sequences and conformations are improved with respect to our previous, polar-hydrogen/CASA study: For several designed complexes, the AspAMP carboxylate forms three interactions with a conserved arginine and a designed lysine, as in the active site of the AspRS:AspAMP complex. The conformations and interactions are well maintained in molecular dynamics simulations and the sequences have an inverted specificity, favoring AspAMP over AsnAMP. The method is not fully successful, since experimental measurements with the seven most promising sequences show that they do not catalyze at a detectable level the adenylation of Asp (or Asn) with ATP. This may be due to weak AspAMP binding and/or disruption of transition-state stabilization.

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

Computational design of protein-ligand interfaces finds optimal amino acid sequences within a small-molecule binding site of a protein for tight binding of a specific small molecule. It requires a search algorithm that can rapidly sample the vast sequence and conformational space, and a scoring function that can identify low energy designs. This review focuses on recent advances in computational design methods and their application to protein-small molecule binding sites. Strategies for increasing affinity, altering specificity, creating broad-spectrum binding, and building novel enzymes from scratch are described. Future prospects for applications in drug development are discussed, including limitations that will need to be overcome to achieve computational design of protein therapeutics with novel modes of action.

The measured and calculated affinity of methyl- and methoxy-substituted benzoquinones for the Q(A) site of bacterial reaction centers

Quinones play important roles in mitochondrial and photosynthetic energy conversion acting as intramembrane, mobile electron, and proton carriers between catalytic sites in various electron transfer proteins. They display different affinity, selectivity, functionality, and exchange dynamics in different binding sites. The computational analysis of quinone binding sheds light on the requirements for quinone affinity and specificity. The affinities of 10 oxidized, neutral benzoquinones were measured for the high affinity Q(A) site in the detergent-solubilized Rhodobacter sphaeroides bacterial photosynthetic reaction center. Multiconformation Continuum Electrostatics was then used to calculate their relative binding free energies by grand canonical Monte Carlo sampling with a rigid protein backbone, flexible ligand, and side chain positions and protonation states. Van der Waals and torsion energies, Poisson-Boltzmann continuum electrostatics, and accessible surface area-dependent ligand-solvent interactions are considered. An initial, single cycle of GROMACS backbone optimization improves the match with experiment as do coupled-ligand and side-chain motions. The calculations match experiment with an root mean square deviation (RMSD) of 2.29 and a slope of 1.28. The affinities are dominated by favorable protein-ligand van der Waals rather than electrostatic interactions. Each quinone appears in a closely clustered set of positions. Methyl and methoxy groups move into the same positions as found for the native quinone. Difficulties putting methyls into methoxy sites are observed. Calculations using a solvent-accessible surface area-dependent implicit van der Waals interaction smoothed out small clashes, providing a better match to experiment with a RMSD of 0.77 and a slope of 0.97.

Experimental library screening demonstrates the successful application of computational protein design to large structural ensembles

The stability, activity, and solubility of a protein sequence are determined by a delicate balance of molecular interactions in a variety of conformational states. Even so, most computational protein design methods model sequences in the context of a single native conformation. Simulations that model the native state as an ensemble have been mostly neglected due to the lack of sufficiently powerful optimization algorithms for multistate design. Here, we have applied our multistate design algorithm to study the potential utility of various forms of input structural data for design. To facilitate a more thorough analysis, we developed new methods for the design and high-throughput stability determination of combinatorial mutation libraries based on protein design calculations. The application of these methods to the core design of a small model system produced many variants with improved thermodynamic stability and showed that multistate design methods can be readily applied to large structural ensembles. We found that exhaustive screening of our designed libraries helped to clarify several sources of simulation error that would have otherwise been difficult to ascertain. Interestingly, the lack of correlation between our simulated and experimentally measured stability values shows clearly that a design procedure need not reproduce experimental data exactly to achieve success. This surprising result suggests potentially fruitful directions for the improvement of computational protein design technology.

Conformational diversity and computational enzyme design

The application of computational protein design methods to the design of enzyme active sites offers potential routes to new catalysts and new reaction specificities. Computational design methods have typically treated the protein backbone as a rigid structure for the sake of computational tractability. However, this fixed-backbone approximation introduces its own special challenges for enzyme design and it contrasts with an emerging picture of natural enzymes as dynamic ensembles with multiple conformations and motions throughout a reaction cycle. This review considers the impact of conformational variation and dynamics on computational enzyme design and it highlights new approaches to addressing protein conformational diversity in enzyme design including recent advances in multi-state design, backbone flexibility, and computational library design.

Specificity in computational protein design

A long-standing goal of computational protein design is to create proteins similar to those found in Nature. One motivation is to harness the exquisite functional capabilities of proteins for our own purposes. The extent of similarity between designed and natural proteins also reports on how faithfully our models represent the selective pressures that determine protein sequences. As the field of protein design shifts emphasis from reproducing native-like protein structure to function, it has become important that these models treat the notion of specificity in molecular interactions. Although specificity may, in some cases, be achieved by optimization of a desired protein in isolation, methods have been developed to address directly the desire for proteins that exhibit specific functions and interactions.

Computational design of protein-ligand binding: modifying the specificity of asparaginyl-tRNA synthetase

A method for computational design of protein-ligand interactions is implemented and tested on the asparaginyl- and aspartyl-tRNA synthetase enzymes (AsnRS, AspRS). The substrate specificity of these enzymes is crucial for the accurate translation of the genetic code. The method relies on a molecular mechanics energy function and a simple, continuum electrostatic, implicit solvent model. As test calculations, we first compute AspRS-substrate binding free energy changes due to nine point mutations, for which experimental data are available; we also perform large-scale redesign of the entire active site of each enzyme (40 amino acids) and compare to experimental sequences. We then apply the method to engineer an increased binding of aspartyl-adenylate (AspAMP) into AsnRS. Mutants are obtained using several directed evolution protocols, where four or five amino acid positions in the active site are randomized. Promising mutants are subjected to molecular dynamics simulations; Poisson-Boltzmann calculations provide an estimate of the corresponding, AspAMP, binding free energy changes, relative to the native AsnRS. Several of the mutants are predicted to have an inverted binding specificity, preferring to bind AspAMP rather than the natural substrate, AsnAMP. The computed binding affinities are significantly weaker than the native, AsnRS:AsnAMP affinity, and in most cases, the active site structure is significantly changed, compared to the native complex. This almost certainly precludes catalytic activity. One of the designed sequences has a higher affinity and more native-like structure and may represent a valid candidate for Asp activity.

An efficient algorithm for multistate protein design based on FASTER

Most of the methods that have been developed for computational protein design involve the selection of side-chain conformations in the context of a single, fixed main-chain structure. In contrast, multistate design (MSD) methods allow sequence selection to be driven by the energetic contributions of multiple structural or chemical states simultaneously. This methodology is expected to be useful when the design target is an ensemble of related states rather than a single structure, or when a protein sequence must assume several distinct conformations to function. MSD can also be used with explicit negative design to suggest sequences with altered structural, binding, or catalytic specificity. We report implementation details of an efficient multistate design optimization algorithm based on FASTER (MSD-FASTER). We subjected the algorithm to a battery of computational tests and found it to be generally applicable to various multistate design problems; designs with a large number of states and many designed positions are completely feasible. A direct comparison of MSD-FASTER and multistate design Monte Carlo indicated that MSD-FASTER discovers low-energy sequences much more consistently. MSD-FASTER likely performs better because amino acid substitutions are chosen on an energetic basis rather than randomly, and because multiple substitutions are applied together. Through its greater efficiency, MSD-FASTER should allow protein designers to test experimentally better-scoring sequences, and thus accelerate progress in the development of improved scoring functions and models for computational protein design.

Evolution: a guide to perturb protein function and networks

Protein interactions give rise to networks that control cell fate in health and disease; selective means to probe these interactions are therefore of wide interest. We discuss here Evolutionary Tracing (ET), a comparative method to identify protein functional sites and to guide experiments that selectively block, recode, or mimic their amino acid determinants. These studies suggest, in principle, a scalable approach to perturb individual links in protein networks.

Computational design of ligand binding is not a solved problem

Computational design has been very successful in recent years: multiple novel ligand binding proteins as well as enzymes have been reported. We wanted to know in molecular detail how precise the predictions of the interactions of protein and ligands are. Therefore, we performed a structural analysis of a number of published receptors designed onto the periplasmic binding protein scaffold that were reported to bind to the new ligands with nano- to micromolar affinities. It turned out that most of these designed proteins are not suitable for structural studies due to instability and aggregation. However, we were able to solve the crystal structure of an arabinose binding protein designed to bind serotonin to 2.2 A resolution. While crystallized in the presence of an excess of serotonin, the protein is in an open conformation with no serotonin bound, although the side-chain conformations in the empty binding pocket are very similar to the conformations predicted. During subsequent characterization using isothermal titration calorimetry, CD, and NMR spectroscopy, no indication of binding could be detected for any of the tested designed receptors, whereas wild-type proteins bound their ligands as expected. We conclude that although the computational prediction of side-chain conformations appears to be working, it does not necessarily confer binding as expected. Hence, the computational design of ligand binding is not a solved problem and needs to be revisited.

SEST: Simulated Electronic Structure Theory

A novel approach to empirically modeling the electronic structure of molecules is introduced. The theory is based on relationships between molecular orbital energy components and the average distance between electrons and electrons and nuclei. The electron-electron and electron-nucleus distances are subsequently related to interatomic distances which provides a means for modeling the electronic structure of molecules. The general energy expression for a simulated electronic structure theory is defined, along with the functional form of the interatomic distance dependent energy functions. The theory is used to model the hydrogen molecule, the first-row hydrides, and ethane. The models, which have the correct RHF/6-31G(d) optimized geometries, also fit the RHF/6-31G(d) energy at equilibrium and the UHF/6-31G(d) energy at the bond dissociation limit as well as some vibrational frequencies.