Abstract:Mutation in calcium-binding proteins (CBPs) can significantly influence Ca binding affinity (BA), resulting in substantial impairment in the signaling process and leading to several lethal diseases. The knowledge behind the changes in the binding affinity can help in understanding the signaling process and designing inhibitors for therapeutic usage. However, accurate prediction of BA for a large number of mutations has been elusive. In this work, for an important calcium binding protein, cardiac Troponin-C, we have developed an integrative modeling approach that combines molecular dynamics (MD) based binding free energy calculations, prediction of plausible mutants using evolutionary information, and an interpretable machine learning model to predict Ca BA for a large number of mutations (seventy-six in all). For the binding free energy calculation, we have used a charge-scaling based MD simulation that considers the polarization in the system, which is critical for divalent ion binding with proteins. The well-known molecular mechanics Poisson-Boltzmann surface area (MMPBSA) method was used for the binding free energy calculations. The calculated results for twenty-four disease mutants, which are associated with different cardiomyopathies and have experimental binding affinity, are in close agreement with the experimental results. To study other plausible mutations, we have probed the evolutionary landscape of cardiac Troponin-C and used the EVmutation method of Hopf et al.(Nature biotechnology 2017, 35, 128-135) to generate sixty-one additional mutants. Finally, a Support vector regression model was developed for both observed and plausible mutations. Our machine learning model used simple structure and sequence-based descriptors along with MD-based descriptors and gave a mean squared error (MSE) of only 0.16 kcal/mol. Assessment of the contribution of each descriptor shows that the number of water molecules within the Ca binding site, type of amino acid substitution (e.g. polar to hydrophobic reduces the binding affinity), and the distance of mutation with Ca are the most important factors in determining the binding affinity. This integrative modeling can be used for other CBPs and can lay the path for modeling the complex and astronomically large mutational landscape of Calcium binding proteins.

Prediction of Ca binding site in proteins with a fast and accurate method based on statistical mechanics and analysis of crystal structures

Prediction of Ca2+ Binding Site in Proteins With a Fast and Accurate Method Based on Statistical Mechanics and Analysis of Crystal Structures

Inferring metal binding sites in flexible regions of proteins

Calcium Ion Binding to the Mutants of Calmodulin: A Structure-Based Computational Predictive Model of Binding Affinity Using a Charge Scaling Approach in Molecular Dynamics Simulation

Quintessence with two energy scales

Machine Learning Meets Physics-based Modeling: A Mass-spring System to Predict Protein-ligand Binding Affinity

Accurate Prediction for Protein-Peptide Binding Based on High-Temperature Molecular Dynamics Simulations

Predicting Ca2+ and Mg2+ ligand binding sites by deep neural network algorithm

Predicting metal-binding proteins and structures through integration of evolutionary-scale and physics-based modeling

Phosphate Binding Sites Prediction in Phosphorylation-Dependent Protein-Protein Interactions

A machine learning approach to predicting protein–ligand binding affinity with applications to molecular docking

Are predicted structures good enough to preserve functional sites?

Computational prediction of protein–protein binding affinities

Machine learning accelerates MD-based binding pose prediction between ligands and proteins

Accurate prediction of dynamic protein–ligand binding using P‐score ranking

Predicting metal-protein interactions using cofolding methods: Status quo

Colabind: A Cloud-Based Approach for Prediction of Binding Sites Using Coarse-Grained Simulations with Molecular Probes

M-Ionic: prediction of metal-ion-binding sites from sequence using residue embeddings

Calcium Binding Affinity in the Mutational Landscape of Troponin-C: Free Energy Calculation, Coevolution modeling and Machine Learning

Template‐guided method for protein–ligand complex structure prediction: Application to CASP15 protein–ligand studies

Reliable and Practical Computational Prediction of Molecular Crystal Polymorphs