Publications

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.

Cell-surface polysaccharides are essential to many aspects of physiology, serving as a highly conserved evolutionary feature of life and as an important part of the innate immune system in mammals. Here, as simplified biophysical models of these sugar coatings, we present results of molecular dynamics simulations of hyaluronic acid and heparin brushes that show important effects of ion pairing, water dielectric decrease, and coion exclusion. As in prior studies of macromolecular crowding under physiologically relevant salt concentrations, our results show equilibria with electroneutrality attained through screening and pairing of brush anionic charges by monovalent cations at the atomistic detail. Most surprising is the reversal of the Donnan potential obtained from both nonpolarizable and Drude polarizable force fields, in contrast to what would be expected based on electrostatic Boltzmann partitioning alone. Water dielectric decrement within the brush domain is also associated with Born hydration-driven cation exclusion from the brush. We observe that the primary partition energy attracting cations to attain brush electroneutrality is the ion pairing or salt-bridge energy. Potassium and sodium pairings to glycosaminoglycan carboxylates and sulfates show similar abundance of contact-pairing and solvent-separated pairing. We conclude that in these crowded macromolecular brushes, ion-pairing, Born-hydration, and electrostatic potential energies all contribute to attain electroneutrality and should therefore contribute in mean-field models to accurately represent brush electrostatics.

Angiotensin II type 1 receptor (AT1R) blockers (ARBs) are among the most prescribed drugs. However, ARB effectiveness varies widely, which may be due to non-synonymous single nucleotide polymorphisms (nsSNPs) within the AT1R gene. The AT1R coding sequence contains over 100 nsSNPs; therefore, this study embarked on determining which nsSNPs may abrogate the binding of selective ARBs. The crystal structure of olmesartan-bound human AT1R (PDB:4ZUD) served as a template to create an inactive apo-AT1R via molecular dynamics simulation (n = 3). All simulations resulted in a water accessible ligand-binding pocket that lacked sodium ions. The model remained inactive displaying little movement in the receptor core; however, helix 8 showed considerable flexibility. A single frame representing the average stable AT1R was used as a template to dock Olmesartan via AutoDock 4.2, MOE, and AutoDock Vina to obtain predicted binding poses and mean Boltzmann weighted average affinity. The docking results did not match the known pose and affinity of Olmesartan. Thus, an optimization protocol was initiated using AutoDock 4.2 that provided more accurate poses and affinity for Olmesartan (n = 6). Atomic models of 103 of the known human AT1R polymorphisms were constructed using the molecular dynamics equilibrated apo-AT1R. Each of the eight ARBs was then docked, using ARB-optimized parameters, to each polymorphic AT1R (n = 6). Although each nsSNP has a negligible effect on the global AT1R structure, most nsSNPs drastically alter a sub-set of ARBs affinity to the AT1R. Alterations within N298 -L314 strongly effected predicted ARB affinity, which aligns with early mutagenesis studies. The current study demonstrates the potential of utilizing in silico approaches towards personalized ARB therapy. The results presented here will guide further biochemical studies and refinement of the model to increase the accuracy of the prediction of ARB resistance in order to increase overall ARB effectiveness.

Protein dynamic network analysis provides a powerful tool for investigating protein allosteric regulation. We recently developed a current-flow betweenness scheme for protein network analysis and demonstrated that this method, that is, using current-flow betweenness as edge weights, is more appropriate and more robust for investigating the signal transmission between two predefined protein residues or domains as compared with direct usage of correlation scores as edge weights. Here we seek to expand the current-flow scheme to study allosteric regulations involving protein-protein binding. Specifically, we investigated three gain-of-function mutations located at the binding interface of ALK2 (also known as ACVR1) kinase and its inhibitory protein FKBP12. We first searched for the optimal smoothing function for contact network construction and then calculated the subnetwork between FKBP12 protein and the kinase ATP binding site using current-flow betweenness. By comparing the networks between the wild-type and three mutants, we have identified statistically significant changes in the protein-protein networks that are common among all three mutants that allosterically shift the kinase toward a catalytically competent configuration. (c) 2019 Wiley Periodicals, Inc.