Browse
Submit Data
Databases
API
Help

Dataset Information

14 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Theoretical exploration of the binding selectivity of inhibitors to BRD7 and BRD9 with multiple short molecular dynamics simulations.

ABSTRACT: Bromodomain-containing proteins 7 and 9 (BRD7 and BRD9) have been considered as potential targets of clinical drug design toward treatment of human cancers and other diseases. Multiple short molecular dynamics simulations and binding free energy predictions were carried out to decipher the binding selectivity of three inhibitors 4L2, 5U6, and 6KT toward BRD7 and BRD9. The results show that 4L2 has more favorable binding ability to BRD7 over BRD9 compared to 5U6 and 6KT, while 5U6 and 6KT possess more favorable associations with BRD9 than BRD7. Furthermore, estimations of residue-based free energy decompositions further identify that four common residue pairs, including (F155, F44), (V160, V49), (Y168, Y57) and (Y217, Y106) in (BRD7, BRD9) generate obvious binding differences with 4L2, 5U6, and 6KT, which mostly drives the binding selectivity of 4L2, 5U6, and 6KT to BRD7 and BRD9. Dynamic information arising from trajectory analysis also suggests that inhibitor bindings affect structural flexibility and motion modes, which is responsible for the partial selectivity of 4L2, 5U6, and 6KT toward BRD7 and BRD9. As per our expectation, this study theoretically provides useful hints for design of dual inhibitors with high selectivity on BRD7 and BRD9.

SUBMITTER: Wang L

PROVIDER: S-EPMC9169554 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Similar Datasets

Deciphering Selectivity Mechanism of BRD9 and TAF1(2) toward Inhibitors Based on Multiple Short Molecular Dynamics Simulations and MM-GBSA Calculations.

Project description:BRD9 and TAF1(2) have been regarded as significant targets of drug design for clinically treating acute myeloid leukemia, malignancies, and inflammatory diseases. In this study, multiple short molecular dynamics simulations combined with the molecular mechanics generalized Born surface area method were employed to investigate the binding selectivity of three ligands, 67B, 67C, and 69G, to BRD9/TAF1(2) with IC50 values of 230/59 nM, 1400/46 nM, and 160/410 nM, respectively. The computed binding free energies from the MM-GBSA method displayed good correlations with that provided by the experimental data. The results indicate that the enthalpic contributions played a critical factor in the selectivity recognition of inhibitors toward BRD9 and TAF1(2), indicating that 67B and 67C could more favorably bind to TAF1(2) than BRD9, while 69G had better selectivity toward BRD9 over TAF1(2). In addition, the residue-based free energy decomposition approach was adopted to calculate the inhibitor-residue interaction spectrum, and the results determined the gatekeeper (Y106 in BRD9 and Y1589 in TAF1(2)) and lipophilic shelf (G43, F44, and F45 in BRD9 and W1526, P1527, and F1528 in TAF1(2)), which could be identified as hotspots for designing efficient selective inhibitors toward BRD9 and TAF1(2). This work is also expected to provide significant theoretical guidance and insightful molecular mechanisms for the rational designs of efficient selective inhibitors targeting BRD9 and TAF1(2).

| S-EPMC10052767 | biostudies-literature

Binding selectivity of inhibitors toward the first over the second bromodomain of BRD4: theoretical insights from free energy calculations and multiple short molecular dynamics simulations.

Project description:Bromodomain-containing protein 4 (BRD4) plays an important role in mediating gene transcription involved in cancers and non-cancer diseases such as acute heart failure and inflammatory diseases. In this work, multiple short molecular dynamics (MSMD) simulations are integrated with a molecular mechanics generalized Born surface area (MM-GBSA) approach to decipher binding selectivity of three inhibitors 8NS, 82Y, and 837 toward two domains BD1 and BD2 of BRD4. The results demonstrate that the enthalpy effects play critical roles in selectivity identification of inhibitors toward BD1 and BD2, determining that 8NS has better selectivity toward BD2 than BD1, while 82Y and 837 more favorably bind to BD1 than BD2. A residue-based free-energy decomposition method was used to calculate an inhibitor-residue interaction spectrum and unveil contributions of separate residues to binding selectivity. The results identify six common residues, containing (P82, P375), (V87, V380), (L92, L385), (L94, L387), (N140, N433), and (I146, V439) individually belonging to (BD1, BD2) of BRD4, and yield a considerable binding difference of inhibitors to BD1 and BD2, suggesting that these residues play key roles in binding selectivity of inhibitors toward BD1 and BD2 of BRD4. Therefore, these results provide useful dynamics information and a structure affinity relationship for the development of highly selective inhibitors targeting BD1 and BD2 of BRD4.

| S-EPMC8693360 | biostudies-literature

Exploration of Catalytic Selectivity for Aminotransferase (BtrR) Based on Multiple Molecular Dynamics Simulations.

Project description:The aminotransferase from Bacillus circulans (BtrR), which is involved in the biosynthesis of butirosin, catalyzes the pyridoxal phosphate (PLP)-dependent transamination reaction to convert valienone to ?-valienamine (a new ?-glycosidase inhibitor for the treatment of lysosomal storage diseases) with an optical purity enantiomeric excess value. To explore the stereoselective mechanism of valienamine generated by BtrR, multiple molecular dynamics (MD) simulations were performed for the BtrR/PLP/valienamine and BtrR/PLP/?-valienamine complexes. The theoretical results showed that ?-valienamine could make BtrR more stable and dense than valienamine. ?-valienamine could increase the hydrogen bond probability and decrease the binding free energy between coenzyme PLP and BtrR by regulating the protein structure of BtrR, which was conducive to the catalytic reaction. ?-valienamine maintained the formation of cation-p interactions between basic and aromatic amino acids in BtrR, thus enhancing its stability and catalytic activity. In addition, CAVER 3.0 analysis revealed that ?-valienamine could make the tunnel of BtrR wider and straight, which was propitious to the removal of products from BtrR. Steered MD simulation results showed that valienamine interacted with more residues in the tunnel during dissociation compared with ?-valienamine, resulting in the need for a stronger force to be acquired from BtrR. Taken together, BtrR was more inclined to catalyze the substrates to form ?-valienamine, either from the point of view of the catalytic reaction or product removal.

| S-EPMC6429434 | biostudies-literature

Theoretical Studies on Selectivity of HPK1/JAK1 Inhibitors by Molecular Dynamics Simulations and Free Energy Calculations.

Project description:Hematopoietic progenitor kinase 1 (HPK1) is a negative regulator of T cell receptor, which has been regarded as a potential target for immunotherapy. Yu et al. observed the off-target effect of the high-throughput screening HPK1 kinase inhibitor hits on JAK1 kinase. The off-target effect is usually due to the lack of specificity of the drug, resulting in toxic side effects. Therefore, exploring the mechanisms to selectively inhibit HPK1 is critical for developing effective and safe inhibitors. In this study, two indazole compounds as HPK1 inhibitors with different selectivity towards JAK1 were used to investigate the selectivity mechanism using multiple computational methods, including conventional molecular dynamics simulations, binding free energy calculations and umbrella sampling simulations. The results indicate that the salt bridge between the inhibitor and residue Asp101 of HPK1 favors their selectivity towards HPK1 over JAK1. Information obtained from this study can be used to discover and design more potent and selective HPK1 inhibitors for immunotherapy.

| S-EPMC9916865 | biostudies-literature

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

Project description: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.

| S-EPMC8154142 | biostudies-literature

Binding selectivity-dependent molecular mechanism of inhibitors towards CDK2 and CDK6 investigated by multiple short molecular dynamics and free energy landscapes.

Project description:Understanding selectivity-dependent molecular mechanism of inhibitors towards CDK2 over CDK6 is prominent for improving drug design towards the CDK family. Multiple short molecular dynamics (MD) simulations combined with MM-GBSA approach are adopted to investigate molecular mechanism on binding selectivity of inhibitors X64, X3A, and 4 AU to CDK2 and CDK6. The RMSF analysis and calculations of molecular surface areas indicate that local structural and global flexibility of CDK6 are stronger than that of CDK2. Based on dynamics cross-correlation maps (DCCMs), motion modes of CDK2 and CDK6 produce difference due to associations of X64, X3A, and 4 AU. The calculated binding free energies (BFEs) demonstrate that the compensation between binding enthalpy and entropy of X64, X34, and 4 AU is a key force driving selectivity of inhibitors towards CDK2 over CDK6. This work provides valuable information for designing highly selective inhibitors towards CDK2 and CDK6 and further promotes identification of efficient anticancer drugs in the future.

| S-EPMC9645286 | biostudies-literature

Exploration of virtual candidates for human HMG-CoA reductase inhibitors using pharmacophore modeling and molecular dynamics simulations.

Project description:3-hydroxy-3-methylglutaryl coenzyme A reductase (HMGR) is a rate-controlling enzyme in the mevalonate pathway which involved in biosynthesis of cholesterol and other isoprenoids. This enzyme catalyzes the conversion of HMG-CoA to mevalonate and is regarded as a drug target to treat hypercholesterolemia. In this study, ten qualitative pharmacophore models were generated based on chemical features in active inhibitors of HMGR. The generated models were validated using a test set. In a validation process, the best hypothesis was selected based on the statistical parameters and used for virtual screening of chemical databases to find novel lead candidates. The screened compounds were sorted by applying drug-like properties. The compounds that satisfied all drug-like properties were used for molecular docking study to identify their binding conformations at active site of HMGR. The final hit compounds were selected based on docking score and binding orientation. The HMGR structures in complex with the hit compounds were subjected to 10 ns molecular dynamics simulations to refine the binding orientation as well as to check the stability of the hits. After simulation, binding modes including hydrogen bonding patterns and molecular interactions with the active site residues were analyzed. In conclusion, four hit compounds with new structural scaffold were suggested as novel and potent HMGR inhibitors.

| S-EPMC3875450 | biostudies-literature

Molecular Dynamics Exploration of Selectivity of Dual Inhibitors 5M7, 65X, and 65Z toward Fatty Acid Binding Proteins 4 and 5.

Project description:Designing highly selective inhibitors of fatty acid binding proteins 4 and 5 (FABP4 and FABP5) is of importance for treatment of some diseases related with inflammation, metabolism, and tumor growth. In this study, molecular dynamics (MD) simulations combined with molecular mechanics generalized Born surface area (MM-GBSA) method were performed to probe binding selectivity of three inhibitors (5M7, 65X, and 65Z) to FABP4/FABP5 with Ki values of 0.022/0.50 μM, 0.011/0.086 μM, and 0.016/0.12 μM, respectively. The results not only suggest that all inhibitors associate more tightly with FABP4 than FABP5, but also prove that the main forces driving the selective bindings of inhibitors to FABP4 and FABP5 stem from the difference in the van der Waals interactions and polar interactions of inhibitors with two proteins. Meanwhile, a residue-based free energy decomposition method was applied to reveal molecular basis that inhibitors selectively interact with individual residues of two different proteins. The calculated results show that the binding difference of inhibitors to the residues (Phe16, Phe19), (Ala33, Gly36), (Phe57, Leu60), (Ala75, Ala78), (Arg126, Arg129), and (Tyr128, Tyr131) in (FABP4, FABP5) drive the selectivity of inhibitors toward FABP4 and FABP5. This study will provide great help for further design of effective drugs to protect against a series of metabolic diseases, arteriosclerosis, and inflammation.

| S-EPMC6164837 | biostudies-literature

Exploration of plant-derived natural polyphenols toward COVID-19 main protease inhibitors: DFT, molecular docking approach, and molecular dynamics simulations.

Project description:Recent outbreaks of coronavirus have brought serious challenges to public health around the world, and it is essential to find effective treatments. In this study, the 3C-like proteinase (3CLpro) of SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) has been considered as an important drug target because of its role in viral replication. We initially optimized 251 compounds at the PM7 level of theory for docking with 3CLpro, and then we selected the top 12 compounds for further optimization with the B3LYP-D3/6-311G** method and obtained the top four compounds by further molecular docking. Quantum chemistry calculations were performed to predict molecular properties, such as the electrostatic potential and some CDFT descriptors. We also performed molecular dynamics simulations and free energy calculations to determine the relative stability of the selected four potential compounds. We have identified key residues controlling the 3CLpro/ligand binding from per-residue based decomposition of the binding free energy. Convincingly, the comprehensive results support the conclusion that the compounds have the potential to become a candidate for anti-coronavirus treatment.

| S-EPMC8981245 | biostudies-literature

Dynamics of the EAG1 K<sup>+</sup> channel selectivity filter assessed by molecular dynamics simulations.

Project description:EAG1 channels belong to the KCNH family of voltage gated potassium channels. They are expressed in several brain regions and increased expression is linked to certain cancer types. Recent cryo-EM structure determination finally revealed the structure of these channels in atomic detail, allowing computational investigations. In this study, we performed molecular dynamics simulations to investigate the ion binding sites and the dynamical behavior of the selectivity filter. Our simulations suggest that sites S2 and S4 form stable ion binding sites, while ions placed at sites S1 and S3 rapidly switched to sites S2 and S4. Further, ions tended to dissociate away from S0 within less than 20 ns, due to increased filter flexibility. This was followed by water influx from the extracellular side, leading to a widening of the filter in this region, and likely non-conductive filter configurations. Simulations with the inactivation-enhancing mutant Y464A or Na+ ions lead to trapped water molecules behind the SF, suggesting that these simulations captured early conformational changes linked to C-type inactivation.

| S-EPMC6531291 | 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