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.

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.

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.

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 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.

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.

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.

CFAP45 deficiency causes situs abnormalities and asthenospermia by disrupting an axonemal adenine nucleotide homeostasis module.

Axonemal dynein ATPases direct ciliary and flagellar beating via adenosine triphosphate (ATP) hydrolysis. The modulatory effect of adenosine monophosphate (AMP) and adenosine diphosphate (ADP) on flagellar beating is not fully understood. Here, we describe a deficiency of cilia and flagella associated protein 45 (CFAP45) in humans and mice that presents a motile ciliopathy featuring situs inversus totalis and asthenospermia. CFAP45-deficient cilia and flagella show normal morphology and axonemal ultrastructure. Proteomic profiling links CFAP45 to an axonemal module including dynein ATPases and adenylate kinase as well as CFAP52, whose mutations cause a similar ciliopathy. CFAP45 binds AMP in vitro, consistent with structural modelling that identifies an AMP-binding interface between CFAP45 and AK8. Microtubule sliding of dyskinetic sperm from Cfap45-/- mice is rescued with the addition of either AMP or ADP with ATP, compared to ATP alone. We propose that CFAP45 supports mammalian ciliary and flagellar beating via an adenine nucleotide homeostasis module.

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.

Remote Activation of a Latent Epitope in an Autoantigen Decoded With Simulated B-Factors.

Mutants of a catalytically inactive variant of Proteinase 3 (PR3)-iPR3-Val103 possessing a Ser195Ala mutation relative to wild-type PR3-Val103-offer insights into how autoantigen PR3 interacts with antineutrophil cytoplasmic antibodies (ANCAs) in granulomatosis with polyangiitis (GPA) and whether such interactions can be interrupted. Here we report that iHm5-Val103, a triple mutant of iPR3-Val103, bound a monoclonal antibody (moANCA518) from a GPA patient on an epitope remote from the mutation sites, whereas the corresponding epitope of iPR3-Val103 was latent to moANCA518. Simulated B-factor analysis revealed that the binding of moANCA518 to iHm5-Val103 was due to increased main-chain flexibility of the latent epitope caused by remote mutations, suggesting rigidification of epitopes with therapeutics to alter pathogenic PR3·ANCA interactions as new GPA treatments.

A pilot clinical trial testing topical resiquimod and a xenopeptide as immune adjuvants for a melanoma vaccine targeting MART-1.

A vaccine that could expand melanoma-specific T cells might reduce the risk of recurrence of resected melanoma and could provide an alternative or adjunct to standard immunotherapy options. We tested the safety and immunogenicity of a vaccine coupling a melanoma-associated peptide with a xenogenic peptide (to promote epitope spreading) and/or resiquimod (to activate antigen-presenting cells). HLA-A2-positive patients with resected stage II, III, and IV melanoma were assigned to treatment on one of three schedules. All patients received three subcutaneous doses of the peptide MART-1a mixed with Montanide. In addition, patients on schedule 1 received the xenoantigen peptide Gag267-274, patients on schedule 2 received topical resiquimod, and patients on schedule 3 received both Gag267-274 and resiquimod. Blood samples were tested for the frequency of antigen-specific T cells by tetramer assay, as well as immune cell subtypes and plasma cytokine levels. Patients enrolled from October 2012 to December 2014, with 10 patients enrolling to each schedule. The most common adverse events were injection site reaction (26 patients) and fatigue (15 patients). Tetramer analysis revealed antigen-specific responses (defined as doubling of MART-1a-specific T cells from pretreatment to post-treatment) in 20, 60, and 40% of patients treated on schedules 1, 2, and 3, respectively. Vaccine treatment consisting of MART-1a peptide, Gag267-274, Montanide, and topical resiquimod was well-tolerated. The addition of the Gag267-274 xenoantigen was not associated with an increase in the response to MART-1a, whereas use of topical resiquimod was associated with a higher frequency of MART-1a-specific T-cell responses that did not meet statistical significance.

Peptide-Binding Groove Contraction Linked to the Lack of T Cell Response: Using Complex Structure and Energy To Identify Neoantigens.

Using personalized peptide vaccines (PPVs) to target tumor-specific nonself-antigens (neoantigens) is a promising approach to cancer treatment. However, the development of PPVs is hindered by the challenge of identifying tumor-specific neoantigens, in part because current in silico methods for identifying such neoantigens have limited effectiveness. In this article, we report the results of molecular dynamics simulations of 12 oligopeptides bound with an HLA, revealing a previously unrecognized association between the inability of an oligopeptide to elicit a T cell response and the contraction of the peptide-binding groove upon binding of the oligopeptide to the HLA. Our conformational analysis showed that this association was due to incompatibility at the interface between the contracted groove and its αß-T cell Ag receptor. This structural demonstration that having the capability to bind HLA does not guarantee immunogenicity prompted us to develop an atom-based method (SEFF12MC) to predict immunogenicity through using the structure and energy of a peptide·HLA complex to assess the propensity of the complex for further complexation with its TCR. In predicting the immunogenicities of the 12 oligopeptides, SEFF12MC achieved a 100% success rate, compared with success rates of 25-50% for 11 publicly available residue-based methods including NetMHC-4.0. Although further validation and refinements of SEFF12MC are required, our results suggest a need to develop in silico methods that assess peptide characteristics beyond their capability to form stable binary complexes with HLAs to help remove hurdles in using the patient tumor DNA information to develop PPVs for personalized cancer immunotherapy.

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.

Use of multiple picosecond high-mass molecular dynamics simulations to predict crystallographic B-factors of folded globular proteins.

Predicting crystallographic B-factors of a protein from a conventional molecular dynamics simulation is challenging, in part because the B-factors calculated through sampling the atomic positional fluctuations in a picosecond molecular dynamics simulation are unreliable, and the sampling of a longer simulation yields overly large root mean square deviations between calculated and experimental B-factors. This article reports improved B-factor prediction achieved by sampling the atomic positional fluctuations in multiple picosecond molecular dynamics simulations that use uniformly increased atomic masses by 100-fold to increase time resolution. Using the third immunoglobulin-binding domain of protein G, bovine pancreatic trypsin inhibitor, ubiquitin, and lysozyme as model systems, the B-factor root mean square deviations (mean ± standard error) of these proteins were 3.1 ± 0.2-9 ± 1 Å2 for Cα and 7.3 ± 0.9-9.6 ± 0.2 Å2 for Cγ, when the sampling was done for each of these proteins over 20 distinct, independent, and 50-picosecond high-mass molecular dynamics simulations with AMBER forcefield FF12MC or FF14SB. These results suggest that sampling the atomic positional fluctuations in multiple picosecond high-mass molecular dynamics simulations may be conducive to a priori prediction of crystallographic B-factors of a folded globular protein.

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.

Physiological roles for butyrylcholinesterase: A BChE-ghrelin axis.

Butyrylcholinesterase (BChE) has long been regarded as an "orphan enzyme" with no specific physiological role other than to metabolize exogenous bioactive esters in the diet or in medicines. Human beings with genetic mutations that eliminate all BChE activity appear completely normal, and BChE-knockout mice have been described as "lacking a phenotype" except for faster weight gain on high-fat diets. However, our recent studies with viral gene transfer of BChE in mice reveal that BChE hydrolyzes the so-called "hunger hormone," ghrelin, at a rate which strongly affects the circulating levels of this peptide hormone. This action has important consequences for weight gain and fat metabolism. Surprisingly, it also impacts emotional behaviors such as aggression. Overexpression of BChE leads to low ghrelin levels in the blood stream and reduces aggression and social stress in mice. Under certain circumstances these combined effects contribute to increased life-span in group-housed animals. These findings may generalize to humans, as recent clinical studies by multiple investigators indicate that, among patients with severe cardiovascular disease, longevity correlates with increasing levels of plasma BChE activity.

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.