Optimizing Coarse-Grained Models for Large-Scale Membrane Protein Simulation.

Coarse-grained (CG) models have been developed for studying membrane proteins at physiologically relevant scales. Such methods, including popular CG lipid models, exhibit stability and efficiency at moderate scales, but they can become impractical or even unusable beyond a critical size due to various technical issues. Here, we report that these scale-dependent issues can arise from progressively slower relaxation dynamics and become confounded by unforeseen instabilities observed only at larger scales. To address these issues, we systemically optimized a 4-site solvent-free CG lipid model that is suitable for conducting micron-scale molecular dynamics simulations of membrane proteins under various membrane properties. We applied this lipid model to explore the long-range membrane deformation induced by a large mechanosensitive ion channel, PIEZO. We show that the optimized CG models are powerful in elucidating the structural and dynamic interplay between PIEZO and the membrane. Furthermore, we anticipate that our methodological insights can prove useful for resolving issues stemming from scale-dependent limitations of similar CG methodologies.

Identification and neuroprotective properties of NA-184, a calpain-2 inhibitor.

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.

Simulations of Functional Motions of Super Large Biomolecules with a Mixed-Resolution Model.

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.

Targeting PRMT9-mediated arginine methylation suppresses cancer stem cell maintenance and elicits cGAS-mediated anticancer immunity.

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.

Calpain-2 Inhibitors as Therapy for Traumatic Brain Injury.

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.

Importance of molecular dynamics equilibrium protocol on protein-lipid interaction near channel pore.

Multiscale molecular dynamics simulations using Martini coarse-grained (CG) and all-atom (AA) force fields are commonly used in membrane protein studies. In particular, reverse mapping an equilibrated CG model to an AA model offers an efficient way for preparing large membrane protein systems with complex protein shapes and lipid compositions. Here, we report that this hybrid CG-equilibrium-AA-production protocol may artificially increase lipid density and decrease hydration in ion channel pores walled with transmembrane gaps. To understand the origin of this conundrum, we conducted replicas of CG, AA, and CG reverse-mapped AA simulations of the pore domain of the mechanosensitive Piezo1 channel in a nonconducting conformation. Lipid/water density analysis and free energy calculations reveal that the lack of initial pore hydration allows excessive lipids to enter the upper pore lumen through gaps between pore helices during CG simulation. Due to the mismatch between CG and AA lipid kinetics, these pore lipids remain trapped in the subsequent AA simulations, despite unfavorable binding free energy. We tested several CG equilibrium protocols and found that a protocol restraining the whole lipid produces pore hydration consistent with AA results, thus eliminating this artifact for further studies of lipid gating and protein-lipid interactions.

Unifying Single-Channel Permeability From Rare-Event Sampling and Steady-State Flux.

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.

Mechanism-Based and Computational-Driven Covalent Drug Design.

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.

Free energy and kinetics of cAMP permeation through connexin26 via applied voltage and milestoning.

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 µs 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 µs 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.

Mechanical properties of anionic asymmetric bilayers from atomistic simulations.

Mechanotransduction, the biological response to mechanical stress, is often initiated by activation of mechanosensitive (MS) proteins upon mechanically induced deformations of the cell membrane. A current challenge in fully understanding this process is in predicting how lipid bilayers deform upon the application of mechanical stress. In this context, it is now well established that anionic lipids influence the function of many proteins. Here, we test the hypothesis that anionic lipids could indirectly modulate MS proteins by alteration of the lipid bilayer mechanical properties. Using all-atom molecular dynamics simulations, we computed the bilayer bending rigidity (KC), the area compressibility (KA), and the surface shear viscosity (ηm) of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (PC) lipid bilayers with and without phosphatidylserine (PS) or phosphatidylinositol bisphosphate (PIP2) at physiological concentrations in the lower leaflet. Tensionless leaflets were first checked for each asymmetric bilayer model, and a formula for embedding an asymmetric channel in an asymmetric bilayer is proposed. Results from two different sized bilayers show consistently that the addition of 20% surface charge in the lower leaflet of the PC bilayer with PIP2 has minimal impact on its mechanical properties, while PS reduced the bilayer bending rigidity by 22%. As a comparison, supplementing the PIP2-enriched PC membrane with 30% cholesterol, a known rigidifying steroid lipid, produces a significant increase in all three mechanical constants. Analysis of pairwise splay moduli suggests that the effect of anionic lipids on bilayer bending rigidity largely depends on the number of anionic lipid pairs formed during simulations. The potential implication of bilayer bending rigidity is discussed in the framework of MS piezo channels.

Concepts, Practices, and Interactive Tutorial for Allosteric Network Analysis of Molecular Dynamics Simulations.

Over the past decade, concepts of network theory in combination with dynamical information from conformational ensembles have been widely applied to gain insights in understanding allosteric regulation in biomolecules. In this chapter, we introduce the basic theories and protocols used in protein dynamics network analysis through a series of interactive python Jupyter notebook scripts. While various network analysis methods exist in the literature, here we focus on the two popular methods based on correlated atomic motions and pairwise interaction energies. While the tutorial is based on a small prototypic protein, the workflow and protocol introduced here are optimized to handle large membrane proteins.

Crowding-induced opening of the mechanosensitive Piezo1 channel in silico.

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.

Ion Channels in Biophysics and Physiology: Methods & Challenges to Study Mechanosensitive Ion Channels.

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.