Abstract:The BAX protein is a pro-apoptotic member of the Bcl-2 family, which triggers apoptosis by causing permeabilization of the mitochondrial outer membrane. However, the activation mechanism of BAX is far from being understood. Although a few small-molecule BAX activators have been reported in the literature, their crystal structures in complex with BAX have not been resolved. So far, their binding modes were modeled at most by simple molecular docking efforts. Lack of an in-depth understanding of the activation mechanism of BAX hinders the development of more effective BAX activators. In this work, we employed cosolvent molecular dynamics simulation to detect the potential binding sites on the surface of BAX and performed a long-time molecular dynamics simulation (50 μs in total) to derive the possible binding modes of three BAX activators (i.e., BAM7, BTC-8, and BTSA1) reported in the literature. Our results indicate that the trigger, S184, and vMIA sites are the three major binding sites on the full-length BAX structure. Moreover, the canonical hydrophobic groove is clearly detected on the α9-truncated BAX structure, which is consistent with the outcomes of relevant experimental studies. Interestingly, it is observed that solvent probes bind to the trigger bottom pocket more stably than the PPI trigger site. Each activator was subjected to unbiased molecular dynamics simulations started at the three major binding sites in five parallel jobs. Our MD results indicate that all three activators tend to stay at the trigger site with favorable MM-GB/SA binding energies. BAM7 and BTSA1 can enter the trigger bottom pocket and thereby enhance the movement of the α1−α2 loop, which may be a key factor at the early stage of BAX activation. Our molecular modeling results may provide useful guidance for designing smart biological experiments to further explore BAX activation and directing structure-based efforts toward discovering more effective BAX activators.The Supporting Information is available free of charge at <a class="ext-link" href="/doi/10.1021/acs.jcim.0c01420?goto=supporting-info">https://pubs.acs.org/doi/10.1021/acs.jcim.0c01420</a>.The complex models shown in <a class="internalNav" href="#fig5">Figure </a><a class="internalNav" href="#fig5">5</a>, <a class="internalNav" href="#fig7">Figure </a><a class="internalNav" href="#fig7">7</a>, and all relevant figures in the above material are also provided as PDB-format files (<a class="ext-link" href="/doi/suppl/10.1021/acs.jcim.0c01420/suppl_file/ci0c01420_si_001.zip">ZIP</a>)(Part 1) Parameters and RDF of phenol and methyl-pyrazolone; (Part 2) cosolvent MD simulation and definition of binding sites; (Part 3) MD simulation of the binding of three BAX activators; (Part 4) cluster analysis and MM-GB/SA calculation; (Part 5) dihedral distributions of M137 and F93; (Part 6) analysis of the α1−α2 loop; (Part 7) hydrogen bonds formed by BAX activators; (Part 8) cluster analysis of BAX conformations (<a class="ext-link" href="/doi/suppl/10.1021/acs.jcim.0c01420/suppl_file/ci0c01420_si_002.pdf">PDF</a>)This article has not yet been cited by other publications.

Insights Into Unbinding Mechanisms Upon Two Mutations Investigated By Molecular Dynamics Study Of Gsk3 Beta-Axin Complex: Role Of Packing Hydrophobic Residues

Mechanism of Kinase Inactivation and Nonbinding of FRATide to GSK3β Due to K85M Mutation: Molecular Dynamics Simulation and Normal Mode Analysis.

Effect of Double Mutations K214/A–E215/Q of FRATide on GSK3β: Insights from Molecular Dynamics Simulation and Normal Mode Analysis

Structural basis for the complete loss of GSK3beta catalytic activity due to R96 mutation investigated by molecular dynamics study.

Effect of Mutation K85R on GSK-3beta: Molecular Dynamics Simulation.

Unraveling the role of Arg4 and Arg6 in the auto-inhibition mechanism of GSK3β from molecular dynamics simulation.

Mutation of the Conserved G66 Residue in GS Region Decreased Structural Stability and Activity of Arginine Kinase

Axin and GSK3- Control Smad3 Protein Stability and Modulate TGF- Signaling.

Differential Molecular Assemblies Underlie the Dual Function of Axin in Modulating the Wnt and JNK Pathways

Insights into the Folding and Unfolding Processes of Wild-Type and Mutated SH3 Domain by Molecular Dynamics and Replica Exchange Molecular Dynamics Simulations

Mutational analysis of GSK3β protein kinase together with kinome-wide binding and stability studies suggest context-dependent recognition of kinases by the chaperone HSP90.

Some Insights into Mechanism for Binding and Drug Resistance of Wild Type and I50V V82A and I84V Mutations in HIV-1 Protease with GRL-98065 Inhibitor from Molecular Dynamic Simulations.

Molecular modeling and molecular dynamics simulation studies of the GSK3β/ATP/substrate complex: understanding the unique P+4 primed phosphorylation specificity for GSK3β substrates.

Molecular Dynamics Simulation Reveals Insights into the Mechanism of Unfolding by the A130T/V Mutations within the MID1 Zinc-Binding Bbox1 Domain

Homology modeling and molecular dynamics study of GSK3/SHAGGY-like kinase.

Capturing the Effects of Single Atom Substitutions on the Inhibition Efficiency of Glycogen Synthase Kinase-3β Inhibitors via Markov State Modeling and Experiments

Structural Basis of the Interaction between Human Axin2 and SIAH1 in the Wnt/β-Catenin Signaling Pathway

Glycine Substitution in SH3-SH2 Connector of Hck Tyrosine Kinase Causes Population Shift from Assembled to Disassembled State

Structural and functional effects of phosphopriming and scaffolding in the kinase GSK-3β

Mutation-Induced Impacts on the Switch Transformations of the GDP- and GTP-Bound K-Ras: Insights from Multiple Replica Gaussian Accelerated Molecular Dynamics and Free Energy Analysis.

Analysis of the Binding Sites on BAX and the Mechanism of BAX Activators through Extensive Molecular Dynamics Simulations