Multiscale Simulations of a Two-Pore Potassium Channel

Multiple-Basin Go̅-Martini for Investigating Conformational Transitions and Environmental Interactions of Proteins

Proteins are inherently dynamic molecules, and their conformational transitions among various states are essential for numerous biological processes, which are often modulated by their interactions with surrounding environments. Although molecular dynamics (MD) simulations are widely used to investigate these transitions, all-atom (AA) methods are often limited by short time scales and high computational costs, and coarse-grained (CG) implicit-solvent Go̅-like models are usually incapable of studying the interactions between proteins and their environments. Here, we present an approach called Multiple-basin Go̅-Martini, which combines the recent Go̅-Martini model with an exponential mixing scheme to facilitate the simulation of spontaneous protein conformational transitions in explicit environments. We demonstrate the versatility of our method through five diverse case studies: GlnBP, Arc, Hinge, SemiSWEET, and TRAAK, representing ligand-binding proteins, fold-switching proteins, de novo designed proteins, transporters, and mechanosensitive ion channels, respectively. Multiple-basin Go̅-Martini offers a new computational tool for investigating protein conformational transitions, identifying key intermediate states, and elucidating essential interactions between proteins and their environments, particularly protein-membrane interactions. In addition, this approach can efficiently generate thermodynamically meaningful data sets of protein conformational space, which may enhance deep learning-based models for predicting protein conformation distributions.

Ion Conduction Mechanism as a Fingerprint of Potassium Channels

K⁺-channels are membrane proteins that regulate the selective conduction of potassium ions across cell membranes. Although the atomic mechanisms of K⁺ permeation have been extensively investigated, previous work focused on characterizing the selectivity and occupancy of the binding sites, the role of water molecules in the conduction process, or the identification of the minimum energy pathways enabling permeation. Here, we exploit molecular dynamics simulations and the analytical power of Markov state models to perform a comparative study of ion conduction in three distinct channel models. Significant differences emerged in terms of permeation mechanisms and binding site occupancy by potassium ions and/or water molecules from 100 μs cumulative trajectories. We found that, at odds with the current paradigm, each system displays a characteristic permeation mechanism, and thus, there is not a unique way by which potassium ions move through K⁺-channels. The high functional diversity of K⁺-channels can be attributed in part to the differences in conduction features that have emerged from this work. This study provides crucial information and further inspiration for wet-lab chemists designing new synthetic strategies to produce versatile artificial ion channels that emulate membrane transport for their applications in diagnosis, sensors, the next generation of water treatment technologies, etc., as the ability of synthetic channels to transport molecular ions across a bilayer in a controlled way is usually governed through the choice of metal ions, their oxidation states, or their coordination geometries.

Quantum Calculations on Ion Channels: Why Are They More Useful Than Classical Calculations, and for Which Processes Are They Essential?

There are reasons to consider quantum calculations to be necessary for ion channels, for two types of reasons. The calculations must account for charge transfer, and the possible switching of hydrogen bonds, which are very difficult with classical force fields. Without understanding charge transfer and hydrogen bonding in detail, the channel cannot be understood. Thus, although classical approximations to the correct force fields are possible, they are unable to reproduce at least some details of the behavior of a system that has atomic scale. However, there is a second class of effects that is essentially quantum mechanical. There are two types of such phenomena: exchange and correlation energies, which have no classical analogues, and tunneling. Tunneling, an intrinsically quantum phenomenon, may well play a critical role in initiating a proton cascade critical to gating. As there is no classical analogue of tunneling, this cannot be approximated classically. Finally, there are energy terms, exchange and correlation energy, whose values can be approximated classically, but these approximations must be subsumed within classical terms, and as a result, will not have the correct dependence on interatomic distances. Charge transfer, and tunneling, require quantum calculations for ion channels. Some results of quantum calculations are shown.

Biophysics and Modeling of Mechanotransduction in Neurons: A Review

Mechanosensing is a key feature through which organisms can receive inputs from the environment and convert them into specific functional and behavioral outputs. Mechanosensation occurs in many cells and tissues, regulating a plethora of molecular processes based on the distribution of forces and stresses both at the cell membrane and at the intracellular organelles levels, through complex interactions between cells’ microstructures, cytoskeleton, and extracellular matrix. Although several primary and secondary mechanisms have been shown to contribute to mechanosensation, a fundamental pathway in simple organisms and mammals involves the presence of specialized sensory neurons and the presence of different types of mechanosensitive ion channels on the neuronal cell membrane. In this contribution, we present a review of the main ion channels which have been proven to be significantly involved in mechanotransduction in neurons. Further, we discuss recent studies focused on the biological mechanisms and modeling of mechanosensitive ion channels’ gating, and on mechanotransduction modeling at different scales and levels of details.

Conduction and Gating Properties of the TRAAK Channel from Molecular Dynamics Simulations with Different Force Fields

In recent years, the K2P family of potassium channels has been the subject of intense research activity. Owing to the complex function and regulation of this family of ion channels, it is common practice to complement experimental findings with the atomistic description provided by computational approaches such as molecular dynamics (MD) simulations, especially, in light of the unprecedented timescales accessible at present. However, despite recent substantial improvements, the accuracy of MD simulations is still undermined by the intrinsic limitations of force fields. Here, we systematically assessed the performance of the most popular force fields employed to study ion channels at timescales that are orders of magnitude greater than the ones accessible when these energy functions were first developed. Using 32 μs of trajectories, we investigated the dynamics of a member of the K2P ion channel family, the TRAAK channel, using two established force fields in simulations of biological systems: AMBER and CHARMM. We found that while results are comparable on the nanosecond timescales, significant inconsistencies arise at microsecond timescales.

Interfacial Binding Sites for Cholesterol on Kir, Kv, K2P, and Related Potassium Channels

Inwardly rectifying, voltage-gated, two-pore domain, and related K+ channels are located in eukaryotic membranes rich in cholesterol. Here, molecular docking is used to detect specific binding sites ("hot spots") for cholesterol on K+ channels with characteristics that match those of known cholesterol binding sites. The transmembrane surfaces of all available high-resolution structures for K+ channels were swept for potential binding sites. Cholesterol poses were found to be located largely in hollows between protein ridges. A comparison between cholesterol poses and resolved phospholipids suggests that not all cholesterol molecules binding to the transmembrane surface of a K+ channel will result in displacement of a phospholipid molecule from the surface. Competition between cholesterol binding and binding of anionic phospholipids essential for activity could explain some of the effects of cholesterol on channel function.

Interaction between mechanosensitive channels embedded in lipid membrane

The study of the gating mechanism of mechanosensitive channels opens a window to the exploration of how different mechanical stimuli induce adaptive cellular behaviors of both the protein and the lipid, across different time and length scales. In this work, through a molecular dynamics-decorated finite element method (MDeFEM), the gating behavior of mechanosensitive channels of small conductance (MscS) in Escherichia coli (E. coli) is studied upon membrane stretch or global bending. The local membrane curvature around MscS is incorporated, as well as multiple MscL (mechanosensitive channels of large conductance) molecules in proximity to MscS. The local membrane curvature is found to delay MscS opening and diminishes moderately upon membrane stretching. Mimicking the insertion of lysophosphatidylcholine (LPC) molecules into the lipid, both downward and upward bending can active MscS, as long as the global membrane curvature radius reaches 34 nm. Based on the different MscS pore evolutions observed with the presence of one or more MscLs nearby, we propose that when coreconstituted, multiple MscL molecules tend to be located at the local membrane curvature zone around MscS. In another word, as MscL "swims around" in the lipid bilayer, it can be trapped by the membrane's local curvature. Collectively, the current study provides valuable insights into the interplay between mechanosensitive channels and lipid membrane at structural and physical levels, and specific predictions are proposed for further experimental investigations.

Brownian Dynamics Simulations of Biological Molecules

Brownian dynamics (BD) is a technique for carrying out computer simulations of physical systems that are driven by thermal fluctuations. Biological systems at the macromolecular and cellular level, while falling in the gap between well-established atomic-level models and continuum models, are especially suitable for such simulations. We present a brief history, examples of important biological processes that are driven by thermal motion, and those that have been profitably studied by BD. We also present some of the challenges facing developers of algorithms and software, especially in the attempt to simulate larger systems more accurately and for longer times.

Kinetics of Drug Binding and Residence Time

The kinetics of drug binding and unbinding is assuming an increasingly crucial role in the long, costly process of bringing a new medicine to patients. For example, the time a drug spends in contact with its biological target is known as residence time (the inverse of the kinetic constant of the drug-target unbinding, 1/ koff). Recent reports suggest that residence time could predict drug efficacy in vivo, perhaps even more effectively than conventional thermodynamic parameters (free energy, enthalpy, entropy). There are many experimental and computational methods for predicting drug-target residence time at an early stage of drug discovery programs. Here, we review and discuss the methodological approaches to estimating drug binding kinetics and residence time. We first introduce the theoretical background of drug binding kinetics from a physicochemical standpoint. We then analyze the recent literature in the field, starting from the experimental methodologies and applications thereof and moving to theoretical and computational approaches to the kinetics of drug binding and unbinding. We acknowledge the central role of molecular dynamics and related methods, which comprise a great number of the computational methods and applications reviewed here. However, we also consider kinetic Monte Carlo. We conclude with the outlook that drug (un)binding kinetics may soon become a go/no go step in the discovery and development of new medicines.

Gating and inactivation of mechanosensitive channels of small conductance: A continuum mechanics study

Mechanosensitive channels of small conductance (MscS) in Escherichia coli (E. coli) serve as a paradigm for understanding the gating behaviors of the MscS family of ion channels. In this work, we develop a continuum mechanics framework to explore the conformational states of MscS during the gating transition. A complete gating transition trajectory from the closed to the open state along with partially open intermediates is obtained, and the open structure is close to the available structural model from crystallographic studies. The computational efficiency of the modeling framework makes it possible to explore the roles of various structural elements (e.g., loops that connect transmembrane helices) and specific interactions in the gating transition. It is observed that removing either the Asp62-Arg131 salt bridge or the Phe68-Leu111 non-polar interaction leads to essentially non-conducting structures even with a membrane tension close to the lysis limit. The loop connecting TM2 (the second transmembrane helix) and TM3 is found to be essential for force transmission during gating, while the loop connecting TM1 and TM2 does not make any major contribution. Based on the different structural evolutions observed when the TM3 kink is treated as a loop or a helical segment, we propose that the helical propensity of the kink plays a central role in inactivation; i.e., under prolonged sub-threshold membrane tension, transition of the initially flexible loop to a helical segment in TM3 may lead to MscS inactivation. Finally, the gating transition of MscS under different transmembrane voltages is explored and found to be essentially voltage independent. Collectively, results from the current continuum mechanics analysis provide further insights into the gating transition of MscS at structural and physical levels, and specific predictions are proposed for further experimental investigations.

Mutations in KCNK4 that Affect Gating Cause a Recognizable Neurodevelopmental Syndrome

Aberrant activation or inhibition of potassium (K⁺) currents across the plasma membrane of cells has been causally linked to altered neurotransmission, cardiac arrhythmias, endocrine dysfunction, and (more rarely) perturbed developmental processes. The K⁺ channel subfamily K member 4 (KCNK4), also known as TRAAK (TWIK-related arachidonic acid-stimulated K⁺ channel), belongs to the mechano-gated ion channels of the TRAAK/TREK subfamily of two-pore-domain (K2P) K⁺ channels. While K2P channels are well known to contribute to the resting membrane potential and cellular excitability, their involvement in pathophysiological processes remains largely uncharacterized. We report that de novo missense mutations in KCNK4 cause a recognizable syndrome with a distinctive facial gestalt, for which we propose the acronym FHEIG (facial dysmorphism, hypertrichosis, epilepsy, intellectual disability/developmental delay, and gingival overgrowth). Patch-clamp analyses documented a significant gain of function of the identified KCNK4 channel mutants basally and impaired sensitivity to mechanical stimulation and arachidonic acid. Co-expression experiments indicated a dominant behavior of the disease-causing mutations. Molecular dynamics simulations consistently indicated that mutations favor sealing of the lateral intramembrane fenestration that has been proposed to negatively control K⁺ flow by allowing lipid access to the central cavity of the channel. Overall, our findings illustrate the pleiotropic effect of dysregulated KCNK4 function and provide support to the hypothesis of a gating mechanism based on the lateral fenestrations of K2P channels.

Microscopic view of lipids and their diverse biological functions

Biological membranes and their diverse lipid constituents play key roles in a broad spectrum of cellular and physiological processes. Characterization of membrane-associated phenomena at a microscopic level is therefore essential to our fundamental understanding of such processes. Due to the semi-fluid and dynamic nature of lipid bilayers, and their complex compositions, detailed characterization of biological membranes at an atomic scale has been refractory to experimental approaches. Computational modeling and simulation offer a highly complementary toolset with sufficient spatial and temporal resolutions to fill this gap. Here, we review recent molecular dynamics studies focusing on the diversity of lipid composition of biological membranes, or aiming at the characterization of lipid-protein interaction, with the overall goal of dissecting how lipids impact biological roles of the cellular membranes.

Energetic differences between non-domain-swapped and domain-swapped chain connectivities in the K2P potassium channel TRAAK

Two-pore domain (K2P) channels are twofold symmetric K⁺ channels which control cell excitability by enabling the leak of potassium ions from cells in response to physicochemical stimuli. Crystallization of K2P channels revealed the presence of several structural features, which include an external cap. In the available crystallographic structures, the cap is present as non-domain-swapped (NDS) and domain-swapped (DS) chain conformations, where DS chain conformation exchanges two opposing outer helices 180° around the channel. In this work, energy differences between the residues located at the highest point of the cap in NDS and DS conformations were evaluated for TRAAK, a K2P channel that was crystallized in both conformations. Results indicated a preference for DS conformation, but this result is not extensible to TASK K2P channels.

In the available crystallographic structures of K2P channels, the cap is present as non-domain-swapped (NDS) and domain-swapped (DS) chain conformations.