scholarly journals Exploring the Binding Mechanism between Human Profilin (PFN1) and Polyproline-10 through Binding Mode Screening

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.

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 26 (21) ◽  
pp. 3890-3910 ◽  
Author(s):  
Branislava Gemovic ◽  
Neven Sumonja ◽  
Radoslav Davidovic ◽  
Vladimir Perovic ◽  
Nevena Veljkovic

Background: The significant number of protein-protein interactions (PPIs) discovered by harnessing concomitant advances in the fields of sequencing, crystallography, spectrometry and two-hybrid screening suggests astonishing prospects for remodelling drug discovery. The PPI space which includes up to 650 000 entities is a remarkable reservoir of potential therapeutic targets for every human disease. In order to allow modern drug discovery programs to leverage this, we should be able to discern complete PPI maps associated with a specific disorder and corresponding normal physiology. Objective: Here, we will review community available computational programs for predicting PPIs and web-based resources for storing experimentally annotated interactions. Methods: We compared the capacities of prediction tools: iLoops, Struck2Net, HOMCOS, COTH, PrePPI, InterPreTS and PRISM to predict recently discovered protein interactions. Results: We described sequence-based and structure-based PPI prediction tools and addressed their peculiarities. Additionally, since the usefulness of prediction algorithms critically depends on the quality and quantity of the experimental data they are built on; we extensively discussed community resources for protein interactions. We focused on the active and recently updated primary and secondary PPI databases, repositories specialized to the subject or species, as well as databases that include both experimental and predicted PPIs. Conclusion: PPI complexes are the basis of important physiological processes and therefore, possible targets for cell-penetrating ligands. Reliable computational PPI predictions can speed up new target discoveries through prioritization of therapeutically relevant protein–protein complexes for experimental studies.


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.


Polymers ◽  
2021 ◽  
Vol 13 (8) ◽  
pp. 1243
Author(s):  
Fan Zhang ◽  
Yufei Cao ◽  
Xuan Liu ◽  
Huan Xu ◽  
Diannan Lu ◽  
...  

Understanding the aging mechanism of polypropylene (PP) is fundamental for the fabrication and application of PP-based materials. In this paper, we present our study in which we first used reactive molecular dynamics (RMD) simulations to explore the thermo-oxidative aging of PP in the presence of acetic acid or acetone. We studied the effects of temperature and oxygen on the aging process and discussed the formation pathways of typical small molecule products (H2, CO, CO2, CH4, C2H4, and C2H6). The effect of two infection agents, acetic acid and acetone, on the aging reaction was analyzed emphatically. The simulation results showed that acetone has a weak impact on accelerating the aging process, while acetic acid has a significant effect, consistent with previous experimental studies. By tracking the simulation trajectories, both acetic acid and acetone produced small active free radicals to further react with other fragment products, thus accelerating the aging process. The first reaction step of acetic acid is often the shedding of the H atom on the hydroxyl group, while the reaction of acetone is often the shedding of the H atom or the methyl. The latter requires higher energy at lower temperatures. This is why the acceleration effect of acetone for the thermo-oxidative aging of PP was not so significant compared to acetic acid in the experimental temperature (383.15 K).


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Renu Wadhwa ◽  
Neetu Singh Yadav ◽  
Shashank P. Katiyar ◽  
Tomoko Yaguchi ◽  
Chohee Lee ◽  
...  

AbstractPoor bioavailability due to the inability to cross the cell membrane is one of the major reasons for the failure of a drug in clinical trials. We have used molecular dynamics simulations to predict the membrane permeability of natural drugs—withanolides (withaferin-A and withanone) that have similar structures but remarkably differ in their cytotoxicity. We found that whereas withaferin-A, could proficiently transverse through the model membrane, withanone showed weak permeability. The free energy profiles for the interaction of withanolides with the model bilayer membrane revealed that whereas the polar head group of the membrane caused high resistance for the passage of withanone, the interior of the membrane behaves similarly for both withanolides. The solvation analysis further revealed that the high solvation of terminal O5 oxygen of withaferin-A was the major driving force for its high permeability; it interacted with the phosphate group of the membrane that led to its smooth passage across the bilayer. The computational predictions were tested by raising and recruiting unique antibodies that react to withaferin-A and withanone. The time-lapsed analyses of control and treated cells demonstrated higher permeation of withaferin-A as compared to withanone. The concurrence between the computation and experimental results thus re-emphasised the use of computational methods for predicting permeability and hence bioavailability of natural drug compounds in the drug development process.


1991 ◽  
Vol 46 (7) ◽  
pp. 616-620 ◽  
Author(s):  
Junko Habasaki

MD simulation has been performed to learn the microscopic mechanism of diffusion of ions in the Li2SiO3 system. The motion of lithium ions can be explained by the trapping model, where lithium is trapped in the polyhedron and moves with fluctuation of the coordination number. The mean square displacement of lithium was found to correlate well with the net changes in coordination number.


2018 ◽  
Vol 19 (11) ◽  
pp. 3524 ◽  
Author(s):  
Guodong Hu ◽  
Xiu Yu ◽  
Yunqiang Bian ◽  
Zanxia Cao ◽  
Shicai Xu ◽  
...  

ToxIN is a triangular structure formed by three protein toxins (ToxNs) and three specific noncoding RNA antitoxins (ToxIs). To respond to stimuli, ToxI is preferentially degraded, releasing the ToxN. Thus, the dynamic character is essential in the normal function interactions between ToxN and ToxI. Here, equilibrated molecular dynamics (MD) simulations were performed to study the stability of ToxN and ToxI. The results indicate that ToxI adjusts the conformation of 3′ and 5′ termini to bind to ToxN. Steered molecular dynamics (SMD) simulations combined with the recently developed thermodynamic integration in 3nD (TI3nD) method were carried out to investigate ToxN unbinding from the ToxIN complex. The potentials of mean force (PMFs) and atomistic pictures suggest the unbinding mechanism as follows: (1) dissociation of the 5′ terminus from ToxN, (2) missing the interactions involved in the 3′ terminus of ToxI without three nucleotides (G31, A32, and A33), (3) starting to unfold for ToxI, (4) leaving the binding package of ToxN for three nucleotides of ToxI, (5) unfolding of ToxI. This work provides information on the structure-function relationship at the atomistic level, which is helpful for designing new potent antibacterial drugs in the future.


2021 ◽  
Author(s):  
Fréderic Célerse ◽  
Theo Jaffrelot-Inizan ◽  
Louis Lagardère ◽  
Olivier Adjoua ◽  
Pierre Monmarché ◽  
...  

We detail a novel multi-level enhanced sampling strategy grounded on Gaussian accelerated Molecular Dynamics (GaMD). First, we propose a GaMD multi-GPUs-accelerated implementation within the Tinker-HP molecular dynamics package. We then introduce the new "dual-water" mode and its use with the flexible AMOEBA polarizable force field. By adding harmonic boosts to the water stretching and bonding terms, it accelerates the solvent-solute interactions while enabling speedups thanks to the use of fast multiple--timestep integrators. To further reduce time-to-solution, we couple GaMD to Umbrella Sampling (US). The GaMD—US/dual-water approach is tested on the 1D Potential of Mean Force (PMF) of the CD2-CD58 system (168000 atoms) allowing the AMOEBA PMF to converge within 1 kcal/mol of the experimental value. Finally, Adaptive Sampling (AS) is added enabling AS-GaMD capabilities but also the introduction of the new Adaptive Sampling--US--GaMD (ASUS--GaMD) scheme. The highly parallel ASUS--GaMD setup decreases time to convergence by respectively 10 and 20 compared to GaMD--US and US.


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.


Sign in / Sign up

Export Citation Format

Share Document