Abstract:Abstract The protein design problem is to identify an amino acid sequence which folds to a desired structure. Given Anfinsen’s thermodynamic hypothesis of folding, this can be recast as finding an amino acid sequence for which the lowest energy conformation is that structure. As this calculation involves not only all possible amino acid sequences but also all possible structures, most current approaches focus instead on the more tractable problem of finding the lowest energy amino acid sequence for the desired structure, often checking by protein structure prediction in a second step that the desired structure is indeed the lowest energy conformation for the designed sequence, and discarding the in many cases large fraction of designed sequences for which this is not the case. Here we show that by backpropagating gradients through the trRosetta structure prediction network from the desired structure to the input amino acid sequence, we can directly optimize over all possible amino acid sequences and all possible structures, and in one calculation explicitly design amino acid sequences predicted to fold into the desired structure and not any other. We find that trRosetta calculations, which consider the full conformational landscape, can be more effective than Rosetta single point energy estimations in predicting folding and stability of de novo designed proteins. We compare sequence design by landscape optimization to the standard fixed backbone sequence design methodology in Rosetta, and show that the results of the former, but not the latter, are sensitive to the presence of competing low-lying states. We show further that more funneled energy landscapes can be designed by combining the strengths of the two approaches: the low resolution trRosetta model serves to disfavor alternative states, and the high resolution Rosetta model, to create a deep energy minimum at the design target structure. Significance Computational protein design has primarily focused on finding sequences which have very low energy in the target designed structure. However, what is most relevant during folding is not the absolute energy of the folded state, but the energy difference between the folded state and the lowest lying alternative states. We describe a deep learning approach which captures the entire folding landscape, and show that it can enhance current protein design methods.

How to Solve Quantum Nonlinear Abelian Gauge Theory in Two Dimension in the Heisenberg Picture

A Discontinuous Potential Model for Protein-Protein Interactions.

Folding Lattice HP Model of Proteins Using the Bond-Fluctuation Model

Learning effective amino acid interactions through iterative stochastic techniques

A Coarse-Grained Approach to Protein Design: Learning from Design to Understand Folding

Transferable coarse-grained potential for $\textit{de novo}$ protein folding and design

The cavity method to protein design problem

Protein sequence design by explicit energy landscape optimization

Protein Design Based on the Relative Entropy

Protein Design by Integrating Machine Learning with Quantum Annealing and Quantum-inspired Optimization

Variational Theory for Site Resolved Protein Folding Free Energy Surfaces

A New-Protein Folding Recognition Potential Function

Monte Carlo Sampling of Protein Folding by Combining an All-Atom Physics-Based Model with a Native State Bias

Inverse Protein Folding Using Deep Bayesian Optimization

Lattice protein design using Bayesian learning

Protein structures and optimal folding emerging from a geometrical variational principle

Variational Methods for Biomolecular Modeling

Quasi-Physical Algorithm of an Off-Lattice Model for Protein Folding Problem

Energy Optimization for Off-Lattice Protein Folding

Protein Design Using Physics Informed Neural Networks

Variational Scheme to Compute Protein Reaction Pathways using Atomistic Force Fields with Explicit Solvent