Computational investigation of the binding mode of bis(hydroxylphenyl)arenes in 17β-HSD1: molecular dynamics simulations, MM-PBSA free energy calculations, and molecular electrostatic potential maps

Anticancer potential of β-sitosterol and oleanolic acid as through inhibition of human estrogenic 17beta-hydroxysteroid dehydrogenase type-1 based on an in silico approach

The human estrogenic enzyme 17beta-hydroxysteroid dehydrogenase type-1 (HSD17B1) provides biosynthesis regulation of active estrogen in stimulating the development of breast cancer through cell proliferation. The β-sitosterol is classified as a steroid compound and is actually a type of triterpenoid compound that has a similar structure to a steroid. This similarity provides a great opportunity for the inhibitor candidate to bind to the HDS17B1 enzyme because of the template similarity on the active site. Several in silico approaches have been applied in this study to examine the potential of these two inhibitor candidates. Pharmacokinetic studies showed positive results by meeting several drug candidate criteria, such as drug-likeness, bioavailability, and ADMET properties. A combination of molecular docking and MD simulation showed good conformational interaction of the inhibitors and HSD17B1. Prediction of binding free energy (ΔG bind) using the Molecular Mechanics-Generalized Born Surface Area (MM-GBSA) approach shows ΔG bind (kcal mol-1) of C1-HSD17B1: -49.31 ± 0.23 and C2-HSD17B1: -33.54 ± 0.34. Meanwhile, decomposition energy analysis (ΔG residue bind) suggested several key residues that were also responsible for the interaction with inhibitors, such as C1-HSD17B1 (six residues: Leu96, Leu149, Pro187, Met193, Val225, and Phe226) and C2-HSD17B1 (four residues: Ile14, Gly94, Pro187, and Val188). Hopefully, the obtained results from this research could be considered for the mechanistic inhibition of the HSDS17B1 enzyme at molecular and atomistic levels.

Identification of therapeutic drugs against COVID-19 through computational investigation on drug repurposing and structural modification

Global prevalence of coronavirus disease 2019 (COVID-19) calls for an urgent development of anti-viral regime. Compared with the development of new drugs, drug repurposing can significantly reduce the cost, time, and safety risks. Given the fact that coronavirus harnesses spike protein to invade host cells through angiotensin-converting enzyme 2 (ACE2), hence we see if any previous anti-virtual compounds can block spike-ACE2 interaction and inhibit the virus entry. The results of molecular docking and molecular dynamic simulations revealed that remdesivir exhibits better than expected anti-viral invasion potential against COVID-19 among the three types of compounds including remdesivir, tenofovir and lopinavir. In addition, a positive correlation between the surface area occupied by remdesivir and anti-viral invasion potential was also found. As such, the structure of remdesivir was modified by linking an N-benzyl substituted diamidine derivative to its hydroxyl group through an ester bond. It was found that this compound has a higher anti-viral invasion potential and greater specificity.

The investigation of the G-quadruplex aptamer selectivity to Pb2+ ion: a joint molecular dynamics simulation and density functional theory study

The aptamers with the ability to form a G-quadruplex structure can be stable in the presence of some ions. Hence, study of the interactions between such aptamers and ions can be beneficial to determine the highest selective aptamer toward an ion. In this article, molecular dynamics (MD) simulations and quantum mechanics (QM) calculations have been applied to investigate the selectivity of the T30695 aptamer toward Pb²⁺ in comparison with some ions. The Free Energy Landscape (FEL) analysis indicates that Pb²⁺ has remained inside the aptamer during the MD simulation, while the other ions have left it. The Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) binding energies prove that the conformational stability of the aptamer is the highest in the presence of Pb²⁺. According to the compaction parameters, the greatest compressed ion-aptamer complex, and hence, the highest ion-aptamer interaction have been induced in the presence of Pb²⁺. The contact maps clarify the closer contacts between the nucleotides of the aptamer in the presence of Pb²⁺. The density functional theory (DFT) results show that Pb²⁺ forms the most stable complex with the aptamer, which is consistent with the MD results. The QM calculations reveal that the N-H bonds and the O…H distances are the longest and the shortest, respectively, in the presence of Pb²⁺. The obtained results verify that the strongest hydrogen bonds (HBs), and hence, the most compressed aptamer structure are induced by Pb²⁺. Besides, atoms in molecules (AIM) and natural bond orbital (NBO) analyses confirm the results.Communicated by Ramaswamy H. Sarma.

Investigation of the binding mode of a novel cruzain inhibitor by docking, molecular dynamics, ab initio and MM/PBSA calculations

Chagas disease remains a major health problem in South America, and throughout the world. The two drugs clinically available for its treatment have limited efficacy and cause serious adverse effects. Cruzain is an established therapeutic target of Trypanosoma cruzi, the protozoan that causes Chagas disease. Our group recently identified a competitive cruzain inhibitor (compound 1) with an IC₅₀ = 15 µM that is also more synthetically accessible than the previously reported lead, compound 2. Prior studies, however, did not propose a binding mode for compound 1, hindering understanding of the structure-activity relationship and optimization. Here, the cruzain binding mode of compound 1 was investigated using docking, molecular dynamics (MD) simulations with ab initio derived parameters, ab initio calculations, and MM/PBSA. Two ligand protonation states and four binding poses were evaluated. A careful ligand parameterization method was employed to derive more physically meaningful parameters than those obtained by automated tools. The poses of unprotonated 1 were unstable in MD, showing large conformational changes and diffusing away from the binding site, whereas the protonated form showed higher stability and interaction with negatively charged residues Asp161 and Cys25. MM/PBSA also suggested that these two residues contribute favorably to binding of compound 1. By combining results from MD, ab initio calculations, and MM/PBSA, a binding mode of 1 is proposed. The results also provide insights for further optimization of 1, an interesting lead compound for the development of new cruzain inhibitors.

Exploration of the binding affinities between ecdysone agonists and EcR/USP by docking and MM-PB/GBSA approaches

Ecdysone receptor (EcR) is a significant target in the identification of new environmentally friendly pesticides. There are two types of ecdysone agonists: steroidal ecdysone agonists and dibenzoylhydrazines (DBHs). In this study, various modeling methods (homology modeling, molecular docking, MD simulation, binding free energy calculation, and per-residue binding free energy decomposition) were utilized to study the different binding mechanisms of two types of ecdysone agonists. Our theoretical results indicated that the relative binding potencies of DBHs can be ranked sufficiently accurately using the MOE docking method. However, MM/PBSA calculations more accurately predicted the binding affinities between steroidal ecdysone agonists and EcR-LBD. To identify the key residues involved in ecdysone agonist binding, the binding free energy (ΔG _Bind) was decomposed into the energy contributions of individual residues. The results revealed that nine residues-Ile339, Thr343, Met380, Met381, Tyr403, Tyr408, Asp419, Gln503, and Asn504-determined the binding affinities of the DBHs. Glu309, Met342, Arg383, Arg387, and Leu396 were important influences on the binding affinities of the steroidal ecdysone agonists. Graphical abstract The ecdysone receptor (EcR) is related to insect growth and has been shown to be a useful target for insecticides.

Inhibition of acetylcholinesterase by two genistein derivatives: kinetic analysis, molecular docking and molecular dynamics simulation

In this study two genistein derivatives (G1 and G2) are reported as inhibitors of acetylcholinesterase (AChE) and butyrylcholinesterase (BuChE), and differences in the inhibition of AChE are described. Although they differ in structure by a single methyl group, the inhibitory effect of G1 (IC₅₀=264 nmol/L) on AChE was 80 times stronger than that of G2 (IC₅₀=21,210 nmol/L). Enzyme-kinetic analysis, molecular docking and molecular dynamics (MD) simulations were conducted to better understand the molecular basis for this difference. The results obtained by kinetic analysis demonstrated that G1 can interact with both the catalytic active site and peripheral anionic site of AChE. The predicted binding free energies of two complexes calculated by the molecular mechanics/generalized born surface area (MM/GBSA) method were consistent with the experimental data. The analysis of the individual energy terms suggested that a difference between the net electrostatic contributions (ΔE_ele+ΔG_GB) was responsible for the binding affinities of these two inhibitors. Additionally, analysis of the molecular mechanics and MM/GBSA free energy decomposition revealed that the difference between G1 and G2 originated from interactions with Tyr124, Glu292, Val294 and Phe338 of AChE. In conclusion, the results reveal significant differences at the molecular level in the mechanism of inhibition of AChE by these structurally related compounds.

Molecular basis of HHQ biosynthesis: molecular dynamics simulations, enzyme kinetic and surface plasmon resonance studies

BACKGROUND

PQS (PseudomonasQuinolone Signal) and its precursor HHQ are signal molecules of the P. aeruginosa quorum sensing system. They explicate their role in mammalian pathogenicity by binding to the receptor PqsR that induces virulence factor production and biofilm formation. The enzyme PqsD catalyses the biosynthesis of HHQ.

RESULTS

Enzyme kinetic analysis and surface plasmon resonance (SPR) biosensor experiments were used to determine mechanism and substrate order of the biosynthesis. Comparative analysis led to the identification of domains involved in functionality of PqsD. A kinetic cycle was set up and molecular dynamics (MD) simulations were used to study the molecular bases of the kinetics of PqsD. Trajectory analysis, pocket volume measurements, binding energy estimations and decompositions ensured insights into the binding mode of the substrates anthraniloyl-CoA and β-ketodecanoic acid.

CONCLUSIONS

Enzyme kinetics and SPR experiments hint at a ping-pong mechanism for PqsD with ACoA as first substrate. Trajectory analysis of different PqsD complexes evidenced ligand-dependent induced-fit motions affecting the modified ACoA funnel access to the exposure of a secondary channel. A tunnel-network is formed in which Ser317 plays an important role by binding to both substrates. Mutagenesis experiments resulting in the inactive S317F mutant confirmed the importance of this residue. Two binding modes for β-ketodecanoic acid were identified with distinct catalytic mechanism preferences.

Exploring the PDE5 H-pocket by ensemble docking and structure-based design and synthesis of novel β-carboline derivatives

By studying the co-crystal information of interactions between PDE5 and its inhibitors, forty new tetrahydro-β-carbolines based-analogues were synthesized, and tested for their PDE5 inhibition. Some compounds were as active as tadalafil in inhibiting PDE5 and of better selectivity profile particularly versus PDE11A, the nature of the terminal ring and its nitrogen substituent are the main determinants of selectivity. Ensemble docking confirmed the role of H-loop closed conformer in activity versus its occluded and open forms. Conformational studies showed the effect of bulkiness of the terminal ring N-alkyl substituent on the formation of stable enzyme ligands conformers. The difference in potencies of hydantoin and piperazinedione analogues, together with the necessity of C-5/C-6 R-absolute configuration has been revealed through molecular docking.

Structural basis for species specific inhibition of 17β-hydroxysteroid dehydrogenase type 1 (17β-HSD1): computational study and biological validation

17β-Hydroxysteroid dehydrogenase type 1 (17β-HSD1) catalyzes the reduction of estrone to estradiol, which is the most potent estrogen in humans. Inhibition of 17β-HSD1 and thereby reducing the intracellular estradiol concentration is thus a promising approach for the treatment of estrogen dependent diseases. In the past, several steroidal and non-steroidal inhibitors of 17β-HSD1 have been described but so far there is no cocrystal structure of the latter in complex with 17β-HSD1. However, a distinct knowledge of active site topologies and protein-ligand interactions is a prerequisite for structure-based drug design and optimization. An elegant strategy to enhance this knowledge is to compare inhibition values obtained for one compound toward ortholog proteins from various species, which are highly conserved in sequence and differ only in few residues. In this study the inhibitory potencies of selected members of different non-steroidal inhibitor classes toward marmoset 17β-HSD1 were determined and the data were compared with the values obtained for the human enzyme. A species specific inhibition profile was observed in the class of the (hydroxyphenyl)naphthols. Using a combination of computational methods, including homology modelling, molecular docking, MD simulation, and binding energy calculation, a reasonable model of the three-dimensional structure of marmoset 17β-HSD1 was developed and inhibition data were rationalized on the structural basis. In marmoset 17β-HSD1, residues 190 to 196 form a small α-helix, which induces conformational changes compared to the human enzyme. The docking poses suggest these conformational changes as determinants for species specificity and energy decomposition analysis highlighted the outstanding role of Asn152 as interaction partner for inhibitor binding. In summary, this strategy of comparing the biological activities of inhibitors toward highly conserved ortholog proteins might be an alternative to laborious x-ray or site-directed mutagenesis experiments in certain cases. Additionally, it facilitates inhibitor design and optimization by offering new information on protein-ligand interactions.