Subunit Symmetry at the Extracellular Domain-Transmembrane Domain Interface in Acetylcholine Receptor Channel Gating

Transmitter molecules bind to synaptic acetylcholine receptor channels (AChRs) to promote a global channel-opening conformational change. Although the detailed mechanism that links ligand binding and channel gating is uncertain, the energy changes caused by mutations appear to be more symmetrical between subunits in the transmembrane domain compared with the extracellular domain. The only covalent connection between these domains is the pre-M1 linker, a stretch of five amino acids that joins strand β10 with the M1 helix. In each subunit, this linker has a central Arg (Arg^3′), which only in the non-α-subunits is flanked by positively charged residues. Previous studies showed that mutations of Arg^3′ in the α-subunit alter the gating equilibrium constant and reduce channel expression. We recorded single-channel currents and estimated the gating rate and equilibrium constants of adult mouse AChRs with mutations at the pre-M1 linker and the nearby residue Glu⁴⁵ in non-α-subunits. In all subunits, mutations of Arg^3′ had similar effects as in the α-subunit. In the ϵ-subunit, mutations of the flanking residues and Glu⁴⁵ had only small effects, and there was no energy coupling between ϵGlu⁴⁵ and ϵArg^3′. The non-α-subunit Arg^3′ residues had Φ-values that were similar to those for the α-subunit. The results suggest that there is a general symmetry between the AChR subunits during gating isomerization in this linker and that the central Arg is involved in expression more so than gating. The energy transfer through the AChR during gating appears to mainly involve Glu⁴⁵, but only in the α-subunits.     

Acetylcholine receptor gating: movement in the alpha-subunit extracellular domain

Acetylcholine receptor channel gating is a brownian conformational cascade in which nanometer-sized domains ("Phi blocks") move in staggering sequence to link an affinity change at the transmitter binding sites with a conductance change in the pore. In the alpha-subunit, the first Phi-block to move during channel opening is comprised of residues near the transmitter binding site and the second is comprised of residues near the base of the extracellular domain. We used the rate constants estimated from single-channel currents to infer the gating dynamics of Y127 and K145, in the inner and outer sheet of the beta-core of the alpha-subunit. Y127 is at the boundary between the first and second Phi blocks, at a subunit interface. alphaY127 mutations cause large changes in the gating equilibrium constant and with a characteristic Phi-value (Phi = 0.77) that places this residue in the second Phi-block. We also examined the effect on gating of mutations in neighboring residues deltaI43 (Phi = 0.86), epsilonN39 (complex kinetics), alphaI49 (no effect) and in residues that are homologous to alphaY127 on the epsilon, beta, and delta subunits (no effect). The extent to which alphaY127 gating motions are coupled to its neighbors was estimated by measuring the kinetic and equilibrium constants of constructs having mutations in alphaY127 (in both alpha subunits) plus residues alphaD97 or deltaI43. The magnitude of the coupling between alphaD97 and alphaY127 depended on the alphaY127 side chain and was small for both H (0.53 kcal/mol) and C (-0.37 kcal/mol) substitutions. The coupling across the single alpha-delta subunit boundary was larger (0.84 kcal/mol). The Phi-value for K145 (0.96) indicates that its gating motion is correlated temporally with the motions of residues in the first Phi-block and is not synchronous with those of alphaY127. This suggests that the inner and outer sheets of the alpha-subunit beta-core do not rotate as a rigid body.     

Loss-of-function mutations identified in the Helical domain of the G protein α-subunit,Gα2, of Dictyostelium discoideum

The guanine nucleotide binding regulatory proteins (G proteins) play essential roles in a wide variety of physiological processes, such as vision, hormone responses, olfaction, immune response, and development. The heterotrimeric G proteins consist of α-, β-, and γ-subunits and act as molecular switches to relay information from transmembrane receptors to intracellular effectors. The switch mechanism is a function of the inherent GTPase activity of the α-subunit. The α-subunit is comprised of two domains, the GTPase domain and the Helical domain. The GTPase domain performs all of the known α-subunit functions while little is know about the role of the Helical domain. To gain a better understanding of α-subunit function, we performed a screen for loss-of-function mutations, using the Gα2-subunit of Dictyostelium. Gα2 is essential for the developmental life cycle of Dictyostelium. It is known that the loss of Gα2 function results in a failure of cells to enter the developmental phase, producing a visibly abnormal phenotype. This allows the easy identification of amino acids essential to Gα2 function. A library of random point mutations in the gα2 cDNA was constructed using low fidelity polymerase chain reaction (PCR). The library was then expressed in a gα2 null cell line and screened for loss-of-function mutations. Mutations were identified in isolated clones by sequencing the gα2 insert. To date, sixteen single amino acids changes have been identified in Gα2 which result in loss-of-function. Of particular interest are seven mutations found in the Helical domain of the α-subunit. These loss-of-function mutations in the α-subunit Helical domain may provide important insight into its function.     

Energy and Structure of the M2 Helix in Acetylcholine Receptor-Channel Gating

We studied single-channel currents from neuromuscular acetylcholine receptor-channels with mutations in the pore-lining, M2 helix of the ɛ-subunit. Three parameters were quantified: 1), the diliganded gating equilibrium constant (E₂), which reflects the energy difference between C(losed) and O(pen) conformations; 2), the correlation between the opening rate constant and E₂ on a log-log scale (Φ), which illuminates the energy character of the residue (C- versus O-like) within the C↔O isomerization process; and 3), the open-channel current amplitude (i₀), which reports whether a mutation alters the energetics of ion permeation. The largest E₂ changes were observed in the cytoplasmic half of ɛM2 (5′, 9′, 12′, 13′, and 16′), with smaller changes apparent for residues ≥17′. Φ was ∼0.54 for most ɛM2 residues, but was ∼0.32 at the positions that had largest E₂ changes. An arginine substitution reduced i₀ significantly at six positions, with the magnitude of the reduction increasing, 16′→2′. The measurements suggest that the 9′, 12′, and 13′ residues experience large and late free-energy changes in the channel-opening process. We speculate that in the gating isomerization the pore-facing residues >6′ and <16′ experience multiple energy perturbations associated with changes in protein structure and, perhaps, hydration.     

Catch-and-Hold Activation of Muscle Acetylcholine Receptors Having Transmitter Binding Site Mutations

Agonists turn on receptors because their target sites have a higher affinity in the active versus resting conformation of the protein. We used single-channel electrophysiology to measure the lower-affinity (LA) and higher-affinity (HA) equilibrium dissociation constants for acetylcholine in adult-type muscle mouse nicotinic receptors (AChRs) having mutations of agonist binding site amino acids. For a series of agonists and for all mutations of αY93, αG147, αW149, αY190, αY198, εW55, and δW57, the change in LA binding energy was approximately half that in HA binding energy. The results were analyzed as a linear free energy relationship between LA and HA agonist binding, the slope of which (κ) gives the fraction of the overall binding chemical potential where the LA complex is established. The linear correlation between LA and HA binding energies suggests that the overall binding process is by an integrated mechanism (catch-and-hold). For the agonist and the above mutations, κ ∼ 0.5, but side-chain substitutions of two residues had a slope that was significantly higher (0.90; αG153) or lower (0.25; εP121). The results suggest that backbone rearrangements in loop B, loop C, and the non-α surface participate in both LA binding and the LA ↔ HA affinity switch. It appears that all of the intermediate steps in AChR activation comprise a single, energetically coupled process.     

Second Transmembrane Helix (M2) and Long Range Coupling in Ca2+-ATPase

The actuator (A) domain of sarco(endo)plasmic reticulum Ca²⁺-ATPase not only plays a catalytic role but also undergoes large rotational movements that influence the distant transport sites through connections with transmembrane helices M1 and M2. Here we explore the importance of long helix M2 and its junction with the A domain by disrupting the helix structure and elongating with insertions of five glycine residues. Insertions into the membrane region of M2 and the top junctional segment impair Ca²⁺ transport despite reasonable ATPase activity, indicating that they are uncoupled. These mutants fail to occlude Ca²⁺. Those at the top segment also exhibited accelerated phosphoenzyme isomerization E1P → E2P. Insertions into the middle of M2 markedly accelerate E2P hydrolysis and cause strong resistance to inhibition by luminal Ca²⁺. Insertions along almost the entire M2 region inhibit the dephosphorylated enzyme transition E2 → E1. The results pinpoint which parts of M2 control cytoplasm gating and which are critical for luminal gating at each stage in the transport cycle and suggest that proper gate function requires appropriate interactions, tension, and/or rigidity in the M2 region at appropriate times for coupling with A domain movements and catalysis.     

Role of conserved Fα-helix residues in the native fold and stability of the kinase-inhibited oxy state of the oxygen-sensing FixL protein from Sinorhizobium meliloti

The oxygen-sensing FixL protein from Sinorhizobium meliloti is part of the heme-PAS family of gas sensors that regulate many important signal transduction pathways in a wide variety of organisms. We examined the role of the conserved F_α-9 arginine 200 and several other conserved residues on the proximal F_α-helix in the heme domain of SmFixL* using site-directed mutagenesis in conjunction with UV-visible, EPR, and resonance Raman spectroscopy. The F_α-helix variants R200A, E, Q, H, Y197A, and D195A were expressed at reasonable levels and purified to homogeneity. The R200I and Y201A variants did not express in observable quantities. Tyrosine 201 is crucial for forming the native protein fold of SmFixL* while Y197 and R200 are important for stabilizing the kinase-inhibited oxy state. Our results show a clear correlation between H-bond donor ability of the F_α-9 side chain and the rate of heme autoxidation. This trend in conjunction with crystal structures of liganded BjFixL heme domains, show that H-bonding between the conserved F_α-9 arginine and the heme-6-propionate group contributes to the kinetic stability of the kinase-inactivated, oxy state of SmFixL*.     

A residue at the cytoplasmic entrance of BK-type channels regulating single-channel opening by its hydrophobicity

Single large-conductance calcium-activated K⁺ (BK) channels encoded by the mSlo gene usually have synchronous gating, but a Drosophila dSlo (A2/C2/E2/G5/10) splice variant (dSlo1A) exhibits very flickery openings. To probe this difference in gating, we constructed a mutant I323T. This channel exhibits four subconductance levels similar to those of dSlo1A. Rectification of the single-channel current-voltage relation of I323T decreased as [Ca²⁺ ]_in increased from 10 to 300 μM. Mutagenesis suggests that the hydrophobicity of the residue at the position is important for the wild-type gating; i.e., increasing hydrophobicity prolongs open duration. Molecular dynamics simulation suggests that four hydrophobic pore-lining residues at position 323 of mSlo act cooperatively in a “shutter-like” mechanism gating the permeation of K⁺ ions. Rate-equilibrium free energy relations analysis shows that the four I323 residues in an mSlo channel have a conformation 65% similar to the closed conformation during gating. Based on these observations, we suggest that the appearance of rectification and substates of BK-type channels arise from a reduction of the cooperativity among these four residues and a lower probability of being open.     

Modeling Neuronal Nicotinic and GABA Receptors: Important Interface Salt-Links and Protein Dynamics

Protein motions in the Cys-loop ligand-gated ion receptors that govern the gating mechanism are still not well understood. The details as to how motions in the ligand-binding domain are translated to the transmembrane domain and how subunit rotations are linked to bring about the cooperative movements involved in gating are under investigation. Homology models of the α4β2 nicotinic acetylcholine (nACh) and β2α1γ2 GABA receptors were constructed based on the torpedo neuromuscular-like nicotinic receptor structure. The template constructed for the full electron microscopy structure must be considered more reliable for structure-function studies due to the preservation of the E45-R209 salt-link. Many other salt-links are seen to transiently form, including switching off of the E45-R209 link, within a network of potential salt-links at the binding domain to the transmembrane domain interface region. Several potentially important intersubunit salt-links form in both the nAChR and GABAR structures during the simulation and appear conserved across many subunit combinations, such as the salt-link between α4.E262 and β2.K255 in nAChR (β2.E262 and α1.K263 in GABAR), at the top of the pore-lining M2 helices, and the intersubunit link of R210 on the M1-linker to E168 on the β8-sheet of the adjacent subunit in the GABA receptor (E175-K46 being the structurally equivalent link in the nAChR, with reversed polarity). A network of other salt-links may be vital for transmitting the cooperative gating motions between subunits that become biased upon ligand binding. The changes seen in the simulations suggest that this network of salt-links helps to set limits and specific states for the conformational changes involved in gating of the receptor. We hope that these hypotheses will be tested experimentally in the near future.     

Mutations in a Conserved Domain of E. coli MscS to the Most Conserved Superfamily Residue Leads to Kinetic Changes

In Escherichia coli (E. coli) the mechanosensitive channel of small conductance, MscS, gates in response to membrane tension created from acute external hypoosmotic shock, thus rescuing the bacterium from cell lysis. E. coli MscS is the most well studied member of the MscS superfamily of channels, whose members are found throughout the bacterial and plant kingdoms. Homology to the pore lining helix and upper vestibule domain of E. coli MscS is required for inclusion into the superfamily. Although highly conserved, in the second half of the pore lining helix (TM3B), E. coli MscS has five residues significantly different from other members of the superfamily. In superfamilies such as this, it remains unclear why variations within such a homologous region occur: is it tolerance of alternate residues, or does it define functional variance within the superfamily? Point mutations (S114I/T, L118F, A120S, L123F, F127E/K/T) and patch clamp electrophysiology were used to study the effect of changing these residues in E. coli MscS on sensitivity and gating. The data indicate that variation at these locations do not consistently lead to wildtype channel phenotypes, nor do they define large changes in mechanosensation, but often appear to effect changes in the E. coli MscS channel gating kinetics.     

Long-Range Coupling in an Allosteric Receptor Revealed by Mutant Cycle Analysis

The functional coupling of residues that are far apart in space is the quintessential property of allosteric proteins. For example, in Cys-loop receptors, the gating of an intrinsic ion channel is allosterically regulated by the binding of small molecule neurotransmitters 50-60 Å from the channel gate. Some residues near the binding site must have as their primary function the communication of the binding event to the gating region. These gating pathway residues are essential to function, but their identification and characterization can be challenging. This work introduces a simple strategy, derived from mutant cycle analysis, for identifying gating pathway residues using macroscopic measurements alone. In the exemplar Cys-loop receptor, the nicotinic acetylcholine receptor, a well-characterized reporter mutation (βL9′S) known to impact gating, was combined with mutations of target residues in the ligand-binding domain hypothesized or previously found to be functionally significant. A mutant cycle analysis of the macroscopic EC₅₀ measurements can then provide insights into the role of the target residue. This new method, elucidating long-range functional coupling in allosteric receptors, can be applied to several reporter mutations in a wide variety of receptors to identify previously characterized and novel mutations that impact the gating pathway. We support our interpretation of macroscopic data with single-channel studies. Elucidating long-range functional coupling in allosteric receptors should be broadly applicable to determining functional roles of residues in allosteric receptors.     

Extracellular intersubunit interactions modulate epithelial Na+ channel gating

Epithelial Na⁺ channels (ENaCs) and related channels have large extracellular domains where specific factors interact and induce conformational changes, leading to altered channel activity. However, extracellular structural transitions associated with changes in ENaC activity are not well defined. Using crosslinking and two-electrode voltage clamp in Xenopus oocytes, we identified several pairs of functional intersubunit contacts where mouse ENaC activity was modulated by inducing or breaking a disulfide bond between introduced Cys residues. Specifically, crosslinking E499C in the β-subunit palm domain and N510C in the α-subunit palm domain activated ENaC, whereas crosslinking βE499C with αQ441C in the α-subunit thumb domain inhibited ENaC. We determined that bridging βE499C to αN510C or αQ441C altered the Na⁺ self-inhibition response via distinct mechanisms. Similar to bridging βE499C and αQ441C, we found that crosslinking palm domain αE557C with thumb domain γQ398C strongly inhibited ENaC activity. In conclusion, we propose that certain residues at specific subunit interfaces form microswitches that convey a conformational wave during ENaC gating and its regulation.     

Mapping heat exchange in an allosteric protein

Nicotinic acetylcholine receptors (AChRs) are synaptic ion channels that spontaneously isomerize (i.e., gate) between resting and active conformations. We used single-molecule electrophysiology to measure the temperature dependencies of mouse neuromuscular AChR gating rate and equilibrium constants. From these we estimated free energy, enthalpy, and entropy changes caused by mutations of amino acids located between the transmitter binding sites and the middle of the membrane domain. The range of equilibrium enthalpy change (13.4 kcal/mol) was larger than for free energy change (5.5 kcal/mol at 25°C). For two residues, the slope of the rate-equilibrium free energy relationship (Φ) was approximately constant with temperature. Mutant cycle analysis showed that both free energies and enthalpies are additive for energetically independent mutations. We hypothesize that changes in energy associated with changes in structure mainly occur close to the site of the mutation, and, hence, that it is possible to make a residue-by-residue map of heat exchange in the AChR gating isomerization. The structural correlates of enthalpy changes are discussed for 12 different mutations in the protein.     

Assignments of human integrin α1I domain in the apo and Mg2+ bound states

The α1β1 integrin receptor binds to its main extracellular ligand, collagen, through an inserted domain in its α-subunit called the αI domain (αI). αI contains a metal binding site that allows collagen to coordinate to the domain through a divalent metal ion. Here we report the backbone assignments of the apo and Mg²⁺ bound state of the isolated human α1I and the chemical shift changes resulting from metal coordination.     

Role of the DELSEED Loop in Torque Transmission of F1-ATPase

F₁-ATPase is an ATP-driven rotary motor that generates torque at the interface between the catalytic β-subunits and the rotor γ-subunit. The β-subunit inwardly rotates the C-terminal domain upon nucleotide binding/dissociation; hence, the region of the C-terminal domain that is in direct contact with γ—termed the DELSEED loop—is thought to play a critical role in torque transmission. We substituted all the DELSEED loop residues with alanine to diminish specific DELSEED loop-γ interactions and with glycine to disrupt the loop structure. All the mutants rotated unidirectionally with kinetic parameters comparable to those of the wild-type F₁, suggesting that the specific interactions between DELSEED loop and γ is not involved in cooperative interplays between the catalytic β-subunits. Glycine substitution mutants generated half the torque of the wild-type F₁, whereas the alanine mutant generated comparable torque. Fluctuation analyses of the glycine/alanine mutants revealed that the γ-subunit was less tightly held in the α₃β₃-stator ring of the glycine mutant than in the wild-type F₁ and the alanine mutant. Molecular dynamics simulation showed that the DELSEED loop was disordered by the glycine substitution, whereas it formed an α-helix in the alanine mutant. Our results emphasize the importance of loop rigidity for efficient torque transmissions.     

On the Activation of Integrin αIIbβ3: Outside-in and Inside-out Pathways

Integrin αIIbβ3 is a member of the integrin family of transmembrane proteins present on the plasma membrane of platelets. Integrin αIIbβ3 is widely known to regulate the process of thrombosis via activation at its cytoplasmic side by talin and interaction with the soluble fibrinogen. It is also reported that three groups of interactions restrain integrin family members in the inactive state, including a set of salt bridges on the cytoplasmic side of the transmembrane domain of the integrin α- and β-subunits known as the inner membrane clasp, hydrophobic packing of a few transmembrane residues on the extracellular side between the α- and β-subunits that is known as the outer membrane clasp, and the key interaction group of the βA domain (located on the β-subunit head domain) with the βTD (proximal to the plasma membrane on the β-subunit). However, molecular details of this key interaction group as well as events that lead to detachment of the βTD and βA domains have remained ambiguous. In this study, we use molecular dynamics models to take a comprehensive outside-in and inside-out approach at exploring how integrin αIIbβ3 is activated. First, we show that talin’s interaction with the membrane-proximal and membrane-distal regions of integrin cytoplasmic-transmembrane domains significantly loosens the inner membrane clasp. Talin also interacts with an additional salt bridge (R734-E1006), which facilitates integrin activation through the separation of the integrin’s α- and β-subunits. The second part of our study classifies three types of interactions between RGD peptides and the extracellular domains of integrin αIIbβ3. Finally, we show that the interaction of the Arg of the RGD sequence may activate integrin via disrupting the key interaction group between K350 on the βA domain and S673/S674 on the βTD.     

Multiple Roles of the Extracellular Vestibule Amino Acid Residues in the Function of the Rat P2X4 Receptor

The binding of ATP to trimeric P2X receptors (P2XR) causes an enlargement of the receptor extracellular vestibule, leading to opening of the cation-selective transmembrane pore, but specific roles of vestibule amino acid residues in receptor activation have not been evaluated systematically. In this study, alanine or cysteine scanning mutagenesis of V47–V61 and F324–N338 sequences of rat P2X4R revealed that V49, Y54, Q55, F324, and G325 mutants were poorly responsive to ATP and trafficking was only affected by the V49 mutation. The Y54F and Y54W mutations, but not the Y54L mutation, rescued receptor function, suggesting that an aromatic residue is important at this position. Furthermore, the Y54A and Y54C receptor function was partially rescued by ivermectin, a positive allosteric modulator of P2X4R, suggesting a rightward shift in the potency of ATP to activate P2X4R. The Q55T, Q55N, Q55E, and Q55K mutations resulted in non-responsive receptors and only the Q55E mutant was ivermectin-sensitive. The F324L, F324Y, and F324W mutations also rescued receptor function partially or completely, ivermectin action on channel gating was preserved in all mutants, and changes in ATP responsiveness correlated with the hydrophobicity and side chain volume of the substituent. The G325P mutant had a normal response to ATP, suggesting that G325 is a flexible hinge. A topological analysis revealed that the G325 and F324 residues disrupt a β-sheet upon ATP binding. These results indicate multiple roles of the extracellular vestibule amino acid residues in the P2X4R function: the V49 residue is important for receptor trafficking to plasma membrane, the Y54 and Q55 residues play a critical role in channel gating and the F324 and G325 residues are critical for vestibule widening.     

Conformation of the transmembrane domain of the c-erbB-2 oncogene-encoded protein in its monomeric and dimeric states

The human c-erbB-2 oncogene is homologous to the ratneu oncogene, both encoding transmembrane growth factor receptors. Overexpression and point mutations in the transmembrane domain of the encoded proteins in both cases have been implicated in cell transformation and carcinogenesis. In the case of theneu protein, it has been proposed that these effects are mediated by conformational preferences for anα-helix in the transmembrane domain, which facilitates receptor dimerization, an important step in the signal transduction process. To examine whether this is the case for c-erbB-2 as well, we have used conformational energy analysis to determine the preferred three-dimensional structures for the transmembrane domain of the c-erbB-2 protein from residues 650 to 668 with Val (nontransforming) and Glu (transforming) at position 659. The global minimum energy conformation for the Val-659 peptide from the normal, nontransforming protein was found to contain several bends, whereas the global minimum energy conformation for Glu-659 peptide from the mutant, transforming protein was found to beα-helical. Thus, the difference in conformational preferences for these transmembrane domains may explain the difference in transforming ability of these proteins. The presence of higher-energyα-helical conformations for the transmembrane domain from the normal Val-659 protein may provide an explanation for the presence of a transforming effect from overexpression of c-erbB-2. In addition, docking of the oncogenic sequences in theirα-helical and bend conformations shows that the all-α-helical dimer is clearly favored energetically over the bend dimer.     

The site for the allosteric activator GTP of Escherichia coli UMP kinase

UMP kinase (UMPK), a key bacterial pyrimidine nucleotide biosynthesis enzyme, is UTP-inhibited and GTP-activated. We delineate the GTP site of Escherichia coli UMPK by alanine mutagenesis of R92, H96, R103, W119 or R130, abolishing GTP activation; of S124 and R127, decreasing affinity for GTP; and of N111 and D115, with little detrimental effect. We exclude the correspondence with the modulatory ATP site of Bacillus anthracis UMPK, confirming the functionality of the GTP site found by Evrin. Mutants R92A, H96A and R127A are constitutively activated, suggesting key roles of these residues in allosteric signal transduction and of positive charge neutralization in triggering activation. No mutation hampered UTP inhibition, excluding overlapping of the UTP and GTP sites.     

A potential interaction between the SARS-CoV-2 spike protein and nicotinic acetylcholine receptors

Changeux et al. (Changeux et al. C. R. Biol. 343:33–39.) recently suggested that the SARS-CoV-2 spike protein may interact with nicotinic acetylcholine receptors (nAChRs) and that such interactions may be involved in pathology and infectivity. This hypothesis is based on the fact that the SARS-CoV-2 spike protein contains a sequence motif similar to known nAChR antagonists. Here, we use molecular simulations of validated atomically detailed structures of nAChRs and of the spike to investigate the possible binding of the Y674-R685 region of the spike to nAChRs. We examine the binding of the Y674-R685 loop to three nAChRs, namely the human α4β2 and α7 subtypes and the muscle-like αβγδ receptor from Tetronarce californica. Our results predict that Y674-R685 has affinity for nAChRs. The region of the spike responsible for binding contains a PRRA motif, a four-residue insertion not found in other SARS-like coronaviruses. The conformational behavior of the bound Y674-R685 is highly dependent on the receptor subtype; it adopts extended conformations in the α4β2 and α7 complexes but is more compact when bound to the muscle-like receptor. In the α4β2 and αβγδ complexes, the interaction of Y674-R685 with the receptors forces the loop C region to adopt an open conformation, similar to other known nAChR antagonists. In contrast, in the α7 complex, Y674-R685 penetrates deeply into the binding pocket in which it forms interactions with the residues lining the aromatic box, namely with TrpB, TyrC1, and TyrC2. Estimates of binding energy suggest that Y674-R685 forms stable complexes with all three nAChR subtypes. Analyses of simulations of the glycosylated spike show that the Y674-R685 region is accessible for binding. We suggest a potential binding orientation of the spike protein with nAChRs, in which they are in a nonparallel arrangement to one another.