The interaction between the anti-apoptotic Bcl-2 protein and its antagonist Bax is essential to the regulation of the mitochondrial pathway of apoptosis. For this work, we built models by homology of Bcl-2 full-sequence length in monomeric form (apo-Bcl-2) and in complex with the BH3 domain of Bax (holo-Bcl-2). The Bcl-2 protein was analyzed with its transmembrane domain anchored to a lipidic bilayer of DPPC, imitating physiological conditions. We performed molecular dynamics (MD) simulations using the GROMACS program. Conformational changes showed that the flexible loop domain (FLD) tends to fold on itself and move towards the main core. Furthermore, the BH3 peptide of pro-apoptotic protein Bax, showed an allosteric stabilizing effect on FLD upon being bound to the hydrophobic cleft of the anti-apoptotic protein Bcl-2, causing a reduction in its structural flexibility. However, FLD is distal from the main core of Bcl-2. Principal component analysis (PCA) showed a weak correlation between FLD residues and BH3 peptide from Bax. Upon MD simulations, several new contacts appeared between FLD and some α-helices of the core of Bcl-2, which contribute to maintaining the stability of Bcl-2. This knowledge sheds light on the behavior of Bcl-2 in the cell's native environment.Communicated by Ramaswamy H. Sarma.
Keywords: Bcl-2; apo/holo Bcl2 conformation; apoptosis regulation; flexible loop domain; membrane protein; molecular dynamics simulation; transmembrane domain.