Chelator-Based Parameterization of the 12-6-4 Lennard-Jones Molecular Mechanics Potential for More Realistic Metal Ion-Protein Interactions

Discovery of Potent Selective HDAC6 Inhibitors with 5-Phenyl-1H-indole Fragment: Virtual Screening, Rational Design, and Biological Evaluation

Among the HDACs family, histone deacetylase 6 (HDAC6) has attracted extensive attention due to its unique structure and biological functions. Numerous studies have shown that compared with broad-spectrum HDACs inhibitors, selective HDAC6 inhibitors exert ideal efficacy in tumor treatment with insignificant toxic and side effects, demonstrating promising clinical application prospect. Herein, we carried out rational drug design by integrating a deep learning model, molecular docking, and molecular dynamics simulation technology to construct a virtual screening process. The designed derivatives with 5-phenyl-1H-indole fragment as Cap showed desirable cytotoxicity to the various tumor cell lines, all of which were within 15 μM (ranging from 0.35 to 14.87 μM), among which compound 5i had the best antiproliferative activities against HL-60 (IC50 = 0.35 ± 0.07 μM) and arrested HL-60 cells in the G0/G1 phase. In addition, 5i exhibited better isotype selective inhibitory activities due to the potent potency against HDAC6 (IC50 = 5.16 ± 0.25 nM) and the reduced inhibitory activities against HDAC1 (selective index ≈ 124), which was further verified by immunoblotting results. Moreover, the representative binding conformation of 5i on HDAC6 was revealed and the key residues contributing 5i's binding were also identified via decomposition free-energy analysis. The discovery of lead compound 5i also indicates that virtual screening is still a beneficial tool in drug discovery and can provide more molecular skeletons with research potential for drug design, which is worthy of widespread application.

Modeling Lanthanide Ions in Solution: A Versatile Force Field in Aqueous and Organic Solvents

In this paper, we propose a new nonpolarizable force field for describing the Ln3+ (Ln = lanthanide) series based on a 12-6-4 Lennard-Jones potential. The development of the force field was performed in pure water by adjusting both the ion-oxygen distance and the hydration free energy. This force field accurately reproduces the Ln3+ hydration properties through the series, especially the coordination number that is hardly accessible using a nonpolarizable force field. Then, the validity and the transferability of the current force field were evaluated for two different systems containing Ln3+ in various solvents, namely, 0.1 mol L-1 La(NO3)3 salts in methanol and Eu(NO3)3 salts in solvent organic phases composed of DMDOHEMA molecules in n-heptane. The good agreement between our simulations and the data available in the literature confirms the accuracy of the force field for describing the lanthanide cations in both aqueous and nonaqueous media.

Comparison of force fields to study the zinc-finger containing protein NPL4, a target for disulfiram in cancer therapy

Molecular dynamics (MD) simulations are a powerful approach to studying the structure and dynamics of proteins related to health and disease. Advances in the MD field allow modeling proteins with high accuracy. However, modeling metal ions and their interactions with proteins is still challenging. NPL4 is a zinc-binding protein and works as a cofactor for p97 to regulate protein homeostasis. NPL4 is of biomedical importance and has been proposed as the target of disulfiram, a drug recently repurposed for cancer treatment. Experimental studies proposed that the disulfiram metabolites, bis-(diethyldithiocarbamate)‑copper and cupric ions, induce NPL4 misfolding and aggregation. However, the molecular details of their interactions with NPL4 and consequent structural effects are still elusive. Here, biomolecular simulations can help to shed light on the related structural details. To apply MD simulations to NPL4 and its interaction with copper the first important step is identifying a suitable force field to describe the protein in its zinc-bound states. We examined different sets of non-bonded parameters because we want to study the misfolding mechanism and cannot rule out that the zinc may detach from the protein during the process and copper replaces it. We investigated the force-field ability to model the coordination geometry of the metal ions by comparing the results from MD simulations with optimized geometries from quantum mechanics (QM) calculations using model systems of NPL4. Furthermore, we investigated the performance of a force field including bonded parameters to treat copper ions in NPL4 that we obtained based on QM calculations.

Lanmodulin peptides – unravelling the binding of the EF-Hand loop sequences stripped from the structural corset

Lanmodulin (LanM), a naturally lanthanide (Ln)-binding protein with a remarkable selectivity for Lns over Ca(ii) and affinities in the picomolar range, is an attractive target to address challenges in Ln separation. Why LanM has such a high selectivity is currently not entirely understood; both specific amino acid sequences of the EF-Hand loops and cooperativity effects have been suggested. Here, we removed the effect of cooperativity and synthesised all four 12-amino acid EF-Hand loop peptides, and investigated their affinity for two Lns (Eu(iii) and Tb(iii)), the actinide Cm(iii) and Ca(ii). Using isothermal titration calorimetry and time-resolved laser fluorescence spectroscopy (TRLFS) combined with parallel factor analysis, we show that the four short peptides behave very similarly, having affinities in the micromolar range for Eu(iii) and Tb(iii). Ca(ii) was shown not to bind to the peptides, which was verified with circular dichroism spectroscopy. This technique also revealed an increase in structural organisation upon Eu(iii) addition, which was supported by molecular dynamics simulations. Lastly, we put Eu(iii) and Cm(iii) in direct competition using TRLFS. Remarkably, a slightly higher affinity for Cm(iii) was found. Our results demonstrate that the picomolar affinities in LanM are largely an effect of pre-structuring and therefore a reduction of flexibility in combination with cooperative effects, and that all EF-Hand loops possess similar affinities when detached from the protein backbone, albeit still retaining the high selectivity for lanthanides and actinides over calcium.

Taking a closer look at Lanmodulin’s remarkable selectivity for lanthanides (Ln) over Ca(ii) and high Ln/actinide affinities on the amino acid level by investigating the four binding-loops as peptides with Ca(ii), Eu(iii), Tb(iii) and Cm(iii).