Undiagnosed Disease Network collaborative approach in diagnosing rare disease in a patient with a mosaic CACNA1D variant

The Undiagnosed Disease Network (UDN) is comprised of clinical and research experts collaborating to diagnose rare disease. The UDN is funded by the National Institutes of Health and includes 12 different clinical sites (About Us, 2022). Here we highlight the success of collaborative efforts within the UDN Clinical Site at Vanderbilt University Medical Center (VUMC) in utilizing a cohort of experts in bioinformatics, structural biology, and genetics specialists in diagnosing rare disease. Our UDN team identified a de novo mosaic CACNA1D variant c.2299T>C in a 5-year-old female with a history of global developmental delay, dystonia, dyskinesis, and seizures. Using a collaborative multidisciplinary approach, our VUMC UDN team diagnosed the participant with Primary Aldosteronism, Seizures, and Neurologic abnormalities (PASNA) OMIM: 615474 due to a rare mosaic CACNA1D variant (O'Neill, 2013). Interestingly, this patient was mosaic, a phenotypic trait previously unreported in PASNA cases. This report highlights the importance of a multidisciplinary approach in diagnosing rare disease.

The novel non-hallucinogenic compound DM506 (3-methyl-1,2,3,4,5,6-hexahydroazepino[4,5-b]indole) induces sedative- and anxiolytic-like activity in mice by a mechanism involving 5-HT2A receptor activation

The anxiolytic and sedative-like effects of 3-methyl-1,2,3,4,5,6-hexahydroazepino[4,5-b]indole (DM506), a non-hallucinogenic compound derived from ibogamine, were studied in mice. The behavioral effects were examined using Elevated O-maze and novelty suppressed feeding (NSFT) tests, open field test, and loss of righting reflex (LORR) test. The results showed that 15 mg/kg DM506 induced acute and long-lasting anxiolytic-like activity in naive and stressed/anxious mice, respectively. Repeated administration of 5 mg/kg DM506 did not cause cumulative anxiolytic activity or any side effects. Higher doses of DM506 (40 mg/kg) induced sedative-like activity, which was inhibited by a selective 5-HT2A receptor antagonist, volinanserin. Electroencephalography results showed that 15 mg/kg DM506 fumarate increased the transition from a highly alert state (fast γ wavelength) to a more synchronized deep-sleeping activity (δ wavelength), which is reflected in the sedative/anxiolytic activity in mice but without the head-twitch response observed in hallucinogens. The functional, radioligand binding, and molecular docking results showed that DM506 binds to the agonist sites of human 5-HT2A (Ki = 24 nM) and 5-HT2B (Ki = 16 nM) receptors and activates them with a potency (EC50) of 9 nM and 3 nM, respectively. DM506 was relatively less potent and behaved as a partial agonist (efficacy <80%) for both receptor subtypes compared to the full agonist DOI (2,5-dimethoxy-4-iodoamphetamine). Our study showed for the first time that the non-hallucinogenic compound DM506 induces anxiolytic- and sedative-like activities in naïve and stressed/anxious mice in a dose-, time-, and volinanserin-sensitive manner, likely through mechanisms involving 5-HT2A receptor activation.

Probable digenic inheritance of Diamond-Blackfan anemia

A 26-year-old female proband with a clinical diagnosis and consistent phenotype of Diamond-Blackfan anemia (DBA, OMIM 105650) without an identified genotype was referred to the Undiagnosed Diseases Network. DBA is classically associated with monoallelic variants that have an autosomal-dominant or -recessive mode of inheritance. Intriguingly, her case was solved by a detection of a digenic interaction between non-allelic RPS19 and RPL27 variants. This was confirmed with a machine learning structural model, co-segregation analysis, and RNA sequencing. This is the first report of DBA caused by a digenic effect of two non-allelic variants demonstrated by machine learning structural model. This case suggests that atypical phenotypic presentations of DBA may be caused by digenic inheritance in some individuals. We also conclude that a machine learning structural model can be useful in detecting digenic models of possible interactions between products encoded by alleles of different genes inherited from non-affected carrier parents that can result in DBA with an unrealized 25% recurrence risk.

Template-free prediction of a new monotopic membrane protein fold and assembly by AlphaFold2

AlphaFold2 (AF2) has revolutionized the field of protein structural prediction. Here, we test its ability to predict the tertiary and quaternary structure of a previously undescribed scaffold with new folds and unusual architecture, the monotopic membrane protein caveolin-1 (CAV1). CAV1 assembles into a disc-shaped oligomer composed of 11 symmetrically arranged protomers, each assuming an identical new fold, and contains the largest parallel β-barrel known to exist in nature. Remarkably, AF2 predicts both the fold of the protomers and the interfaces between them. It also assembles between seven and 15 copies of CAV1 into disc-shaped complexes. However, the predicted multimers are energetically strained, especially the parallel β-barrel. These findings highlight the ability of AF2 to correctly predict new protein folds and oligomeric assemblies at a granular level while missing some elements of higher-order complexes, thus positing a new direction for the continued development of deep-learning protein structure prediction approaches.

Structural analysis of the P132L disease mutation in caveolin-1 reveals its role in the assembly of oligomeric complexes

Caveolin-1 (CAV1) is a membrane-sculpting protein that oligomerizes to generate flask-shaped invaginations of the plasma membrane known as caveolae. Mutations in CAV1 have been linked to multiple diseases in humans. Such mutations often interfere with oligomerization and the intracellular trafficking processes required for successful caveolae assembly, but the molecular mechanisms underlying these defects have not been structurally explained. Here, we investigate how a disease-associated mutation in one of the most highly conserved residues in CAV1, P132L, affects CAV1 structure and oligomerization. We show that P132 is positioned at a major site of protomer-protomer interactions within the CAV1 complex, providing a structural explanation for why the mutant protein fails to homo-oligomerize correctly. Using a combination of computational, structural, biochemical, and cell biological approaches, we find that despite its homo-oligomerization defects P132L is capable of forming mixed hetero-oligomeric complexes with WT CAV1 and that these complexes can be incorporated into caveolae. These findings provide insights into the fundamental mechanisms that control the formation of homo- and hetero-oligomers of caveolins that are essential for caveolae biogenesis, as well as how these processes are disrupted in human disease.

Silent agonists for α7 nicotinic acetylcholine receptors

We discuss models for the activation and desensitization of α7 nicotinic acetylcholine receptors (nAChRs) and the effects of efficacious type II positive allosteric modulators (PAMs) that destabilize α7 desensitized states. Type II PAMs such as PNU-120596 can be used to distinguish inactive compounds from silent agonists, compounds that produce little or no channel activation but stabilize the non-conducting conformations associated with desensitization. We discuss the effects of α7 nAChRs in cells of the immune system and their roles in modulating inflammation and pain through what has come to be known as the cholinergic anti-inflammatory system (CAS). Cells controlling CAS do not generate ion channel currents but rather respond to α7 drugs by modulating intracellular signaling pathways analogous to the effects of metabotropic receptors. Metabotropic signaling by α7 receptors appears to be mediated by receptors in nonconducting conformations and can be accomplished by silent agonists. We discuss electrophysiological structure-activity relationships for α7 silent agonists and their use in cell-based and in vivo assays for CAS regulation. We discuss the strongly desensitizing partial agonist GTS-21 and its effectiveness in modulation of CAS. We also review the properties of the silent agonist NS6740, which is remarkably effective at maintaining α7 receptors in PAM-sensitive desensitized states. Most silent agonists bind to sites overlapping those for orthosteric agonists, but some appear to bind to allosteric sites. Finally, we discuss α9^⁎ nAChRs and their potential role in CAS, and ligands that will be useful in defining and distinguishing the specific roles of α7 and α9 in CAS.

Benchmarking AlphaFold2 on peptide structure prediction

Recent advancements in computational tools have allowed protein structure prediction with high accuracy. Computational prediction methods have been used for modeling many soluble and membrane proteins, but the performance of these methods in modeling peptide structures has not yet been systematically investigated. We benchmarked the accuracy of AlphaFold2 in predicting 588 peptide structures between 10 and 40 amino acids using experimentally determined NMR structures as reference. Our results showed AlphaFold2 predicts α-helical, β-hairpin, and disulfide-rich peptides with high accuracy. AlphaFold2 performed at least as well if not better than alternative methods developed specifically for peptide structure prediction. AlphaFold2 showed several shortcomings in predicting Φ/Ψ angles, disulfide bond patterns, and the lowest RMSD structures failed to correlate with lowest pLDDT ranked structures. In summary, computation can be a powerful tool to predict peptide structures, but additional steps may be necessary to analyze and validate the results.

Molecular architecture of the human caveolin-1 complex

Membrane-sculpting proteins shape the morphology of cell membranes and facilitate remodeling in response to physiological and environmental cues. Complexes of the monotopic membrane protein caveolin function as essential curvature-generating components of caveolae, flask-shaped invaginations that sense and respond to plasma membrane tension. However, the structural basis for caveolin's membrane remodeling activity is currently unknown. Here, we show that, using cryo-electron microscopy, the human caveolin-1 complex is composed of 11 protomers organized into a tightly packed disc with a flat membrane-embedded surface. The structural insights suggest a previously unrecognized mechanism for how membrane-sculpting proteins interact with membranes and reveal how key regions of caveolin-1, including its scaffolding, oligomerization, and intramembrane domains, contribute to its function.

Veratridine Can Bind to a Site at the Mouth of the Channel Pore at Human Cardiac Sodium Channel NaV1.5

The cardiac sodium ion channel (Na_V1.5) is a protein with four domains (DI-DIV), each with six transmembrane segments. Its opening and subsequent inactivation results in the brief rapid influx of Na⁺ ions resulting in the depolarization of cardiomyocytes. The neurotoxin veratridine (VTD) inhibits Na_V1.5 inactivation resulting in longer channel opening times, and potentially fatal action potential prolongation. VTD is predicted to bind at the channel pore, but alternative binding sites have not been ruled out. To determine the binding site of VTD on Na_V1.5, we perform docking calculations and high-throughput electrophysiology experiments in the present study. The docking calculations identified two distinct binding regions. The first site was in the pore, close to the binding site of Na_V1.4 and Na_V1.5 blocking drugs in experimental structures. The second site was at the “mouth” of the pore at the cytosolic side, partly solvent-exposed. Mutations at this site (L409, E417, and I1466) had large effects on VTD binding, while residues deeper in the pore had no effect, consistent with VTD binding at the mouth site. Overall, our results suggest a VTD binding site close to the cytoplasmic mouth of the channel pore. Binding at this alternative site might indicate an allosteric inactivation mechanism for VTD at Na_V1.5.

The Allosteric Activation of α7 nAChR by α-Conotoxin MrIC Is Modified by Mutations at the Vestibular Site

α-conotoxins are 13–19 amino acid toxin peptides that bind various nicotinic acetylcholine receptor (nAChR) subtypes. α-conotoxin Mr1.7c (MrIC) is a 17 amino acid peptide that targets α7 nAChR. Although MrIC has no activating effect on α7 nAChR when applied by itself, it evokes a large response when co-applied with the type II positive allosteric modulator PNU-120596, which potentiates the α7 nAChR response by recovering it from a desensitized state. A lack of standalone activity, despite activation upon co-application with a positive allosteric modulator, was previously observed for molecules that bind to an extracellular domain allosteric activation (AA) site at the vestibule of the receptor. We hypothesized that MrIC may activate α7 nAChR allosterically through this site. We ran voltage-clamp electrophysiology experiments and in silico peptide docking calculations in order to gather evidence in support of α7 nAChR activation by MrIC through the AA site. The experiments with the wild-type α7 nAChR supported an allosteric mode of action, which was confirmed by the significantly increased MrIC + PNU-120596 responses of three α7 nAChR AA site mutants that were designed in silico to improve MrIC binding. Overall, our results shed light on the allosteric activation of α7 nAChR by MrIC and suggest the involvement of the AA site.

Prediction of amphipathic helix-membrane interactions with Rosetta

Amphipathic helices have hydrophobic and hydrophilic/charged residues situated on opposite faces of the helix. They can anchor peripheral membrane proteins to the membrane, be attached to integral membrane proteins, or exist as independent peptides. Despite the widespread presence of membrane-interacting amphipathic helices, there is no computational tool within Rosetta to model their interactions with membranes. In order to address this need, we developed the AmphiScan protocol with PyRosetta, which runs a grid search to find the most favorable position of an amphipathic helix with respect to the membrane. The performance of the algorithm was tested in benchmarks with the RosettaMembrane, ref2015_memb, and franklin2019 score functions on six engineered and 44 naturally-occurring amphipathic helices using membrane coordinates from the OPM and PDBTM databases, OREMPRO server, and MD simulations for comparison. The AmphiScan protocol predicted the coordinates of amphipathic helices within less than 3Å of the reference structures and identified membrane-embedded residues with a Matthews Correlation Constant (MCC) of up to 0.57. Overall, AmphiScan stands as fast, accurate, and highly-customizable protocol that can be pipelined with other Rosetta and Python applications.

Amphipathic helices are important targets as antibacterial peptides and as domains of membrane proteins that play a role in sensing the membrane environment. Understanding how amphipathic helices interact with membrane enables us to design better peptides and understand how membrane proteins use them to interact with their environment. However, there is a limited number of tools available for the modeling of amphipathic helices in membranes. Implicit membrane models can be used for this purpose as simplistic representations of the membrane environment. In this work, we developed the AmphiScan protocol that can be used to predict membrane coordinates of amphipathic helices starting with a helix structure in an implicit membrane environment. We benchmarked the performance of AmphiScan on engineered LK peptides, naturally-occurring amphipathic helices, and hydrophobic and hydrophilic peptides. Our approach provides a reliable and customizable tool to model amphipathic helix–membrane interactions, and pose a platform for the screening of amphipathic helix properties in silico.

Modeling Immunity with Rosetta: Methods for Antibody and Antigen Design

Structure-based antibody and antigen design has advanced greatly in recent years, due not only to the increasing availability of experimentally determined structures but also to improved computational methods for both prediction and design. Constant improvements in performance within the Rosetta software suite for biomolecular modeling have given rise to a greater breadth of structure prediction, including docking and design application cases for antibody and antigen modeling. Here, we present an overview of current protocols for antibody and antigen modeling using Rosetta and exemplify those by detailed tutorials originally developed for a Rosetta workshop at Vanderbilt University. These tutorials cover antibody structure prediction, docking, and design and antigen design strategies, including the addition of glycans in Rosetta. We expect that these materials will allow novice users to apply Rosetta in their own projects for modeling antibodies and antigens.

A Computational Analysis of the Factors Governing the Dynamics of α7 nAChR and Its Homologs

The α7 nicotinic acetylcholine receptor is a homopentameric ion channel from the Cys-loop receptor superfamily targeted for psychiatric indications and inflammatory pain. Molecular dynamics studies of the receptor have focused on residue mobility and global conformational changes to address receptor function. However, a comparative analysis of α7 with its homologs that cannot trigger channel opening has not been made so far. To identify the residues involved in α7 activation, we ran triplicate 500-ns molecular dynamics simulations with an α7 extracellular domain homology model and two acetylcholine-binding protein homologs. We tested the effect of ligand binding and amino acid sequence on the structure and dynamics of the three proteins. We found that mobile regions identified based on root mean-square deviation and root mean-square fluctuation values are not always consistent among the individual α7 extracellular domain simulations. Comparison of the replica-average properties of the three proteins based on dynamic cross-correlation maps showed that ligand binding affects the coupling between the C-loop and the Cys-loop, vestibular loop, and β1-β2 loops. In addition, the main-immunogenic-region-like domain of α7 went through correlated motions with multiple domains of the receptor. These correlated motions were absent or diminished in α7 homologs, suggesting a unique role in α7 activation.

An Investigation of Three-Finger Toxin-nAChR Interactions through Rosetta Protein Docking

Three-finger toxins (3FTX) are a group of peptides that affect multiple receptor types. One group of proteins affected by 3FTX are nicotinic acetylcholine receptors (nAChR). Structural information on how neurotoxins interact with nAChR is limited and is confined to a small group of neurotoxins. Therefore, in silico methods are valuable in understanding the interactions between 3FTX and different nAChR subtypes, but there are no established protocols to model 3FTX-nAChR interactions. We followed a homology modeling and protein docking protocol to address this issue and tested its success on three different systems. First, neurotoxin peptides co-crystallized with acetylcholine binding protein (AChBP) were re-docked to assess whether Rosetta protein-protein docking can reproduce the native poses. Second, experimental data on peptide binding to AChBP was used to test whether the docking protocol can qualitatively distinguish AChBP-binders from non-binders. Finally, we docked eight peptides with known α7 and muscle-type nAChR binding properties to test whether the protocol can explain the differential activities of the peptides at the two receptor subtypes. Overall, the docking protocol predicted the qualitative and some specific aspects of 3FTX binding to nAChR with reasonable success and shed light on unknown aspects of 3FTX binding to different receptor subtypes.

Nicotinic receptor pharmacology in silico: Insights and challenges

Nicotinic acetylcholine receptors (nAChR) are homo- or hetero-pentameric ligand-gated ion channels of the Cys-loop superfamily and play important roles in the nervous system and muscles. Studies on nAChR benefit from in silico modeling due to the lack of high-resolution structures for most receptor subtypes and challenges in experiments addressing the complex mechanism of activation involving allosteric sites. Although there is myriad of computational modeling studies on nAChR, the multitude of the methods and parameters used in these studies makes modeling nAChR a daunting task, particularly for the non-experts in the field. To address this problem, the modeling literature on Torpedo nAChR and α7 nAChR were focused on as examples of heteromeric and homomeric nAChR, and the key in silico modeling studies between the years 1995-2019 were concisely reviewed. This was followed by a critical analysis of these studies by comparing the findings with each other and with the emerging experimental and computational data on nAChR. Based on these critical analyses, suggestions were made to guide the future researchers in the field of in silico modeling of nAChR. This article is part of the special issue on 'Contemporary Advances in Nicotine Neuropharmacology'.

In Silico Modeling of the α7 Nicotinic Acetylcholine Receptor: New Pharmacological Challenges Associated with Multiple Modes of Signaling

The α7 nicotinic acetylcholine receptor is a homopentameric ion-channel of the Cys-loop superfamily characterized by its low probability of opening, high calcium permeability, and rapid desensitization. The α7 receptor has been targeted for the treatment of the cognitive symptoms of schizophrenia, depression, and Alzheimer's disease, but it is also involved in inflammatory modulation as a part of the cholinergic anti-inflammatory pathway. Despite its functional importance, in silico studies of the α7 receptor cannot produce a general model explaining the structural features of receptor activation, nor predict the mode of action for various ligand classes. Two particular problems in modeling the α7 nAChR are the absence of a high-resolution structure and the presence of five potentially nonequivalent orthosteric ligand binding sites. There is wide variability regarding the templates used for homology modeling, types of ligands investigated, simulation methods, and simulation times. However, a systematic survey focusing on the methodological similarities and differences in modeling α7 has not been done. In this work, we make a critical analysis of the modeling literature of α7 nAChR by comparing the findings of computational studies with each other and with experimental studies under the main topics of structural studies, ligand binding studies, and comparisons with other nAChR. In light of our findings, we also summarize current problems in the field and make suggestions for future studies concerning modeling of the α7 receptor.

Allosteric Agonism of α7 Nicotinic Acetylcholine Receptors: Receptor Modulation Outside the Orthosteric Site

Nicotinic acetylcholine receptors (nAChRs) are members of the Cys-loop superfamily of ligand-gated ion channels. Typically, channel activation follows the binding of agonists to the orthosteric binding sites of the receptor. α7 nAChRs have a very low probability of channel activation, which can be reversed by the binding of α7 selective positive allosteric modulators (PAMs) to putative sites within the transmembrane domains. Although typical PAMs, like PNU-120596, require coapplication of an orthosteric agonist to produce large channel activations, some, like GAT107 and B-973B [(S)-3-(3,4-difluorophenyl)-N-(1-(6-(4-(pyridin-2-yl)piperazin-1-yl)pyrazin-2-yl)ethyl)propanamide], are characterized as allosteric activating PAMs, which also bind to an allosteric activation (AA) site in the extracellular domain and activate the α7 ion channel by themselves. We had previously characterized N,N-diethyl-N'-phenylpiperazine analogs with various functions. In this work, we docked members of this family to a homology model of the α7 receptor extracellular domain. The compound 1,1-diethyl-4(naphthalene-2-yl)piperazin-1-ium (2NDEP) a weak partial agonist, showed particularly favorable docking and binding energies at the putative AA site of the receptor. We hypothesized that 2NDEP could couple with PAMs through the AA site. This hypothesis was tested with the α7 mutant C190A, which is not activated by orthosteric agonists but is effectively activated by GAT107. The results showed that 2NDEP acts as an allosteric agonist of α7C190A when coapplied with the PAM PNU-120596. Also, the allosteric activity was nearly abolished upon coapplication with the AA site-selective antagonist 2,3,5,6MP-TQS (cis-trans-4-(2,3,5,6-tetramethylphenyl)-3a,4,5,9b-tetrahydro-3H-cyclopenta[c]quinoline-8-sulfonamide), consistent with AA site involvement. Overall, our findings show a novel mode of agonism through an allosteric site in the extracellular domain of α7 nAChR.

Macroscopic and Microscopic Activation of α7 Nicotinic Acetylcholine Receptors by the Structurally Unrelated Allosteric Agonist-Positive Allosteric Modulators (ago-PAMs) B-973B and GAT107

B-973 is an efficacious type II positive allosteric modulator (PAM) of α7 nicotinic acetylcholine receptors that, like 4BP-TQS and its active isomer GAT107, can produce direct allosteric activation in addition to potentiation of orthosteric agonist activity, which identifies it as an allosteric activating (ago)-PAM. We compared the properties of B-973B, the active enantiomer of B-973, with those of GAT107 regarding the separation of allosteric potentiation and activation. Both ago-PAMs can strongly activate mutants of α7 that are insensitive to standard orthosteric agonists like acetylcholine. Likewise, the activity of both ago-PAMs is largely eliminated by the M254L mutation in the putative transmembrane PAM-binding site. Allosteric activation by B-973B appeared more protracted than that produced by GAT107, and B-973B responses were relatively insensitive to the noncompetitive antagonist mecamylamine compared with GAT107 responses. Similar differences are also seen in the single-channel currents. The two agents generate unique profiles of full-conductance and subconductance states, with B-973B producing protracted bursts, even in the presence of mecamylamine. Modeling and docking studies suggest that the molecular basis for these effects depends on specific interactions in both the extracellular and transmembrane domains of the receptor.

Sulfonium as a Surrogate for Ammonium: A New α7 Nicotinic Acetylcholine Receptor Partial Agonist with Desensitizing Activity

Weak partial agonists that promote a desensitized state of the α7 nicotinic acetylcholine receptor (nAChR) have been associated with anti-inflammatory effects. Exemplar compounds feature a tertiary or quaternary ammonium group. We report the synthesis, structure, and electrophysiological evaluation of 1-ethyl-4-phenylthiomorpholin-1-ium triflate, a weak partial agonist with a sulfonium isostere of the ammonium pharmacophore. These results offer new insights in understanding nAChR-ligand interactions and provide a new chemical space to target the α7 nAChR.