The role of loop dynamics in the prediction of ligand-protein binding enthalpy

The enthalpic and entropic components of ligand-protein binding free energy reflect the interactions and dynamics between ligand and protein. Despite decades of study, our understanding and hence our ability to predict these individual components remains poor. In recent years, there has been substantial effort and success in the prediction of relative and absolute binding free energies, but the prediction of the enthalpic (and entropic) contributions in biomolecular systems remains challenging. Indeed, it is not even clear what kind of performance in terms of accuracy could currently be obtained for such systems. It is, however, relatively straight-forward to compute the enthalpy of binding. We thus evaluated the performance of absolute enthalpy of binding calculations using molecular dynamics simulation for ten inhibitors against a member of the bromodomain family, BRD4-1, against isothermal titration calorimetry data. Initial calculations, with the AMBER force-field showed good agreement with experiment (R² = 0.60) and surprisingly good accuracy with an average of root-mean-square error (RMSE) = 2.49 kcal mol^-1. Of the ten predictions, three were obvious outliers that were all over-predicted compared to experiment. Analysis of various simulation factors, including parameterization, buffer concentration and conformational dynamics, revealed that the behaviour of a loop (the ZA loop on the periphery of the binding site) strongly dictates the enthalpic prediction. Consistent with previous observations, the loop exists in two distinct conformational states and by considering one or the other or both states, the prediction for the three outliers can be improved dramatically to the point where the R² = 0.95 and the accuracy in terms of RMSE improves to 0.90 kcal mol^-1. However, performance across force-fields is not consistent: if OPLS and CHARMM are used, different outliers are observed and the correlation with the ZA loop behaviour is not recapitulated, likely reflecting parameterization as a confounding problem. The results provide a benchmark standard for future study and comparison.

Unraveling the Abnormal Molecular Mechanism of Suicide Inhibition of Cytochrome P450 3A4

Suicide inhibition of the CYP3A4 enzyme by a drug inactivates the enzyme in the drug biotransformation process and often shows safety concerns about the drug. Despite extensive experimental studies, the abnormal molecular mechanism of a suicide inhibitor that forms a covalent bond with the residue far away from the catalytically active center of CYP3A4 inactivating the enzyme remains elusive. Here, the authors used molecular simulation approaches to study in detail how diquinone methide (DQR), the metabolite product of raloxifene, unbinds from CYP3A4 and inactivates the enzyme at the atomistic level. The results clearly indicate that in one of the intermediate states formed in its unbinding process, DQR covalently binds to Cys239, a residue far away from the catalytically active center of CYP3A4, and hinders the substrate from entering or leaving the enzyme. This work therefore provides an unprecedented way of clarifying the abnormal mechanism of suicide inhibition of the CYP3A4 enzyme.

Exploration of Somatostatin Binding Mechanism to Somatostatin Receptor Subtype 4

Somatostatin (also named as growth hormone-inhibiting hormone or somatotropin release-inhibiting factor) is a regulatory peptide important for the proper functioning of the endocrine system, local inflammatory reactions, mood and motor coordination, and behavioral responses to stress. Somatostatin exerts its effects via binding to G-protein-coupled somatostatin receptors of which the fourth subtype (SSTR4) is a particularly important receptor mediating analgesic, anti-inflammatory, and anti-depressant effects without endocrine actions. Thus, SSTR4 agonists are promising drug candidates. Although the knowledge of the atomic resolution-binding modes of SST would be essential for drug development, experimental elucidation of the structures of SSTR4 and its complexes is still awaiting. In the present study, structures of the somatostatin–SSTR4 complex were produced using an unbiased, blind docking approach. Beyond the static structures, the binding mechanism of SST was also elucidated in the explicit water molecular dynamics (MD) calculations, and key binding modes (external, intermediate, and internal) were distinguished. The most important residues on both receptor and SST sides were identified. An energetic comparison of SST binding to SSTR4 and 2 offered a residue-level explanation of receptor subtype selectivity. The calculated structures show good agreement with available experimental results and indicate that somatostatin binding is realized via prerequisite binding modes and an induced fit mechanism. The identified binding modes and the corresponding key residues provide useful information for future drug design targeting SSTR4.

Dynamics of the Histone Acetyltransferase Lysine-Rich Loop in the Catalytic Core of the CREB-Binding Protein

The tight control of transcriptional coactivators is a fundamental aspect of gene expression in cells. The regulation of the CREB-binding protein (CBP) and p300 coactivators, two paralog multidomain proteins, involves an autoinhibitory loop (AIL) of the histone acetyltransferase (HAT) domain. There is experimental evidence for the AIL engaging with the HAT binding site, thus interrupting the acetylation of histone tails or other proteins. Both CBP and p300 contain a domain of about 110 residues (called the bromodomain) that recognizes histone tails with one or more acetylated lysine side chains. Here, we investigate by molecular dynamics simulations whether the AIL of CBP (residues 1556-1618) acetylated at the side chain of Lys1595 can bind to the bromodomain. The structural instability and fast unbinding kinetics of the AIL from the bromodomain pocket suggest that the AIL is not a ligand of the bromodomain on the same protein chain. This is further supported by the absence of strong and persistent contacts at the binding interface. Furthermore, the simulations of unbinding show an initial fast detachment of the acetylated lysine and a slower phase necessary for complete AIL dissociation. We provide further evidence for the instability of the AIL intramolecular binding by comparison with a natural ligand, the histone peptide H3K56ac, which shows higher stability in the pocket.

ATAD2 interacts with C/EBPβ to promote esophageal squamous cell carcinoma metastasis via TGF-β1/Smad3 signaling

Background

Distant metastasis is the leading cause of death for esophageal squamous cell carcinoma (ESCC) with limited treatment options and unsatisfactory effectiveness. Bromodomain (BRD) containing proteins are emerging targets for cancer therapy with promising effects. As a unique member of BRD family, the function and molecular mechanism of ATAD2 in cancer development is seldomly investigated.

Methods

The clinical impact of ATAD2 was assessed both at RNA and protein level in 75 and 112 ESCC patients separately. The biological function of ATAD2 was investigated in vitro and in vivo. Signaling pathway and downstream effectors of ATAD2 were identified by RNA sequencing, luciferase reporter, co-immunoprecipitation, chromatin immunoprecipitation, immunofluorescence and western blot assay.

Results

We found that elevated ATAD2 expression was significantly associated with lymph node metastasis, advanced clinical stage as well as poor survival of ESCC patients. Silencing ATAD2 significantly suppressed ESCC cell migration and invasion in vitro, and inhibited tumor growth and lung metastasis in vivo. Mechanically, we identified a new cofactor, C/EBPβ. ATAD2 directly interacted with C/EBPβ and promoted its nuclear translocation, which directly bound to the promoter region of TGF-β1 and activated its expression. Further, we demonstrated that TGF-β1 activated its downstream effectors in a Smad3 dependent manner. In addition, we further found that ATAD2 promoted ESCC metastasis through TGF-β signaling induced Snail expression and the subsequent epithelial-mesenchymal transition.

Conclusion

Our findings demonstrated the pro-metastatic function of ATAD2 and uncovered the new molecular mechanism by regulating C/EBPβ/TGF-β1/Smad3/Snail signaling pathway, thus providing a potential target for the treatment of ESCC metastasis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13046-021-01905-x.

Emerging oncogene ATAD2: Signaling cascades and therapeutic initiatives

ATAD2 is a promising oncoprotein with tumor-promoting functions in many cancers. It is a valid cancer drug-target and a potential cancer-biomarker for multiple malignancies. As a cancer/testis antigen (CTA), ATAD2 could also be a probable candidate for immunotherapy. It is a unique CTA that belongs to both AAA+ ATPase and bromodomain family proteins. Since 2007, several research groups have been reported on the pleiotropic oncogenic functions of ATAD2 in diverse signaling pathways, including Rb/E2F-cMyc pathway, steroid hormone signaling pathway, p53 and p38-MAPK-mediated apoptotic pathway, AKT pathway, hedgehog signaling pathway, HIF1α signaling pathway, and Epithelial to Mesenchymal Transition (EMT) pathway in various cancers. In all these pathways, ATAD2 participates in chromatin dynamics, DNA replication, and gene transcription, demonstrating its role as an epigenetic reader and transcription factor or coactivator to promote tumorigenesis. However, despite the progress, an overall mechanism of ATAD2-mediated oncogenesis in diverse origin is elusive. In this review, we summarize the accumulated evidence to envision the overall ATAD2 signaling networks during carcinogenesis and highlight the area where missing links await further research. Besides, the structure-function aspect of ATAD2 is also discussed. Since the efforts have already been initiated to explore targeted drug molecules and RNA-based therapeutic alternatives against ATAD2, their potency and prospects have been elucidated. Together, we believe this is a well-rounded review on ATAD2, facilitating a new drift in ATAD2 research, essential for its clinical implication as a biomarker and/or cancer drug-target.

Discovery of a hidden transient state in all bromodomain families

Bromodomains (BDs) are small protein modules that interact with acetylated marks in histones. These posttranslational modifications are pivotal to regulate gene expression, making BDs promising targets to treat several diseases. While the general structure of BDs is well known, their dynamical features and their interplay with other macromolecules are poorly understood, hampering the rational design of potent and selective inhibitors. Here, we combine extensive molecular dynamics simulations, Markov state modeling, and available structural data to reveal a transiently formed state that is conserved across all BD families. It involves the breaking of two backbone hydrogen bonds that anchor the ZA-loop with the αA helix, opening a cryptic pocket that partially occludes the one associated to histone binding. By analyzing more than 1,900 experimental structures, we unveil just two adopting the hidden state, explaining why it has been previously unnoticed and providing direct structural evidence for its existence. Our results suggest that this state is an allosteric regulatory switch for BDs, potentially related to a recently unveiled BD-DNA-binding mode.

Dissecting the Structural Plasticity and Dynamics of Cytochrome P450 2B4 by Molecular Dynamics Simulations

The plasticity of cytochromes P450 (P450s) is known to contribute significantly to their catalytic capacity of metabolizing various substrates. Although numerous studies have been performed, factors governing the plasticity and dynamics of P450s are still not fully understood. In this study, taking CYP2B4 as an example, we dissect the protein plasticity and dynamics in different environments. CYP2B4 is featured by a high degree of plasticity, which exhibits open, closed, and intermediate states. By analyzing the CYP2B4 crystal structures, we identified the structural features for the closed, open, and intermediate states. Interestingly, formation of the dimer structure was found in the open and intermediate states. The subsequent molecular dynamics (MD) simulations of the open structure in water confirmed the importance of the dimer form in stabilizing the open conformations. MD simulations of the closed and open structures in the membrane environment and the free energies for opening the F-G cassette obtained from the umbrella sampling calculations indicate that the membrane environment is important for stabilizing the F-G cassette. The dynamical network analysis indicates that Asp105 on the B-C loop plays an important role in transiting the structure from the open to the intermediate state. Our results thus unveil the mechanisms of dimer formation and open-to-intermediate transition for CYP2B4 in the water and membrane environments.