The binding of heparin to spike glycoprotein inhibits SARS-CoV-2 infection by three mechanisms

Heparin, a naturally occurring glycosaminoglycan, has been found to have antiviral activity against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative virus of COVID-19. To elucidate the mechanistic basis for the antiviral activity of heparin, we investigated the binding of heparin to the SARS-CoV-2 spike glycoprotein by means of sliding window docking, molecular dynamics simulations, and biochemical assays. Our simulations show that heparin binds at long, positively charged patches on the spike glycoprotein, thereby masking basic residues of both the receptor-binding domain (RBD) and the multifunctional S1/S2 site. Biochemical experiments corroborated the simulation results, showing that heparin inhibits the furin-mediated cleavage of spike by binding to the S1/S2 site. Our simulations showed that heparin can act on the hinge region responsible for motion of the RBD between the inactive closed and active open conformations of the spike glycoprotein. In simulations of the closed spike homotrimer, heparin binds the RBD and the N-terminal domain of two adjacent spike subunits and hinders opening. In simulations of open spike conformations, heparin induces stabilization of the hinge region and a change in RBD motion. Our results indicate that heparin can inhibit SARS-CoV-2 infection by three mechanisms: by allosterically hindering binding to the host cell receptor, by directly competing with binding to host heparan sulfate proteoglycan coreceptors, and by preventing spike cleavage by furin. Furthermore, these simulations provide insights into how host heparan sulfate proteoglycans can facilitate viral infection. Our results will aid the rational optimization of heparin derivatives for SARS-CoV-2 antiviral therapy.     

Elucidation of interactions regulating conformational stability and dynamics of SARS-CoV-2 S-protein

The ongoing COVID-19 pandemic caused by the new coronavirus, SARS-CoV-2, calls for urgent developments of vaccines and antiviral drugs. The spike protein of SARS-CoV-2 (S-protein), which consists of trimeric polypeptide chains with glycosylated residues on the surface, triggers the virus entry into a host cell. Extensive structural and functional studies on this protein have rapidly advanced our understanding of the S-protein structure at atomic resolutions, although most of these structural studies overlook the effect of glycans attached to the S-protein on the conformational stability and functional motions between the inactive down and active up forms. Here, we performed all-atom molecular dynamics simulations of both down and up forms of a fully glycosylated S-protein in solution as well as targeted molecular dynamics simulations between them to elucidate key interdomain interactions for stabilizing each form and inducing the large-scale conformational transitions. The residue-level interaction analysis of the simulation trajectories detects distinct amino acid residues and N-glycans as determinants on conformational stability of each form. During the conformational transitions between them, interdomain interactions mediated by glycosylated residues are switched to play key roles on the stabilization of another form. Electrostatic interactions, as well as hydrogen bonds between the three receptor binding domains, work as driving forces to initiate the conformational transitions toward the active form. This study sheds light on the mechanisms underlying conformational stability and functional motions of the S-protein, which are relevant for vaccine and antiviral drug developments.     

Interplay between structural rigidity and electrostatic interactions in the ligand binding domain of GluR2

Using molecular dynamics (MD) simulations, computational protein modifications, and a novel theoretical methodology that determines structural rigidity/flexibility (the FIRST algorithm), we investigate how molecular structure and dynamics of the glutamate receptor ligand binding domain (GluR2 S1S2) facilitate its conformational transition. S1S2 is a two-lobe protein, which undergoes a cleft closure conformational transition upon binding an agonist in the cleft between the two lobes; hence it is expected that the mechanism of this conformational transition can be characterized as a hinge-type. However, in the rigidity analysis one lobe of the protein is identified as a single rigid cluster while the other one is structurally flexible, inconsistent with a presumed mechanical hinge mechanism. Instead, we characterize the cleft-closing transition as a load and lock mechanism. We find that when two cross-cleft hydrogen bonds are disrupted the protein undergoes a rapid cleft opening transition. At the same time, the dynamical behavior of the cleft in the presence of the glutamate ligand is only weakly affected by the S652 peptide bond in its flipped conformation observed in the crystal structure. The residue E705 plays significant role in stabilization of the closed conformation via electrostatic interactions. The presence of the E705-K730 salt bridge seems to correlate strongly withthe cleft opening transition in the MD simulations.     

Population Dynamics and Structural Effects at Short and Long Range Support the Hypothesis of the Selective Advantage of the G614 SARS-CoV-2 Spike Variant

SARS-CoV-2 epidemics quickly propagated worldwide, sorting virus genomic variants in newly established propagules of infections. Stochasticity in transmission within and between countries or an actual selective advantage could explain the global high frequency reached by some genomic variants. Using statistical analyses, demographic reconstructions, and molecular dynamics simulations, we show that the globally invasive G614 spike variant 1) underwent a significant demographic expansion in most countries explained neither by stochastic effects nor by overrepresentation in clinical samples, 2) increases the spike S1/S2 furin-like site conformational plasticity (short-range effect), and 3) modifies the internal motion of the receptor-binding domain affecting its cross-connection with other functional domains (long-range effect). Our results support the hypothesis of a selective advantage at the basis of the spread of the G614 variant, which we suggest may be due to structural modification of the spike protein at the S1/S2 proteolytic site, and provide structural information to guide the design of variant-specific drugs.     

How the headpiece hinge angle is opened: New insights into the dynamics of integrin activation

How the integrin head transitions to the high-affinity conformation is debated. Although experiments link activation with the opening of the hinge angle between the betaA and hybrid domains in the ligand-binding headpiece, this hinge is closed in the liganded alpha(v)beta3 integrin crystal structure. We replaced the RGD peptide ligand of this structure with the 10th type III fibronectin module (FnIII10) and discovered through molecular dynamics (MD) equilibrations that when the conformational constraints of the leg domains are lifted, the betaA/hybrid hinge opens spontaneously. Together with additional equilibrations on the same nanosecond timescale in which small structural variations impeded hinge-angle opening, these simulations allowed us to identify the allosteric pathway along which ligand-induced strain propagates via elastic distortions of the alpha1 helix to the betaA/hybrid domain hinge. Finally, we show with steered MD how force accelerates hinge-angle opening along the same allosteric pathway. Together with available experimental data, these predictions provide a novel framework for understanding integrin activation.     

The flexibility of ACE2 in the context of SARS-CoV-2 infection

The coronavirus disease 2019 (COVID-19) pandemic has swept over the world in the past months, causing significant loss of life and consequences to human health. Although numerous drug and vaccine development efforts are underway, there are many outstanding questions on the mechanism of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) viral association to angiotensin-converting enzyme 2 (ACE2), its main host receptor, and host cell entry. Structural and biophysical studies indicate some degree of flexibility in the viral extracellular spike glycoprotein and at the receptor-binding domain (RBD)-receptor interface, suggesting a role in infection. Here, we perform explicitly solvated, all-atom, molecular dynamics simulations of the glycosylated, full-length, membrane-bound ACE2 receptor in both an apo and spike RBD-bound state to probe the intrinsic dynamics of the ACE2 receptor in the context of the cell surface. A large degree of fluctuation in the full-length structure is observed, indicating hinge bending motions at the linker region connecting the head to the transmembrane helix while still not disrupting the ACE2 homodimer or ACE2-RBD interfaces. This flexibility translates into an ensemble of ACE2 homodimer conformations that could sterically accommodate binding of the spike trimer to more than one ACE2 homodimer and suggests a mechanical contribution of the host receptor toward the large spike conformational changes required for cell fusion. This work presents further structural and functional insights into the role of ACE2 in viral infection that can potentially be exploited for the rational design of effective SARS-CoV-2 therapeutics.     

Engineered disulfide reveals structural dynamics of locked SARS-CoV-2 spike

The spike (S) protein of SARS-CoV-2 has been observed in three distinct pre-fusion conformations: locked, closed and open. Of these, the function of the locked conformation remains poorly understood. Here we engineered a SARS-CoV-2 S protein construct “S-R/x3” to arrest SARS-CoV-2 spikes in the locked conformation by a disulfide bond. Using this construct we determined high-resolution structures confirming that the x3 disulfide bond has the ability to stabilize the otherwise transient locked conformations. Structural analyses reveal that wild-type SARS-CoV-2 spike can adopt two distinct locked-1 and locked-2 conformations. For the D614G spike, based on which all variants of concern were evolved, only the locked-2 conformation was observed. Analysis of the structures suggests that rigidified domain D in the locked conformations interacts with the hinge to domain C and thereby restrains RBD movement. Structural change in domain D correlates with spike conformational change. We propose that the locked-1 and locked-2 conformations of S are present in the acidic high-lipid cellular compartments during virus assembly and egress. In this model, release of the virion into the neutral pH extracellular space would favour transition to the closed or open conformations. The dynamics of this transition can be altered by mutations that modulate domain D structure, as is the case for the D614G mutation, leading to changes in viral fitness. The S-R/x3 construct provides a tool for the further structural and functional characterization of the locked conformations of S, as well as how sequence changes might alter S assembly and regulation of receptor binding domain dynamics.     

Ca2+-dependent mechanism of membrane insertion and destabilization by the SARS-CoV-2 fusion peptide

Cell penetration after recognition of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus by the ACE2 receptor and the fusion of its viral envelope membrane with cellular membranes are the early steps of infectivity. A region of the Spike protein of the virus, identified as the “fusion peptide” (FP), is liberated at its N-terminal site by a specific cleavage occurring in concert with the interaction of the receptor-binding domain of the Spike. Studies have shown that penetration is enhanced by the required binding of Ca²⁺ ions to the FPs of coronaviruses, but the mechanisms of membrane insertion and destabilization remain unclear. We have predicted the preferred positions of Ca²⁺ binding to the SARS-CoV-2-FP, the role of Ca²⁺ ions in mediating peptide-membrane interactions, the preferred mode of insertion of the Ca²⁺-bound SARS-CoV-2-FP, and consequent effects on the lipid bilayer from extensive atomistic molecular dynamics simulations and trajectory analyses. In a systematic sampling of the interactions of the Ca²⁺-bound peptide models with lipid membranes, SARS-CoV-2-FP penetrated the bilayer and disrupted its organization only in two modes involving different structural domains. In one, the hydrophobic residues F833/I834 from the middle region of the peptide are inserted. In the other, more prevalent mode, the penetration involves residues L822/F823 from the LLF motif, which is conserved in CoV-2-like viruses, and is achieved by the binding of Ca²⁺ ions to the D830/D839 and E819/D820 residue pairs. FP penetration is shown to modify the molecular organization in specific areas of the bilayer, and the extent of membrane binding of the SARS-CoV-2 FP is significantly reduced in the absence of Ca²⁺ ions. These findings provide novel mechanistic insights regarding the role of Ca²⁺ in mediating SARS-CoV-2 fusion and provide a detailed structural platform to aid the ongoing efforts in rational design of compounds to inhibit SARS-CoV-2 cell entry.     

Modeling coronavirus spike protein dynamics: implications for immunogenicity and immune escape

The ongoing COVID-19 pandemic is a global public health emergency requiring urgent development of efficacious vaccines. While concentrated research efforts have focused primarily on antibody-based vaccines that neutralize SARS-CoV-2, and several first-generation vaccines have either been approved or received emergency use authorization, it is forecasted that COVID-19 will become an endemic disease requiring updated second-generation vaccines. The SARS-CoV-2 surface spike (S) glycoprotein represents a prime target for vaccine development because antibodies that block viral attachment and entry, i.e., neutralizing antibodies, bind almost exclusively to the receptor-binding domain. Here, we develop computational models for a large subset of S proteins associated with SARS-CoV-2, implemented through coarse-grained elastic network models and normal mode analysis. We then analyze local protein domain dynamics of the S protein systems and their thermal stability to characterize structural and dynamical variability among them. These results are compared against existing experimental data and used to elucidate the impact and mechanisms of SARS-CoV-2 S protein mutations and their associated antibody binding behavior. We construct a SARS-CoV-2 antigenic map and offer predictions about the neutralization capabilities of antibody and S mutant combinations based on protein dynamic signatures. We then compare SARS-CoV-2 S protein dynamics to SARS-CoV and MERS-CoV S proteins to investigate differing antibody binding and cellular fusion mechanisms that may explain the high transmissibility of SARS-CoV-2. The outbreaks associated with SARS-CoV, MERS-CoV, and SARS-CoV-2 over the last two decades suggest that the threat presented by coronaviruses is ever-changing and long term. Our results provide insights into the dynamics-driven mechanisms of immunogenicity associated with coronavirus S proteins and present a new, to our knowledge, approach to characterize and screen potential mutant candidates for immunogen design, as well as to characterize emerging natural variants that may escape vaccine-induced antibody responses.     

Pores formed by the nicotinic receptor m2delta Peptide: a molecular dynamics simulation study

The M2delta peptide self-assembles to form a pentameric bundle of transmembrane alpha-helices that is a model of the pore-lining region of the nicotinic acetylcholine receptor. Long (>15 ns) molecular dynamics simulations of a model of the M2delta(5) bundle in a POPC bilayer have been used to explore the conformational dynamics of the channel assembly. On the timescale of the simulation, the bundle remains relatively stable, with the polar pore-lining side chains remaining exposed to the lumen of the channel. Fluctuations at the helix termini, and in the helix curvature, result in closing/opening transitions at both mouths of the channel, on a timescale of approximately 10 ns. On average, water within the pore lumen diffuses approximately 4x more slowly than water outside the channel. Examination of pore water trajectories reveals both single-file and path-crossing regimes to occur at different times within the simulation.     

Effect of clinical isolate or cleavage site mutations in the SARS-CoV-2 spike protein on protein stability,cleavage, and cell–cell fusion

The trimeric severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike protein (S) is the sole viral protein responsible for both viral binding to a host cell and the membrane fusion event needed for cell entry. In addition to facilitating fusion needed for viral entry, S can also drive cell–cell fusion, a pathogenic effect observed in the lungs of SARS-CoV-2–infected patients. While several studies have investigated S requirements involved in viral particle entry, examination of S stability and factors involved in S cell–cell fusion remain limited. A furin cleavage site at the border between the S1 and S2 subunits (S1/S2) has been identified, along with putative cathepsin L and transmembrane serine protease 2 cleavage sites within S2. We demonstrate that S must be processed at the S1/S2 border in order to mediate cell–cell fusion and that mutations at potential cleavage sites within the S2 subunit alter S processing at the S1/S2 border, thus preventing cell–cell fusion. We also identify residues within the internal fusion peptide and the cytoplasmic tail that modulate S-mediated cell–cell fusion. In addition, we examined S stability and protein cleavage kinetics in a variety of mammalian cell lines, including a bat cell line related to the likely reservoir species for SARS-CoV-2, and provide evidence that proteolytic processing alters the stability of the S trimer. This work therefore offers insight into S stability, proteolytic processing, and factors that mediate S cell–cell fusion, all of which help give a more comprehensive understanding of this high-profile therapeutic target.     

Prefusion spike protein conformational changes are slower in SARS-CoV-2 than in SARS-CoV-1

Within the last 2 decades, severe acute respiratory syndrome coronaviruses 1 and 2 (SARS-CoV-1 and SARS-CoV-2) have caused two major outbreaks; yet, for reasons not fully understood, the coronavirus disease 2019 pandemic caused by SARS-CoV-2 has been significantly more widespread than the 2003 SARS epidemic caused by SARS-CoV-1, despite striking similarities between these two viruses. The SARS-CoV-1 and SARS-CoV-2 spike proteins, both of which bind to host cell angiotensin-converting enzyme 2, have been implied to be a potential source of their differential transmissibility. However, the mechanistic details of prefusion spike protein binding to angiotensin-converting enzyme 2 remain elusive at the molecular level. Here, we performed an extensive set of equilibrium and nonequilibrium microsecond-level all-atom molecular dynamics simulations of SARS-CoV-1 and SARS-CoV-2 prefusion spike proteins to determine their differential dynamic behavior. Our results indicate that the active form of the SARS-CoV-2 spike protein is more stable than that of SARS-CoV-1 and the energy barrier associated with the activation is higher in SARS-CoV-2. These results suggest that not only the receptor-binding domain but also other domains such as the N-terminal domain could play a crucial role in the differential binding behavior of SARS-CoV-1 and SARS-CoV-2 spike proteins.     

The polybasic insert,the RBD of the SARS-CoV-2 spike protein,and the feline coronavirus – evolved or yet to evolve

Recent research on the SARS-CoV-2 pandemic has exploded around the furin-cleavable polybasic insert PRRAR↓S, found within the spike protein. The insert and the receptor-binding domain, (RBD), are vital clues in the Sherlock Holmes-like investigation into the origin of the virus and in its zoonotic crossover. Based on comparative analysis of the whole genome and the sequence features of the insert and the RBD domain, the bat and the pangolin have been proposed as very likely intermediary hosts. In this study, using the various databases, in-house developed tools, sequence comparisons, structure-guided docking, and molecular dynamics simulations, we cautiously present a fresh, theoretical perspective on the SARS-CoV-2 virus activation and its intermediary host. They are a) the SARS-CoV-2 has not yet acquired a fully optimal furin binding site or this seemingly less optimal sequence, PRRARS, has been selected for survival; b) in structural models of furin complexed with peptides, PRRAR↓S binds less well and with distinct differences as compared to the all basic RRKRR↓S; c) these differences may be exploited for the design of virus-specific inhibitors; d) the novel polybasic insert of SARS-CoV-2 may be promiscuous enough to be cleaved by multiple enzymes of the human airway epithelium and tissues which may explain its unexpected broad tropism; e) the RBD domain of the feline coronavirus spike protein carries residues that are responsible for high-affinity binding of the SARS-CoV-2 to the ACE 2 receptor; f) en route zoonotic transfer, the virus may have passed through the domestic cat whose very human-like ACE2 receptor and furin may have played some role in optimizing the traits required for zoonotic transfer.     

Elastase-mediated Activation of the Severe Acute Respiratory Syndrome Coronavirus Spike Protein at Discrete Sites within the S2 Domain

Proteolytic priming is a common method of controlling the activation of membrane fusion mediated by viral glycoproteins. The severe acute respiratory syndrome coronavirus spike protein (SARS-CoV S) can be primed by a variety of host cell proteases, with proteolytic cleavage occurring both as the S1/S2 boundary and adjacent to a fusion peptide in the S2 domain. Here, we studied the priming of SARS-CoV S by elastase and show an important role for residue Thr⁷⁹⁵ in the S2 domain. A series of alanine mutants were generated in the vicinity of the S2 cleavage site, with the goal of examining elastase-mediated cleavage within S2. Both proteolytic cleavage and fusion activation were modulated by altering the cleavage site position. We propose a novel mechanism whereby SARS-CoV fusion protein function can be controlled by spatial regulation of the proteolytic priming site, with important implications for viral pathogenesis.     

A second SARS-CoV S2 glycoprotein internal membrane-active peptide. Biophysical characterization and membrane interaction

The severe acute respiratory syndrome coronavirus (SARS-CoV) envelope spike (S) glycoprotein, a class I viral fusion protein, is responsible for the fusion between the membranes of the virus and the target cell. The S2 domain of protein S has been suggested to have two fusion peptides, one located at its N-terminus, downstream of the furin cleavage, and another, more internal, located immediately upstream of the HR1. Therefore, we have carried out a study of the binding and interaction with model membranes of a peptide corresponding to segment 873-888 of the SARS-CoV S glycoprotein, peptide SARS IFP, as well as the structural changes taking place in both the phospholipid and the peptide induced by the binding of the peptide to the membrane. We demonstrate that SARS IFP peptide binds to and interacts with phospholipid model membranes and shows a higher affinity for negatively charged phospholipids than for zwitterionic ones. SARS IFP peptide specifically decreases the mobility of the phospholipid acyl chains of negatively charged phospholipids and adopts different conformations in the membrane depending upon their composition. These data support its role in SARS-mediated membrane fusion and suggest that the regions where this peptide resides might assist the fusion peptide and/or the pretransmembrane segment of the SARS-CoV spike glycoprotein in the fusion process.     

Markov state modeling of membrane transport proteins

The flux of ions and molecules in and out of the cell is vital for maintaining the basis of various biological processes. The permeation of substrates across the cellular membrane is mediated through the function of specialized integral membrane proteins commonly known as membrane transporters. These proteins undergo a series of structural rearrangements that allow a primary substrate binding site to be accessed from either side of the membrane at a given time. Structural insights provided by experimentally resolved structures of membrane transporters have aided in the biophysical characterization of these important molecular drug targets. However, characterizing the transitions between conformational states remains challenging to achieve both experimentally and computationally. Though molecular dynamics simulations are a powerful approach to provide atomistic resolution of protein dynamics, a recurring challenge is its ability to efficiently obtain relevant timescales of large conformational transitions as exhibited in transporters. One approach to overcome this difficulty is to adaptively guide the simulation to favor exploration of the conformational landscape, otherwise known as adaptive sampling. Furthermore, such sampling is greatly benefited by the statistical analysis of Markov state models. Historically, the use of Markov state models has been effective in quantifying slow dynamics or long timescale behaviors such as protein folding. Here, we review recent implementations of adaptive sampling and Markov state models to not only address current limitations of molecular dynamics simulations, but to also highlight how Markov state modeling can be applied to investigate the structure–function mechanisms of large, complex membrane transporters.     

Molecular adaptive evolution of SARS-COV-2 spike protein in Saudi Arabia

The sequences of SARS-CoV-2 spike (S) from Saudi Arabia along with SARS-CoV and bat SARS-like CoVs were obtained. Positive selection analysis and secondary structure investigation of spike sequences were performed. Adaptive molecular evolution was observed in SARS-CoV-2 displayed by positive selection pressure at N-terminal domain (NTD; codons 41, 163, 174 and 218), Receptor binding domain (RBD; codons 378 and 404) and S1/S2 Cleavage site (codon 690). Furthermore, the spike protein secondary structure depicted by the homo-trimer structure showed a high similarity between Saudi SARS-CoV-2 isolate and the parental strain (bat SL-COVZC45). Despite the high similarity depicted in the spike sequence model alignment, it displayed a significant difference when each chain was treated solely owing to 7 motif differences in the three composing chains. In addition, SARS-CoV-2 S trimer model uncovered the presence of N-acetyl glucosamine ligands. Eventually, 3C-like proteinase cleavage site was observed in S2 domain could be used as a site for drug discovery. Genetics and molecular evolutionary facts are useful for assessment of evolution, host adaptation and epidemic patterns ultimately helpful for adaptation of control strategies.     

SARS-CoV-2 variants with mutations at the S1/S2 cleavage site are generated in vitro during propagation in TMPRSS2-deficient cells

The spike (S) protein of Severe Acute Respiratory Syndrome-Coronavirus-2 (SARS-CoV-2) binds to a host cell receptor which facilitates viral entry. A polybasic motif detected at the cleavage site of the S protein has been shown to broaden the cell tropism and transmissibility of the virus. Here we examine the properties of SARS-CoV-2 variants with mutations at the S protein cleavage site that undergo inefficient proteolytic cleavage. Virus variants with S gene mutations generated smaller plaques and exhibited a more limited range of cell tropism compared to the wild-type strain. These alterations were shown to result from their inability to utilize the entry pathway involving direct fusion mediated by the host type II transmembrane serine protease, TMPRSS2. Notably, viruses with S gene mutations emerged rapidly and became the dominant SARS-CoV-2 variants in TMPRSS2-deficient cells including Vero cells. Our study demonstrated that the S protein polybasic cleavage motif is a critical factor underlying SARS-CoV-2 entry and cell tropism. As such, researchers should be alert to the possibility of de novo S gene mutations emerging in tissue-culture propagated virus strains.     

Binding of the peptide deformylase on the ribosome surface modulates the exit tunnel interior

Proteosynthesis on ribosomes is regulated at many levels. Conformational changes of the ribosome, possibly induced by external factors, may transfer over large distances and contribute to the regulation. The molecular principles of this long-distance allostery within the ribosome remain poorly understood. Here, we use structural analysis and atomistic molecular dynamics simulations to investigate peptide deformylase (PDF), an enzyme that binds to the ribosome surface near the ribosomal protein uL22 during translation and chemically modifies the emerging nascent peptide. Our simulations of the entire ribosome-PDF complex reveal that the PDF undergoes a swaying motion on the ribosome surface at the submicrosecond timescale. We show that the PDF affects the conformational dynamics of parts of the ribosome over distances of more than 5 nm. Using a supervised-learning algorithm, we demonstrate that the exit tunnel is influenced by the presence or absence of PDF. Our findings suggest a possible effect of the PDF on the nascent peptide translocation through the ribosome exit tunnel.     

A gate in the selectivity filter of potassium channels

The selectivity filter of potassium channels is the structural element directly responsible for the selective and rapid conduction of K+, whereas other parts of the protein are thought to function as a molecular gate that either permits or blocks the passage of ions. However, whether the selectivity filter itself also possesses the ability to play the role of a gate is an unresolved question. Using free energy molecular dynamics simulations, it is shown that the reorientation of two peptide linkages in the selectivity filter of the KcsA K+ channel can lead to a stable nonconducting conformational state. Two microscopic factors influence the transition toward such a conformational state: the occupancy of one specific cation binding site in the selectivity filter (S2), and the strength of intersubunit interactions involving the GYG signature sequence. These results suggest that such conformational transitions occurring in the selectivity filter might be related to different K+ channel gating events, including C-type (slow) inactivation.