Modeling the human Nav1.5 sodium channel: structural and mechanistic insights of ion permeation and drug blockade

Identification of Lauric Acid as a Potent Sodium Channel NaV1.5 Blocker from Compound Chinese Medicine Wenxin Keli

Purpose

The major cardiac voltage-gated sodium channel NaV1.5 (INa) is essential for cardiac action potential initiation and subsequent propagation. Compound Chinese medicine Wenxin Keli (WXKL) has been shown to suppress arrhythmias and heart failure. However, its active components have not been fully elucidated. This study focused on identifying the active inhibitor of INa in WXKL and exploring their mode of action in electrophysiological conduction.

Methods

A chemical fraction library was constructed from an aqueous extract of WXKL and screened using an automated patch-clamping system in cells stably expressing the NaV1.5 gene SCN5A. Candidate fractions with INa-inhibition activity were analyzed by HPLC-ESI-IT-TOF-MS and GC-MS to identify the ingredients. NaV1.5 blocker molecules identified by single-cell electrocardiogram were tested in hiPSC-derived cardiomyocytes. We evaluated the SCN5A inhibitory potential of Wenxin Keli effective monomer employing molecular docking and molecular dynamics simulation approaches.

Results

A primary screen of the WXKL chemical library identified five fractions that significantly inhibited the NaV1.5 channel, with one of them rich in poly-saturated fatty acids. Molecular structural characterization revealed the presence of lauric acid, myristic acid, palmitic acid, and stearic acid in the active subfraction. Electrophysiological characterization demonstrated lauric acid (LA) as the most effective monomer for INa-inhibition with an IC50 at 27.40 ± 12.78 μM. LA shifted the steady-state inactivation of INa to more negative potentials and decreased the amplitude of extracellular field potential in hiPSC-derived cardiomyocytes. We demonstrate for the first time that naturally poly-saturated fatty acid, lauric acid, as a potential novel INa blocker. Molecular docking and molecular dynamics simulation suggested that LA binds to the NaV1.5 protein, with a significant binding affinity forming interactions with functionally essential residues and blocks the inward flow of Na+. Mechanistically, lauric acid acts on the fast inactivation of NaV1.5 alter electrophysiology conduction of hiPSC-derived cardiomyocytes and contribute to the antiarrhythmic effect of WXKL.

Conclusion

Lauric acid is a potent blocker for sodium channel NaV1.5 and alleviates arrhythmia via inhibiting INa.

Role of Hydration and Amino Acid Interactions on the Ion Permeation Mechanism in the hNaV1.5 Channel

This study explores the ion selectivity and conduction mechanisms of the hNaV1.5 sodium channel using classical molecular dynamics simulations under an externally applied electric field. Our findings reveal distinct conduction mechanisms for Na+ and K+, primarily driven by differences in their hydration states when they diffuse close to the channel's selective filter (DEKA) and extracellular ring (EEDD). The Na+ ions undergo partial dehydration in the EEDD region, followed by a rehydration step in the DEKA ring, resulting in longer retention times and a deeper free energy minimum compared to K+. Conversely, the K+ ions exhibit a continuous dehydration process, facilitating a smoother passage through these key regions. These results indicate that ion selectivity and conductance are primarily governed by solvation dynamics, which, in turn, depend on the interactions with key charged residues within the channel. Additionally, we show that the delicate energetic balance between the interactions of the ions with the protein residues and with their solvation shells during the dehydration and rehydration processes is not properly captured by the force field. As a consequence, the selectivity of the channel is not well described, indicating that more accurate computational models must be applied to simulate ion conduction through eukaryotic NaV channels.

Structural Insight on the Selectivity of Calyx[4]Arene-Based Inhibitors of Mg2+-Dependent Atp-Hydrolases

Located in plasma membranes, ATP hydrolases are involved in several dynamic transport processes, helping to control the movement of ions across cell membranes. ATP hydrolase acts as a transport protein, converting energy from ATP hydrolysis into transport molecules against their concentration gradients. In addition to energy metabolism and active transport, ATP hydrolase is essential for maintaining cellular homeostasis and cell function. This study focused on the domain architecture model of P-type ATPases, which participate in the reaction cycles of ATP hydrolysis carried out by membrane transport systems - Na+, K+-ATPase and Ca2+, Mg2+-ATPase. Targeted modulation of Na+, K+-ATPase and Ca2+, Mg2+-ATPase by unnatural drugs is of greatest interest due to the lack of known effectors. This new discovery presents a convenient model based on our recent experimental studies of the membrane structures and myocytes of the uterine smooth muscle, the myometrium. This current study strongly supports the fact that nanosized calix[4]arenes functionalised on the upper rings of the macrocycle with biologically active phosphonic acid fragments can serve as selective and potent inhibitors of cation-transporting electroenzymes. This is how we discovered that calix[4]arene of methylenebisphosphonic acid C-97 and calix[4]arene of bis-aminophosphonic acid C-107 selectively and effectively (I0.5 <100 nM) inhibit the activity of Mg2+, ATP-dependent electrogenic Na+ K+ plasma membrane pump. As drug discovery in the field of Mg2+-ATPase inhibitors is uncharted territory, basic research holds the key to explaining and predicting the mechanism of interaction and action of different classes of compounds. In light of the presented results, new calix[4]arene compounds can be used as potent inhibitors of Mg2+, ATP-dependent electrogenic ion pumps.

The sodium/glucose cotransporter 2 inhibitor Empagliflozin inhibits long QT 3 late sodium currents in a mutation specific manner

BACKGROUND

Sodium/glucose cotransporter 2 inhibitors (SGLT2is) like empagliflozin have demonstrated cardioprotective effects in patients with or without diabetes. SGLT2is have been shown to selectively inhibit the late component of cardiac sodium current (late INa). Induction of late INa is the primary mechanism in the pathophysiology of congenital long QT syndrome type 3 (LQT3) gain-of-function mutations in the SCN5A gene encoding Nav1.5. We investigated empagliflozin's effect on late INa in thirteen known LQT3 mutations located in distinct regions of the channel.

METHODS

The whole-cell patch-clamp technique was used to investigate the effect of empagliflozin on late INa in recombinantly expressed Nav1.5 channels containing different LQT3 mutations. Molecular modeling of human Nav1.5 and simulations in a mathematical model of human ventricular myocytes were used to extrapolate our experimental results to excitation-contraction coupling.

RESULTS

Empagliflozin selectively inhibited late INa in LQT3 mutations in the inactivation gate region of Nav1.5, without affecting peak current or channel kinetics. In contrast, empagliflozin inhibited both peak and late INa in mutations in the S4 voltage-sensing regions, altered channel gating, and slowed recovery from inactivation. Empagliflozin had no effect on late/peak INa or channel kinetics in channels with mutations in the putative empagliflozin binding region. Simulation results predict that empagliflozin may have a desirable therapeutic effect in LQT3 mutations in the inactivation gate region.

CONCLUSIONS

Empagliflozin selectively inhibits late INa, without affecting channel kinetics, in LQT3 mutations in the inactivation gate region. Empagliflozin may thus be a promising precision medicine approach for patients with specific LQT3 mutations.

An overview of drug-induced sodium channel blockade and changes in cardiac conduction: Implications for drug safety

The human voltage-gated sodium channel Nav1.5 (hNav1.5/SCN5A) plays a critical role in the initiation and propagation of action potentials in cardiac myocytes, and its modulation by various drugs has significant implications for cardiac safety. Drug-dependent block of Nav1.5 current (INa) can lead to significant alterations in cardiac electrophysiology, potentially resulting in conduction slowing and an increased risk of proarrhythmic events. This review aims to provide a comprehensive overview of the mechanisms by which various pharmacological agents interact with Nav1.5, focusing on the molecular determinants of drug binding and the resultant electrophysiological effects. We discuss the structural features of Nav1.5 that influence drug affinity and specificity. Special attention is given to the concept of state-dependent block, where drug binding is influenced by the conformational state of the channel, and its relevance to therapeutic efficacy and safety. The review also examines the clinical implications of INa block, highlighting case studies of drugs that have been associated with adverse cardiac events, and how the Vaughan-Williams Classification system has been employed to qualify "unsafe" sodium channel block. Furthermore, we explore the methodologies currently used to assess INa block in nonclinical and clinical settings, with the hope of providing a weight of evidence approach including in silico modeling, in vitro electrophysiological assays and in vivo cardiac safety studies for mitigating proarrhythmic risk early in drug discovery. This review underscores the importance of understanding Nav1.5 pharmacology in the context of drug development and cardiac risk assessment.

Structural insights into SARS-CoV-2 main protease conformational plasticity

The spread of different severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants underscores the need for insights into the structural properties of its structural and non-structural proteins. The highly conserved homo-dimeric chymotrypsin-like protease (3CL M^PRO ), belonging to the class of cysteine hydrolases, plays an indispensable role in processing viral polyproteins that are involved in viral replication and transcription. Studies have successfully demonstrated the role of M^PRO as an attractive drug target for designing antiviral treatments because of its importance in the viral life cycle. Herein, we report the structural dynamics of six experimentally solved structures of M^PRO (i.e., 6LU7, 6M03, 6WQF, 6Y2E, 6Y84, and 7BUY including both ligand-free and ligand-bound states) at different resolutions. We have employed a structure-based balanced forcefield, CHARMM36m through state-of-the-art all-atoms molecular dynamics simulations at µ-seconds scale at room temperature (303K) and pH 7.0 to explore their structure-function relationship. The helical domain-III responsible for dimerization mostly contributes to the altered conformational states and destabilization of M^PRO . A keen observation of the high degree of flexibility in the P5 binding pocket adjoining domain II-III highlights the reason for observation of conformational heterogeneity among the structural ensembles of M^PRO . We also observe a differential dynamics of the catalytic pocket residues His41, Cys145, and Asp187, which may lead to catalytic impairment of the monomeric proteases. Among the highly populated conformational states of the six systems, 6LU7 and 7M03 forms the most stable and compact M^PRO conformation with intact catalytic site and structural integrity. Altogether, our findings from this extensive study provides a benchmark to identify physiologically relevant structures of such promising drug targets for structure-based drug design and discovery of potent drug-like compounds having clinical potential.

Relationship among surface electric double layer of cardiomyocyte membrane and toxicology of digoxin and opening of ion channels

We applied a new idea that the potential effect can change the ion adsorption structure on the cell surface to explore the mechanism of digoxin poisoning and the regulation of ion channels. The effects of digoxin on the electrophoretic mobility and behaviors (non-contraction or contraction or autorhythmicity) of cardiomyocytes were observed by single-cell electrophoresis technique (imitate the opening method of in vivo channel) and the method of decomposing surface potential components on the cells. As well as affect the association with electrical activity. The results suggested that the increase of cardiomyocytes transmembrane potential and the Na⁺-K⁺ exchange on the cell surface of the action potential phase 4 caused by the poisoning dose of digoxin, leading to the oscillation of adsorbed ions on the cell surface and the incomplete channel structure, which were the mechanism of cardiac ectopic beats. The results revealed that the opening of ion channels is regulated by the surface electric double layer of the cell membrane.

In Silico Evaluation of Hexamethylene Amiloride Derivatives as Potential Luminal Inhibitors of SARS-CoV-2 E Protein

The coronavirus E proteins are small membrane proteins found in the virus envelope of alpha and beta coronaviruses that have a high degree of overlap in their biochemical and functional properties despite minor sequence variations. The SARS-CoV-2 E is a 75-amino acid transmembrane protein capable of acting as an ion channel when assembled in a pentameric fashion. Various studies have found that hexamethylene amiloride (HMA) can inhibit the ion channel activity of the E protein in bilayers and also inhibit viral replication in cultured cells. Here, we use the available structural data in conjunction with homology modelling to build a comprehensive model of the E protein to assess potential binding sites and molecular interactions of HMA derivatives. Furthermore, we employed an iterative cycle of molecular modelling, extensive docking simulations, molecular dynamics and leveraging steered molecular dynamics to better understand the pore characteristics and quantify the affinity of the bound ligands. Results from this work highlight the potential of acylguanidines as blockers of the E protein and guide the development of subsequent small molecule inhibitors.

"Virtual patch clamp analysis" for predicting the functional significance of pathogenic variants in sodium channels

OBJECTIVE

Opening of voltage-gated sodium channels is crucial for neuronal depolarization. Proper channel opening and influx of Na⁺ through the ion pore, is dependent upon binding of Na⁺ ion to a specific amino-acid motif (DEKA) within the pore. In this study we used molecular dynamic simulations, an advanced bioinformatic tool, to research the dysfunction caused by pathogenic variants in SCN1a, SCN2a and SCN8a genes.

METHOD

Molecular dynamic simulations were performed in six patients: three patients with Dravet syndrome (p.Gly177Ala,p.Ser259Arg and p.Met1267Ile, SCN1a), two patients with early onset drug resistant epilepsy(p.Ala263Val, SCN2a and p.Ile251Arg, SCN8a), and a patient with autism (p.Thr155Ala, SCN2a). After predicting the 3D-structure of mutated proteins by homology modeling, time dependent molecular dynamic simulations were performed, using the Schrödinger algorithm. The opening of the sodium channel, including the detachment of the sodium ion to the DEKA motif and pore diameter were assessed. Results were compared to the existent patch clamp analysis in four patients, and consistency with clinical phenotype was noted.

RESULTS

The Na⁺ ion remained attached to DEKA filter longer when compared to wild type in the p.Gly177Ala, p.Ser259Arg,SCN1a, and p.Thr155Ala, SCN2a variants, consistent with loss-of-function. In contrast, it detached quicker from DEKA than wild type in the p.Ala263Val,SCN2a variant, consistent with gain-of-function. In the p.Met1267Ile,SCN1a variant, detachment from DEKA was quicker, but pore diameter decreased, suggesting partial loss-of-function. In the p.Leu251Arg,SCN8a variant, the pore remained opened longer when compared to wild type, consistent with a gain-of-function. The molecular dynamic simulation results were consistent with the existing patch-clamp analysis studies, as well as the clinical phenotype.

SIGNIFICANCE

Molecular dynamic simulation can be useful in predicting pathogenicity of variants and the disease phenotype, and selecting targeted treatment based on channel dysfunction. Further development of these bioinformatic tools may lead to "virtual patch-clamp analysis".

Simulation and Machine Learning Methods for Ion-Channel Structure Determination, Mechanistic Studies and Drug Design

Ion channels are expressed in almost all living cells, controlling the in-and-out communications, making them ideal drug targets, especially for central nervous system diseases. However, owing to their dynamic nature and the presence of a membrane environment, ion channels remain difficult targets for the past decades. Recent advancement in cryo-electron microscopy and computational methods has shed light on this issue. An explosion in high-resolution ion channel structures paved way for structure-based rational drug design and the state-of-the-art simulation and machine learning techniques dramatically improved the efficiency and effectiveness of computer-aided drug design. Here we present an overview of how simulation and machine learning-based methods fundamentally changed the ion channel-related drug design at different levels, as well as the emerging trends in the field.

Structural modeling of Nav1.5 pore domain in closed state

The voltage-dependent cardiac sodium channel plays a key role in cardiac excitability and conduction and it is the drug target of medically important. However, its atomic- resolution structure is still lack. Here, we report a modeled structure of Na_v1.5 pore domain in closed state. The structure was constructed by Rosetta-membrane homology modeling method based on the template of eukaryotic Na_v channel Na_vPaS and selected by energy and direct coupling analysis (DCA). Moreover, this structure was optimized through molecular dynamical simulation in the lipid membrane bilayer. Finally, to validate the constructed model, the binding energy and binding sites of closed-state local anesthetics (LAs) in the modeled structure were computed by the MM-GBSA method and the results are in agreement with experiments. The modeled structure of Na_v1.5 pore domain in closed state may be useful to explore molecular mechanism of a state-dependent drug binding and helpful for new drug development.

Antiepileptic Drug Tiagabine Does Not Directly Target Key Cardiac Ion Channels Kv11.1, Nav1.5 and Cav1.2

Tiagabine is an antiepileptic drug used for the treatment of partial seizures in humans. Recently, this drug has been found useful in several non-epileptic conditions, including anxiety, chronic pain and sleep disorders. Since tachycardia-an impairment of cardiac rhythm due to cardiac ion channel dysfunction-is one of the most commonly reported non-neurological adverse effects of this drug, in the present paper we have undertaken pharmacological and numerical studies to assess a potential cardiovascular risk associated with the use of tiagabine. A chemical interaction of tiagabine with a model of human voltage-gated ion channels (VGICs) is described using the molecular docking method. The obtained in silico results imply that the adverse effects reported so far in the clinical cardiological of tiagabine could not be directly attributed to its interactions with VGICs. This is also confirmed by the results from the isolated organ studies (i.e., calcium entry blocking properties test) and in vivo (electrocardiogram study) assays of the present research. It was found that tachycardia and other tiagabine-induced cardiac complications are not due to a direct effect of this drug on ventricular depolarization and repolarization.

Cardiac Late Sodium Channel Current Is a Molecular Target for the Sodium/Glucose Cotransporter 2 Inhibitor Empagliflozin

BACKGROUND

SGLT2 (sodium/glucose cotransporter 2) inhibitors exert robust cardioprotective effects against heart failure in patients with diabetes, and there is intense interest to identify the underlying molecular mechanisms that afford this protection. Because the induction of the late component of the cardiac sodium channel current (late-INa) is involved in the etiology of heart failure, we investigated whether these drugs inhibit late-INa.

METHODS

Electrophysiological, in silico molecular docking, molecular, calcium imaging, and whole heart perfusion techniques were used to address this question.

RESULTS

The SGLT2 inhibitor empagliflozin reduced late-INa in cardiomyocytes from mice with heart failure and in cardiac Nav1.5 sodium channels containing the long QT syndrome 3 mutations R1623Q or ΔKPQ. Empagliflozin, dapagliflozin, and canagliflozin are all potent and selective inhibitors of H2O2-induced late-INa (half maximal inhibitory concentration = 0.79, 0.58, and 1.26 µM, respectively) with little effect on peak sodium current. In mouse cardiomyocytes, empagliflozin reduced the incidence of spontaneous calcium transients induced by the late-INa activator veratridine in a similar manner to tetrodotoxin, ranolazine, and lidocaine. The putative binding sites for empagliflozin within Nav1.5 were investigated by simulations of empagliflozin docking to a three-dimensional homology model of human Nav1.5 and point mutagenic approaches. Our results indicate that empagliflozin binds to Nav1.5 in the same region as local anesthetics and ranolazine. In an acute model of myocardial injury, perfusion of isolated mouse hearts with empagliflozin or tetrodotoxin prevented activation of the cardiac NLRP3 (nuclear-binding domain-like receptor 3) inflammasome and improved functional recovery after ischemia.

CONCLUSIONS

Our results provide evidence that late-INa may be an important molecular target in the heart for the SGLT2 inhibitors, contributing to their unexpected cardioprotective effects.

Modelling of an autonomous Nav1.5 channel system as a part of in silico pharmacology study

A homology model of Nav1.5, based mainly on the crystal structures of Nav1.2/1.5 was built, optimized and successfully inserted into the membrane bilayer. We applied steered and free MD simulation protocols for the visualization of the mechanism of Nav1.5 activation. We constrained dihedrals of S4 trigger to introduce a structural tension with further rearrangement and movement of secondary structure elements. From these, we observed an intracellular gate opening and movement of the Lys1419 residue caused by a gradual displacement of the distal S6 α-helix with the extended S4 3-10 helix of voltage-sensing domains (VSD). A construction containing the Lys1419 residue in P-loop also changed its position due to the extension of this helix and subsequent induction of the pore-forming helixes motion. From this point, a double membrane system was generated, implying a free of ligand Nav1.5 protein and on the opposite side its copy containing a docked bupivacaine molecule inside the pore channel. The system can be used for the design of selective inhibitors against the Nav1.7 channel, instead of mixed effect on both channels.

GPCR_LigandClassify.py; a rigorous machine learning classifier for GPCR targeting compounds

The current study describes the construction of various ligand-based machine learning models to be used for drug-repurposing against the family of G-Protein Coupled Receptors (GPCRs). In building these models, we collected > 500,000 data points, encompassing experimentally measured molecular association data of > 160,000 unique ligands against > 250 GPCRs. These data points were retrieved from the GPCR-Ligand Association (GLASS) database. We have used diverse molecular featurization methods to describe the input molecules. Multiple supervised ML algorithms were developed, tested and compared for their accuracy, F scores, as well as for their Matthews' correlation coefficient scores (MCC). Our data suggest that combined with molecular fingerprinting, ensemble decision trees and gradient boosted trees ML algorithms are on the accuracy border of the rather sophisticated deep neural nets (DNNs)-based algorithms. On a test dataset, these models displayed an excellent performance, reaching a ~ 90% classification accuracy. Additionally, we showcase a few examples where our models were able to identify interesting connections between known drugs from the Drug-Bank database and members of the GPCR family of receptors. Our findings are in excellent agreement with previously reported experimental observations in the literature. We hope the models presented in this paper synergize with the currently ongoing interest of applying machine learning modeling in the field of drug repurposing and computational drug discovery in general.

Molecular dynamics and in silico mutagenesis on the reversible inhibitor-bound SARS-CoV-2 main protease complexes reveal the role of lateral pocket in enhancing the ligand affinity

The 2019 novel coronavirus pandemic caused by SARS-CoV-2 remains a serious health threat to humans and there is an urgent need to develop therapeutics against this deadly virus. Recent scientific evidences have suggested that the main protease (Mpro) enzyme in SARS-CoV-2 can be an ideal drug target due to its crucial role in the viral replication and transcription processes. Therefore, there are ongoing research efforts to identify drug candidates against SARS-CoV-2 Mpro that resulted in hundreds of X-ray crystal structures of ligand-bound Mpro complexes in the Protein Data Bank (PDB) describing the interactions of different fragment chemotypes within different sites of the Mpro. In this work, we performed rigorous molecular dynamics (MD) simulation of 62 reversible ligand-Mpro complexes in the PDB to gain mechanistic insights about their interactions at the atomic level. Using a total of over 3 µs long MD trajectories, we characterized different pockets in the apo Mpro structure, and analyzed the dynamic interactions and binding affinity of ligands within those pockets. Our results identified the key residues that stabilize the ligands in the catalytic sites and other pockets of Mpro. Our analyses unraveled the role of a lateral pocket in the catalytic site in Mpro that is critical for enhancing the ligand binding to the enzyme. We also highlighted the important contribution from HIS163 in the lateral pocket towards ligand binding and affinity against Mpro through computational mutation analyses. Further, we revealed the effects of explicit water molecules and Mpro dimerization in the ligand association with the target. Thus, comprehensive molecular-level insights gained from this work can be useful to identify or design potent small molecule inhibitors against SARS-CoV-2 Mpro.

Structure-based virtual screening, molecular dynamics and binding affinity calculations of some potential phytocompounds against SARS-CoV-2

COVID-19 caused by a positive-sense single stranded RNA virus named as severe acute respiratory syndrome-Coronavirus-2 (SARS-CoV-2) triggered the global pandemic. This virus has infected about 10.37 Crores and taken lives of 2.24 Crores people of 213 countries to date. To cope-up this emergency clinical trials are undergoing with some existing drugs like remdesivir, flavipiravir, lopinavir-ritonavir, nafamostat, doxycycline, hydroxy-chloroquine, dexamethasone, etc., despite their severe toxicity and health hazards among diabetics, hypertensive, cardiac patients or normal individuals. The lack of safe and approved treatment for COVID-19 has forced the scientific community to find novel and safe compounds with potential efficacy. This study evaluates a few selective herbal compounds like glucoraphanin, vitexin, niazinin, etc., as a potential inhibitor of the spike protein and 3-chymotrypsin-like protease (3CLpro) or main protease (Mpro) of SARS-COV-2 through in-silico virtual studies such as molecular docking, target analysis, toxicity prediction and ADME prediction and supported by a Molecular-Dynamic simulation. Selective phytocompounds were docked successfully in the binding site of spike glycoprotein and 3CLpro (Mpro) of SARS-CoV-2. In-silico approaches also predict this molecule to have good solubility, pharmacodynamic property and target accuracy through MD simulation and ADME studies. These hit molecules niazinin, vitexin, glucoraphanin also obey Lipinski's rule along with their stable binding towards target protein of the virus, which makes them suitable for further biochemical and cell-based assays followed by clinical investigations to highlight their potential use in COVID-19 treatment.Communicated by Ramaswamy H. Sarma.

Machine learning-based QSAR models to predict sodium ion channel (Nav 1.5) blockers

Aim: Conventional experimental approaches used for the evaluation of the proarrhythmic potential of compounds in the drug discovery process are expensive and time consuming but an integral element in the safety profile required for a new drug to be approved. The voltage-gated sodium ion channel 1.5 (Na_v 1.5), a target known for arrhythmic drugs, causes adverse cardiac complications when the channel is blocked. Results: Machine learning classification and regression models were built to predict the possibility of blocking these channels by small molecules. The finalized models tested with balanced accuracies of 0.88, 0.93 and 0.94 at three thresholds (1, 10 and 30 µmol, respectively). The regression model built to predict the pIC₅₀ of compounds had q² of 0.84 (root-mean-square error = 0.46). Conclusion: The machine learning models that have been built can act as effective filters to screen out the potentially toxic compounds in the early stages of drug discovery.

Atomistic modeling and molecular dynamics analysis of human CaV1.2 channel using external electric field and ion pulling simulations

BACKGROUND

Human Ca_V1.2 (hCav1.2), a calcium selective voltage-gated channel, plays important roles in normal cardiac and neuronal functions. Calcium influx and gating mechanisms leading to the activation of hCa_V1.2 are critical for its functionalities. Lack of an experimentally resolved structure of hCa_V1.2 remains a significant impediment in molecular-level understanding of this channel. This work focuses on building atomistic hCa_V1.2 model and studying calcium influx using computational approaches.

METHODS

We employed homology modeling and molecular dynamics (MD) to build the structure of hCa_V1.2. Subsequently, we employed steered molecular dynamics (SMD) to understand calcium ion permeation in hCa_V1.2.

RESULTS

We report a comprehensive three-dimensional model of a closed state hCa_V1.2 refined under physiological membrane-bound conditions using MD simulations. Our SMD simulations on the model revealed four important barriers for ion permeation: this includes three calcium binding sites formed by the EEEE- and TTTT- rings within the selectivity filter region and a large barrier rendered by the hydrophobic internal gate. Our results also revealed that the first hydration shell of calcium remained intact throughout the simulations, thus playing an important role in ion permeation in hCa_V1.2.

CONCLUSIONS

Our results have provided some important mechanistic insights into the structure, dynamics and ion permeation in hCa_V1.2. The significant barriers for ion permeation formed by the four phenylalanine residues at the internal gate region suggest that this site is important for channel activation.

Vasodilation Elicited by Isoxsuprine, Identified by High-Throughput Virtual Screening of Compound Libraries, Involves Activation of the NO/cGMP and H₂S/KATP Pathways and Blockade of α₁-Adrenoceptors and Calcium Channels

Recently, our research group demonstrated that uvaol and ursolic acid increase NO and H₂S production in aortic tissue. Molecular docking studies showed that both compounds bind with high affinity to endothelial NO synthase (eNOS) and cystathionine gamma-lyase (CSE). The aim of this study was to identify hits with high binding affinity for the triterpene binding-allosteric sites of eNOS and CSE and to evaluate their vasodilator effect. Additionally, the mechanism of action of the most potent compound was explored. A high-throughput virtual screening (HTVS) of 107,373 compounds, obtained from four ZINC database libraries, was performed employing the crystallographic structures of eNOS and CSE. Among the nine top-scoring ligands, isoxsuprine showed the most potent vasodilator effect. Pharmacological evaluation, employing the rat aorta model, indicated that the vasodilation produced by this compound involved activation of the NO/cGMP and H₂S/K_ATP signaling pathways and blockade of α₁-adrenoceptors and L-type voltage-dependent Ca²⁺ channels. Incubation of aorta homogenates in the presence of isoxsuprine caused 2-fold greater levels of H₂S, which supported our preliminary in silico data. This study provides evidence to propose that the vasodilator effect of isoxsuprine involves various mechanisms, which highlights its potential to treat a wide variety of cardiovascular diseases.

Challenges and advances in atomistic simulations of potassium and sodium ion channel gating and permeation

Ion channels are implicated in many essential physiological events such as electrical signal propagation and cellular communication. The advent of K⁺ and Na⁺ ion channel structure determination has facilitated numerous investigations of molecular determinants of their behaviour. At the same time, rapid development of computer hardware and molecular simulation methodologies has made computational studies of large biological molecules in all-atom representation tractable. The concurrent evolution of experimental structural biology with biomolecular computer modelling has yielded mechanistic details of fundamental processes unavailable through experiments alone, such as ion conduction and ion channel gating. This review is a short survey of the atomistic computational investigations of K⁺ and Na⁺ ion channels, focusing on KcsA and several voltage-gated channels from the K_V and Na_V families, which have garnered many successes and engendered several long-standing controversies regarding the nature of their structure-function relationship. We review the latest advancements and challenges facing the field of molecular modelling and simulation regarding the structural and energetic determinants of ion channel function and their agreement with experimental observations.

Attenuation of quorum sensing regulated virulence and biofilm development in Pseudomonas aeruginosa PAO1 by Diaporthe phaseolorum SSP12

In recent years, Pseudomonas aeruginosa PAO1 emerged as the significant pathogenic microorganism in majority of the hospital-acquired infections due to its resistance to the conventional antibiotics by virtue of its highly organized quorum sensing and associated biofilm formation. In the present study, quorum sensing attenuation potential of Diaporthe phaseolorum SSP12 extract was investigated against P. aeruginosa PAO1 amply supported by molecular docking studies. D. phaseolorum SSP12 extract significantly inhibited the production of LasI/R mediated LasA protease, LasB elastase and chitinase with 66.52 ± 5.41, 71.26 ± 4.58 and 61.16 ± 4.28% of inhibition respectively at a concentration of 750 μg mL^-1. In addition, RhlI/R mediated production of pyocyanin, exopolysaccharides and rhamnolipids were also down-regulated by 74.71 ± 3.97, 66.41 ± 3.62 and 63.75 ± 3.76% respectively on treatment with sub-MIC concentration of D. phaseolorum SSP12. The light, fluorescence and confocal laser scanning microscopic (CLSM) analysis confirmed the significant disruption in biofilm formation. The presence of bioactive constituents such as phenyl ethylalcohol, 2, 4-di-tert-butylphenol, fenaclon, 1, 4-phenylenediacetic acid, and benzyl hydrazine in D. phaseolorum SSP12 extract was evident from Gas chromatography-mass spectrophotometric (GC-MS) analysis. From the in silico molecular docking studies, fenaclon and 2, 4-di-tert-butylphenol competitively binds to QS receptors LasR and RhlR and alters the binding of its cognate ligands and modulates the expression of virulence phenotypes. The promising anti quorum sensing efficacy of D. phaseolorum SSP12 extract suggested new avenues for development of anti-infective drugs from fungal derived metabolites to counteract the problems associated with conventional antibiotic therapies.

Effects of protein-protein interactions and ligand binding on the ion permeation in KCNQ1 potassium channel

The voltage-gated KCNQ1 potassium ion channel interacts with the type I transmembrane protein minK (KCNE1) to generate the slow delayed rectifier (I_Ks) current in the heart. Mutations in these transmembrane proteins have been linked with several heart-related issues, including long QT syndromes (LQTS), congenital atrial fibrillation, and short QT syndrome. Off-target interactions of several drugs with that of KCNQ1/KCNE1 ion channel complex have been known to cause fatal cardiac irregularities. Thus, KCNQ1/KCNE1 remains an important avenue for drug-design and discovery research. In this work, we present the structural and mechanistic details of potassium ion permeation through an open KCNQ1 structural model using the combined molecular dynamics and steered molecular dynamics simulations. We discuss the processes and key residues involved in the permeation of a potassium ion through the KCNQ1 ion channel, and how the ion permeation is affected by (i) the KCNQ1-KCNE1 interactions and (ii) the binding of chromanol 293B ligand and its derivatives into the complex. The results reveal that interactions between KCNQ1 with KCNE1 causes a pore constriction in the former, which in-turn forms small energetic barriers in the ion-permeation pathway. These findings correlate with the previous experimental reports that interactions of KCNE1 dramatically slows the activation of KCNQ1. Upon ligand-binding onto the complex, the energy-barriers along ion permeation path are more pronounced, as expected, therefore, requiring higher force in our steered-MD simulations. Nevertheless, pulling the ion when a weak blocker is bound to the channel does not necessitate high force in SMD. This indicates that our SMD simulations have been able to discern between strong and week blockers and reveal their influence on potassium ion permeation. The findings presented here will have some implications in understanding the potential off-target interactions of the drugs with the KCNQ1/KCNE1 channel that lead to cardiotoxic effects.

Guanidinium Toxins and Their Interactions with Voltage-Gated Sodium Ion Channels

Guanidinium toxins, such as saxitoxin (STX), tetrodotoxin (TTX) and their analogs, are naturally occurring alkaloids with divergent evolutionary origins and biogeographical distribution, but which share the common chemical feature of guanidinium moieties. These guanidinium groups confer high biological activity with high affinity and ion flux blockage capacity for voltage-gated sodium channels (NaV). Members of the STX group, known collectively as paralytic shellfish toxins (PSTs), are produced among three genera of marine dinoflagellates and about a dozen genera of primarily freshwater or brackish water cyanobacteria. In contrast, toxins of the TTX group occur mainly in macrozoa, particularly among puffer fish, several species of marine invertebrates and a few terrestrial amphibians. In the case of TTX and analogs, most evidence suggests that symbiotic bacteria are the origin of the toxins, although endogenous biosynthesis independent from bacteria has not been excluded. The evolutionary origin of the biosynthetic genes for STX and analogs in dinoflagellates and cyanobacteria remains elusive. These highly potent molecules have been the subject of intensive research since the latter half of the past century; first to study the mode of action of their toxigenicity, and later as tools to characterize the role and structure of NaV channels, and finally as therapeutics. Their pharmacological activities have provided encouragement for their use as therapeutants for ion channel-related pathologies, such as pain control. The functional role in aquatic and terrestrial ecosystems for both groups of toxins is unproven, although plausible mechanisms of ion channel regulation and chemical defense are often invoked. Molecular approaches and the development of improved detection methods will yield deeper understanding of their physiological and ecological roles. This knowledge will facilitate their further biotechnological exploitation and point the way towards development of pharmaceuticals and therapeutic applications.

A comprehensive structural model for the human KCNQ1/KCNE1 ion channel

The voltage-gated KCNQ1/KCNE1 potassium ion channel complex, forms the slow delayed rectifier (I_Ks) current in the heart, which plays an important role in heart signaling. The importance of KCNQ1/KCNE1 channel's function is further implicated by the linkage between loss-of-function and gain-of-function mutations in KCNQ1 or KCNE1, and long QT syndromes, congenital atrial fibrillation, and short QT syndrome. Also, KCNQ1/KCNE1 channels are an off-target for many non-cardiovascular drugs, leading to fatal cardiac irregularities. One solution to address and study the mentioned aspects of KCNQ1/KNCE1 channel would be the structural studies using a validated and accurate model. Along the same line in this study, we have used several top-notch modeling approaches to build a structural model for the open state of KCNQ1 protein, which is both accurate and compatible with available experimental data. Next, we included the KCNE1 protein components using data-driven protein-protein docking simulations, encompassing a 4:2 stoichiometry to complete the picture of the channel complex formed by these two proteins. All the protein systems generated through these processes were refined by long Molecular Dynamics simulations. The refined models were analyzed extensively to infer data about the interaction of KCNQ1 channel with its accessory KCNE1 beta subunits.