Recent Advances on Small-Molecule Bromodomain-Containing Histone Acetyltransferase Inhibitors

In recent years, substantial research has been conducted on molecular mechanisms and inhibitors targeting bromodomains (BRDs) and extra-terminal (BET) family proteins. On this basis, non-BET BRD is gradually becoming a research hot spot. BRDs are abundant in histone acetyltransferase (HAT)-associated activating transcription factors, and BRD-containing HATs have been linked to cancer, inflammation, and viral replication. Therefore, the development of BRD-containing HATs as chemical probes is useful for understanding the specific biological roles of BRDs in diseases and drug discovery. Several types of BRD-containing HATs, including CBP/P300, PCAF/GCN5, and TAF1, are discussed in this context in terms of their structures, functions, and small-molecule inhibitors. Additionally, progress in BRD inhibitors/chemical probes and proteolysis targeting chimeras in terms of drug design, biological activity, and disease application are summarized. These findings provide insights into the development of BRD inhibitors as potential drug candidates for various diseases.

Atomic Mechanisms of Transthyretin Tetramer Dissociation Studied by Molecular Dynamics Simulations

The dissociation of the transthyretin (TTR) tetramer into a monomer is closely related to various TTR amyloidoses in humans. While the tetramer dissociation has been reported to be the rate-limiting step for TTR aggregation, few details are known about the mechanism. Here, molecular dynamics (MD) simulations were performed by combining conventional MD and biased metadynamics to investigate the mechanism for the wild-type (WT) and mutant (T119M) structures. Both were found to have a great deal in common. Conventional MD simulations reveal that interfacial hydrophobic interactions contribute significantly to stabilize the tetramer. Interfacial residues including L17, V20, L110, and V121 with close contacts form a hydrophobic channel. Metadynamics simulations indicate that the mouth opening of the hydrophobic channel is the first and the most difficult step for dissociation. Interactions of V20 between opposing dimers lock four monomers into the tetramer, and disruption of the interactions is found to be involved in the final step. During the dissociation, an increasing extent of solvation was observed by calculating the radial distribution functions of water around interfacial hydrophobic residues, suggesting that water plays a role in driving the tetramer dissociation. Moreover, compared to T119, residue M119 has a longer side chain that extends into the hydrophobic channel, making solvation more difficult, consistent with a higher energy barrier for dissociation of the T119M tetramer. This result provides a good explanation for the protective role of the T119M mutation. Overall, this study can provide atomic-level insights to better understand the pathogenesis of TTR amyloidosis and guide rational drug design in the future.

Exploring the Binding Mechanism of a Supramolecular Tweezer CLR01 to 14-3-3σ Protein via Well-Tempered Metadynamics

Using supramolecules for protein function regulation is an effective strategy in chemical biology and drug discovery. However, due to the presence of multiple binding sites on protein surfaces, protein function regulation via selective binding of supramolecules is challenging. Recently, the functions of 14-3-3 proteins, which play an important role in regulating intracellular signaling pathways via protein–protein interactions, have been modulated using a supramolecular tweezer, CLR01. However, the binding mechanisms of the tweezer molecule to 14-3-3 proteins are still unclear, which has hindered the development of novel supramolecules targeting the 14-3-3 proteins. Herein, the binding mechanisms of the tweezer to the lysine residues on 14-3-3σ (an isoform in 14-3-3 protein family) were explored by well-tempered metadynamics. The results indicated that the inclusion complex formed between the protein and supramolecule is affected by both kinetic and thermodynamic factors. In particular, simulations confirmed that K214 could form a strong binding complex with the tweezer; the binding free energy was calculated to be −10.5 kcal·mol⁻¹ with an association barrier height of 3.7 kcal·mol⁻¹. In addition, several other lysine residues on 14-3-3σ were identified as being well-recognized by the tweezer, which agrees with experimental results, although only K214/tweezer was co-crystallized. Additionally, the binding mechanisms of the tweezer to all lysine residues were analyzed by exploring the representative conformations during the formation of the inclusion complex. This could be helpful for the development of new inhibitors based on tweezers with more functions against 14-3-3 proteins via modifications of CLR01. We also believe that the proposed computational strategies can be extended to understand the binding mechanism of multi-binding sites proteins with supramolecules and will, thus, be useful toward drug design.

Mechanisms Underlying Gastrodin Alleviating Vincristine-Induced Peripheral Neuropathic Pain

Gastrodin (GAS) is the main bioactive ingredient of Gastrodia, a famous Chinese herbal medicine widely used as an analgesic, but the underlying analgesic mechanism is still unclear. In this study, we first observed the effects of GAS on the vincristine-induced peripheral neuropathic pain by alleviating the mechanical and thermal hyperalgesia. Further studies showed that GAS could inhibit the current density of Na_V1.7 and Na_V1.8 channels and accelerate the inactivation process of Na_V1.7 and Na_V1.8 channel, thereby inhibiting the hyperexcitability of neurons. Additionally, GAS could significantly reduce the over-expression of Na_V1.7 and Na_V1.8 on DRG neurons from vincristine-treated rats according to the analysis of Western blot and immunofluorescence results. Moreover, based on the molecular docking and molecular dynamic simulation, the binding free energies of the constructed systems were calculated, and the binding sites of GAS on the sodium channels (Na_V1.7 and Na_V1.8) were preliminarily determined. This study has shown that modulation of Na_V1.7 and Na_V1.8 sodium channels by GAS contributing to the alleviation of vincristine-induced peripheral neuropathic pain, thus expanding the understanding of complex action of GAS as a neuromodulator.

The prediction of protein-ligand unbinding for modern drug discovery

INTRODUCTION

Drug-target thermodynamic and kinetic information have perennially important roles in drug design. The prediction of protein-ligand unbinding, which can provide important kinetic information, in experiments continues to face great challenges. Uncovering protein-ligand unbinding through molecular dynamics simulations has become efficient and inexpensive with the progress and enhancement of computing power and sampling methods.

AREAS COVERED

In this review, various sampling methods for protein-ligand unbinding and their basic principles are firstly briefly introduced. Then, their applications in predicting aspects of protein-ligand unbinding, including unbinding pathways, dissociation rate constants, residence time and binding affinity, are discussed.

EXPERT OPINION

Although various sampling methods have been successfully applied in numerous systems, they still have shortcomings and deficiencies. Most enhanced sampling methods require researchers to possess a wealth of prior knowledge of collective variables or reaction coordinates. In addition, most systems studied at present are relatively simple, and the study of complex systems in real drug research remains greatly challenging. Through the combination of machine learning and enhanced sampling methods, prediction accuracy can be further improved, and some problems encountered in complex systems also may be solved.

Targeting protein arginine methyltransferase 5 in cancers: Roles, inhibitors and mechanisms

The protein arginine methyltransferase 5 (PRMT5) as the major type II arginine methyltransferase catalyzes the mono- and symmetric dimethylation of arginine residues in both histone and non-histone proteins. Recently, increasing evidence has demonstrated that PRMT5 plays an indispensable role in the occurrence and development of various human cancers by promoting the cell proliferation, invasion, and migration. It has become a promising and valuable target in the cancer epigenetic therapy. This review is to summarize the clinical significance of PRMT5 in the cancers such as lung cancer, breast cancer and colorectal cancer, and the drug discovery targeting PRMT5. Importantly, the existing PRMT5 inhibitors representing different molecular mechanisms, and their pharmacological effect, mechanism of action and biological affinity are analyzed. Clinical status, current problems and future perspective of PRMT5 inhibitors for the treatment of cancers are also discussed, all of which provides crucial help for the future discovery of PRMT5 targeted drugs for cancer treatment.

Exploring the thermodynamic, kinetic and inhibitory mechanisms of 5-iTU targeting mitotic kinase haspin by integrated molecular dynamics

As a human mitotic kinase, haspin is considered as a promising target for various diseases including cancers. However, no inhibitors targeting haspin have entered clinical trials presently. 5-iTU (5-iodotubercidin) is a useful and classical chemical probe for the investigation of haspin activity, but its inhibitory mechanism remains unclear. In this study, integrated molecular dynamics (MD) of conventional MD, extended adaptive biasing force (eABF), random acceleration MD and well-tempered metadynamics were applied to investigate the thermodynamic and kinetic features of 5-iTU and three derivatives targeting haspin. To emphasize the importance of gatekeeper Phe605, two haspin mutants (F605Y and F605T) were also built. The results showed that the binding affinity of 5-iTU and haspin was highest in all wild type (WT) systems, relying on the strong halogen aromatic π interaction between 5-iTU and gatekeeper Phe605. Gatekeeper mutations, because of damage to this interaction, led to the rearrangement of water distributions at the binding site and the decrease of 5-iTU residence times. Additionally, compared with the smaller 5-fTU, 5-iTU dissociated from WT haspin with more difficulty through distinct unbinding pathways. These findings will provide crucial guidance for the design and development of novel haspin inhibitors and the rational modification of existing inhibitors.

Selectively targeting individual bromodomain: Drug discovery and molecular mechanisms

Bromodomain-containing proteins include bromodomain and extra-terminal (BET) and non-BET families. Due to the conserved bromodomain (BD) module between BD-containing proteins, and especially BETs with each member having two BDs (BD1 and BD2), the high degree of structural similarity makes BD-selective inhibitors much difficult to be designed. However, increasing evidences emphasized that individual BDs had distinct functions and different cellular phenotypes after pharmacological inhibition, and selectively targeting one of the BDs could result in a different efficacy and tolerability profile. This review is to summarize the pioneering progress of BD-selective inhibitors targeting BET and non-BET proteins, focusing on their structural features, biological activity, therapeutic application and experimental/theoretical mechanisms. The present proteolysis targeting chimeras (PROTAC) degraders targeting BDs, and clinical status of BD-selective inhibitors were also analyzed, providing a new insight into future direction of bromodomain-selective drug discovery.

Insights into interaction mechanism of inhibitors E3T, E3H and E3B with CREB binding protein by using molecular dynamics simulations and MM-GBSA calculations

CREB binding protein (CBP) and its paralog E1A binding protein (p300) are related to the development of inflammatory diseases, cancers and other diseases, and have been potential targets for the treatment of human diseases. In this work, interaction mechanism of three small molecules E3T, E3H, and E3B with CBP was investigated by employing molecular dynamics (MD) simulations, principal component analysis (PCA), and molecular mechanics/generalized born surface area (MM-GBSA) method. The results indicate that inhibitor bindings cause the changes of movement modes and structural flexibility of CBP, and van der Waals interactions mostly drive associations of inhibitors with CBP. In the meantime, the results based on inhibitor-residue interactions not only show that eight residues of CBP can strongly interact with E3T, E3H and E3B but also verify that the CH-CH, CH-π, and π-π interactions are responsible for vital contributions in associations of E3T, E3H and E3B with CBP. In addition, the H-O radial distribution functions (RDFs) were computed to assess the stability of hydrogen bonding interactions between inhibitors and CBP, and the obtained information identifies several key hydrogen bonds playing key roles in bindings of E3T, E3H and E3B to CBP. Potential new inhibitors have been proposed.

Structure-based drug optimization and biological evaluation of tetrahydroquinolin derivatives as selective and potent CBP bromodomain inhibitors

CBP bromodomain could recognize acetylated lysine and function as transcription coactivator to regulate transcription and downstream gene expression. Furthermore, CBP has been shown to be related to many human malignancies including acute myeloid leukemia. Herein, we identified DC-CPin734 as a potent CBP bromodomain inhibitor with a TR-FRET IC₅₀ value of 19.5 ± 1.1 nM and over 400-fold of selectivity against BRD4 bromodomains through structure based rational drug design guided iterative chemical modification endeavoring to discover optimal tail-substituted tetrahydroquinolin derivatives. Moreover, DC-CPin734 showed potent inhibitory activity to AML cell line MV4-11 with an IC₅₀ value of 0.55 ± 0.04 μM, and its cellular on-target effects were further evidenced by c-Myc downregulation results. In summary, DC-CPin734 showing good potency, selectivity and anti AML activity could serve as a potent and selective in vitro and in vivo probe of CBP bromodomain and a promising lead compound for future drug development.

Design, synthesis, and biological evaluation of tetrahydroquinolin derivatives as potent inhibitors of CBP bromodomain

CREB-binding protein (CBP) is a large multi-domain protein containing a HAT domain catalyzing transacetylation and a bromodomain responsible for acetylated lysine recognition. CBPs could act as transcription co-activators to regulate gene expression and have been shown to play a significant role in the development and progression of many cancers. Herein, through in silico screening two hit compounds with tetrahydroquinolin methyl carbamate scaffold were discovered, among which DC-CPin7 showed an in vitro inhibitory activity with the TR-FRET IC₅₀ value of 2.5 ± 0.3 μM. We obtained a high-resolution co-crystal structure of the CBP bromodomain in complex with DC-CPin7 to guide following structure-based rational drug design, which yielded over ten DC-CPin7 derivatives with much higher potency, among which DC-CPin711 showed approximately 40-fold potency compared with hit compound DC-CPin7 with an in vitro TR-FRET IC₅₀ value of 63.3 ± 4.0 nM. Notably, DC-CPin711 showed over 150-fold selectivity against BRD4 bromodomains. Moreover, DC-CPin711 showed micromolar level of anti-leukemia proliferation through G1 phase cell cycle arrest and cell apoptosis. In summary, through a combination of computational and crystal-based structure optimization, DC-CPin711 showed potent in vitro inhibitory activities to CBP bromodomain with a decent selectivity towards BRD4 bromodomains and good cellular activity to leukemia cells, which could further be applied to related biological and translational studies as well as serve as a lead compound for future development of potent and selective CBP bromodomain inhibitors.

Decoding molecular mechanism of inhibitor bindings to CDK2 using molecular dynamics simulations and binding free energy calculations

CDK2 can be used as an attractive target for development of efficient inhibitors curing multiple disease relating with CDK2. In this work, molecular dynamics (MD) simulations and binding free energy calculations were coupled to probe conformational changes of CDK2 due to inhibitor associations and binding mechanisms of inhibitors PM1, FMD and X64 to CDK2. The results suggest that the binding strength of FMD and X64 to CDK2 is stronger than that of PM1. Principal component (PC) analysis and cross-correlation map calculations based on the equilibrated MD trajectories demonstrate that the structural difference in inhibitors exerts important impact on motion modes and dynamics behavior of CDK2. Residue-based free energy decomposition method was adopted to estimate the inhibitor-residue spectrum. The results not only efficiently identify the hot interaction spot of inhibitors with CDK2 but also show that the hydrophobic rings R1, R2 and R3 as well as polar groups of three inhibitors play key roles in favorably binding of inhibitors to CDK2. This work is expected to contribute energetic basis and dynamics information to development of promising inhibitors toward CDK2.Communicated by Ramaswamy H. Sarma.

Understanding the selectivity of inhibitors toward PI4KIIIα and PI4KIIIβ based molecular modeling

Molecular dynamics simulations and binding free energy calculations are combined to investigate the selectivity of inhibitors toward type III phosphatidylinositol 4 kinases.