scholarly journals Rapid constriction of the selectivity filter underlies C-type inactivation in the KcsA potassium channel

2018 ◽  
Vol 150 (10) ◽  
pp. 1408-1420 ◽  
Author(s):  
Jing Li ◽  
Jared Ostmeyer ◽  
Luis G. Cuello ◽  
Eduardo Perozo ◽  
Benoît Roux

C-type inactivation is a time-dependent process observed in many K+ channels whereby prolonged activation by an external stimulus leads to a reduction in ionic conduction. While C-type inactivation is thought to be a result of a constriction of the selectivity filter, the local dynamics of the process remain elusive. Here, we use molecular dynamics (MD) simulations of the KcsA channel to elucidate the nature of kinetically delayed activation/inactivation gating coupling. Microsecond-scale MD simulations based on the truncated form of the KcsA channel (C-terminal domain deleted) provide a first glimpse of the onset of C-type inactivation. We observe over multiple trajectories that the selectivity filter consistently undergoes a spontaneous and rapid (within 1–2 µs) transition to a constricted conformation when the intracellular activation gate is fully open, but remains in the conductive conformation when the activation gate is closed or partially open. Multidimensional umbrella sampling potential of mean force calculations and nonequilibrium voltage-driven simulations further confirm these observations. Electrophysiological measurements show that the truncated form of the KcsA channel inactivates faster and greater than full-length KcsA, which is consistent with truncated KcsA opening to a greater degree because of the absence of the C-terminal domain restraint. Together, these results imply that the observed kinetics underlying activation/inactivation gating reflect a rapid conductive-to-constricted transition of the selectivity filter that is allosterically controlled by the slow opening of the intracellular gate.

BIBECHANA ◽  
2020 ◽  
Vol 17 ◽  
pp. 1-12
Author(s):  
Bikash Panthi ◽  
Nurapati Pantha

Molecular Dynamics (MD) simulations of propane dimer in different solvents (water, acetonitrile and methanol) were performed by using CHARMM platform for modeling the solute and solvents. A series of Umbrella sampling MD simulations were carried out in each solvent separately and potential of mean force (PMFs) were calculated by using Weighted Histogram Analysis Method. Results show that two minima (contact minima and solvent separated minima) characterize the PMF of propane dimer in all three solvent environments. The contact minima are deeper and less sensitive to solvent environment for its position. However, significant effect in the position of second minima, solvent separated minima, was observed. Our study reveals that the interaction between propane dimer is softer in methanol and acetonitrile than in water. BIBECHANA 17 (2020) 1-12  


2018 ◽  
Vol 115 (21) ◽  
pp. 5426-5431 ◽  
Author(s):  
Alain J. Labro ◽  
D. Marien Cortes ◽  
Cholpon Tilegenova ◽  
Luis G. Cuello

The selectivity filter and the activation gate in potassium channels are functionally and structurally coupled. An allosteric coupling underlies C-type inactivation coupled to activation gating in this ion-channel family (i.e., opening of the activation gate triggers the collapse of the channel’s selectivity filter). We have identified the second Threonine residue within the TTVGYGD signature sequence of K+ channels as a crucial residue for this allosteric communication. A Threonine to Alanine substitution at this position was studied in three representative members of the K+-channel family. Interestingly, all of the mutant channels exhibited lack of C-type inactivation gating and an inversion of their allosteric coupling (i.e., closing of the activation gate collapses the channel’s selectivity filter). A state-dependent crystallographic study of KcsA-T75A proves that, on activation, the selectivity filter transitions from a nonconductive and deep C-type inactivated conformation to a conductive one. Finally, we provide a crystallographic demonstration that closed-state inactivation can be achieved by the structural collapse of the channel’s selectivity filter.


2017 ◽  
Vol 19 (18) ◽  
pp. 11017-11025 ◽  
Author(s):  
Parisa Akhshi ◽  
Gang Wu

We have applied the umbrella sampling (US) method in all-atom molecular dynamics (MD) simulations to obtain potential of mean force (PMF) profiles for ion transport through three representative G-quadruplex DNA channels: [d(TG4T)]4, [d(G3T4G4)]2, and d[G4(T4G4)3].


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Yogeeshwar Ajjugal ◽  
Kripi Tomar ◽  
D. Krishna Rao ◽  
Thenmalarchelvi Rathinavelan

AbstractBase pair mismatches in DNA can erroneously be incorporated during replication, recombination, etc. Here, the influence of A…A mismatch in the context of 5′CAA·5′TAG sequence is explored using molecular dynamics (MD) simulation, umbrella sampling MD, circular dichroism (CD), microscale thermophoresis (MST) and NMR techniques. MD simulations reveal that the A…A mismatch experiences several transient events such as base flipping, base extrusion, etc. facilitating B–Z junction formation. A…A mismatch may assume such conformational transitions to circumvent the effect of nonisostericity with the flanking canonical base pairs so as to get accommodated in the DNA. CD and 1D proton NMR experiments further reveal that the extent of B–Z junction increases when the number of A…A mismatch in d(CAA)·d(T(A/T)G) increases (1–5). CD titration studies of d(CAA)·d(TAG)n=5 with the hZαADAR1 show the passive binding between the two, wherein, the binding of protein commences with B–Z junction recognition. Umbrella sampling simulation indicates that the mismatch samples anti…+ syn/+ syn…anti, anti…anti & + syn…+ syn glycosyl conformations. The concomitant spontaneous transitions are: a variety of hydrogen bonding patterns, stacking and minor or major groove extrahelical movements (with and without the engagement of hydrogen bonds) involving the mismatch adenines. These transitions frequently happen in anti…anti conformational region compared with the other three regions as revealed from the lifetime of these states. Further, 2D-NOESY experiments indicate that the number of cross-peaks diminishes with the increasing number of A…A mismatches implicating its dynamic nature. The spontaneous extrahelical movement seen in A…A mismatch may be a key pre-trapping event in the mismatch repair due to the accessibility of the base(s) to the sophisticated mismatch repair machinery.


2021 ◽  
Vol 0 (0) ◽  
Author(s):  
Rajendra P. Koirala ◽  
Rudramani Pokhrel ◽  
Prabin Baral ◽  
Purushottam B. Tiwari ◽  
Prem P. Chapagain ◽  
...  

Abstract Methylation induced DNA base-pairing damage is one of the major causes of cancer. O6-alkylguanine-DNA alkyltransferase (AGT) is considered a demethylation agent of the methylated DNA. Structural investigations with thermodynamic properties of the AGT-DNA complex are still lacking. In this report, we modeled two catalytic states of AGT-DNA interactions and an AGT-DNA covalent complex and explored structural features using molecular dynamics (MD) simulations. We utilized the umbrella sampling method to investigate the changes in the free energy of the interactions in two different AGT-DNA catalytic states, one with methylated GUA in DNA and the other with methylated CYS145 in AGT. These non-covalent complexes represent the pre- and post-repair complexes. Therefore, our study encompasses the process of recognition, complex formation, and separation of the AGT and the damaged (methylated) DNA base. We believe that the use of parameters for the amino acid and nucleotide modifications and for the protein-DNA covalent bond will allow investigations of the DNA repair mechanism as well as the exploration of cancer therapeutics targeting the AGT-DNA complexes at various functional states as well as explorations via stabilization of the complex.


2007 ◽  
Vol 129 (5) ◽  
pp. 403-418 ◽  
Author(s):  
Gyorgy Panyi ◽  
Carol Deutsch

Slow inactivation involves a local rearrangement of the outer mouth of voltage-gated potassium channels, but nothing is known regarding rearrangements in the cavity between the activation gate and the selectivity filter. We now report that the cavity undergoes a conformational change in the slow-inactivated state. This change is manifest as altered accessibility of residues facing the aqueous cavity and as a marked decrease in the affinity of tetraethylammonium for its internal binding site. These findings have implications for global alterations of the channel during slow inactivation and putative coupling between activation and slow-inactivation gates.


PLoS ONE ◽  
2021 ◽  
Vol 16 (2) ◽  
pp. e0246187
Author(s):  
Tianyi Jin ◽  
Samarthaben J. Patel ◽  
Reid C. Van Lehn

Quorum sensing (QS) is a bacterial communication process mediated by both native and non-native small-molecule quorum sensing modulators (QSMs), many of which have been synthesized to disrupt QS pathways. While structure-activity relationships have been developed to relate QSM structure to the activation or inhibition of QS receptors, less is known about the transport mechanisms that enable QSMs to cross the lipid membrane and access intracellular receptors. In this study, we used atomistic MD simulations and an implicit solvent model, called COSMOmic, to analyze the partitioning and translocation of QSMs across lipid bilayers. We performed umbrella sampling at atomistic resolution to calculate partitioning and translocation free energies for a set of naturally occurring QSMs, then used COSMOmic to screen the water-membrane partition and translocation free energies for 50 native and non-native QSMs that target LasR, one of the LuxR family of quorum-sensing receptors. This screening procedure revealed the influence of systematic changes to head and tail group structures on membrane partitioning and translocation free energies at a significantly reduced computational cost compared to atomistic MD simulations. Comparisons with previously determined QSM activities suggest that QSMs that are least likely to partition into the bilayer are also less active. This work thus demonstrates the ability of the computational protocol to interrogate QSM-bilayer interactions which may help guide the design of new QSMs with engineered membrane interactions.


2021 ◽  
Author(s):  
Fréderic Célerse ◽  
Theo Jaffrelot-Inizan ◽  
Louis Lagardère ◽  
Olivier Adjoua ◽  
Pierre Monmarché ◽  
...  

We detail a novel multi-level enhanced sampling strategy grounded on Gaussian accelerated Molecular Dynamics (GaMD). First, we propose a GaMD multi-GPUs-accelerated implementation within the Tinker-HP molecular dynamics package. We then introduce the new "dual-water" mode and its use with the flexible AMOEBA polarizable force field. By adding harmonic boosts to the water stretching and bonding terms, it accelerates the solvent-solute interactions while enabling speedups thanks to the use of fast multiple--timestep integrators. To further reduce time-to-solution, we couple GaMD to Umbrella Sampling (US). The GaMD—US/dual-water approach is tested on the 1D Potential of Mean Force (PMF) of the CD2-CD58 system (168000 atoms) allowing the AMOEBA PMF to converge within 1 kcal/mol of the experimental value. Finally, Adaptive Sampling (AS) is added enabling AS-GaMD capabilities but also the introduction of the new Adaptive Sampling--US--GaMD (ASUS--GaMD) scheme. The highly parallel ASUS--GaMD setup decreases time to convergence by respectively 10 and 20 compared to GaMD--US and US.


2021 ◽  
Author(s):  
Fréderic Célerse ◽  
Theo Jaffrelot-Inizan ◽  
Louis Lagardère ◽  
Olivier Adjoua ◽  
Pierre Monmarché ◽  
...  

We introduce a novel multi-level enhanced sampling strategy grounded on Gaussian accelerated Molecular Dynamics (GaMD). First, we propose a GaMD multi-GPUs -accelerated implementation within Tinker-HP. For the specific use with the flexible AMOEBA polarizable force field (PFF), we introduce the new "dual–water" GaMD mode. By adding harmonic boosts to the water stretching and bonding terms, it accelerates the solvent-solute interactions while enabling speedups with fast multiple–timestep integrators. To further reduce time-to-solution, we couple GaMD to Umbrella Sampling (US). The GaMD—US/dual–water approach is tested on the 1D Potential of Mean Force (PMF) of the CD2–CD58 system (168000 atoms) allowing the AMOEBA PMF to converge within 1 kcal/mol of the experimental value. Finally, Adaptive Sampling (AS) is added enabling AS–GaMD capabilities but also the introduction of the new Adaptive Sampling–US–GaMD (ASUS–GaMD) scheme. The highly parallel ASUS–GaMD setup decreases time to convergence by respectively 10 and 20 compared to GaMD–US and US.


Sign in / Sign up

Export Citation Format

Share Document