Publications

Our laboratory has shown that calpain-2 activation in the brain following acute injury is directly related to neuronal damage and the long-term functional consequences of the injury, while calpain-1 activation is generally neuroprotective and calpain-1 deletion exacerbates neuronal injury. We have also shown that a relatively selective calpain-2 inhibitor, referred to as C2I, enhanced long-term potentiation and learning and memory, and provided neuroprotection in the controlled cortical impact (CCI) model of traumatic brain injury (TBI) in mice. Using molecular dynamic simulation and Site Identification by Ligand Competitive Saturation (SILCS) software, we generated about 130 analogs of C2I and tested them in a number of in vitro and in vivo assays. These led to the identification of two interesting compounds, NA-112 and NA-184. Further analyses indicated that NA-184, (S)-2-(3-benzylureido)-N-((R,S)-1-((3-chloro-2-methoxybenzyl)amino)-1,2-dioxopentan-3-yl)-4-methylpentanamide, selectively and dose-dependent inhibited calpain-2 activity without evident inhibition of calpain-1 at the tested concentrations in mouse brain tissues and human cell lines. Like NA-112, NA-184 inhibited TBI-induced calpain-2 activation and cell death in mice and rats, both male and females. Pharmacokinetic and pharmacodynamic analyses indicated that NA-184 exhibited properties, including stability in plasma and liver and blood-brain barrier permeability, that make it a good clinical candidate for the treatment of TBI.

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.

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.