The nervous system plays an important role in the regulation of epithelial homeostasis and has also been postulated to play a role in tumorigenesis. We provide evidence that proper innervation is critical at all stages of gastric tumorigenesis. In three separate mouse models of gastric cancer, surgical or pharmacological denervation of the stomach (bilateral or unilateral truncal vagotomy, or local injection of botulinum toxin type A) markedly reduced tumor incidence and progression, but only in the denervated portion of the stomach. Vagotomy or botulinum toxin type A treatment also enhanced the therapeutic effects of systemic chemotherapy and prolonged survival. Denervation-induced suppression of tumorigenesis was associated with inhibition of Wnt signaling and suppression of stem cell expansion. In gastric organoid cultures, neurons stimulated growth in a Wnt-mediated fashion through cholinergic signaling. Furthermore, pharmacological inhibition or genetic knockout of the muscarinic acetylcholine M3 receptor suppressed gastric tumorigenesis. In gastric cancer patients, tumor stage correlated with neural density and activated Wnt signaling, whereas vagotomy reduced the risk of gastric cancer. Together, our findings suggest that vagal innervation contributes to gastric tumorigenesis via M3 receptor–mediated Wnt signaling in the stem cells, and that denervation might represent a feasible strategy for the control of gastric cancer.

No abstract

Summary Within the gastrointestinal stem cell niche, nerves help to regulate both normal and neoplastic stem cell dynamics. Here, we reveal the mechanisms underlying the cancer-nerve partnership. We find that Dclk1+ tuft cells and nerves are the main sources of acetylcholine (ACh) within the gastric mucosa. Cholinergic stimulation of the gastric epithelium induced nerve growth factor (NGF) expression, and in turn NGF overexpression within gastric epithelium expanded enteric nerves and promoted carcinogenesis. Ablation of Dclk1+ cells or blockade of NGF/Trk signaling inhibited epithelial proliferation and tumorigenesis in a muscarinic acetylcholine receptor-3 (M3R)-dependent manner, in part through suppression of Yes-Associated Protein (YAP) function. This feed-forward ACh-NGF axis activates the gastric cancer niche and offers a compelling target for tumor treatment and prevention.

The Poisson-Boltzmann equation (PBE) is an established model for the electrostatic analysis of biomolecules.The development of advanced computational techniques for the solution of the PBE has been an important topic in the past two decades. This article presents a matched interface and boundary (MIB)-based PBE software package, the MIBPB solver, for electrostatic analysis. The MIBPB has a unique feature that it is the first interface technique-based PBE solver that rigorously enforces the solution and flux continuity conditions at the dielectric interface between the biomolecule and the solvent. For protein molecular surfaces, which may possess troublesome geometrical singularities, the MIB scheme makes the MIBPB by far the only existing PBE solver that is able to deliver the second-order convergence, that is, the accuracy increases four times when the mesh size is halved. The MIBPB method is also equipped with a Dirichlet-toNeumann mapping technique that builds a Green's function approach to analytically resolve the singular charge distribution in biomolecules in order to obtain reliable solutions at meshes as coarse as 1 Å -whereas it usually takes other traditional PB solvers 0.25 Å to reach similar level of reliability. This work further accelerates the rate of convergence of linear equation systems resulting from the MIBPB by using the Krylov subspace (KS) techniques. Condition numbers of the MIBPB matrices are significantly reduced by using appropriate KS solver and preconditioner combinations. Both linear and nonlinear PBE solvers in the MIBPB package are tested by protein-solvent solvation energy calculations and analysis of salt effects on protein-protein binding energies, respectively.

The synthetic channel [acetyl-(LeuSerSerLeuLeuSerLeu)3-CONH2]6 (pore diameter approximately 8 A, length approximately 30 A) is a bundle of six alpha-helices with blocked termini. This simple channel has complex properties, which are difficult to explain, even qualitatively, by traditional theories: its single-channel currents rectify in symmetrical solutions and its selectivity (defined by reversal potential) is a sensitive function of bathing solution. These complex properties can be fit quantitatively if the channel has fixed charge at its ends, forming a kind of macrodipole, bracketing a central charged region, and the shielding of the fixed charges is described by the Poisson-Nernst-Planck (PNP) equations. PNP fits current voltage relations measured in 15 solutions with an r.m.s. error of 3.6% using four adjustable parameters: the diffusion coefficients in the channel's pore DK = 2.1 x 10(-6) and DCl = 2.6 x 10(-7) cm2/s; and the fixed charge at the ends of the channel of +/- 0.12e (with unequal densities 0.71 M = 0.021e/A on the N-side and -1.9 M = -0.058e/A on the C-side). The fixed charge in the central region is 0.31e (with density P2 = 0.47 M = 0.014e/A). In contrast to traditional theories, PNP computes the electric field in the open channel from all of the charges in the system, by a rapid and accurate numerical procedure. In essence, PNP is a theory of the shielding of fixed (i.e., permanent) charge of the channel by mobile charge and by the ionic atmosphere in and near the channel's pore. The theory fits a wide range of data because the ionic contents and potential profile in the channel change significantly with experimental conditions, as they must, if the channel simultaneously satisfies the Poisson and Nernst-Planck equations and boundary conditions. Qualitatively speaking, the theory shows that small changes in the ionic atmosphere of the channel (i.e., shielding) make big changes in the potential profile and even bigger changes in flux, because potential is a sensitive function of charge and shielding, and flux is an exponential function of potential.