Exploring the conformational and binding properties of unphosphorylated/phosphorylated monomeric and trimeric Bcl-2 through docking and molecular dynamics simulations

B-cell lymphoma (Bcl-2) is commonly associated with the progression and preservation of cancer and certain lymphomas; therefore, it is considered as a biological target against cancer. Nevertheless, evidence of all its structural binding sites has been hidden because of the lack of a complete Bcl-2 model, given the presence of a flexible loop domain (FLD), which is responsible for its complex behavior. FLD region has been implicated in phosphorylation, homotrimerization, and heterodimerization associated with Bcl-2 antiapoptotic function. In this contribution, homology modeling, molecular dynamics (MD) simulations in the microsecond (μs) time-scale and docking calculations were combined to explore the conformational complexity of unphosphorylated/phosphorylated monomeric and trimeric Bcl-2 systems. Conformational ensembles generated through MD simulations allowed for identifying the most populated unphosphorylated/phosphorylated monomeric conformations, which were used as starting models to obtain trimeric complexes through protein-protein docking calculations, also submitted to μs MD simulations. Principal component analysis showed that FLD represents the main contributor to total Bcl-2 mobility, and is affected by phosphorylation and oligomerization. Subsequently, based on the most representative unphosphorylated/phosphorylated monomeric and trimeric Bcl-2 conformations, docking studies were initiated to identify the ligand binding site of several known Bcl-2 inhibitors to explain their influence in homo-complex formation and phosphorylation. Docking studies showed that the different conformational states experienced by FLD, such as phosphorylation and oligomerization, play an essential role in the ability to make homo and hetero-complexes.