Recessively Inherited Deficiency of Secreted WFDC2 (HE4) Causes Nasal Polyposis and Bronchiectasis.

Rationale: Bronchiectasis is a pathological dilatation of the bronchi in the respiratory airways associated with environmental or genetic causes (e.g., cystic fibrosis, primary ciliary dyskinesia, and primary immunodeficiency disorders), but most cases remain idiopathic. Objectives: To identify novel genetic defects in unsolved cases of bronchiectasis presenting with severe rhinosinusitis, nasal polyposis, and pulmonary Pseudomonas aeruginosa infection. Methods: DNA was analyzed by next-generation or targeted Sanger sequencing. RNA was analyzed by quantitative PCR and single-cell RNA sequencing. Patient-derived cells, cell cultures, and secretions (mucus, saliva, seminal fluid) were analyzed by Western blotting and immunofluorescence microscopy, and mucociliary activity was measured. Blood serum was analyzed by electrochemiluminescence immunoassay. Protein structure and proteomic analyses were used to assess the impact of a disease-causing founder variant. Measurements and Main Results: We identified biallelic pathogenic variants in WAP four-disulfide core domain 2 (WFDC2) in 11 individuals from 10 unrelated families originating from the United States, Europe, Asia, and Africa. Expression of WFDC2 was detected predominantly in secretory cells of control airway epithelium and also in submucosal glands. We demonstrate that WFDC2 is below the limit of detection in blood serum and hardly detectable in samples of saliva, seminal fluid, and airway surface liquid from WFDC2-deficient individuals. Computer simulations and deglycosylation assays indicate that the disease-causing founder variant p.Cys49Arg structurally hampers glycosylation and, thus, secretion of mature WFDC2. Conclusions: WFDC2 dysfunction defines a novel molecular etiology of bronchiectasis characterized by the deficiency of a secreted component of the airways. A commercially available blood test combined with genetic testing allows its diagnosis.

Parkin Regulates Mitosis and Genomic Stability through Cdc20/Cdh1.

Mutations in the E3 ubiquitin ligase Parkin have been linked to familial Parkinson's disease. Parkin has also been implicated in mitosis through mechanisms that are unclear. Here we show that Parkin interacts with anaphase promoting complex/cyclosome (APC/C) coactivators Cdc20 and Cdh1 to mediate the degradation of several key mitotic regulators independent of APC/C. We demonstrate that ordered progression through mitosis is orchestrated by two distinct E3 ligases through the shared use of Cdc20 and Cdh1. Furthermore, Parkin is phosphorylated and activated by polo-like kinase 1 (Plk1) during mitosis. Parkin deficiency results in overexpression of its substrates, mitotic defects, genomic instability, and tumorigenesis. These results suggest that the Parkin-Cdc20/Cdh1 complex is an important regulator of mitosis.

Plasma butyrylcholinesterase regulates ghrelin to control aggression.

Ongoing mouse studies of a proposed therapy for cocaine abuse based on viral gene transfer of butyrylcholinesterase (BChE) mutated for accelerated cocaine hydrolysis have yielded surprising effects on aggression. Further investigation has linked these effects to a reduction in circulating ghrelin, driven by BChE at levels ∼ 100-fold above normal. Tests with human BChE showed ready ghrelin hydrolysis at physiologic concentrations, and multiple low-mass molecular dynamics simulations revealed that ghrelin's first five residues fit sterically and electrostatically into BChE's active site. Consistent with in vitro results, male BALB/c mice with high plasma BChE after gene transfer exhibited sharply reduced plasma ghrelin. Unexpectedly, such animals fought less, both spontaneously and in a resident/intruder provocation model. One mutant BChE was found to be deficient in ghrelin hydrolysis. BALB/c mice transduced with this variant retained normal plasma ghrelin levels and did not differ from untreated controls in the aggression model. In contrast, C57BL/6 mice with BChE gene deletion exhibited increased ghrelin and fought more readily than wild-type animals. Collectively, these findings indicate that BChE-catalyzed ghrelin hydrolysis influences mouse aggression and social stress, with potential implications for humans.

How fast fast-folding proteins fold in silico.

In reported microcanonical molecular dynamics simulations, fast-folding proteins CLN025 and Trp-cage autonomously folded to experimentally determined native conformations. However, the folding times of these proteins derived from the simulations were more than 4-10 times longer than their experimental values. This article reports autonomous folding of CLN025 and Trp-cage in isobaric-isothermal molecular dynamics simulations with agreements within factors of 0.69-1.75 between simulated and experimental folding times at different temperatures. These results show that CLN025 and Trp-cage can now autonomously fold in silico as fast as in experiments, and suggest that the accuracy of folding simulations for fast-folding proteins begins to overlap with the accuracy of folding experiments. This opens new prospects of developing computer algorithms that can predict both ensembles of conformations and their interconversion rates for a protein from its sequence for artificial intelligence on how and when a protein acts as a receiver, switch, and relay to facilitate various subcellular-to-tissue communications. Then the genetic information that encodes proteins can be better read in the context of intricate biological functions.

On the use of the experimentally determined enzyme inhibition constant as a measure of absolute binding affinity.

Defined as a state function representing an inhibitor's absolute affinity for its target enzyme, the experimentally determined enzyme inhibition constant (Ki) is widely used to rank order binding affinities of different inhibitors for a common enzyme or different enzymes for a common inhibitor and to benchmark computational approaches to predicting binding affinity. Herein, we report that adsorption of bis(7)-tacrine to the glass container surface increased its Ki against Electrophorus electricus acetylcholinesterase (eeAChE) to 3.2 ± 0.1 nM (n = 5) compared to 2.9 ± 0.4 pM (n = 5) that was determined using plastic containers with other assay conditions kept the same. We also report that, due to binding or "adsorption" of bis(7)-tacrine to the inactive eeAChE, the bis(7)-tacrine Ki increased from 2.9 ± 0.4 pM (n = 5) to 734 ± 70 pM (n = 5) as the specific eeAChE activity decreased from 342 U/mg to 26 U/mg while other assay conditions were kept the same. These results caution against using Kis to rank order binding potencies, define selectivity, or benchmark computational methods without knowing detailed assay conditions.

Prime, Shock, and Kill: Priming CD4 T Cells from HIV Patients with a BCL-2 Antagonist before HIV Reactivation Reduces HIV Reservoir Size.

UNLABELLED: Understanding how some HIV-infected cells resist the cytotoxicity of HIV replication is crucial to enabling HIV cure efforts. HIV killing of CD4 T cells that replicate HIV can involve HIV protease-mediated cleavage of procaspase 8 to generate a fragment (Casp8p41) that directly binds and activates the mitochondrial proapoptotic protein BAK. Here, we demonstrate that Casp8p41 also binds with nanomolar affinity to the antiapoptotic protein Bcl-2, which sequesters Casp8p41 and prevents apoptosis. Further, we show that central memory CD4 T cells (TCM) from HIV-infected individuals have heightened expression of BCL-2 relative to procaspase 8, possibly explaining the persistence of HIV-infected TCMdespite generation of Casp8p41. Consistent with this hypothesis, the selective BCL-2 antagonist venetoclax induced minimal killing of uninfected CD4 T cells but markedly increased the death of CD4 T cells and diminished cell-associated HIV DNA when CD4 T cells from antiretroviral therapy (ART)-suppressed HIV patients were induced with αCD3/αCD28 to reactivate HIVex vivo Thus, priming CD4 T cells from ART suppressed HIV patients with a BCL-2 antagonist, followed by HIV reactivation, achieves reductions in cell-associated HIV DNA, whereas HIV reactivation alone does not. IMPORTANCE: HIV infection is incurable due to a long-lived reservoir of HIV(+)memory CD4 T cells, and no clinically relevant interventions have been identified that reduce the number of these HIV DNA-containing cells. Since postintegration HIV replication can result in HIV protease generation of Casp8p41, which activates BAK, causing infected CD4 T cell death, we sought to determine whether this occurs in memory CD4 T cells. Here, we demonstrate that memory CD4 T cells can generate Casp8p41 and yet are intrinsically resistant to death induced by diverse stimuli, including Casp8p41. Furthermore, BCL-2 expression is relatively increased in these cells and directly binds and inhibits Casp8p41's proapoptotic effects. Antagonizing BCL-2 with venetoclax derepresses this antagonism, resulting in death, preferentially in HIV DNA containing cells, since only these cells generate Casp8p41. Thus, BCL-2 antagonism is a clinically relevant intervention with the potential to reduce HIV reservoir size in patients.

FF12MC: A revised AMBER forcefield and new protein simulation protocol.

Specialized to simulate proteins in molecular dynamics (MD) simulations with explicit solvation, FF12MC is a combination of a new protein simulation protocol employing uniformly reduced atomic masses by tenfold and a revised AMBER forcefield FF99 with (i) shortened CH bonds, (ii) removal of torsions involving a nonperipheral sp(3) atom, and (iii) reduced 1-4 interaction scaling factors of torsions ϕ and ψ. This article reports that in multiple, distinct, independent, unrestricted, unbiased, isobaric-isothermal, and classical MD simulations FF12MC can (i) simulate the experimentally observed flipping between left- and right-handed configurations for C14-C38 of BPTI in solution, (ii) autonomously fold chignolin, CLN025, and Trp-cage with folding times that agree with the experimental values, (iii) simulate subsequent unfolding and refolding of these miniproteins, and (iv) achieve a robust Z score of 1.33 for refining protein models TMR01, TMR04, and TMR07. By comparison, the latest general-purpose AMBER forcefield FF14SB locks the C14-C38 bond to the right-handed configuration in solution under the same protein simulation conditions. Statistical survival analysis shows that FF12MC folds chignolin and CLN025 in isobaric-isothermal MD simulations 2-4 times faster than FF14SB under the same protein simulation conditions. These results suggest that FF12MC may be used for protein simulations to study kinetics and thermodynamics of miniprotein folding as well as protein structure and dynamics. Proteins 2016; 84:1490-1516. © 2016 The Authors Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.

Evaluation of the BH3-only protein Puma as a direct Bak activator.

Interactions among Bcl-2 family proteins play critical roles in cellular life and death decisions. Previous studies have established the BH3-only proteins Bim, tBid, and Noxa as "direct activators" that are able to directly initiate the oligomerization and activation of Bak and/or Bax. Earlier studies of Puma have yielded equivocal results, with some concluding that it also acts as a direct activator and other studies suggesting that it acts solely as a sensitizer BH3-only protein. In the present study we examined the interaction of Puma BH3 domain or full-length protein with Bak by surface plasmon resonance, assessed Bak oligomerization status by cross-linking followed by immunoblotting, evaluated the ability of the Puma BH3 domain to induce Bak-mediated permeabilization of liposomes and mitochondria, and determined the effect of wild type and mutant Puma on cell viability in a variety of cellular contexts. Results of this analysis demonstrate high affinity (KD = 26 ± 5 nM) binding of the Puma BH3 domain to purified Bak ex vivo, leading to Bak homo-oligomerization and membrane permeabilization. Mutations in Puma that inhibit (L141E/M144E/L148E) or enhance (M144I/A145G) Puma BH3 binding to Bak also produce corresponding alterations in Bak oligomerization, Bak-mediated membrane permeabilization and, in a cellular context, Bak-mediated killing. Collectively, these results provide strong evidence that Puma, like Bim, Noxa, and tBid, is able to act as a direct Bak activator.

At least 10% shorter C-H bonds in cryogenic protein crystal structures than in current AMBER forcefields.

High resolution protein crystal structures resolved with X-ray diffraction data at cryogenic temperature are commonly used as experimental data to refine forcefields and evaluate protein folding simulations. However, it has been unclear hitherto whether the C-H bond lengths in cryogenic protein structures are significantly different from those defined in forcefields to affect protein folding simulations. This article reports the finding that the C-H bonds in high resolution cryogenic protein structures are 10-14% shorter than those defined in current AMBER forcefields, according to 3709 C-H bonds in the cryogenic protein structures with resolutions of 0.62-0.79 Å. Also, 20 all-atom, isothermal-isobaric, 0.5-µs molecular dynamics simulations showed that chignolin folded from a fully-extended backbone formation to the native ß-hairpin conformation in the simulations using AMBER forcefield FF12SB at 300 K with an aggregated native state population including standard error of 10 ± 4%. However, the aggregated native state population with standard error reduced to 3 ± 2% in the same simulations except that C-H bonds were shortened by 10-14%. Furthermore, the aggregated native state populations with standard errors increased to 35 ± 3% and 26 ± 3% when using FF12MC, which is based on AMBER forcefield FF99, with and without the shortened C-H bonds, respectively. These results show that the 10-14% bond length differences can significantly affect protein folding simulations and suggest that re-parameterization of C-H bonds according to the cryogenic structures could improve the ability of a forcefield to fold proteins in molecular dynamics simulations.

Use of 1-4 interaction scaling factors to control the conformational equilibrium between α-helix and ß-strand.

1-4 interaction scaling factors are used in AMBER forcefields to reduce the exaggeration of short-range repulsion caused by the 6-12 Lennard-Jones potential and a nonpolarizable charge model and to obtain better agreements of small-molecule conformational energies with experimental data. However, the effects of these scaling factors on protein secondary structure conformations have not been investigated until now. This article reports the finding that the 1-4 interactions among the protein backbone atoms separated by three consecutive covalent bonds are more repulsive in the α-helix conformation than in two ß-strand conformations. Therefore, the 1-4 interaction scaling factors of protein backbone torsions ϕ and ψ control the conformational equilibrium between α-helix and ß-strand. Molecular dynamics simulations confirm that reducing the ϕ and ψ scaling factors readily converts the α-helix conformation of AcO-(AAQAA)3-NH2 to a ß-strand conformation, and the reverse occurs when these scaling factors are increased. These results suggest that the ϕ and ψ scaling factors can be used to generate the α-helix or ß-strand conformation in situ and to control the propensities of a forcefield for adopting secondary structure elements.

Regulation of error-prone translesion synthesis by Spartan/C1orf124.

Translesion synthesis (TLS) employs low fidelity polymerases to replicate past damaged DNA in a potentially error-prone process. Regulatory mechanisms that prevent TLS-associated mutagenesis are unknown; however, our recent studies suggest that the PCNA-binding protein Spartan plays a role in suppression of damage-induced mutagenesis. Here, we show that Spartan negatively regulates error-prone TLS that is dependent on POLD3, the accessory subunit of the replicative DNA polymerase Pol δ. We demonstrate that the putative zinc metalloprotease domain SprT in Spartan directly interacts with POLD3 and contributes to suppression of damage-induced mutagenesis. Depletion of Spartan induces complex formation of POLD3 with Rev1 and the error-prone TLS polymerase Pol ζ, and elevates mutagenesis that relies on POLD3, Rev1 and Pol ζ. These results suggest that Spartan negatively regulates POLD3 function in Rev1/Pol ζ-dependent TLS, revealing a previously unrecognized regulatory step in error-prone TLS.

Low-mass molecular dynamics simulation: a simple and generic technique to enhance configurational sampling.

CLN025 is one of the smallest fast-folding proteins. Until now it has not been reported that CLN025 can autonomously fold to its native conformation in a classical, all-atom, and isothermal-isobaric molecular dynamics (MD) simulation. This article reports the autonomous and repeated folding of CLN025 from a fully extended backbone conformation to its native conformation in explicit solvent in multiple 500-ns MD simulations at 277K and 1atm with the first folding event occurring as early as 66.1ns. These simulations were accomplished by using AMBER forcefield derivatives with atomic masses reduced by 10-fold on Apple Mac Pros. By contrast, no folding event was observed when the simulations were repeated using the original AMBER forcefields of FF12SB and FF14SB. The results demonstrate that low-mass MD simulation is a simple and generic technique to enhance configurational sampling. This technique may propel autonomous folding of a wide range of miniature proteins in classical, all-atom, and isothermal-isobaric MD simulations performed on commodity computers-an important step forward in quantitative biology.

ABCC9 mutations identified in human dilated cardiomyopathy disrupt catalytic KATP channel gating.

Stress tolerance of the heart requires high-fidelity metabolic sensing by ATP-sensitive potassium (K(ATP)) channels that adjust membrane potential-dependent functions to match cellular energetic demand. Scanning of genomic DNA from individuals with heart failure and rhythm disturbances due to idiopathic dilated cardiomyopathy identified two mutations in ABCC9, which encodes the regulatory SUR2A subunit of the cardiac K(ATP) channel. These missense and frameshift mutations mapped to evolutionarily conserved domains adjacent to the catalytic ATPase pocket within SUR2A. Mutant SUR2A proteins showed aberrant redistribution of conformations in the intrinsic ATP hydrolytic cycle, translating into abnormal K(ATP) channel phenotypes with compromised metabolic signal decoding. Defective catalysis-mediated pore regulation is thus a mechanism for channel dysfunction and susceptibility to dilated cardiomyopathy.

A BAK subdomain that binds mitochondrial lipids selectively and releases cytochrome C.

How BAK and BAX induce mitochondrial outer membrane (MOM) permeabilization (MOMP) during apoptosis is incompletely understood. Here we have used molecular dynamics simulations, surface plasmon resonance, and assays for membrane permeabilization in vitro and in vivo to assess the structure and function of selected BAK subdomains and their derivatives. Results of these studies demonstrate that BAK helical regions α5 and α6 bind the MOM lipid cardiolipin. While individual peptides corresponding to these helical regions lack the full biological activity of BAK, tandem peptides corresponding to α4-α5, α5-α6, or α6-α7/8 can localize exogenous proteins to mitochondria, permeabilize liposomes composed of MOM lipids, and cause MOMP in the absence of the remainder of the BAK protein. Importantly, the ability of these tandem helices to induce MOMP under cell-free conditions is diminished by mutations that disrupt the U-shaped helix-turn-helix structure of the tandem peptides or decrease their lipid binding. Likewise, BAK-induced apoptosis in intact cells is diminished by CLS1 gene interruption, which decreases mitochondrial cardiolipin content, or by BAK mutations that disrupt the U-shaped tandem peptide structure or diminish lipid binding. Collectively, these results suggest that BAK structural rearrangements during apoptosis might mobilize helices involved in specific protein-lipid interactions that are critical for MOMP.

Functional and Clinical Characterization of Variants of Uncertain Significance Identifies a Hotspot for Inactivating Missense Variants in RAD51C.

Pathogenic protein-truncating variants of RAD51C, which plays an integral role in promoting DNA damage repair, increase the risk of breast and ovarian cancer. A large number of RAD51C missense variants of uncertain significance (VUS) have been identified, but the effects of the majority of these variants on RAD51C function and cancer predisposition have not been established. Here, analysis of 173 missense variants by a homology-directed repair (HDR) assay in reconstituted RAD51C-/- cells identified 30 nonfunctional (deleterious) variants, including 18 in a hotspot within the ATP-binding region. The deleterious variants conferred sensitivity to cisplatin and olaparib and disrupted formation of RAD51C/XRCC3 and RAD51B/RAD51C/RAD51D/XRCC2 complexes. Computational analysis indicated the deleterious variant effects were consistent with structural effects on ATP-binding to RAD51C. A subset of the variants displayed similar effects on RAD51C activity in reconstituted human RAD51C-depleted cancer cells. Case-control association studies of deleterious variants in women with breast and ovarian cancer and noncancer controls showed associations with moderate breast cancer risk [OR, 3.92; 95% confidence interval (95% CI), 2.18-7.59] and high ovarian cancer risk (OR, 14.8; 95% CI, 7.71-30.36), similar to protein-truncating variants. This functional data supports the clinical classification of inactivating RAD51C missense variants as pathogenic or likely pathogenic, which may improve the clinical management of variant carriers. SIGNIFICANCE: Functional analysis of the impact of a large number of missense variants on RAD51C function provides insight into RAD51C activity and information for classification of the cancer relevance of RAD51C variants.

How neocarcerand Octacid4 self-assembles with guests into irreversible noncovalent complexes and what accelerates the assembly.

Cram's supramolecular capsule Octacid4 can irreversibly and noncovalently self-assemble with small-molecule guests at room temperature, but how they self-assemble and what accelerates their assembly remain poorly understood. This article reports 81 distinct Octacid4•guest self-assembly pathways captured in unrestricted, unbiased molecular dynamics simulations. These pathways reveal that the self-assembly was initiated by the guest interaction with the cavity portal exterior of Octacid4 to increase the portal collisions that led to the portal expansion for guest ingress, and completed by the portal contraction caused by the guest docking inside the cavity to impede guest egress. The pathways also reveal that the self-assembly was accelerated by engaging populated host and guest conformations for the exterior interaction to increase the portal collision frequency. These revelations may help explain why the presence of an exterior binding site at the rim of the enzyme active site is a fundamental feature of fast enzymes such as acetylcholinesterase and why small molecules adopt local minimum conformations when binding to proteins. Further, these revelations suggest that irreversible noncovalent complexes with fast assembly rates could be developed-by engaging populated host and guest conformations for the exterior interactions-for materials technology, data storage and processing, molecular sensing and tagging, and drug therapy.

TREM2 interacts with TDP-43 and mediates microglial neuroprotection against TDP-43-related neurodegeneration.

Triggering receptor expressed on myeloid cell 2 (TREM2) is linked to risk of neurodegenerative disease. However, the function of TREM2 in neurodegeneration is still not fully understood. Here, we investigated the role of microglial TREM2 in TAR DNA-binding protein 43 (TDP-43)-related neurodegeneration using virus-mediated and transgenic mouse models. We found that TREM2 deficiency impaired phagocytic clearance of pathological TDP-43 by microglia and enhanced neuronal damage and motor impairments. Mass cytometry analysis revealed that human TDP-43 (hTDP-43) induced a TREM2-dependent subpopulation of microglia with high CD11c expression and phagocytic ability. Using mass spectrometry (MS) and surface plasmon resonance (SPR) analysis, we further demonstrated an interaction between TDP-43 and TREM2 in vitro and in vivo as well as in human tissues from individuals with amyotrophic lateral sclerosis (ALS). We computationally identified regions within hTDP-43 that interact with TREM2. Our data highlight that TDP-43 is a possible ligand for microglial TREM2 and that this interaction mediates neuroprotection of microglia in TDP-43-related neurodegeneration.

How the water-soluble hemicarcerand incarcerates guests at room temperature decoded with modular simulations.

Molecular dynamics simulations of hemicarcerands and related variants allow the study of constrictive binding and offer insight into the rules of molecular complexation, but are limited because three-dimensional models of hemicarcerands are tedious to build and their atomic charges are complicated to derive. There have been no molecular dynamics simulations of the reported water-soluble hemicarcerand (Octacid4) that explain how Octacid4 encapsulates guests at 298 K and keeps them encapsulated at 298 K in NMR experiments. Herein we report a modular approach to hemicarcerand simulations that simplifies the model building and charge derivation in a manner reminiscent of the approach to protein simulations with truncated amino acids as building blocks. We also report that in aqueous molecular dynamics simulations at 298 K apo Octacid4 adopts two clusters of conformations one of which has an equatorial portal open but the guest-bound Octacid4 adopts one cluster of conformations with all portals closed. These results explain how Octacid4 incarcerates guests at room temperature and suggest that the guest-induced host conformational change that impedes decomplexation is a previously unrecognized conformational characteristic that promotes strong molecular complexation.

Enzymatic activation of pyruvate kinase increases cytosolic oxaloacetate to inhibit the Warburg effect.

Pharmacological activation of the glycolytic enzyme PKM2 or expression of the constitutively active PKM1 isoform in cancer cells results in decreased lactate production, a phenomenon known as the PKM2 paradox in the Warburg effect. Here we show that oxaloacetate (OAA) is a competitive inhibitor of human lactate dehydrogenase A (LDHA) and that elevated PKM2 activity increases de novo synthesis of OAA through glutaminolysis, thereby inhibiting LDHA in cancer cells. We also show that replacement of human LDHA with rabbit LDHA, which is relatively resistant to OAA inhibition, eliminated the paradoxical correlation between the elevated PKM2 activity and the decreased lactate concentration in cancer cells treated with a PKM2 activator. Furthermore, rabbit LDHA-expressing tumours, compared to human LDHA-expressing tumours in mice, displayed resistance to the PKM2 activator. These findings describe a mechanistic explanation for the PKM2 paradox by showing that OAA accumulates and inhibits LDHA following PKM2 activation.

Characterization of an alternative BAK-binding site for BH3 peptides.

Many cellular stresses are transduced into apoptotic signals through modification or up-regulation of the BH3-only subfamily of BCL2 proteins. Through direct or indirect mechanisms, these proteins activate BAK and BAX to permeabilize the mitochondrial outer membrane. While the BH3-only proteins BIM, PUMA, and tBID have been confirmed to directly activate BAK through its canonical BH3 binding groove, whether the BH3-only proteins BMF, HRK or BIK can directly activate BAK is less clear. Here we show that BMF and HRK bind and directly activate BAK. Through NMR studies, site-directed mutagenesis, and advanced molecular dynamics simulations, we also find that BAK activation by BMF and possibly HRK involves a previously unrecognized binding groove formed by BAK α4, α6, and α7 helices. Alterations in this groove decrease the ability of BMF and HRK to bind BAK, permeabilize membranes and induce apoptosis, suggesting a potential role for this BH3-binding site in BAK activation.