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.

Promising tacrine/huperzine A-based dimeric acetylcholinesterase inhibitors for neurodegenerative disorders: From relieving symptoms to modifying diseases through multitarget.

Neurodegenerative disorders, such as Alzheimer's disease and Parkinson's disease, are devastating diseases in the elderly world, which are closely associated with progressive neuronal loss induced by a variety of genetic and/or environmental factors. Unfortunately, currently available treatments for neurodegenerative disorders can only relieve the symptoms but not modify the pathological processes. Over the past decades, our group by collaborating with Profs. Yuan-Ping Pang and Paul R. Carlier has developed three series of homo/hetero dimeric acetylcholinesterase inhibitors derived from tacrine and/or huperzine A. The representative dimers bis(3)-Cognitin (B3C), bis(12)-hupyridone, and tacrine(10)-hupyridone might possess disease-modifying effects through the modulation of N-methyl-d-aspartic acid receptors, the activation of myocyte enhancer factor 2D gene transcription, and the promotion of neurotrophic factor secretion. In this review, we summarize that the representative dimers, such as B3C, provide neuroprotection against a variety of neurotoxins via multiple targets, including the inhibitions of N-methyl-d-aspartic acid receptor with pathological-activated potential, neuronal nitric oxide synthase, and ß-amyloid cascades synergistically. More importantly, B3C might offer disease-modifying potentials by activating myocyte enhancer factor 2D transcription, inducing neuritogenesis, and promoting the expressions of neurotrophic factors in vitro and in vivo. Taken together, the novel dimers might offer synergistic disease-modifying effects, proving that dimerization might serve as one of the strategies to develop new generation of therapeutics for neurodegenerative disorders.

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.

Significant combination of Aß aggregation inhibitory and neuroprotective properties in silico, in vitro and in vivo by bis(propyl)-cognitin, a multifunctional anti-Alzheimer's agent.

Inhibition of Aß aggregation and neurotoxicity has been developed as an attractive therapeutic strategy to combat Alzheimer's disease (AD). Bis(propyl)-cognitin (B3C) is a multifunctional dimer derived from tacrine. Herein, the anti-aggregation and disassembly effects of B3C on Aß, together with the neuroprotective effects and underlying mechanisms of B3C against Aß-induced neurotoxicity were investigated in silico, in vitro and in vivo. Data from Thioflavin-T fluorescence and atomic force microscopy assays indicated that B3C (1-10 µM), but not its monomer tacrine, greatly inhibited the formation of Aß fibrils and disaggregated pre-formed mature Aß fibrils. Comparative molecular dynamics simulation results revealed a possible binding mode that prevented Aß fibrils formation, showing that B3C favorably bound to Aß via hydrophobic interactions. Additionally, B3C was able to block the neurotoxicity caused by Aß fibrils in cultured PC12 cells. Very encouragingly, B3C (0.3 and 0.45 mg/kg) markedly alleviated the cognitive impairments in rats insulted by intra-hippocampal injection of Aß1-42 fibrils, more potently than tacrine (1 and 2 mg/kg). Furthermore, mechanistic studies demonstrated that B3C reversed the inhibition of phospho-GSK3ß at Ser9 site in vitro and in vivo caused by Aß, suggesting the neuroprotection of B3C was achieved through the inhibition of GSK3ß pathway. These findings indicate that B3C could serve as an effective inhibitor of Aß aggregation and neurotoxicity, and provide novel molecular insights into the potential application of B3C in AD prevention and treatment.

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.

Bis(propyl)-cognitin potentiates rehabilitation of treadmill exercise after a transient focal cerebral ischemia, possibly via inhibiting NMDA receptor and regulating VEGF expression.

Combination therapies may have greater efficacy compared with monotherapy in treating stroke. We investigated the molecular mechanisms by which the combination of bis(propyl)-cognitin, an uncompetitive antagonist of NMDA receptor, and treadmill exercise promote rehabilitation after ischemic stroke. Rats were distributed into 3 treatment groups: infarct/bis(propyl)-cognitin(drug only group, DO); infarct/treadmill exercise(exercise only group, EO); infarct/bis(propyl)-cognitin + treadmill exercise (drug + exercise group, DE). The DE group had further separated to 3 sub-groups to investigate the effects achieved by different time for drug administration (60 min before stroke (DE-60 m), 15 min (DE+15 m) and 60 min (DE+60 m) after stroke). Although all infarct groups improved over time, the combination of bis(propyl)-cognitin and treadmill exercise effectively enhanced motor recovery during 14-day intervention. Early drug intervention has a best recovery result, the DE+15 m group with drug intervention at 15-min after stroke had better motor recovery than DE+60 m, DO, EO and control groups. Both bis(propyl)-cognitin and treadmill exercise significantly elevated brain VEGF expression and decreased brain infarct volume at 14 day post-ischemia. Our study reveals that bis(propyl)-cognitin potentiated rehabilitation of treadmill exercise after ischemic stroke, possibly via regulating brain VEGF expression, indicating that the combination of NMDA receptor antagonists and exercise might be useful for stroke rehabilitation.

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.