Deorphanizing Peptides Using Structure Prediction.

Many endogenous peptides rely on signaling pathways to exert their function, but identifying their cognate receptors remains a challenging problem. We investigate the use of AlphaFold-Multimer complex structure prediction together with transmembrane topology prediction for peptide deorphanization. We find that AlphaFold's confidence metrics have strong performance for prioritizing true peptide-receptor interactions. In a library of 1112 human receptors, the method ranks true receptors in the top percentile on average for 11 benchmark peptide-receptor pairs.

An open state of a voltage-gated sodium channel involving a π-helix and conserved pore-facing asparagine.

Voltage-gated sodium (Nav) channels play critical roles in propagating action potentials and otherwise manipulating ionic gradients in excitable cells. These channels open in response to membrane depolarization, selectively permeating sodium ions until rapidly inactivating. Structural characterization of the gating cycle in this channel family has proved challenging, particularly due to the transient nature of the open state. A structure from the bacterium Magnetococcus marinus Nav (NavMs) was initially proposed to be open, based on its pore diameter and voltage-sensor conformation. However, the functional annotation of this model, and the structural details of the open state, remain disputed. In this work, we used molecular modeling and simulations to test possible open-state models of NavMs. The full-length experimental structure, termed here the α-model, was consistently dehydrated at the activation gate, indicating an inability to conduct ions. Based on a spontaneous transition observed in extended simulations, and sequence/structure comparison to other Nav channels, we built an alternative π-model featuring a helix transition and the rotation of a conserved asparagine residue into the activation gate. Pore hydration, ion permeation, and state-dependent drug binding in this model were consistent with an open functional state. This work thus offers both a functional annotation of the full-length NavMs structure and a detailed model for a stable Nav open state, with potential conservation in diverse ion-channel families.

Dual regulation of TRPV1 channels by phosphatidylinositol via functionally distinct binding sites.

Regulation of the heat- and capsaicin-activated transient receptor potential vanilloid 1 (TRPV1) channel by phosphoinositides is complex and controversial. In the most recent TRPV1 cryo-EM structure, endogenous phosphatidylinositol (PtdIns) was detected in the vanilloid binding site, and phosphoinositides were proposed to act as competitive vanilloid antagonists. This model is difficult to reconcile with phosphatidylinositol 4,5-bisphosphate [PtdIns(4,5)P2] being a well-established positive regulator of TRPV1. Here we show that in the presence of PtdIns(4,5)P2 in excised patches, PtdIns, but not PtdIns(4)P, partially inhibited TRPV1 activity at low, but not at high capsaicin concentrations. This is consistent with PtdIns acting as a competitive vanilloid antagonist. However, in the absence of PtdIns(4,5)P2, PtdIns partially stimulated TRPV1 activity. We computationally identified residues, which are in contact with PtdIns, but not with capsaicin in the vanilloid binding site. The I703A mutant of TRPV1 showed increased sensitivity to capsaicin, as expected when removing the effect of an endogenous competitive antagonist. I703A was not inhibited by PtdIns in the presence of PtdIns(4,5)P2, but it was still activated by PtdIns in the absence of PtdIns(4,5)P2 indicating that inhibition, but not activation by PtdIns proceeds via the vanilloid binding site. In molecular dynamics simulations, PtdIns was more stable than PtdIns(4,5)P2 in this inhibitory site, whereas PtdIns(4,5)P2 was more stable than PtdIns in a previously identified, nonoverlapping, putative activating binding site. Our data indicate that phosphoinositides regulate channel activity via functionally distinct binding sites, which may explain some of the complexities of the effects of these lipids on TRPV1.

Pulsed Electric Fields Can Create Pores in the Voltage Sensors of Voltage-Gated Ion Channels.

Pulsed electric fields are increasingly used in medicine to transiently increase the cell membrane permeability via electroporation to deliver therapeutic molecules into the cell. One type of event that contributes to this increase in membrane permeability is the formation of pores in the membrane lipid bilayer. However, electrophysiological measurements suggest that membrane proteins are affected as well, particularly voltage-gated ion channels (VGICs). The molecular mechanisms by which the electric field could affects these molecules remain unidentified. In this study, we used molecular dynamics simulations to unravel the molecular events that take place in different VGICs when exposing them to electric fields mimicking electroporation conditions. We show that electric fields can induce pores in the voltage-sensor domains (VSDs) of different VGICs and that these pores form more easily in some channels than in others. We demonstrate that poration is more likely in VSDs that are more hydrated and are electrostatically more favorable for the entry of ions. We further show that pores in VSDs can expand into so-called complex pores, which become stabilized by lipid headgroups. Our results suggest that such complex pores are considerably more stable than conventional lipid pores, and their formation can lead to severe unfolding of VSDs from the channel. We anticipate that such VSDs become dysfunctional and unable to respond to changes in transmembrane voltage, which is in agreement with previous electrophysiological measurements showing a decrease in the voltage-dependent transmembrane ionic currents after pulse treatment. Finally, we discuss the possibility of activation of VGICs by submicrosecond-duration pulses. Overall, our study reveals a new, to our knowledge, mechanism of electroporation through membranes containing VGICs.

Molecular Insights from Conformational Ensembles via Machine Learning.

Biomolecular simulations are intrinsically high dimensional and generate noisy data sets of ever-increasing size. Extracting important features from the data is crucial for understanding the biophysical properties of molecular processes, but remains a big challenge. Machine learning (ML) provides powerful dimensionality reduction tools. However, such methods are often criticized as resembling black boxes with limited human-interpretable insight. We use methods from supervised and unsupervised ML to efficiently create interpretable maps of important features from molecular simulations. We benchmark the performance of several methods, including neural networks, random forests, and principal component analysis, using a toy model with properties reminiscent of macromolecular behavior. We then analyze three diverse biological processes: conformational changes within the soluble protein calmodulin, ligand binding to a G protein-coupled receptor, and activation of an ion channel voltage-sensor domain, unraveling features critical for signal transduction, ligand binding, and voltage sensing. This work demonstrates the usefulness of ML in understanding biomolecular states and demystifying complex simulations.

Evolutionarily Conserved Interactions within the Pore Domain of Acid-Sensing Ion Channels.

Despite the sequence homology between acid-sensing ion channels (ASICs) and epithelial sodium channel (ENaCs), these channel families display very different functional characteristics. Whereas ASICs are gated by protons and show a relatively low degree of selectivity for sodium over potassium, ENaCs are constitutively active and display a remarkably high degree of sodium selectivity. To decipher if some of the functional diversity originates from differences within the transmembrane helices (M1 and M2) of both channel families, we turned to a combination of computational and functional interrogations, using statistical coupling analysis and mutational studies on mouse ASIC1a. The coupling analysis suggests that the relative position of M1 and M2 in the upper part of the pore domain is likely to remain constant during the ASIC gating cycle, whereas they may undergo relative movements in the lower part. Interestingly, our data suggest that to account for coupled residue pairs being in close structural proximity, both domain-swapped and nondomain-swapped ASIC M2 conformations need to be considered. Such conformational flexibility is consistent with structural work, which suggested that the lower part of M2 can adopt both domain-swapped and nondomain-swapped conformations. Overall, mutations to residues in the middle and lower pore were more likely to affect gating and/or ion selectivity than those in the upper pore. Indeed, disrupting the putative interaction between a highly conserved Trp/Glu residue pair in the lower pore is detrimental to gating and selectivity, although this interaction might occur in both domain-swapped and nonswapped conformations. Finally, our results suggest that the greater number of larger, aromatic side chains in the ENaC M2 helix may contribute to the constitutive activity of these channels at a resting pH. Together, the data highlight differences in the transmembrane domains of these closely related ion channels that may help explain some of their distinct functional properties.

A Structural Basis for 129 Xe Hyper-CEST Signal in TEM-1 ß-Lactamase.

Genetically encoded (GE) contrast agents detectable by magnetic resonance imaging (MRI) enable non-invasive visualization of gene expression and cell proliferation at virtually unlimited penetration depths. Using hyperpolarized 129 Xe in combination with chemical exchange saturation transfer, an MR contrast approach known as hyper-CEST, enables ultrasensitive protein detection and biomolecular imaging. GE MRI contrast agents developed to date include nanoscale proteinaceous gas vesicles as well as the monomeric bacterial proteins TEM-1 ß-lactamase (bla) and maltose binding protein (MBP). To improve understanding of hyper-CEST NMR with proteins, structural and computational studies were performed to further characterize the Xe-bla interaction. X-ray crystallography validated the location of a high-occupancy Xe binding site predicted by MD simulations, and mutagenesis experiments confirmed this Xe site as the origin of the observed CEST contrast. Structural studies and MD simulations with representative bla mutants offered additional insight regarding the relationship between local protein structure and CEST contrast.

Exploring the Viral Channel KcvPBCV-1 Function via Computation.

Viral potassium channels (Kcv) are homologous to the pore module of complex [Formula: see text]-selective ion channels of cellular organisms. Due to their relative simplicity, they have attracted interest towards understanding the principles of [Formula: see text] conduction and channel gating. In this work, we construct a homology model of the [Formula: see text] open state, which we validate by studying the binding of known blockers and by monitoring ion conduction through the channel. Molecular dynamics simulations of this model reveal that the re-orientation of selectivity filter carbonyl groups coincides with the transport of potassium ions, suggesting a possible mechanism for fast gating. In addition, we show that the voltage sensitivity of this mechanism can originate from the relocation of potassium ions inside the selectivity filter. We also explore the interaction of [Formula: see text] with the surrounding bilayer and observe the binding of lipids in the area between two adjacent subunits. The model is available to the scientific community to further explore the structure/function relationship of Kcv channels.

Free-energy landscape of ion-channel voltage-sensor-domain activation.

Voltage sensor domains (VSDs) are membrane-bound protein modules that confer voltage sensitivity to membrane proteins. VSDs sense changes in the transmembrane voltage and convert the electrical signal into a conformational change called activation. Activation involves a reorganization of the membrane protein charges that is detected experimentally as transient currents. These so-called gating currents have been investigated extensively within the theoretical framework of so-called discrete-state Markov models (DMMs), whereby activation is conceptualized as a series of transitions across a discrete set of states. Historically, the interpretation of DMM transition rates in terms of transition state theory has been instrumental in shaping our view of the activation process, whose free-energy profile is currently envisioned as composed of a few local minima separated by steep barriers. Here we use atomistic level modeling and well-tempered metadynamics to calculate the configurational free energy along a single transition from first principles. We show that this transition is intrinsically multidimensional and described by a rough free-energy landscape. Remarkably, a coarse-grained description of the system, based on the use of the gating charge as reaction coordinate, reveals a smooth profile with a single barrier, consistent with phenomenological models. Our results bridge the gap between microscopic and macroscopic descriptions of activation dynamics and show that choosing the gating charge as reaction coordinate masks the topological complexity of the network of microstates participating in the transition. Importantly, full characterization of the latter is a prerequisite to rationalize modulation of this process by lipids, toxins, drugs, and genetic mutations.

Voltage-gated ion channel modulation by lipids: insights from molecular dynamics simulations.

Cells commonly use lipids to modulate the function of ion channels. The lipid content influences the amplitude of the ionic current and changes the probability of voltage-gated ion channels being in the active or in the resting states. Experimental findings inferred from a variety of techniques and molecular dynamics studies have revealed a direct interaction between the lipid headgroups and the ion channel residues, suggesting an influence on the ion channel function. On the other hand the alteration of the lipids may in principle modify the overall electrostatic environment of the channel, and hence the transmembrane potential, leading to an indirect modulation, i.e. a global effect. Here we have investigated the structural and dynamical properties of the voltage-gated potassium channel Kv1.2 embedded in bilayers with modified upper or lower leaflet compositions corresponding to realistic biological scenarios: the first relates to the effects of sphingomyelinase, an enzyme that modifies the composition of lipids of the outer membrane leaflets, and the second to the effect of the presence of a small fraction of PIP2, a highly negatively charged lipid known to modulate voltage-gated channel function. Our molecular dynamics simulations do not enable to exclude the global effect mechanism in the former case. For the latter, however, it is shown that local interactions between the ion channel and the lipid headgroups are key-elements of the modulation.

Membrane Protein Structure, Function, and Dynamics: a Perspective from Experiments and Theory.

Membrane proteins mediate processes that are fundamental for the flourishing of biological cells. Membrane-embedded transporters move ions and larger solutes across membranes; receptors mediate communication between the cell and its environment and membrane-embedded enzymes catalyze chemical reactions. Understanding these mechanisms of action requires knowledge of how the proteins couple to their fluid, hydrated lipid membrane environment. We present here current studies in computational and experimental membrane protein biophysics, and show how they address outstanding challenges in understanding the complex environmental effects on the structure, function, and dynamics of membrane proteins.

Dual effect of phosphatidylinositol (4,5)-bisphosphate PIP(2) on Shaker K(+) [corrected] channels.

Phosphatidylinositol (4,5)-bisphosphate (PIP(2)) is a phospholipid of the plasma membrane that has been shown to be a key regulator of several ion channels. Functional studies and more recently structural studies of Kir channels have revealed the major impact of PIP(2) on the open state stabilization. A similar effect of PIP(2) on the delayed rectifiers Kv7.1 and Kv11.1, two voltage-gated K(+) channels, has been suggested, but the molecular mechanism remains elusive and nothing is known on PIP(2) effect on other Kv such as those of the Shaker family. By combining giant-patch ionic and gating current recordings in COS-7 cells, and voltage-clamp fluorimetry in Xenopus oocytes, both heterologously expressing the voltage-dependent Shaker channel, we show that PIP(2) exerts 1) a gain-of-function effect on the maximal current amplitude, consistent with a stabilization of the open state and 2) a loss-of-function effect by positive-shifting the activation voltage dependence, most likely through a direct effect on the voltage sensor movement, as illustrated by molecular dynamics simulations.

Modulation of the IKS channel by PIP2 requires two binding sites per monomer.

The phosphatidyl-inositol-4,5-bisphosphate (PIP2) lipid has been shown to be crucial for the coupling between the voltage sensor and the pore of the potassium voltage-gated KV7 channel family, especially the KV7.1 channel. Expressed in the myocardium membrane, KV7.1 forms a complex with KCNE1 auxiliary subunits to generate the IKS current. Here we present molecular models of the transmembrane region of this complex in its three known states, namely the Resting/Closed (RC), the Intermediate/Closed (IC), and the Activated/Open (AO), robustness of which is assessed by agreement with a range of biophysical data. Molecular Dynamics (MD) simulations of these models embedded in a lipid bilayer including phosphatidyl-inositol-4,5-bisphosphate (PIP2) lipids show that in presence of KCNE1, two PIP2 lipids are necessary to stabilize each state. The simulations also show that KCNE1 interacts with both PIP2 binding sites, forming a tourniquet around the pore and preventing its opening. The present investigation provides therefore key molecular elements that govern the role of PIP2 in KCNE1 modulation of IKS channels, possibly a common mechanism by which auxiliary KCNE subunits might modulate a variety of other ion channels.

An integrated platform approach enables discovery of potent, selective and ligand-competitive cyclic peptides targeting the GIP receptor.

In any drug discovery effort, the identification of hits for further optimisation is of crucial importance. For peptide therapeutics, display technologies such as mRNA display have emerged as powerful methodologies to identify these desired de novo hit ligands against targets of interest. The diverse peptide libraries are genetically encoded in these technologies, allowing for next-generation sequencing to be used to efficiently identify the binding ligands. Despite the vast datasets that can be generated, current downstream methodologies, however, are limited by low throughput validation processes, including hit prioritisation, peptide synthesis, biochemical and biophysical assays. In this work we report a highly efficient strategy that combines bioinformatic analysis with state-of-the-art high throughput peptide synthesis to identify nanomolar cyclic peptide (CP) ligands of the human glucose-dependent insulinotropic peptide receptor (hGIP-R). Furthermore, our workflow is able to discriminate between functional and remote binding non-functional ligands. Efficient structure-activity relationship analysis (SAR) combined with advanced in silico structural studies allow deduction of a thorough and holistic binding model which informs further chemical optimisation, including efficient half-life extension. We report the identification and design of the first de novo, GIP-competitive, incretin receptor family-selective CPs, which exhibit an in vivo half-life up to 10.7 h in rats. The workflow should be generally applicable to any selection target, improving and accelerating hit identification, validation, characterisation, and prioritisation for therapeutic development.

Two-stage electro-mechanical coupling of a KV channel in voltage-dependent activation.

In voltage-gated potassium (KV) channels, the voltage-sensing domain (VSD) undergoes sequential activation from the resting state to the intermediate state and activated state to trigger pore opening via electro-mechanical (E-M) coupling. However, the spatial and temporal details underlying E-M coupling remain elusive. Here, utilizing KV7.1's unique two open states, we report a two-stage E-M coupling mechanism in voltage-dependent gating of KV7.1 as triggered by VSD activations to the intermediate and then activated state. When the S4 segment transitions to the intermediate state, the hand-like C-terminus of the VSD-pore linker (S4-S5L) interacts with the pore in the same subunit. When S4 then proceeds to the fully-activated state, the elbow-like hinge between S4 and S4-S5L engages with the pore of the neighboring subunit to activate conductance. This two-stage hand-and-elbow gating mechanism elucidates distinct tissue-specific modulations, pharmacology, and disease pathogenesis of KV7.1, and likely applies to numerous domain-swapped KV channels.

The contribution of lipid peroxidation to membrane permeability in electropermeabilization: A molecular dynamics study.

Electroporation or electropermeabilization is a technique that enables transient increase in the cell membrane permeability by exposing cells to pulsed electric field. However, the molecular mechanisms of the long-lived cell membrane permeability, which persists on the minutes time scale after the pulse treatment, remain elusive. Experimental studies have suggested that lipid peroxidation could present a mechanism of this prolonged membrane permeabilization. In this study we make the first important step in quantifying the possible contribution of lipid peroxidation to electropermeabilization. We use free energy calculations to quantify the permeability and conductance of bilayers, containing an increasing percentage of hydroperoxide lipid derivatives, to sodium and chloride ions. We then compare our calculations to experimental measurements on electropermeabilized cells. Our results show that the permeability and conductance increase dramatically by several orders of magnitude in peroxidized bilayers. Yet this increase is not sufficient to reasonably account for the entire range of experimental measurements. Nevertheless, lipid peroxidation might be considered an important mechanism of prolonged membrane permeabilization, if exposure of cells to high voltage electric pulses leads to secondary lipid peroxidation products. Our analysis calls for experimental studies, which will determine the type and amount of lipid peroxidation products in electropermeabilized cell membranes.

Helix breaking transition in the S4 of HCN channel is critical for hyperpolarization-dependent gating.

In contrast to most voltage-gated ion channels, hyperpolarization- and cAMP gated (HCN) ion channels open on hyperpolarization. Structure-function studies show that the voltage-sensor of HCN channels are unique but the mechanisms that determine gating polarity remain poorly understood. All-atom molecular dynamics simulations (~20 µs) of HCN1 channel under hyperpolarization reveals an initial downward movement of the S4 voltage-sensor but following the transfer of last gating charge, the S4 breaks into two sub-helices with the lower sub-helix becoming parallel to the membrane. Functional studies on bipolar channels show that the gating polarity strongly correlates with helical turn propensity of the substituents at the breakpoint. Remarkably, in a proto-HCN background, the replacement of breakpoint serine with a bulky hydrophobic amino acid is sufficient to completely flip the gating polarity from inward to outward-rectifying. Our studies reveal an unexpected mechanism of inward rectification involving a linker sub-helix emerging from HCN S4 during hyperpolarization.

Determining the molecular basis of voltage sensitivity in membrane proteins.

Voltage-sensitive membrane proteins are united by their ability to transform changes in membrane potential into mechanical work. They are responsible for a spectrum of physiological processes in living organisms, including electrical signaling and cell-cycle progression. Although the mechanism of voltage-sensing has been well characterized for some membrane proteins, including voltage-gated ion channels, even the location of the voltage-sensing elements remains unknown for others. Moreover, the detection of these elements by using experimental techniques is challenging because of the diversity of membrane proteins. Here, we provide a computational approach to predict voltage-sensing elements in any membrane protein, independent of its structure or function. It relies on an estimation of the propensity of a protein to respond to changes in membrane potential. We first show that this property correlates well with voltage sensitivity by applying our approach to a set of voltage-sensitive and voltage-insensitive membrane proteins. We further show that it correctly identifies authentic voltage-sensitive residues in the voltage-sensor domain of voltage-gated ion channels. Finally, we investigate six membrane proteins for which the voltage-sensing elements have not yet been characterized and identify residues and ions that might be involved in the response to voltage. The suggested approach is fast and simple and enables a characterization of voltage sensitivity that goes beyond mere identification of charges. We anticipate that its application before mutagenesis experiments will significantly reduce the number of potential voltage-sensitive elements to be tested.

Studying Kv Channels Function using Computational Methods.

In recent years, molecular modeling techniques, combined with MD simulations, provided significant insights on voltage-gated (Kv) potassium channels intrinsic properties. Among the success stories are the highlight of molecular level details of the effects of mutations, the unraveling of several metastable intermediate states, and the influence of a particular lipid, PIP2, in the stability and the modulation of Kv channel function. These computational studies offered a detailed view that could not have been reached through experimental studies alone. With the increase of cross disciplinary studies, numerous experiments provided validation of these computational results, which endows an increase in the reliability of molecular modeling for the study of Kv channels. This chapter offers a description of the main techniques used to model Kv channels at the atomistic level.

A hypothetical molecular mechanism for TRPV1 activation that invokes rotation of an S6 asparagine.

The transient receptor potential channel vanilloid type 1 (TRPV1) is activated by a variety of endogenous and exogenous stimuli and is involved in nociception and body temperature regulation. Although the structure of TRPV1 has been experimentally determined in both the closed and open states, very little is known about its activation mechanism. In particular, the conformational changes that occur in the pore domain and result in ionic conduction have not yet been identified. Here we suggest a hypothetical molecular mechanism for TRPV1 activation, which involves rotation of a conserved asparagine in S6 from a position facing the S4-S5 linker toward the pore. This rotation is associated with hydration of the pore and dehydration of the four peripheral cavities located between each S6 and S4-S5 linker. In light of our hypothesis, we perform bioinformatics analyses of TRP and other evolutionary related ion channels, evaluate newly available structures, and reexamine previously reported water accessibility and mutagenesis experiments. These analyses provide several independent lines of evidence to support our hypothesis. Finally, we show that our proposed molecular mechanism is compatible with the prevailing theory that the selectivity filter acts as a secondary gate in TRPV1.