Taxifolin ameliorates cerebral ischemia-reperfusion injury in rats through its anti-oxidative effect and modulation of NF-kappa B activation

Summary Infarction in adult rat brain was induced by middle cerebral arterial occlusion (MCAO) followed by reperfusion to examine whether taxifolin could reduce cerebral ischemic reperfusion (CI/R) injury. Taxifolin administration (0.1 and 1.0 μg/kg, i.v.) 60 min after MCAO ameliorated infarction (by 42%±7% and 62%±6%, respectively), which was accompanied by a dramatic reduction in malondialdehyde and nitrotyrosine adduct formation, two markers for oxidative tissue damage. Overproduction of reactive oxygen species (ROS) and nitric oxide (NO) via oxidative enzymes (e.g., COX-2 and iNOS) was responsible for this oxidative damage. Taxifolin inhibited leukocyte infiltration, and COX-2 and iNOS expressions in CI/R-injured brain. Taxifolin also prevented Mac-1 and ICAM-1 expression, two key counter-receptors involved in firm adhesion/transmigration of leukocytes to the endothelium, which partially accounted for the limited leukocyte infiltration. ROS, generated by leukocytes and microglial cells, activated nuclear factor-kappa B (NF-κB) that in turn signaled up-regulation of inflammatory proteins. NF-κB activity in CI/R was enhanced 2.5-fold over that of sham group and was inhibited by taxifolin. Production of both ROS and NO by leukocytes and microglial cells was significantly antagonized by taxifolin. These data suggest that amelioration of CI/R injury by taxifolin may be attributed to its anti-oxidative effect, which in turn modulates NF-κB activation that mediates CI/R injury. Yea-Hwey Wang, Wen-Yen Wang, Chia-Che Chang, and Kuo-Tong Liou contributed equally to this work.     

Structural biology of the Bcl-2 family of proteins

The proteins of the Bcl-2 family are important regulators of programmed cell death. Structural studies of Bcl-2 family members have provided many important insights into their molecular mechanism of action and how members of this family interact with one another. To date, structural studies have been performed on six Bcl-2 family members encompassing both anti- (Bcl-x(L), Bcl-2, KSHV-Bcl-2, Bcl-w) and pro-apoptotic (Bax, Bid) members. They all show a remarkably similar fold despite an overall divergence in amino acid sequence and function (pro-apoptotic versus anti-apoptotic). The three-dimensional structures of Bcl-2 family members consist of two central, predominantly hydrophobic alpha-helices surrounded by six or seven amphipathic alpha-helices of varying lengths. A long, unstructured loop is present between the first two alpha-helices. The structures of the Bcl-2 proteins show a striking similarity to the overall fold of the pore-forming domains of bacterial toxins. This finding led to experiments which demonstrated that Bcl-x(L), Bcl-2, and Bax all form pores in artificial membranes. A prominent hydrophobic groove is present on the surface of the anti-apoptotic proteins. This groove is the binding site for peptides that mimic the BH3 region of various pro-apoptotic proteins such as Bak and Bad. Structures of Bcl-x(L) in complex with these BH3 peptides showed that they bind as an amphipathic alpha-helix and make extensive hydrophobic contacts with the protein. These data have not only helped to elucidate the interactions important for hetero-dimerization of Bcl-2 family members but have also been used to guide the discovery of small molecules that block Bcl-x(L) and Bcl-2 function. In the recently determined structure of the anti-apoptotic Bcl-w protein, the protein was also found to have a hydrophobic groove on its surface capable of binding BH3-containing proteins and peptides. However, in the native protein an additional carboxy-terminal alpha-helix interacts with the hydrophobic groove. This is reminiscent of how the carboxy-terminal alpha-helix of the pro-apoptotic protein Bax binds into its hydrophobic groove. This interaction may play a regulatory role and for Bax may explain why it is found predominately in the cytoplasm prior to activation. The hydrophobic groove of the pro-apoptotic protein, Bid protein, is neither as long nor as deep as that found in Bcl-x(L), Bcl-2, or Bax. In addition, Bid contains an extra alpha-helix, which is located between alpha1 and alpha2 with respect to Bcl-x(L), Bcl-2, and Bax. Although there are still many unanswered questions regarding the exact mechanism by which the Bcl-2 family of proteins modulates apoptosis, structural studies of these proteins have deepened our understanding of apoptosis on the molecular level.     

Biophysical characterization of recombinant human Bcl-2 and its interactions with an inhibitory ligand, antimycin A

Apoptosis is an essential physiological process, regulated by the family of Bcl-2-related proteins. However, the molecular mechanism by which Bcl-2 regulates apoptosis still remains elusive. Here we report the functional studies of recombinant human Bcl-2 with the deletion of 22 residues at the C-terminal membrane-anchoring region (rhBcl-2Delta22). Characterization of rhBcl-2Delta22 showed that the recombinant protein is homogeneous and monodisperse in nondenaturing solutions, stable at room temperature in the presence of a metal chelator, and an alpha-helical protein with unfolding of secondary structure at a T(m) of 62.8 degrees C. Optimal membrane pore formation by rhBcl-2Delta22 required negatively charged phospholipids. The existence of a hydrophobic groove in rhBcl-2Delta22 was demonstrated by the fluorescence enhancement of the hydrophobic ANS probe with which a pro-apoptotic Bak BH3 peptide competed. The respiratory inhibitor antimycin A also bound to the hydrophobic groove of rhBcl-2Delta22 with a K(d) of 0.82 microM. The optimal binding conformation of antimycin A was predicted from molecular docking of antimycin A with the hBcl-2 model created by homology modeling. Antimycin A selectively induces apoptosis in cells overexpressing Bcl-2, suggesting that hydrophobic groove-binding compounds may act as selective apoptotic triggers in tumor cells.     

Structure, dynamics and hydration of the nogalamycin-d(ATGCAT)2Complex determined by NMR and molecular dynamics simulations in solution.

The structure of the 1:1 nogalamycin:d(ATGCAT)2 complex has been determined in solution from high-resolution NMR data and restrained molecular dynamics (rMD) simulations using an explicit solvation model. The antibiotic intercalates at the 5'-TpG step with the nogalose lying along the minor groove towards the centre of the duplex. Many drug-DNA nuclear Overhauser enhancements (NOEs) in the minor groove are indicative of hydrophobic interactions over the TGCA sequence. Steric occlusion prevents a second nogalamycin molecule from binding at the symmetry-related 5'-CpA site, leading to the conclusion that the observed binding orientation in this complex is the preferred orientation free of the complication of end-effects (drug molecules occupy terminal intercalation sites in all X-ray structures) or steric interactions between drug molecules (other NMR structures have two drug molecules bound in close proximity), as previously suggested. Fluctuations in key structural parameters such as rise, helical twist, slide, shift, buckle and sugar pucker have been examined from an analysis of the final 500 ps of a 1 ns rMD simulation, and reveal that many sequence-dependent structural features previously identified by comparison of different X-ray structures lie within the range of dynamic fluctuations observed in the MD simulations. Water density calculations on MD simulation data reveal a time-averaged pattern of hydration in both the major and minor groove, in good agreement with the extensive hydration observed in two related X-ray structures in which nogalamycin is bound at terminal 5'-TpG sites. However, the pattern of hydration determined from the sign and magnitude of NOE and ROE cross-peaks to water identified in 2D NOESY and ROESY experiments identifies only a few "bound" water molecules with long residence times. These solvate the charged bicycloaminoglucose sugar ring, suggesting an important role for water molecules in mediating drug-DNA electrostatic interactions within the major groove. The high density of water molecules found in the minor groove in X-ray structures and MD simulations is found to be associated with only weakly bound solvent in solution.     

Dual role of JNK1-mediated phosphorylation of Bcl-2 in autophagy and apoptosis regulation

Autophagy and apoptosis are fundamental cellular pathways that are both regulated by JNK-mediated Bcl-2 phosphorylation. Several years ago, JNK-mediated Bcl-2 phosphorylation was shown to interfere with its binding to proapoptotic BH3 domain-containing proteins such as Bax and recently, our laboratory demonstrated that JNK1-mediated Bcl-2 phosphorylation interferes with its binding to the proautophagy BH3 domain-containing protein Beclin 1. Here, we examined the kinetic relationship between Bcl-2 phosphorylation, Bcl-2-Beclin 1 interactions, Bcl-2-Bax interactions, and caspase 3 activation during nutrient starvation. We found that after a short period of nutrient deprivation (4 hours), a small amount of Bcl-2 phosphorylation dissociates Bcl-2 from the Bcl-2-Beclin 1 complex but not from the Bcl-2-Bax complex. After 16 hours of nutrient deprivation, Bcl-2 phosphorylation reaches maximal levels, the Bcl-2-Bax complex is disrupted, and active caspase 3 is detected, indicating the initiation of apoptosis. Based on this result, we propose a speculative model for understanding the interrelationship between autophagy and apoptosis regulated by JNK1-mediated Bcl-2 phosphorylation. According to this model, rapid Bcl-2 phosphorylation may occur initially to promote cell survival by disrupting the Bcl-2-Beclin 1 complex and activating autophagy. At a certain point when autophagy is no longer able to keep the cell alive, Bcl-2 phosphorylation might then serve to inactivate its antiapoptotic function.     

Interaction of taxifolin (dihydroquercetin) with dimyristoylphosphatidylcholine multilamellar liposomes

Differential scanning calorimetry was used to study the influence of the flavonoid taxifolin (dihydroquercetin) on the temperature-dependent phase transition of dimyristoylphosphatidylcholine multilamellar liposomes. Taxifolin was added to organic solution of the lipid during the procedure of liposomes preparation (addition from-within) or to a suspension of prepared liposomes (addition from-without). In the first case, liposomes contained from 2 to 50 mol% of taxifolin added from-within; in the second case, lyposomes were treated with 0.001% or 0.01% taxifolin. In both cases, the effect was similar. When the concentration of taxifolin increased, the temperature of lipid melting decreased while the width of transition considerably enlarged. Freeze-fracture electron microscopy revealed that taxifolin did not rupture multilamellar liposomes, while the formation of ripple-phase was retarded in all bilayers even when the liposomes were treated from without. This suggested the ability of taxifolin to penetrate through numerous bilayers of multilamellar liposomes.     

Dual Role of JNK1-mediated phosphorylation of Bcl-2 in autophagy and apoptosis regulation

Autophagy and apoptosis are fundamental cellular pathways that are both regulated by JNK-mediated Bcl-2 phosphorylation. Several years ago, JNK-mediated Bcl-2 phosphorylation was shown to interfere with its binding to pro-apoptotic BH3 domain-containing proteins such as Bax and recently, our laboratory demonstrated that JNK1-mediated Bcl-2 phosphorylation interferes with its binding to the pro-autophagy BH3 domain-containing protein Beclin 1. Here, we examined the kinetic relationship between Bcl-2 phosphorylation, Bcl-2-Beclin 1 interactions, Bcl-2-Bax interactions and caspase 3 activation during nutrient starvation. We found that after a short period of nutrient deprivation (4 hours), a small amount of Bcl-2 phosphorylation dissociates Bcl-2 from the Bcl-2-Beclin 1 complex but not from the Bcl-2-Bax complex. After 16 hours of nutrient deprivation, Bcl-2 phosphorylation reaches maximal levels, the Bcl-2-Bax complex is disrupted, and active caspase 3 is detected, indicating the initiation of apoptosis. Based on this result, we propose a speculative model for understanding the interrelationship between autophagy and apoptosis regulated by JNK1-mediated Bcl-2 phosphorylation. According to this model, rapid Bcl-2 phosphorylation may occur initially to promote cell survival by disrupting the Bcl-2-Beclin 1 complex and activating autophagy. At a certain point when autophagy is no longer able to keep the cell alive, Bcl-2 phosphorylation might then serve to inactivate its anti-apoptotic function.

Addendum to: Wei Y, Pattingre S, Sinha S, Bassik M, Levine B. JNK1-mediated phosphorylation of Bcl-2 regulates starvation-induced autophagy. Mol Cell 2008; 30:678-88.     

Structure and molecular dynamics simulation of archaeal prefoldin: the molecular mechanism for binding and recognition of nonnative substrate proteins

Prefoldin (PFD) is a heterohexameric molecular chaperone complex in the eukaryotic cytosol and archaea with a jellyfish-like structure containing six long coiled-coil tentacles. PFDs capture protein folding intermediates or unfolded polypeptides and transfer them to group II chaperonins for facilitated folding. Although detailed studies on the mechanisms for interaction with unfolded proteins or cooperation with chaperonins of archaeal PFD have been performed, it is still unclear how PFD captures the unfolded protein. In this study, we determined the X-ray structure of Pyrococcus horikoshii OT3 PFD (PhPFD) at 3.0 Å resolution and examined the molecular mechanism for binding and recognition of nonnative substrate proteins by molecular dynamics (MD) simulation and mutation analyses. PhPFD has a jellyfish-like structure with six long coiled-coil tentacles and a large central cavity. Each subunit has a hydrophobic groove at the distal region where an unfolded substrate protein is bound. During MD simulation at 330 K, each coiled coil was highly flexible, enabling it to widen its central cavity and capture various nonnative proteins. Docking MD simulation of PhPFD with unfolded insulin showed that the β subunit is essentially involved in substrate binding and that the α subunit modulates the shape and width of the central cavity. Analyses of mutant PhPFDs with amino acid replacement of the hydrophobic residues of the β subunit in the hydrophobic groove have shown that βIle107 has a critical role in forming the hydrophobic groove.     

The structure of Bcl-w reveals a role for the C-terminal residues in modulating biological activity

Pro-survival Bcl-2-related proteins, critical regulators of apoptosis, contain a hydrophobic groove targeted for binding by the BH3 domain of the pro-apoptotic BH3-only proteins. The solution structure of the pro-survival protein Bcl-w, presented here, reveals that the binding groove is not freely accessible as predicted by previous structures of pro-survival Bcl-2-like molecules. Unexpectedly, the groove appears to be occluded by the C-terminal residues. Binding and kinetic data suggest that the C-terminal residues of Bcl-w and Bcl-x(L) modulate pro-survival activity by regulating ligand access to the groove. Binding of the BH3-only proteins, critical for cell death initiation, is likely to displace the hydrophobic C-terminal region of Bcl-w and Bcl-x(L). Moreover, Bcl-w does not act only by sequestering the BH3-only proteins. There fore, pro-survival Bcl-2-like molecules probably control the activation of downstream effectors by a mechanism that remains to be elucidated.     

Hydrophobic Interactions Are a Key to MDM2 Inhibition by Polyphenols as Revealed by Molecular Dynamics Simulations and MM/PBSA Free Energy Calculations

p53, a tumor suppressor protein, has been proven to regulate the cell cycle, apoptosis, and DNA repair to prevent malignant transformation. MDM2 regulates activity of p53 and inhibits its binding to DNA. In the present study, we elucidated the MDM2 inhibition potential of polyphenols (Apigenin, Fisetin, Galangin and Luteolin) by MD simulation and MM/PBSA free energy calculations. All polyphenols bind to hydrophobic groove of MDM2 and the binding was found to be stable throughout MD simulation. Luteolin showed the highest negative binding free energy value of -173.80 kJ/mol followed by Fisetin with value of -172.25 kJ/mol. It was found by free energy calculations, that hydrophobic interactions (vdW energy) have major contribution in binding free energy.     

Natural polyphenolic inhibitors against the antiapoptotic BCL-2

The apoptotic mechanism is regulated by the BCL-2 family of proteins, such as BCL-2 or Bcl-xL, which block apoptosis while Bad, Bak, Bax, Bid, Bim or Hrk induce apoptosis. The overexpression of BCL-2 was found to be related to the progression of cancer and also providing resistance towards chemotherapeutic treatments. In the present study, we found that all polyphenols (apigenin, fisetin, galangin and luteolin) bind to the hydrophobic groove of BCL-2 and the interaction is stable throughout MD simulation run. Luteolin was found to bind with highest negative binding energy and thus, claimed highest potency towards BCL-2 inhibition followed by fisetin. The hydrophobic interactions were found to be critical for stable complex formation as revealed by the vdW energy and ligplot analysis. Finally, on the basis of data obtained during the study, it can be concluded that these polyphenols have the potential to be used as lead molecules for BCL-2 inhibition.     

Effect of taxifolin on cyclophosphamide-induced oxidative and inflammatory bladder injury in rats

The role of oxidative stress and inflammation in the pathogenesis of cyclophosphamide-related side effects has been demonstrated in previous studies. This study aimed to investigate the effect of taxifolin, due to its antioxidant and anti-inflammatory properties, on cyclophosphamide-induced oxidative and inflammatory bladder injury in albino Wistar rats. The taxifolin+cyclophosphamide (TCYC) group was given 50 mg/kg of taxifolin orally by gavage. Normal saline was used as a solvent for the cyclophosphamide (CYC) group and the healthy control (HC) group. One hour after taxifolin administration, 75 mg/kg of cyclophosphamide was intraperitoneally injected in the TCYC and CYC groups. This procedure was repeated once a day for 30 days. At the end of this period, biochemical markers were studied in the excised bladder tissues and histopathological evaluations were conducted. In the histopathological evaluation of the CYC group, severe epithelial irregularity, dilatation, congestion, and polymorphonuclear leukocyte accumulation in the vascular structures were observed. Additionally, the malondialdehyde (MDA), tumor necrosis factor-α (TNF-α), interleukin-1β (IL-1β), and interleukin-6 (IL-6) levels, the total oxidant status (TOS), and the oxidative stress index (OSI) values were significantly higher, and the total glutathione (tGSH) levels and total antioxidant status (TAS) were significantly lower in the CYC group in comparison to the HC group (P<0.001). Taxifolin reduced the cyclophosphamide-induced increases in the MDA, TNF-α, IL-1β, and IL-6 levels and the TOS and OSI values; it decreased the tGSH and TAS levels and reduced histopathological damage (P<0.001). Taxifolin may be useful in the treatment of cyclophosphamide-induced bladder damage.     

Structure-based discovery of a new class of Bcl-xL antagonists

Apoptosis, or programmed cell death, plays a key role in normal tissue homeostasis ensuring a proper balance between cell production and cell loss. Anti-apoptotic Bcl-2-family proteins are central regulators of the apoptotic pathway and due to their ability to confer tumor resistance to chemotherapy or radiation, have been recently validated as targets for cancer drug discovery. Since the crucial interaction between pro- and anti-apoptotic members occurs via a conserved region located on the surface of the protein, a viable way to inhibit the anti-death activity of Bcl-2 proteins is to design small molecule inhibitors that occupy this cavity. Here, we describe a structure-based approach that led to the identification of four small molecule inhibitors directed at the hydrophobic groove on the surface of the Bcl-2 family protein Bcl-xL. The compounds were characterized in a number of assays including in vitro binding using 15N-labeled protein, a displacement DELFIA assay, and a cell-based viability assay with human cancer cells.     

A surface groove essential for viral Bcl-2 function during chronic infection in vivo

Antiapoptotic Bcl-2 family proteins inhibit apoptosis in cultured cells by binding BH3 domains of proapoptotic Bcl-2 family members via a hydrophobic BH3 binding groove on the protein surface. We investigated the physiological importance of the BH3 binding groove of an antiapoptotic Bcl-2 protein in mammals in vivo by analyzing a viral Bcl-2 family protein. We show that the gamma-herpesvirus 68 (gammaHV68) Bcl-2 family protein (gammaHV68 v-Bcl-2), which is known to inhibit apoptosis in cultured cells, inhibits both apoptosis in primary lymphocytes and Bax toxicity in yeast. Nuclear magnetic resonance determination of the gammaHV68 v-Bcl-2 structure revealed a BH3 binding groove that binds BH3 domain peptides from proapoptotic Bcl-2 family members Bax and Bak via a molecular mechanism shared with host Bcl-2 family proteins, involving a conserved arginine in the BH3 peptide binding groove. Mutations of this conserved arginine and two adjacent amino acids to alanine (SGR to AAA) within the BH3 binding groove resulted in a properly folded protein that lacked the capacity of the wild-type gammaHV68 v-Bcl-2 to bind Bax BH3 peptide and to block Bax toxicity in yeast. We tested the physiological importance of this v-Bcl-2 domain during viral infection by engineering viral mutants encoding a v-Bcl-2 containing the SGR to AAA mutation. This mutation resulted in a virus defective for both efficient reactivation of gammaHV68 from latency and efficient persistent gammaHV68 replication. These studies demonstrate an essential functional role for amino acids in the BH3 peptide binding groove of a viral Bcl-2 family member during chronic infection.     

Molecular dynamics investigations of structural and functional changes in Bcl-2 induced by the novel antagonist BDA-366

Apoptosis is a fundamental biological phenomenon, in which anti- or proapoptotic proteins of the Bcl-2 family regulate a committed step. Overexpression of Bcl-2, the prototypical antiapoptotic protein in this family, is associated with therapy resistance in various human cancers. Accordingly, Bcl-2 inhibitors intended for cancer therapy have been developed, typically against the BH3 domain. Recent experimental evidences have shown that the antiapoptotic function of Bcl-2 is not immutable, and that BDA-366, a novel antagonist of the BH4 domain, converts Bcl-2 from a survival molecule to an inducer of cell death. In this study, the underlying mechanisms of this functional conversion were investigated by accelerated molecular dynamics simulation. Results revealed that Pro127 and Trp30 in the BH4 domain rotate to stabilize BDA-366 via π-π interactions, and trigger a series of significant conformational changes of the α3 helix. This rearrangement blocks the hydrophobic binding site (HBS) in the BH3 domain and further prevents binding of BH3-only proteins, which consequently allows the BH3-only proteins to activate the proapoptotic proteins. Analysis of binding free energy confirmed that BDA-366 cross-inhibits BH3-only proteins, implying negative cooperative effects across separate binding sites. The newly identified blocked conformation of the HBS along with the open to closed transition pathway revealed by this study advances the understanding of the Bcl-2 transition from antiapoptotic to proapoptotic function, and yielded new structural insights for novel drug design against the BH4 domain.

Communicated by Ramaswamy H. Sarma     

Bax dimerizes via a symmetric BH3:groove interface during apoptosis

During apoptotic cell death, Bax and Bak change conformation and homo-oligomerize to permeabilize mitochondria. We recently reported that Bak homodimerizes via an interaction between the BH3 domain and hydrophobic surface groove, that this BH3:groove interaction is symmetric, and that symmetric dimers can be linked via the α6-helices to form the high order oligomers thought responsible for pore formation. We now show that Bax also dimerizes via a BH3:groove interaction after apoptotic signaling in cells and in mitochondrial fractions. BH3:groove dimers of Bax were symmetric as dimers but not higher order oligomers could be linked by cysteine residues placed in both the BH3 and groove. The BH3:groove interaction was evident in the majority of mitochondrial Bax after apoptotic signaling, and correlated strongly with cytochrome c release, supporting its central role in Bax function. A second interface between the Bax α6-helices was implicated by cysteine linkage studies, and could link dimers to higher order oligomers. We also found that a population of Bax:Bak heterodimers generated during apoptosis formed via a BH3:groove interaction, further demonstrating that Bax and Bak oligomerize via similar mechanisms. These findings highlight the importance of BH3:groove interactions in apoptosis regulation by the Bcl-2 protein family.     

Structural insights into the effect of hydration and ions on A-tract DNA: a molecular dynamics study

DNA structure is known to be sensitive to hydration and ionic environment. To explore the dynamics, hydration, and ion binding features of A-tract sequences, a 7-ns Molecular dynamics (MD) study has been performed on the dodecamer d(CGCAAATTTGCG)(2). The results suggest that the intrusion of Na(+) ion into the minor groove is a rare event and the structure of this dodecamer is not very sensitive to the location of the sodium ions. The prolonged MD simulation successfully leads to the formation of sequence dependent hydration patterns in the minor groove, often called spine of hydration near the A-rich region and ribbon of hydration near the GC regions. Such sequence dependent differences in the hydration patterns have been seen earlier in the high resolution crystal structure of the Drew-Dickerson sequence, but not reported for the medium resolution structures (2.0 approximately 3.0 A). Several water molecules are also seen in the major groove of the MD simulated structure, though they are not highly ordered over the extended MD. The characteristic narrowing of the minor groove in the A-tract region is seen to precede the formation of the spine of hydration. Finally, the occurrence of cross-strand C2-H2.O2 hydrogen bonds in the minor groove of A-tract sequences is confirmed. These are found to occur even before the narrowing of the minor groove, indicating that such interactions are an intrinsic feature of A-tract sequences.     

Molecular dynamics simulation and steered molecular dynamics simulation on irisin dimers

Irisin is found closely associated with promoting the browning of beige fat cells in white adipose tissue. The crystal structure reveals that irisin forms a continuous inter-subunit β-sheet dimer. Here, molecular dynamics (MD) simulation and steered molecular dynamics (SMD) simulation were performed to investigate the dissociation process and the intricate interactions between the two irisin monomers. In the process of MD, the interactions between the monomers were roughly analyzed through the average numbers of both hydrophobic contacts and H-bonds. Then, SMD was performed to investigate the accurate interaction energy between the monomers. By the analysis of dissociation energy, the van der Waals (vdW) force was identified as the major energy to maintain the dimer structure, which also verified the results of MD simulation. Meanwhile, 11 essential residues were discovered by the magnitude of rupture force during dissociation. Among them, residues Arg75, Glu79, Ile77, Ala88, and Trp90 were reported in a previous study using the method of mutagenesis and size exclusion chromatography, and several new important residues (Arg72, Leu74, Phe76, Gln78, Val80, and Asp91) were also identified. Interestingly, the new important residues that we discovered and the important residues that were reported are located in the opposite side of the β-sheet of the dimer.     

Simulation of the substrate cavity dynamics of quercetinase

Molecular dynamics (MD) simulations have been performed on quercetin 2,3 dioxygenase (2,3QD) to study the mobility and flexibility of the substrate cavity. 2,3QD is the only firmly established Cu-containing dioxygenase known so far. It catalyses the breakage of the O-heterocycle of flavonols. The substrates occupy a shallow and overall hydrophobic cavity proximal to the metal centre of the homo-dimeric enzyme. The linker connecting the C-terminal and N-terminal domains in the monomer is partly disordered in the crystal structure and part of it forms a flexible lid at the entrance of the substrate cavity. This loop has been tentatively assigned a role in the enzyme mechanism: it helps lock the substrate into place. The dynamics of this loop has been investigated by MD simulation. The initial coordinates were taken from the crystal structure of 2,3QD in the presence of the substrate kaempferol (KMP). After equilibration and simulation over 7.2ns the substrate was removed and another equilibration and simulation of 7.2ns was performed. The results show that the structures of the free enzyme as well as of the enzyme-substrate complex are stable in MD simulation. The linker shows strongly enhanced mobility in the loop region that is close to the entrance to the substrate cavity (residues 154-169). Movement of the loop takes place on a timescale of 5-10ns. To confirm the conclusions about the loop dynamics drawn from the 7.2ns simulation, the simulation was extended with another 8ns. When substrate binds into the cavity the loop orders remarkably, although mobility is retained by residues 155-158. Some regions of the loop (residues 154-160 and 164-176) move over a considerable distance and approach the substrate closely, reinforcing the idea that they lock the substrate in the substrate cavity. The enthalpic component of the interaction of the loop with the protein and the KMP appears to favour the locking of the substrate. Two water molecules were found immobilised in the cavity, one of which exhibited rotation on the picosecond timescale. When the substrate is removed, the empty cavity fills up with water within 200ps.     

Molecular dynamics study of small molecule inhibitors of the Bcl-2 family

We carried out docking and molecular dynamics simulations on ABT-737 and obatoclax, which are inhibitors of the Bcl-2 family of proteins. We modeled the binding mode of ABT-737 with Bcl-x(L) , Bcl-2, and Mcl-1 and examined their dynamical behavior. We found that the binding of the chlorobiphenyl end of ABT-737 was quite stable across all three proteins. However, the phenylpiperazine linker group was dramatically more mobile in Mcl-1 compared to either Bcl-x(L) or Bcl-2. The S-phenyl group at the p4 binding site was well-anchored in Bcl-x(L) and Bcl-2 but was somewhat more mobile in Mcl-1 although the phenyl ring itself on average stayed close to the p4 binding site in Mcl-1. This greater mobility is likely due to the greater openness of the p3 and p4 binding sites on Mcl-1. The calculated binding free energies were consistent with the much weaker binding affinity of ABT-737 for Mcl-1. Obatoclax was predicted to bind at the p1 and p2 binding sites of Mcl-1 and the binding mode was quite stable during the molecular dynamics simulation with Mcl-1 wrapping around the molecule. The modeled binding mode suggests that obatoclax is able to inhibit all three proteins because it makes use of the p1 and p2 binding sites alone, which is a fairly narrow groove in all three proteins unlike the p4 binding site, which is much broader in Mcl-1.