Bases of Bacterial Sodium Channel Selectivity Among Organic Cations

Molecular Dynamics Simulations of Ion Permeation in Human Voltage-Gated Sodium Channels

The recent determination of cryo-EM structures of voltage-gated sodium (Na_v) channels has revealed many details of these proteins. However, knowledge of ionic permeation through the Na_v pore remains limited. In this work, we performed atomistic molecular dynamics (MD) simulations to study the structural features of various neuronal Na_v channels based on homology modeling of the cryo-EM structure of the human Na_v1.4 channel and, in addition, on the recently resolved configuration for Na_v1.2. In particular, single Na⁺ permeation events during standard MD runs suggest that the ion resides in the inner part of the Na_v selectivity filter (SF). On-the-fly free energy parametrization (OTFP) temperature-accelerated molecular dynamics (TAMD) was also used to calculate two-dimensional free energy surfaces (FESs) related to single/double Na⁺ translocation through the SF of the homology-based Na_v1.2 model and the cryo-EM Na_v1.2 structure, with different realizations of the DEKA filter domain. These additional simulations revealed distinct mechanisms for single and double Na⁺ permeation through the wild-type SF, which has a charged lysine in the DEKA ring. Moreover, the configurations of the ions in the SF corresponding to the metastable states of the FESs are specific for each SF motif. Overall, the description of these mechanisms gives us new insights into ion conduction in human Na_v cryo-EM-based and cryo-EM configurations that could advance understanding of these systems and how they differ from potassium and bacterial Na_v channels.

Opening of glutamate receptor channel to subconductance levels

Ionotropic glutamate receptors (iGluRs) are tetrameric ligand-gated ion channels that open their pores in response to binding of the agonist glutamate1-3. An ionic current through a single iGluR channel shows up to four discrete conductance levels (O1-O4)4-6. Higher conductance levels have been associated with an increased number of agonist molecules bound to four individual ligand-binding domains (LBDs)6-10. Here we determine structures of a synaptic complex of AMPA-subtype iGluR and the auxiliary subunit γ2 in non-desensitizing conditions with various occupancy of the LBDs by glutamate. We show that glutamate binds to LBDs of subunits B and D only after it is already bound to at least the same number of LBDs that belong to subunits A and C. Our structures combined with single-channel recordings, molecular dynamics simulations and machine-learning analysis suggest that channel opening requires agonist binding to at least two LBDs. Conversely, agonist binding to all four LBDs does not guarantee maximal channel conductance and favours subconductance states O1 and O2, with O3 and O4 being rare and not captured structurally. The lack of subunit independence and low efficiency coupling of glutamate binding to channel opening underlie the gating of synaptic complexes to submaximal conductance levels, which provide a potential for upregulation of synaptic activity.

Dissecting the Role of N-Glycan at N413 in Toll-like Receptor 3 via Molecular Dynamics Simulations

Toll-like receptor 3 (TLR3) is an endosomal receptor involved in initiating immune responses upon viral infection by directly recognizing double-stranded RNA (dsRNA). As one of the most heavily glycosylated TLR family members, the role of glycan at N413 of TLR3 in ligand recognition has been in debate for decades. Herein, to investigate the role of glycans in TLR3, specifically at amino acid residue N413, molecular dynamic simulations were performed. The loop region of LRR12 (residues 323-355), which protrudes from the dsRNA binding TLR3 lateral surface was found to be vital for interacting with dsRNA via the formation of hydrogen bonds. The glycan at N413 not only prevented dsRNA from being exposed to the bulk water during the binding process but further stabilized dsRNA in the TLR3 binding site. When N413 was in the glycosylated form, the binding free energy of TLR3 interacting with dsRNA was significantly lower than that of TLR3 in the N413 unglycosylated form. Additionally, as the glycan at N413 functioned to alter the dynamics of the dsRNA binding process, its flexibility was meanwhile influenced by dsRNA. In all, these results demonstrate that the size, length, and branch of glycan at N413 affect the thermodynamics and dynamics of TLR3 recognition with dsRNA. This study further extends our understanding of the biological role of glycans in the innate immune recognition of dsRNA by TLR3 and provides a new perspective for modulating TLR3 function.

Structural Mechanism of ω-Currents in a Mutated Kv7.2 Voltage Sensor Domain from Molecular Dynamics Simulations

Activation of voltage-gated ion channels is regulated by conformational changes of the voltage sensor domains (VSDs), four water- and ion-impermeable modules peripheral to the central, permeable pore domain. Anomalous currents, defined as ω-currents, have been recorded in response to mutations of residues on the VSD S4 helix and associated with ion fluxes through the VSDs. In humans, gene defects in the potassium channel Kv7.2 result in a broad range of epileptic disorders, from benign neonatal seizures to severe epileptic encephalopathies. Experimental evidence suggests that the R207Q mutation in S4, associated with peripheral nerve hyperexcitability, induces ω-currents at depolarized potentials, but the fine structural details are still elusive. In this work, we use atom-detailed molecular dynamics simulations and a refined model structure of the Kv7.2 VSD in the active conformation in a membrane/water environment to study the effect of R207Q and four additional mutations of proven clinical importance. Our results demonstrate that the R207Q mutant shows the most pronounced increase of hydration in the internal VSD cavity, a feature favoring the occurrence of ω-currents. Free energy and kinetics calculations of sodium permeation through the native and mutated VSD indicate as more favorable the formation of a cationic current in the latter. Overall, our simulations establish a mechanistic linkage between genetic variations and their physiological outcome, by providing a computational description that includes both thermodynamic and kinetic features of ion permeation associated with ω-currents.

Exploring the Thermodynamics of 7-Amino Actinomycin D-Induced Single-Stranded DNA Hairpin by Spectroscopic Techniques and Computational Simulations

NMR studies have indicated that the anti-tumor therapeutic agent actinomycin D (ACTD) can induce seemingly single-stranded DNA (ssDNA) oligomer 5'-CCGTT₃GTGG-3' to form a hairpin structure with tandem GT mismatches at the stem region next to a loop of three stacked thymine bases. In an effort to uncover the preference of binding sequence and to elucidate the thermodynamics properties of the binding, a combination of spectroscopic techniques and computational simulation studies was performed with d(CCGTT_nGTGG) and d(CCGAA_nGAGG) (denoted as GTT_n and GAA_n, respectively; n = 3, 5, and 7) sequences. In the presence of 7-amino actinomycin D (7AACTD), all the six oligomers formed stable hairpin structures. The GTT₅-7AACTD/GAA₅-7AACTD hairpin structure was more stable than the corresponding GTT_n-7AACTD and GAA_n-7AACTD (n = 3, 7). No significant ΔG difference was observed between GTT_n-7AACTD and GAA_n-7AACTD complexes with the same loop length. In agreement with the 7AACTD-induced hairpin stability results, the binding affinity of GTT_n and GAA_n with 7AACTD increased from n = 3 to n = 5 and then decreased when n is 7. Moreover, GTT_n and GAA_n with the same loop length showed comparable binding affinities to 7AACTD. Furthermore, molecular dynamics simulations found that van der Waals interactions between GTT_n/GAA_n and 7AACTD were the primary attractive forces for 7AACTD binding, and the electrostatic interactions between the carbonyl groups of 7AACTD and bases in the hairpin were the major unfavorable forces. These findings furthered our understanding that 7AACTD is sensitive to the loop size and sequence as well as tandem GT/GA mismatches of their deoxyribonucleic acid (DNA) targets. A deep understanding of the thermodynamics and the molecular recognition mechanism of 7AACTD with ssDNAs would further the development of ACTD-like antitumor agents.

Canales de calcio como blanco de interés farmacológico

Los canales de calcio son proteínas de membrana que constituyen la vía más importante para el ingreso del ion calcio (Ca2+) a la célula. Al abrirse, permiten el ingreso selectivo del ion, iniciando una variedad de procesos como contracción muscular, secreción endocrina y liberación de neurotransmisores, entre otros. Estas proteínas se agrupan en tres categorías de acuerdo con sus propiedades estructurales y funcionales: (i) Canales de Ca2+ operados por interacción receptor-ligando (ROCC), (ii) Canales activados por parámetros físicos (Transient Receptor Potencial, TRP) y (iii) Canales de Calcio dependientes de voltaje (VDCCs), siendo estos últimos los más estudiados debido a su presencia en células excitables. Dada la importancia de Ca2+ en la fisiología celular, los canales de Ca2+ constituyen un punto de acción farmacológica importante para múltiples tratamientos y, por tanto, son objeto de estudio para el desarrollo de nuevos fármacos. El objetivo de esta revisión es explicar la importancia de los canales de Ca2+ desde una proyección farmacológica, a partir de la exploración documental de artículos publicados hasta la fecha teniendo en cuenta temas relacionados con la estructura de los canales Ca2+, sus propiedades biofísicas, localización celular, funcionamiento y su interacción farmacológica.