Intrinsic Dynamics of the Binding Rail and Its Allosteric Effect in the Class I Histone Deacetylases

Density functional theory (DFT) studies in HDAC-based chemotherapeutics: Current findings, case studies and future perspectives

Density Functional Theory (DFT) is a quantum chemical computational method used to predict and analyze the electronic properties of atoms, molecules, and solids based on the density of electrons rather than wavefunctions. It provides insights into the structure, bonding, and behavior of different molecules, including those involved in the development of chemotherapeutic agents, such as histone deacetylase inhibitors (HDACis). HDACs are a wide group of metalloenzymes that facilitate the removal of acetyl groups from acetyl-lysine residues situated in the N-terminal tail of histones. Abnormal HDAC recruitment has been linked to several human diseases, especially cancer. Therefore, it has been recognized as a prospective target for accelerating the development of anticancer therapies. Researchers have studied HDACs and its inhibitors extensively using a combination of experimental methods and diverse in-silico approaches such as machine learning and quantitative structure-activity relationship (QSAR) methods, molecular docking, molecular dynamics, pharmacophore mapping, and more. In this context, DFT studies can make significant contribution by shedding light on the molecular properties, interactions, reaction pathways, transition states, reactivity and mechanisms involved in the development of HDACis. This review attempted to elucidate the scope in which DFT methodologies may be used to enhance our comprehension of the molecular aspects of HDAC inhibitors, aiding in the rational design and optimization of these compounds for therapeutic applications in cancer and other ailments. The insights gained can guide experimental efforts toward developing more potent and selective HDAC inhibitors.

Mono-phosphorylation at Ser4 of barrier-to-autointegration factor (Banf1) significantly reduces its DNA binding capability by inducing critical changes in its local conformation and DNA binding surface

Barrier-to-autointegration factor (Banf1) is a small DNA-bridging protein. The binding status of Banf1 to DNA is regulated by its N-terminal phosphorylation and dephosphorylation, which plays a critical role in cell proliferation. Banf1 can be phosphorylated at Ser4 into mono-phosphorylated Banf1, which is further phosphorylated at Thr3 to form di-phosphorylated Banf1. It was observed decades ago that mono-phosphorylated Banf1 cannot bind to DNA. However, the underlying molecular- and atomic-level mechanisms remain unclear. A clear understanding of these mechanisms will aid in interfering with the cell proliferation process for better global health. Herein, we explored the detailed atomic bases of unphosphorylated Banf1-DNA binding and how mono- and di-phosphorylation of Banf1 impair these atomic bases to eliminate its DNA-binding capability, followed by exploring the DNA-binding capability of mono- and di-phosphorylation Banf1, using comprehensive and systematic molecular modelling and molecular dynamics simulations. This work presented in detail the residue-level binding energies, hydrogen bonds and water bridges between Banf1 and DNA, some of which have not been reported. Moreover, we revealed that mono-phosphorylation of Banf1 causes its N-terminal secondary structure changes, which in turn induce significant changes in Banf1's DNA binding surface, thus eliminating its DNA-binding capability. At the atomic level, we also uncovered the alterations in interactions due to the induction of mono-phosphorylation that result in the N-terminal secondary structure changes of Banf1. Additionally, our modelling showed that phosphorylated Banf1 with their dominant N-terminal secondary structures bind to DNA with a significantly lower affinity and the docked binding pose are not stable in MD simulations. These findings help future studies in predicting effect of mutations in Banf1 on its DNA-binding capability and open a novel avenue for the development of therapeutics such as cancer drugs, targeting cell proliferation by inducing conformational changes in Banf1's N-terminal domain.

Crystal structures of Schistosoma mansoni histone deacetylase 8 reveal a novel binding site for allosteric inhibitors

Parasitic diseases cause significant global morbidity and mortality particularly in the poorest regions of the world. Schistosomiasis, one of the most widespread neglected tropical diseases, affects more than 200 million people worldwide. Histone deacetylase (HDAC) inhibitors are prominent epigenetic drugs that are being investigated in the treatment of several diseases, including cancers and parasitic diseases. Schistosoma mansoni HDAC8 (SmHDAC8) is highly expressed in all life cycle stages of the parasite and selective inhibition is required in order to avoid undesirable off-target effects in the host. Herein, by X-ray crystal structures of SmHDAC8-inhibitor complexes, biochemical and phenotypic studies, we found two schistosomicidal spiroindoline-derivatives binding a novel site, next to Trp198, on the enzyme surface. We determined that by acting on this site, either by mutation of the Trp198 or by compound binding, a decrease in the activity of the enzyme is achieved. Remarkably, this allosteric site differs from the human counterpart; rather, it is conserved in all Schistosoma spp., as well as Rhabidoptera and Trematoda classes, thus paving the way for the design of HDAC8-selective allosteric inhibitors with improved properties.

Structure-Based Inhibitor Discovery of Class I Histone Deacetylases (HDACs)

Class I histone deacetylases (HDACs) are promising targets for epigenetic therapies for a range of diseases such as cancers, inflammations, infections and neurological diseases. Although six HDAC inhibitors are now licensed for clinical treatments, they are all pan-inhibitors with little or no HDAC isoform selectivity, exhibiting undesirable side effects. A major issue with the currently available HDAC inhibitors is that they have limited specificity and target multiple deacetylases. Except for HDAC8, Class I HDACs (1, 2 and 3) are recruited to large multiprotein complexes to function. Therefore, there are rising needs to develop new, hopefully, therapeutically efficacious HDAC inhibitors with isoform or complex selectivity. Here, upon the introduction of the structures of Class I HDACs and their complexes, we provide an up-to-date overview of the structure-based discovery of Class I HDAC inhibitors, including pan-, isoform-selective and complex-specific inhibitors, aiming to provide an insight into the discovery of additional HDAC inhibitors with greater selectivity, specificity and therapeutic utility.

Allostery of multidomain proteins with disordered linkers

Intrinsically disordered regions are often involved in allosteric regulation of multidomain proteins. They can act as disordered linkers to connect and interact with domains, resulting in rather complex allosteric mechanism and novel protein behavior. Therefore, it is necessary to analyze the diverse functions of disordered linkers in order to better understand allostery and relevant regulation process. Here we summarize recent advances in understanding the function of linkers and the advantages of adopting mutlidomain architecture with disorder linkers. It was shown that linkers between domains enhance the local domain concentration and make the allosteric regulation of weakly interacting partners possible, while linkers with only one tethered end cause an entropy effect to reduce binding affinity and prevent aggregation.

Conformational dynamics and allosteric effect modulated by the unique zinc-binding motif in class IIa HDACs

Class IIa histone deacetylases (HDACs) have been considered as potential targets for the treatment of several diseases. Compared to other HDACs, class IIa HDACs have an additional second zinc binding motif. So far, the function of the unique zinc-binding motif is still not very clear. In this work, extensive classical molecular dynamics (MD) simulations were employed to illuminate the conformational change modulated by the unique zinc-binding motif. It has been revealed that the unique zinc-binding motif is a crucial structural stabilization factor in retaining the catalytic activity of the enzyme and the stability of the multi-protein complex, by remotely modulating the active site pocket in a "closed" conformation. Moreover, it is also revealed that the Loop2 motion in HDAC4 is less flexible than that in HDAC7, which opens a new avenue to design selective inhibitors by utilizing the local conformational dynamics difference between the structurally highly similar HDAC4 and HDAC7. Finally, by comparative studies with class I HDACs (HDAC1-3), it is found that the reversible "in-out" conformational transformation of the binding rail (highly conserved both in class I and IIa HDACs) occurs spontaneously in HDAC1-3, whereas the binding rail is trapped in an "in" conformation owing to the strong metal coordination interaction of the unique CCHC zinc-binding motif in class IIa HDACs. Thus, the CCHC zinc-binding motif may be a feasible allosteric site for the development of class IIa-selective inhibitors.

Fluctuation correlations as major determinants of structure- and dynamics-driven allosteric effects

Both structure- and dynamics-driven allosteric effects are determined by the correlation of distance fluctuations in proteins.