Exploration of the selective binding mechanism of GSK3β via molecular modeling and molecular dynamics simulation studies

Molecular dynamics simulations on interactions of five antibiotics with luciferase of Vibrio Qinghaiensis sp.-Q67

A large number of antibiotics have been used in the medical industry, agriculture, and animal husbandry industry in recent years. It may cause pollution to the aquatic environment and ultimately threaten to human health due to their prolonged exposure to the environment. We aim to study the toxicity mechanism of enrofloxacin (ENR), chlortetracycline hydrochloride (CTC), trimethoprim (TMP), chloramphenicol (CMP), and erythromycin (ETM) to luciferase of Vibrio Qinghaiensis sp.-Q67 (Q67) by using toxicity testing combined with molecular docking, molecular dynamics, and binding free energy analysis. The curve categories for ENR were different from the other four antibiotics, with ENR being J-type and the rest being S-type, and the toxicity of these five antibiotics (pEC50) followed the order of ENR (7.281) > ETM (6.814) > CMP (6.672) > CTC (6.400) > TMP (6.123), the order of toxicity value is consistent with the the magnitude of the binding free energy (ENR (-47.759 kcal/mol), ETM (-46.821 kcal/mol), CMP (-42.905 kcal/mol), CTC (-40.946 kcal/mol), TMP (-28.251 kcal/mol)). The van der Waals force provided the most important contribution to the binding free energy of the five antibiotics in the binding system with Q67 luciferase. Therefore, the dominant factor for the binding of antibiotics to luciferase was shape compensation. The face-to-face π-π stacking interaction between the diazohexane structure outside the active pocket region and the indoles structure of Phe194 and Phe250 in the molecular structure was the main reason for the highest toxicity value of antibiotic ENR. The hormesis effect of ENR has a competitive binding relationship with the α and β subunits of luciferase. Homology modeling, molecular docking, molecular dynamics simulations and binding free energy calculations were used to derive the toxicity magnitude of different antibiotics against Q67, and insights at the molecular level. The conclusion of toxicological experiments verified the correctness of the simulation results. This study contributes to the understanding of toxicity mechanisms of five antibiotics and facilitates risk assessment of antibiotic contaminants in the aquatic environment.

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

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).

In Silico Studies to Develop New GSK3β Inhibitors Effective in Alzheimer's Disease

Background:

Alzheimer's disease affects a large part of the world population by prolonging the human life span and becoming an economic burden in the health system. Therefore, its treatment becomes more and more important every day. With the insufficiency of existing drug molecules, new drug targets started to be emphasized. The most important of these is the Glycogen Synthase Kinase 3β enzyme, thought to be of key importance in Tau hyperphosphorylation and Amyloid β accumulation mechanisms.

Objective:

In this research, computational studies were conducted to develop a new GSK3β enzyme inhibitor.

Method:

Leading compounds suitable for pharmacophore models obtained by the 3D QSAR method were scanned in databases. In silico ADME/Tox analyses were performed on the obtained molecules.

Results:

Although the three molecules (ENA99104, CNR13756, TIM405938) had strong Dock Scores (42.869, 53.344, and 41.119, respectively) in molecular docking calculations, only the CNR13756 molecule was found successful according to molecular dynamics simulations.

Conclusion:

All computational studies have revealed that the CNR13756 molecule can exhibit a therapeutic scaffold property, thus obtaining a selective GSK3β inhibitor with minimal side effects.