scholarly journals Ranking protein-protein docking results using steered molecular dynamics and potential of mean force calculations

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


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.


2018 ◽  
Author(s):  
Leili Zhang ◽  
David R. Bell ◽  
Binquan Luan ◽  
Ruhong Zhou

AbstractThe large magnitude of protein-protein interaction (PPI) pairs within the human interactome necessitates the development of predictive models and screening tools to better understand this fundamental molecular communication. However, despite enormous efforts from various groups to develop predictive techniques in the last decade, PPI complex structures are in general still very challenging to predict due to the large number of degrees of freedom. In this study, we use the binding complex of human profilin (PFN1) and polyproline-10 (P10) as a model system to examine various approaches, with the aim of going beyond normal protein docking for PPI prediction and evaluation. The potential of mean force (PMF) was first obtained from the timeconsuming umbrella sampling, which confirmed that the most stable binding structure identified by the maximal PMF difference is indeed the crystallographic binding structure. Moreover, crucial residues previously identified in experimental studies, W3, H133 and S137 of PFN1, were found to form favorable hydrogen bonds with P10, suggesting a zipping process during the binding between PFN1 and P10. We then explored both regular molecular dynamics (MD) and steered molecular dynamics (SMD) simulations, seeking for better criteria of ranking the PPI prediction. Despite valuable information obtained from conventional MD simulations, neither the commonly used interaction energy between the two binding parties nor the long-term root mean square displacement (RMSD) correlates well with the PMF results. On the other hand, with a sizable collection of trajectories, we demonstrated that the average rupture work calculated from SMD simulations correlates fairly well with the PMFs (R2 = 0.67), making it a promising PPI screening method.


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.


2021 ◽  
Vol 22 (10) ◽  
pp. 5220
Author(s):  
Jarosław J. Panek ◽  
Joanna Zasada ◽  
Bartłomiej M. Szyja ◽  
Beata Kizior ◽  
Aneta Jezierska

The O-H...N and O-H...O hydrogen bonds were investigated in 10-hydroxybenzo[h]quinoline (HBQ) and benzo[h]quinoline-2-methylresorcinol complex in vacuo, solvent and crystalline phases. The chosen systems contain analogous donor and acceptor moieties but differently coupled (intra- versus intermolecularly). Car–Parrinello molecular dynamics (CPMD) was employed to shed light onto principle components of interactions responsible for the self-assembly. It was applied to study the dynamics of the hydrogen bonds and vibrational features as well as to provide initial geometries for incorporation of quantum effects and electronic structure studies. The vibrational features were revealed using Fourier transformation of the autocorrelation function of atomic velocity and by inclusion of nuclear quantum effects on the O-H stretching solving vibrational Schrödinger equation a posteriori. The potential of mean force (Pmf) was computed for the whole trajectory to derive the probability density distribution and for the O-H stretching mode from the proton vibrational eigenfunctions and eigenvalues incorporating statistical sampling and nuclear quantum effects. The electronic structure changes of the benzo[h]quinoline-2-methylresorcinol dimer and trimers were studied based on Constrained Density Functional Theory (CDFT) whereas the Electron Localization Function (ELF) method was applied for all systems. It was found that the bridged proton is localized on the donor side in both investigated systems in vacuo. The crystalline phase simulations indicated bridged proton-sharing and transfer events in HBQ. These effects are even more pronounced when nuclear quantization is taken into account, and the quantized Pmf allows the proton to sample the acceptor area more efficiently. The CDFT indicated the charge depletion at the bridged proton for the analyzed dimer and trimers in solvent. The ELF analysis showed the presence of the isolated proton (a signature of the strongest hydrogen bonds) only in some parts of the HBQ crystal simulation. The collected data underline the importance of the intramolecular coupling between the donor and acceptor moieties.


Sign in / Sign up

Export Citation Format

Share Document