Browse
Submit Data
Databases
API
Help

Dataset Information

18 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Binding Selectivity of Inhibitors toward Bromodomains BAZ2A and BAZ2B Uncovered by Multiple Short Molecular Dynamics Simulations and MM-GBSA Calculations.

ABSTRACT: Two Bromodomain-Containing proteins BAZ2A and BAZ2B are responsible for remodeling chromatin and regulating noncoding RNAs. As for our current studies, integration of multiple short molecular dynamics simulations (MSMDSs) with molecular mechanics generalized Born surface area (MM-GBSA) method is adopted for insights into binding selectivity of three small molecules D8Q, D9T and UO1 to BAZ2A against BAZ2B. The calculations of MM-GBSA unveil that selectivity of inhibitors toward BAZ2A and BAZ2B highly depends on the enthalpy changes and the details uncover that D8Q has better selectivity toward BAZ2A than BAZ2B, D9T more favorably bind to BAZ2B than BAZ2A, and UO1 does not show obvious selectivity toward these two proteins. The analysis of interaction network between residues and inhibitors indicates that seven residues are mainly responsible for the selectivity of D8Q, six residues for D9T and four residues provide significant contributions to associations of UO1 with two proteins. Moreover the analysis of interaction network not only reveals warm spots of inhibitor bindings to BAZ2A and BAZ2B but also unveils that common residue pairs, including (W1816, W1887), (P1817, P1888), (F1818, F1889), (V1822, V1893), (N1823, N1894),(L1826, L1897), (V1827, V1898), (F1872, F1943), (N1873, N1944) and (V1879, I1950) belonging to (BAZ2A, BAZ2B), induce mainly binding differences of inhibitors to BAZ2A and BAZ2B. Hence, insights from our current studies offer useful dynamics information relating with conformational alterations and structure-affinity relationship at atomistic levels for novel therapeutic strategies toward BAZ2A and BAZ2B.

SUBMITTER: Wang L

PROVIDER: S-EPMC8154142 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Similar Datasets

Discovery and Characterization of GSK2801, a Selective Chemical Probe for the Bromodomains BAZ2A and BAZ2B.

Project description:Bromodomains are acetyl-lysine specific protein interaction domains that have recently emerged as a new target class for the development of inhibitors that modulate gene transcription. The two closely related bromodomain containing proteins BAZ2A and BAZ2B constitute the central scaffolding protein of the nucleolar remodeling complex (NoRC) that regulates the expression of noncoding RNAs. However, BAZ2 bromodomains have low predicted druggability and so far no selective inhibitors have been published. Here we report the development of GSK2801, a potent, selective and cell active acetyl-lysine competitive inhibitor of BAZ2A and BAZ2B bromodomains as well as the inactive control compound GSK8573. GSK2801 binds to BAZ2 bromodomains with dissociation constants (KD) of 136 and 257 nM for BAZ2B and BAZ2A, respectively. Crystal structures demonstrated a canonical acetyl-lysine competitive binding mode. Cellular activity was demonstrated using fluorescent recovery after photobleaching (FRAP) monitoring displacement of GFP-BAZ2A from acetylated chromatin. A pharmacokinetic study in mice showed that GSK2801 had reasonable in vivo exposure after oral dosing, with modest clearance and reasonable plasma stability. Thus, GSK2801 represents a versatile tool compound for cellular and in vivo studies to understand the role of BAZ2 bromodomains in chromatin biology.

| S-EPMC4770311 | biostudies-literature

Structure enabled design of BAZ2-ICR, a chemical probe targeting the bromodomains of BAZ2A and BAZ2B.

Project description:The bromodomain containing proteins BAZ2A/B play essential roles in chromatin remodeling and regulation of noncoding RNAs. We present the structure based discovery of a potent, selective, and cell active inhibitor 13 (BAZ2-ICR) of the BAZ2A/B bromodomains through rapid optimization of a weakly potent starting point. A key feature of the presented inhibitors is an intramolecular aromatic stacking interaction that efficiently occupies the shallow bromodomain pockets. 13 represents an excellent chemical probe for functional studies of the BAZ2 bromodomains in vitro and in vivo.

| S-EPMC4441536 | biostudies-literature

Assessing the performance of the MM/PBSA and MM/GBSA methods. 1. The accuracy of binding free energy calculations based on molecular dynamics simulations.

Project description:The Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) and the Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) methods calculate binding free energies for macromolecules by combining molecular mechanics calculations and continuum solvation models. To systematically evaluate the performance of these methods, we report here an extensive study of 59 ligands interacting with six different proteins. First, we explored the effects of the length of the molecular dynamics (MD) simulation, ranging from 400 to 4800 ps, and the solute dielectric constant (1, 2, or 4) on the binding free energies predicted by MM/PBSA. The following three important conclusions could be observed: (1) MD simulation length has an obvious impact on the predictions, and longer MD simulation is not always necessary to achieve better predictions. (2) The predictions are quite sensitive to the solute dielectric constant, and this parameter should be carefully determined according to the characteristics of the protein/ligand binding interface. (3) Conformational entropy often show large fluctuations in MD trajectories, and a large number of snapshots are necessary to achieve stable predictions. Next, we evaluated the accuracy of the binding free energies calculated by three Generalized Born (GB) models. We found that the GB model developed by Onufriev and Case was the most successful model in ranking the binding affinities of the studied inhibitors. Finally, we evaluated the performance of MM/GBSA and MM/PBSA in predicting binding free energies. Our results showed that MM/PBSA performed better in calculating absolute, but not necessarily relative, binding free energies than MM/GBSA. Considering its computational efficiency, MM/GBSA can serve as a powerful tool in drug design, where correct ranking of inhibitors is often emphasized.

| S-EPMC3029230 | biostudies-literature

Binding-affinity predictions of HSP90 in the D3R Grand Challenge 2015 with docking, MM/GBSA, QM/MM, and free-energy simulations.

Project description:We have estimated the binding affinity of three sets of ligands of the heat-shock protein 90 in the D3R grand challenge blind test competition. We have employed four different methods, based on five different crystal structures: first, we docked the ligands to the proteins with induced-fit docking with the Glide software and calculated binding affinities with three energy functions. Second, the docked structures were minimised in a continuum solvent and binding affinities were calculated with the MM/GBSA method (molecular mechanics combined with generalised Born and solvent-accessible surface area solvation). Third, the docked structures were re-optimised by combined quantum mechanics and molecular mechanics (QM/MM) calculations. Then, interaction energies were calculated with quantum mechanical calculations employing 970-1160 atoms in a continuum solvent, combined with energy corrections for dispersion, zero-point energy and entropy, ligand distortion, ligand solvation, and an increase of the basis set to quadruple-zeta quality. Fourth, relative binding affinities were estimated by free-energy simulations, using the multi-state Bennett acceptance-ratio approach. Unfortunately, the results were varying and rather poor, with only one calculation giving a correlation to the experimental affinities larger than 0.7, and with no consistent difference in the quality of the predictions from the various methods. For one set of ligands, the results could be strongly improved (after experimental data were revealed) if it was recognised that one of the ligands displaced one or two water molecules. For the other two sets, the problem is probably that the ligands bind in different modes than in the crystal structures employed or that the conformation of the ligand-binding site or the whole protein changes.

| S-EPMC5078160 | biostudies-literature

Assessing the performance of docking scoring function, FEP, MM-GBSA, and QM/MM-GBSA approaches on a series of PLK1 inhibitors.

Project description:Over-expressed polo-like kinases 1, a key regulator of cell mitosis, is associated with carcinogenesis and poor prognosis. It is very necessary to develop a reliable computational affinity prediction protocol targeting PLK1. In this study, the performance of different docking scoring function, free energy perturbation, MM-GBSA and QM/MM-GBSA were evaluated. The ranking capability of FEP is the best with rs = 0.854. However, the rs obtained from MM-GBSA can reach 0.767, which requires only about one-eighth of the simulation time of FEP. As for the sampling method, single long molecular dynamics (SLMD) surpass the multiple short molecular dynamics (MSMD) in ranking of the 20 congeneric compounds by about 0.1 in rs. In addition, ligands treated by QM can significantly improve the ranking performance. As for the docking scoring functions, a force field-based scoring function is more suitable for ranking congeneric compounds.

| S-EPMC6072424 | biostudies-literature

Insight into the LFA-1/SARS-CoV-2 Orf7a Complex by Protein-Protein Docking, Molecular Dynamics, and MM-GBSA Calculations.

Project description:In the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) genome, open reading frames (ORFs) encode for viral accessory proteins. Among these, Orf7a structurally resembles the members of the immunoglobulin (Ig) superfamily and intracellular adhesion molecules (ICAMs), in particular. ICAMs are involved in integrin binding through lymphocyte function-associated antigen 1 (LFA-1). Based on such considerations and on previous findings on SARS-CoV, it has been postulated that the formation of the LFA-1/Orf7a complex could contribute to SARS-CoV-2 infectivity and pathogenicity. With the current work, we aim at providing insight into this macromolecular assembly, taking advantage of the recently reported SARS-CoV-2 Orf7a structure. Protein-protein docking, molecular dynamics (MD) simulations, and a Molecular Mechanical-Generalized Born Surface Area (MM-GBSA)-based stage were enrolled to provide refined models.

| S-EPMC8189040 | biostudies-literature

Targeting low-druggability bromodomains: fragment based screening and inhibitor design against the BAZ2B bromodomain.

Project description:Bromodomains are epigenetic reader domains that have recently become popular targets. In contrast to BET bromodomains, which have proven druggable, bromodomains from other regions of the phylogenetic tree have shallower pockets. We describe successful targeting of the challenging BAZ2B bromodomain using biophysical fragment screening and structure-based optimization of high ligand-efficiency fragments into a novel series of low-micromolar inhibitors. Our results provide attractive leads for development of BAZ2B chemical probes and indicate the whole family may be tractable.

| S-EPMC3905694 | biostudies-literature

Comparison of radii sets, entropy, QM methods, and sampling on MM-PBSA, MM-GBSA, and QM/MM-GBSA ligand binding energies of F. tularensis enoyl-ACP reductase (FabI).

Project description:To validate a method for predicting the binding affinities of FabI inhibitors, three implicit solvent methods, MM-PBSA, MM-GBSA, and QM/MM-GBSA were carefully compared using 16 benzimidazole inhibitors in complex with Francisella tularensis FabI. The data suggests that the prediction results are sensitive to radii sets, GB methods, QM Hamiltonians, sampling protocols, and simulation length, if only one simulation trajectory is used for each ligand. In this case, QM/MM-GBSA using 6 ns MD simulation trajectories together with GB(neck2) , PM3, and the mbondi2 radii set, generate the closest agreement with experimental values (r(2) = 0.88). However, if the three implicit solvent methods are averaged from six 1 ns MD simulations for each ligand (called "multiple independent sampling"), the prediction results are relatively insensitive to all the tested parameters. Moreover, MM/GBSA together with GB(HCT) and mbondi, using 600 frames extracted evenly from six 0.25 ns MD simulations, can also provide accurate prediction to experimental values (r(2) = 0.84). Therefore, the multiple independent sampling method can be more efficient than a single, long simulation method. Since future scaffold expansions may significantly change the benzimidazole's physiochemical properties (charges, etc.) and possibly binding modes, which may affect the sensitivities of various parameters, the relatively insensitive "multiple independent sampling method" may avoid the need of an entirely new validation study. Moreover, due to large fluctuating entropy values, (QM/)MM-P(G)BSA were limited to inhibitors' relative affinity prediction, but not the absolute affinity. The developed protocol will support an ongoing benzimidazole lead optimization program.

| S-EPMC4688044 | biostudies-literature

Modeling Peptide Nucleic Acid Binding Enthalpies Using MM-GBSA.

Project description:The binding enthalpies of peptide nucleic acid (PNA) homoduplexes were predicted using a molecular mechanics generalized Born surface area approach. Using the nucleic acid nearest-neighbor model, these were decomposed into sequence parameters which could replicate the enthalpies from thermal melting experiments with a mean error of 8.7%. These results present the first systematic computational investigation into the relationship between sequence and binding energy for PNA homoduplexes and identified a stabilizing helix initiation enthalpy not observed for nucleic acids with phosphoribose backbones.

| S-EPMC9706560 | biostudies-literature

QM/MM calculations with deMon2k.

Project description:The density functional code deMon2k employs a fitted density throughout (Auxiliary Density Functional Theory), which offers a great speed advantage without sacrificing necessary accuracy. Powerful Quantum Mechanical/Molecular Mechanical (QM/MM) approaches are reviewed. Following an overview of the basic features of deMon2k that make it efficient while retaining accuracy, three QM/MM implementations are compared and contrasted. In the first, deMon2k is interfaced with the CHARMM MM code (CHARMM-deMon2k); in the second MM is coded directly within the deMon2k software; and in the third the Chemistry in Ruby (Cuby) wrapper is used to drive the calculations. Cuby is also used in the context of constrained-DFT/MM calculations. Each of these implementations is described briefly; pros and cons are discussed and a few recent applications are described briefly. Applications include solvated ions and biomolecules, polyglutamine peptides important in polyQ neurodegenerative diseases, copper monooxygenases and ultra-rapid electron transfer in cryptochromes.

| S-EPMC6272552 | biostudies-literature

OmicsDI is part of the ELIXIR infrastructure

OmicsDI is an Elixir interoperability service. Learn more ›

Tweets

OmicsDI Databases

PRIDE
PeptideAtlas
MassIVE
JPOST Repository
Physiome Model Repository

EGA
EVA
ENA
LINCS
PAXDB
Cell Collective

MetaboLights
Metabolomics Workbench
MetabolomeExpress
GNPS
BioModels
FAIRDOMHub

ArrayExpress
dbGaP
ExpressionAtlas
GEO
NODE

Information

Databases
Help
API
Contact us
Code on GitHub
Terms of use
Submit Data