Reduced state transition barrier of CDK6 from open to closed state induced by Thr177 phosphorylation and its implication in binding modes of inhibitors

Unveiling Conformational States of CDK6 Caused by Binding of Vcyclin Protein and Inhibitor by Combining Gaussian Accelerated Molecular Dynamics and Deep Learning

CDK6 plays a key role in the regulation of the cell cycle and is considered a crucial target for cancer therapy. In this work, conformational transitions of CDK6 were identified by using Gaussian accelerated molecular dynamics (GaMD), deep learning (DL), and free energy landscapes (FELs). DL finds that the binding pocket as well as the T-loop binding to the Vcyclin protein are involved in obvious differences of conformation contacts. This result suggests that the binding pocket of inhibitors (LQQ and AP9) and the binding interface of CDK6 to the Vcyclin protein play a key role in the function of CDK6. The analyses of FELs reveal that the binding pocket and the T-loop of CDK6 have disordered states. The results from principal component analysis (PCA) indicate that the binding of the Vcyclin protein affects the fluctuation behavior of the T-loop in CDK6. Our QM/MM-GBSA calculations suggest that the binding ability of LQQ to CDK6 is stronger than AP9 with or without the binding of the Vcyclin protein. Interaction networks of inhibitors with CDK6 were analyzed and the results reveal that LQQ contributes more hydrogen binding interactions (HBIs) and hot interaction spots with CDK6. In addition, the binding pocket endures flexibility changes from opening to closing states and the Vcyclin protein plays an important role in the stabilizing conformation of the T-loop. We anticipate that this work could provide useful information for further understanding the function of CDK6 and developing new promising inhibitors targeting CDK6.

Roles of Accelerated Molecular Dynamics Simulations in Predictions of Binding Kinetic Parameters

Rational predictions on binding kinetics parameters of drugs to targets play significant roles in future drug designs. Full conformational samplings of targets are requisite for accurate predictions of binding kinetic parameters. In this review, we mainly focus on the applications of enhanced sampling technologies in calculations of binding kinetics parameters and residence time of drugs. The methods involved in molecular dynamics simulations are applied to not only probe conformational changes of targets but also reveal calculations of residence time that is significant for drug efficiency. For this review, special attention are paid to accelerated molecular dynamics (aMD) and Gaussian aMD (GaMD) simulations that have been adopted to predict the association or disassociation rate constant. We also expect that this review can provide useful information for future drug design.

Exploring the structural-activity relationship of hexahydropyrazino[1,2-d]pyrido[3,2-b][1,4]oxazine derivatives as potent and orally-bioavailable PARP7 inhibitors

PARP7 has emerged as a promising anti-tumor target due to its crucial roles in nucleic acid sensing and immune regulation. Herein, we explored the structural-activity relationship of tricyclic PARP7 inhibitors containing a hexahydropyrazino[1,2-d]pyrido[3,2-b][1,4]oxazine motif. The effects of the chirality of the fused rings, the group conjugated to the fused rings, and the size of the linker on PARP7 inhibition were fully investigated. Our work leads to the discovery of an extremely potent and orally-bioavailable PARP7 inhibitor, namely 18 (PARP7 inhibition IC50 = 0.56 nM), for efficacious treatment of lung cancer in vivo. Notably, 18 showed acceptable bioavailability in ICR mice (F = 33.9%) and Beagle dogs (F = 45.2%). Further investigation of ADME-T properties suggested that 18 has the potential to be developed as a candidate drug molecule for PARP7-sensitive tumors.

Inclusion complex of 20(S)-protopanaxatriol with modified β-cyclodextrin: Characterization, solubility, and interaction with bovine serum albumin

20(S)-protopanaxatriol (PPT) is one of the ginsenosides isolated from Panax ginseng which have many pharmaceutical activities. However, the poor water solubility of PPT restrict its applications. Herein, a novel bridged-bis-[6-(3,3'-(ethylenedioxy) bis (propylamine))-6-deoxy-β-cyclodextrin] (EDBA-bis-β-CD) was designed and synthesized, and the inclusion complex (IC) of EDBA-bis-β-CD with PPT was successfully prepared in the solid state, and characterized by UV, ¹H NMR, 2D ROESY, FT-IR, XRD and SEM and molecular modelling methods. The continuous variation method analysis indicated that the stoichiometry of the IC was 1:1. UV-vis spectral analysis demonstrated the binding constant K_s was 995.94 M^-1, and the solubility study showed that the solubility of PPT improved 290 times. The interaction of the IC with bovine serum albumin (BSA) was investigated via fluorescence spectroscopy. The results indicated that fluorescence quenching of BSA by IC was static quenching. Thermodynamic studies showed that van der Waals forces and hydrogen bonding play significant roles in interaction. The esterase-like activity of BSA in the presence of IC showed that it reduce the esterase activity of BSA in a competitive manner. Furthermore, molecular docking and molecular dynamics simulations for EDBA-bis-β-CD/PPT and BSA/IC systems were generated to provide information on the stability and the forces in the binding.

Unveiling conformational dynamics changes of H-Ras induced by mutations based on accelerated molecular dynamics

Uncovering molecular basis with regard to the conformational change of two switches I and II in the GppNHp (GNP)-bound H-Ras is highly significant for the understanding of Ras signaling. For this purpose, accelerated molecular dynamics (aMD) simulations and principal component (PC) analysis are integrated to probe the effect of mutations G12V, T35S and Q61K on conformational transformation between two switches of the GNP-bound H-Ras. The RMSF and cross-correlation analyses suggest that three mutations exert a vital effect on the flexibility and internal dynamics of two switches in the GNP-bound H-Ras. The results stemming from PC analysis indicate that two switches in the GNP-bound WT H-Ras tend to form a closed state in most conformations, while those in the GNP-bound mutated H-Ras display transformation between different states. This conclusion is further supported by free energy landscapes constructed by using the distances of residues 12 away from 35 and 35 away from 61 as reaction coordinates and different experimental studies. Interaction scanning is performed on aMD trajectories and the information shows that conformational transformations of two switches I and II induced by mutations extremely affect the GNP-residue interactions. Meanwhile, the scanning results also signify that residues G15, A18, F28, K117, A146 and K147 form stable contacts with GNP, while residues D30, E31, Y32, D33, P34 and E62 in two switches I and II produce unstable contacts with GNP. This study not only reveals dynamic behavior changes of two switches in H-Ras induced by mutations, but also unveils general principles and mechanisms with regard to functional conformational changes of H-Ras.

Synthesis of a novel p-hydroxycinnamic amide with anticancer capability and its interaction with human serum albumin

In the present study, a novel p-hydroxycinnamic amide (E)-3-(4-hydroxyphenyl)-N-(4-(N-(5-meth oxypyrimidin-2-yl)-sulfamoyl)phenyl)acrylamide (HMSP) was synthesized and confirmed. In vitro cytotoxic assays indicated that HMSP was able to inhibit the proliferation of various cancer cell lines. The interaction between HMSP and human serum albumin (HSA) was examined by fluorescence, UV-Vis and circular dichroism (CD) spectra, in addition to molecular simulation. The fluorescence and UV-Vis spectra data indicated that the binding of HMSP with HSA was a static process. According to the fluorescence quenching calculation, the corresponding thermodynamic parameters, bimolecular quenching rate constant and apparent quenching constants were calculated. Van der Walls forces and hydrogen bonds were vital in the binding of HMSP on HSA. The distances between HSA and its derivatives were obtained. Furthermore, competitive experiments and molecular modeling results suggested that the binding of the compound on HSA mainly occurred in site I (sub-domain IIA). Changes in HSA conformation were observed from synchronous fluorescence and CD spectra, which were further investigated by molecular dynamic simulations.