Protein redesign

The Curious Case of A31P, a Topology-Switching Mutant of the Repressor of Primer Protein: A Molecular Dynamics Study of Its Folding and Misfolding

The effect of mutations on protein structures is usually rather localized and minor. Finding a mutation that can single-handedly change the fold and/or topology of a protein structure is a rare exception. The A31P mutant of the homodimeric Repressor of primer (Rop) protein is one such exception: This single mutation ─and as demonstrated by two independent crystal structure determinations─ can convert the canonical (left-handed/all-antiparallel) 4-α-helical bundle of Rop to a new form (right-handed/mixed parallel and antiparallel bundle) displaying a previously unobserved "bisecting U" topology. The main problem with understanding the dramatic effect of this mutation on the folding of Rop is to understand its very existence: Most computational methods appear to agree that the mutation should have had no appreciable effect, with the majority of energy minimization methods and protein structure prediction protocols indicating that this mutation is fully consistent with the native Rop structure, requiring only a local and minor change at the mutation site. Here we use two long (10 μs each) molecular dynamics simulations to compare the stability and dynamics of the native Rop versus a hypothetical structure that is identical with the native Rop but is carrying this single Alanine31 to Proline mutation. Comparative analysis of the two trajectories convincingly shows that, in contrast to the indications from energy minimization ─but in agreement with the experimental data─, this hypothetical native-like A31P structure is unstable, with its turn regions almost completely unfolding, even under the relatively mild 320 K NpT simulations that we have used for this study. We discuss the implication of these findings for the folding of the A31P mutant, especially with respect to the proposed model of a double-funneled energy landscape.

PLUG (Pruning of Local Unrealistic Geometries) removes restrictions on biophysical modeling for protein design

Protein design algorithms must search an enormous conformational space to identify favorable conformations. As a result, those that perform this search with guarantees of accuracy generally start with a conformational pruning step, such as dead-end elimination (DEE). However, the mathematical assumptions of DEE-based pruning algorithms have up to now severely restricted the biophysical model that can feasibly be used in protein design. To lift these restrictions, I propose to prune local unrealistic geometries (PLUG) using a linear programming-based method. PLUG's biophysical model consists only of well-known lower bounds on interatomic distances. PLUG is intended as preprocessing for energy-based protein design calculations, whose biophysical model need not support DEE pruning. Based on 96 test cases, PLUG is at least as effective at pruning as DEE for larger protein designs-the type that most require pruning. When combined with the LUTE protein design algorithm, PLUG greatly facilitates designs that account for continuous entropy, large multistate designs with continuous flexibility, and designs with extensive continuous backbone flexibility and advanced nonpairwise energy functions. Many of these designs are tractable only with PLUG, either for empirical reasons (LUTE's machine learning step achieves an accurate fit only after PLUG pruning), or for theoretical reasons (many energy functions are fundamentally incompatible with DEE).

CATS (Coordinates of Atoms by Taylor Series): protein design with backbone flexibility in all locally feasible directions

Motivation

When proteins mutate or bind to ligands, their backbones often move significantly, especially in loop regions. Computational protein design algorithms must model these motions in order to accurately optimize protein stability and binding affinity. However, methods for backbone conformational search in design have been much more limited than for sidechain conformational search. This is especially true for combinatorial protein design algorithms, which aim to search a large sequence space efficiently and thus cannot rely on temporal simulation of each candidate sequence.

Results

We alleviate this difficulty with a new parameterization of backbone conformational space, which represents all degrees of freedom of a specified segment of protein chain that maintain valid bonding geometry (by maintaining the original bond lengths and angles and ω dihedrals). In order to search this space, we present an efficient algorithm, CATS, for computing atomic coordinates as a function of our new continuous backbone internal coordinates. CATS generalizes the iMinDEE and EPIC protein design algorithms, which model continuous flexibility in sidechain dihedrals, to model continuous, appropriately localized flexibility in the backbone dihedrals ϕ and ψ as well. We show using 81 test cases based on 29 different protein structures that CATS finds sequences and conformations that are significantly lower in energy than methods with less or no backbone flexibility do. In particular, we show that CATS can model the viability of an antibody mutation known experimentally to increase affinity, but that appears sterically infeasible when modeled with less or no backbone flexibility.

Availability and implementation

Our code is available as free software at https://github.com/donaldlab/OSPREY_refactor.

Supplementary information

Supplementary data are available at Bioinformatics online.

From DARPins to LoopDARPins: novel LoopDARPin design allows the selection of low picomolar binders in a single round of ribosome display

Antibodies are the most versatile binding proteins in nature with six loops creating a flexible continuous interaction surface. However, in some molecular formats, antibodies are aggregation prone. Designed ankyrin repeat proteins (DARPins) were successfully created as alternative design solutions. Nevertheless, their concave shape, rigidity and incompletely randomized binding surface may limit the epitopes that can be targeted by this extremely stable scaffold. Combining conformational diversity and a continuous convex paratope found in many antibodies with the beneficial biophysical properties of DARPins, we created LoopDARPins, a next generation of DARPins with extended epitope binding properties. We employed X-ray structure determination of a LoopDARPin for design validation. Biophysical characterizations show that the introduction of an elongated loop through consensus design does not decrease the stability of the scaffold,consistent with molecular dynamics simulations. Ribosome-display selections against extracellular signal-regulated kinase 2 (ERK2) and four members of the BCL-2 family (BCL-2, BCL-XL, BCL-W and MCL-1) of anti-apoptotic regulators yielded LoopDARPins with affinities in the mid-picomolar to low nanomol arrange against all targets. The BCL-2 family binders block the interaction with their natural interaction partner and will be valuable reagents to test the apoptotic response in functional assays. With the LoopDARPin scaffold, binders for BCL-2 with an affinity of 30 pM were isolated with only a single round of ribosome display,an enrichment that has not been described for any scaffold. Identical stringent one-round selections with conventional DARPins without loop yielded no binders. The LoopDARPin scaffold may become a highly valuable tool for biotechnological high-throughput applications.

Site-Specific Protein Bioconjugation via a Pyridoxal 5'-Phosphate-Mediated N-Terminal Transamination Reaction

The covalent attachment of chemical groups to proteins is a critically important tool for the study of protein function and the creation of protein-based materials. Methods of site-specific protein modification are necessary for the generation of well defined bioconjugates possessing a new functional group in a single position in the amino acid sequence. This article describes a pyridoxal 5'-phosphate (PLP)-mediated transamination reaction that is specific for the N-terminus of a protein. The reaction oxidizes the N-terminal amine to a ketone or an aldehyde, which can form a stable oxime linkage with an alkoxyamine reagent of choice. Screening studies have identified the most reactive N-terminal residues, facilitating the use of site-directed mutagenesis to achieve high levels of conversion. Additionally, this reaction has been shown to be effective for a number of targets that are not easily accessed through heterologous expression, such as monoclonal antibodies. Curr. Protoc. Chem. Biol. 2:125-134 © 2010 by John Wiley & Sons, Inc.

Dead-end elimination with perturbations (DEEPer): a provable protein design algorithm with continuous sidechain and backbone flexibility

Computational protein and drug design generally require accurate modeling of protein conformations. This modeling typically starts with an experimentally determined protein structure and considers possible conformational changes due to mutations or new ligands. The DEE/A* algorithm provably finds the global minimum-energy conformation (GMEC) of a protein assuming that the backbone does not move and the sidechains take on conformations from a set of discrete, experimentally observed conformations called rotamers. DEE/A* can efficiently find the overall GMEC for exponentially many mutant sequences. Previous improvements to DEE/A* include modeling ensembles of sidechain conformations and either continuous sidechain or backbone flexibility. We present a new algorithm, DEEPer (Dead-End Elimination with Perturbations), that combines these advantages and can also handle much more extensive backbone flexibility and backbone ensembles. DEEPer provably finds the GMEC or, if desired by the user, all conformations and sequences within a specified energy window of the GMEC. It includes the new abilities to handle arbitrarily large backbone perturbations and to generate ensembles of backbone conformations. It also incorporates the shear, an experimentally observed local backbone motion never before used in design. Additionally, we derive a new method to accelerate DEE/A*-based calculations, indirect pruning, that is particularly useful for DEEPer. In 67 benchmark tests on 64 proteins, DEEPer consistently identified lower-energy conformations than previous methods did, indicating more accurate modeling. Additional tests demonstrated its ability to incorporate larger, experimentally observed backbone conformational changes and to model realistic conformational ensembles. These capabilities provide significant advantages for modeling protein mutations and protein-ligand interactions.

Chimeric glutathione S-transferases containing inserts of kininogen peptides: potential novel protein therapeutics

The study of synthetic peptides corresponding to discrete regions of proteins has facilitated the understanding of protein structure-activity relationships. Short peptides can also be used as powerful therapeutic agents. However, in many instances, small peptides are prone to rapid degradation or aggregation and may lack the conformation required to mimic the functional motifs of the protein. For peptides to function as pharmacologically active agents, efficient production or expression, high solubility, and retention of biological activity through purification and storage steps are required. We report here the design, expression, and functional analysis of eight engineered GST proteins (denoted GSHKTs) in which peptides ranging in size from 8 to 16 amino acids and derived from human high molecular weight kininogen (HK) domain 5 were inserted into GST (between Gly-49 and Leu-50). Peptides derived from HK are known to inhibit cell proliferation, angiogenesis, and tumor metastasis, and the biological activity of the HK peptides was dramatically (>50-fold) enhanced following insertion into GST. GSHKTs are soluble and easily purified from Escherichia coli by affinity chromatography. Functionally, these hybrid proteins cause inhibition of endothelial cell proliferation. Crystallographic analysis of GSHKT10 and GSHKT13 (harboring 10- and 13-residue HK peptides, respectively) showed that the overall GST structure was not perturbed. These results suggest that the therapeutic efficacy of short peptides can be enhanced by insertion into larger proteins that are easily expressed and purified and that GST may potentially be used as such a carrier.

A rewired green fluorescent protein: folding and function in a nonsequential, noncircular GFP permutant

The sequential order of secondary structural elements in proteins affects the folding and activity to an unknown extent. To test the dependence on sequential connectivity, we reconnected secondary structural elements by their solvent-exposed ends, permuting their sequential order, called "rewiring". This new protein design strategy changes the topology of the backbone without changing the core side chain packing arrangement. While circular and noncircular permutations have been observed in protein structures that are not related by sequence homology, to date no one has attempted to rationally design and construct a protein with a sequence that is noncircularly permuted while conserving three-dimensional structure. Herein, we show that green fluorescent protein can be rewired, still functionally fold, and exhibit wild-type fluorescence excitation and emission spectra.

Nanoscale protein assemblies from a circular permutant of the tobacco mosaic virus

The protein coat of the tobacco mosaic virus (TMV) has been explored extensively for the construction of nanoscale architectures. In previous work, we have reported efficient TMV-based light harvesting systems bearing chromophores in a hollow channel of the assembled protein. We have also reported an N-terminal transamination/oximation method that could be used to attach electrodes and catalytic groups to the exterior surface of the rods. To complement these techniques, we report herein a new circular permutant of the TMV capsid protein that repositions the N- and C-termini to the center of the assemblies. This protein can be produced in very high yield through E. coli expression and self-assembles into light harvesting rods that are much like those assembled from the wild-type protein. However, the disks formed from the permutant structure are stable over a significantly wider pH range, greatly improving the practicality of this assembled form for materials applications. The new position of the N-terminus allows functional groups to be installed in the inner pore of the disks, affording geometries reminiscent of natural photosynthetic systems. The permutant also shows the ability to coassemble with regular monomers, allowing the future generation of multicomponent rod structures that are modified on the exterior and interior surfaces, as well as in the internal RNA channel.

Laboratory-evolved vanillyl-alcohol oxidase produces natural vanillin

The flavoenzyme vanillyl-alcohol oxidase was subjected to random mutagenesis to generate mutants with enhanced reactivity to creosol (2-methoxy-4-methylphenol). The vanillyl-alcohol oxidase-mediated conversion of creosol proceeds via a two-step process in which the initially formed vanillyl alcohol (4-hydroxy-3-methoxybenzyl alcohol) is oxidized to the widely used flavor compound vanillin (4-hydroxy-3-methoxybenzaldehyde). The first step of this reaction is extremely slow due to the formation of a covalent FAD N-5-creosol adduct. After a single round of error-prone PCR, seven mutants were generated with increased reactivity to creosol. The single-point mutants I238T, F454Y, E502G, and T505S showed an up to 40-fold increase in catalytic efficiency (kcat/Km) with creosol compared with the wild-type enzyme. This enhanced reactivity was due to a lower stability of the covalent flavin-substrate adduct, thereby promoting vanillin formation. The catalytic efficiencies of the mutants were also enhanced for other ortho-substituted 4-methylphenols, but not for p-cresol (4-methylphenol). The replaced amino acid residues are not located within a distance of direct interaction with the substrate, and the determined three-dimensional structures of the mutant enzymes are highly similar to that of the wild-type enzyme. These results clearly show the importance of remote residues, not readily predicted by rational design, for the substrate specificity of enzymes.

Combinatorial approaches to protein stability and structure

Why do proteins adopt the conformations that they do, and what determines their stabilities? While we have come to some understanding of the forces that underlie protein architecture, a precise, predictive, physicochemical explanation is still elusive. Two obstacles to addressing these questions are the unfathomable vastness of protein sequence space, and the difficulty in making direct physical measurements on large numbers of protein variants. Here, we review combinatorial methods that have been applied to problems in protein biophysics over the last 15 years. The effects of hydrophobic core composition, the most important determinant of structure and stability, are still poorly understood. Particular attention is given to core composition as addressed by library methods. Increasingly useful screens and selections, in combination with modern high-throughput approaches borrowed from genomics and proteomics efforts, are making the empirical, statistical correlation between sequence and structure a tractable problem for the coming years.

Mapping of unit boundaries of a protein: exhaustive search for permissive sites for duplication by complementation analysis of random fragment libraries of tryptophan synthase alpha subunit

To identify peptide units that make up a single-domain protein, we searched for possible combinations of N and C-fragments that exhibit functional complementation, and attempted an exhaustive evaluation of peptide unit boundaries. The tryptophan synthase alpha subunit was used as a model enzyme, which has a single TIM (beta8/alpha8) barrel domain. Libraries comprising randomly digested N and C-fragments were constructed, and clones expressing enzymatic activity were selected by the ability to confer growth of the Escherichia coli trpA mutant on a medium lacking tryptophan. More than 50 clones were obtained, and two cleavable positions were found on the loops after extra-helix 2' and helix 5. Half of the clones harbored N and C-fragments having an overlap between two fragments. The remaining clones harbored one fragment made by the fusion of N and C-fragments with insertional sequence duplication. Mapping the frequency of occurrence of fragment overlap and insertional duplication showed significant peaks at two loops, which coincide with the cleavable sites. These results suggest that the boundaries of unit fragments are located at these positions, and that bisection, fragment overlap and insertion are all possible at the unit boundaries.

In vitro disassembly of a parvovirus capsid and effect on capsid stability of heterologous peptide insertions in surface loops

We have analyzed the in vitro disassembly of the capsid of the minute virus of mice, and the stability of capsid chimeras carrying heterologous epitope insertions. Upon heating in a physiological buffer, empty capsids formed by 60 copies of protein VP2 underwent first a reversible conformational change with a small enthalpy change detected by fluorescence. This change was associated with, but not limited to, externalization of the VP2 N terminus. Irreversible capsid dissociation as detected by changes in fluorescence, hemagglutination activity, and electrophoretic mobility occurred at much higher temperatures. Differential scanning calorimetry in the same conditions indicated that the dissociation/denaturation transition involved a high enthalpy change and proceeded through one or more intermediates. In contrast, in the presence of 1.5 M guanidinium chloride, heat-induced disassembly fitted a two-state irreversible process. Both thermally and chemically induced dissociation/denaturation yielded a form that had lost a part of the tertiary structure, but still retained the native secondary structure. Data from chemical dissociation indicates this form may correspond to a molten globule-like monomeric state of the capsid protein. All five antigenic peptide insertions attempted in exposed loops, despite being perhaps among the least disruptive, led to defects in folding/assembly of the capsid and, in most cases, to reduced capsid stability against thermal dissociation. The results with one of the simplest viral capsids reveal a complex pathway for disassembly, and a reduction in capsid assembly and stability upon insertion of peptides, even within the most exposed capsid loops.

Design of protein struts for self-assembling nanoconstructs

Bacteriophage T4 tail fibers have a quaternary structure of bent rigid rods, 3 x 160 nm in size. The four proteins which make up these organelles are able to self-assemble in an essentially irreversible manner. To use the self-assembly domains of these proteins as elements in construction of mesoscale structures, we must be able to rearrange these domains without affecting the self-assembly properties and add internal binding sites for other functional elements. Here we present results on several alterations of the P37 component of the T4 tail fiber that change its length and add novel protein sequences into the protein. One of these sequences is an antibody binding site that is used to inactivate phage carrying the modified gene.

Denatured state thermodynamics: residual structure, chain stiffness and scaling factors

A set of nine variants of yeast iso-1-cytochrome c with zero or one surface histidine have been engineered such that the N-terminal amino group is acetylated in vivo. N-terminal acetylation has been confirmed by mass spectral analysis of intact and proteolytically digested protein. The histidine-heme loop-forming equilibrium, under denaturing conditions (3 M guanidine hydrochloride), has been measured by pH titration providing an observed pK(a), pK(a)(obs), for each variant. N-terminal acetylation prevents the N-terminal amino group-heme binding equilibrium from interfering with measurements of histidine-heme affinity. Significant deviation is observed from the linear dependence of pK(a)(obs) on the log of the number of monomers in the loop formed, expected for a random coil denatured state. The maximum histidine-heme affinity occurs for a loop size of 37 monomers. For loop sizes of 37-83 monomers, histidine-heme pK(a)(obs) values are consistent with a scaling factor of -4.2+/-0.3. This value is much larger than the scaling factor of -1.5 for a freely jointed random coil, which is commonly used to represent the conformational properties of protein denatured states. For loop sizes of nine to 22 monomers, chain stiffness is likely responsible for the decreases in histidine-heme affinity relative to a loop size of 37. The results are discussed in terms of residual structure and sequence composition effects on the conformational properties of the denatured states of proteins.

Designing metal-peptide models for protein structure and function

Recent progress in the rational design of metal sites within peptide model systems shows increasing control in the placement of metals within helical bundles and inclusion of sophisticated elements such as second-sphere ligand interactions. A crystallographically characterized two-metal peptide model for diiron proteins represents a major achievement in de novo design methodologies. Increasingly complex and robust models for electron transfer through and between helices, and electrode-supported electron-transfer peptides, have been constructed. Design elements for peptide-supported ferredoxins and mononuclear Fe(II) and Zn(II) sites have been refined.

Altering dimerization specificity by changes in surface electrostatics

Arc repressor forms a homodimer in which the subunits intertwine to create a single globular domain. To obtain Arc sequences that fold preferentially as heterodimers, variants with surface patches of excess positive or negative charge were designed. Several but not all oppositely charged sequence pairs showed preferential heterodimer formation. In the most successful design pair, alpha helix B of one subunit contained glutamic acids at positions 43, 46, 47, 48, and 50, whereas the other subunit contained lysines or arginines at these positions. A continuum electrostatic model captures many features of the experimental results and suggests that the most successful designs include elements of both positive and negative design.

Modeling, mutagenesis, and structural studies on the fully conserved phosphate-binding loop (loop 8) of triosephosphate isomerase: toward a new substrate specificity

Loop 8 (residues 232-242) in triosephosphate isomerase (TIM) is a highly conserved loop that forms a tight binding pocket for the phosphate moiety of the substrate. Its sequence includes the fully conserved, solvent-exposed Leu238. The tight phosphate-binding pocket explains the high substrate specificity of TIM being limited to the in vivo substrates dihydroxyacetone-phosphate and D-glyceraldehyde-3-phosphate. Here we use the monomeric variant of trypanosomal TIM for exploring the structural consequences of shortening this loop. The mutagenesis, guided by extensive modeling calculations and followed up by crystallographic characterization, is aimed at widening the phosphate-binding pocket and, consequently, changing the substrate specificity. Two new variants were characterized. The crystal structures of these variants indicate that in monomeric forms of TIM, the Leu238 side-chain is nicely buried in a hydrophobic cluster. Monomeric forms of wild-type dimeric TIM are known to exist transiently as folding intermediates; our structural analysis suggests that in this monomeric form, Leu238 of loop 8 also adopts this completely buried conformation, which explains its full conservation across the evolution. The much wider phosphate-binding pocket of the new variant allows for the development of a new TIM variant with a different substrate specificity.

beta-hairpin stability and folding: molecular dynamics studies of the first beta-hairpin of tendamistat

The stability and (un)folding of the 19-residue peptide, SCVTLYQSWRYSQADNGCA, corresponding to the first beta-hairpin (residues 10 to 28) of the alpha-amylase inhibitor tendamistat (PDB entry 3AIT) has been studied by molecular dynamics simulations in explicit water under periodic boundary conditions at several temperatures (300 K, 360 K and 400 K), starting from various conformations for simulation lengths, ranging from 10 to 30 ns. Comparison of trajectories of the reduced and oxidized native peptides reveals the importance of the disulphide bridge closing the beta-hairpin in maintaining a proper turn conformation, thereby insuring a proper side-chain arrangement of the conserved turn residues. This allows rationalization of the conservation of those cysteine residues among the family of alpha-amylase inhibitors. High temperature simulations starting from widely different initial configurations (native beta-hairpin, alpha and left-handed helical and extended conformations) begin sampling similar regions of the conformational space within tens of nanoseconds, and both native and non-native beta-hairpin conformations are recovered. Transitions between conformational clusters are accompanied by an increase in energy fluctuations, which is consistent with the increase in heat capacity measured experimentally upon protein folding. The folding events observed in the various simulations support a model for beta-hairpin formation in which the turn is formed first, followed by hydrogen bond formation closing the hairpin, and subsequent stabilization by side-chain hydrophobic interactions.