RESUMEN
To design safe, selective, and effective new therapies, there must be a deep understanding of the structure and function of the drug target. One of the most difficult problems to solve has been resolution of discrete conformational states of transmembrane ion channel proteins. An example is KV11.1 (hERG), comprising the primary cardiac repolarizing current, IKr. hERG is a notorious drug anti-target against which all promising drugs are screened to determine potential for arrhythmia. Drug interactions with the hERG inactivated state are linked to elevated arrhythmia risk, and drugs may become trapped during channel closure. However, the structural details of multiple conformational states have remained elusive. Here, we guided AlphaFold2 to predict plausible hERG inactivated and closed conformations, obtaining results consistent with myriad available experimental data. Drug docking simulations demonstrated hERG state-specific drug interactions aligning well with experimental results, revealing that most drugs bind more effectively in the inactivated state and are trapped in the closed state. Molecular dynamics simulations demonstrated ion conduction that aligned with earlier studies. Finally, we identified key molecular determinants of state transitions by analyzing interaction networks across closed, open, and inactivated states in agreement with earlier mutagenesis studies. Here, we demonstrate a readily generalizable application of AlphaFold2 as a novel method to predict discrete protein conformations and novel linkages from structure to function.
RESUMEN
The authors demonstrate the feasibility of technological innovation for personalized medicine in the context of drug-induced arrhythmia. The authors use atomistic-scale structural models to predict rates of drug interaction with ion channels and make predictions of their effects in digital twins of induced pluripotent stem cell-derived cardiac myocytes. The authors construct a simplified multilayer, 1-dimensional ring model with sufficient path length to enable the prediction of arrhythmogenic dispersion of repolarization. Finally, the authors validate the computational pipeline prediction of drug effects with data and quantify drug-induced propensity to repolarization abnormalities in cardiac tissue. The technology is high throughput, computationally efficient, and low cost toward personalized pharmacologic prediction.
Asunto(s)
Arritmias Cardíacas , Células Madre Pluripotentes Inducidas , Humanos , Canales Iónicos , Miocitos Cardíacos , TecnologíaRESUMEN
Human voltage-gated sodium (hNaV) channels are responsible for initiating and propagating action potentials in excitable cells, and mutations have been associated with numerous cardiac and neurological disorders. hNaV1.7 channels are expressed in peripheral neurons and are promising targets for pain therapy. The tarantula venom peptide protoxin-II (PTx2) has high selectivity for hNaV1.7 and is a valuable scaffold for designing novel therapeutics to treat pain. Here, we used computational modeling to study the molecular mechanisms of the state-dependent binding of PTx2 to hNaV1.7 voltage-sensing domains (VSDs). Using Rosetta structural modeling methods, we constructed atomistic models of the hNaV1.7 VSD II and IV in the activated and deactivated states with docked PTx2. We then performed microsecond-long all-atom molecular dynamics (MD) simulations of the systems in hydrated lipid bilayers. Our simulations revealed that PTx2 binds most favorably to the deactivated VSD II and activated VSD IV. These state-specific interactions are mediated primarily by PTx2's residues R22, K26, K27, K28, and W30 with VSD and the surrounding membrane lipids. Our work revealed important protein-protein and protein-lipid contacts that contribute to high-affinity state-dependent toxin interaction with the channel. The workflow presented will prove useful for designing novel peptides with improved selectivity and potency for more effective and safe treatment of pain.
Asunto(s)
Canal de Sodio Activado por Voltaje NAV1.7 , Péptidos , Venenos de Araña , Humanos , Potenciales de Acción , Interneuronas , Simulación de Dinámica Molecular , Dolor , Canal de Sodio Activado por Voltaje NAV1.7/metabolismo , Venenos de Araña/metabolismo , Péptidos/metabolismoRESUMEN
The human ether-a-go-go-related gene (hERG) not only encodes a potassium-selective voltage-gated ion channel essential for normal electrical activity in the heart but is also a major drug anti-target. Genetic hERG mutations and blockage of the channel pore by drugs can cause long QT syndrome, which predisposes individuals to potentially deadly arrhythmias. However, not all hERG-blocking drugs are proarrhythmic, and their differential affinities to discrete channel conformational states have been suggested to contribute to arrhythmogenicity. We used Rosetta electron density refinement and homology modeling to build structural models of open-state hERG channel wild-type and mutant variants (Y652A, F656A, and Y652A/F656 A) and a closed-state wild-type channel based on cryo-electron microscopy structures of hERG and EAG1 channels. These models were used as protein targets for molecular docking of charged and neutral forms of amiodarone, nifekalant, dofetilide, d/l-sotalol, flecainide, and moxifloxacin. We selected these drugs based on their different arrhythmogenic potentials and abilities to facilitate hERG current. Our docking studies and clustering provided atomistic structural insights into state-dependent drug-channel interactions that play a key role in differentiating safe and harmful hERG blockers and can explain hERG channel facilitation through drug interactions with its open-state hydrophobic pockets.
RESUMEN
Cardiac function is tightly regulated by the autonomic nervous system (ANS). Activation of the sympathetic nervous system increases cardiac output by increasing heart rate and stroke volume, while parasympathetic nerve stimulation instantly slows heart rate. Importantly, imbalance in autonomic control of the heart has been implicated in the development of arrhythmias and heart failure. Understanding of the mechanisms and effects of autonomic stimulation is a major challenge because synapses in different regions of the heart result in multiple changes to heart function. For example, nerve synapses on the sinoatrial node (SAN) impact pacemaking, while synapses on contractile cells alter contraction and arrhythmia vulnerability. Here, we present a multiscale neurocardiac modelling and simulator tool that predicts the effect of efferent stimulation of the sympathetic and parasympathetic branches of the ANS on the cardiac SAN and ventricular myocardium. The model includes a layered representation of the ANS and reproduces firing properties measured experimentally. Model parameters are derived from experiments and atomistic simulations. The model is a first prototype of a digital twin that is applied to make predictions across all system scales, from subcellular signalling to pacemaker frequency to tissue level responses. We predict conditions under which autonomic imbalance induces proarrhythmia and can be modified to prevent or inhibit arrhythmia. In summary, the multiscale model constitutes a predictive digital twin framework to test and guide high-throughput prediction of novel neuromodulatory therapy. KEY POINTS: A multi-layered model representation of the autonomic nervous system that includes sympathetic and parasympathetic branches, each with sparse random intralayer connectivity, synaptic dynamics and conductance based integrate-and-fire neurons generates firing patterns in close agreement with experiment. A key feature of the neurocardiac computational model is the connection between the autonomic nervous system and both pacemaker and contractile cells, where modification to pacemaker frequency drives initiation of electrical signals in the contractile cells. We utilized atomic-scale molecular dynamics simulations to predict the association and dissociation rates of noradrenaline with the ß-adrenergic receptor. Multiscale predictions demonstrate how autonomic imbalance may increase proclivity to arrhythmias or be used to terminate arrhythmias. The model serves as a first step towards a digital twin for predicting neuromodulation to prevent or reduce disease.
Asunto(s)
Sistema Nervioso Autónomo , Corazón , Humanos , Sistema Nervioso Autónomo/fisiología , Arritmias Cardíacas , Sistema Nervioso Parasimpático , Sistema Nervioso Simpático , Frecuencia Cardíaca/fisiología , Nodo SinoatrialRESUMEN
Human voltage-gated sodium (hNaV) channels are responsible for initiating and propagating action potentials in excitable cells and mutations have been associated with numerous cardiac and neurological disorders. hNaV1.7 channels are expressed in peripheral neurons and are promising targets for pain therapy. The tarantula venom peptide protoxin-2 (PTx2) has high selectivity for hNaV1.7 and serves as a valuable scaffold to design novel therapeutics to treat pain. Here, we used computational modeling to study the molecular mechanisms of the state-dependent binding of PTx2 to hNaV1.7 voltage-sensing domains (VSDs). Using Rosetta structural modeling methods, we constructed atomistic models of the hNaV1.7 VSD II and IV in the activated and deactivated states with docked PTx2. We then performed microsecond-long all-atom molecular dynamics (MD) simulations of the systems in hydrated lipid bilayers. Our simulations revealed that PTx2 binds most favorably to the deactivated VSD II and activated VSD IV. These state-specific interactions are mediated primarily by PTx2's residues R22, K26, K27, K28, and W30 with VSD as well as the surrounding membrane lipids. Our work revealed important protein-protein and protein-lipid contacts that contribute to high-affinity state-dependent toxin interaction with the channel. The workflow presented will prove useful for designing novel peptides with improved selectivity and potency for more effective and safe treatment of pain.
RESUMEN
G protein-coupled receptors (GPCRs) represent the largest group of membrane receptors for transmembrane signal transduction. Ligand-induced activation of GPCRs triggers G protein activation followed by various signaling cascades. Understanding the structural and energetic determinants of ligand binding to GPCRs and GPCRs to G proteins is crucial to the design of pharmacological treatments targeting specific conformations of these proteins to precisely control their signaling properties. In this study, we focused on interactions of a prototypical GPCR, beta-2 adrenergic receptor (ß2AR), with its endogenous agonist, norepinephrine (NE), and the stimulatory G protein (Gs). Using molecular dynamics (MD) simulations, we demonstrated the stabilization of cationic NE, NE(+), binding to ß2AR by Gs protein recruitment, in line with experimental observations. We also captured the partial dissociation of the ligand from ß2AR and the conformational interconversions of Gs between closed and open conformations in the NE(+)-ß2AR-Gs ternary complex while it is still bound to the receptor. The variation of NE(+) binding poses was found to alter Gs α subunit (Gsα) conformational transitions. Our simulations showed that the interdomain movement and the stacking of Gsα α1 and α5 helices are significant for increasing the distance between the Gsα and ß2AR, which may indicate a partial dissociation of Gsα The distance increase commences when Gsα is predominantly in an open state and can be triggered by the intracellular loop 3 (ICL3) of ß2AR interacting with Gsα, causing conformational changes of the α5 helix. Our results help explain molecular mechanisms of ligand and GPCR-mediated modulation of G protein activation.
Asunto(s)
Subunidades alfa de la Proteína de Unión al GTP Gs , Receptores Adrenérgicos beta 2 , Ligandos , Transducción de Señal , Simulación de Dinámica Molecular , NorepinefrinaRESUMEN
The voltage-gated potassium channel, KV11.1, encoded by the human Ether-à-go-go-Related Gene (hERG), is expressed in cardiac myocytes, where it is crucial for the membrane repolarization of the action potential. Gating of the hERG channel is characterized by rapid, voltage-dependent, C-type inactivation, which blocks ion conduction and is suggested to involve constriction of the selectivity filter. Mutations S620T and S641A/T within the selectivity filter region of hERG have been shown to alter the voltage dependence of channel inactivation. Because hERG channel blockade is implicated in drug-induced arrhythmias associated with both the open and inactivated states, we used Rosetta to simulate the effects of hERG S620T and S641A/T mutations to elucidate conformational changes associated with hERG channel inactivation and differences in drug binding between the two states. Rosetta modeling of the S641A fast-inactivating mutation revealed a lateral shift of the F627 side chain in the selectivity filter into the central channel axis along the ion conduction pathway and the formation of four lateral fenestrations in the pore. Rosetta modeling of the non-inactivating mutations S620T and S641T suggested a potential molecular mechanism preventing F627 side chain from shifting into the ion conduction pathway during the proposed inactivation process. Furthermore, we used Rosetta docking to explore the binding mechanism of highly selective and potent hERG blockers - dofetilide, terfenadine, and E4031. Our structural modeling correlates well with much, but not all, existing experimental evidence involving interactions of hERG blockers with key residues in hERG pore and reveals potential molecular mechanisms of ligand interactions with hERG in an inactivated state.
RESUMEN
We report two structures of the human voltage-gated potassium channel (Kv) Kv1.3 in immune cells alone (apo-Kv1.3) and bound to an immunomodulatory drug called dalazatide (dalazatide-Kv1.3). Both the apo-Kv1.3 and dalazatide-Kv1.3 structures are in an activated state based on their depolarized voltage sensor and open inner gate. In apo-Kv1.3, the aromatic residue in the signature sequence (Y447) adopts a position that diverges 11 Å from other K+ channels. The outer pore is significantly rearranged, causing widening of the selectivity filter and perturbation of ion binding within the filter. This conformation is stabilized by a network of intrasubunit hydrogen bonds. In dalazatide-Kv1.3, binding of dalazatide to the channel's outer vestibule narrows the selectivity filter, Y447 occupies a position seen in other K+ channels, and this conformation is stabilized by a network of intersubunit hydrogen bonds. These remarkable rearrangements in the selectivity filter underlie Kv1.3's transition into the drug-blocked state.
Asunto(s)
Canal de Potasio Kv1.3/metabolismo , Canal de Potasio Kv1.3/ultraestructura , Secuencia de Aminoácidos/genética , Sitios de Unión/fisiología , Humanos , Activación del Canal Iónico/fisiología , Canal de Potasio Kv1.3/efectos de los fármacos , Potenciales de la Membrana , Microscopía Electrónica/métodos , Modelos Moleculares , Conformación Molecular , Potasio/metabolismo , Canales de Potasio/metabolismo , Canales de Potasio/ultraestructura , Canales de Potasio con Entrada de Voltaje/metabolismo , Canales de Potasio con Entrada de Voltaje/ultraestructura , Alineación de Secuencia/métodosRESUMEN
The development of induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) has been a critical in vitro advance in the study of patient-specific physiology, pathophysiology, and pharmacology. We designed a new deep learning multitask network approach intended to address the low throughput, high variability, and immature phenotype of the iPSC-CM platform. The rationale for combining translation and classification tasks is because the most likely application of the deep learning technology we describe here is to translate iPSC-CMs following application of a perturbation. The deep learning network was trained using simulated action potential (AP) data and applied to classify cells into the drug-free and drugged categories and to predict the impact of electrophysiological perturbation across the continuum of aging from the immature iPSC-CMs to the adult ventricular myocytes. The phase of the AP extremely sensitive to perturbation due to a steep rise of the membrane resistance was found to contain the key information required for successful network multitasking. We also demonstrated successful translation of both experimental and simulated iPSC-CM AP data validating our network by prediction of experimental drug-induced effects on adult cardiomyocyte APs by the latter.
Asunto(s)
Algoritmos , Aprendizaje Profundo , Técnicas Electrofisiológicas Cardíacas , Miocitos Cardíacos/fisiología , Potenciales de Acción/fisiología , Diferenciación Celular/fisiología , Simulación por Computador , Canal de Potasio ERG1/genética , Canal de Potasio ERG1/metabolismo , Fenómenos Electrofisiológicos/fisiología , Regulación de la Expresión Génica/efectos de los fármacos , Humanos , Células Madre Pluripotentes Inducidas/fisiología , Modelos Biológicos , Fenetilaminas/farmacología , Sulfonamidas/farmacologíaRESUMEN
Drug isomers may differ in their proarrhythmia risk. An interesting example is the drug sotalol, an antiarrhythmic drug comprising d- and l- enantiomers that both block the hERG cardiac potassium channel and confer differing degrees of proarrhythmic risk. We developed a multi-scale in silico pipeline focusing on hERG channel - drug interactions and used it to probe and predict the mechanisms of pro-arrhythmia risks of the two enantiomers of sotalol. Molecular dynamics (MD) simulations predicted comparable hERG channel binding affinities for d- and l-sotalol, which were validated with electrophysiology experiments. MD derived thermodynamic and kinetic parameters were used to build multi-scale functional computational models of cardiac electrophysiology at the cell and tissue scales. Functional models were used to predict inactivated state binding affinities to recapitulate electrocardiogram (ECG) QT interval prolongation observed in clinical data. Our study demonstrates how modeling and simulation can be applied to predict drug effects from the atom to the rhythm for dl-sotalol and also increased proarrhythmia proclivity of d- vs. l-sotalol when accounting for stereospecific beta-adrenergic receptor blocking.
Asunto(s)
Antagonistas Adrenérgicos beta/química , Antagonistas Adrenérgicos beta/metabolismo , Antiarrítmicos/química , Antiarrítmicos/metabolismo , Canales de Potasio Éter-A-Go-Go/metabolismo , Síndrome de QT Prolongado/metabolismo , Bloqueadores de los Canales de Potasio/química , Bloqueadores de los Canales de Potasio/metabolismo , Transducción de Señal/efectos de los fármacos , Sotalol/química , Sotalol/metabolismo , Antagonistas Adrenérgicos beta/farmacología , Antiarrítmicos/farmacología , Microscopía por Crioelectrón/métodos , Canales de Potasio Éter-A-Go-Go/antagonistas & inhibidores , Canales de Potasio Éter-A-Go-Go/química , Células HEK293 , Humanos , Simulación de Dinámica Molecular , Miocitos Cardíacos/efectos de los fármacos , Miocitos Cardíacos/metabolismo , Bloqueadores de los Canales de Potasio/farmacología , Unión Proteica/efectos de los fármacos , Sotalol/farmacología , EstereoisomerismoRESUMEN
Different mechanisms have been proposed to explain the permeation of charged compounds through lipid membranes. Overall, it is expected that an ion-induced defect permeation mechanism, where substantial membrane deformations accompany ion movement, should be dominant in thin membranes but that a solubility-diffusion mechanism, where ions partition into the membrane core with large associated dehydration energy costs, becomes dominant in thicker membranes. However, while this physical picture is intuitively reasonable, capturing the interconversion between these two permeation mechanisms in molecular dynamics (MD) simulations based on atomic models is challenging. In particular, simulations relying on nonpolarizable force fields are artificially unfavorable to the solubility-diffusion mechanism, as induced polarization of the nonpolar hydrocarbon is ignored, causing overestimated free energy costs for charged molecules to enter into this region of the membrane. In this study, all-atom MD simulations based on nonpolarizable and polarizable force fields are used to quantitatively characterize the permeation process for the arginine side chain analog methyl-guanidinium through bilayer membranes of mono-unsaturated phosphatidylcholine lipids with and without cholesterol, resulting in thicknesses spanning from â¼24 to â¼42 Å. With simulations based on a nonpolarizable force field, ion translocation can take place solely through an ion-induced defect mechanism, with free energy barriers increasing linearly from 14 to 40 kcal/mol, depending on the thickness. However, with simulations based on a polarizable force field, ion translocation is predominantly dominated by an ion-induced defect mechanism in thin membranes, which progressively converts to a solubility-diffusion mechanism as the membranes get thicker. The transition between the two mechanisms occurs at a thickness of â¼29 Å, with lipid tails of 22 or more carbon atoms. This situation appears to represent the upper limit for ion-induced defect permeation within the current polarizable models. Beyond this thickness, it becomes energetically preferable for the ion to dehydrate and partition into the membrane core-a phenomenon that cannot be captured using the nonpolarizable models. Induced electronic polarizability therefore leads not just to a shift in permeation energetics but to an interconversion between two strikingly different physical mechanisms. The result highlights the importance of induced polarizability in modeling lipid membranes.
Asunto(s)
Membrana Dobles de Lípidos , Simulación de Dinámica Molecular , Difusión , Entropía , Guanidina , Iones , TermodinámicaRESUMEN
RATIONALE: Drug-induced proarrhythmia is so tightly associated with prolongation of the QT interval that QT prolongation is an accepted surrogate marker for arrhythmia. But QT interval is too sensitive a marker and not selective, resulting in many useful drugs eliminated in drug discovery. OBJECTIVE: To predict the impact of a drug from the drug chemistry on the cardiac rhythm. METHODS AND RESULTS: In a new linkage, we connected atomistic scale information to protein, cell, and tissue scales by predicting drug-binding affinities and rates from simulation of ion channel and drug structure interactions and then used these values to model drug effects on the hERG channel. Model components were integrated into predictive models at the cell and tissue scales to expose fundamental arrhythmia vulnerability mechanisms and complex interactions underlying emergent behaviors. Human clinical data were used for model framework validation and showed excellent agreement, demonstrating feasibility of a new approach for cardiotoxicity prediction. CONCLUSIONS: We present a multiscale model framework to predict electrotoxicity in the heart from the atom to the rhythm. Novel mechanistic insights emerged at all scales of the system, from the specific nature of proarrhythmic drug interaction with the hERG channel, to the fundamental cellular and tissue-level arrhythmia mechanisms. Applications of machine learning indicate necessary and sufficient parameters that predict arrhythmia vulnerability. We expect that the model framework may be expanded to make an impact in drug discovery, drug safety screening for a variety of compounds and targets, and in a variety of regulatory processes.
Asunto(s)
Antiarrítmicos/química , Arritmias Cardíacas/tratamiento farmacológico , Cardiotoxinas/química , Simulación por Computador , Descubrimiento de Drogas/métodos , Canal de Potasio ERG1/química , Antiarrítmicos/metabolismo , Antiarrítmicos/uso terapéutico , Arritmias Cardíacas/metabolismo , Cardiotoxicidad/metabolismo , Cardiotoxicidad/prevención & control , Cardiotoxinas/efectos adversos , Cardiotoxinas/metabolismo , Descubrimiento de Drogas/tendencias , Canal de Potasio ERG1/metabolismo , Femenino , Humanos , Síndrome de QT Prolongado/tratamiento farmacológico , Síndrome de QT Prolongado/metabolismo , Aprendizaje Automático , Masculino , Moxifloxacino/química , Moxifloxacino/metabolismo , Moxifloxacino/uso terapéutico , Miocitos Cardíacos/efectos de los fármacos , Miocitos Cardíacos/fisiología , Fenetilaminas/química , Fenetilaminas/metabolismo , Fenetilaminas/uso terapéutico , Estructura Secundaria de Proteína , Sulfonamidas/química , Sulfonamidas/metabolismo , Sulfonamidas/uso terapéutico , Inhibidores de Topoisomerasa II/química , Inhibidores de Topoisomerasa II/metabolismo , Inhibidores de Topoisomerasa II/uso terapéuticoRESUMEN
The human ether-á-go-go-related gene (hERG1) channel conducts small outward K+ currents that are critical for cardiomyocyte membrane repolarization. The gain-of-function mutation N629D at the outer mouth of the selectivity filter (SF) disrupts inactivation and K+-selective transport in hERG1, leading to arrhythmogenic phenotypes associated with long-QT syndrome. Here, we combined computational electrophysiology with Markov state model analysis to investigate how SF-level gating modalities control selective cation transport in wild-type (WT) and mutant (N629D) hERG1 variants. Starting from the recently reported cryogenic electron microscopy (cryo-EM) open-state channel structure, multiple microseconds-long molecular-dynamics (MD) trajectories were generated using different cation configurations at the filter, voltages, electrolyte concentrations, and force-field parameters. Most of the K+ permeation events observed in hERG1-WT simulations occurred at microsecond timescales, influenced by the spontaneous dehydration/rehydration dynamics at the filter. The SF region displayed conductive, constricted, occluded, and dilated states, in qualitative agreement with the well-documented flickering conductance of hERG1. In line with mutagenesis studies, these gating modalities resulted from dynamic interaction networks involving residues from the SF, outer-mouth vestibule, P-helices, and S5-P segments. We found that N629D mutation significantly stabilizes the SF in a state that is permeable to both K+ and Na+, which is reminiscent of the SF in the nonselective bacterial NaK channel. Increasing the external K+ concentration induced "WT-like" SF dynamics in N629D, in qualitative agreement with the recovery of flickering currents in experiments. Overall, our findings provide an understanding of the molecular mechanisms controlling selective transport in K+ channels with a nonconventional SF sequence.
Asunto(s)
Canal de Potasio ERG1/química , Canal de Potasio ERG1/metabolismo , Secuencias de Aminoácidos , Canal de Potasio ERG1/genética , Mutación con Ganancia de Función , Humanos , Cinética , Síndrome de QT Prolongado/genética , Síndrome de QT Prolongado/metabolismo , Mutación Missense , Potasio/metabolismo , Dominios Proteicos , Estructura Secundaria de ProteínaRESUMEN
Membrane ion channels are the fundamental electrical components in the nervous system. Recent developments in X-ray crystallography and cryo-EM microscopy have revealed what these proteins look like in atomic detail but do not tell us how they function. Molecular dynamics simulations have progressed to the point that we can now simulate realistic molecular assemblies to produce quantitative calculations of the thermodynamic and kinetic quantities that control function. In this review, we summarize the state of atomistic simulation methods for ion channels to understand their conduction, activation, and drug modulation mechanisms. We are at a crossroads in atomistic simulation, where long time scale observation can provide unbiased exploration of mechanisms, supplemented by biased free energy methodologies. We illustrate the use of these approaches to describe ion conduction and selectivity in voltage-gated sodium and acid-sensing ion channels. Studies of channel gating present a significant challenge, as activation occurs on longer time scales. Enhanced sampling approaches can ensure convergence on minimum free energy pathways for activation, as illustrated here for pentameric ligand-gated ion channels that are principal to nervous system function and the actions of general anesthetics. We also examine recent studies of local anesthetic and antiepileptic drug binding to a sodium channel, revealing sites and pathways that may offer new targets for drug development. Modern simulations thus offer a range of molecular-level insights into ion channel function and modulation as a learning platform for mechanistic discovery and drug development.
Asunto(s)
Activación del Canal Iónico , Canales Iónicos/química , Membrana Celular/química , Membrana Celular/metabolismo , Humanos , Interacciones Hidrofóbicas e Hidrofílicas , Canales Iónicos/metabolismo , Modelos Químicos , Modelos Moleculares , Simulación de Dinámica Molecular , Preparaciones Farmacéuticas/química , Preparaciones Farmacéuticas/metabolismo , Conformación Proteica , TermodinámicaRESUMEN
Multi-scale computational modeling is a major branch of computational biology as evidenced by the US federal interagency Multi-Scale Modeling Consortium and major international projects. It invariably involves specific and detailed sequences of data analysis and simulation, often with multiple tools and datasets, and the community recognizes improved modularity, reuse, reproducibility, portability and scalability as critical unmet needs in this area. Scientific workflows are a well-recognized strategy for addressing these needs in scientific computing. While there are good examples if the use of scientific workflows in bioinformatics, medical informatics, biomedical imaging and data analysis, there are fewer examples in multi-scale computational modeling in general and cardiac electrophysiology in particular. Cardiac electrophysiology simulation is a mature area of multi-scale computational biology that serves as an excellent use case for developing and testing new scientific workflows. In this article, we develop, describe and test a computational workflow that serves as a proof of concept of a platform for the robust integration and implementation of a reusable and reproducible multi-scale cardiac cell and tissue model that is expandable, modular and portable. The workflow described leverages Python and Kepler-Python actor for plotting and pre/post-processing. During all stages of the workflow design, we rely on freely available open-source tools, to make our workflow freely usable by scientists.
Asunto(s)
Corazón/fisiología , Modelos Cardiovasculares , Flujo de Trabajo , Simulación por Computador , Humanos , Prueba de Estudio Conceptual , Reproducibilidad de los ResultadosRESUMEN
The human voltage-gated sodium channel, hNaV1.5, is responsible for the rapid upstroke of the cardiac action potential and is target for antiarrhythmic therapy. Despite the clinical relevance of hNaV1.5-targeting drugs, structure-based molecular mechanisms of promising or problematic drugs have not been investigated at atomic scale to inform drug design. Here, we used Rosetta structural modeling and docking as well as molecular dynamics simulations to study the interactions of antiarrhythmic and local anesthetic drugs with hNaV1.5. These calculations revealed several key drug binding sites formed within the pore lumen that can simultaneously accommodate up to two drug molecules. Molecular dynamics simulations identified a hydrophilic access pathway through the intracellular gate and a hydrophobic access pathway through a fenestration between DIII and DIV. Our results advance the understanding of molecular mechanisms of antiarrhythmic and local anesthetic drug interactions with hNaV1.5 and will be useful for rational design of novel therapeutics.
Asunto(s)
Antiarrítmicos/química , Simulación de Dinámica Molecular , Canal de Sodio Activado por Voltaje NAV1.5/química , Canales de Sodio/química , Secuencia de Aminoácidos/genética , Antiarrítmicos/uso terapéutico , Sitios de Unión , Interacciones Farmacológicas , Flecainida/química , Humanos , Lidocaína/química , Modelos Moleculares , Simulación del Acoplamiento Molecular , Canal de Sodio Activado por Voltaje NAV1.5/genética , Unión Proteica , Conformación Proteica/efectos de los fármacos , Sodio/química , Canales de Sodio/genéticaRESUMEN
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 KV and NaV 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.
Asunto(s)
Activación del Canal Iónico/fisiología , Canales Iónicos/metabolismo , Potasio/metabolismo , Sodio/metabolismo , Sitios de Unión/fisiología , Humanos , Simulación de Dinámica MolecularRESUMEN
G-protein-coupled receptors (GPCRs) are a large group of membrane-bound receptor proteins that are involved in a plethora of diverse processes (e.g., vision, hormone response). In mammals, and particularly in humans, GPCRs are involved in many signal transduction pathways and, as such, are heavily studied for their immense pharmaceutical potential. Indeed, a large fraction of drugs target various GPCRs, and drug-development is often aimed at GPCRs. Therefore, understanding the activation of GPCRs is a challenge of major importance both from fundamental and practical considerations. And yet, despite the remarkable progress in structural understanding, we still do not have a translation of the structural information to an energy-based picture. Here we use coarse-grained (CG) modeling to chart the free-energy landscape of the activation process of the ß-2 adrenergic receptor (ß2AR) as a representative GPCR. The landscape provides the needed tool for analyzing the processes that lead to activation of the receptor upon binding of the ligand (adrenaline) while limiting constitutive activation. Our results pave the way to better understand the biological mechanisms of action of the ß2AR and GPCRs, from a physical chemistry point of view rather than simply by observing the receptor's behavior physiologically.