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.
In a seminal work published in 1950, Sir B. Katz showed that the electrical response of the frog muscle spindle varies directly with the rate and amplitude of muscle stretch. This observation led him to propose the existence of a piezoelectric substance in this organ, setting the stage for the field of mechanobiology (Katz, J Physiol 111, 261-282, 1950). Despite this early work, the identity of the molecules responsible for the conversion of mechanical stimuli into biological signals has remained hidden for decades. This delay is often attributed to the inherent difficulty to precisely quantify the mechanical deformations of biological samples. In contrast to other forms of stimuli such as ligand concentration and membrane potential, quantifying mechanical deformations of cell membranes is not trivial. Mechanical forces produce a complex array of membrane deformations including bending, thinning, compression, expansion, and shear, and thus, have components in many strain dimensions. In addition, due to the viscoelastic nature of cells, these deformations may have linear and nonlinear components. In spite of these experimental challenges, Sukharev et al. cloned the first mechanosensitive ion channel from the bacteria E. coli in the mid-1990s (Sukharev et al. Nature, 265-268, 1994). Two decades later, several protein families encompassing dozens of eukaryotic mechanosensitive ion channels have been identified, depicting an astonishing diversity of force-activated molecular machines. In this chapter, we intend to provide an overview of the current state of knowledge and technical challenges to study how cell membranes deform upon mechanical stress and how ion channel proteins detect these deformations to engage homeostatic cellular responses.
Covalent drugs offer higher efficacy and longer duration of action than their noncovalent counterparts. Significant advances in computational methods for modeling covalent drugs are poised to shift the paradigm of small molecule therapeutics within the next decade. This viewpoint discusses the advantages of a two-state model for ranking reversible and irreversible covalent ligands and of more complex models for dissecting reaction mechanisms. The relation between these models highlights the complexity and diversity of covalent drug binding and provides opportunities for mechanism-based rational design.
