A simple model of backbone flexibility improves modeling of side-chain conformational variability

The Open Force Field Initiative: Open Software and Open Science for Molecular Modeling

Force fields are a key component of physics-based molecular modeling, describing the energies and forces in a molecular system as a function of the positions of the atoms and molecules involved. Here, we provide a review and scientific status report on the work of the Open Force Field (OpenFF) Initiative, which focuses on the science, infrastructure and data required to build the next generation of biomolecular force fields. We introduce the OpenFF Initiative and the related OpenFF Consortium, describe its approach to force field development and software, and discuss accomplishments to date as well as future plans. OpenFF releases both software and data under open and permissive licensing agreements to enable rapid application, validation, extension, and modification of its force fields and software tools. We discuss lessons learned to date in this new approach to force field development. We also highlight ways that other force field researchers can get involved, as well as some recent successes of outside researchers taking advantage of OpenFF tools and data.

Uncovering Protein Ensembles: Automated Multiconformer Model Building for X-ray Crystallography and Cryo-EM

In their folded state, biomolecules exchange between multiple conformational states that are crucial for their function. Traditional structural biology methods, such as X-ray crystallography and cryogenic electron microscopy (cryo-EM), produce density maps that are ensemble averages, reflecting molecules in various conformations. Yet, most models derived from these maps explicitly represent only a single conformation, overlooking the complexity of biomolecular structures. To accurately reflect the diversity of biomolecular forms, there is a pressing need to shift towards modeling structural ensembles that mirror the experimental data. However, the challenge of distinguishing signal from noise complicates manual efforts to create these models. In response, we introduce the latest enhancements to qFit, an automated computational strategy designed to incorporate protein conformational heterogeneity into models built into density maps. These algorithmic improvements in qFit are substantiated by superior R f r e e and geometry metrics across a wide range of proteins. Importantly, unlike more complex multicopy ensemble models, the multiconformer models produced by qFit can be manually modified in most major model building software (e.g. Coot) and fit can be further improved by refinement using standard pipelines (e.g. Phenix, Refmac, Buster). By reducing the barrier of creating multiconformer models, qFit can foster the development of new hypotheses about the relationship between macromolecular conformational dynamics and function.

Analyses of Mutation Displacements from Homology Models

Evaluation of the structural perturbations introduced by a single amino acid mutation is the main issue for protein structural biology. We propose here to present some recent advances in methods, allowing the splitting of distortion between the actual substitution effect and the contribution of the local flexibility of the position where the mutation occurs. Its main drawback is the need of many structures with a single mutation in each of them. To bypass this difficulty, we propose to use molecular modeling tools, with several software enabling us to build a model from a template, given the sequence. As a proof of concept, we rely on a gold standard, the human lysozyme. Both wild-type and three mutant structures are available in the PDB. Two of these mutations result in amyloid fibril formation, and the last one is neutral. As a conclusion, irrespective of the algorithm used for modeling, side chain conformations at the site of mutation are reliable, although long-range effects are out of reach of these tools.

Computation-Aided Engineering of Cytochrome P450 for the Production of Pravastatin

CYP105AS1 is a cytochrome P450 from Amycolatopsis orientalis that catalyzes monooxygenation of compactin to 6-epi-pravastatin. For fermentative production of the cholesterol-lowering drug pravastatin, the stereoselectivity of the enzyme needs to be inverted, which has been partially achieved by error-prone PCR mutagenesis and screening. In the current study, we report further optimization of the stereoselectivity by a computationally aided approach. Using the CoupledMoves protocol of Rosetta, a virtual library of mutants was designed to bind compactin in a pro-pravastatin orientation. By examining the frequency of occurrence of beneficial substitutions and rational inspection of their interactions, a small set of eight mutants was predicted to show the desired selectivity and these variants were tested experimentally. The best CYP105AS1 variant gave >99% stereoselective hydroxylation of compactin to pravastatin, with complete elimination of the unwanted 6-epi-pravastatin diastereomer. The enzyme-substrate complexes were also examined by ultrashort molecular dynamics simulations of 50 × 100 ps and 5 × 22 ns, which revealed that the frequency of occurrence of near-attack conformations agreed with the experimentally observed stereoselectivity. These results show that a combination of computational methods and rational inspection could improve CYP105AS1 stereoselectivity beyond what was obtained by directed evolution. Moreover, the work lays out a general in silico framework for specificity engineering of enzymes of known structure.

Molecular insights into the regulation of constitutive activity by RNA editing of 5HT2C serotonin receptors

RNA editing is a process by which post-transcriptional changes of mRNA nucleotides alter protein function through modification of the amino acid content. The 5HT_2C serotonin receptor, which undergoes 32 distinct RNA-editing events leading to 24 protein isoforms, is a notable example of this process. These 5HT_2C isoforms display differences in constitutive activity, agonist/inverse agonist potencies, and efficacies. To elucidate the molecular mechanisms responsible for these effects of RNA editing, we present four active-state 5HT_2C-transducer-coupled structures of three representative isoforms (INI, VGV, and VSV) with the selective drug lorcaserin (Belviq) and the classic psychedelic psilocin. We also provide a comprehensive analysis of agonist activation and constitutive activity across all 24 protein isoforms. Collectively, these findings reveal a unique hydrogen-bonding network located on intracellular loop 2 that is subject to RNA editing, which differentially affects GPCR constitutive and agonist signaling activities.

It is established that the 5HT_2C receptor undergoes RNA editing leading to 24 isoforms. Several isoforms exhibit changes in basal activity and are linked to pathologies. Gumpper et al. have done a systematic structure-function characterization of all the isoforms revealing the underlying mechanisms that govern basal activity of the 5HT_2C receptor.

Key aspects of the past 30 years of protein design

Proteins are the workhorse of life. They are the building infrastructure of living systems; they are the most efficient molecular machines known, and their enzymatic activity is still unmatched in versatility by any artificial system. Perhaps proteins' most remarkable feature is their modularity. The large amount of information required to specify each protein's function is analogically encoded with an alphabet of just ∼20 letters. The protein folding problem is how to encode all such information in a sequence of 20 letters. In this review, we go through the last 30 years of research to summarize the state of the art and highlight some applications related to fundamental problems of protein evolution.

Modular peptide binders - development of a predictive technology as alternative for reagent antibodies

Current biomedical research and diagnostics critically depend on detection agents for specific recognition and quantification of protein molecules. Monoclonal antibodies have been used for this purpose over decades and facilitated numerous biological and biomedical investigations. Recently, however, it has become apparent that many commercial reagent antibodies lack specificity or do not recognize their target at all. Thus, synthetic alternatives are needed whose complex designs are facilitated by multidisciplinary approaches incorporating experimental protein engineering with computational modeling. Here, we review the status of such an engineering endeavor based on the modular armadillo repeat protein scaffold and discuss challenges in its implementation.

A Collection of Designed Peptides to Target SARS-CoV-2 Spike RBD-ACE2 Interaction

The angiotensin-converting enzyme 2 (ACE2) is the receptor used by SARS-CoV and SARS-CoV-2 coronaviruses to attach to cells via the receptor-binding domain (RBD) of their viral spike protein. Since the start of the COVID-19 pandemic, several structures of protein complexes involving ACE2 and RBD as well as monoclonal antibodies and nanobodies have become available. We have leveraged the structural data to design peptides to target the interaction between the RBD of SARS-CoV-2 and ACE2 and SARS-CoV and ACE2, as contrasting exemplar, as well as the dimerization surface of ACE2 monomers. The peptides were modelled using our original method: PiPreD that uses native elements of the interaction between the targeted protein and cognate partner(s) that are subsequently included in the designed peptides. These peptides recapitulate stretches of residues present in the native interface plus novel and highly diverse conformations surrogating key interactions at the interface. To facilitate the access to this information we have created a freely available and dedicated web-based repository, PepI-Covid19 database, providing convenient access to this wealth of information to the scientific community with the view of maximizing its potential impact in the development of novel therapeutic and diagnostic agents.

Recent advances in de novo protein design: Principles, methods, and applications

The computational de novo protein design is increasingly applied to address a number of key challenges in biomedicine and biological engineering. Successes in expanding applications are driven by advances in design principles and methods over several decades. Here, we review recent innovations in major aspects of the de novo protein design and include how these advances were informed by principles of protein architecture and interactions derived from the wealth of structures in the Protein Data Bank. We describe developments in de novo generation of designable backbone structures, optimization of sequences, design scoring functions, and the design of the function. The advances not only highlight design goals reachable now but also point to the challenges and opportunities for the future of the field.

Ensemble-based enzyme design can recapitulate the effects of laboratory directed evolution in silico

The creation of artificial enzymes is a key objective of computational protein design. Although de novo enzymes have been successfully designed, these exhibit low catalytic efficiencies, requiring directed evolution to improve activity. Here, we use room-temperature X-ray crystallography to study changes in the conformational ensemble during evolution of the designed Kemp eliminase HG3 (k_cat/K_M 146 M⁻¹s⁻¹). We observe that catalytic residues are increasingly rigidified, the active site becomes better pre-organized, and its entrance is widened. Based on these observations, we engineer HG4, an efficient biocatalyst (k_cat/K_M 103,000 M⁻¹s⁻¹) containing key first and second-shell mutations found during evolution. HG4 structures reveal that its active site is pre-organized and rigidified for efficient catalysis. Our results show how directed evolution circumvents challenges inherent to enzyme design by shifting conformational ensembles to favor catalytically-productive sub-states, and suggest improvements to the design methodology that incorporate ensemble modeling of crystallographic data.

Kemp eliminases are artificial enzymes that catalyze the concerted deprotonation and ring-opening of benzisoxazoles. Here, the authors use room-temperature X-ray crystallography to investigate changes to the conformational ensemble of the Kemp eliminase HG3 along a directed evolutionary trajectory, and develop an experimentally guided, ensemble-based computational enzyme design procedure.

Macromolecular modeling and design in Rosetta: recent methods and frameworks

The Rosetta software for macromolecular modeling, docking and design is extensively used in laboratories worldwide. During two decades of development by a community of laboratories at more than 60 institutions, Rosetta has been continuously refactored and extended. Its advantages are its performance and interoperability between broad modeling capabilities. Here we review tools developed in the last 5 years, including over 80 methods. We discuss improvements to the score function, user interfaces and usability. Rosetta is available at http://www.rosettacommons.org.

Affinity-matured variants derived from nimotuzumab keep the original fine specificity and exhibit superior biological activity

Nimotuzumab is a humanized monoclonal antibody against the Epidermal Growth Factor Receptor with a long history of therapeutic use, recognizing an epitope different from the ones targeted by other antibodies against the same antigen. It is also distinguished by much less toxicity resulting in a better safety profile, which has been attributed to its lower affinity compared to these other antibodies. Nevertheless, the ideal affinity window for optimizing the balance between anti-tumor activity and toxic effects has not been determined. In the current work, the paratope of the phage-displayed nimotuzumab Fab fragment was evolved in vitro to obtain affinity-matured variants. Soft-randomization of heavy chain variable region CDRs and phage selection resulted in mutated variants with improved binding ability. Two recombinant antibodies were constructed using these variable regions, which kept the original fine epitope specificity and showed moderate affinity increases against the target (3-4-fold). Such differences were translated into a greatly enhanced inhibitory capacity upon ligand-induced receptor phosphorylation on tumor cells. The new antibodies, named K4 and K5, are valuable tools to explore the role of affinity in nimotuzumab biological properties, and could be used for applications requiring a fine-tuning of the balance between binding to tumor cells and healthy tissues.

Comparison of Rosetta flexible-backbone computational protein design methods on binding interactions

Computational design of binding sites in proteins remains difficult, in part due to limitations in our current ability to sample backbone conformations that enable precise and accurate geometric positioning of side chains during sequence design. Here we present a benchmark framework for comparison between flexible-backbone design methods applied to binding interactions. We quantify the ability of different flexible backbone design methods in the widely used protein design software Rosetta to recapitulate observed protein sequence profiles assumed to represent functional protein/protein and protein/small molecule binding interactions. The CoupledMoves method, which combines backbone flexibility and sequence exploration into a single acceptance step during the sampling trajectory, better recapitulates observed sequence profiles than the BackrubEnsemble and FastDesign methods, which separate backbone flexibility and sequence design into separate acceptance steps during the sampling trajectory. Flexible-backbone design with the CoupledMoves method is a powerful strategy for reducing sequence space to generate targeted libraries for experimental screening and selection.

Molecular basis for heat desensitization of TRPV1 ion channels

The transient receptor potential vanilloid 1 (TRPV1) ion channel is a prototypical molecular sensor for noxious heat in mammals. Its role in sustained heat response remains poorly understood, because rapid heat-induced desensitization (Dh) follows tightly heat-induced activation (Ah). To understand the physiological role and structural basis of Dh, we carried out a comparative study of TRPV1 channels in mouse (mV1) and those in platypus (pV1), which naturally lacks Dh. Here we show that a temperature-sensitive interaction between the N- and C-terminal domains of mV1 but not pV1 drives a conformational rearrangement in the pore leading to Dh. We further show that knock-in mice expressing pV1 sensed heat normally but suffered scald damages in a hot environment. Our findings suggest that Dh evolved late during evolution as a protective mechanism and a delicate balance between Ah and Dh is crucial for mammals to sense and respond to noxious heat.

The conformational wave in capsaicin activation of transient receptor potential vanilloid 1 ion channel

The capsaicin receptor TRPV1 has been intensively studied by cryo-electron microscopy and functional tests. However, though the apo and capsaicin-bound structural models are available, the dynamic process of capsaicin activation remains intangible, largely due to the lack of a capsaicin-induced open structural model and the low occupancy of the transition states. Here we report that reducing temperature toward the freezing point substantially increased channel closure events even in the presence of saturating capsaicin. We further used a combination of fluorescent unnatural amino acid (fUAA) incorporation, computational modeling, and rate-equilibrium linear free-energy relationships analysis (Φ-analysis) to derive the fully open capsaicin-bound state model, and reveal how the channel transits from the apo to the open state. We observed that capsaicin initiates a conformational wave that propagates through the S4–S5 linker towards the S6 bundle and finally reaching the selectivity filter. Our study provides a temporal mechanism for capsaicin activation of TRPV1.

The capsaicin receptor TRPV1 has been structurally characterized, but the capsaicin activation dynamics remain elusive. Here authors use fluorescent unnatural amino acid incorporation, computational modeling and Φ-analysis to derive the capsaicin-bound open state model and reveal the capsaicin induced conformational changes.

Chimeric Fatty Acyl-Acyl Carrier Protein Thioesterases Provide Mechanistic Insight into Enzyme Specificity and Expression

Medium-chain fatty acids are commodity chemicals. Increasing and modifying the activity of thioesterases (TEs) on medium-chain fatty acyl-acyl carrier protein (acyl-ACP) esters may enable a high-yield microbial production of these molecules. The plant Cuphea palustris harbors two distinct TEs: C. palustris FatB1 (CpFatB1) (C₈ specificity, lower activity) and CpFatB2 (C₁₄ specificity, higher activity) with 78% sequence identity. We combined structural features from these two enzymes to create several chimeric TEs, some of which showed nonnatural fatty acid production as measured by an enzymatic assay and gas chromatography-mass spectrometry (GC-MS). Notably, chimera 4 exhibited an increased C₈ fatty acid production in correlation with improved microbial expression. This chimera led us to identify CpFatB2-specific amino acids between positions 219 and 272 that lead to higher protein levels. Chimera 7 produced a broad range of fatty acids and appeared to combine a fatty acid binding pocket with long-chain specificity and an ACP interaction site that may activate fatty acid extrusion. Using homology modeling and in silico docking with ACP, we identified a "positive patch" within amino acids 162 to 218, which may direct the ACP interaction and regulate access to short-chain fatty acids. On the basis of this modeling, we transplanted putative ACP interaction sequences from CpFatB1 into CpFatB2 and created a chimeric thioesterase that produced medium-chain as well as long-chain fatty acids. Thus, the engineering of chimeric enzymes and characterizing their microbial activity and chain-length specificity suggested mechanistic insights into TE functions and also generated thioesterases with potentially useful properties. These observations may inform a rational engineering of TEs to allow alkyl chain length control.IMPORTANCE Medium-chain fatty acids are important commodity chemicals. These molecules are used as plastic precursors and in shampoos and other detergents and could be used as biofuel precursors if production economics were favorable. Hydrocarbon-based liquid fuels must be optimized to have a desired boiling point, low freezing point, low viscosity, and other physical characteristics. Similarly, the solubility and harshness of detergents and the flexibility of plastic polymers can be modulated. The length and distribution of the carbon chains in the hydrophobic tails determine these properties. The biological synthesis of cell membranes and fatty acids produces chains of primarily 16 to 18 carbons, which give rise to current biofuels. The ultimate goal of the work presented here is to engineer metabolic pathways to produce designer molecules with the correct number of carbons in a chain, so that such molecules could be used directly as specialty commodity chemicals or as fuels after minimal processing.

Flex ddG: Rosetta Ensemble-Based Estimation of Changes in Protein-Protein Binding Affinity upon Mutation

Computationally modeling changes in binding free energies upon mutation (interface ΔΔ G) allows large-scale prediction and perturbation of protein-protein interactions. Additionally, methods that consider and sample relevant conformational plasticity should be able to achieve higher prediction accuracy over methods that do not. To test this hypothesis, we developed a method within the Rosetta macromolecular modeling suite (flex ddG) that samples conformational diversity using "backrub" to generate an ensemble of models and then applies torsion minimization, side chain repacking, and averaging across this ensemble to estimate interface ΔΔ G values. We tested our method on a curated benchmark set of 1240 mutants, and found the method outperformed existing methods that sampled conformational space to a lesser degree. We observed considerable improvements with flex ddG over existing methods on the subset of small side chain to large side chain mutations, as well as for multiple simultaneous non-alanine mutations, stabilizing mutations, and mutations in antibody-antigen interfaces. Finally, we applied a generalized additive model (GAM) approach to the Rosetta energy function; the resulting nonlinear reweighting model improved the agreement with experimentally determined interface ΔΔ G values but also highlighted the necessity of future energy function improvements.

Data driven flexible backbone protein design

Protein design remains an important problem in computational structural biology. Current computational protein design methods largely use physics-based methods, which make use of information from a single protein structure. This is despite the fact that multiple structures of many protein folds are now readily available in the PDB. While ensemble protein design methods can use multiple protein structures, they treat each structure independently. Here, we introduce a flexible backbone strategy, FlexiBaL-GP, which learns global protein backbone movements directly from multiple protein structures. FlexiBaL-GP uses the machine learning method of Gaussian Process Latent Variable Models to learn a lower dimensional representation of the protein coordinates that best represent backbone movements. These learned backbone movements are used to explore alternative protein backbones, while engineering a protein within a parallel tempered MCMC framework. Using the human ubiquitin-USP21 complex as a model we demonstrate that our design strategy outperforms current strategies for the interface design task of identifying tight binding ubiquitin variants for USP21.

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

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

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

Flexible Backbone Methods for Predicting and Designing Peptide Specificity

Protein-protein interactions play critical roles in essentially every cellular process. These interactions are often mediated by protein interaction domains that enable proteins to recognize their interaction partners, often by binding to short peptide motifs. For example, PDZ domains, which are among the most common protein interaction domains in the human proteome, recognize specific linear peptide sequences that are often at the C-terminus of other proteins. Determining the set of peptide sequences that a protein interaction domain binds, or it's "peptide specificity," is crucial for understanding its cellular function, and predicting how mutations impact peptide specificity is important for elucidating the mechanisms underlying human diseases. Moreover, engineering novel cellular functions for synthetic biology applications, such as the biosynthesis of biofuels or drugs, requires the design of protein interaction specificity to avoid crosstalk with native metabolic and signaling pathways. The ability to accurately predict and design protein-peptide interaction specificity is therefore critical for understanding and engineering biological function. One approach that has recently been employed toward accomplishing this goal is computational protein design. This chapter provides an overview of recent methodological advances in computational protein design and highlights examples of how these advances can enable increased accuracy in predicting and designing peptide specificity.

Finding Our Way in the Dark Proteome

The traditional structure-function paradigm has provided significant insights for well-folded proteins in which structures can be easily and rapidly revealed by X-ray crystallography beamlines. However, approximately one-third of the human proteome is comprised of intrinsically disordered proteins and regions (IDPs/IDRs) that do not adopt a dominant well-folded structure, and therefore remain "unseen" by traditional structural biology methods. This Perspective considers the challenges raised by the "Dark Proteome", in which determining the diverse conformational substates of IDPs in their free states, in encounter complexes of bound states, and in complexes retaining significant disorder requires an unprecedented level of integration of multiple and complementary solution-based experiments that are analyzed with state-of-the art molecular simulation, Bayesian probabilistic models, and high-throughput computation. We envision how these diverse experimental and computational tools can work together through formation of a "computational beamline" that will allow key functional features to be identified in IDP structural ensembles.

Protein design: Past, present, and future

Building on the pioneering work of Ho and DeGrado (J Am Chem Soc 1987, 109, 6751-6758) in the late 1980s, protein design approaches have revealed many fundamental features of protein structure and stability. We are now in the era that the early work presaged - the design of new proteins with practical applications and uses. Here we briefly survey some past milestones in protein design, in addition to highlighting recent progress and future aspirations.

Novel Ricin Subunit Antigens With Enhanced Capacity to Elicit Toxin-Neutralizing Antibody Responses in Mice

RiVax is a candidate ricin toxin subunit vaccine antigen that has proven to be safe in human phase I clinical trials. In this study, we introduced double and triple cavity-filling point mutations into the RiVax antigen with the expectation that stability-enhancing modifications would have a beneficial effect on overall immunogenicity of the recombinant proteins. We demonstrate that 2 RiVax triple mutant derivatives, RB (V81L/C171L/V204I) and RC (V81I/C171L/V204I), when adsorbed to aluminum salts adjuvant and tested in a mouse prime-boost-boost regimen were 5- to 10-fold more effective than RiVax at eliciting toxin-neutralizing serum IgG antibody titers. Increased toxin neutralizing antibody values and seroconversion rates were evident at different antigen dosages and within 7 days after the first booster. Quantitative stability/flexibility relationships analysis revealed that the RB and RC mutations affect rigidification of regions spanning residues 98-103, which constitutes a known immunodominant neutralizing B-cell epitope. A more detailed understanding of the immunogenic nature of RB and RC may provide insight into the fundamental relationship between local protein stability and antibody reactivity.

Principles and Overview of Sampling Methods for Modeling Macromolecular Structure and Dynamics

Investigation of macromolecular structure and dynamics is fundamental to understanding how macromolecules carry out their functions in the cell. Significant advances have been made toward this end in silico, with a growing number of computational methods proposed yearly to study and simulate various aspects of macromolecular structure and dynamics. This review aims to provide an overview of recent advances, focusing primarily on methods proposed for exploring the structure space of macromolecules in isolation and in assemblies for the purpose of characterizing equilibrium structure and dynamics. In addition to surveying recent applications that showcase current capabilities of computational methods, this review highlights state-of-the-art algorithmic techniques proposed to overcome challenges posed in silico by the disparate spatial and time scales accessed by dynamic macromolecules. This review is not meant to be exhaustive, as such an endeavor is impossible, but rather aims to balance breadth and depth of strategies for modeling macromolecular structure and dynamics for a broad audience of novices and experts.

Protein backbone ensemble generation explores the local structural space of unseen natural homologs

MOTIVATION

Mutations in homologous proteins affect changes in the backbone conformation that involve a complex interplay of forces which are difficult to predict. Protein design algorithms need to anticipate these backbone changes in order to accurately calculate the energy of the structure given an amino acid sequence, without knowledge of the final, designed sequence. This is related to the problem of predicting small changes in the backbone between highly similar sequences.

RESULTS

We explored the ability of the Rosetta suite of protein design tools to move the backbone from its position in one structure (template) to its position in a close homologous structure (target) as a function of the diversity of a backbone ensemble constructed using the template structure, the percent sequence identity between the template and target, and the size of local zone being considered in the ensemble. We describe a pareto front in the likelihood of moving the backbone toward the target as a function of ensemble diversity and zone size. The equations and protocols presented here will be useful for protein design.

AVAILABILITY AND IMPLEMENTATION

PyRosetta scripts available at www.bioinfo.rpi.edu/bystrc/downloads.html#ensemble

CONTACT

bystrc@rpi.edu.

Identification of protein scaffolds for enzyme design using scaffold selection

The identification of suitable protein structures that can serve as scaffolds for the introduction of catalytic residues is crucial for the design of new enzymes. Here we describe how the automated and rapid scaffold search program ScaffoldSelection can be used to find the best starting points, namely protein structures that are most likely to tolerate the introduction and promote the proper formation of a specific catalytic motif.

Correlation between local structural dynamics of proteins inferred from NMR ensembles and evolutionary dynamics of homologues of known structure

Conformational changes in proteins are extremely important for their biochemical functions. Correlation between inherent conformational variations in a protein and conformational differences in its homologues of known structure is still unclear. In this study, we have used a structural alphabet called Protein Blocks (PBs). PBs are used to perform abstraction of protein 3-D structures into a 1-D strings of 16 alphabets (a-p) based on dihedral angles of overlapping pentapeptides. We have analyzed the variations in local conformations in terms of PBs represented in the ensembles of 801 protein structures determined using NMR spectroscopy. In the analysis of concatenated data over all the residues in all the NMR ensembles, we observe that the overall nature of inherent local structural variations in NMR ensembles is similar to the nature of local structural differences in homologous proteins with a high correlation coefficient of .94. High correlation at the alignment positions corresponding to helical and β-sheet regions is only expected. However, the correlation coefficient by considering only the loop regions is also quite high (.91). Surprisingly, segregated position-wise analysis shows that this high correlation does not hold true to loop regions at the structurally equivalent positions in NMR ensembles and their homologues of known structure. This suggests that the general nature of local structural changes is unique; however most of the local structural variations in loop regions of NMR ensembles do not correlate to their local structural differences at structurally equivalent positions in homologues.

Effect of single-point mutations on the stability and immunogenicity of a recombinant ricin A chain subunit vaccine antigen

There is great interest in the design and development of highly thermostable and immunogenic protein subunit vaccines for biodefense. In this study, we used two orthogonal and complementary computational protein design approaches to generate a series of single-point mutants of RiVax, an attenuated recombinant ricin A chain (RTA) protein subunit vaccine antigen. As assessed by differential scanning calorimetry, the conformational stabilities of the designed mutants ranged from 4°C less stable to 4.5°C more stable than RiVax, depending on solution pH. Two more thermostable (V18P, C171L) and two less thermostable (T13V, S89T) mutants that displayed native-like secondary and tertiary structures (as determined by circular dichroism and fluorescence spectral analysis, respectively) were tested for their capacity to elicit RTA-specific antibodies and toxin-neutralizing activity. Following a prime-boost regimen, we found qualitative differences with respect to specific antibody titers and toxin neutralizing antibody levels induced by the different mutants. Upon a second boost with the more thermostable mutant C171L, a statistically significant increase in RTA-specific antibody titers was observed when compared with RiVax-immunized mice. Notably, the results indicate that single residue changes can be made to the RiVax antigen that increase its thermal stability without adversely impacting the efficacy of the vaccine.

Flexible backbone sampling methods to model and design protein alternative conformations

Sampling alternative conformations is key to understanding how proteins work and engineering them for new functions. However, accurately characterizing and modeling protein conformational ensembles remain experimentally and computationally challenging. These challenges must be met before protein conformational heterogeneity can be exploited in protein engineering and design. Here, as a stepping stone, we describe methods to detect alternative conformations in proteins and strategies to model these near-native conformational changes based on backrub-type Monte Carlo moves in Rosetta. We illustrate how Rosetta simulations that apply backrub moves improve modeling of point mutant side-chain conformations, native side-chain conformational heterogeneity, functional conformational changes, tolerated sequence space, protein interaction specificity, and amino acid covariation across protein-protein interfaces. We include relevant Rosetta command lines and RosettaScripts to encourage the application of these types of simulations to other systems. Our work highlights that critical scoring and sampling improvements will be necessary to approximate conformational landscapes. Challenges for the future development of these methods include modeling conformational changes that propagate away from designed mutation sites and modulating backbone flexibility to predictively design functionally important conformational heterogeneity.

GRID: a high-resolution protein structure refinement algorithm

The energy-based refinement of protein structures generated by fold prediction algorithms to atomic-level accuracy remains a major challenge in structural biology. Energy-based refinement is mainly dependent on two components: (1) sufficiently accurate force fields, and (2) efficient conformational space search algorithms. Focusing on the latter, we developed a high-resolution refinement algorithm called GRID. It takes a three-dimensional protein structure as input and, using an all-atom force field, attempts to improve the energy of the structure by systematically perturbing backbone dihedrals and side-chain rotamer conformations. We compare GRID to Backrub, a stochastic algorithm that has been shown to predict a significant fraction of the conformational changes that occur with point mutations. We applied GRID and Backrub to 10 high-resolution (≤ 2.8 Å) crystal structures from the Protein Data Bank and measured the energy improvements obtained and the computation times required to achieve them. GRID resulted in energy improvements that were significantly better than those attained by Backrub while expending about the same amount of computational resources. GRID resulted in relaxed structures that had slightly higher backbone RMSDs compared to Backrub relative to the starting crystal structures. The average RMSD was 0.25 ± 0.02 Å for GRID versus 0.14 ± 0.04 Å for Backrub. These relatively minor deviations indicate that both algorithms generate structures that retain their original topologies, as expected given the nature of the algorithms.

Prediction of mutational tolerance in HIV-1 protease and reverse transcriptase using flexible backbone protein design

Predicting which mutations proteins tolerate while maintaining their structure and function has important applications for modeling fundamental properties of proteins and their evolution; it also drives progress in protein design. Here we develop a computational model to predict the tolerated sequence space of HIV-1 protease reachable by single mutations. We assess the model by comparison to the observed variability in more than 50,000 HIV-1 protease sequences, one of the most comprehensive datasets on tolerated sequence space. We then extend the model to a second protein, reverse transcriptase. The model integrates multiple structural and functional constraints acting on a protein and uses ensembles of protein conformations. We find the model correctly captures a considerable fraction of protease and reverse-transcriptase mutational tolerance and shows comparable accuracy using either experimentally determined or computationally generated structural ensembles. Predictions of tolerated sequence space afforded by the model provide insights into stability-function tradeoffs in the emergence of resistance mutations and into strengths and limitations of the computational model.

Many related protein sequences can be consistent with the structure and function of a given protein, suggesting that proteins may be quite robust to mutations. This tolerance to mutations is frequently exploited by pathogens. In particular, pathogens can rapidly evolve mutated proteins that have a new function - resistance against a therapeutic inhibitor - without abandoning other functions essential for the pathogen. This principle may also hold more generally: Proteins tolerant to mutational changes can more easily acquire new functions while maintaining their existing properties. The ability to predict the tolerance of proteins to mutation could thus help both to analyze the emergence of resistance mutations in pathogens and to engineer proteins with new functions. Here we develop a computational model to predict protein mutational tolerance towards point mutations accessible by single nucleotide changes, and validate it using two important pathogenic proteins and therapeutic targets: the protease and reverse transcriptase from HIV-1. The model provides insights into how resistance emerges and makes testable predictions on mutations that have not been seen yet. Similar models of mutational tolerance should be useful for characterizing and reengineering the functions of other proteins for which a three-dimensional structure is available.

Multistate approaches in computational protein design

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

The role of local backrub motions in evolved and designed mutations

Amino acid substitutions in protein structures often require subtle backbone adjustments that are difficult to model in atomic detail. An improved ability to predict realistic backbone changes in response to engineered mutations would be of great utility for the blossoming field of rational protein design. One model that has recently grown in acceptance is the backrub motion, a low-energy dipeptide rotation with single-peptide counter-rotations, that is coupled to dynamic two-state sidechain rotamer jumps, as evidenced by alternate conformations in very high-resolution crystal structures. It has been speculated that backrubs may facilitate sequence changes equally well as rotamer changes. However, backrub-induced shifts and experimental uncertainty are of similar magnitude for backbone atoms in even high-resolution structures, so comparison of wildtype-vs.-mutant crystal structure pairs is not sufficient to directly link backrubs to mutations. In this study, we use two alternative approaches that bypass this limitation. First, we use a quality-filtered structure database to aggregate many examples for precisely defined motifs with single amino acid differences, and find that the effectively amplified backbone differences closely resemble backrubs. Second, we directly apply a provably-accurate, backrub-enabled protein design algorithm to idealized versions of these motifs, and discover that the lowest-energy computed models match the average-coordinate experimental structures. These results support the hypothesis that backrubs participate in natural protein evolution and validate their continued use for design of synthetic proteins.

Protein design has the potential to generate useful molecules for medicine and chemistry, including sensors, drugs, and catalysts for arbitrary reactions. When protein design is carried out starting from an experimentally determined structure, as is often the case, one important aspect to consider is backbone flexibility, because in response to a mutation the backbone often must shift slightly to reconcile the new sidechain with its environment. In principle, one may model the backbone in many ways, but not all are physically realistic or experimentally validated. Here we study the "backrub" motion, which has been previously documented in atomic detail, but only for sidechain movements within single structures. By a twopronged approach involving both structural bioinformatics and computation with a principled design algorithm, we demonstrate that backrubs are sufficient to explain the backbone differences between mutation-related sets of very precisely defined motifs from the protein structure database. Our findings illustrate that backrubs are useful for describing evolutionary sequence change and, by extension, suggest that they are also appropriate for rational protein design calculations.

Increasing sequence diversity with flexible backbone protein design: the complete redesign of a protein hydrophobic core

Protein design tests our understanding of protein stability and structure. Successful design methods should allow the exploration of sequence space not found in nature. However, when redesigning naturally occurring protein structures, most fixed backbone design algorithms return amino acid sequences that share strong sequence identity with wild-type sequences, especially in the protein core. This behavior places a restriction on functional space that can be explored and is not consistent with observations from nature, where sequences of low identity have similar structures. Here, we allow backbone flexibility during design to mutate every position in the core (38 residues) of a four-helix bundle protein. Only small perturbations to the backbone, 1-2 Å, were needed to entirely mutate the core. The redesigned protein, DRNN, is exceptionally stable (melting point >140°C). An NMR and X-ray crystal structure show that the side chains and backbone were accurately modeled (all-atom RMSD = 1.3 Å).

A smoothed backbone-dependent rotamer library for proteins derived from adaptive kernel density estimates and regressions

Rotamer libraries are used in protein structure determination, prediction, and design. The backbone-dependent rotamer library consists of rotamer frequencies, mean dihedral angles, and variances as a function of the backbone dihedral angles. Structure prediction and design methods that employ backbone flexibility would strongly benefit from smoothly varying probabilities and angles. A new version of the backbone-dependent rotamer library has been developed using adaptive kernel density estimates for the rotamer frequencies and adaptive kernel regression for the mean dihedral angles and variances. This formulation allows for evaluation of the rotamer probabilities, mean angles, and variances as a smooth and continuous function of phi and psi. Continuous probability density estimates for the nonrotameric degrees of freedom of amides, carboxylates, and aromatic side chains have been modeled as a function of the backbone dihedrals and rotamers of the remaining degrees of freedom. New backbone-dependent rotamer libraries at varying levels of smoothing are available from http://dunbrack.fccc.edu.

Predicting the tolerated sequences for proteins and protein interfaces using RosettaBackrub flexible backbone design

Predicting the set of sequences that are tolerated by a protein or protein interface, while maintaining a desired function, is useful for characterizing protein interaction specificity and for computationally designing sequence libraries to engineer proteins with new functions. Here we provide a general method, a detailed set of protocols, and several benchmarks and analyses for estimating tolerated sequences using flexible backbone protein design implemented in the Rosetta molecular modeling software suite. The input to the method is at least one experimentally determined three-dimensional protein structure or high-quality model. The starting structure(s) are expanded or refined into a conformational ensemble using Monte Carlo simulations consisting of backrub backbone and side chain moves in Rosetta. The method then uses a combination of simulated annealing and genetic algorithm optimization methods to enrich for low-energy sequences for the individual members of the ensemble. To emphasize certain functional requirements (e.g. forming a binding interface), interactions between and within parts of the structure (e.g. domains) can be reweighted in the scoring function. Results from each backbone structure are merged together to create a single estimate for the tolerated sequence space. We provide an extensive description of the protocol and its parameters, all source code, example analysis scripts and three tests applying this method to finding sequences predicted to stabilize proteins or protein interfaces. The generality of this method makes many other applications possible, for example stabilizing interactions with small molecules, DNA, or RNA. Through the use of within-domain reweighting and/or multistate design, it may also be possible to use this method to find sequences that stabilize particular protein conformations or binding interactions over others.

Assessment of flexible backbone protein design methods for sequence library prediction in the therapeutic antibody Herceptin-HER2 interface

Computational protein design methods can complement experimental screening and selection techniques by predicting libraries of low-energy sequences compatible with a desired structure and function. Incorporating backbone flexibility in computational design allows conformational adjustments that should broaden the range of predicted low-energy sequences. Here, we evaluate computational predictions of sequence libraries from different protocols for modeling backbone flexibility using the complex between the therapeutic antibody Herceptin and its target human epidermal growth factor receptor 2 (HER2) as a model system. Within the program RosettaDesign, three methods are compared: The first two use ensembles of structures generated by Monte Carlo protocols for near-native conformational sampling: kinematic closure (KIC) and backrub, and the third method uses snapshots from molecular dynamics (MD) simulations. KIC or backrub methods were better able to identify the amino acid residues experimentally observed by phage display in the Herceptin-HER2 interface than MD snapshots, which generated much larger conformational and sequence diversity. KIC and backrub, as well as fixed backbone simulations, captured the key mutation Asp98Trp in Herceptin, which leads to a further threefold affinity improvement of the already subnanomolar parental Herceptin-HER2 interface. Modeling subtle backbone conformational changes may assist in the design of sequence libraries for improving the affinity of antibody-antigen interfaces and could be suitable for other protein complexes for which structural information is available.

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

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

Rapid flexible docking using a stochastic rotamer library of ligands

Existing flexible docking approaches model the ligand and receptor flexibility either separately or in a loosely coupled manner, which captures the conformational changes inefficiently. Here, we propose a flexible docking approach, MedusaDock, which models both ligand and receptor flexibility simultaneously with sets of discrete rotamers. We developed an algorithm to build the ligand rotamer library "on-the-fly" during docking simulations. MedusaDock benchmarks demonstrate a rapid sampling efficiency and high prediction accuracy in both self- (to the cocrystallized state) and cross-docking (to a state cocrystallized with a different ligand), the latter of which mimics the virtual screening procedure in computational drug discovery. We also perform a virtual screening test of four flexible kinase targets, including cyclin-dependent kinase 2, vascular endothelial growth factor receptor 2, HIV reverse transcriptase, and HIV protease. We find significant improvements of virtual screening enrichments when compared to rigid-receptor methods. The predictive power of MedusaDock in cross-docking and preliminary virtual-screening benchmarks highlights the importance to model both ligand and receptor flexibility simultaneously in computational docking.

Assessing the native state conformational distribution of ubiquitin by peptide acidity

At equilibrium, every energetically feasible conformation of a protein occurs with a non-zero probability. Quantitative analysis of protein flexibility is thus synonymous with determining the proper Boltzmann-weighting of this conformational distribution. The exchange reactivity of solvent-exposed amide hydrogens greatly varies with conformation, while the short-lived peptide anion intermediate implies an insensitivity to the dynamics of conformational motion. Amides that are well-exposed in model conformational ensembles of ubiquitin vary a million-fold in exchange rates which continuum dielectric methods can predict with an rmsd of 3. However, the exchange rates for many of the more rarely exposed amides are markedly overestimated in the PDB-deposited 2K39 and 2KN5 ubiquitin ensembles, while the 2NR2 ensemble predictions are largely consistent with those of the Boltzmann-weighted conformational distribution sampled at the level of 1%. The correlation between the fraction of solvent-accessible conformations for a given amide hydrogen and the exchange rate constant for that residue provides a useful monitor of the degree of completeness with which a given ensemble has sampled the energetically accessible conformational space. These exchange predictions correlate with the degree to which each ensemble deviates from a set of 46 ubiquitin X-ray structures. Kolmogorov-Smirnov analysis for the distribution of intra- and inter-ensemble pairwise structural rmsd values assisted the identification of a subensemble of 2K39 that eliminates the overestimations of hydrogen exchange rates observed for the full ensemble. The relative merits of incorporating experimental restraints into the conformational sampling process are compared to using these restraints as filters to select subpopulations consistent with the experimental data.