Directed Evolution of a Selective and Sensitive Serotonin Sensor via Machine Learning.

Serotonin plays a central role in cognition and is the target of most pharmaceuticals for psychiatric disorders. Existing drugs have limited efficacy; creation of improved versions will require better understanding of serotonergic circuitry, which has been hampered by our inability to monitor serotonin release and transport with high spatial and temporal resolution. We developed and applied a binding-pocket redesign strategy, guided by machine learning, to create a high-performance, soluble, fluorescent serotonin sensor (iSeroSnFR), enabling optical detection of millisecond-scale serotonin transients. We demonstrate that iSeroSnFR can be used to detect serotonin release in freely behaving mice during fear conditioning, social interaction, and sleep/wake transitions. We also developed a robust assay of serotonin transporter function and modulation by drugs. We expect that both machine-learning-guided binding-pocket redesign and iSeroSnFR will have broad utility for the development of other sensors and in vitro and in vivo serotonin detection, respectively.

H3K4me1 recruits DNA repair proteins in plants.

DNA repair proteins can be recruited by their histone reader domains to specific epigenomic features, with consequences on intragenomic mutation rate variation. Here, we investigated H3K4me1-associated hypomutation in plants. We first examined 2 proteins which, in plants, contain Tudor histone reader domains: PRECOCIOUS DISSOCIATION OF SISTERS 5 (PDS5C), involved in homology-directed repair, and MUTS HOMOLOG 6 (MSH6), a mismatch repair protein. The MSH6 Tudor domain of Arabidopsis (Arabidopsis thaliana) binds to H3K4me1 as previously demonstrated for PDS5C, which localizes to H3K4me1-rich gene bodies and essential genes. Mutations revealed by ultradeep sequencing of wild-type and msh6 knockout lines in Arabidopsis show that functional MSH6 is critical for the reduced rate of single-base substitution (SBS) mutations in gene bodies and H3K4me1-rich regions. We explored the breadth of these mechanisms among plants by examining a large rice (Oryza sativa) mutation data set. H3K4me1-associated hypomutation is conserved in rice as are the H3K4me1-binding residues of MSH6 and PDS5C Tudor domains. Recruitment of DNA repair proteins by H3K4me1 in plants reveals convergent, but distinct, epigenome-recruited DNA repair mechanisms from those well described in humans. The emergent model of H3K4me1-recruited repair in plants is consistent with evolutionary theory regarding mutation modifier systems and offers mechanistic insight into intragenomic mutation rate variation in plants.

Elucidation of a dynamic interplay between a beta-2 adrenergic receptor, its agonist, and stimulatory G protein.

G protein-coupled receptors (GPCRs) represent the largest group of membrane receptors for transmembrane signal transduction. Ligand-induced activation of GPCRs triggers G protein activation followed by various signaling cascades. Understanding the structural and energetic determinants of ligand binding to GPCRs and GPCRs to G proteins is crucial to the design of pharmacological treatments targeting specific conformations of these proteins to precisely control their signaling properties. In this study, we focused on interactions of a prototypical GPCR, beta-2 adrenergic receptor (ß2AR), with its endogenous agonist, norepinephrine (NE), and the stimulatory G protein (Gs). Using molecular dynamics (MD) simulations, we demonstrated the stabilization of cationic NE, NE(+), binding to ß2AR by Gs protein recruitment, in line with experimental observations. We also captured the partial dissociation of the ligand from ß2AR and the conformational interconversions of Gs between closed and open conformations in the NE(+)-ß2AR-Gs ternary complex while it is still bound to the receptor. The variation of NE(+) binding poses was found to alter Gs α subunit (Gsα) conformational transitions. Our simulations showed that the interdomain movement and the stacking of Gsα α1 and α5 helices are significant for increasing the distance between the Gsα and ß2AR, which may indicate a partial dissociation of Gsα The distance increase commences when Gsα is predominantly in an open state and can be triggered by the intracellular loop 3 (ICL3) of ß2AR interacting with Gsα, causing conformational changes of the α5 helix. Our results help explain molecular mechanisms of ligand and GPCR-mediated modulation of G protein activation.

Energetic landscape of polycystin channel gating.

Members of the polycystin family (PKD2 and PKD2L1) of transient receptor potential (TRP) channels conduct Ca2+ and depolarizing monovalent cations. Variants in PKD2 cause autosomal dominant polycystic kidney disease (ADPKD) in humans, whereas loss of PKD2L1 expression causes seizure susceptibility in mice. Understanding structural and functional regulation of these channels will provide the basis for interpreting their molecular dysregulation in disease states. However, the complete structures of polycystins are unresolved, as are the conformational changes regulating their conductive states. To provide a holistic understanding of the polycystin gating cycle, we use computational prediction tools to model missing PKD2L1 structural motifs and evaluate more than 150 mutations in an unbiased mutagenic functional screen of the entire pore module. Our results provide an energetic landscape of the polycystin pore, which enumerates gating sensitive sites and interactions required for opening, inactivation, and subsequent desensitization. These findings identify the external pore helices and specific cross-domain interactions as critical structural regulators controlling the polycystin ion channel conductive and nonconductive states.

Opening of capsaicin receptor TRPV1 is stabilized equally by its four subunits.

Capsaicin receptor TRPV1 is a nociceptor for vanilloid molecules, such as capsaicin and resiniferatoxin (RTX). Even though cryo-EM structures of TRPV1 in complex with these molecules are available, how their binding energetically favors the open conformation is not known. Here, we report an approach to control the number of bound RTX molecules (0-4) in functional rat TRPV1. The approach allowed direct measurements of each of the intermediate open states under equilibrium conditions at both macroscopic and single-molecule levels. We found that RTX binding to each of the four subunits contributes virtually the same activation energy, which we estimated to be 1.70 to 1.86 kcal/mol and found to arise predominately from destabilizing the closed conformation. We further showed that sequential bindings of RTX increase open probability without altering single-channel conductance, confirming that there is likely a single open-pore conformation for TRPV1 activated by RTX.

α-Actinin-1 promotes activity of the L-type Ca2+ channel Cav 1.2.

The L-type Ca2+ channel CaV 1.2 governs gene expression, cardiac contraction, and neuronal activity. Binding of α-actinin to the IQ motif of CaV 1.2 supports its surface localization and postsynaptic targeting in neurons. We report a bi-functional mechanism that restricts CaV 1.2 activity to its target sites. We solved separate NMR structures of the IQ motif (residues 1,646-1,664) bound to α-actinin-1 and to apo-calmodulin (apoCaM). The CaV 1.2 K1647A and Y1649A mutations, which impair α-actinin-1 but not apoCaM binding, but not the F1658A and K1662E mutations, which impair apoCaM but not α-actinin-1 binding, decreased single-channel open probability, gating charge movement, and its coupling to channel opening. Thus, α-actinin recruits CaV 1.2 to defined surface regions and simultaneously boosts its open probability so that CaV 1.2 is mostly active when appropriately localized.

A Novel Peptide Prevents Enterotoxin- and Inflammation-Induced Intestinal Fluid Secretion by Stimulating Sodium-Hydrogen Exchanger 3 Activity.

BACKGROUND & AIMS: Acute diarrheal diseases are the second most common cause of infant mortality in developing countries. This is contributed to by lack of effective drug therapy that shortens the duration or lessens the volume of diarrhea. The epithelial brush border sodium (Na+)/hydrogen (H+) exchanger 3 (NHE3) accounts for a major component of intestinal Na+ absorption and is inhibited in most diarrheas. Because increased intestinal Na+ absorption can rehydrate patients with diarrhea, NHE3 has been suggested as a potential druggable target for drug therapy for diarrhea. METHODS: A peptide (sodium-hydrogen exchanger 3 stimulatory peptide [N3SP]) was synthesized to mimic the part of the NHE3 C-terminus that forms a multiprotein complex that inhibits NHE3 activity. The effect of N3SP on NHE3 activity was evaluated in NHE3-transfected fibroblasts null for other plasma membrane NHEs, a human colon cancer cell line that models intestinal absorptive enterocytes (Caco-2/BBe), human enteroids, and mouse intestine in vitro and in vivo. N3SP was delivered into cells via a hydrophobic fluorescent maleimide or nanoparticles. RESULTS: N3SP uptake stimulated NHE3 activity at nmol/L concentrations under basal conditions and partially reversed the reduced NHE3 activity caused by elevated adenosine 3',5'-cyclic monophosphate, guanosine 3',5'-cyclic monophosphate, and Ca2+ in cell lines and in in vitro mouse intestine. N3SP also stimulated intestinal fluid absorption in the mouse small intestine in vivo and prevented cholera toxin-, Escherichia coli heat-stable enterotoxin-, and cluster of differentiation 3 inflammation-induced fluid secretion in a live mouse intestinal loop model. CONCLUSIONS: These findings suggest pharmacologic stimulation of NHE3 activity as an efficacious approach for the treatment of moderate/severe diarrheal diseases.

Structural modeling of peptide toxin-ion channel interactions using RosettaDock.

Voltage-gated ion channels play essential physiological roles in action potential generation and propagation. Peptidic toxins from animal venoms target ion channels and provide useful scaffolds for the rational design of novel channel modulators with enhanced potency and subtype selectivity. Despite recent progress in obtaining experimental structures of peptide toxin-ion channel complexes, structural determination of peptide toxins bound to ion channels in physiologically important states remains challenging. Here we describe an application of RosettaDock approach to the structural modeling of peptide toxins interactions with ion channels. We tested this approach on 10 structures of peptide toxin-ion channel complexes and demonstrated that it can sample near-native structures in all tested cases. Our approach will be useful for improving the understanding of the molecular mechanism of natural peptide toxin modulation of ion channel gating and for the structural modeling of novel peptide-based ion channel modulators.

Toggle switch residues control allosteric transitions in bacterial adhesins by participating in a concerted repacking of the protein core.

Critical molecular events that control conformational transitions in most allosteric proteins are ill-defined. The mannose-specific FimH protein of Escherichia coli is a prototypic bacterial adhesin that switches from an 'inactive' low-affinity state (LAS) to an 'active' high-affinity state (HAS) conformation allosterically upon mannose binding and mediates shear-dependent catch bond adhesion. Here we identify a novel type of antibody that acts as a kinetic trap and prevents the transition between conformations in both directions. Disruption of the allosteric transitions significantly slows FimH's ability to associate with mannose and blocks bacterial adhesion under dynamic conditions. FimH residues critical for antibody binding form a compact epitope that is located away from the mannose-binding pocket and is structurally conserved in both states. A larger antibody-FimH contact area is identified by NMR and contains residues Leu-34 and Val-35 that move between core-buried and surface-exposed orientations in opposing directions during the transition. Replacement of Leu-34 with a charged glutamic acid stabilizes FimH in the LAS conformation and replacement of Val-35 with glutamic acid traps FimH in the HAS conformation. The antibody is unable to trap the conformations if Leu-34 and Val-35 are replaced with a less bulky alanine. We propose that these residues act as molecular toggle switches and that the bound antibody imposes a steric block to their reorientation in either direction, thereby restricting concerted repacking of side chains that must occur to enable the conformational transition. Residues homologous to the FimH toggle switches are highly conserved across a diverse family of fimbrial adhesins. Replacement of predicted switch residues reveals that another E. coli adhesin, galactose-specific FmlH, is allosteric and can shift from an inactive to an active state. Our study shows that allosteric transitions in bacterial adhesins depend on toggle switch residues and that an antibody that blocks the switch effectively disables adhesive protein function.

A Computational Pipeline to Predict Cardiotoxicity: From the Atom to the Rhythm.

RATIONALE: Drug-induced proarrhythmia is so tightly associated with prolongation of the QT interval that QT prolongation is an accepted surrogate marker for arrhythmia. But QT interval is too sensitive a marker and not selective, resulting in many useful drugs eliminated in drug discovery. OBJECTIVE: To predict the impact of a drug from the drug chemistry on the cardiac rhythm. METHODS AND RESULTS: In a new linkage, we connected atomistic scale information to protein, cell, and tissue scales by predicting drug-binding affinities and rates from simulation of ion channel and drug structure interactions and then used these values to model drug effects on the hERG channel. Model components were integrated into predictive models at the cell and tissue scales to expose fundamental arrhythmia vulnerability mechanisms and complex interactions underlying emergent behaviors. Human clinical data were used for model framework validation and showed excellent agreement, demonstrating feasibility of a new approach for cardiotoxicity prediction. CONCLUSIONS: We present a multiscale model framework to predict electrotoxicity in the heart from the atom to the rhythm. Novel mechanistic insights emerged at all scales of the system, from the specific nature of proarrhythmic drug interaction with the hERG channel, to the fundamental cellular and tissue-level arrhythmia mechanisms. Applications of machine learning indicate necessary and sufficient parameters that predict arrhythmia vulnerability. We expect that the model framework may be expanded to make an impact in drug discovery, drug safety screening for a variety of compounds and targets, and in a variety of regulatory processes.

Opening TRPP2 (PKD2L1) requires the transfer of gating charges.

The opening of voltage-gated ion channels is initiated by transfer of gating charges that sense the electric field across the membrane. Although transient receptor potential ion channels (TRP) are members of this family, their opening is not intrinsically linked to membrane potential, and they are generally not considered voltage gated. Here we demonstrate that TRPP2, a member of the polycystin subfamily of TRP channels encoded by the PKD2L1 gene, is an exception to this rule. TRPP2 borrows a biophysical riff from canonical voltage-gated ion channels, using 2 gating charges found in its fourth transmembrane segment (S4) to control its conductive state. Rosetta structural prediction demonstrates that the S4 undergoes ∼3- to 5-Å transitional and lateral movements during depolarization, which are coupled to opening of the channel pore. Here both gating charges form state-dependent cation-π interactions within the voltage sensor domain (VSD) during membrane depolarization. Our data demonstrate that the transfer of a single gating charge per channel subunit is requisite for voltage, temperature, and osmotic swell polymodal gating of TRPP2. Taken together, we find that irrespective of stimuli, TRPP2 channel opening is dependent on activation of its VSDs.

Structural basis for antiarrhythmic drug interactions with the human cardiac sodium channel.

The human voltage-gated sodium channel, hNaV1.5, is responsible for the rapid upstroke of the cardiac action potential and is target for antiarrhythmic therapy. Despite the clinical relevance of hNaV1.5-targeting drugs, structure-based molecular mechanisms of promising or problematic drugs have not been investigated at atomic scale to inform drug design. Here, we used Rosetta structural modeling and docking as well as molecular dynamics simulations to study the interactions of antiarrhythmic and local anesthetic drugs with hNaV1.5. These calculations revealed several key drug binding sites formed within the pore lumen that can simultaneously accommodate up to two drug molecules. Molecular dynamics simulations identified a hydrophilic access pathway through the intracellular gate and a hydrophobic access pathway through a fenestration between DIII and DIV. Our results advance the understanding of molecular mechanisms of antiarrhythmic and local anesthetic drug interactions with hNaV1.5 and will be useful for rational design of novel therapeutics.

Molecular determinants of pro-arrhythmia proclivity of d- and l-sotalol via a multi-scale modeling pipeline.

Drug isomers may differ in their proarrhythmia risk. An interesting example is the drug sotalol, an antiarrhythmic drug comprising d- and l- enantiomers that both block the hERG cardiac potassium channel and confer differing degrees of proarrhythmic risk. We developed a multi-scale in silico pipeline focusing on hERG channel - drug interactions and used it to probe and predict the mechanisms of pro-arrhythmia risks of the two enantiomers of sotalol. Molecular dynamics (MD) simulations predicted comparable hERG channel binding affinities for d- and l-sotalol, which were validated with electrophysiology experiments. MD derived thermodynamic and kinetic parameters were used to build multi-scale functional computational models of cardiac electrophysiology at the cell and tissue scales. Functional models were used to predict inactivated state binding affinities to recapitulate electrocardiogram (ECG) QT interval prolongation observed in clinical data. Our study demonstrates how modeling and simulation can be applied to predict drug effects from the atom to the rhythm for dl-sotalol and also increased proarrhythmia proclivity of d- vs. l-sotalol when accounting for stereospecific beta-adrenergic receptor blocking.

Pathogenic effects of agrin V1727F mutation are isoform specific and decrease its expression and affinity for HSPGs and LRP4.

Agrin is a large extracellular matrix protein whose isoforms differ in their tissue distribution and function. Motoneuron-derived y+z+ agrin regulates the formation of the neuromuscular junction (NMJ), while y-z- agrin is widely expressed and has diverse functions. Previously we identified a missense mutation (V1727F) in the second laminin globular (LG2) domain of agrin that causes severe congenital myasthenic syndrome. Here, we define pathogenic effects of the agrin V1727F mutation that account for the profound dysfunction of the NMJ. First, by expressing agrin variants in heterologous cells, we show that the V1727F mutation reduces the secretion of y+z+ agrin compared to wild type, whereas it has no effect on the secretion of y-z- agrin. Second, we find that the V1727F mutation significantly impairs binding of y+z+ agrin to both heparin and the low-density lipoprotein receptor-related protein 4 (LRP4) coreceptor. Third, molecular modeling of the LG2 domain suggests that the V1727F mutation primarily disrupts the y splice insert, and consistent with this we find that it partially occludes the contribution of the y splice insert to agrin binding to heparin and LRP4. Together, these findings identify several pathogenic effects of the V1727F mutation that reduce its expression and ability to bind heparan sulfate proteoglycan and LRP4 coreceptors involved in the muscle-specific kinase signaling pathway. These defects primarily impair the function of neural y+z+ agrin and combine to cause a severe CMS phenotype, whereas y-z- agrin function in other tissues appears preserved.

In vivo optophysiology reveals that G-protein activation triggers osmotic swelling and increased light scattering of rod photoreceptors.

The light responses of rod and cone photoreceptors have been studied electrophysiologically for decades, largely with ex vivo approaches that disrupt the photoreceptors' subretinal microenvironment. Here we report the use of optical coherence tomography (OCT) to measure light-driven signals of rod photoreceptors in vivo. Visible light stimulation over a 200-fold intensity range caused correlated rod outer segment (OS) elongation and increased light scattering in wild-type mice, but not in mice lacking the rod G-protein alpha subunit, transducin (Gαt), revealing these responses to be triggered by phototransduction. For stimuli that photoactivated one rhodopsin per Gαt the rod OS swelling response reached a saturated elongation of 10.0 ± 2.1%, at a maximum rate of 0.11% s-1 Analyzing swelling as osmotically driven water influx, we find the H2O membrane permeability of the rod OS to be (2.6 ± 0.4) × 10-5 cm⋅s-1, comparable to that of other cells lacking aquaporin expression. Application of Van't Hoff's law reveals that complete activation of phototransduction generates a potentially harmful 20% increase in OS osmotic pressure. The increased backscattering from the base of the OS is explained by a model combining cytoplasmic swelling, translocation of dissociated G-protein subunits from the disc membranes into the cytoplasm, and a relatively higher H2O permeability of nascent discs in the basal rod OS. Translocation of phototransduction components out of the OS may protect rods from osmotic stress, which could be especially harmful in disease conditions that affect rod OS structural integrity.

A distinct structural mechanism underlies TRPV1 activation by piperine.

Piperine, the principle pungent compound in black peppers, is known to activate the capsaicin receptor TRPV1 ion channel. How piperine interacts with the channel protein, however, remains unclear. Here we show that piperine binds to the same ligand-binding pocket as capsaicin but in different poses. There was no detectable detrimental effect when T551 and E571, two major sites known to form hydrogen bond with capsaicin, were mutated to a hydrophobic amino acid. Computational structural modeling suggested that piperine makes interactions with multiple amino acids within the ligand binding pocket, including T671 on the pore-forming S6 segment. Mutations of this residue could substantially reduce or even eliminate piperine-induced activation, confirming that T671 is an important site. Our results suggest that the bound piperine may directly interact with the pore-forming S6 segment to induce channel opening. These findings help to explain why piperine is a weak agonist, and may guide future efforts to develop novel pharmaceutical reagents targeting TRPV1.

Rational design and validation of a vanilloid-sensitive TRPV2 ion channel.

Vanilloids activation of TRPV1 represents an excellent model system of ligand-gated ion channels. Recent studies using cryo-electron microcopy (cryo-EM), computational analysis, and functional quantification revealed the location of capsaicin-binding site and critical residues mediating ligand-binding and channel activation. Based on these new findings, here we have successfully introduced high-affinity binding of capsaicin and resiniferatoxin to the vanilloid-insensitive TRPV2 channel, using a rationally designed minimal set of four point mutations (F467S-S498F-L505T-Q525E, termed TRPV2_Quad). We found that binding of resiniferatoxin activates TRPV2_Quad but the ligand-induced open state is relatively unstable, whereas binding of capsaicin to TRPV2_Quad antagonizes resiniferatoxin-induced activation likely through competition for the same binding sites. Using Rosetta-based molecular docking, we observed a common structural mechanism underlying vanilloids activation of TRPV1 and TRPV2_Quad, where the ligand serves as molecular "glue" that bridges the S4-S5 linker to the S1-S4 domain to open these channels. Our analysis revealed that capsaicin failed to activate TRPV2_Quad likely due to structural constraints preventing such bridge formation. These results not only validate our current working model for capsaicin activation of TRPV1 but also should help guide the design of drug candidate compounds for this important pain sensor.

Gain-of-function mutation of a voltage-gated sodium channel NaV1.7 associated with peripheral pain and impaired limb development.

Dominant mutations in voltage-gated sodium channel NaV1.7 cause inherited erythromelalgia, a debilitating pain disorder characterized by severe burning pain and redness of the distal extremities. NaV1.7 is preferentially expressed within peripheral sensory and sympathetic neurons. Here, we describe a novel NaV1.7 mutation in an 11-year-old male with underdevelopment of the limbs, recurrent attacks of burning pain with erythema, and swelling in his feet and hands. Frequency and duration of the episodes gradually increased with age, and relief by cooling became less effective. The patient's sister had short stature and reported similar complaints of erythema and burning pain, but with less intensity. Genetic analysis revealed a novel missense mutation in NaV1.7 (2567G>C; p.Gly856Arg) in both siblings. The G856R mutation, located within the DII/S4-S5 linker of the channel, substitutes a highly conserved non-polar glycine by a positively charged arginine. Voltage-clamp analysis of G856R currents revealed that the mutation hyperpolarized (-11.2 mV) voltage dependence of activation and slowed deactivation but did not affect fast inactivation, compared with wild-type channels. A mutation of Gly-856 to aspartic acid was previously found in a family with limb pain and limb underdevelopment, and its functional assessment showed hyperpolarized activation, depolarized fast inactivation, and increased ramp current. Structural modeling using the Rosetta computational modeling suite provided structural clues to the divergent effects of the substitution of Gly-856 by arginine and aspartic acid. Although the proexcitatory changes in gating properties of G856R contribute to the pathophysiology of inherited erythromelalgia, the link to limb underdevelopment is not well understood.

A novel tarantula toxin stabilizes the deactivated voltage sensor of bacterial sodium channel.

Voltage-gated sodium channels (NaVs) are activated by transiting the voltage sensor from the deactivated to the activated state. The crystal structures of several bacterial NaVs have captured the voltage sensor module (VSM) in an activated state, but structure of the deactivated voltage sensor remains elusive. In this study, we sought to identify peptide toxins stabilizing the deactivated VSM of bacterial NaVs. We screened fractions from several venoms and characterized a cystine knot toxin called JZTx-27 from the venom of tarantula Chilobrachys jingzhao as a high-affinity antagonist of the prokaryotic NaVs NsVBa (nonselective voltage-gated Bacillus alcalophilus) and NaChBac (bacterial sodium channel from Bacillus halodurans) (IC50 = 112 nM and 30 nM, respectively). JZTx-27 was more efficacious at weaker depolarizing voltages and significantly slowed the activation but accelerated the deactivation of NsVBa, whereas the local anesthetic drug lidocaine was shown to antagonize NsVBa without affecting channel gating. Mutation analysis confirmed that JZTx-27 bound to S3-4 linker of NsVBa, with F98 being the critical residue in determining toxin affinity. All electrophysiological data and in silico analysis suggested that JZTx-27 trapped VSM of NsVBa in one of the deactivated states. In mammalian NaVs, JZTx-27 preferably inhibited the inactivation of NaV1.5 by targeting the fourth transmembrane domain. To our knowledge, this is the first report of peptide antagonist for prokaryotic NaVs. More important, we proposed that JZTx-27 stabilized the NsVBa VSM in the deactivated state and may be used as a probe to determine the structure of the deactivated VSM of NaVs.-Tang, C., Zhou, X., Nguyen, P. T., Zhang, Y., Hu, Z., Zhang, C., Yarov-Yarovoy, V., DeCaen, P. G., Liang, S., Liu, Z. A novel tarantula toxin stabilizes the deactivated voltage sensor of bacterial sodium channel.