Menthol in mints elicits coolness sensation by selectively activating TRPM8 channel. Although structures of TRPM8 were determined in the apo and liganded states, the menthol-bounded state is unresolved. To understand how menthol activates the channel, we docked menthol to the channel and systematically validated our menthol binding models with thermodynamic mutant cycle analysis. We observed that menthol uses its hydroxyl group as a hand to specifically grab with R842, and its isopropyl group as legs to stand on I846 and L843. By imaging with fluorescent unnatural amino acid, we found that menthol binding induces wide-spread conformational rearrangements within the transmembrane domains. By Φ analysis based on single-channel recordings, we observed a temporal sequence of conformational changes in the S6 bundle crossing and the selectivity filter leading to channel activation. Therefore, our study suggested a ‘grab and stand’ mechanism of menthol binding and how menthol activates TRPM8 at the atomic level.

Blocking the association between the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike protein receptor-binding domain (RBD) and the human angiotensin-converting enzyme 2 (ACE2) is an attractive therapeutic approach to prevent the virus from entering human cells. While antibodies and other modalities have been developed to this end, d -amino acid peptides offer unique advantages, including serum stability, low immunogenicity, and low cost of production. Here, we designed potent novel D-peptide inhibitors that mimic the ACE2 α1-binding helix by searching a mirror-image version of the PDB. The two best designs bound the RBD with affinities of 29 and 31 nM and blocked the infection of Vero cells by SARS-CoV-2 with IC 50 values of 5.76 and 6.56 μM, respectively. Notably, both D-peptides neutralized with a similar potency the infection of two variants of concern: B.1.1.7 and B.1.351 in vitro . These potent D-peptide inhibitors are promising lead candidates for developing SARS-CoV-2 prophylactic or therapeutic treatments.

N-Methyl-D-aspartate (NMDA) receptors are glutamate-gated excitatory channels that play essential roles in brain functions. High-resolution structures have been solved for an allosterically inhibited and agonist-bound form of a functional NMDA receptor; however, other key functional states (particularly the active open-channel state) were only resolved at moderate resolutions by cryo-electron microscopy (cryo-EM). To decrypt the mechanism of the NMDA receptor activation, structural modeling is essential to provide presently missing information about structural dynamics. We performed systematic coarse-grained modeling using an elastic network model and related modeling/analysis tools (e.g., normal mode analysis, flexibility and hotspot analysis, cryo-EM flexible fitting, and transition pathway modeling) based on an active-state cryo-EM map. We observed extensive conformational changes that allosterically couple the extracellular regulatory and agonist-binding domains to the pore-forming trans-membrane domain (TMD), and validated these, to our knowledge, new observations against known mutational and functional studies. Our results predict two key modes of collective motions featuring shearing/twisting of the extracellular domains relative to the TMD, reveal subunit-specific flexibility profiles, and identify functional hotspot residues at key domain-domain interfaces. Finally, by examining the conformational transition pathway between the allosterically inhibited form and the active form, we predict a discrete sequence of domain motions, which propagate from the extracellular domains to the TMD. In summary, our results offer rich structural and dynamic information, which is consistent with the literature on structure-function relationships in NMDA receptors, and will guide in-depth studies on the activation dynamics of this important neurotransmitter receptor.

No abstract

As a prototype cellular sensor, the TRPV1 cation channel undergoes a closed-to-open gating transition in response to various physical and chemical stimuli including noxious heat. Despite recent progress, the molecular mechanism of heat activation of TRPV1 gating remains enigmatic. Toward decrypting the structural basis of TRPV1 heat activation, we performed extensive molecular dynamics simulations (with cumulative simulation time of ∼11 μs) for the wild-type channel and a constitutively active double mutant at different temperatures (30, 60, and 72°C), starting from a high-resolution closed-channel structure of TRPV1 solved by cryo-electron microscopy. In the wild-type simulations, we observed heat-activated conformational changes (e.g., expansion or contraction) in various key domains of TRPV1 (e.g., the S2-S3 and S4-S5 linkers) to prime the channel for gating. These conformational changes involve a number of dynamic hydrogen-bond interactions that were validated with previous mutational studies. Next, our mutant simulations observed channel opening after a series of conformational changes that propagate from the channel periphery to the channel pore via key intermediate domains (including the S2-S3 and S4-S5 linkers). The gating transition is accompanied by a large increase in the protein-water electrostatic interaction energy, which supports the contribution of desolvation of polar/charged residues to the temperature-sensitive TRPV1 gating. Taken together, our molecular dynamics simulations and analyses offered, to our knowledge, new structural, dynamic, and energetic information to guide future mutagenesis and functional studies of the TRPV1 channels and development of TRPV1-targeting drugs.

Protein-peptide interactions play a fundamental role in many cellular processes, but remain underexplored experimentally and difficult to model computationally. Here, we present PepNN-Struct and PepNN-Seq, structure and sequence-based approaches for the prediction of peptide binding sites on a protein. A main difficulty for the prediction of peptide-protein interactions is the flexibility of peptides and their tendency to undergo conformational changes upon binding. Motivated by this, we developed reciprocal attention to simultaneously update the encodings of peptide and protein residues while enforcing symmetry, allowing for information flow between the two inputs. PepNN integrates this module with modern graph neural network layers and a series of transfer learning steps are used during training to compensate for the scarcity of peptide-protein complex information. We show that PepNN-Struct achieves consistently high performance across different benchmark datasets. We also show that PepNN makes reasonable peptide-agnostic predictions, allowing for the identification of novel peptide binding proteins.

Transient receptor potential vanilloid 6 (TRPV6), a calcium-selective channel possessing six transmembrane domains (S1-S6) and intracellular N-/C-termini, plays crucial roles in calcium absorption in epithelia and bone and is involved in human diseases including vitamin-D deficiency, osteoporosis and cancer. The TRPV6 function and regulation remain poorly understood. Here by electrophysiology and Xenopus oocyte expression we found that Arg510 in the S4-S5 linker and Trp633 in the TRP helix are functionally critical. By coimmunoprecipitation and in vitro binding we found that the Arg510:Trp633 pair mediates binding between the linker (L) and C-terminal TRP helix (C) and that disruption of the linker/TRP helix (L/C) interaction substantially increases the channel activity. We also found that Trp361 in the N-terminal pre-S1 helix (N) and Ile637 mediate the pre-S1/TRP helix (N/C) binding which is functionally critical as well. Calcium imaging in human embryonic kidney cells further supported functional importance of Trp361 and Trp633. Besides, disruption of either interaction by blocking peptides activated TRPV6. The L/C interaction was required for the N/C interaction but not reversely. Together, we showed that the TRPV6 intramolecular S4-S5 linker to TRP helix and pre-S1 helix to TRP helix interactions, mediated by Arg510:Trp633 and Trp361:Ile637 bonding, respectively, are autoinhibitory and are required for maintaining TRPV6 at basal states. This study reveals a regulatory mechanism of TRPV6 activation-autoinhibition, which would help elucidating the corresponding mechanisms in other TRP channels. Supported by AIGSS (to RC), NSERC and KFoC (to XZC).

NMDA receptors are excitatory channels with critical functions in the physiology of central synapses. Their activation reaction proceeds as a series of kinetically distinguishable, reversible steps, whose structural bases are currently under investigation. Very likely, the earliest steps include glutamate binding to glycine-bound receptors and subsequent constriction of the ligand-binding domain. Later, three short linkers transduce this movement to open the gate by mechanical pulling on transmembrane helices. Here, we used molecular and kinetic simulations and double-mutant cycle analyses to show that a direct chemical interaction between GluN1-I642 (on M3 helix) and GluN2A-L550 (on L1-M1 linker) stabilizes receptors after they have opened and thus represents one of the structural changes that occur late in the activation reaction. This native interaction extends the current decay, and its absence causes deficits in charge transfer by GluN1-I642L, a pathogenic human variant.