scholarly journals Protein docking and steered molecular dynamics suggest alternative phospholamban-binding sites on the SERCA calcium transporter

2020 ◽  
Vol 295 (32) ◽  
pp. 11262-11274
Author(s):  
Rebecca F. Alford ◽  
Nikolai Smolin ◽  
Howard S. Young ◽  
Jeffrey J. Gray ◽  
Seth L. Robia

The transport activity of the sarco(endo)plasmic reticulum calcium ATPase (SERCA) in cardiac myocytes is modulated by an inhibitory interaction with a transmembrane peptide, phospholamban (PLB). Previous biochemical studies have revealed that PLB interacts with a specific inhibitory site on SERCA, and low-resolution structural evidence suggests that PLB interacts with distinct alternative sites on SERCA. High-resolution details of the structural determinants of SERCA regulation have been elusive because of the dynamic nature of the regulatory complex. In this study, we used computational approaches to develop a structural model of SERCA–PLB interactions to gain a mechanistic understanding of PLB-mediated SERCA transport regulation. We combined steered molecular dynamics and membrane protein–protein docking experiments to achieve both a global search and all-atom force calculations to determine the relative affinities of PLB for candidate sites on SERCA. We modeled the binding of PLB to several SERCA conformations, representing different enzymatic states sampled during the calcium transport catalytic cycle. The results of the steered molecular dynamics and docking experiments indicated that the canonical PLB-binding site (comprising transmembrane helices M2, M4, and M9) is the preferred site. This preference was even more stringent for a superinhibitory PLB variant. Interestingly, PLB-binding specificity became more ambivalent for other SERCA conformers. These results provide evidence for polymorphic PLB interactions with novel sites on M3 and with the outside of the SERCA helix M9. Our findings are compatible with previous physical measurements that suggest that PLB interacts with multiple binding sites, conferring dynamic responsiveness to changing physiological conditions.

2019 ◽  
Author(s):  
Rebecca F. Alford ◽  
Nikolai Smolin ◽  
Howard S. Young ◽  
Jeffrey J. Gray ◽  
Seth L. Robia

AbstractThe transport activity of the calcium ATPase SERCA is modulated by an inhibitory interaction with a 52-residue transmembrane peptide, phospholamban (PLB). Biochemical and structural studies have revealed the primary inhibitory site on SERCA, but PLB has been hypothesized to interact with alternative sites on SERCA that are distinct from the inhibitory site. The present study was undertaken to test these hypotheses and explore structural determinants of SERCA regulation by PLB. Steered molecular dynamics (SMD) and membrane protein-protein docking experiments were performed to investigate the apparent affinity of PLB interactions with candidate sites on SERCA. We modeled the relative binding of PLB to several different conformations of SERCA, representing different enzymatic states sampled during the calcium transport catalytic cycle. Overall, the SMD and docking experiments suggest that the canonical binding site is preferred, but also provide evidence for alternative sites that are favorable for certain conformational states of SERCA.


2016 ◽  
Vol 37 (20) ◽  
pp. 1861-1865 ◽  
Author(s):  
Laura J. Kingsley ◽  
Juan Esquivel-Rodríguez ◽  
Ying Yang ◽  
Daisuke Kihara ◽  
Markus A. Lill

2019 ◽  
Vol 116 (42) ◽  
pp. 21001-21011 ◽  
Author(s):  
Benjamin W. Parker ◽  
Edward J. Goncz ◽  
David T. Krist ◽  
Alexander V. Statsyuk ◽  
Alexey I. Nesvizhskii ◽  
...  

Short linear peptide motifs that are intracellular ligands of folded proteins are a modular, incompletely understood molecular interaction language in signaling systems. Such motifs, which frequently occur in intrinsically disordered protein regions, often bind partner proteins with modest affinity and are difficult to study with conventional structural biology methods. We developed LiF-MS (ligand-footprinting mass spectrometry), a method to map peptide binding sites on folded protein domains that allows consideration of their dynamic disorder, and used it to analyze a set of D-motif peptide–mitogen-activated protein kinase (MAPK) associations to validate the approach and define unknown binding structures. LiF-MS peptide ligands carry a short-lived, indiscriminately reactive cleavable crosslinker that marks contacts close to ligand binding sites with high specificity. Each marked amino acid provides an independent constraint for a set of directed peptide–protein docking simulations, which are analyzed by agglomerative hierarchical clustering. We found that LiF-MS provides accurate ab initio identification of ligand binding surfaces and a view of potential binding ensembles of a set of D-motif peptide–MAPK associations. Our analysis provides an MKK4–JNK1 structural model, which has thus far been crystallographically unattainable, a potential alternate binding mode for part of the NFAT4–JNK interaction, and evidence of bidirectional association of MKK4 peptide with ERK2. Overall, we find that LiF-MS is an effective noncrystallographic way to understand how short linear motifs associate with specific sites on folded protein domains at the level of individual amino acids.


2019 ◽  
Author(s):  
Frédéric Célerse ◽  
Louis Lagardere ◽  
Étienne Derat ◽  
Jean-Philip Piquemal

This paper is dedicated to the massively parallel implementation of Steered Molecular Dynamics in the Tinker-HP softwtare. It allows for direct comparisons of polarizable and non-polarizable simulations of realistic systems.


2019 ◽  
Author(s):  
Frédéric Célerse ◽  
Louis Lagardere ◽  
Étienne Derat ◽  
Jean-Philip Piquemal

This paper is dedicated to the massively parallel implementation of Steered Molecular Dynamics in the Tinker-HP softwtare. It allows for direct comparisons of polarizable and non-polarizable simulations of realistic systems.


Life ◽  
2021 ◽  
Vol 11 (2) ◽  
pp. 74
Author(s):  
Justin Spiriti ◽  
Chung F. Wong

Most early-stage drug discovery projects focus on equilibrium binding affinity to the target alongside selectivity and other pharmaceutical properties. Since many approved drugs have nonequilibrium binding characteristics, there has been increasing interest in optimizing binding kinetics early in the drug discovery process. As focal adhesion kinase (FAK) is an important drug target, we examine whether steered molecular dynamics (SMD) can be useful for identifying drug candidates with the desired drug-binding kinetics. In simulating the dissociation of 14 ligands from FAK, we find an empirical power–law relationship between the simulated time needed for ligand unbinding and the experimental rate constant for dissociation, with a strong correlation depending on the SMD force used. To improve predictions, we further develop regression models connecting experimental dissociation rate with various structural and energetic quantities derived from the simulations. These models can be used to predict dissociation rates from FAK for related compounds.


Molecules ◽  
2021 ◽  
Vol 26 (4) ◽  
pp. 1051
Author(s):  
Edgardo Becerra ◽  
Giovanny Aguilera-Durán ◽  
Laura Berumen ◽  
Antonio Romo-Mancillas ◽  
Guadalupe García-Alcocer

Multidrug resistance protein-4 (MRP4) belongs to the ABC transporter superfamily and promotes the transport of xenobiotics including drugs. A non-synonymous single nucleotide polymorphisms (nsSNPs) in the ABCC4 gene can promote changes in the structure and function of MRP4. In this work, the interaction of certain endogen substrates, drug substrates, and inhibitors with wild type-MRP4 (WT-MRP4) and its variants G187W and Y556C were studied to determine differences in the intermolecular interactions and affinity related to SNPs using protein threading modeling, molecular docking, all-atom, coarse grained, and umbrella sampling molecular dynamics simulations (AA-MDS and CG-MDS, respectively). The results showed that the three MRP4 structures had significantly different conformations at given sites, leading to differences in the docking scores (DS) and binding sites of three different groups of molecules. Folic acid (FA) had the highest variation in DS on G187W concerning WT-MRP4. WT-MRP4, G187W, Y556C, and FA had different conformations through 25 ns AA-MD. Umbrella sampling simulations indicated that the Y556C-FA complex was the most stable one with or without ATP. In Y556C, the cyclic adenosine monophosphate (cAMP) and ceefourin-1 binding sites are located out of the entrance of the inner cavity, which suggests that both cAMP and ceefourin-1 may not be transported. The binding site for cAMP and ceefourin-1 is quite similar and the affinity (binding energy) of ceefourin-1 to WT-MRP4, G187W, and Y556C is greater than the affinity of cAMP, which may suggest that ceefourin-1 works as a competitive inhibitor. In conclusion, the nsSNPs G187W and Y556C lead to changes in protein conformation, which modifies the ligand binding site, DS, and binding energy.


Sign in / Sign up

Export Citation Format

Share Document