Binding affinity of pro-apoptotic BH3 peptides for the anti-apoptotic Mcl-1 and A1 proteins: Molecular dynamics simulations of Mcl-1 and A1 in complex with six different BH3 peptides

Abstract

The anti-apoptotic members of Bcl-2 family of proteins bind to their pro-apoptotic counterparts to induce or prevent cell death.Based on the distinct binding profiles for specific pro-apoptotic BH3 peptides, the anti-apoptotic Bcl-2 proteins can be divided into at least two subclasses. The subclass that includes Bcl-XL binds strongly to Bad BH3 peptide while it has weak binding affinity for the second subclass of Bcl-2 proteins such as Mcl-1 and A1. Anti-apoptotic Bcl-2 proteins are considered to be attractive drug targets for anti-cancer drugs. BH3-mimetic inhibitors such as ABT-737 have been shown to be specific to Bcl-XL subclass while Mcl-1 and A1 show resistance to the same drug. An efficacious inhibitor should target all the anti-apoptotic Bcl-2 proteins. Hence, development of inhibitors selective to Mcl-1 and A1 is of prime importance for targeted cancer therapeutics. The first step to achieve this goal is to understand the molecular basis of high binding affinities of specific pro-apoptotic BH3 peptides for Mcl-1 and A1. To understand the interactions between the BH3 peptides and Mcl-1/A1, we performed multi-nanosecond molecular dynamics (MD) simulations of six complex structures of Mcl-1 and A1. With the exception of Bad, all complex structures were experimentally determined. Bad complex structures were modeled. Our simulation studies identified specific pattern of polar interactions between Mcl-1/A1 and high-affinity binding BH3 peptides. The lack of such polar interactions in Bad peptide complex is attributed to specific basic residues present before and after the highly conserved Leu residue. The close approach of basic residues in Bad and Mcl-1/A1 is hypothesized to be the cause of weak binding affinity. To test this hypothesis, we generated in silico mutants of these basic residues in Bad peptide and Mcl-1/A1 proteins. MD simulations of the mutant systems established the pattern of stable polar interactions observed in high-affinity binding BH3 peptides. We have thus identified specific residue positions in Bad and Mcl-1/A1 responsible for the weak binding affinity. Results from these simulation studies will aid in the development of inhibitors specific to Mcl-1 and A1 proteins.