Rational design of new class of BH3-mimetics as inhibitors of the Bcl-xL protein

Characterizing the consensus residue specificity and surface of BCL-2 binding to BH3 ligands using the Knob-Socket model

Cancer cells bypass cell death by changing the expression of the BCL-2 family of proteins, which are apoptotic pathway regulators. Upregulation of pro-survival BCL-2 proteins or downregulation of cell death effectors BAX and BAK interferes with the initiation of the intrinsic apoptotic pathway. In normal cells, apoptosis can occur through pro-apoptotic BH3-only proteins interacting and inhibiting pro-survival BCL-2 proteins. When cancer cells over-express pro-survival BCL-2 proteins, a potential remedy is the sequestration of these pro-survival proteins through a class of anti-cancer drugs called BH3 mimetics that bind in the hydrophobic groove of pro-survival BCL-2 proteins. To improve the design of these BH3 mimetics, the packing interface between BH3 domain ligands and pro-survival BCL-2 proteins was analyzed using the Knob-Socket model to identify the amino acid residues responsible for interaction affinity and specificity. A Knob-Socket analysis organizes all the residues in a binding interface into simple 4 residue units: 3-residue sockets defining surfaces on a protein that pack a 4th residue knob from the other protein. In this way, the position and composition of the knobs packing into sockets across the BH3/BCL-2 interface can be classified. A Knob-Socket analysis of 19 BCL-2 protein and BH3 helix co-crystals reveal multiple conserved binding patterns across protein paralogs. Conserved knob residues such as a Gly, Leu, Ala and Glu most likely define binding specificity in the BH3/BCL-2 interface, whereas other residues such as Asp, Asn, and Val are important for forming surface sockets that bind these knobs. These findings can be used to inform the design of BH3 mimetics that are specific to pro-survival BCL-2 proteins for cancer therapeutics.

Rational Drug Design of Antineoplastic Agents Using 3D-QSAR, Cheminformatic, and Virtual Screening Approaches

Background:

Computer-Aided Drug Design has strongly accelerated the development of novel antineoplastic agents by helping in the hit identification, optimization, and evaluation.

Results:

Computational approaches such as cheminformatic search, virtual screening, pharmacophore modeling, molecular docking and dynamics have been developed and applied to explain the activity of bioactive molecules, design novel agents, increase the success rate of drug research, and decrease the total costs of drug discovery. Similarity, searches and virtual screening are used to identify molecules with an increased probability to interact with drug targets of interest, while the other computational approaches are applied for the design and evaluation of molecules with enhanced activity and improved safety profile.

Conclusion:

In this review are described the main in silico techniques used in rational drug design of antineoplastic agents and presented optimal combinations of computational methods for design of more efficient antineoplastic drugs.

Understanding the microscopic binding mechanism of hydroxylated and sulfated polybrominated diphenyl ethers with transthyretin by molecular docking, molecular dynamics simulations and binding free energy calculations

Polybrominated diphenyl ethers (PBDEs), one typical type of persistent environmental contaminant, have toxicological effects such as disrupting thyroid homeostasis in the human body. The high binding affinities of hydroxylated metabolites of PBDEs (OH-PBDEs) with transthyretin (TTR) were considered to be one major reason for their extraordinary capacity of passing through the blood-brain barrier via competitive thyroid hormone (T4) transport protein binding. Recent findings showed that sulfated PBDEs can be formed in human liver cytosol as phase-II metabolites. However, experimentally determined data for the TTR binding potential of the sulfated PBDEs are still not available. Therefore, molecular docking and molecular dynamics (MD) simulations were employed in the present study to probe the molecular basis of TTR interacting with hydroxylated and sulfated PBDEs at the atomic level. The docking scores of LeDock were used to construct the structure-based predictive model. The calculated results showed that the sulfated PBDEs have stronger affinity for TTR than the corresponding OH-PBDEs. Further analysis of structural characteristics based on MD simulations indicated that upon the binding of PBDE metabolites, the stability of TTR was enhanced and the dissociation rate of the tetrameric protein structure was potentially decreased. Subsequent binding free energy calculations implied that van der Waals interactions are the dominant forces for the binding of these metabolites of PBDEs at the T4 site of TTR. The residues Ser117/Ser117' and Lys15/Lys15' were identified, by both residue energy decomposition and computational alanine-scanning mutagenesis methods, as key residues which play an important role in determining the binding orientations of the -OSO₃^- group of sulfated PBDEs by formation of either hydrogen bonds or electrostatic interactions, respectively. In general, the combination of docking calculations with MD simulations provided a theoretically toxicological assessment for the metabolites of PBDEs, deep insight into the recognition mechanism of TTR for these compounds, and thus more comprehensive understanding of the thyroid-related toxic effects of PBDEs as well.

Innovative computer-aided methods for the discovery of new kinase ligands

Recent evidence points to significant roles played by protein kinases in cell signaling and cellular proliferation. Faulty protein kinases are involved in cancer, diabetes and chronic inflammation. Efforts are continuously carried out to discover new inhibitors for selected protein kinases. In this review, we discuss two new computer-aided methodologies we developed to mine virtual databases for new bioactive compounds. One method is ligand-based exploration of the pharmacophoric space of inhibitors of any particular biotarget followed by quantitative structure-activity relationship-based selection of the best pharmacophore(s). The second approach is structure-based assuming that potent ligands come into contact with binding site spots distinct from those contacted by weakly potent ligands. Both approaches yield pharmacophores useful as 3D search queries for the discovery of new bioactive (kinase) inhibitors.

HybridDock: A Hybrid Protein-Ligand Docking Protocol Integrating Protein- and Ligand-Based Approaches

Structure-based molecular docking and ligand-based similarity search are two commonly used computational methods in computer-aided drug design. Structure-based docking tries to utilize the structural information on a drug target like protein, and ligand-based screening takes advantage of the information on known ligands for a target. Given their different advantages, it would be desirable to use both protein- and ligand-based approaches in drug discovery when information for both the protein and known ligands is available. Here, we have presented a general hybrid docking protocol, referred to as HybridDock, to utilize both the protein structures and known ligands by combining the molecular docking program MDock and the ligand-based similarity search method SHAFTS, and evaluated our hybrid docking protocol on the CSAR 2013 and 2014 exercises. The results showed that overall our hybrid docking protocol significantly improved the performance in both binding affinity and binding mode predictions, compared to the sole MDock program. The efficacy of the hybrid docking protocol was further confirmed using the combination of DOCK and SHAFTS, suggesting an alternative docking approach for modern drug design/discovery.

A study of the structural properties and thermal stability of human Bcl-2 by molecular dynamics simulations

The anti-apoptotic B-cell lymphoma 2 (Bcl-2) protein interacts with several proteins that regulate the apoptotic properties of cells. In this research, we conduct several all-atom molecular dynamics (MD) simulations under high-temperature unfolding conditions, from 400 to 800 K, for 25 ns. These simulations were performed using a model of an engineered Bcl-2 human protein (Bcl-2-Δ22Σ3), which lacks 22 C-terminal residues of the transmembrane domain. The aim of this study is to gain insight into the structural behavior of Bcl-2-Δ22Σ3 by mapping the conformational movements involved in Bcl-2 stability and its biological function. To build a Bcl-2-Δ22Σ3 three-dimensional model, the protein core was built by homology modeling and the flexible loop domain (FLD, residues 33-91) by ab initio methods. Further, the entire protein model was refined by MD simulations. Afterwards, the production MD simulations showed that the FLD at 400 and 500 K has several conformations reaching into the protein core, whereas at 600 K some of the alpha-helices were lost. At 800 K, the Bcl-2 core is destabilized suggesting a possible mechanism for protein unfolding, where the alpha helices 1 and 6 were the most stable, and a reduction in the number of hydrogen bonds initially occurs. In conclusion, the structural changes and the internal protein interactions suggest that the core and the FLD are crucial components of Bcl-2 in its function of regulate ng access to the recognition sites of kinases and caspases.

Combination of ligand- and structure-based methods in virtual screening

The combination of ligand- and structure-based molecular modelling methods has become a common approach in virtual screening. This review describes different strategies for integration of ligand- and structure-based methods which can be divided into sequential, parallel or hybrid approaches. Although no thorough performance comparisons between combined approaches are available, examples of successful applications in prospective and retrospective virtual screening are discussed. Most published studies use a sequential approach, utilising well-documented single methods successfully.

BH3 mimetics: status of the field and new developments

Targeting apoptosis is an attractive approach in cancer therapy. The BH3-only proteins of the BCL-2 family (having only the BCL-2 homology domain BH3) can trigger apoptosis by binding to the prosurvival members of this family and neutralizing their functional activity (sequestration of the proapoptotic Bcl-2 family members). The "BH3 mimetic" concept has prompted the development of small molecules capable of mimicking BH3-only proteins and thus inducing apoptosis. The prototype BH3 mimetic ABT-737 selectively targets the three prosurvival proteins BCL-XL, BCL-2, and BCL-W (but not MCL-1 or A1) and its oral derivative ABT-263 has proved promising in clinical trials. Some putative BH3 mimetics are also tested clinically while others are still being characterized. This article recapitulates the various known BH3 mimetics and presents the recent developments in the field. The latter include (i) the identification of molecular determinants responsible for the specific interactions between BH3 motifs and the binding grooves of prosurvival proteins and (ii) the characterization of new compounds and particularly BH3 mimetics that antagonize either selectively MCL-1 or BCL-2 or a broad range of prosurvival proteins. These data are critical advances toward the discovery of novel anticancer agents.

Screening efficient BH3-mimetics to hBcl-B by means of peptidodynmimetic method

The crucial residues of hBaxBH3 peptide for interaction with hBcl-B, an anti-apoptotic protein, were identified using molecular docking studies on the polypeptides and temperature-specific molecular dynamic simulations performed for the protein-peptide complex at near-physiological conditions (pH 7.0, 1 atmospheric pressure and 0.1 M NaCl). The data from the methods were examined by a 'strong residue contacts' filter strategy and the data analyses of the former and latter methods identified 10 (Q52, K57, S60, L63, K64, R65, G67, D68, D71 & S72) and 3 (S60, E61 & K64) crucial residues of the hBaxBH3 peptide for interacting with the protein, respectively. We have herein demonstrated that BH3-chemical mimetics screened using the pharmacophoric residues of hBaxBH3 obtained from the 'peptidodynmimetic method' were superior in terms of ligand efficiencies, bioavailability and pharmacokinetic properties vis-à-vis that of small molecule BH3-mimetics retrieved using the conventional 'peptidomimetic method'. The unique advantages of the 'peptidodynmimetic method' to identify efficient BH3-mimetics for modulating interfaces (composed of a large number of amino acids) of other anti-apoptotic proteins-BH3-only peptides have also been discussed in detail.

Behavior of solvent-exposed hydrophobic groove in the anti-apoptotic Bcl-XL protein: clues for its ability to bind diverse BH3 ligands from MD simulations

Bcl-XL is a member of Bcl-2 family of proteins involved in the regulation of intrinsic pathway of apoptosis. Its overexpression in many human cancers makes it an important target for anti-cancer drugs. Bcl-XL interacts with the BH3 domain of several pro-apoptotic Bcl-2 partners. This helical bundle protein has a pronounced hydrophobic groove which acts as a binding region for the BH3 domains. Eight independent molecular dynamics simulations of the apo/holo forms of Bcl-XL were carried out to investigate the behavior of solvent-exposed hydrophobic groove. The simulations used either a twin-range cut-off or particle mesh Ewald (PME) scheme to treat long-range interactions. Destabilization of the BH3 domain-containing helix H2 was observed in all four twin-range cut-off simulations. Most of the other major helices remained stable. The unwinding of H2 can be related to the ability of Bcl-XL to bind diverse BH3 ligands. The loss of helical character can also be linked to the formation of homo- or hetero-dimers in Bcl-2 proteins. Several experimental studies have suggested that exposure of BH3 domain is a crucial event before they form dimers. Thus unwinding of H2 seems to be functionally very important. The four PME simulations, however, revealed a stable helix H2. It is possible that the H2 unfolding might occur in PME simulations at longer time scales. Hydrophobic residues in the hydrophobic groove are involved in stable interactions among themselves. The solvent accessible surface areas of bulky hydrophobic residues in the groove are significantly buried by the loop LB connecting the helix H2 and subsequent helix. These observations help to understand how the hydrophobic patch in Bcl-XL remains stable in the solvent-exposed state. We suggest that both the destabilization of helix H2 and the conformational heterogeneity of loop LB are important factors for binding of diverse ligands in the hydrophobic groove of Bcl-XL.

MicroRNA miR-491-5p targeting both TP53 and Bcl-XL induces cell apoptosis in SW1990 pancreatic cancer cells through mitochondria mediated pathway

MicroRNA (miRNA) actively participates in a broad range of cellular processes such as proliferation, differentiation, cell survival and apoptosis. Deregulated expression of miRNA may affect cell growth and eventually lead to cancer. In this study, we found that hsa-miR491-5p (miR491-5p) displays a significantly high level of expression in normal human pancreas tissue versus pancreatic cancer cells. Targeted site prediction indicated that both Bcl-XL and TP53 contain miR-491-5p recognizing sites in their 3' UTRs. Overexpression of miR-491-5p in the pancreatic cancer cell line SW1990 effectively inhibited both endogenous Bcl-XL and TP53 gene expressions. Mutagenesis at the seed match region of both targeted genes further confirmed the specificity of miR491-5p recognition. Cell proliferation rate was inversely related to the increased doses of miR-491-5p. Flow cytometric analysis showed that the proportions of total apoptotic and early apoptotic cells were significantly induced as the dose of miR491-5p increased. Moreover, a mechanistic study indicated that miR-R491-5p-mediated cell apoptosis was associated with the activation of intrinsic mitochondria mediated pathways. miR491-5p also markedly inhibited mitogenic signaling pathways such as STAT3 and PI-3K/Akt, but not Ras/MAPK. Thus, our results demonstrated that miR491-5p could effectively target both Bcl-xL and TP53 and induce cell apoptosis independent of TP53.