scholarly journals Mode of inhibitory binding of epigallocatechin gallate to the ubiquitin-activating enzyme Uba1 via accelerated molecular dynamics

RSC Advances ◽  
2021 ◽  
Vol 11 (14) ◽  
pp. 8264-8276
Author(s):  
Paras Gaur ◽  
Gabriel Fenteany ◽  
Chetna Tyagi

The hinge-like movement of the SCCH domain upon ligand binding closes the ubiquitin binding site and disrupts the interfaces crucial for thioester bond formation.

2015 ◽  
Vol 48 (4) ◽  
pp. 479-487 ◽  
Author(s):  
Kalli Kappel ◽  
Yinglong Miao ◽  
J. Andrew McCammon

AbstractElucidating the detailed process of ligand binding to a receptor is pharmaceutically important for identifying druggable binding sites. With the ability to provide atomistic detail, computational methods are well poised to study these processes. Here, accelerated molecular dynamics (aMD) is proposed to simulate processes of ligand binding to a G-protein-coupled receptor (GPCR), in this case the M3 muscarinic receptor, which is a target for treating many human diseases, including cancer, diabetes and obesity. Long-timescale aMD simulations were performed to observe the binding of three chemically diverse ligand molecules: antagonist tiotropium (TTP), partial agonist arecoline (ARc) and full agonist acetylcholine (ACh). In comparison with earlier microsecond-timescale conventional MD simulations, aMD greatly accelerated the binding of ACh to the receptor orthosteric ligand-binding site and the binding of TTP to an extracellular vestibule. Further aMD simulations also captured binding of ARc to the receptor orthosteric site. Additionally, all three ligands were observed to bind in the extracellular vestibule during their binding pathways, suggesting that it is a metastable binding site. This study demonstrates the applicability of aMD to protein–ligand binding, especially the drug recognition of GPCRs.


Molecules ◽  
2021 ◽  
Vol 26 (5) ◽  
pp. 1250
Author(s):  
Hien T. T. Lai ◽  
Alejandro Giorgetti ◽  
Giulia Rossetti ◽  
Toan T. Nguyen ◽  
Paolo Carloni ◽  
...  

The translocator protein (TSPO) is a 18kDa transmembrane protein, ubiquitously present in human mitochondria. It is overexpressed in tumor cells and at the sites of neuroinflammation, thus representing an important biomarker, as well as a promising drug target. In mammalian TSPO, there are cholesterol–binding motifs, as well as a binding cavity able to accommodate different chemical compounds. Given the lack of structural information for the human protein, we built a model of human (h) TSPO in the apo state and in complex with PK11195, a molecule routinely used in positron emission tomography (PET) for imaging of neuroinflammatory sites. To better understand the interactions of PK11195 and cholesterol with this pharmacologically relevant protein, we ran molecular dynamics simulations of the apo and holo proteins embedded in a model membrane. We found that: (i) PK11195 stabilizes hTSPO structural fold; (ii) PK11195 might enter in the binding site through transmembrane helices I and II of hTSPO; (iii) PK11195 reduces the frequency of cholesterol binding to the lower, N–terminal part of hTSPO in the inner membrane leaflet, while this impact is less pronounced for the upper, C–terminal part in the outer membrane leaflet, where the ligand binding site is located; (iv) very interestingly, cholesterol most frequently binds simultaneously to the so-called CRAC and CARC regions in TM V in the free form (residues L150–X–Y152–X(3)–R156 and R135–X(2)–Y138–X(2)–L141, respectively). However, when the protein is in complex with PK11195, cholesterol binds equally frequently to the CRAC–resembling motif that we observed in TM I (residues L17–X(2)–F20–X(3)–R24) and to CRAC in TM V. We expect that the CRAC–like motif in TM I will be of interest in future experimental investigations. Thus, our MD simulations provide insight into the structural features of hTSPO and the previously unknown interplay between PK11195 and cholesterol interactions with this pharmacologically relevant protein.


2020 ◽  
Vol 12 (13) ◽  
pp. 1213-1225 ◽  
Author(s):  
Shristi Pawnikar ◽  
Yinglong Miao

Background: Chemokine GPCRs play key roles in biology and medicine. Particularly, CXCR4 promotes cancer metastasis and facilitate HIV entry into host cells. Plerixafor (PLX) is a CXCR4 drug, but the pathway and binding site of PLX in CXCR4 remain unknown. Results & methodology: We have performed molecular docking and all-atom simulations using Gaussian accelerated molecular dynamics (GaMD), which are consistent with previous mutation experiments, suggesting that PLX binds to the orthosteric site of CXCR4 as an antagonist. The GaMD simulations further revealed an intermediate allosteric binding site at the extracellular mouth of CXCR4. Conclusion: The newly identified allosteric site can be targeted for novel drug design targeting CXCR4 and other chemokine receptors.


2020 ◽  
Author(s):  
Yinglong Miao ◽  
Apurba Bhattarai ◽  
Jinan Wang

AbstractCalculations of ligand binding free energies and kinetic rates are important for drug design. However, such tasks have proven challenging in computational chemistry and biophysics. To address this challenge, we have developed a new computational method “LiGaMD”, which selectively boosts the ligand non-bonded interaction potential energy based on the Gaussian accelerated molecular dynamics (GaMD) enhanced sampling technique. Another boost potential could be applied to the remaining potential energy of the entire system in a dual-boost algorithm (LiGaMD_Dual) to facilitate ligand binding. LiGaMD has been demonstrated on host-guest and protein-ligand binding model systems. Repetitive guest binding and unbinding in the β-cyclodextrin host were observed in hundreds-of-nanosecond LiGaMD simulations. The calculated binding free energies of guest molecules with sufficient sampling agreed excellently with experimental data (< 1.0 kcal/mol error). In comparison with previous microsecond-timescale conventional molecular dynamics simulations, accelerations of ligand kinetic rate constants in LiGaMD simulations were properly estimated using Kramers’ rate theory. Furthermore, LiGaMD allowed us to capture repetitive dissociation and binding of the benzamidine inhibitor in trypsin within 1 μs simulations. The calculated ligand binding free energy and kinetic rate constants compared well with the experimental data. In summary, LiGaMD provides a promising approach for characterizing ligand binding thermodynamics and kinetics simultaneously, which is expected to facilitate computer-aided drug design.


Blood ◽  
2006 ◽  
Vol 108 (11) ◽  
pp. 1805-1805
Author(s):  
Marta Murcia ◽  
Marketa Jirouskova ◽  
Jihong Li ◽  
Barry S. Coller ◽  
Filizola Marta

Abstract Although the role of the β3 MIDAS metal ion in ligand binding to αIIbβ3 is well established, serving as the site of interaction of the ligand Asp residue, the role of the nearby LIMBS metal ion is less well defined. Previous studies suggested a role for the LIMBS in ligand binding. We confirmed this by showing that HEK293 cells expressing normal αIIbβ3 adhered to both immobilized fibrinogen and the RGD-containing venom echistatin in the presence of either Mg++/Ca++ or Mn++, whereas two different αIIbβ3 LIMBS mutants (β3 N215A and D217A) failed to adhere to either protein. In addition, we found that both mutations also increased the binding of mAb AP5, which recognizes a ligand-induced binding site (LIBS) in the β3 PSI domain (normal 7±4% vs N215A 46±12% and D217A 41±20% of mAb anti-αIIb (HIP8) binding; mean±SD, n=6, p<0.05 for both), indicating that the mutations caused allosteric changes in the conformation of the receptor. To define the mechanism(s) by which the LIMBS mutants affect ligand binding, we carried out equilibrium and non-equilibrium (steered) molecular dynamics (MD) simulations of the cyclic peptide ligand eptifibatide in complex with either the fully hydrated normal αIIbβ3 integrin headpiece (PDB 1TY6) or the equivalent β3 D217A mutant, with and without the LIMBS metal ion. Simulations were carried out using the GROMACS package with the OPLS all-atom force-field. During the simulation, the hybrid domain of the D217A mutant demonstrated greater structural fluctuations than the normal αIIbβ3. Although Craig et al. have reported the appearance of a new contact between the RGD peptide ligand Asp carboxyl and the LIMBS metal ion in αVβ3 after 10 ps of a 1 ns simulation, we did not observe the appearance of such an interaction between the eptifibatide carboxyl and the normal αIIbβ3 LIMBS metal ion even after 20 ns. We did, however, observe such an interaction with the LIMBS metal ion in the D217A mutant. This interaction was facilitated by the movement of the LIMBS ~ 2 Å closer to the MIDAS, and was accompanied by rearrangements of the LIMBS coordinating residues D158 and N215. When the D217A mutant simulation was performed in the absence of the LIMBS metal ion, changes in the orientation of E220 were also observed. The D217A mutant demonstrated increased fluctuations in the C177–C184 specificity-determining loop (SDL), which has been implicated in ligand binding, and decreased fluctuations in K209. Steered MD were used to investigate the pulling forces required to unbind eptifibatide from its binding site. Notably, although the unbinding force decreased modestly when the LIMBS metal ion was removed, it required removal of both the LIMBS and MIDAS metal ions to effect a marked reduction in unbinding force. The binding free energies of the association of the αIIb and β3 subunits were also calculated, and the D217A mutant in the presence of the LIMBS metal ion demonstrated much tighter binding than normal integrin αIIbβ3 (ΔGb −162±6 vs −119±6 Kcal/mol; mean±SD; n=500). We conclude that the LIMBS plays a crucial role in ligand binding to αIIbβ3, perhaps by virtue of its effects on the coordination of the MIDAS, the accentuated mobility of specific domains (e.g., the SDL and the hybrid domains), and/or the number and strength of contacts between αIIb and β3.


2014 ◽  
Vol 28 (4) ◽  
pp. 534-545 ◽  
Author(s):  
P.C.T. Souza ◽  
A.C. Puhl ◽  
L. Martínez ◽  
R. Aparício ◽  
A.S. Nascimento ◽  
...  

Abstract Thyroid hormone receptors (TRs) are members of the nuclear receptor superfamily of ligand-activated transcription factors involved in cell differentiation, growth, and homeostasis. Although X-ray structures of many nuclear receptor ligand-binding domains (LBDs) reveal that the ligand binds within the hydrophobic core of the ligand-binding pocket, a few studies suggest the possibility of ligands binding to other sites. Here, we report a new x-ray crystallographic structure of TR-LBD that shows a second binding site for T3 and T4 located between H9, H10, and H11 of the TRα LBD surface. Statistical multiple sequence analysis, site-directed mutagenesis, and cell transactivation assays indicate that residues of the second binding site could be important for the TR function. We also conducted molecular dynamics simulations to investigate ligand mobility and ligand-protein interaction for T3 and T4 bound to this new TR surface-binding site. Extensive molecular dynamics simulations designed to compute ligand-protein dissociation constant indicate that the binding affinities to this surface site are of the order of the plasma and intracellular concentrations of the thyroid hormones, suggesting that ligands may bind to this new binding site under physiological conditions. Therefore, the second binding site could be useful as a new target site for drug design and could modulate selectively TR functions.


Sign in / Sign up

Export Citation Format

Share Document