Discovery of Novel Pyridin-2-yl Urea Inhibitors Targeting ASK1 Kinase and Its Binding Mode by Absolute Protein–Ligand Binding Free Energy Calculations


4.3. Binding Free Energy (BFE) Calculation by Molecular Dynamics

Molecular dynamics (MD) simulation was performed using the Particle Mesh Ewald (PME) module of AMBER software (v.2020) [48]. Based on the docking poses of ligand with ASK1 catalytic domain, the complex was solvated in a cubic water box, maintaining a minimal distance of 8.0 Å from system boundaries. Necessary counter ions (Na+ or Cl ions) were added to neutralize the system for simulation. The ff14SB force field and TIP3P water model were applied to the protein and water molecules, respectively, while the GAFF2 force field and RESP charge for the ligand were assigned based on geometry optimization at the HF/6-31G(d) level of theory using Gaussain16 [49].
The stability of complex was initially examined via standard MD simulation, as described in previous papers [25,50]. Briefly, the complex was optimized in 1000 cycles with half of steepest descent and half of conjugate gradient minimization. Subsequently, the complex box was gradually heated to 298 K over 200 ps using an NVT ensemble simulation (constant number of atoms, volume and temperature) and changed to NPT ensemble (constant number of atoms, pressure and temperature) to optimize the solvent environment for 100 ps simulation. During the NVT and NPT processes, the complex was restrained with 0.5 kcal/mol/Å2. After that, 500 ps NPT equilibrium simulation was performed without any constraint. Finally, a 5 ns production run was conducted to obtain an equilibrated complex.
For each complex, three parallel independent production runs were performed. The root-mean-square deviations (RMSD) of the heavy atoms in ligand were calculated. If the RMSD exceeded the threshold value of 2.0 Å, the ligand was deemed unstable in binding pocket [51,52]. Otherwise, the final snapshot with the smallest RMSD of ligand was selected as the starting geometry for the subsequent binding free energy (BFE) calculation.
For the BFE calculation, a force constraint was introduced to enhance sampling [33]. Three reference atoms from the protein and ligand were selected based on their fluctuations during MD simulation. The detailed selection of reference atoms is described in previous studies [25]. In the BFE calculation, the contributions were divided into electrostatic (ELE) and van der Waals (vdW) components using a series of discrete windows (λ) based on thermodynamic integration (TI) method. This process was similar to the above general MD simulation, but with different coupling λ values. Firstly, the final snapshot from previous MD run was further optimized using 500 steepest descent cycles followed by 500 conjugate gradient steps. It should be noted that only 1000 steepest descent cycles were performed in the minimization of the vdW term. Next, a 200 ps heating procedure was conducted in NVT ensemble, followed by 500 ps equilibration in the NPT ensemble. Finally, a 10 ns production run was carried out. For the ELE calculation, 21 windows with equal intervals were used (λ = 0.00, 0.05, 0.10, 0.15, 0.20, 0.25, 0.30, 0.35, 0.40, 0.45, 0.50, 0.55, 0.60, 0.65, 0.70, 0.75, 0.80, 0.85, 0.90, 0.95, 1.00), whereas 26 windows were used for vdW calculation (λ = 0.00, 0.04, 0.08, 0.12, 0.16, 0.20, 0.24, 0.28, 0.32, 0.36, 0.40, 0.44, 0.48, 0.52, 0.56, 0.60, 0.64, 0.68, 0.72, 0.76, 0.80, 0.84, 0.88, 0.92, 0.96, 1.00). The final binding affinity was estimated using the TI method from the alchemlyb toolkit [53].

4.4. General Route of Chemical Synthesis

Compounds were synthesized based on our efficient and environmentally friendly protocol [23]. Initially, we selected the commercially available 6-(4-isopropyl-4H-1,2,4-triazol-3-yl) pyridin-2-amine (IPTPA) as the starting reactant. Subsequently, 1.0 equivalent of phenyl chloroformate was mixed and agitated overnight. Following the addition of aromatic/non-aromatic amines together with N,N-diisopropylethylamine, the desired pyridin-2-yl ureas were obtained with high yield via a concerted mechanism [23].
Aside from compounds 1, 6 and 7, which were synthesized previously [23], we successfully synthesized four new pyridin-2-yl urea compounds in this study.

1H and 13C spectra were collected on 500 MHz NMR spectrometer (Bruker AVANCE). Chemical shifts for protons are reported in parts per million (ppm) downfield and are referenced to residual protium in the NMR solvent (CDCl3 = δ 7.26). Chemical shifts for carbon are reported in parts per million downfield and are referenced to the carbon resonances of solvent (CDCl3 = δ 77.16). Data are represented as follows: chemical shift, multiplicity (br = broad, s = singlet, d = double, t = triplet, q = quartet, m = multiplet, heptet using full name), coupling constants in Hertz (Hz) and integration.

High-resolution mass spectra (HRMS) were collected using a Waters XEVO G2-S Q-TOF-MS mass spectrometer (Waters Corporation, Boston, MA, USA) in positive mode using MeCN/H2O; I Class liquid chromatograph; ACQUITY UPLC HSS T3 C18 (2.1 × 100 mm, 1.8 μm) chromatographic column; and MassLynxV4.2 analysis software.

Analytical thin layer chromatography (TLC) was performed on precoated silica gel GF254 HPTLC plates (5 × 10 cm2) purchased from Yantai Jiangyou Chemical Co., Ltd. (Yantai, China). The developed chromatogram was analyzed by UV lamp (254 and 365 nm). The non-UV active compounds were generally visualized through placing the plates in a sealed TLC tank containing iodine and silica gel (the 160–200 mesh mentioned above), if necessary, with the aid of a heating gun. The progress of reaction was also monitored by TLC stained with 5% v/v ethanol aqueous solution of concentrated H2SO4 and estimated in a concentration- and time-dependent manner.

Pyrrolidine (99%) was supplied by Shanghai Macklin Biochemical Technology Co., Ltd. (Shanghai, China) Phenyl chloroformate (98%) was purchased from Shanghai Aladdin Biochemical Technology Co., Ltd. (Shanghai, China) 6-(Trifluoromethoxy)indoline (98%), 6-methoxyindoline (98%), 5-(trifluoromethyl)indoline (97%), 1,2,3,4-tetrahydroquinoline (99%), 1,2,3,4-tetrahydroisoquinoline (98%), 6-(4-isopropyl-4H-1,2,4-triazol-3-yl)pyridin-2-amine (IPTPA, 95%) and N,N-diisopropylethylamine (DIPEA, 99%) were obtained from Shanghai Bidepharm. 5-Methoxyindoline (95%) was provided by Shanghai Leyan Co., Ltd., Shanghai, China.

4.5. Expression and Purification of ASK1

The catalytic domain of ASK1 (encoding residues 659–951) was modified to include an N-terminal 6XHis tag and cloned into expression vector, pMCSG7. The recombinant human ASK1 protein was expressed in E. coli Rosetta2 (DE3) in LB medium, containing 100 µg/mL ampicillin and 34 µg/mL chloramphenicol. A total of 2 mL overnight culture was diluted 100-fold into 200 mL LB and grown at 37 °C with shaking at 200 rpm until the OD600 reached 0.8. Subsequently, 0.5 mM IPTG was added, and protein expression was induced overnight at 20 °C.

The pellets were collected by centrifugation and resuspended in HEPES buffer (50 mM HEPES, 300 mM NaCl). After adding ddH2O, the suspension was placed on ice and subjected to ultrasonication (3 s pulses with 5 s intervals) for 15 min. The lysate was then separated by centrifugation at 18,000× g for 15 min, and the supernatant was purified using Ni-affinity chromatography columns with a wash buffer (50 mM HEPES, 300 mM NaCl and 25 mM imidazole) and an elution buffer of the same composition but containing 250 mM imidazole. The concentration of the purified protein was estimated using a BCA protein assay kit (Solarbio, Beijing, China).

4.6. ASK1 Enzyme Inhibition Assay

The IC50 value for each inhibitor was determined using the ADP-Glo™ protocol. Inhibitors were prepared as a series of 1:1 serial dilutions (final concentration of 10 μM) in 40 mM HEPES buffer (pH 7.5) containing 20 mM MgCl2, 0.1 mg/mL BSA, 50 μM DTT and 5% DMSO.

The kinase assays were conducted in a 5 μL volume containing the following final concentrations: 6.25 ng/μL active ASK1, 25 μM ATP and 0.1 μg/μL MBP. The reaction was terminated by adding 5 μL of ADP-Glo reagent and incubated at room temperature for 40 min. Then, 10 μL of kinase detection reagent was added, and the mixture was incubated for an additional 30 min. Luminescence was measured using a VICTOR Nivo microplate reader (PerkinElmer, Waltham, MA, USA), and the concentration–response curve and half-maximal inhibitory concentration (IC50) values were obtained by fitting the data using GraphPad Prism (v.9.0) software.



Source link

Lingzhi Wang www.mdpi.com