CoMFA, CoMSIA, molecular docking and MOLCAD studies of pyrimidinone derivatives to design novel and selective tankyrase inhibitors

Role of Thylakoid Lipids in Protochlorophyllide Oxidoreductase Activation: Allosteric Mechanism Elucidated by a Computational Study

Light-dependent protochlorophyllide oxidoreductase (LPOR) is a chlorophyll synthetase that catalyzes the reduction of protochlorophyllide (Pchlide) to chlorophyllide (Chlide) with indispensable roles in regulating photosynthesis processes. A recent study confirmed that thylakoid lipids (TL) were able to allosterically enhance modulator-induced LPOR activation. However, the allosteric modulation mechanism of LPOR by these compounds remains unclear. Herein, we integrated multiple computational approaches to explore the potential cavities in the Arabidopsis thaliana LPOR and an allosteric site around the helix-G region where high affinity for phosphatidyl glycerol (PG) was identified. Adopting accelerated molecular dynamics simulation for different LPOR states, we rigorously analyzed binary LPOR/PG and ternary LPOR/NADPH/PG complexes in terms of their dynamics, energetics, and attainable allosteric regulation. Our findings clarify the experimental observation of increased NADPH binding affinity for LPOR with PGs. Moreover, the simulations indicated that allosteric regulators targeting LPOR favor a mechanism involving lid opening upon binding to an allosteric hinge pocket mechanism. This understanding paves the way for designing novel LPOR activators and expanding the applications of LPOR.

Unravelling the Structural Mechanism of Action of 5-methyl-5-[4-(4-oxo-3H-quinazolin-2-yl)phenyl]imidazolidine-2,4-dione in Dual-Targeting Tankyrase 1 and 2: A Novel Avenue in Cancer Therapy

Tankyrase (TNKS) belonging to the poly(ADPribose) polymerase family, are known for their multi-functioning capabilities, and play an essential role in the Wnt β-catenin pathway and various other cellular processes. Although showing inhibitory potential at a nanomolar level, the structural dual-inhibitory mechanism of the novel TNKS inhibitor, 5-methyl-5-[4-(4-oxo-3H-quinazolin-2-yl)phenyl]imidazolidine-2,4-dione, remains unexplored. By employing advanced molecular modeling, this study provides these insights. Results of sequence alignments of binding site residues identified conserved residues; GLY1185 and ILE1224 in TNKS-1 and PHE1035 and PRO1034 in TNKS-2 as crucial mediators of the dual binding mechanism of 5-methyl-5-[4-(4-oxo-3H-quinazolin-2-yl)phenyl]imidazolidine-2,4-dione, corroborated by high per-residue energy contributions and consistent high-affinity interactions of these residues. Estimation of the binding free energy of 5-methyl-5-[4-(4-oxo-3H-quinazolin-2-yl)phenyl]imidazolidine-2,4-dione showed estimated total energy of -43.88 kcal/mol and -30.79 kcal/mol towards TNKS-1 and 2, respectively, indicating favorable analogous dual binding as previously reported. Assessment of the conformational dynamics of TNKS-1 and 2 upon the binding of 5-methyl-5-[4-(4-oxo-3H-quinazolin-2-yl)phenyl]imidazolidine-2,4-dione revealed similar structural changes characterized by increased flexibility and solvent assessible surface area of the residues inferring an analogous structural binding mechanism. Insights from this study show that peculiar, conserved residues are the driving force behind the dual inhibitory mechanism of 5-methyl-5-[4-(4-oxo-3H-quinazolin-2-yl)phenyl]imidazolidine-2,4-dione and could aid in the design of novel dual inhibitors of TNKS-1 and 2 with improved therapeutic properties.

Developing a variation of 3D-QSAR/MD method in drug design

In continuation of the previous reports on a combination of 3D-quantitative structure-activity relationships (QSAR) with computational molecular dynamics (MD) studies, a new variation of 3D-QSAR/MD method has been employed for drug-design as an alternative or supplementary for the typical experimental methods. The presented method is more cost-effective and less time-consuming than the previous methods and avoids several restrictions of experimental methods, such as validity estimation, and predictability. For this purpose, seven inhibitors for bromodomain (BRD)-containing protein, as an important protein in the development of different types of cancer and responsible for oncogenic rearrangements, have been selected to study of their interactions by docking and MD simulations using molecular mechanics/generalized born surface area (MM/GBSA) method. To build the proposed model, a common variant of 3D-QSAR methods, comparative molecular field analysis has been employed using a dataset of 100 MD-extracted ligand conformations and their corresponding MM/GBSA BRD4-binding energies. The results showed excellent predictability of the generated model for both the training set and test groups. Finally, two new inhibitors were selected among total 4000 designed derivatives (generated through evolutionary techniques) using the proposed 3D-QSAR-MD model. The potentials of these inhibitors were assessed by MD simulations, which showed the higher inhibitory of these compounds than the previous inhibitors. Therefore, this method showed high potentials for acceleration of the procedure of drug design and a basis for joining researchers in computational biology and pharmaceutical sciences.