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.

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

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

Binding Energy and Free Energy of Calcium Ion to Calmodulin EF-Hands with the Drude Polarizable Force Field

The Effect of Macromolecular Crowding, Ionic Strength and Calcium Binding on Calmodulin Dynamics

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

Investigating the high affinity and low sequence specificity of calmodulin binding to its targets

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

Mathematical modeling and biochemical analysis support partially ordered calmodulin-myosin light chain kinase binding

Prediction of Three Dimensional Structure of Calmodulin

Homology Modeling Identifies C-Terminal Residues That Contribute to the Ca2+ Sensitivity of A Bkca Channel

Molecular Mechanisms of Calmodulin's Functional Versatility

Computational insights into the binding pattern of mitochondrial calcium uniporter inhibitor through homology modeling, molecular dynamics simulation, binding free energy prediction and density functional theory calculation

Mobility capillary electrophoresis-native mass spectrometry reveals dynamic conformational equilibrium of calmodulin and its complexes

Modulation of Calmodulin Lobes by Different Targets: An Allosteric Model with Hemiconcerted Conformational Transitions

Molecular Insights into the MLCK Activation by CaM

Molecular Mechanism Of Multispecific Recognition Of Calmodulin Through Conformational Changes

Determining the atomic charge of calcium ion in a calcium-binding protein requires the information of its dynamic coordination geometry in aqueous solution

Effect of Ca2+ Binding States of Calmodulin on the Conformational Dynamics and Force Responses of Myosin Lever Arm

Divers models of divalent cation interaction to calcium-binding proteins: techniques and anthology

Decreased Interactions between Calmodulin and a Mutant Huntingtin Model Might Reduce the Cytotoxic Level of Intracellular Ca2+: A Molecular Dynamics Study

Coarse-grained molecular simulations of allosteric cooperativity