Abstract:Proteins play a pivotal role in many biological processes, and changes in their amino acid sequences can lead to dysfunction and disease. These changes can affect protein folding or interaction with other biomolecules, such as preventing antibodies from inhibiting a viral infection or causing proteins to misfold. The ability to predict the effects of mutations in proteins is crucial. Although experimental techniques can accurately quantify the effect of mutations on protein folding free energies and protein- protein binding free energies, they are often time-consuming and costly. By contrast, computational techniques offer fast and cost-effective alternatives for estimating free energies, but they typically suffer from lower accuracy. Enhancing the accuracy of computational predictions is therefore of high importance, with the potential to greatly impact fields ranging from drug design to understanding disease mechanisms. One such widely used computational method, FoldX, is capable of rapidly predicting the relative folding stability (ΔΔG ) for a protein as well as the relative binding affinity (ΔΔG ) between proteins using a single protein structure as input. However, it can suffer from low accuracy, especially for antibody-antigen systems. In this work, we trained a neural network on FoldX output to enhance its prediction accuracy. We first performed FoldX calculations on the largest datasets available for mutations that affect binding (SKEMPIv2) and folding (ProTherm4) with experimentally measured ΔΔG. Features were then extracted from the FoldX output files including its prediction for ΔΔG. We then developed and optimized a neural network framework to predict the difference between FoldX estimated ΔΔG and the experimental data, creating a model capable of producing a correction factor. Our approach showed significant improvements in Pearson correlation performance. For single mutations affecting folding, the correlation improved from a baseline of 0.3 to 0.66. In terms of binding, performance increased from 0.37 to 0.61 for single mutations and from 0.52 to 0.81 for double mutations. For epistasis, the correlation for binding affinity (both singles and doubles) improved from 0.19 to 0.59. Our results also indicated that models trained on double mutations enhanced accuracy when predicting higher-order mutations (such as triple or quadruple mutations), whereas models trained on singles did not. This suggests that interaction energy and epistasis effects present in the FoldX output are not fully utilized by FoldX itself. Once trained, these models add minimal computational time but provide a substantial increase in performance, especially for higher-order mutations and epistasis. This makes them a valuable addition to any free energy prediction pipeline using FoldX. Furthermore, we believe this technique can be further optimized and tested for predicting antibody escape, aiding in the efficient development of watch lists.

Improvements of Network Approach for Analysis of the Folding Free‐energy Surface of Peptides and Proteins

Network Approach to Identify the Folding Transition States of Peptides and Proteins.

Folding rate prediction based on neural network model

Simplification of complex networks by folding

Investigation of the Proteins Folding Rates and Their Properties of Amino Acid Networks

Prediction of Protein Folding Rates from Structural Topology and Complex Network Properties

A WAVELET APPROACH FOR THE ANALYSIS OF FOLDING TRAJECTORY OF PROTEIN TRP-CAGE

Folding with a protein's native shortcut network

Networks of electrostatic and hydrophobic interactions modulate the complex folding free energy surface of a designed βα protein

Improvement on A Simplified Model for Protein Folding Simulation

Characterization of folding mechanisms of Trp-cage and WW-domain by network analysis of simulations with a hybrid-resolution model.

Protein folding optimization based on 3D off-lattice model via an improved artificial bee colony algorithm

Exploration Of Two-Dimensional Hydrophobic-Polar Lattice Model By Combining Local Search With Elastic Net Algorithm

Folding of Small Helical Proteins Assisted by Small-Angle X-ray Scattering Profiles.

Energy Optimization for Off-Lattice Protein Folding

Static and dynamic characteristics of protein contact networks

Ab Initio Folding of Terminal Segments with Secondary Structures Reveals the Fine Difference Between Two Closely Related All‐atom Statistical Energy Functions

A New-Protein Folding Recognition Potential Function

Leveraging neural networks to correct FoldX free energy estimates

Protein folding simulations: from coarse-grained model to all-atom model.

Multi-Dimensional Theory of Protein Folding