Permeation and block of the Kv1.2 channel examined using brownian and molecular dynamics

Targeting thermoTRP ion channels: in silico preclinical approaches and opportunities

INTRODUCTION

A myriad of cellular pathophysiological responses are mediated by polymodal ion channels that respond to chemical and physical stimuli such as thermoTRP channels. Intriguingly, these channels are pivotal therapeutic targets with limited clinical pharmacology. In silico methods offer an unprecedented opportunity for discovering new lead compounds targeting thermoTRP channels with improved pharmacological activity and therapeutic index.

AREAS COVERED

This article reviews the progress on thermoTRP channel pharmacology because of (i) advances in solving their atomic structure using cryo-electron microscopy and, (ii) progress on computational techniques including homology modeling, molecular docking, virtual screening, molecular dynamics, ADME/Tox and artificial intelligence. Together, they have increased the number of lead compounds with clinical potential to treat a variety of pathologies. We used original and review articles from Pubmed (1997-2020), as well as the clinicaltrials.gov database, containing the terms thermoTRP, artificial intelligence, docking, and molecular dynamics.

EXPERT OPINION

The atomic structure of thermoTRP channels along with computational methods constitute a realistic first line strategy for designing drug candidates with improved pharmacology and clinical translation. In silico approaches can also help predict potential side-effects that can limit clinical development of drug candidates. Together, they should provide drug candidates with upgraded therapeutic properties.

Queueing arrival and release mechanism for K+ permeation through a potassium channel

The mechanism underlying ion permeation through potassium channels still remains controversial. K+ ions permeate across a narrow selectivity filter (SF) in a single file. Conventional scenarios assume that K+ ions are tightly bound in the SF, and, thus, they are displaced from their energy well by ion-ion repulsion with an incoming ion. This tight coupling between entering and exiting ions has been called the "knock-on" mechanism. However, this paradigm is contradicted by experimental data measuring the water-ion flux coupling ratio, demonstrating fewer ion occupancies. Here, the results of molecular dynamics simulations of permeation through the KcsA potassium channel revealed an alternative mechanism. In the aligned ions in the SF (an ion queue), the outermost K+ was readily and spontaneously released toward the extracellular space, and the affinity of the relevant ion was ~ 50 mM. Based on this low-affinity regime, a simple queueing mechanism described by loose coupling of entering and exiting ions is proposed.

Conductance simulation of the purinergic P2X2, P2X4, and P2X7 ionic channels using a combined Brownian dynamics and molecular dynamics approach

This paper investigates the application of an original combined approach of molecular and Brownian dynamic methods with quantum chemistry calculations for modeling the process of conductance of ion channels using purinergic P2X family receptors P2X2, P2X4, and P2X7 as a case study. A simplified model of the ionic channel in the lipid bilayer has been developed. A high level of conductance (30 pS) of P2X2 ionic channel together with the key role of Asp349 in forming the selectivity filter of P2X2 has been shown by using this approach. Calculated P2X2 permeability to monovalent cations Li(+), Na(+), and K(+) conforms to the free diffusion coefficient of these ions, which shows the low selectivity of P2X2 ionic channel.

Computational methods of studying the binding of toxins from venomous animals to biological ion channels: theory and applications

The discovery of new drugs that selectively block or modulate ion channels has great potential to provide new treatments for a host of conditions. One promising avenue revolves around modifying or mimicking certain naturally occurring ion channel modulator toxins. This strategy appears to offer the prospect of designing drugs that are both potent and specific. The use of computational modeling is crucial to this endeavor, as it has the potential to provide lower cost alternatives for exploring the effects of new compounds on ion channels. In addition, computational modeling can provide structural information and theoretical understanding that is not easily derivable from experimental results. In this review, we look at the theory and computational methods that are applicable to the study of ion channel modulators. The first section provides an introduction to various theoretical concepts, including force-fields and the statistical mechanics of binding. We then look at various computational techniques available to the researcher, including molecular dynamics, brownian dynamics, and molecular docking systems. The latter section of the review explores applications of these techniques, concentrating on pore blocker and gating modifier toxins of potassium and sodium channels. After first discussing the structural features of these channels, and their modes of block, we provide an in-depth review of past computational work that has been carried out. Finally, we discuss prospects for future developments in the field.

Redox artifacts in electrophysiological recordings

Electrophysiological techniques make use of Ag/AgCl electrodes that are in direct contact with cells or bath. In the bath, electrodes are exposed to numerous experimental conditions and chemical reagents that can modify electrode voltage. We examined voltage offsets created in Ag/AgCl electrodes by exposure to redox reagents used in electrophysiological studies. Voltage offsets were measured in reference to an electrode separated from the solution by an agar bridge. The reducing reagents Tris-2-carboxyethly-phosphine, dithiothreitol (DTT), and glutathione, as well as the oxidizing agent H(2)O(2) used at experimentally relevant concentrations reacted with Ag in the electrodes to produce voltage offsets. Chloride ions and strong acids and bases produced offsets at millimolar concentrations. Electrolytic depletion of the AgCl layer, to replicate voltage clamp and sustained use, resulted in increased sensitivity to flow and DTT. Offsets were sensitive to electrode silver purity and to the amount and method of chloride deposition. For example, exposure to 10 μM DTT produced a voltage offset between 10 and 284 mV depending on the chloride deposition method. Currents generated by these offsets are significant and dependent on membrane conductance and by extension the expression of ion channels and may therefore appear to be biological in origin. These data demonstrate a new source of artifacts in electrophysiological recordings that can affect measurements obtained from a variety of experimental techniques from patch clamp to two-electrode voltage clamp.

Extension of Brownian dynamics for studying blockers of ion channels

We present new Brownian dynamics techniques for studying blockers of ion channels. By treating the channel as a fixed body, simulating the blocker molecules using rigid bodies, and using an implicit water force field with explicit ions, we are able to carry out fast simulations that can be used to investigate the dynamics of block and unblock, deduce binding modes, and calculate binding affinities. We test our program using the NavAb bacterial sodium channel, whose structure was recently solved (Payandeh et al. Nature, 2011, 475, 353-358) in conjunction with the μ-conotoxin PIIIA blocker. We derive an ohmic current-voltage relationship for channel permeation, calculate potentials of mean force for blocker unbinding, and deduce multiple binding modes for the blocker. Our results are shown to be compatible with other computational and experimental results. Finally, we discuss future improvements such as the inclusion of flexible side chains. After these improvements are carried out, we anticipate our program will be an extremely useful new tool that could be used to help develop new drugs to treat a range of ion-channel related diseases.

Conductance properties of the inwardly rectifying channel, Kir3.2: molecular and Brownian dynamics study

Using the recently unveiled crystal structure, and molecular and Brownian dynamics simulations, we elucidate several conductance properties of the inwardly rectifying potassium channel, Kir3.2, which is implicated in cardiac and neurological disorders. We show that the pore is closed by a hydrophobic gating mechanism similar to that observed in Kv1.2. Once open, potassium ions move into, but not out of, the cell. The asymmetrical current-voltage relationship arises from the lack of negatively charged residues at the narrow intracellular mouth of the channel. When four phenylalanine residues guarding the intracellular gate are mutated to glutamate residues, the channel no longer shows inward rectification. Inward rectification is restored in the mutant Kir3.2 when it becomes blocked by intracellular Mg(2+). Tertiapin, a polypeptide toxin isolated from the honey bee, is known to block several subtypes of the inwardly rectifying channels with differing affinities. We identify critical residues in the toxin and Kir3.2 for the formation of the stable complex. A lysine residue of tertiapin protrudes into the selectivity filter of Kir3.2, while two other basic residues of the toxin form hydrogen bonds with acidic residues located just outside the channel entrance. The depth of the potential of mean force encountered by tertiapin is -16.1kT, thus indicating that the channel will be half-blocked by 0.4μM of the toxin.