Modeling and fitting protein-protein complexes to predict change of binding energy

Insights into Transient Dimerisation of Carnitine/Acylcarnitine Carrier (SLC25A20) from Sarkosyl/PAGE, Cross-Linking Reagents, and Comparative Modelling Analysis

The carnitine/acylcarnitine carrier (CAC) is a crucial protein for cellular energy metabolism, facilitating the exchange of acylcarnitines and free carnitine across the mitochondrial membrane, thereby enabling fatty acid β-oxidation and oxidative phosphorylation (OXPHOS). Although CAC has not been crystallised, structural insights are derived from the mitochondrial ADP/ATP carrier (AAC) structures in both cytosolic and matrix conformations. These structures underpin a single binding centre-gated pore mechanism, a common feature among mitochondrial carrier (MC) family members. The functional implications of this mechanism are well-supported, yet the structural organization of the CAC, particularly the formation of dimeric or oligomeric assemblies, remains contentious. Recent investigations employing biochemical techniques on purified and reconstituted CAC, alongside molecular modelling based on crystallographic AAC dimeric structures, suggest that CAC can indeed form dimers. Importantly, this dimerization does not alter the transport mechanism, a phenomenon observed in various other membrane transporters across different protein families. This observation aligns with the ping-pong kinetic model, where the dimeric form potentially facilitates efficient substrate translocation without necessitating mechanistic alterations. The presented findings thus contribute to a deeper understanding of CAC's functional dynamics and its structural parallels with other MC family members.

Mining the Protein Data Bank to improve prediction of changes in protein-protein binding

Predicting the effect of mutations on protein-protein interactions is important for relating structure to function, as well as for in silico affinity maturation. The effect of mutations on protein-protein binding energy (ΔΔG) can be predicted by a variety of atomic simulation methods involving full or limited flexibility, and explicit or implicit solvent. Methods which consider only limited flexibility are naturally more economical, and many of them are quite accurate, however results are dependent on the atomic coordinate set used. In this work we perform a sequence and structure based search of the Protein Data Bank to find additional coordinate sets and repeat the calculation on each. The method increases precision and Positive Predictive Value, and decreases Root Mean Square Error, compared to using single structures. Given the ongoing growth of near-redundant structures in the Protein Data Bank, our method will only increase in applicability and accuracy.

SSIPe: accurately estimating protein-protein binding affinity change upon mutations using evolutionary profiles in combination with an optimized physical energy function

MOTIVATION

Most proteins perform their biological functions through interactions with other proteins in cells. Amino acid mutations, especially those occurring at protein interfaces, can change the stability of protein-protein interactions (PPIs) and impact their functions, which may cause various human diseases. Quantitative estimation of the binding affinity changes (ΔΔGbind) caused by mutations can provide critical information for protein function annotation and genetic disease diagnoses.

RESULTS

We present SSIPe, which combines protein interface profiles, collected from structural and sequence homology searches, with a physics-based energy function for accurate ΔΔGbind estimation. To offset the statistical limits of the PPI structure and sequence databases, amino acid-specific pseudocounts were introduced to enhance the profile accuracy. SSIPe was evaluated on large-scale experimental data containing 2204 mutations from 177 proteins, where training and test datasets were stringently separated with the sequence identity between proteins from the two datasets below 30%. The Pearson correlation coefficient between estimated and experimental ΔΔGbind was 0.61 with a root-mean-square-error of 1.93 kcal/mol, which was significantly better than the other methods. Detailed data analyses revealed that the major advantage of SSIPe over other traditional approaches lies in the novel combination of the physical energy function with the new knowledge-based interface profile. SSIPe also considerably outperformed a former profile-based method (BindProfX) due to the newly introduced sequence profiles and optimized pseudocount technique that allows for consideration of amino acid-specific prior mutation probabilities.

AVAILABILITY AND IMPLEMENTATION

Web-server/standalone program, source code and datasets are freely available at https://zhanglab.ccmb.med.umich.edu/SSIPe and https://github.com/tommyhuangthu/SSIPe.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

ProAffiMuSeq: sequence-based method to predict the binding free energy change of protein-protein complexes upon mutation using functional classification

MOTIVATION

Protein-protein interactions are essential for the cell and mediate various functions. However, mutations can disrupt these interactions and may cause diseases. Currently available computational methods require a complex structure as input for predicting the change in binding affinity. Further, they have not included the functional class information for the protein-protein complex. To address this, we have developed a method, ProAffiMuSeq, which predicts the change in binding free energy using sequence-based features and functional class.

RESULTS

Our method shows an average correlation between predicted and experimentally determined ΔΔG of 0.73 and mean absolute error (MAE) of 0.86 kcal/mol in 10-fold cross-validation and correlation of 0.75 with MAE of 0.94 kcal/mol in the test dataset. ProAffiMuSeq was also tested on an external validation set and showed results comparable to structure-based methods. Our method can be used for large-scale analysis of disease-causing mutations in protein-protein complexes without structural information.

AVAILABILITY AND IMPLEMENTATION

Users can access the method at https://web.iitm.ac.in/bioinfo2/proaffimuseq/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

In Silico Analysis of Homologous Heterodimers of Cruzipain-Chagasin from Structural Models Built by Homology

The present study gives an overview of the binding energetics of the homologous heterodimers of cruzipain−chagasin based on the binding energy (ΔG_b) prediction obtained with FoldX. This analysis involves a total of 70 homologous models of the cruzipain−chagasin complex which were constructed by homology from the combinatory variation of nine papain-like cysteine peptidase structures and seven cysteine protease inhibitor structures (as chagasin-like and cystatin-like inhibitors). Only 32 systems have been evaluated experimentally, ΔG_b^experimental values previously reported. Therefore, the result of the multiple analysis in terms of the thermodynamic parameters, are shown as relative energy |ΔΔG| = |ΔG_b^fromFoldX − ΔG_b^experimental|. Nine models were identified that recorded |ΔΔG| < 1.3, five models to 2.8 > |ΔΔG| > 1.3 and the other 18 models, values of |ΔΔG| > 2.8. The energetic analysis of the contribution of ΔH and ΔS to ΔG_b to the 14-molecular model presents a ΔG_b mostly ΔH-driven at neutral pH and at an ionic strength (I) of 0.15 M. The dependence of ΔG_b(I,pH) at 298 K to the cruzipain−chagasin complex predicts a linear dependence of ΔG_b(I). The computational protocol allowed the identification and prediction of thermodynamics binding energy parameters for cruzipain−chagasin-like heterodimers.

Tools for the cryo-EM gold rush: going from the cryo-EM map to the atomistic model

As cryo-electron microscopy (cryo-EM) enters mainstream structural biology, the demand for fitting methods is high. Here, we review existing flexible fitting methods for cryo-EM. We discuss their importance, potential concerns and assessment strategies. We aim to give readers concrete descriptions of cryo-EM flexible fitting methods with corresponding examples.

BindProfX: Assessing Mutation-Induced Binding Affinity Change by Protein Interface Profiles with Pseudo-Counts

Understanding how gene-level mutations affect the binding affinity of protein-protein interactions is a key issue of protein engineering. Due to the complexity of the problem, using physical force field to predict the mutation-induced binding free-energy change remains challenging. In this work, we present a renewed approach to calculate the impact of gene mutations on the binding affinity through the structure-based profiling of protein-protein interfaces, where the binding free-energy change (ΔΔG) is counted as the logarithm of relative probability of mutant amino acids over wild-type ones in the interface alignment matrix; three pseudo-counts are introduced to alleviate the limit of the current interface library. Compared with a previous profile score that was based on the log-odds likelihood calculation, the correlation between predicted and experimental ΔΔG of single-site mutations is increased in this approach from 0.33 to 0.68. The structure-based profile score is found complementary to the physical potentials, where a linear combination of the profile score with the FoldX potential could increase the ΔΔG correlation from 0.46 to 0.74. It is also shown that the profile score is robust for counting the coupling effect of multiple individual mutations. For the mutations involving more than two mutation sites where the correlation between FoldX and experimental data vanishes, the profile-based calculation retains a strong correlation with the experimental measurements.