Publications

Current anticancer therapies cannot eliminate all cancer cells, which hijack normal arginine methylation as a means to promote their maintenance via unknown mechanisms. Here we show that targeting protein arginine N-methyltransferase 9 (PRMT9), whose activities are elevated in blasts and leukemia stem cells (LSCs) from patients with acute myeloid leukemia (AML), eliminates disease via cancer-intrinsic mechanisms and cancer-extrinsic type I interferon (IFN)-associated immunity. PRMT9 ablation in AML cells decreased the arginine methylation of regulators of RNA translation and the DNA damage response, suppressing cell survival. Notably, PRMT9 inhibition promoted DNA damage and activated cyclic GMP-AMP synthase, which underlies the type I IFN response. Genetically activating cyclic GMP-AMP synthase in AML cells blocked leukemogenesis. We also report synergy of a PRMT9 inhibitor with anti-programmed cell death protein 1 in eradicating AML. Overall, we conclude that PRMT9 functions in survival and immune evasion of both LSCs and non-LSCs; targeting PRMT9 may represent a potential anticancer strategy.

Connexin hemichannels were identified as the first members of the eukaryotic large-pore channel family that mediate permeation of both atomic ions and small molecules between the intracellular and extracellular environments. The conventional view is that their pore is a large passive conduit through which both ions and molecules diffuse in a similar manner. In stark contrast to this notion, we demonstrate that the permeation of ions and of molecules in connexin hemichannels can be uncoupled and differentially regulated. We find that human connexin mutations that produce pathologies and were previously thought to be loss-of-function mutations due to the lack of ionic currents are still capable of mediating the passive transport of molecules with kinetics close to those of wild-type channels. This molecular transport displays saturability in the micromolar range, selectivity, and competitive inhibition, properties that are tuned by specific interactions between the permeating molecules and the N-terminal domain that lies within the pore-a general feature of large-pore channels. We propose that connexin hemichannels and, likely, other large-pore channels, are hybrid channel/transporter-like proteins that might switch between these two modes to promote selective ion conduction or autocrine/paracrine molecular signaling in health and disease processes.

Many large protein machines function through an interplay between large-scale movements and intricate conformational changes. Understanding functional motions of these proteins through simulations becomes challenging for both all-atom and coarse-grained (CG) modeling techniques because neither approach alone can readily capture the full details of these motions. In this study, we develop a multiscale model by employing the popular MARTINI CG model to represent a heterogeneous environment and structurally stable proteins and using the united-atom (UA) model PACE to describe proteins undergoing subtle conformational changes. PACE was previously developed to be compatible with the MARTINI solvent and membrane. Here, we couple the protein descriptions of the two models by directly mixing UA and CG interaction parameters to greatly simplify parameter determination. Through extensive validations with diverse protein systems in solution or membrane, we demonstrate that only additional parameter rescaling is needed to enable the resulting model to recover the stability of native structures of proteins under mixed representation. Moreover, we identify the optimal scaling factors that can be applied to various protein systems, rendering the model potentially transferable. To further demonstrate its applicability for realistic systems, we apply the model to a mechanosensitive ion channel Piezo1 that has peripheral arms for sensing membrane tension and a central pore for ion conductance. The model can reproduce the coupling between Piezo1's large-scale arm movement and subtle pore opening in response to membrane stress while consuming much less computational costs than all-atom models. Therefore, our model shows promise for studying functional motions of large protein machines.

We present a checklist to improve the reliability and reproducibility of molecular dynamics simulations and related methods.

While calpains have long been implicated in neurodegeneration, no calpain inhibitor has been developed for the treatment of neurodegeneration. This is partly due to the lack of understanding of the specific functions of most of the 15 members of the calpain family. Work from our laboratory over the last 5-10 years has revealed that calpain-1 and calpain-2, two of the major calpain isoforms in the brain, play opposite roles in both synaptic plasticity/learning and memory and neuroprotection/neurodegeneration. Thus, calpain-1 activation is required for triggering certain forms of synaptic plasticity and for learning some types of information and is neuroprotective. In contrast, calpain-2 activation limits the extent of synaptic plasticity and of learning and is neurodegenerative. These results have been validated with the use of calpain-1 knock-out mice and mice with a selective calpain-2 deletion in excitatory neurons of the forebrain. Through a medicinal chemistry campaign, we have identified a number of selective calpain-2 inhibitors and shown that these inhibitors do facilitate learning of certain tasks and are neuroprotective in a number of animal models of acute neurodegeneration. One of these inhibitors, NA-184, is currently being developed for the treatment of traumatic brain injury, and clinical trials are being planned.

Mechanosensitive PIEZO channels constitute potential pharmacological targets for multiple clinical conditions, spurring the search for potent chemical PIEZO modulators. Among them is Yoda1, a widely used synthetic small molecule PIEZO1 activator discovered through cell-based high-throughput screening. Yoda1 is thought to bind to PIEZO1's mechanosensory arm domain, sandwiched between two transmembrane regions near the channel pore. However, how the binding of Yoda1 to this region promotes channel activation remains elusive. Here, we first demonstrate that cross-linking PIEZO1 repeats A and B with disulfide bridges reduces the effects of Yoda1 in a redox-dependent manner, suggesting that Yoda1 acts by perturbing the contact between these repeats. Using molecular dynamics-based absolute binding free energy simulations, we next show that Yoda1 preferentially occupies a deeper, amphipathic binding site with higher affinity in PIEZO1 open state. Using Yoda1's binding poses in open and closed states, relative binding free energy simulations were conducted in the membrane environment, recapitulating structure-activity relationships of known Yoda1 analogs. Through virtual screening of an 8 million-compound library using computed fragment maps of the Yoda1 binding site, we subsequently identified two chemical scaffolds with agonist activity toward PIEZO1. This study supports a pharmacological model in which Yoda1 activates PIEZO1 by wedging repeats A and B, providing a structural and thermodynamic framework for the rational design of PIEZO1 modulators. Beyond PIEZO channels, the three orthogonal computational approaches employed here represent a promising path toward drug discovery in highly heterogeneous membrane protein systems.

Various all-atom molecular dynamics (MD) simulation methods have been developed to compute free energies and crossing rates of ions and small molecules through ion channels. However, a systemic comparison across different methods is scarce. Using a carbon nanotube as a model of small conductance ion channel, we computed the single-channel permeability for potassium ion using umbrella sampling, Markovian milestoning, and steady-state flux under applied voltage. We show that a slightly modified inhomogeneous solubility-diffusion equation yields a single-channel permeability consistent with the mean first passage time (MFPT) based method. For milestoning, applying cylindrical and spherical bulk boundary conditions yield consistent MFPT if factoring in the effective bulk concentration. The sensitivity of the MFPT to the output frequency of collective variables is highlighted using the convergence and symmetricity of the inward and outward MFPT profiles. The consistent transport kinetic results from all three methods demonstrated the robustness of MD-based methods in computing ion channel permeation. The advantages and disadvantages of each technique are discussed, focusing on the future applications of milestoning in more complex systems.

The connexin family is a diverse group of highly regulated wide-pore channels permeable to biological signaling molecules. Despite the critical roles of connexins in mediating selective molecular signaling in health and disease, the basis of molecular permeation through these pores remains unclear. Here, we report the thermodynamics and kinetics of binding and transport of a second messenger, adenosine-3',5'-cyclophosphate (cAMP), through a connexin26 hemichannel (Cx26). First, inward and outward fluxes of cAMP molecules solvated in KCl solution were obtained from 4 mus of +/- 200 mV simulations. These fluxes data yielded a single-channel permeability of cAMP and cAMP/K(+) permeability ratio consistent with experimentally measured values. The results from voltage simulations were then compared with the potential of mean force (PMF) and the mean first passage times (MFPTs) of a single cAMP without voltage, obtained from a total of 16.5 mus of Voronoi-tessellated Markovian milestoning simulations. Both the voltage simulations and the milestoning simulations revealed two cAMP-binding sites, for which the binding constants KD and dissociation rates koff were computed from PMF and MFPTs. The protein dipole inside the pore produces an asymmetric PMF, reflected in unequal cAMP MFPTs in each direction once within the pore. The free energy profiles under opposite voltages were derived from the milestoning PMF and revealed the interplay between voltage and channel polarity on the total free energy. In addition, we show how these factors influence the cAMP dipole vector during permeation, and how cAMP affects the local and nonlocal pore diameter in a position-dependent manner.

Mechanosensitive Piezo1 channels are essential mechanotransduction proteins in eukaryotes. Their curved transmembrane domains, called arms, create a convex membrane deformation, or footprint, which is predicted to flatten in response to increased membrane tension. Here, using a hyperbolic tangent model, we show that, due to the intrinsic bending rigidity of the membrane, the overlap of neighboring Piezo1 footprints produces a flattening of the Piezo1 footprints and arms. Multiple all-atom molecular dynamics simulations of Piezo1 further reveal that this tension-independent flattening is accompanied by gating motions that open an activation gate in the pore. This open state recapitulates experimentally obtained ionic selectivity, unitary conductance, and mutant phenotypes. Tracking ion permeation along the open pore reveals the presence of intracellular and extracellular fenestrations acting as cation-selective sites. Simulations also reveal multiple potential binding sites for phosphatidylinositol 4,5-bisphosphate. We propose that the overlap of Piezo channel footprints may act as a cooperative mechanism to regulate channel activity.