Identification of Cryptic Binding Sites Using MixMD with Standard and Accelerated Molecular Dynamics

Author(s):  
Richard D. Smith ◽  
Heather A. Carlson
2020 ◽  
Vol 53 (3) ◽  
pp. 654-661 ◽  
Author(s):  
Antonija Kuzmanic ◽  
Gregory R. Bowman ◽  
Jordi Juarez-Jimenez ◽  
Julien Michel ◽  
Francesco L. Gervasio

PeerJ ◽  
2020 ◽  
Vol 8 ◽  
pp. e9241
Author(s):  
Hongye Hu ◽  
Buran Chen ◽  
Danni Zheng ◽  
Guanli Huang

Background Research has shown that Poly-ADP-ribose polymerases 1 (PARP-1) is a potential therapeutic target in the clinical treatment of breast cancer. An increasing number of studies have focused on the development of highly selective inhibitors that target PARP-1 over PARP-2, its closest isoform, to mitigate potential side effects. However, due to the highly conserved and similar binding sites of PARP-1 and PARP-2, there is a huge challenge for the discovery and design of PARP-1 inhibitors. Recently, it was reported that a potent PARP-1 inhibitor named NMS-P118 exhibited greater selectivity to PARP-1 over PARP-2 compared with a previously reported drug (Niraparib). However, the mechanisms underlying the effect of this inhibitor remains unclear. Methods In the present study, classical molecular dynamics (MD) simulations and accelerated molecular dynamics (aMD) simulations combined with structural and energetic analysis were used to investigate the structural dynamics and selective mechanisms of PARP-1 and PARP-2 that are bound to NMS-P118 and Niraparib with distinct selectivity. Results The results from classical MD simulations indicated that the selectivity of inhibitors may be controlled by electrostatic interactions, which were mainly due to the residues of Gln-322, Ser-328, Glu-335, and Tyr-455 in helix αF. The energetic differences were corroborated by the results from aMD simulations. Conclusion This study provides new insights about how inhibitors specifically bind to PARP-1 over PARP-2, which may help facilitate the design of highly selective PARP-1 inhibitors in the future.


2017 ◽  
Vol 57 (6) ◽  
pp. 1388-1401 ◽  
Author(s):  
S. Roy Kimura ◽  
Hai Peng Hu ◽  
Anatoly M. Ruvinsky ◽  
Woody Sherman ◽  
Angelo D. Favia

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.


Biomolecules ◽  
2021 ◽  
Vol 11 (2) ◽  
pp. 334
Author(s):  
Shih-Ting Hong ◽  
Yu-Cheng Su ◽  
Yu-Jen Wang ◽  
Tian-Lu Cheng ◽  
Yeng-Tseng Wang

Humira is a monoclonal antibody that binds to TNF alpha, inactivates TNF alpha receptors, and inhibits inflammation. Neonatal Fc receptors can mediate the transcytosis of Humira–TNF alpha complex structures and process them toward degradation pathways, which reduces the therapeutic effect of Humira. Allowing the Humira–TNF alpha complex structures to dissociate to Humira and soluble TNF alpha in the early endosome to enable Humira recycling is crucial. We used the cytoplasmic pH (7.4), the early endosomal pH (6.0), and pKa of histidine side chains (6.0–6.4) to mutate the residues of complementarity-determining regions with histidine. Our engineered Humira (W1-Humira) can bind to TNF alpha in plasma at neutral pH and dissociate from the TNF alpha in the endosome at acidic pH. We used the constant-pH molecular dynamics, Gaussian accelerated molecular dynamics, two-dimensional potential mean force profiles, and in vitro methods to investigate the characteristics of W1-Humira. Our results revealed that the proposed Humira can bind TNF alpha with pH-dependent affinity in vitro. The W1-Humira was weaker than wild-type Humira at neutral pH in vitro, and our prediction results were close to the in vitro results. Furthermore, our approach displayed a high accuracy in antibody pH-dependent binding characteristics prediction, which may facilitate antibody drug design. Advancements in computational methods and computing power may further aid in addressing the challenges in antibody drug design.


Author(s):  
Nikolay Kondratyuk ◽  
Vsevolod Nikolskiy ◽  
Daniil Pavlov ◽  
Vladimir Stegailov

Classical molecular dynamics (MD) calculations represent a significant part of the utilization time of high-performance computing systems. As usual, the efficiency of such calculations is based on an interplay of software and hardware that are nowadays moving to hybrid GPU-based technologies. Several well-developed open-source MD codes focused on GPUs differ both in their data management capabilities and in performance. In this work, we analyze the performance of LAMMPS, GROMACS and OpenMM MD packages with different GPU backends on Nvidia Volta and AMD Vega20 GPUs. We consider the efficiency of solving two identical MD models (generic for material science and biomolecular studies) using different software and hardware combinations. We describe our experience in porting the CUDA backend of LAMMPS to ROCm HIP that shows considerable benefits for AMD GPUs comparatively to the OpenCL backend.


Sign in / Sign up

Export Citation Format

Share Document