Engineering protein stability with atomic precision in a monomeric miniprotein

Methyl-π Interactions. Nature of Bonding and Limits of Strength

Both methyl groups and benzene rings are exceedingly common, and they lie near one another in many chemical situations. DFT calculations are used to gauge the strength of the attractive forces between them, and to better understand the phenomena that underlie this attraction. Methane and benzene are taken as the starting point, and substituents of both electron-withdrawing and donating types are added to each. The interaction energy varies between 1.4 and 5.0 kcal/mol, depending upon the substituents placed on the two groups. The nature of the binding is analyzed via Atoms in Molecules (AIM), Natural Bond Orbital (NBO), Symmetry-Adapted Perturbation Theory (SAPT), nuclear magnetic resonance (NMR) chemical shifts, and electron density shift diagrams. While there is a sizable electrostatic component, it is dispersion that dominates these interactions, particularly the weaker ones. As such, these interactions cannot be categorized unambiguously as either H-bonds or tetrel bonds.

Stabilization of a miniprotein fold by an unpuckered proline surrogate

The unique role of proline in modulating protein folding and recognition makes it an attractive target for substitution to generate new proteomimetics. The design, synthesis, and conformational analysis of non-canonical surrogates can also aid in parsing the role of prolyl stereoelectronic effects on structure. We recently described the synthesis and conformational analysis of dehydro-δ-azaproline (ΔaPro), a novel unsaturated analogue of proline featuring a planar dehydropyrazine ring. When incorporated into host sequences, this backbone N-aminated proline surrogate forms an acylhydrazone bond with an unusually high trans rotamer bias and low isomerization barrier. Here, we used CD, NMR spectroscopy, and MD simulations to evaluate the impact of ΔaPro substitution within the polyproline II (PPII) and loop regions of the avian pancreatic polypeptide (aPP). The ΔaPro residue strongly favors PPII conformation and stabilizes the aPP tertiary fold when incorporated at select positions within the miniprotein. A variant featuring three ΔaPro substitutions was found to significantly enhance the thermal stability of wild-type aPP despite compromising protein dimerization. Our results suggest that the stability of proline-rich folds relies more on backbone torsional preferences than ring puckering and informs strategies for the incorporation of ΔaPro into thermally stable and functional proteomimetics.

Nanoscale Evolution of Charge Transport Through C-H···π Interactions

C-H···π interactions, a prevalent intermolecular force, play a pivotal role in chemistry, materials science, and life sciences. Despite extensive studies of their influence on intermolecular binding configurations and energetics, their impact on intermolecular coupling and charge transport remains unexplored. Here, we investigate the charge transport within supramolecular junctions connected by C-H···π and π-π interactions, respectively, and find that C-H···π interactions exhibit conductances that are 3.5 times those of π-π interactions. Angstrom-scale distance-dependent experiments indicate that the conductance of C-H···π supramolecular junctions experiences initial decay under stretching, followed by gradual convergence, in contrast with the periodic fluctuations in π-π stacked supramolecular junctions. Theoretical calculations show that charge transport within C-H···π interactions transitions from destructive to constructive quantum interference under stretching, with a larger range of constructive quantum interference compared with π-π stacking. This study establishes that C-H···π interactions facilitate efficient intermolecular charge transport and elucidates the evolution of quantum interference effects with assembly configuration, offering critical insights for the design of supramolecular materials and devices.

Rationally seeded computational protein design of ɑ-helical barrels

Computational protein design is advancing rapidly. Here we describe efficient routes starting from validated parallel and antiparallel peptide assemblies to design two families of α-helical barrel proteins with central channels that bind small molecules. Computational designs are seeded by the sequences and structures of defined de novo oligomeric barrel-forming peptides, and adjacent helices are connected by loop building. For targets with antiparallel helices, short loops are sufficient. However, targets with parallel helices require longer connectors; namely, an outer layer of helix-turn-helix-turn-helix motifs that are packed onto the barrels. Throughout these computational pipelines, residues that define open states of the barrels are maintained. This minimizes sequence sampling, accelerating the design process. For each of six targets, just two to six synthetic genes are made for expression in Escherichia coli. On average, 70% of these genes express to give soluble monomeric proteins that are fully characterized, including high-resolution structures for most targets that match the design models with high accuracy.

Architectonic principles of polyproline II helix bundle protein domains

Glycine rich polyproline II helix assemblies are an emerging class of natural domains found in several proteins with different functions and diverse origins. The distinct properties of these domains relative to those composed of α-helices and β-sheets could make glycine-rich polyproline II helix assemblies a useful building block for protein design. Whereas the high population of polyproline II conformers in disordered state ensembles could facilitate glycine-rich polyproline II helix folding, the architectonic bases of these structures are not well known. Here, we compare and analyze their structures to uncover common features. These protein domains are found to be highly tolerant of distinct flanking sequences. This speaks to the robustness of this fold and strongly suggests that glycine rich polyproline II assemblies could be grafted with other protein domains to engineer new structures and functions. These domains are also well packed with few or no cavities. Moreover, a significant trend towards antiparallel helix configuration is observed in all these domains and could provide stabilizing interactions among macrodipoles. Finally, extensive networks of Cα-H···OC hydrogen bonds are detected in these domains. Despite their diverse evolutionary origins and activities, glycine-rich polyproline II helix assemblies share architectonic features which could help design novel proteins.

Tunable CH/π Interactions within a Tryptophan Zipper Motif to Stabilize the Fold of Long β-Hairpin Peptides

The tryptophan zipper (Trpzip) is an iconic folding motif of β-hairpin peptides capitalizing on two pairs of cross-strand tryptophans, each stabilized by an aromatic-aromatic stacking in an edge-to-face (EtF) geometry. Yet, the origins and the contribution of this EtF packing to the unique Trpzip stability remain poorly understood. To address this question of structure-stability relationship, a library of Trpzip hairpins was developed by incorporating readily accessible nonproteinogenic tryptophans of varying electron densities. We found that each EtF geometry was, in fact, stabilized by an intricate combination of XH/π interactions. By tuning the π-electron density of Trpface rings, CH/π interactions are strengthened to gain additional stability. On the contrary, our DFT calculations support the notion that Trpedge modulations are challenging due to their simultaneous paradoxical engagement as H-bond donors in CH/π and acceptors in NH/π interactions.

The Energetic Origins of Pi-Pi Contacts in Proteins

Accurate potential energy models of proteins must describe the many different types of noncovalent interactions that contribute to a protein's stability and structure. Pi-pi contacts are ubiquitous structural motifs in all proteins, occurring between aromatic and nonaromatic residues and play a nontrivial role in protein folding and in the formation of biomolecular condensates. Guided by a geometric criterion for isolating pi-pi contacts from classical molecular dynamics simulations of proteins, we use quantum mechanical energy decomposition analysis to determine the molecular interactions that stabilize different pi-pi contact motifs. We find that neutral pi-pi interactions in proteins are dominated by Pauli repulsion and London dispersion rather than repulsive quadrupole electrostatics, which is central to the textbook Hunter-Sanders model. This results in a notable lack of variability in the interaction profiles of neutral pi-pi contacts even with extreme changes in the dielectric medium, explaining the prevalence of pi-stacked arrangements in and between proteins. We also find interactions involving pi-containing anions and cations to be extremely malleable, interacting like neutral pi-pi contacts in polar media and like typical ion-pi interactions in nonpolar environments. Like-charged pairs such as arginine-arginine contacts are particularly sensitive to the polarity of their immediate surroundings and exhibit canonical pi-pi stacking behavior only if the interaction is mediated by environmental effects, such as aqueous solvation.

AQcalc: A web server that identifies weak molecular interactions in protein structures

Weak molecular interactions play an important role in protein structure and function. Computational tools that identify weak molecular interactions are, therefore, valuable for the study of proteins. Here, we present AQcalc, a web server (https://aqcalcbiocomputing.com/) that can be used to identify anion-quadrupole (AQ) interactions, which are weak interactions involving aromatic residue (Trp, Tyr, and Phe) ring edges and anions (Asp, Glu, and phosphate ion) both within proteins and at their interfaces (protein-protein, protein-nucleic acids, and protein-lipid bilayer). AQcalc identifies AQ interactions as well as clusters involving AQ, cation-π, and salt bridges, among others. Utilizing AQcalc we analyzed weak interactions in protein models, even in the absence of experimental structures, to understand the contributions of weak interactions to deleterious structural changes, including those associated with oncogenic and germline disease variants. We identified several deleterious variants with disrupted AQ interactions (comparable in frequency to cation-π disruptions). Amyloid fibrils utilize AQ to bury anions at frequencies that far exceed those observed for globular proteins. AQ interactions were detected three and five times more frequently than the hydrogen-bonded AQ (HBAQ) in fibril structures and protein-lipid bilayer interfaces, respectively. By contrast, AQ and HBAQ interactions were detected with similar frequencies in globular proteins. Collectively, these findings suggest AQcalc will be effective in facilitating fine structural analysis. As other web utilities designed to identify protein residue interaction networks do not report AQ interactions, wide use of AQcalc will enrich our understanding of residue interaction networks and facilitate hypothesis testing by identifying and experimentally characterizing these comparably weak but important interactions.

Visualization Tools for Detecting Microcystin-LR in the Biological System via Near-Infrared Fluorescent Probes

Numerous toxicological and epidemiological studies have shown that microcystin-LR (MC-LR) could cause a variety of toxicity to humans and animals. However, the absence of effective methods to trace MC-LR in biological systems has hindered the in-depth understanding of the mechanism of MC-LR toxicity. Near-infrared (NIR) fluorescent probes are crucial tools for accurate visualization and in-depth study of specific molecules in biological systems. Due to the lack of effective design strategies, NIR fluorescent probes for imaging MC-LR specifically in biological systems have not been reported yet. In order to address this pressing issue, herein, we have introduced a new and facile strategy to improve MC-LR detection and imaging in biological systems, and based on this design strategy, three NIR fluorescence probes (MC-RdTPA1, MC-RdTPA2, and MC-RdTPE1) have been constructed. These probes have several advantages: (i) have long emission wavelength and large Stokes shifts, which have great potential in vivo imaging applications; (ii) could selectively visualize MC-LR in cells; and (iii) showed stable fluorescence intensity in the pH range of 5.0-7.0. This work may provide a new avenue for the detection of MC-LR in biological systems and new tool to advance our knowledge of the mechanism of MC-LR toxicity.

The Realm of Unconventional Noncovalent Interactions in Proteins: Their Significance in Structure and Function

Proteins and their assemblies are fundamental for living cells to function. Their complex three-dimensional architecture and its stability are attributed to the combined effect of various noncovalent interactions. It is critical to scrutinize these noncovalent interactions to understand their role in the energy landscape in folding, catalysis, and molecular recognition. This Review presents a comprehensive summary of unconventional noncovalent interactions, beyond conventional hydrogen bonds and hydrophobic interactions, which have gained prominence over the past decade. The noncovalent interactions discussed include low-barrier hydrogen bonds, C5 hydrogen bonds, C-H···π interactions, sulfur-mediated hydrogen bonds, n → π* interactions, London dispersion interactions, halogen bonds, chalcogen bonds, and tetrel bonds. This Review focuses on their chemical nature, interaction strength, and geometrical parameters obtained from X-ray crystallography, spectroscopy, bioinformatics, and computational chemistry. Also highlighted are their occurrence in proteins or their complexes and recent advances made toward understanding their role in biomolecular structure and function. Probing the chemical diversity of these interactions, we determined that the variable frequency of occurrence in proteins and the ability to synergize with one another are important not only for ab initio structure prediction but also to design proteins with new functionalities. A better understanding of these interactions will promote their utilization in designing and engineering ligands with potential therapeutic value.

Protein Backbone Alteration in Non-Hairpin β-Turns: Impacts on Tertiary Folded Structure and Folded Stability

The importance of β-turns to protein folding has motivated extensive efforts to stabilize the motif with non-canonical backbone connectivity. Prior work has focused almost exclusively on turns between strands in a β-sheet (i. e., hairpins). Turns in other structural contexts are also common in nature and have distinct conformational preferences; however, design principles for their mimicry remain poorly understood. Here, we report strategies that stabilize non-hairpin β-turns through systematic evaluation of the impacts of backbone alteration on the high-resolution folded structure and folded stability of a helix-loop-helix prototype protein. Several well-established hairpin turn mimetics are shown detrimental to folded stability and/or hydrophobic core packing, while less-explored modification schemes that reinforce alternate turn types lead to improved stability and more faithful structural mimicry. Collectively, these results have implications in control over protein folding through chemical modification as well as the design of protein mimetics.

Understanding a protein fold: the physics, chemistry, and biology of α-helical coiled coils

Protein science is being transformed by powerful computational methods for structure prediction and design: AlphaFold2 can predict many natural protein structures from sequence, and other AI methods are enabling the de novo design of new structures. This raises a question: how much do we understand the underlying sequence-to-structure/function relationships being captured by these methods? This perspective presents our current understanding of one class of protein assembly, the α-helical coiled coils. At first sight, these are straightforward: sequence repeats of hydrophobic (h) and polar (p) residues, (hpphppp)_n, direct the folding and assembly of amphipathic α helices into bundles. However, many different bundles are possible: they can have two or more helices (different oligomers); the helices can have parallel, antiparallel or mixed arrangements (different topologies); and the helical sequences can be the same (homomers) or different (heteromers). Thus, sequence-to-structure relationships must be present within the hpphppp repeats to distinguish these states. I discuss the current understanding of this problem at three levels: First, physics gives a parametric framework to generate the many possible coiled-coil backbone structures. Second, chemistry provides a means to explore and deliver sequence-to-structure relationships. Third, biology shows how coiled coils are adapted and functionalized in nature, inspiring applications of coiled coils in synthetic biology. I argue that the chemistry is largely understood; the physics is partly solved, though the considerable challenge of predicting even relative stabilities of different coiled-coil states remains; but there is much more to explore in the biology and synthetic biology of coiled coils.

Dissecting the stability determinants of a challenging de novo protein fold using massively parallel design and experimentation

Designing entirely new protein structures remains challenging because we do not fully understand the biophysical determinants of folding stability. Yet, some protein folds are easier to design than others. Previous work identified the 43-residue ɑββɑ fold as especially challenging: The best designs had only a 2% success rate, compared to 39 to 87% success for other simple folds [G. J. Rocklin et al., Science 357, 168-175 (2017)]. This suggested the ɑββɑ fold would be a useful model system for gaining a deeper understanding of folding stability determinants and for testing new protein design methods. Here, we designed over 10,000 new ɑββɑ proteins and found over 3,000 of them to fold into stable structures using a high-throughput protease-based assay. NMR, hydrogen-deuterium exchange, circular dichroism, deep mutational scanning, and scrambled sequence control experiments indicated that our stable designs fold into their designed ɑββɑ structures with exceptional stability for their small size. Our large dataset enabled us to quantify the influence of universal stability determinants including nonpolar burial, helix capping, and buried unsatisfied polar atoms, as well as stability determinants unique to the ɑββɑ topology. Our work demonstrates how large-scale design and test cycles can solve challenging design problems while illuminating the biophysical determinants of folding.

Miniproteins in medicinal chemistry

Miniproteins exhibit great potential as scaffolds for drug candidates because of their well-defined structure and good synthetic availability. Because of recently described methodologies for their de novo design, the field of miniproteins is emerging and can provide molecules that effectively bind to problematic targets, i.e., those that have been previously considered to be undruggable. This review describes methodologies for the development of miniprotein scaffolds and for the construction of biologically active miniproteins.

Structure-guided design of CPPC-paired disulfide-rich peptide libraries for ligand and drug discovery

Peptides constrained through multiple disulfides (or disulfide-rich peptides, DRPs) have been an emerging frontier for ligand and drug discovery. Such peptides have the potential to combine the binding capability of biologics with the stability and bioavailability of smaller molecules. However, DRPs with stable three-dimensional (3D) structures are usually of natural origin or engineered from natural ones. Here, we report the discovery and identification of CPPC (cysteine-proline-proline-cysteine) motif-directed DRPs with stable 3D structures (i.e., CPPC-DRPs). A range of new CPPC-DRPs were designed or selected from either random or structure-convergent peptide libraries. Thus, for the first time we revealed that the CPPC-DRPs can maintain diverse 3D structures by taking advantage of constraints from unique dimeric CPPC mini-loops, including irregular structures and regular α-helix and β-sheet folds. New CPPC-DRPs that can specifically bind the receptors (CD28) on the cell surface were also successfully discovered and identified using our DRP-discovery platform. Overall, this study provides the basis for accessing an unconventional peptide structure space previously inaccessible by natural DRPs and computational designs, inspiring the development of new peptide ligands and therapeutics.

Transferable and Polarizable Coarse Grained Model for Proteins─ProMPT

The application of classical molecular dynamics (MD) simulations at atomic resolution (fine-grained level, FG), to most biomolecular processes, remains limited because of the associated computational complexity of representing all the atoms. This problem is magnified in the presence of protein-based biomolecular systems that have a very large conformational space, and MD simulations with fine-grained resolution have slow dynamics to explore this space. Current transferable coarse grained (CG) force fields in literature are either limited to only peptides with the environment encoded in an implicit form or cannot capture transitions into secondary/tertiary peptide structures from a primary sequence of amino acids. In this work, we present a transferable CG force field with an explicit representation of the environment for accurate simulations with proteins. The force field consists of a set of pseudoatoms representing different chemical groups that can be joined/associated together to create different biomolecular systems. This preserves the transferability of the force field to multiple environments and simulation conditions. We have added electronic polarization that can respond to environmental heterogeneity/fluctuations and couple it to protein's structural transitions. The nonbonded interactions are parametrized with physics-based features such as solvation and partitioning free energies determined by thermodynamic calculations and matched with experiments and/or atomistic simulations. The bonded potentials are inferred from corresponding distributions in nonredundant protein structure databases. We present validations of the CG model with simulations of well-studied aqueous protein systems with specific protein fold types─Trp-cage, Trpzip4, villin, WW-domain, and β-α-β. We also explore the applications of the force field to study aqueous aggregation of Aβ 16-22 peptides.

Miniature β-Hairpin Mimetic by Intramolecular Hydrogen Bond and C-H···π Interactions

Canonically, protein β-hairpin motifs are stabilized by intramolecular hydrogen bonds. Here, we attempt to develop a rational design recipe for a miniature hairpin structure stabilized by hydrogen bonding as well as C-H···π interaction and try to understand how such a stabilization effect varies with different functional groups at each terminus. Database analysis shows that the α-amino acids with an aromatic side chain will not favor that kind of C-H···π stabilized hairpin structure. However, hybrid tripeptides with an N-terminal Boc-Trp-Aib corner residue and C-terminal aromatic ω-amino acids fold into the hairpin conformation with a central β-turn/open-turn that is reinforced by a C-H···π interaction. The CCDC database analysis further confirms that this C-H···π stabilized hairpin motif is general for Boc-protected tripeptides containing Aib in the middle and aromatic functionality at the C-terminus. The different α-amino acids like Leu/Ala/Phe/Pro/Ser at the N-terminus have a minor influence on the C-H···π interaction and stabilities of the folded structures in solid-state. However, the hybrid peptides exhibit different degrees of conformational heterogeneity both in the solid and solution phase, which is common for this kind of flexible small molecule. Conformational heterogeneity in the solution phase including the C-H···π stabilized β-hairpin structures are characterized by the molecular dynamics (MD) simulations explaining their plausible origin at an atomistic level.

Target highlights in CASP14: Analysis of models by structure providers

The biological and functional significance of selected Critical Assessment of Techniques for Protein Structure Prediction 14 (CASP14) targets are described by the authors of the structures. The authors highlight the most relevant features of the target proteins and discuss how well these features were reproduced in the respective submitted predictions. The overall ability to predict three-dimensional structures of proteins has improved remarkably in CASP14, and many difficult targets were modeled with impressive accuracy. For the first time in the history of CASP, the experimentalists not only highlighted that computational models can accurately reproduce the most critical structural features observed in their targets, but also envisaged that models could serve as a guidance for further studies of biologically-relevant properties of proteins.

Structural basis for the hyperthermostability of an archaeal enzyme induced by succinimide formation

Stability of proteins from hyperthermophiles (organisms existing under boiling water conditions) enabled by a reduction of conformational flexibility is realized through various mechanisms. A succinimide (SNN) arising from the post-translational cyclization of the side chains of aspartyl/asparaginyl residues with the backbone amide -NH of the succeeding residue would restrain the torsion angle Ψ and can serve as a new route for hyperthermostability. However, such a succinimide is typically prone to hydrolysis, transforming to either an aspartyl or β-isoaspartyl residue. Here, we present the crystal structure of Methanocaldococcus jannaschii glutamine amidotransferase and, using enhanced sampling molecular dynamics simulations, address the mechanism of its increased thermostability, up to 100°C, imparted by an unexpectedly stable succinimidyl residue at position 109. The stability of SNN109 to hydrolysis is seen to arise from its electrostatic shielding by the side-chain carboxylate group of its succeeding residue Asp110, as well as through n → π^∗ interactions between SNN109 and its preceding residue Glu108, both of which prevent water access to SNN. The stable succinimidyl residue induces the formation of an α-turn structure involving 13-atom hydrogen bonding, which locks the local conformation, reducing protein flexibility. The destabilization of the protein upon replacement of SNN with a Φ-restricted prolyl residue highlights the specificity of the succinimidyl residue in imparting hyperthermostability to the enzyme. The conservation of the succinimide-forming tripeptide sequence (E(N/D)(E/D)) in several archaeal GATases strongly suggests an adaptation of this otherwise detrimental post-translational modification as a harbinger of thermostability.

A Brief History of De Novo Protein Design: Minimal, Rational, and Computational

Protein design has come of age, but how will it mature? In the 1980s and the 1990s, the primary motivation for de novo protein design was to test our understanding of the informational aspect of the protein-folding problem; i.e., how does protein sequence determine protein structure and function? This necessitated minimal and rational design approaches whereby the placement of each residue in a design was reasoned using chemical principles and/or biochemical knowledge. At that time, though with some notable exceptions, the use of computers to aid design was not widespread. Over the past two decades, the tables have turned and computational protein design is firmly established. Here, I illustrate this progress through a timeline of de novo protein structures that have been solved to atomic resolution and deposited in the Protein Data Bank. From this, it is clear that the impact of rational and computational design has been considerable: More-complex and more-sophisticated designs are being targeted with many being resolved to atomic resolution. Furthermore, our ability to generate and manipulate synthetic proteins has advanced to a point where they are providing realistic alternatives to natural protein functions for applications both in vitro and in cells. Also, and increasingly, computational protein design is becoming accessible to non-specialists. This all begs the questions: Is there still a place for minimal and rational design approaches? And, what challenges lie ahead for the burgeoning field of de novo protein design as a whole?

Computer-aided understanding and engineering of enzymatic selectivity

Enzymes offering chemo-, regio-, and stereoselectivity enable the asymmetric synthesis of high-value chiral molecules. Unfortunately, the drawback that naturally occurring enzymes are often inefficient or have undesired selectivity toward non-native substrates hinders the broadening of biocatalytic applications. To match the demands of specific selectivity in asymmetric synthesis, biochemists have implemented various computer-aided strategies in understanding and engineering enzymatic selectivity, diversifying the available repository of artificial enzymes. Here, given that the entire asymmetric catalytic cycle, involving precise interactions within the active pocket and substrate transport in the enzyme channel, could affect the enzymatic efficiency and selectivity, we presented a comprehensive overview of the computer-aided workflow for enzymatic selectivity. This review includes a mechanistic understanding of enzymatic selectivity based on quantum mechanical calculations, rational design of enzymatic selectivity guided by enzyme-substrate interactions, and enzymatic selectivity regulation via enzyme channel engineering. Finally, we discussed the computational paradigm for designing enzyme selectivity in silico to facilitate the advancement of asymmetric biosynthesis.

Biomolecular Modeling and Simulation: A Prospering Multidisciplinary Field

We reassess progress in the field of biomolecular modeling and simulation, following up on our perspective published in 2011. By reviewing metrics for the field's productivity and providing examples of success, we underscore the productive phase of the field, whose short-term expectations were overestimated and long-term effects underestimated. Such successes include prediction of structures and mechanisms; generation of new insights into biomolecular activity; and thriving collaborations between modeling and experimentation, including experiments driven by modeling. We also discuss the impact of field exercises and web games on the field's progress. Overall, we note tremendous success by the biomolecular modeling community in utilization of computer power; improvement in force fields; and development and application of new algorithms, notably machine learning and artificial intelligence. The combined advances are enhancing the accuracy andscope of modeling and simulation, establishing an exemplary discipline where experiment and theory or simulations are full partners.

Dynamic Community Composition Unravels Allosteric Communication in PDZ3

The third domain of PSD-95 (PDZ3) is a model for investigating allosteric communication in protein and ligand interactions. While motifs contributing to its binding specificity have been scrutinized, a conformational dynamical basis is yet to be established. Despite the miniscule structural changes due to point mutants, the observed significant binding affinity differences have previously been assessed with a focus on two α-helices located at the binding groove (α₂) and the C-terminus (α₃). Here, we employ a new computational approach to develop a generalized view on the molecular basis of PDZ3 binding selectivity and interaction communication for a set of point mutants of the protein (G330T, H372A, G330T-H372A) and its ligand (CRIPT, named L₁, and its T-2F variant, L₂) along with the wild type (WT). To analyze the dynamical aspects hidden in the conformations that are produced by molecular dynamics simulations, we utilize variations in community composition calculated based on the betweenness centrality measure from graph theory. We find that the highly charged N-terminus, which is located far from the ligand, has the propensity to share the same community with the ligand in the biologically functional complexes, indicating a distal segment might mediate the binding dynamics. N- and C-termini of PDZ3 share communities, and α₃ acts as a hub for the whole protein by sustaining the communication with all structural segments, albeit being a trait not unique to the functional complexes. Moreover, α₂ which lines the binding cavity frequently parts communities with the ligand and is not a controller of the binding but is rather a slave to the overall dynamics coordinated by the N-terminus. Thus, ligand binding fate in PDZ3 is traced to the population of community compositions extracted from dynamics despite the lack of significant conformational changes.

CH-π interaction between cross-strand amino acid pairs stabilizes β-hairpins

We identified several CH-π donor-acceptor pairs involving amino acid side chains with less polarized C-H bonds at a solvent-exposed site between the strands of a β-hairpin peptide. Therein, we observe a distance-dependent induction of CH-π interaction within the aliphatic-aromatic amino acid pair. Our results also suggest an interplay of hydrophobicity and CH-π interaction in dictating the stability of β-hairpins, where a leucine-tryptophan pair contributes -1.14 kcal mol-1 to the overall foldedness of the β-hairpin.

More Than π-π-π Stacking: Contribution of Amide-π and CH-π Interactions to Crotonyllysine Binding by the AF9 YEATS Domain

Lysine crotonylation (Kcr) is a histone post-translational modification that is implicated in numerous epigenetic pathways and diseases. Recognition of Kcr by YEATS domains has been proposed to occur through intermolecular amide-π and alkene-π interactions, but little is known about the driving force of these key interactions. Herein, we probed the recognition of lysine crotonylation and acetylation by the AF9 YEATS domain through incorporation of noncanonical Phe analogs with distinct electrostatics at two positions. We found that amide-π interactions between AF9 and acyllysines are electrostatically tunable, with electron-rich rings providing more favorable interactions. This differs from trends in amide-heteroarene interactions and provides insightful information for therapeutic design. Additionally, we report for the first time that CH-π interactions at Phe28 directly contribute to AF9's recognition of acyllysines, illuminating differences among YEATS domains, as this residue is not highly conserved but has been shown to impart selectivity for specific post-translational modification.

Design of novel protein building modules and modular architectures

Nature uses only a limited number of protein topologies and while several folds have evolved independently over time, there are clearly many possible topologies that have not been explored by evolution. With recent advances of protein design concepts, computational modeling tools, high resolution and high-throughput experimental methods it is now possible to design new protein architectures. The collection of building blocks and design principles widened both in size and complexity, offering an expanded toolset for building new modular folds and functional protein structures. Here we review and discuss recent achievements of protein design, focusing in particular on the use and prospects of modular approaches for assembling new protein folds.

Secondary Forces in Protein Folding

A complete inventory of the forces governing protein folding is critical for productive protein modeling, including structure prediction and de novo design, as well as understanding protein misfolding diseases of clinical significance. The dominant contributors to protein folding include the hydrophobic effect and conventional hydrogen bonding, along with Coulombic and van der Waals interactions. Over the past few decades, important additional contributors have been identified, including C-H···O hydrogen bonding, n→π* interactions, C5 hydrogen bonding, chalcogen bonding, and interactions involving aromatic rings (cation-π, X-H···π, π-π, anion-π, and sulfur-arene). These secondary contributions fall into two general classes: (1) weak but abundant interactions of the protein main chain and (2) strong but less frequent interactions involving protein side chains. Though interactions with high individual energies play important roles in specifying nonlocal molecular contacts and ligand binding, we estimate that weak but abundant interactions are likely to make greater overall contributions to protein folding, particularly at the level of secondary structure. Further research is likely to illuminate additional roles of these noncanonical interactions and could also reveal contributions yet unknown.

Stabilizing and Understanding a Miniprotein by Rational Redesign

Miniproteins reduce the complexity of the protein-folding problem allowing systematic studies of contributions to protein folding and stabilization. Here, we describe the rational redesign of a miniprotein, PPα, comprising a polyproline II helix, a loop, and an α helix. The redesign provides a de novo framework for interrogating noncovalent interactions. Optimized PPα has significantly improved thermal stability with a midpoint unfolding temperature (T_M) of 51 °C. Its nuclear magnetic resonance structure indicates a density of stabilizing noncovalent interactions that is higher than that of the parent peptide, specifically an increased number of CH-π interactions. In part, we attribute this to improved long-range electrostatic interactions between the two helical elements. We probe further sequence-stability relationships in the miniprotein through a series of rational mutations.

The Singular NMR Fingerprint of a Polyproline II Helical Bundle

Polyproline II (PPII) helices play vital roles in biochemical recognition events and structures like collagen and form part of the conformational landscapes of intrinsically disordered proteins (IDPs). Nevertheless, this structure is generally hard to detect and quantify. Here, we report the first thorough NMR characterization of a PPII helical bundle protein, the Hypogastrura harveyi "snow flea" antifreeze protein (sfAFP). J-couplings and nuclear Overhauser enhancement spectroscopy confirm a natively folded structure consisting of six PPII helices. NMR spectral analyses reveal quite distinct Hα2 versus Hα3 chemical shifts for 28 Gly residues as well as ¹³Cα, ¹⁵N, and ¹HN conformational chemical shifts (Δδ) unique to PPII helical bundles. The ¹⁵N Δδ and ¹HN Δδ values and small negative ¹HN temperature coefficients evince hydrogen-bond formation. ¹H-¹⁵N relaxation measurements reveal that the backbone structure is generally highly rigid on ps-ns time scales. NMR relaxation parameters and biophysical characterization reveal that sfAFP is chiefly a dimer. For it, a structural model featuring the packing of long, flat hydrophobic faces at the dimer interface is advanced. The conformational stability, measured by amide H/D exchange to be 6.24 ± 0.2 kcal·mol^-1, is elevated. These are extraordinary findings considering the great entropic cost of fixing Gly residues and, together with the remarkable upfield chemical shifts of 28 Gly ¹Hα, evidence significant stabilizing contributions from CαHα ||| O═C hydrogen bonds. These stabilizing interactions are corroborated by density functional theory calculations and natural bonding orbital analysis. The singular conformational chemical shifts, J-couplings, high hNOE ratios, small negative temperature coefficients, and slowed H/D exchange constitute a unique set of fingerprints to identify PPII helical bundles, which may be formed by hundreds of Gly-rich motifs detected in sequence databases. These results should aid the quantification of PPII helices in IDPs, the development of improved antifreeze proteins, and the incorporation of PPII helices into novel designed proteins.

CXC-Mediated Cellular Uptake of Miniproteins: Forsaking "Arginine Magic"

Miniproteins have a size between that of larger biologics and small molecules and presumably possess the advantages of both; they represent an expanding class of promising scaffolds for the design of affinity reagents, enzymes, and therapeutics. Conventional strategies to promote cellular uptake of miniproteins rely on extensive grafting or embedding of arginine residues. However, the requirement of using cationic arginines would cause problems to the modified miniproteins, for example, low solubility in solutions (proneness of aggregation) and potential toxicity, which are open secrets in the peptide and protein communities. In this work, we report that the cell-permeability of cationic miniproteins can be further markedly increased through appending a magic CXC (cysteine- any-cysteine) motif, which takes advantage of thiol-disulfide exchanges on the cell surface. More importantly, we discovered that the high cell permeability of the CXC-appended miniproteins can still be preserved when the embedded arginines are all substituted with lysine residues, indicating that the "arginine magic" essential to almost all cell-permeable peptides and (mini)proteins is not required for the CXC-mediated cellular uptake. This finding provides a new avenue for designing highly cell-permeable miniproteins without compromise of potential toxicity and stability arising from arginine embedding or grafting.

The Intracellular Loop of the Na+/Ca2+ Exchanger Contains an "Awareness Ribbon"-Shaped Two-Helix Bundle Domain

The Na⁺/Ca²⁺ exchanger (NCX) is a ubiquitous single-chain membrane protein that plays a major role in regulating the intracellular Ca²⁺ homeostasis by the counter transport of Na⁺ and Ca²⁺ across the cell membrane. Other than its prokaryotic counterpart, which contains only the transmembrane domain and is self-sufficient as an active ion transporter, the eukaryotic NCX protein possesses in addition a large intracellular loop that senses intracellular calcium signals and controls the activation of ion transport across the membrane. This provides a necessary layer of regulation for the more complex function of eukaryotic cells. The Ca²⁺ sensor in the intracellular loop is known as the Ca²⁺-binding domain (CBD12). However, how the signaling of the allosteric intracellular Ca²⁺ binding propagates and results in transmembrane ion transportation still lacks a detailed explanation. Further structural and dynamics characterization of the intracellular loop flanking both sides of CBD12 is therefore imperative. Here, we report the identification and characterization of another structured domain that is N-terminal to CBD12 in the intracellular loop using solution nuclear magnetic resonance (NMR) spectroscopy. The atomistic structure of this domain reveals that two tandem long α-helices, connected by a short linker, form a stable crossover two-helix bundle (THB), resembling an "awareness ribbon". Considering the highly conserved amino acid sequence of the THB domain, the detailed structural and dynamics properties of the THB domain will be common among NCXs from different species and will contribute toward the understanding of the regulatory mechanism of eukaryotic Na⁺/Ca²⁺ exchangers.

Increasing Enzyme Stability and Activity through Hydrogen Bond-Enhanced Halogen Bonds

The construction of more stable proteins is important in biomolecular engineering, particularly in the design of biologics-based therapeutics. We show here that replacing the tyrosine at position 18 (Y18) of T4 lysozyme with the unnatural amino acid m-chlorotyrosine (^mClY) increases both the thermal stability (increasing the melting temperature by ∼1 °C and the melting enthalpy by 3 kcal/mol) and the enzymatic activity at elevated temperatures (15% higher than that of the parent enzyme at 40 °C) of this classic enzyme. The chlorine of ^mClY forms a halogen bond (XB) to the carbonyl oxygen of the peptide bond at glycine 28 (G28) in a tight loop near the active site. In this case, the XB potential of the typically weak XB donor Cl is shown from quantum chemical calculations to be significantly enhanced by polarization via an intramolecular hydrogen bond (HB) from the adjacent hydroxyl substituent of the tyrosyl side chain, resulting in a distinctive synergistic HB-enhanced XB (or HeX-B for short) interaction. The larger halogens (bromine and iodine) are not well accommodated within this same loop and, consequently, do not exhibit the effects on protein stability or function associated with the HeX-B interaction. Thus, we have for the first time demonstrated that an XB can be engineered to stabilize and increase the activity of an enzyme, with the increased stabilizing potential of the HeX-B further extending the application of halogenated amino acids in the design of more stable protein therapeutics.

Design of a Short Thermally Stable α-Helix Embedded in a Macrocycle

Although helices play key roles in peptide-protein and protein-protein interactions, the helical conformation is generally unstable for short peptides (10-15 residues) in aqueous solution in the absence of their binding partners. Thus, stabilizing the helical conformation of peptides can lead to increases in binding potency, specificity, and stability towards proteolytic degradation. Helices have been successfully stabilized by introducing side chain-to-side chain crosslinks within the central portion of the helix. However, this approach leaves the ends of the helix free, thus leading to fraying and exposure of the non-hydrogen-bonded amide groups to solvent. Here, we develop a "capped-strapped" peptide strategy to stabilize helices by embedding the entire length of the helix within a macrocycle, which also includes a semirigid organic template as well as end-capping interactions. We have designed a ten-residue capped-strapped helical peptide that behaves like a miniprotein, with a cooperative thermal unfolding transition and Tm ≈70 °C, unprecedented for helical peptides of this length. The NMR structure determination confirmed the design, and X-ray crystallography revealed a novel quaternary structure with implications for foldamer design.

A Comprehensive Analysis of Anion-Quadrupole Interactions in Protein Structures

The edgewise interactions of anions with phenylalanine (Phe) aromatic rings in proteins, known as anion-quadrupole interactions, have been well studied. However, the anion-quadrupole interactions of the tyrosine (Tyr) and tryptophan (Trp) rings have been less well studied, probably because these have been considered weaker than interactions of anions hydrogen bonded to Trp/Tyr side chains. Distinguishing such hydrogen bonding interactions, we comprehensively surveyed the edgewise interactions of certain anions (aspartate, glutamate, and phosphate) with Trp, Tyr, and Phe rings in high-resolution, nonredundant protein single chains and interfaces (protein-protein, DNA/RNA-protein, and membrane-protein). Trp/Tyr anion-quadrupole interactions are common, with Trp showing the highest propensity and average interaction energy for this type of interaction. The energy of an anion-quadrupole interaction (-15.0 to 0.0 kcal/mol, based on quantum mechanical calculations) depends not only on the interaction geometry but also on the ring atom. The phosphate anions at DNA/RNA-protein interfaces interact with aromatic residues with energies comparable to that of aspartate/glutamate anion-quadrupole interactions. At DNA-protein interfaces, the frequency of aromatic ring participation in anion-quadrupole interactions is comparable to that of positive charge participation in salt bridges, suggesting an underappreciated role for anion-quadrupole interactions at DNA-protein (or membrane-protein) interfaces. Although less frequent than salt bridges in single-chain proteins, we observed highly conserved anion-quadrupole interactions in the structures of remote homologues, and evolutionary covariance-based residue contact score predictions suggest that conserved anion-quadrupole interacting pairs, like salt bridges, contribute to polypeptide folding, stability, and recognition.

Direct Observation of CH/CH van der Waals Interactions in Proteins by NMR

van der Waals interactions are important to protein stability and function. These interactions are usually identified empirically based on protein 3D structures. In this work, we performed a solution nuclear magnetic resonance (NMR) spectroscopy study of van der Waals interactions by detecting the through-space vdw JCC-coupling between protein aliphatic side chain groups. Specifically, vdw JCC-coupling values up to ∼0.5 Hz were obtained between the methyl and nearby aliphatic groups in protein GB3, providing direct experimental evidence for the van der Waals interactions. Quantum mechanical calculations suggest that the J-coupling is correlated with the exchange-repulsion term of van der Waals interaction. NMR detection of vdw JCC-coupling offers a new tool to characterize such interactions in proteins.

The fourth wave of biocatalysis is approaching

Biocatalysis has undergone a tremendous development in the past few years. A plethora of methods enable the rather rapid tailored-design of an enzyme for a targeted reaction such as asymmetric synthesis of a chiral building block by the combination of information from sequence and structure databases with modern molecular biology methods and high-throughput screening tools. Moreover, novel non-natural reactions could be implemented into protein scaffolds and new enzyme classes are emerging, both broadening the repertoire of reactions now available for organic synthesis. Furthermore, impressive examples of metabolic engineering-the combination of several newly introduced reaction steps in a microbial host-have been developed, paving the way for large-scale processes for both pharmaceuticals and bulk chemicals. This contribution highlights recent developments in this area and points out future challenges.This article is part of a discussion meeting issue 'Providing sustainable catalytic solutions for a rapidly changing world'.

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.

How do miniproteins fold?

A high-throughput study yields libraries of miniproteins that help to explain how proteins are stabilized