Functional and structural insights into activation of TRPV2 by weak acids.

Transient receptor potential (TRP) ion channels are involved in the surveillance or regulation of the acid-base balance. Here, we demonstrate that weak carbonic acids, including acetic acid, lactic acid, and CO2 activate and sensitize TRPV2 through a mechanism requiring permeation through the cell membrane. TRPV2 channels in cell-free inside-out patches maintain weak acid-sensitivity, but protons applied on either side of the membrane do not induce channel activation or sensitization. The involvement of proton modulation sites for weak acid-sensitivity was supported by the identification of titratable extracellular (Glu495, Glu561) and intracellular (His521) residues on a cryo-EM structure of rat TRPV2 (rTRPV2) treated with acetic acid. Molecular dynamics simulations as well as patch clamp experiments on mutant rTRPV2 constructs confirmed that these residues are critical for weak acid-sensitivity. We also demonstrate that the pore residue Glu609 dictates an inhibition of weak acid-induced currents by extracellular calcium. Finally, TRPV2-expression in HEK293 cells is associated with an increased weak acid-induced cytotoxicity. Together, our data provide new insights into weak acids as endogenous modulators of TRPV2.

Subtype-specific modulation of human KV 7 channels by the anticonvulsant cannabidiol through a lipid-exposed pore-domain site.

BACKGROUND AND PURPOSE: Cannabidiol (CBD) is used clinically as an anticonvulsant. Its precise mechanism of action has remained unclear. CBD was recently demonstrated to enhance the activity of the neuronal KV 7.2/7.3 channel, which may be one important contributor to CBD anticonvulsant effect. Curiously, CBD inhibits the closely related cardiac KV 7.1/KCNE1 channel. Whether and how CBD affects other KV 7 subtypes remains uninvestigated and the CBD interaction sites mediating these diverse effects remain unknown. EXPERIMENTAL APPROACH: Here, we used electrophysiology, molecular dynamics simulations, molecular docking and site-directed mutagenesis to address these questions. KEY RESULTS: We found that CBD modulates the activity of all human KV 7 subtypes and that the effects are subtype dependent. CBD enhanced the activity of KV 7.2-7.5 subtypes, seen as a V50 shift towards more negative voltages or increased maximum conductance. In contrast, CBD inhibited the KV 7.1 and KV 7.1/KCNE1 channels, seen as a V50 shift towards more positive voltages and reduced conductance. In KV 7.2 and KV 7.4, we propose a CBD interaction site at the subunit interface in the pore domain that overlaps with the interaction site of other compounds, notably the anticonvulsant retigabine. However, CBD relies on other residues for its effects than the conserved tryptophan that is critical for retigabine effects. We propose a similar, though not identical CBD site in KV 7.1, with a non-conserved phenylalanine being important. CONCLUSIONS AND IMPLICATIONS: We identify novel targets of CBD, contributing to a better understanding of CBD clinical effects and provide mechanistic insights into how CBD modulates different KV 7 subtypes.

Photogeneration of Spin Quintet Triplet-Triplet Excitations in DNA-Assembled Pentacene Stacks.

Singlet fission (SF), an exciton-doubling process observed in certain molecular semiconductors where two triplet excitons are generated from one singlet exciton, requires correctly tuned intermolecular coupling to allow separation of the two triplets to different molecular units. We explore this using DNA-encoded assembly of SF-capable pentacenes into discrete π-stacked constructs of defined size and geometry. Precise structural control is achieved via a combination of the DNA duplex formation between complementary single-stranded DNA and the local molecular geometry that directs the SF chromophores into a stable and predictable slip-stacked configuration, as confirmed by molecular dynamics (MD) modeling. Transient electron spin resonance spectroscopy revealed that within these DNA-assembled pentacene stacks, SF evolves via a bound triplet pair quintet state, which subsequently converts into free triplets. SF evolution via a long-lived quintet state sets specific requirements on intermolecular coupling, rendering the quintet spectrum and its zero-field-splitting parameters highly sensitive to intermolecular geometry. We have found that the experimental spectra and zero-field-splitting parameters are consistent with a slight systematic strain relative to the MD-optimized geometry. Thus, the transient electron spin resonance analysis is a powerful tool to test and refine the MD-derived structure models. DNA-encoded assembly of coupled semiconductor molecules allows controlled construction of electronically functional structures, but brings with it significant dynamic and polar disorders. Our findings here of efficient SF through quintet states demonstrate that these conditions still allow efficient and controlled semiconductor operation and point toward future opportunities for constructing functional optoelectronic systems.

Markov state modelling reveals heterogeneous drug-inhibition mechanism of Calmodulin.

Calmodulin (CaM) is a calcium sensor which binds and regulates a wide range of target-proteins. This implicitly enables the concentration of calcium to influence many downstream physiological responses, including muscle contraction, learning and depression. The antipsychotic drug trifluoperazine (TFP) is a known CaM inhibitor. By binding to various sites, TFP prevents CaM from associating to target-proteins. However, the molecular and state-dependent mechanisms behind CaM inhibition by drugs such as TFP are largely unknown. Here, we build a Markov state model (MSM) from adaptively sampled molecular dynamics simulations and reveal the structural and dynamical features behind the inhibitory mechanism of TFP-binding to the C-terminal domain of CaM. We specifically identify three major TFP binding-modes from the MSM macrostates, and distinguish their effect on CaM conformation by using a systematic analysis protocol based on biophysical descriptors and tools from machine learning. The results show that depending on the binding orientation, TFP effectively stabilizes features of the calcium-unbound CaM, either affecting the CaM hydrophobic binding pocket, the calcium binding sites or the secondary structure content in the bound domain. The conclusions drawn from this work may in the future serve to formulate a complete model of pharmacological modulation of CaM, which furthers our understanding of how these drugs affect signaling pathways as well as associated diseases.

Structural and dynamic mechanisms of GABAA receptor modulators with opposing activities.

γ-Aminobutyric acid type A (GABAA) receptors are pentameric ligand-gated ion channels abundant in the central nervous system and are prolific drug targets for treating anxiety, sleep disorders and epilepsy. Diverse small molecules exert a spectrum of effects on γ-aminobutyric acid type A (GABAA) receptors by acting at the classical benzodiazepine site. They can potentiate the response to GABA, attenuate channel activity, or counteract modulation by other ligands. Structural mechanisms underlying the actions of these drugs are not fully understood. Here we present two high-resolution structures of GABAA receptors in complex with zolpidem, a positive allosteric modulator and heavily prescribed hypnotic, and DMCM, a negative allosteric modulator with convulsant and anxiogenic properties. These two drugs share the extracellular benzodiazepine site at the α/γ subunit interface and two transmembrane sites at ß/α interfaces. Structural analyses reveal a basis for the subtype selectivity of zolpidem that underlies its clinical success. Molecular dynamics simulations provide insight into how DMCM switches from a negative to a positive modulator as a function of binding site occupancy. Together, these findings expand our understanding of how GABAA receptor allosteric modulators acting through a common site can have diverging activities.

Cryo-EM structure of the human Kv3.1 channel reveals gating control by the cytoplasmic T1 domain.

Kv3 channels have distinctive gating kinetics tailored for rapid repolarization in fast-spiking neurons. Malfunction of this process due to genetic variants in the KCNC1 gene causes severe epileptic disorders, yet the structural determinants for the unusual gating properties remain elusive. Here, we present cryo-electron microscopy structures of the human Kv3.1a channel, revealing a unique arrangement of the cytoplasmic tetramerization domain T1 which facilitates interactions with C-terminal axonal targeting motif and key components of the gating machinery. Additional interactions between S1/S2 linker and turret domain strengthen the interface between voltage sensor and pore domain. Supported by molecular dynamics simulations, electrophysiological and mutational analyses, we identify several residues in the S4/S5 linker which influence the gating kinetics and an electrostatic interaction between acidic residues in α6 of T1 and R449 in the pore-flanking S6T helices. These findings provide insights into gating control and disease mechanisms and may guide strategies for the design of pharmaceutical drugs targeting Kv3 channels.

Surface Electrostatics Govern the Emulsion Stability of Biomolecular Condensates.

Liquid-liquid phase separation underlies the formation of biological condensates. Physically, such systems are microemulsions that in general have a propensity to fuse and coalesce; however, many condensates persist as independent droplets in the test tube and inside cells. This stability is crucial for their function, but the physicochemical mechanisms that control the emulsion stability of condensates remain poorly understood. Here, by combining single-condensate zeta potential measurements, optical microscopy, tweezer experiments, and multiscale molecular modeling, we investigate how the nanoscale forces that sustain condensates impact their stability against fusion. By comparing peptide-RNA (PR25:PolyU) and proteinaceous (FUS) condensates, we show that a higher condensate surface charge correlates with a lower fusion propensity. Moreover, measurements of single condensate zeta potentials reveal that such systems can constitute classically stable emulsions. Taken together, these results highlight the role of passive stabilization mechanisms in protecting biomolecular condensates against coalescence.

Deoxyribonucleic Acid Encoded and Size-Defined π-Stacking of Perylene Diimides.

Natural photosystems use protein scaffolds to control intermolecular interactions that enable exciton flow, charge generation, and long-range charge separation. In contrast, there is limited structural control in current organic electronic devices such as OLEDs and solar cells. We report here the DNA-encoded assembly of π-conjugated perylene diimides (PDIs) with deterministic control over the number of electronically coupled molecules. The PDIs are integrated within DNA chains using phosphoramidite coupling chemistry, allowing selection of the DNA sequence to either side, and specification of intermolecular DNA hybridization. In this way, we have developed a "toolbox" for construction of any stacking sequence of these semiconducting molecules. We have discovered that we need to use a full hierarchy of interactions: DNA guides the semiconductors into specified close proximity, hydrophobic-hydrophilic differentiation drives aggregation of the semiconductor moieties, and local geometry and electrostatic interactions define intermolecular positioning. As a result, the PDIs pack to give substantial intermolecular π wave function overlap, leading to an evolution of singlet excited states from localized excitons in the PDI monomer to excimers with wave functions delocalized over all five PDIs in the pentamer. This is accompanied by a change in the dominant triplet forming mechanism from localized spin-orbit charge transfer mediated intersystem crossing for the monomer toward a delocalized excimer process for the pentamer. Our modular DNA-based assembly reveals real opportunities for the rapid development of bespoke semiconductor architectures with molecule-by-molecule precision.

Regulation of a pentameric ligand-gated ion channel by a semiconserved cationic lipid-binding site.

Pentameric ligand-gated ion channels (pLGICs) are crucial mediators of electrochemical signal transduction in various organisms from bacteria to humans. Lipids play an important role in regulating pLGIC function, yet the structural bases for specific pLGIC-lipid interactions remain poorly understood. The bacterial channel ELIC recapitulates several properties of eukaryotic pLGICs, including activation by the neurotransmitter GABA and binding and modulation by lipids, offering a simplified model system for structure-function relationship studies. In this study, functional effects of noncanonical amino acid substitution of a potential lipid-interacting residue (W206) at the top of the M1-helix, combined with detergent interactions observed in recent X-ray structures, are consistent with this region being the location of a lipid-binding site on the outward face of the ELIC transmembrane domain. Coarse-grained and atomistic molecular dynamics simulations revealed preferential binding of lipids containing a positive charge, particularly involving interactions with residue W206, consistent with cation-π binding. Polar contacts from other regions of the protein, particularly M3 residue Q264, further support lipid binding via headgroup ester linkages. Aromatic residues were identified at analogous sites in a handful of eukaryotic family members, including the human GABAA receptor ε subunit, suggesting conservation of relevant interactions in other evolutionary branches. Further mutagenesis experiments indicated that mutations at this site in ε-containing GABAA receptors can change the apparent affinity of the agonist response to GABA, suggesting a potential role of this site in channel gating. In conclusion, this work details type-specific lipid interactions, which adds to our growing understanding of how lipids modulate pLGICs.

Reentrant liquid condensate phase of proteins is stabilized by hydrophobic and non-ionic interactions.

Liquid-liquid phase separation of proteins underpins the formation of membraneless compartments in living cells. Elucidating the molecular driving forces underlying protein phase transitions is therefore a key objective for understanding biological function and malfunction. Here we show that cellular proteins, which form condensates at low salt concentrations, including FUS, TDP-43, Brd4, Sox2, and Annexin A11, can reenter a phase-separated regime at high salt concentrations. By bringing together experiments and simulations, we demonstrate that this reentrant phase transition in the high-salt regime is driven by hydrophobic and non-ionic interactions, and is mechanistically distinct from the low-salt regime, where condensates are additionally stabilized by electrostatic forces. Our work thus sheds light on the cooperation of hydrophobic and non-ionic interactions as general driving forces in the condensation process, with important implications for aberrant function, druggability, and material properties of biomolecular condensates.

Protein disorder-to-order transition enhances the nucleosome-binding affinity of H1.

Intrinsically disordered proteins are crucial elements of chromatin heterogenous organization. While disorder in the histone tails enables a large variation of inter-nucleosome arrangements, disorder within the chromatin-binding proteins facilitates promiscuous binding to a wide range of different molecular targets, consistent with structural heterogeneity. Among the partially disordered chromatin-binding proteins, the H1 linker histone influences a myriad of chromatin characteristics including compaction, nucleosome spacing, transcription regulation, and the recruitment of other chromatin regulating proteins. Although it is now established that the long C-terminal domain (CTD) of H1 remains disordered upon nucleosome binding and that such disorder favours chromatin fluidity, the structural behaviour and thereby the role/function of the N-terminal domain (NTD) within chromatin is yet unresolved. On the basis of microsecond-long parallel-tempering metadynamics and temperature-replica exchange atomistic molecular dynamics simulations of different H1 NTD subtypes, we demonstrate that the NTD is completely unstructured in solution but undergoes an important disorder-to-order transition upon nucleosome binding: it forms a helix that enhances its DNA binding ability. Further, we show that the helical propensity of the H1 NTD is subtype-dependent and correlates with the experimentally observed binding affinity of H1 subtypes, suggesting an important functional implication of this disorder-to-order transition.

Emergence of chromatin hierarchical loops from protein disorder and nucleosome asymmetry.

Protein flexibility and disorder is emerging as a crucial modulator of chromatin structure. Histone tail disorder enables transient binding of different molecules to the nucleosomes, thereby promoting heterogeneous and dynamic internucleosome interactions and making possible recruitment of a wide-range of regulatory and remodeling proteins. On the basis of extensive multiscale modeling we reveal the importance of linker histone H1 protein disorder for chromatin hierarchical looping. Our multiscale approach bridges microsecond-long bias-exchange metadynamics molecular dynamics simulations of atomistic 211-bp nucleosomes with coarse-grained Monte Carlo simulations of 100-nucleosome systems. We show that the long C-terminal domain (CTD) of H1-a ubiquitous nucleosome-binding protein-remains disordered when bound to the nucleosome. Notably, such CTD disorder leads to an asymmetric and dynamical nucleosome conformation that promotes chromatin structural flexibility and establishes long-range hierarchical loops. Furthermore, the degree of condensation and flexibility of H1 can be fine-tuned, explaining chromosomal differences of interphase versus metaphase states that correspond to partial and hyperphosphorylated H1, respectively. This important role of H1 protein disorder in large-scale chromatin organization has a wide range of biological implications.

Waterdock 2.0: Water placement prediction for Holo-structures with a pymol plugin.

Water is often found to mediate interactions between a ligand and a protein. It can play a significant role in orientating the ligand within a binding pocket and contribute to the free energy of binding. It would thus be extremely useful to be able to accurately predict the position and orientation of water molecules within a binding pocket. Recently, we developed the WaterDock protocol that was able to predict 97% of the water molecules in a test set. However, this approach generated false positives at a rate of over 20% in most cases and whilst this might be acceptable for some applications, in high throughput scenarios this is not desirable. Here we tackle this problem via the inclusion of knowledge regarding the solvation structure of ligand functional groups. We call this new protocol WaterDock2 and demonstrate that this protocol maintains a similar true positive rate to the original implementation but is capable of reducing the false-positive rate by over 50%. To improve the usability of the method, we have also developed a plugin for the popular graphics program PyMOL. The plugin also contains an implementation of the original WaterDock.

Multi-scale molecular dynamics study of cholera pentamer binding to a GM1-phospholipid membrane.

The AB5 type toxin produced by the Vibrio cholerae bacterium is the causative agent of the cholera disease. The cholera toxin (CT) has been shown to bind specifically to GM1 glycolipids on the membrane surface. This binding of CT to the membrane is the initial step in its endocytosis and has been postulated to cause significant disruption to the membrane structure. In this work, we have carried out a combination of coarse-grain and atomistic simulations to study the binding of CT to a membrane modelled as an asymmetrical GM1-DPPC bilayer. Simulation results indicate that the toxin binds to the membrane through only three of its five B subunits, in effect resulting in a tilted bound configuration. Additionally, the binding of the CT can increase the area per lipid of GM1 leaflet, which in turn can cause the membrane regions interacting with the bound subunits to experience significant bilayer thinning and lipid tail disorder across both the leaflets.

The solvation structure of alprazolam.

Alprazolam is a benzodiazepine that is commonly prescribed for the treatment of anxiety and other related disorders. Like other benzodiazepines, it is thought to exert its effect through interaction with GABAA receptors. However, it has also been described as a potent and selective protein interaction inhibitor of bromodomain and extra-terminal (BET) proteins. Indeed, the only crystal structure of alprazolam bound to a protein is a complex between alprazolam and the BRD4 bromodomain. The structure shows that the complex also involves many water interactions that mediate contacts between the drug and the protein, a scenario that exists in many drug-protein complexes. How such waters relate to solvation patterns of small molecules may improve our understanding of what dictates their appearance or absence in bridging positions within complexes and thus will be important in terms of future rational drug-design. Here, we use neutron diffraction in conjunction with molecular dynamics simulations to provide a detailed analysis of how water molecules interact with alprazolam in methanol/water mixtures. The agreement between the neutron diffraction and the molecular dynamics is extremely good. We discuss the results in the context of drug design.

Content-based image retrieval of digitized histopathology in boosted spectrally embedded spaces.

CONTEXT: Content-based image retrieval (CBIR) systems allow for retrieval of images from within a database that are similar in visual content to a query image. This is useful for digital pathology, where text-based descriptors alone might be inadequate to accurately describe image content. By representing images via a set of quantitative image descriptors, the similarity between a query image with respect to archived, annotated images in a database can be computed and the most similar images retrieved. Recently, non-linear dimensionality reduction methods have become popular for embedding high-dimensional data into a reduced-dimensional space while preserving local object adjacencies, thereby allowing for object similarity to be determined more accurately in the reduced-dimensional space. However, most dimensionality reduction methods implicitly assume, in computing the reduced-dimensional representation, that all features are equally important. AIMS: In this paper we present boosted spectral embedding(BoSE), which utilizes a boosted distance metric to selectively weight individual features (based on training data) to subsequently map the data into a reduced-dimensional space. SETTINGS AND DESIGN: BoSE is evaluated against spectral embedding (SE) (which employs equal feature weighting) in the context of CBIR of digitized prostate and breast cancer histopathology images. MATERIALS AND METHODS: The following datasets, which were comprised of a total of 154 hematoxylin and eosin stained histopathology images, were used: (1) Prostate cancer histopathology (benign vs. malignant), (2) estrogen receptor (ER) + breast cancer histopathology (low vs. high grade), and (3) HER2+ breast cancer histopathology (low vs. high levels of lymphocytic infiltration). STATISTICAL ANALYSIS USED: We plotted and calculated the area under precision-recall curves (AUPRC) and calculated classification accuracy using the Random Forest classifier. RESULTS: BoSE outperformed SE both in terms of CBIR-based (area under the precision-recall curve) and classifier-based (classification accuracy) on average across all of the dimensions tested for all three datasets: (1) Prostate cancer histopathology (AUPRC: BoSE = 0.79, SE = 0.63; Accuracy: BoSE = 0.93, SE = 0.80), (2) ER + breast cancer histopathology (AUPRC: BoSE = 0.79, SE = 0.68; Accuracy: BoSE = 0.96, SE = 0.96), and (3) HER2+ breast cancer histopathology (AUPRC: BoSE = 0.54, SE = 0.44; Accuracy: BoSE = 0.93, SE = 0.91). CONCLUSION: Our results suggest that BoSE could serve as an important tool for CBIR and classification of high-dimensional biomedical data.

Coarse-grain molecular dynamics study of fullerene transport across a cell membrane.

The study of the ability of drug molecules to enter cells through the membrane is of vital importance in the field of drug delivery. In cases where the transport of the drug molecules through the membrane is not easily accomplishable, other carrier molecules are used. Spherical fullerene molecules have been postulated as potential carriers of highly hydrophilic drugs across the plasma membrane. Here, we report the coarse-grain molecular dynamics study of the translocation of C60 fullerene and its derivatives across a cell membrane modeled as a 1,2-distearoyl-sn-glycero-3-phosphocholine bilayer. Simulation results indicate that pristine fullerene molecules enter the bilayer quickly and reside within it. The addition of polar functionalized groups makes the fullerenes less likely to reside within the bilayer but increases their residence time in bulk water. Addition of polar functional groups to one half of the fullerene surface, in effect creating a Janus particle, offers the most promise in developing fullerene models that can achieve complete translocation through the membrane bilayer.