Search

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the pandemic that broke out in 2020 and continues to be the cause of massive global upheaval. Coronaviruses are positive-strand RNA viruses with a genome of ~30 kb. The genome is replicated and transcribed by RNA-dependent RNA polymerase together with accessory factors. One of the latter is the protein helicase (NSP13), which is essential for viral replication. The recently solved helicase structure revealed a tertiary structure composed of five domains. Here, we investigated NSP13 from a structural point of view, comparing its RNA-free form with the RNA-engaged form by using atomistic molecular dynamics (MD) simulations at the microsecond timescale. Structural analyses revealed conformational changes that provide insights into the contribution of the different domains, identifying the residues responsible for domain–domain interactions in both observed forms. The RNA-free system appears to be more flexible than the RNA-engaged form. This result underlies the stabilizing role of the nucleic acid and the functional core role of these domains.

The COVID-19 pandemic is caused by SARS-CoV-2. Currently, most of the research efforts towards the development of vaccines and antibodies against SARS-CoV-2 were mainly focused on the spike (S) protein, which mediates virus entry into the host cell by binding to ACE2. As the virus SARS-CoV-2 continues to spread globally, variants have emerged, characterized by multiple mutations of the S glycoprotein. Herein, we employed microsecond-long molecular dynamics simulations to study the impact of the mutations of the S glycoprotein in SARS-CoV-2 Variant of Concern 202012/01 (B.1.1.7), termed the “UK variant”, in comparison with the wild type, with the aim to decipher the structural basis of the reported increased infectivity and virulence. The simulations provided insights on the different dynamics of UK and wild-type S glycoprotein, regarding in particular the Receptor Binding Domain (RBD). In addition, we investigated the role of glycans in modulating the conformational transitions of the RBD. The overall results showed that the UK mutant experiences higher flexibility in the RBD with respect to wild type, this behavior might be correlated with the increased transmission reported for this variant. Our work also adds useful structural information on antigenic “hotspots” and epitopes targeted by neutralizing antibodies.

The conformational state of the activation loop (A-loop) is pivotal for the activity of most protein kinases. Hence, the characterization of the conformational dynamics of the A-loop is important to increase our understanding of the molecular processes related to diseases and to support the discovery of small molecule kinase inhibitors. Here, we carry out a combination of molecular dynamics (MD) and essential dynamics (ED) analyses to fully map the effects of phosphorylation, ADP, and conformation disrupting (CD) inhibitors (i.e., CD532 and MLN8054) on the dynamics of the A-loop of Aurora-A. MD revealed that the stability of the A-loop in an open conformation is enhanced by single phospho-Thr-288, while paradoxically, the presence of a second phosphorylation at Thr-287 decreases such stability and renders the A-loop more fluctuant in time and space. Moreover, we found that this post-translational modification has a significant effect on the direction of the A-loop motions. ED analysis suggests that the presence of the phosphate moiety induces the dynamics of Aurora-A to sample two distinct energy minima, instead of a single large minimum, as in unphosphorylated Aurora-A states. This observation indicates that the conformational distributions of Aurora-A with both single and double phospho-threonine modifications are remarkably different from the unphosphorylated state. In the closed states, binding of CD532 and MLN8054 inhibitors has the effect of increasing the distance of the N- and C-lobes of the kinase domain of Aurora-A, and the angle analysis between those two lobes during MD simulations showed that the N- and C-lobes are kept more open in presence of CD532, compared to MLN8054. As the A-loop is a common feature of Aurora protein kinases, our studies provide a general description of the conformational dynamics of this structure upon phosphorylation and different ligands binding.

Given the enormous social and health impact of the pandemic triggered by severe acute respiratory syndrome 2 (SARS-CoV-2), the scientific community made a huge effort to provide an immediate response to the challenges posed by Coronavirus disease 2019 (COVID-19). One of the most important proteins of the virus is an enzyme, called 3CLpro or main protease, already identified as an important pharmacological target also in SARS and Middle East respiratory syndrome virus (MERS) viruses. This protein triggers the production of a whole series of enzymes necessary for the virus to carry out its replicating and infectious activities. Therefore, it is crucial to gain a deeper understanding of 3CLpro structure and function in order to effectively target this enzyme. All-atoms molecular dynamics (MD) simulations were performed to examine the different conformational behaviors of the monomeric and dimeric form of SARS-CoV-2 3CLpro apo structure, as revealed by microsecond time scale MD simulations. Our results also shed light on the conformational dynamics of the loop regions at the entry of the catalytic site. Studying, at atomic level, the characteristics of the active site and obtaining information on how the protein can interact with its substrates will allow the design of molecules able to block the enzymatic function crucial for the virus.

Molecular dynamics (MD) simulation allows one to predict the time evolution of a system of interacting particles. It is widely used in physics, chemistry and biology to address specific questions about the structural properties and dynamical mechanisms of model systems. MD earned a great success in genome research, as it proved to be beneficial in sorting pathogenic from neutral genomic mutations. Considering their computational requirements, simulations are commonly performed on HPC computing devices, which are generally expensive and hard to administer. However, variables like the software tool used for modeling and simulation or the size of the molecule under investigation might make one hardware type or configuration more advantageous than another or even make the commodity hardware definitely suitable for MD studies. This work aims to shed lights on this aspect.

The conformational behavior of the wild-type endonucleases I-DmoI and two of its mutants has been studied in the presence and in the absence of DNA target sequences by means of extended molecular dynamics simulations. Our results show that in the absence of DNA, the three protein forms explore a similar essential conformational space, whereas when bound to the same DNA target sequence of 25 base pairs, they diversify and restrain the subspace explored. In addition, the differences in the essential subspaces explored by the residues near the catalytic site for both the bound and unbound forms are discussed in background of the experimental protein activity.

Voltage-dependent K(+) channels transfer the voltage sensor movement into gate opening or closure through an electromechanical coupling. To test functionally whether an interaction between the S4-S5 linker (L45) and the cytoplasmic end of S6 (S6(T)) constitutes this coupling, the L45 in hKv1.5 was replaced by corresponding hKv2.1 sequence. This exchange was not tolerated but could be rescued by also swapping S6(T). Exchanging both L45 and S6(T) transferred hKv2.1 kinetics to an hKv1.5 background while preserving the voltage dependence. A one-by-one residue substitution scan of L45 and S6(T) in hKv1.5 further shows that S6(T) needs to be alpha-helical and forms a "crevice" in which residues I422 and T426 of L45 reside. These residues transfer the mechanical energy onto the S6(T) crevice, whereas other residues in S6(T) and L45 that are not involved in the interaction maintain the correct structure of the coupling.

The tau protein belongs to the category of intrinsically disordered proteins, which in their native state do not have an average stable structure and fluctuate between many conformations. In its physiological state, tau helps nucleating and stabilising the microtubules in the axons of the neurons. On the other hand, the same tau is involved in the development of Alzheimer disease, when it aggregates in paired helical filaments forming fibrils, which form insoluble tangles. The beginning of the pathological aggregation of tau has been attributed to a local transition of protein portions from random coil to a β-sheet. These structures would very likely be transient; therefore, we performed a molecular dynamics simulation of tau to gather information on the existence of segments of tau endowed with a secondary structure. We combined the results of our simulation with small-angle X-ray scattering experimental data to extract from the dynamics a set of most probable conformations of tau. The analysis of these ...

The conformational behavior of the wild-type endonucleases I-DmoI and two of its mutants has been studied in the presence and in the absence of DNA target sequences by means of extended molecular dynamics simulations. Our results show that in the absence of DNA, the three protein forms explore a similar essential conformational space, whereas when bound to the same DNA target sequence of 25 base pairs, they diversify and restrain the subspace explored. In addition, the differences in the essential subspaces explored by the residues near the catalytic site for both the bound and unbound forms are discussed in background of the experimental protein activity.

KirBac1.1 and 3.1 are bacterial homologues of mammalian inward rectifier K channels. We have performed extended molecular dynamics simulations (five simulations, each of >20 ns duration) of the transmembrane domain of KirBac in two membrane environments, a palmitoyl oleoyl phosphatidylcholine bilayer and an octane slab. Analysis of these simulations has focused on the conformational dynamics of the pore-lining M2 helices, which form the cytoplasmic hydrophobic gate of the channel. Principal components analysis reveals bending of M2, with a molecular hinge at the conserved glycine (Gly134 in KirBac1.1, Gly120 in KirBac3.1). More detailed analysis reveals a dimer-of-dimers type motion. The first two eigenvectors describing the motions of M2 correspond to helix kink and swivel motions. The conformational flexibility of M2 seen in these simulations correlates with differences in M2 conformation between that seen in the X-ray structures of closed channels (KcsA and KirBac) in which the helix is undistorted, and in open channels (e.g. MthK) in which the M2 helix is kinked. Thus, the simulations, albeit on a time scale substantially shorter than that required for channel gating, suggest a gating model in which the intrinsic flexibility of M2 about a molecular hinge is coupled to conformational transitions of an intracellular 'gatekeeper' domain, the latter changing conformation in response to ligand binding.

The 37/67-kDa human laminin receptor (LamR) is a cell surface protein that interacts with molecules located in the extra-cellular matrix. In particular, interactions between LamR and laminins play a major role in mediating changes in the cellular environment that affect cell adhesion, neurite outgrowth, tumor growth and metastasis. The exact interaction mode of laminin-1 and LamR is not fully understood. Laminin-1 is thought to bind to LamR through interaction with the so-called peptide G (residues 161–180) and the C-terminal helix (residues 205–229). Here we performed 100-ns atomistic force field-based molecular dynamics simulations to explore the structure and dynamics of LamR related to laminin-1 interactions. Our main finding is that loop 188–197 in the C-terminal region is highly flexible. It undergoes a major change resulting in a conformational switch that partially solvent exposes the R180 residue in the final part of the G peptide. So, R180 could contribute to laminin-1 binding. Projection of the simulations along the first two principal components also confirms the importance of this conformational switch in the LamR. This may be a basic prerequisite to clarify the key structural determinants of the interaction of LamR with laminin-1.

The bacterial channel KirBac1.1 provides a structural homolog of mammalian inward rectifier potassium (Kir) channels. The conformational dynamics of the selectivity filter of Kir channels are of some interest in the context of possible permeation and gating mechanisms for this channel. Molecular dynamics simulations of KirBac have been performed on a 10-ns timescale, i.e., comparable to that of ion permeation. The results of five simulations (total simulation time 50 ns) based on three different initial ion configurations and two different model membranes are reported. These simulation data provide evidence for limited (

Ion channels are gated, i.e. they can switch conformation between a closed and an open state. Molecular dynamics simulations may be used to study the conformational dynamics of ion channels and of simple channel models. Simulations on model nanopores reveal that a narrow (

In recent years, increased interest in the origin of protein thermal stability has gained attention both for its possible role in understanding the forces governing the folding of a protein and for the design of new highly stable engineered biocatalysts. To study the origin of thermostability, we have performed molecular dynamics simulations of two rubredoxins, from the mesophile Clostridium pasteurianum and from the hyperthermophile Pyrococcus furiosus. The simulations were carried out at two temperatures, 300 and 373 K, for each molecule. The length of the simulations was within the range of 6-7.2 ns. The rubredoxin from the hyperthermophilic organism was more flexible than its mesophilic counterpart at both temperatures; however, the overall flexibility of both molecules at their optimal growth temperature was the same, despite 59% sequence homology. The conformational space sampled by both molecules was larger at 300 K than at 373 K. The essential dynamics analysis showed that the principal overall motions of the two molecules are significantly different. On the contrary, each molecule showed similar directions of motion at both temperatures.

This chapter focuses on homology modeling and molecular dynamics (MD) simulations studies of ion channels for the range of single cell organisms from prokaryotes to eukaryotes. The glutamate receptor channels (GluRs) share some distant homology in their transmembrane (TM) domains with K channels but possess distinct extracellular ligand-binding domains for which several structures, of both bacterial and mammalian homologs, are known. It can be seen that molecular modeling and simulations can contribute to studies of ion channels in two respects. Modeling studies enable extrapolation from experimental structures of prokaryotic ion channels to molecular models of eukaryotic homologs, thus aiding design and interpretation of, for example, mutation experiments for dissecting structure-function relationships. Ion channel structures and ion channel models may also be used as the basis of multinanosecond MD simulations. Finally, it will become increasingly important to run multiple simulations on multiple channels to allow comparative analysis of simulation results, which in turn will enable the formulation of more general hypotheses concerning the relationship between the conformational dynamics of channel proteins and their physiological functions.

From a specific definition of the roto-translational (external) and intramolecular (internal) coordinates, a constrained dynamics algorithm is derived for removing the roto-translational motions during molecular dynamics simulations, within the leap-frog integration scheme. In the paper the theoretical basis of this new method and its statistical mechanical consistency are reported, together with two applications.

The salivary antimicrobial peptide histatin-5 is able to aggregate and fuse negatively charged small unilamellar vesicles, and this fusogenic activity is selectively induced by the presence of zinc ions. Circular dichroism spectroscopy shows that histatin-5, in the presence of negatively charged vesicles and zinc ions, undergoes a conformational change leading to the stabilization of an alpha-helical secondary structure. We attribute the specific action of the zinc ions to the presence of a consensus sequence, HEXXH, located in the C-terminal functional domain of histatin-5, a recognized zinc-binding motif in many proteins. Two-dimensional proton NMR spectroscopy of histatin-5 in a trifluoroethanol/water mixture (a membrane mimetic environment) has been performed and the results analyzed by means of distance geometry and restrained molecular dynamics simulations. Our results reveal that the peptide chain, including the Zn-binding consensus sequence corresponding to residues 15-19, is in a helicoidal conformation. Comparison of the chemical shifts of the individual amino acids in histatin-5 with those recently reported in other solvents indicates that trifluoroethanol/water has a structuring capability somewhere between water and dimethyl sulfoxide. The mechanism of action of this antimicrobial peptide is discussed on the basis of its structural characteristics with particular attention to the Zn-binding motif.

Voltage-gated K+ channels (Kv) are composed of four subunits, each of which contains six trans-membrane domains (TMs), S1 through S6. The S1-S4 segments comprise the voltage-sensing domain (VS), which senses membrane potential and controls the gating of the pore domain (PD). Although still controversial, the voltage-sensing domain undergoes conformational changes within the membrane electric field, upon membrane depolarization, that is mechanically transferred via the S4-S5 linker to the intracellular gate of the channel. MD simulations of portions of PD and VS regions highlighted the importance of their flexibility for proper channel function. Nevertheless, a comprehensive description of the dynamics of both domains at atomic level has not been provided yet. Here we report the analysis of all-atom multiple molecular dynamics simulations (∼200 ns) of the entire Kv1.2/2.1 chimera, consisting of the α and β chain embedded in a 549 monomer POPC bilayer, and immersed in a box of 135K explicit SPC water molecules at 300 K. We used principal components analysis (PCA) of the Cα atomic fluctuations covariance matrix to analyze the essential subspace that characterizes the channel internal dynamics. Briefly, we observed an up to 4.5 A conformational drift of VS from its starting position. The average RMSF of the S3b-S4 domain was between 1.6 and 3.0 A. Relative to the pore region, i) the second principal component shows that all four VS domains fluctuate in a concerted manner and affect the flexibility of the intracellular gates; ii) the first principal component reveals that T1 domain moves approximately 3.5 A downwards, influencing the local structure and dynamics of the neighboring intracellular gate. Protein-lipids interactions are crucial for channel structure/function. Thus, the contributions of H-bonds and salt-bridges between channel atoms and lipid head-groups on global channel dynamics will be illustrated.

The energetic profile of an ion translated along the axis of an ion channel should reveal whether the structure corresponds to a functionally open or closed state of the channel. In this study, we explore the combined use of Poisson-Boltzmann electrostatic calculations and evaluation of van der Waals interactions between ion and pore to provide an initial appraisal of the gating state of a channel. This approach is exemplified by its application to the bacterial inward rectifier potassium channel KirBac3.1, where it reveals the closed gate to be formed by a ring of leucine (L124) side chains. We have extended this analysis to a comparative survey of gating profiles, including model hydrophobic nanopores, the nicotinic acetylcholine receptor, and a number of potassium channel structures and models. This enables us to identify three gating regimes, and to show the limitation of this computationally inexpensive method. For a (closed) gate radius of 0.4 nm < R < 0.8 nm, a hydrophobic gate may be present. For a gate radius of 0.2 nm < R < 0.4 nm, both electrostatic and van der Waals interactions will contribute to the barrier height. Below R = 0.2 nm, repulsive van der Waals interactions are likely to dominate, resulting in a sterically occluded gate. In general, the method is more useful when the channel is wider; for narrower channels, the flexibility of the protein may allow otherwise-unsurmountable energetic barriers to be overcome.

VSTx1 is a tarantula venom toxin which binds to the archaebacterial voltage-gated potassium channel KvAP. VSTx1 is thought to access the voltage sensor domain of the channel via the lipid bilayer phase. In order to understand its mode of action and implications for the mechanism of channel activation, it is important to characterize the interactions of VSTx1 with lipid bilayers. Molecular dynamics (MD) simulations (for a total simulation time in excess of 0.2 micros) have been used to explore VSTx1 localization and interactions with zwitterionic (POPC) and with anionic (POPE/POPG) lipid bilayers. In particular, three series of MD simulations have been used to explore the net drift of VSTx1 relative to the center of a bilayer, starting from different locations of the toxin. The preferred location of the toxin is at the membrane/water interface. Although there are differences between POPC and POPE/POPG bilayers, in both cases the toxin forms favorable interactions at the interface, maximizing H-bonding to lipid headgroups and to water molecules while retaining interactions with the hydrophobic core of the bilayer. A 30 ns unrestrained simulation reveals dynamic partitioning of VSTx1 into the interface of a POPC bilayer. The preferential location of VSTx1 at the interface is discussed in the context of Kv channel gating models and provides support for a mode of action in which the toxin interacts with the Kv voltage sensor "paddle" formed by the S3 and S4 helices.

Molecular modeling and simulations enable extrapolation for the structure of bacterial potassium channels to the function of their mammalian homologues. Molecular dynamics simulations have revealed the concerted single-file motion of potassium ions and water molecules through the selectivity filter of K channels and the role of filter flexibility in ion permeation and in "fast gating." Principal components analysis of extended K channel simulations suggests that hinge-bending of pore-lining M2 (or S6) helices plays a key role in K channel gating. Based on these and other simulations, a molecular model for gating of inward rectifier K channel gating is presented.

Further experimental and computational studies are required before we reach a complete structural understanding of the mechanisms of voltage-sensing and voltage-gating of Kv channels. Progress is likely to be made using a range of indirect techniques to complement the structural data that have come from X-ray crystallography. In particular, spectroscopic [7xMolecular architecture of the KvAP voltage-dependent K+ channel in a lipid bilayer. Cuello, L.G., Cortes, D.M., and Perozo, E. Science. 2004; 306: 491–495Crossref | PubMed | Scopus (173)See all References[7] and chemical modification [8xMolecular mechanism of voltage sensor movements in a potassium channel. Elliott, D.J.S., Neale, E.J., Aziz, Q., Dunham, J.P., Munsey, T.S., Hunter, M., and Sivaprasadarao, A. EMBO J. 2005; in pressSee all References[8] studies of Kv channels in situ in lipid bilayers will help to resolve the structure of the resting (closed state) channel. These methods will then have to be deployed in combination with a transmembrane voltage difference to reveal the change in sensor conformation and orientation during the transition from the closed to the open state.It is likely that computational methods will be used to integrate data from these diverse sources in a molecular mechanism. Given the proposed changes in conformation of the voltage sensor upon activation of Kv channels, it is important to characterize the intrinsic flexibility of this domain. Molecular dynamics simulations offer one possibility for exploring the conformational dynamics of the sensor domain (our unpublished data).A number of other techniques may also yield information on the location of the voltage sensor relative to the membrane. For example, toxins that interact with the sensor, for example Vstx1 [9xA membrane-access mechanism of ion channel inhibition by voltage sensor toxins from spider venom. Lee, S.Y. and MacKinnon, R. Nature. 2004; 430: 232–235Crossref | PubMed | Scopus (169)See all References[9], provide valuable probes. Their location – and hence by extension the location of the voltage-sensor – relative to a lipid bilayer may be established via molecular simulations. Similarly, electron microscopy [10xElectron microscopic analysis of KvAP voltage-dependent K+ channels in an open conformation. Jiang, Q.X., Wang, D.N., and MacKinnon, R. Nature. 2004; 430: 806–810Crossref | PubMed | Scopus (82)See all References[10] may reveal the overall shape of Kv molecules trapped in different conformational states. Thus, by combining information from these disparate sources, a complete mechanism of Kv voltage-sensing and gating may emerge.

In a previous article (Zbilut et al., Biophys J 2003;85:3544-3557), we demonstrated how an aggregation versus folding choice could be approached considering hydrophobicity distribution and charge. In this work, our aim is highlighting the mutual interaction of charge and hydrophobicity distribution in the aggregation process. Use was made of two different peptides, both derived from a transmembrane protein (amyloid precursor protein; APP), namely, Abeta(1-28) and Abeta(1-40). Abeta(1-28) has a much lower aggregation propensity than Abeta(1-40). The results obtained by means of molecular dynamics simulations show that, when submitted to the most "aggregation-prone" environment, corresponding to the isoelectric point and consequently to zero net charge, both peptides acquire their maximum flexibility, but Abeta(1-40) has a definitely higher conformational mobility than Abeta(1-28). The absence of a hydrophobic "tail," which is the most mobile part of the molecule in Abeta(1-40), is the element lacking in Abeta(1-28) for obtaining a "fully aggregating" phenotype. Our results suggest that conformational flexibility, determined by both hydrophobicity and charge effect, is the main mechanistic determinant of aggregation propensity.

Toxins that block voltage-gated potassium (Kv) channels provide a possible template for improved homology models of the Kv pore. In assessing the interactions of Kv channels and their toxins it is important to determine the dynamic flexibility of the toxins. Multiple 10 ns duration molecular dynamics simulations combined with essential dynamics analysis have been used to explore the flexibility of four different Kv channel-blocking toxins. Three toxins (Tc1, AgTx and ChTx) share a common fold. They also share a common pattern of conformational dynamics, as revealed by essential dynamics analysis of the simulation results. This suggests that some aspects of dynamic behaviour are conserved across a single protein fold class. In each of these three toxins, the residue exhibiting minimum flexibility corresponds to a conserved lysine residue that is suggested to interact with the filter domain of the channel. Thus, comparative simulations reveal functionally important conservation of molecular dynamics as well as protein fold across a family of related toxins.

Pseudomycin A is a cyclic lipodepsinonapeptide phytotoxin produced by a strain of the plant pathogenic bacterium Pseudomonas syringae. Like other members of this family of bacterial metabolites, it is characterised by a fatty acylated cyclic peptide with mixed chirality and lactonic closure. Several biological activities of Pseudomycin A are lower than those found for some of its congeners, a difference which might depend on the diverse number and distribution of charged residues in the peptide moiety. Hence, it was of interest to investigate its conformation in solution. After the complete interpretation of the two-dimensional NMR spectra, NOE data were obtained and the structure was determined by computer simulations, applying distance geometry and molecular dynamics procedures. The conformation of the large ring of Pseudomycin A in solution includes three rigid structural regions interrupted by three short flexible regions that act as hinges. The overall three-dimensional structure of the cyclic moiety is similar to that of previously studied bioactive lipodepsinonapeptides produced by other pseudomonads.

Voltage-gated potassium (Kv) channels are tetramers of 6-transmembrane domain (S1-S6) α-subunits. The activation gate that seals off the ion conducting pore is under control of the voltage-sensing domain and locates at the bundle crossing of the S6 segments. After channel opening most Kv channels display slow inactivation, a process that involves rearrangements of the selectivity filter (SF) resulting in a non-conducting channel although the S6-gate is open. Recent evidence argued for a strong coupling between inactivation and activation: after gate opening the S6 segment undergoes structural rearrangements that would be transmitted up to the level of the SF destabilizing its conformation. Substituting in Kv1.5 residue T480 that locates at the bottom of the SF by an alanine generated mutant channels which instead of inactivating displayed a second open state, characterised by slowly increasing current. Several molecular dynamics (MD) simulations showed a MD trajectory that displays a hydrophobic collapse of the central cavity with S6 dynamics decreasing the pore radius of the S6-gate. The simulations further showed that the SF did constrict in WT channels but not in the mutant. Ionic current measurements of the mutant T480A channels showed that after prolonged depolarizations - pushing the channels into the second conducting state - the activation gate could close and reopen with the channel remaining in the second conducting state. To convert the channels back to their original conducting state longer repolarizing times were needed. Thus, using the T480A pore mutant we could directly determine from ionic currents that gate closure in Kv1.5 channels does not depend on the status of the SF which, as suggested by Deutsch et al. implies the existence of a channel state with both a closed gate and inactivated selectivity filter. (Support: FWO-G025608)