Search

Supporting Information: The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/jacsau.2c00358. MD and SMD simulation parameters, data analysis details of all 50 replicas of SMD simulations, all pre-processed simulation trajectories, and data analysis scripts available at https://unmm-my.sharepoint.com/:f:/g/personal/yihe_unm_edu/EqwwWVZFiTJCmhAXnVCCHkwBgmyDGObyy_F60RJyZKgzFQ?e=rnMbfN (PDF). Copyright © 2022 The Authors. The dynamic association and dissociation between proteins are the basis of cellular signal transduction. This process becomes much more complicated if one or both interaction partners are intrinsically disordered because intrinsically disordered proteins can undergo disorder-to-order transitions upon binding to their partners. p53, a transcription factor with disordered regions, plays significant roles in many cellular signaling pathways. It is critical to understand the binding/unbinding mechanism involving these disordered regions of p53 at the residue level to reveal how p53 performs its biological functions. Here, we studied the dissociation process of the intrinsically disordered N-terminal transactivation domain 2 (TAD2) of p53 and the transcriptional adaptor zinc-binding 2 (Taz2) domain of transcriptional coactivator p300 using a combination of classical molecular dynamics, steered molecular dynamics, self-organizing maps, and time-resolved force distribution analysis (TRFDA). We observed two different dissociation pathways with different probabilities. One dissociation pathway starts from the TAD2 N-terminus and propagates to the α-helix and finally the C-terminus. The other dissociation pathway is in the opposite order. Subsequent TRFDA results reveal that key residues in TAD2 play critical roles. Besides the residues in agreement with previous experimental results, we also highlighted some other residues that play important roles in the disassociation process. In the dissociation process, non-native interactions were formed to partially compensate for the energy loss due to the breaking of surrounding native interactions. Moreover, our statistical analysis results of other experimentally determined complex structures involving either Taz2 or TAD2 suggest that the binding of the Taz2-TAD2 complex is mainly governed by the binding site of Taz2, which includes three main binding regions. Therefore, the complexes involving Taz2 may follow similar binding/unbinding behaviors, which could be studied together to generate common principles. National Science Foundation Graduate Research Fellowship Program (grant no. DGE1939267); the National Science Foundation (grant no. 2137558); Leverhulme Trust (RPG-2017-222); Substance Use Disorders Grand Challenge Pilot Research Award; Research Allocations Committee (RAC) Award; the startup fund from the University of New Mexico; the University of New Mexico Office of the Vice President for Research WeR1 Faculty Success Program.

© 2021 The Authors. Allosteric drugs have been attracting increasing interest over the past few years. In this context, it is common practice to use high-throughput screening for the discovery of non-natural allosteric drugs. While the discovery stage is supported by a growing amount of biological information and increasing computing power, major challenges still remain in selecting allosteric ligands and predicting their effect on the target protein’s function. Indeed, allosteric compounds can act both as inhibitors and activators of biological responses. Computational approaches to the problem have focused on variations on the theme of molecular docking coupled to molecular dynamics with the aim of recovering information on the (long-range) modulation typical of allosteric proteins. AIRC IG 2017 - ID. 20019 project; AIRC Fellowship; EC Research Innovation Action H2020 Programme Project HPC-EUROPA3 (INFRAIA-2016-1-730897)

© 2021 The Authors. Several experimental studies indicated that large conformational changes, including partial domain unfolding, have a role in the functional mechanisms of the basic helix loop helix Per/ARNT/SIM (bHLH-PAS) transcription factors. Recently, single-molecule atomic force microscopy (AFM) revealed two distinct pathways for the mechanical unfolding of the ARNT PAS-B. In this work we used steered molecular dynamics simulations to gain new insights into this process at an atomistic level. To reconstruct and classify pathways sampled in multiple simulations, we designed an original approach based on the use of self-organizing maps (SOMs). This led us to identify two types of unfolding pathways for the ARNT PAS-B, which are in good agreement with the AFM findings. Analysis of average forces mapped on the SOM revealed a stable conformation of the PAS-B along one pathway, which represents a possible structural model for the intermediate state detected by AFM. The approach here proposed will facilitate the study of other signal transmission mechanisms involving the folding/unfolding of PAS domains. National Institutes of Environmental Health Sciences (R01-ES007685); Leverhulme Trust (RPG-2017-222).

The aryl hydrocarbon receptor (AhR) is a ligand-dependent transcription factor that is activated by a structurally diverse array of synthetic and natural chemicals, including toxic halogenated aromatic hydrocarbons such as 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). Analysis of the occurring in the AhR ligand binding and activation processes requires structural information on the AhR Per-Arnt-Sim (PAS) B-containing ligand binding domain, for which no experimentally determined structure has been reported. With the availability of extensive structural information on homologous PAS-containing proteins, a reliable model of the mouse AhR PAS B domain was developed by comparative modeling techniques. The PAS domain structures of the functionally related hypoxia-inducible factor 2α (HIF-2α) and AhR nuclear translocator (ARNT) proteins, which exhibit the highest degree of sequence identity and similarity with AhR, were chosen to develop a two-template model. To confirm the features of the modeled domain, the effects of point mutations in selected residue positions on both TCDD binding to the AhR and TCDD-dependent transformation and DNA binding were analyzed. Mutagenesis and functional analysis results are consistent with the proposed model and confirm that the cavity modeled in the interior of the domain is indeed involved in ligand binding. Moreover, the physicochemical characteristics of some residues and of their mutants, along with the effects of mutagenesis on TCDD and DNA binding, also suggest some key features that are required for ligand binding and activation of mAhR at a molecular level, thus providing a framework for further studies. © 2007 American Chemical Society.