Dasatinib-SIK2 Binding Elucidated by Homology Modeling, Molecular Docking, and Dynamics Simulations.
Ontology highlight
ABSTRACT:
SUBMITTER: Shi M
PROVIDER: S-EPMC8153941 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:BackgroundAn important mechanism of endocrine activity is chemicals entering target cells via transport proteins and then interacting with hormone receptors such as the estrogen receptor (ER). α-Fetoprotein (AFP) is a major transport protein in rodent serum that can bind and sequester estrogens, thus preventing entry to the target cell and where they could otherwise induce ER-mediated endocrine activity. Recently, we reported rat AFP binding affinities for a large set of structurally diverse chemicals, including 53 binders and 72 non-binders. However, the lack of three-dimensional (3D) structures of rat AFP hinders further understanding of the structural dependence for binding. Therefore, a 3D structure of rat AFP was built using homology modeling in order to elucidate rat AFP-ligand binding modes through docking analyses and molecular dynamics (MD) simulations.MethodsHomology modeling was first applied to build a 3D structure of rat AFP. Molecular docking and Molecular Mechanics-Generalized Born Surface Area (MM-GBSA) scoring were then used to examine potential rat AFP ligand binding modes. MD simulations and free energy calculations were performed to refine models of binding modes.ResultsA rat AFP tertiary structure was first obtained using homology modeling and MD simulations. The rat AFP-ligand binding modes of 13 structurally diverse, representative binders were calculated using molecular docking, (MM-GBSA) ranking and MD simulations. The key residues for rat AFP-ligand binding were postulated through analyzing the binding modes.ConclusionThe optimized 3D rat AFP structure and associated ligand binding modes shed light on rat AFP-ligand binding interactions that, in turn, provide a means to estimate binding affinity of unknown chemicals. Our results will assist in the evaluation of the endocrine disruption potential of chemicals.
| S-EPMC3851483 | biostudies-literature
Project description:We present a homology based model of the ligand binding domain (LBD) of the homopentameric alpha1 glycine receptor (GlyR). The model is based on multiple sequence alignment with other members of the nicotinicoid ligand gated ion channel superfamily and two homologous acetylcholine binding proteins (AChBP) from the freshwater (Lymnaea stagnalis) and saltwater (Aplysia californica) snails with known high resolution structure. Using two template proteins with known structure to model three dimensional structure of a target protein is especially advantageous for sequences with low homology as in the case presented in this paper. The final model was cross-validated by critical evaluation of experimental and published mutagenesis, functional and other biochemical studies. In addition, a complex structure with strychnine antagonist in the putative binding site is proposed based on docking simulation using Autodock program. Molecular dynamics (MD) simulations with simulated annealing protocol are reported on the proposed LBD of GlyR, which is stable in 5 ns simulation in water, as well as for a deformed LBD structure modeled on the corresponding domain determined in low-resolution cryomicroscopy structure of the alpha subunit of the full-length acetylcholine receptor (AChR). Our simulations demonstrate that the beta-sandwich central core of the protein monomer is fairly rigid in the simulations and resistant to deformations in water.
| S-EPMC2873197 | biostudies-literature
Project description:We have modeled the structure of human lymphotactin (hLpnt), by homology modeling and molecular dynamics simulations. This chemokine is unique in having a single disulfide bond and a long C-terminal tail. Because other structural classes of chemokines have two pairs of Cys residues, compared to one in Lpnt, and because it has been shown that both disulfide bonds are required for stability and function, the question arises how the Lpnt maintains its structural integrity. The initial structure of hLpnt was constructed by homology modeling. The first 63 residues in the monomer of hLpnt were modeled using the structure of the human CC chemokine, RANTES, whose sequence appeared most similar. The structure of the long C-terminal tail, missing in RANTES, was taken from the human muscle fatty-acid binding protein. In a Protein Data Bank search, this protein was found to contain a sequence that was most homologous to the long tail. Consequently, the modeled hLpnt C-terminal tail consisted of both alpha-helical and beta-motifs. The complete model of the hLpnt monomer consisted of two alpha-helices located above the five-stranded beta-sheet. Molecular dynamics simulations of the solvated initial model have indicated that the stability of the predicted fold is related to the geometry of Pro78. The five-stranded beta-sheet appeared to be preserved only when Pro78 was modeled in the cis conformation. Simulations were also performed both for the C-terminal truncated forms of the hLpnt that contained one or two (CC chemokine-like) disulfide bonds, and for the chicken Lpnt (cLpnt). Our MD simulations indicated that the turn region (T30-G34) in hLpnt is important for the interactions with the receptor, and that the long C-terminal region stabilizes both the turn (T30-G34) and the five-stranded beta-sheet. The major conclusion from our theoretical studies is that the lack of one disulfide bond and the extension of the C-terminus in hLptn are mutually complementary. It is very likely that removal of two Cys residues sufficiently destabilizes the structure of a chemokine molecule, particularly the core beta-sheet, to abolish its biological function. However, this situation is rectified by the long C-terminal segment. The role of this long region is most likely to stabilize the first beta-turn region and alpha-helix H1, explaining how this chemokine can function with a single disulfide bond.
| S-EPMC2144500 | biostudies-other
Project description:Peptides mediate up to 40% of known protein-protein interactions in higher eukaryotes and play a key role in cellular signaling, protein trafficking, immunology, and oncology. However, it is challenging to predict peptide-protein binding with conventional computational modeling approaches, due to slow dynamics and high peptide flexibility. Here, we present a prototype of the approach which combines global peptide docking using ClusPro PeptiDock and all-atom enhanced simulations using Gaussian accelerated molecular dynamics (GaMD). For three distinct model peptides, the lowest backbone root-mean-square deviations (RMSDs) of their bound conformations relative to X-ray structures obtained from PeptiDock were 3.3-4.8 Å, being medium quality predictions according to the Critical Assessment of PRediction of Interactions (CAPRI) criteria. GaMD simulations refined the peptide-protein complex structures with significantly reduced peptide backbone RMSDs of 0.6-2.7 Å, yielding two high quality (sub-angstrom) and one medium quality models. Furthermore, the GaMD simulations identified important low-energy conformational states and revealed the mechanism of peptide binding to the target proteins. Therefore, PeptiDock+GaMD is a promising approach for exploring peptide-protein interactions.
| S-EPMC6835073 | biostudies-literature
Project description:Antibody based bio-molecular drugs are an exciting, new avenue of drug development as an alternative to the more traditional small chemical compounds. However, the binding mechanism and the effect on the conformational ensembles of a therapeutic antibody to its peptide or protein antigen have not yet been well studied. We have utilized dynamic docking and path sampling simulations based on all-atom molecular dynamics to study the binding mechanism between the antibody solanezumab and the peptide amyloid-β (Aβ). Our docking simulations reproduced the experimental structure and gave us representative binding pathways, from which we accurately estimated the binding free energy. Not only do our results show why solanezumab has an explicit preference to bind to the monomeric form of Aβ, but that upon binding, both molecules are stabilized towards a specific conformation, suggesting that their complex formation follows a novel, mutual population-shift model, where upon binding, both molecules impact the dynamics of their reciprocal one.
| S-EPMC6989527 | biostudies-literature
Project description:By performing homology modeling, molecular docking, and molecular dynamics (MD) simulations, we have developed three-dimensional (3D) structural models of the M5 muscarinic acetylcholine receptor (mAChR) and two complexes for M5 mAChR binding with antagonists SVT-40776 and solifenacin in the environment of lipid bilayer and solvent water. According to the simulated results, each of the antagonists is oriented horizontally in the binding pocket formed by transmembrane helices 2, 3, and 5-7. The cationic headgroup of each of the antagonists interacts with a negatively charged residue, Asp110, through electrostatic and hydrogen-bonding interactions. The simulated results also reveal some significant difference between the binding modes of SVT-40776 and solifenacin. In particular, SVT-40776 is persistently hydrogen bonded with the side chain of residue Tyr458, whereas solifenacin cannot form a similar hydrogen bond with residues around its carbonyl group. Such significant difference in the binding structures is consistent with the fact that SVT-40776 has a much higher binding affinity (K(d) = 0.4 nM) to M5 mAChR than that of solifenacin (K(d) = 31 nM) with the same reeptor. The calculated binding free energy change (-2.3 ± 0.3 kcal/mol) from solifenacin to SVT-40776 is in good agreement with the experimentally derived binding free energy change (-2.58 kcal/mol), suggesting that our modeled M5 mAChR structure and its complexes with the antagonists are reliable. The new structural insights obtained from this computational study are expected to stimulate further biochemical and pharmacological studies on the detailed structures of M5 and other subtypes of mAChRs.
| S-EPMC3257414 | biostudies-literature
Project description:Epididymal protease inhibitor (EPPIN) that is located on the sperm surface and specific to the male reproductive system is a non-hormonal contraceptive target, since the binding of EPPIN with the seminal plasma protein semenogelin-1 (SEMG1) causes a loss of sperm function. Here, we investigated the binding interactions between EPPIN and SEMG1 by homology modeling, docking and molecular dynamics simulation. Since no crystal structure was reported for EPPIN, its 3D structure was constructed by homology modeling and refined by dynamics simulation, illustrating the C-terminus domain of EPPIN could bind with its N-terminus domain through the residues 30-32 and 113-116. The binding interaction of SEMG110-8 peptide and EPPIN was investigated by Z-DOCK and dynamics simulation. After evaluating the models according to the calculated binding free energies, we demonstrated that C-terminus domain of EPPIN was important for the binding of SEMG1 via residues Tyr107, Gly112, Asn116, Gln118 and Asn122, while residue Arg32 in N-terminus domain also had contribution for their binding interaction. Additionally, the binding pocket of EPPIN was defined according to these key residues and verified by molecular docking with reported inhibitor EP055, suggesting that the pocket formed by Arg32, Asn114, Asn116, Phe117 and Asn122 could be important for the design of new ligands. This study might be helpful for the understanding of biological function of EPPIN and would encourage the discovery of non-hormonal contraceptive leads/drugs in the future.
| S-EPMC6686837 | biostudies-literature
Project description:The intricate complex system of the differentiation 47 (CD47) and the signal-regulatory protein alpha (SIRPα) cluster is a crucial target for cancer immunotherapy. Although the conformational state of the CD47-SIRPα complex has been revealed through crystallographic studies, further characterization is needed to fully understand the binding mechanism and to identify the hot spot residues involved. In this study, molecular dynamics (MD) simulations were carried out for the complexes of CD47 with two SIRPα variants (SIRPαv1, SIRPαv2) and the commercially available anti-CD47 monoclonal antibody (B6H12.2). The calculated binding free energy of CD47-B6H12.2 is lower than that of CD47-SIRPαv1 and CD47-SIRPαv2 in all the three simulations, indicating that CD47-B6H12.2 has a higher binding affinity than the other two complexes. Moreover, the dynamical cross-correlation matrix reveals that the CD47 protein shows more correlated motions when it binds to B6H12.2. Significant effects were observed in the energy and structural analyses of the residues (Glu35, Tyr37, Leu101, Thr102, Arg103) in the C strand and FG region of CD47 when it binds to the SIRPα variants. The critical residues (Leu30, Val33, Gln52, Lys53, Thr67, Arg69, Arg95, and Lys96) were identified in SIRPαv1 and SIRPαv2, which surround the distinctive groove regions formed by the B2C, C'D, DE, and FG loops. Moreover, the crucial groove structures of the SIRPα variants shape into obvious druggable sites. The C'D loops on the binding interfaces undergo notable dynamical changes throughout the simulation. For B6H12.2, the residues Tyr32LC, His92LC, Arg96LC, Tyr32HC, Thr52HC, Ser53HC, Ala101HC, and Gly102HC in its initial half of the light and heavy chains exhibit obvious energetic and structural impacts upon binding with CD47. The elucidation of the binding mechanism of SIRPαv1, SIRPαv2, and B6H12.2 with CD47 could provide novel perspectives for the development of inhibitors targeting CD47-SIRPα.
| S-EPMC10304963 | biostudies-literature
Project description:When a small molecule binds to the androgen receptor (AR), a conformational change can occur which impacts subsequent binding of co-regulator proteins and DNA. In order to accurately study this mechanism, the scientific community needs a crystal structure of the Wild type AR (WT-AR) ligand binding domain, bound with antagonist. To address this open need, we leveraged molecular docking and molecular dynamics (MD) simulations to construct a structure of the WT-AR ligand binding domain bound with antagonist bicalutamide. The structure of mutant AR (Mut-AR) bound with this same antagonist informed this study. After molecular docking analysis pinpointed the suitable binding orientation of a ligand in AR, the model was further optimized through 1 μs of MD simulations. Using this approach, three molecular systems were studied: (1) WT-AR bound with agonist R1881, (2) WT-AR bound with antagonist bicalutamide, and (3) Mut-AR bound with bicalutamide. Our structures were very similar to the experimentally determined structures of both WT-AR with R1881 and Mut-AR with bicalutamide, demonstrating the trustworthiness of this approach. In our model, when WT-AR is bound with bicalutamide, Val716/Lys720/Gln733, or Met734/Gln738/Glu897 move and thus disturb the positive and negative charge clumps of the AF2 site. This disruption of the AF2 site is key for understanding the impact of antagonist binding on subsequent co-regulator binding. In conclusion, the antagonist induced structural changes in WT-AR detailed in this study will enable further AR research and will facilitate AR targeting drug discovery.
| S-EPMC5962723 | biostudies-literature
Project description:Dopamine (DA) receptors, a class of G-protein coupled receptors (GPCRs), have been targeted for drug development for the treatment of neurological, psychiatric and ocular disorders. The lack of structural information about GPCRs and their ligand complexes has prompted the development of homology models of these proteins aimed at structure-based drug design. Crystal structure of human dopamine D(3) (hD(3)) receptor has been recently solved. Based on the hD(3) receptor crystal structure we generated dopamine D(2) and D(3) receptor models and refined them with molecular dynamics (MD) protocol. Refined structures, obtained from the MD simulations in membrane environment, were subsequently used in molecular docking studies in order to investigate potential sites of interaction. The structure of hD(3) and hD(2L) receptors was differentiated by means of MD simulations and D(3) selective ligands were discriminated, in terms of binding energy, by docking calculation. Robust correlation of computed and experimental K(i) was obtained for hD(3) and hD(2L) receptor ligands. In conclusion, the present computational approach seems suitable to build and refine structure models of homologous dopamine receptors that may be of value for structure-based drug discovery of selective dopaminergic ligands.
| S-EPMC3435408 | biostudies-literature