BAlaS: fast, interactive and accessible computational alanine-scanning using BudeAlaScan.

MOTIVATION: In experimental protein engineering, alanine-scanning mutagenesis involves the replacement of selected residues with alanine to determine the energetic contribution of each side chain to forming an interaction. For example, it is often used to study protein-protein interactions. However, such experiments can be time-consuming and costly, which has led to the development of programmes for performing computational alanine-scanning mutagenesis (CASM) to guide experiments. While programmes are available for this, there is a need for a real-time web application that is accessible to non-expert users. RESULTS: Here, we present BAlaS, an interactive web application for performing CASM via BudeAlaScan and visualizing its results. BAlaS is interactive and intuitive to use. Results are displayed directly in the browser for the structure being interrogated enabling their rapid inspection. BAlaS has broad applications in areas, such as drug discovery and protein-interface design. AVAILABILITY AND IMPLEMENTATION: BAlaS works on all modern browsers and is available through the following website: https://balas.app. The project is open source, distributed using an MIT license and is available on GitHub (https://github.com/wells-wood-research/balas).

Applying graph theory to protein structures: an Atlas of coiled coils.

Motivation: To understand protein structure, folding and function fully and to design proteins de novo reliably, we must learn from natural protein structures that have been characterized experimentally. The number of protein structures available is large and growing exponentially, which makes this task challenging. Indeed, computational resources are becoming increasingly important for classifying and analyzing this resource. Here, we use tools from graph theory to define an Atlas classification scheme for automatically categorizing certain protein substructures. Results: Focusing on the α-helical coiled coils, which are ubiquitous protein-structure and protein-protein interaction motifs, we present a suite of computational resources designed for analyzing these assemblies. iSOCKET enables interactive analysis of side-chain packing within proteins to identify coiled coils automatically and with considerable user control. Applying a graph theory-based Atlas classification scheme to structures identified by iSOCKET gives the Atlas of Coiled Coils, a fully automated, updated overview of extant coiled coils. The utility of this approach is illustrated with the first formal classification of an emerging subclass of coiled coils called α-helical barrels. Furthermore, in the Atlas, the known coiled-coil universe is presented alongside a partial enumeration of the 'dark matter' of coiled-coil structures; i.e. those coiled-coil architectures that are theoretically possible but have not been observed to date, and thus present defined targets for protein design. Availability and implementation: iSOCKET is available as part of the open-source GitHub repository associated with this work (https://github.com/woolfson-group/isocket). This repository also contains all the data generated when classifying the protein graphs. The Atlas of Coiled Coils is available at: http://coiledcoils.chm.bris.ac.uk/atlas/app.

Engineering protein stability with atomic precision in a monomeric miniprotein.

Miniproteins simplify the protein-folding problem, allowing the dissection of forces that stabilize protein structures. Here we describe PPα-Tyr, a designed peptide comprising an α-helix buttressed by a polyproline II helix. PPα-Tyr is water soluble and monomeric, and it unfolds cooperatively with a midpoint unfolding temperature (TM) of 39 °C. NMR structures of PPα-Tyr reveal proline residues docked between tyrosine side chains, as designed. The stability of PPα is sensitive to modifications in the aromatic residues: replacing tyrosine with phenylalanine, i.e., changing three solvent-exposed hydroxyl groups to protons, reduces the TM to 20 °C. We attribute this result to the loss of CH-π interactions between the aromatic and proline rings, which we probe by substituting the aromatic residues with nonproteinogenic side chains. In analyses of natural protein structures, we find a preference for proline-tyrosine interactions over other proline-containing pairs, and observe abundant CH-π interactions in biologically important complexes between proline-rich ligands and SH3 and similar domains.

ISAMBARD: an open-source computational environment for biomolecular analysis, modelling and design.

MOTIVATION: The rational design of biomolecules is becoming a reality. However, further computational tools are needed to facilitate and accelerate this, and to make it accessible to more users. RESULTS: Here we introduce ISAMBARD, a tool for structural analysis, model building and rational design of biomolecules. ISAMBARD is open-source, modular, computationally scalable and intuitive to use. These features allow non-experts to explore biomolecular design in silico. ISAMBARD addresses a standing issue in protein design, namely, how to introduce backbone variability in a controlled manner. This is achieved through the generalization of tools for parametric modelling, describing the overall shape of proteins geometrically, and without input from experimentally determined structures. This will allow backbone conformations for entire folds and assemblies not observed in nature to be generated de novo, that is, to access the 'dark matter of protein-fold space'. We anticipate that ISAMBARD will find broad applications in biomolecular design, biotechnology and synthetic biology. AVAILABILITY AND IMPLEMENTATION: A current stable build can be downloaded from the python package index (https://pypi.python.org/pypi/isambard/) with development builds available on GitHub (https://github.com/woolfson-group/) along with documentation, tutorial material and all the scripts used to generate the data described in this paper. CONTACT: d.n.woolfson@bristol.ac.uk or chris.wood@bristol.ac.uk. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Miniprotein Design: Past, Present, and Prospects.

The design and study of miniproteins, that is, polypeptide chains <40 amino acids in length that adopt defined and stable 3D structures, is resurgent. Miniproteins offer possibilities for reducing the complexity of larger proteins and so present new routes to studying sequence-to-structure and sequence-to-stability relationships in proteins generally. They also provide modules for protein design by pieces and, with this, prospects for building more-complex or even entirely new protein structures. In addition, miniproteins are useful scaffolds for templating functional domains, for example, those involved in protein-protein interactions, catalysis, and biomolecular binding, leading to potential applications in biotechnology and medicine. Here we select examples from almost four decades of miniprotein design, development, and dissection. Simply because of the word limit for this Account, we focus on miniproteins that are cooperatively folded monomers in solution and not stabilized by cross-linking or metal binding. In these cases, the optimization of noncovalent interactions is even more critical for the maintenance of the folded states than in larger proteins. Our chronology and catalogue highlights themes in miniproteins, which we explore further and begin to put on a firmer footing through an analysis of the miniprotein structures that have been deposited in the Protein Data Bank (PDB) thus far. Specifically, and compared with larger proteins, miniproteins generally have a lower proportion of residues in regular secondary structure elements (α helices, ß strands, and polyproline-II helices) and, concomitantly, more residues in well-structured loops. This allows distortions of the backbone enabling mini-hydrophobic cores to be made. This also contrasts with larger proteins, which can achieve hydrophobic cores through tertiary contacts between distant regions of sequence. On average, miniproteins have a higher proportion of aromatic residues than larger proteins, and specifically electron-rich Trp and Tyr, which are often found in combination with Pro and Arg to render networks of CH-π or cation-π interactions. Miniproteins also have a higher proportion of the long-chain charged amino acids (Arg, Glu, and Lys), which presumably reflects salt-bridge formation and their greater surface area-to-volume ratio. Together, these amino-acid preferences appear to support greater densities of noncovalent interactions in miniproteins compared with larger proteins. We anticipate that with recent developments such as parametric protein design, it will become increasingly routine to use computation to generate and evaluate models for miniproteins in silico ahead of experimental studies. This could include accessing new structures comprising secondary structure elements linked in previously unseen configurations. The improved understanding of the noncovalent interactions that stabilize the folded states of such miniproteins that we are witnessing through both in-depth bioinformatics analyses and experimental testing will feed these computational protein designs. With this in mind, we can expect a new and exciting era for miniprotein design, study, and application.

Conformational Dynamics of Asparagine at Coiled-Coil Interfaces.

Coiled coils (CCs) are among the best-understood protein folds. Nonetheless, there are gaps in our knowledge of CCs. Notably, CCs are likely to be structurally more dynamic than often considered. Here, we explore this in an abundant class of CCs, parallel dimers, focusing on polar asparagine (Asn) residues in the hydrophobic interface. It is well documented that such inclusions discriminate between different CC oligomers, which has been rationalized in terms of whether the Asn can make side-chain hydrogen bonds. Analysis of parallel CC dimers in the Protein Data Bank reveals a variety of Asn side-chain conformations, but not all of these make the expected inter-side-chain hydrogen bond. We probe the structure and dynamics of a de novo-designed coiled-coil homodimer, CC-Di, by multidimensional nuclear magnetic resonance spectroscopy, including model-free dynamical analysis and relaxation-dispersion experiments. We find dynamic exchange on the millisecond time scale between Asn conformers with the side chains pointing into and out of the core. We perform molecular dynamics simulations that are consistent with this, revealing that the side chains are highly dynamic, exchanging between hydrogen-bonded-paired conformations in picoseconds to nanoseconds. Combined, our data present a more dynamic view for Asn at CC interfaces. Although inter-side-chain hydrogen bonding states are the most abundant, Asn is not always buried or engaged in such interactions. Because interfacial Asn residues are key design features for modulating CC stability and recognition, these further insights into how they are accommodated within CC structures will aid their predictive modeling, engineering, and design.

Toward a Soluble Model System for the Amyloid State.

The formation and deposition of amyloids is associated with many diseases. ß-Sheet secondary structure is a common feature of amyloids, but the packing of sheets against one another is distinctive relative to soluble proteins. Standard methods that rely on perturbing a polypeptide's sequence and evaluating impact on folding can be problematic for amyloid aggregates because a single sequence can adopt multiple conformations and diverse packing arrangements. We describe initial steps toward a minimum-sized, soluble model system for the amyloid state that supports comparisons among sequence variants. Critical to this goal is development of a new linking strategy to enable intersheet association mediated by side chain interactions, which is characteristic of the amyloid state. The linker design we identified should ultimately support exploration of relationships between sequence and amyloid state stability for specific strand-association modes.

Local and macroscopic electrostatic interactions in single α-helices.

The noncovalent forces that stabilize protein structures are not fully understood. One way to address this is to study equilibria between unfolded states and α-helices in peptides. Electrostatic forces-which include interactions between side chains, the backbone and side chains, and side chains and the helix macrodipole-are believed to contribute to these equilibria. Here we probe these interactions experimentally using designed peptides. We find that both terminal backbone-side chain and certain side chain-side chain interactions (which include both local effects between proximal charges and interatomic contacts) contribute much more to helix stability than side chain-helix macrodipole electrostatics, which are believed to operate at larger distances. This has implications for current descriptions of helix stability, the understanding of protein folding and the refinement of force fields for biomolecular modeling and simulations. In addition, this study sheds light on the stability of rod-like structures formed by single α-helices, which are common in natural proteins such as non-muscle myosins.

Carbohydrate-Aromatic Interactions in Proteins.

Protein-carbohydrate interactions play pivotal roles in health and disease. However, defining and manipulating these interactions has been hindered by an incomplete understanding of the underlying fundamental forces. To elucidate common and discriminating features in carbohydrate recognition, we have analyzed quantitatively X-ray crystal structures of proteins with noncovalently bound carbohydrates. Within the carbohydrate-binding pockets, aliphatic hydrophobic residues are disfavored, whereas aromatic side chains are enriched. The greatest preference is for tryptophan with an increased prevalence of 9-fold. Variations in the spatial orientation of amino acids around different monosaccharides indicate specific carbohydrate C-H bonds interact preferentially with aromatic residues. These preferences are consistent with the electronic properties of both the carbohydrate C-H bonds and the aromatic residues. Those carbohydrates that present patches of electropositive saccharide C-H bonds engage more often in CH-π interactions involving electron-rich aromatic partners. These electronic effects are also manifested when carbohydrate-aromatic interactions are monitored in solution: NMR analysis indicates that indole favorably binds to electron-poor C-H bonds of model carbohydrates, and a clear linear free energy relationships with substituted indoles supports the importance of complementary electronic effects in driving protein-carbohydrate interactions. Together, our data indicate that electrostatic and electronic complementarity between carbohydrates and aromatic residues play key roles in driving protein-carbohydrate complexation. Moreover, these weak noncovalent interactions influence which saccharide residues bind to proteins, and how they are positioned within carbohydrate-binding sites.

CCBuilder: an interactive web-based tool for building, designing and assessing coiled-coil protein assemblies.

MOTIVATION: The ability to accurately model protein structures at the atomistic level underpins efforts to understand protein folding, to engineer natural proteins predictably and to design proteins de novo. Homology-based methods are well established and produce impressive results. However, these are limited to structures presented by and resolved for natural proteins. Addressing this problem more widely and deriving truly ab initio models requires mathematical descriptions for protein folds; the means to decorate these with natural, engineered or de novo sequences; and methods to score the resulting models. RESULTS: We present CCBuilder, a web-based application that tackles the problem for a defined but large class of protein structure, the α-helical coiled coils. CCBuilder generates coiled-coil backbones, builds side chains onto these frameworks and provides a range of metrics to measure the quality of the models. Its straightforward graphical user interface provides broad functionality that allows users to build and assess models, in which helix geometry, coiled-coil architecture and topology and protein sequence can be varied rapidly. We demonstrate the utility of CCBuilder by assembling models for 653 coiled-coil structures from the PDB, which cover >96% of the known coiled-coil types, and by generating models for rarer and de novo coiled-coil structures. AVAILABILITY AND IMPLEMENTATION: CCBuilder is freely available, without registration, at http://coiledcoils.chm.bris.ac.uk/app/cc_builder/.

Interplay of hydrogen bonds and n→π* interactions in proteins.

Protein structures are stabilized by multiple weak interactions, including the hydrophobic effect, hydrogen bonds, electrostatic effects, and van der Waals interactions. Among these interactions, the hydrogen bond is distinct in having its origins in electron delocalization. Recently, another type of electron delocalization, the n→π* interaction between carbonyl groups, has been shown to play a role in stabilizing protein structure. Here we examine the interplay between hydrogen bonding and n→π* interactions. To address this issue, we used data available from high-resolution protein crystal structures to interrogate asparagine side-chain oxygen atoms that are both acceptors of a hydrogen bond and donors of an n→π* interaction. Then we employed natural bond orbital analysis to determine the relative energetic contributions of the hydrogen bonds and n→π* interactions in these systems. We found that an n→π* interaction is worth ~5-25% of a hydrogen bond and that stronger hydrogen bonds tend to attenuate or obscure n→π* interactions. Conversely, weaker hydrogen bonds correlate with stronger n→π* interactions and demixing of the orbitals occupied by the oxygen lone pairs. Thus, these two interactions conspire to stabilize local backbone-side-chain contacts, which argues for the inclusion of n→π* interactions in the inventory of non-covalent forces that contribute to protein stability and thus in force fields for biomolecular modeling.

A de novo peptide hexamer with a mutable channel.

The design of new proteins that expand the repertoire of natural protein structures represents a formidable challenge. Success in this area would increase understanding of protein structure and present new scaffolds that could be exploited in biotechnology and synthetic biology. Here we describe the design, characterization and X-ray crystal structure of a new coiled-coil protein. The de novo sequence forms a stand-alone, parallel, six-helix bundle with a channel running through it. Although lined exclusively by hydrophobic leucine and isoleucine side chains, the 6-Å channel is permeable to water. One layer of leucine residues within the channel is mutable, accepting polar aspartic acid and histidine side chains, which leads to subdivision and organization of solvent within the lumen. Moreover, these mutants can be combined to form a stable and unique (Asp-His)(3) heterohexamer. These new structures provide a basis for engineering de novo proteins with new functions.

Designed coiled coils promote folding of a recombinant bacterial collagen.

Collagen triple helices fold slowly and inefficiently, often requiring adjacent globular domains to assist this process. In the Streptococcus pyogenes collagen-like protein Scl2, a V domain predicted to be largely α-helical, occurs N-terminal to the collagen triple helix (CL). Here, we replace this natural trimerization domain with a de novo designed, hyperstable, parallel, three-stranded, α-helical coiled coil (CC), either at the N terminus (CC-CL) or the C terminus (CL-CC) of the collagen domain. CD spectra of the constructs are consistent with additivity of independently and fully folded CC and CL domains, and the proteins retain their distinctive thermal stabilities, CL at ∼37 °C and CC at >90 °C. Heating the hybrid proteins to 50 °C unfolds CL, leaving CC intact, and upon cooling, the rate of CL refolding is somewhat faster for CL-CC than for CC-CL. A construct with coiled coils on both ends, CC-CL-CC, retains the ∼37 °C thermal stability for CL but shows less triple helix at low temperature and less denaturation at 50 °C. Most strikingly however, in CC-CL-CC, the CL refolds slower than in either CC-CL or CL-CC by almost two orders of magnitude. We propose that a single CC promotes folding of the CL domain via nucleation and in-register growth from one end, whereas initiation and growth from both ends in CC-CL-CC results in mismatched registers that frustrate folding. Bioinformatics analysis of natural collagens lends support to this because, where present, there is generally only one coiled-coil domain close to the triple helix, and it is nearly always N-terminal to the collagen repeat.

Strong contributions from vertical triads to helix-partner preferences in parallel coiled coils.

Pairing preferences in heterodimeric coiled coils are determined by complementarities among side chains that pack against one another at the helix-helix interface. However, relationships between dimer stability and interfacial residue identity are not fully understood. In the context of the "knobs-into-holes" (KIH) packing pattern, one can identify two classes of interactions between side chains from different helices: "lateral", in which a line connecting the adjacent side chains is perpendicular to the helix axes, and "vertical", in which the connecting line is parallel to the helix axes. We have previously analyzed vertical interactions in antiparallel coiled coils and found that one type of triad constellation (a'-a-a') exerts a strong effect on pairing preferences, while the other type of triad (d'-d-d') has relatively little impact on pairing tendencies. Here, we ask whether vertical interactions (d'-a-d') influence pairing in parallel coiled-coil dimers. Our results indicate that vertical interactions can exert a substantial impact on pairing specificity, and that the influence of the d'-a-d' triad depends on the lateral a' contact within the local KIH motif. Structure-informed bioinformatic analyses of protein sequences reveal trends consistent with the thermodynamic data derived from our experimental model system in suggesting that heterotriads involving Leu and Ile are preferred over homotriads involving Leu and Ile.

The d'--d--d' vertical triad is less discriminating than the a'--a--a' vertical triad in the antiparallel coiled-coil dimer motif.

Elucidating relationships between the amino-acid sequences of proteins and their three-dimensional structures, and uncovering non-covalent interactions that underlie polypeptide folding, are major goals in protein science. One approach toward these goals is to study interactions between selected residues, or among constellations of residues, in small folding motifs. The α-helical coiled coil has served as a platform for such studies because this folding unit is relatively simple in terms of both sequence and structure. Amino acid side chains at the helix-helix interface of a coiled coil participate in so-called "knobs-into-holes" (KIH) packing whereby a side chain (the knob) on one helix inserts into a space (the hole) generated by four side chains on a partner helix. The vast majority of sequence-stability studies on coiled-coil dimers have focused on lateral interactions within these KIH arrangements, for example, between an a position on one helix and an a' position of the partner in a parallel coiled-coil dimer, or between a--d' pairs in an antiparallel dimer. More recently, it has been shown that vertical triads (specifically, a'--a--a' triads) in antiparallel dimers exert a significant impact on pairing preferences. This observation provides impetus for analysis of other complex networks of side-chain interactions at the helix-helix interface. Here, we describe a combination of experimental and bioinformatics studies that show that d'--d--d' triads have much less impact on pairing preference than do a'--a--a' triads in a small, designed antiparallel coiled-coil dimer. However, the influence of the d'--d--d' triad depends on the lateral a'--d interaction. Taken together, these results strengthen the emerging understanding that simple pairwise interactions are not sufficient to describe side-chain interactions and overall stability in antiparallel coiled-coil dimers; higher-order interactions must be considered as well.

Squaring the circle in peptide assembly: from fibers to discrete nanostructures by de novo design.

The design of bioinspired nanostructures and materials of defined size and shape is challenging as it pushes our understanding of biomolecular assembly to its limits. In such endeavors, DNA is the current building block of choice because of its predictable and programmable self-assembly. The use of peptide- and protein-based systems, however, has potential advantages due to their more-varied chemistries, structures and functions, and the prospects for recombinant production through gene synthesis and expression. Here, we present the design and characterization of two complementary peptides programmed to form a parallel heterodimeric coiled coil, which we use as the building blocks for larger, supramolecular assemblies. To achieve the latter, the two peptides are joined via peptidic linkers of variable lengths to produce a range of assemblies, from flexible fibers of indefinite length, through large colloidal-scale assemblies, down to closed and discrete nanoscale objects of defined stoichiometry. We posit that the different modes of assembly reflect the interplay between steric constraints imposed by short linkers and the bulk of the helices, and entropic factors that favor the formation of many smaller objects as the linker length is increased. This approach, and the resulting linear and proteinogenic polypeptides, represents a new route for constructing complex peptide-based assemblies and biomaterials.

n-->pi* interactions in proteins.

Hydrogen bonds between backbone amides are common in folded proteins. Here, we show that an intimate interaction between backbone amides also arises from the delocalization of a lone pair of electrons (n) from an oxygen atom to the antibonding orbital (pi*) of the subsequent carbonyl group. Natural bond orbital analysis predicted significant n-->pi* interactions in certain regions of the Ramachandran plot. These predictions were validated by a statistical analysis of a large, non-redundant subset of protein structures determined to high resolution. The correlation between these two independent studies is striking. Moreover, the n-->pi* interactions are abundant and especially prevalent in common secondary structures such as alpha-, 3(10)- and polyproline II helices and twisted beta-sheets. In addition to their evident effects on protein structure and stability, n-->pi* interactions could have important roles in protein folding and function, and merit inclusion in computational force fields.

Understanding p300-transcription factor interactions using sequence variation and hybridization.

The hypoxic response is central to cell function and plays a significant role in the growth and survival of solid tumours. HIF-1 regulates the hypoxic response by activating over 100 genes responsible for adaptation to hypoxia, making it a potential target for anticancer drug discovery. Although there is significant structural and mechanistic understanding of the interaction between HIF-1α and p300 alongside negative regulators of HIF-1α such as CITED2, there remains a need to further understand the sequence determinants of binding. In this work we use a combination of protein expression, chemical synthesis, fluorescence anisotropy and isothermal titration calorimetry for HIF-1α sequence variants and a HIF-1α-CITED hybrid sequence which we term CITIF. We show the HIF-1α sequence is highly tolerant to sequence variation through reduced enthalpic and less unfavourable entropic contributions, These data imply backbone as opposed to side chain interactions and ligand folding control the binding interaction and that sequence variations are tolerated as a result of adopting a more disordered bound interaction or "fuzzy" complex.

Query-guided protein-protein interaction inhibitor discovery.

Protein-protein interactions (PPIs) are central to biological mechanisms, and can serve as compelling targets for drug discovery. Yet, the discovery of small molecule inhibitors of PPIs remains challenging given the large and typically shallow topography of the interacting protein surfaces. Here, we describe a general approach to the discovery of orthosteric PPI inhibitors that mimic specific secondary protein structures. Initially, hot residues at protein-protein interfaces are identified in silico or from experimental data, and incorporated into secondary structure-based queries. Virtual libraries of small molecules are then shape-matched against the queries, and promising ligands docked to target proteins. The approach is exemplified experimentally using two unrelated PPIs that are mediated by an α-helix (p53/hDM2) and a ß-strand (GKAP/SHANK1-PDZ). In each case, selective PPI inhibitors are discovered with low µM activity as determined by a combination of fluorescence anisotropy and 1H-15N HSQC experiments. In addition, hit expansion yields a series of PPI inhibitors with defined structure-activity relationships. It is envisaged that the generality of the approach will enable discovery of inhibitors of a wide range of unrelated secondary structure-mediated PPIs.