Evolutionary-Conserved Allosteric Properties of Three Neuronal Calcium Sensor Proteins

Disrupted protein interaction dynamics in a genetic neurodevelopmental disorder revealed by structural bioinformatics and genetic code expansion

Deciphering the structural effects of gene variants is essential for understanding the pathophysiological mechanisms of genetic diseases. Using a neurodevelopmental disorder called Bosch-Boonstra-Schaaf Optic Atrophy Syndrome (BBSOAS) as a genetic disease model, we applied structural bioinformatics and Genetic Code Expansion (GCE) strategies to assess the pathogenic impact of human NR2F1 variants and their binding with known and novel partners. While the computational analyses of the NR2F1 structure delineated the molecular basis of the impact of several variants on the isolated and complexed structures, the GCE enabled covalent and site-specific capture of transient supramolecular interactions in living cells. This revealed the variable quaternary conformations of NR2F1 variants and highlighted the disrupted interplay with dimeric partners and the newly identified co-factor, CRABP2. The disclosed consequence of the pathogenic mutations on the conformation, supramolecular interplay, and alterations in the cell cycle, viability, and sub-cellular localization of the different variants reflect the heterogeneous disease spectrum of BBSOAS and set up novel foundation for unveiling the complexity of neurodevelopmental diseases.

PyInteraph2 and PyInKnife2 to Analyze Networks in Protein Structural Ensembles

Due to the complex nature of noncovalent interactions and their long-range effects, analyzing protein conformations using network theory can be enlightening. Protein Structure Networks (PSNs) provide a convenient formalism to study protein structures in relation to essential properties such as key residues for structural stability, allosteric communication, and the effects of modifications of the protein. PSNs can be defined according to very different principles, and the available tools have limitations in input formats, supported models, and version control. Other outstanding problems are related to the definition of network cutoffs and the assessment of the stability of the network properties. The protein science community could benefit from a common framework to carry out these analyses and make them easier to reproduce, reuse, and evaluate. We here provide two open-source software packages, PyInteraph2 and PyInKnife2, to implement and analyze PSNs in a reproducible and documented manner. PyInteraph2 interfaces with multiple formats for protein ensembles and incorporates different network models with the possibility of integrating them into a macronetwork and performing various downstream analyses, including hubs, connected components, and several other centrality measures, and visualizes the networks or further analyzes them thanks to compatibility with Cytoscape.PyInKnife2 that supports the network models implemented in PyInteraph2. It employs a jackknife resampling approach to estimate the convergence of network properties and streamline the selection of distance cutoffs. We foresee that the modular structure of the code and the supported version control system will promote the transition to a community-driven effort, boost reproducibility, and establish common protocols in the PSN field. As developers, we will guarantee the introduction of new functionalities and maintenance, assistance, and training of new contributors.

Calmodulin variants associated with congenital arrhythmia impair selectivity for ryanodine receptors

Among its many molecular targets, the ubiquitous calcium sensor protein calmodulin (CaM) recognizes and regulates the activity of ryanodine receptors type 1 (RyR1) and 2 (RyR2), mainly expressed in skeletal and cardiac muscle, respectively. Such regulation is essential to achieve controlled contraction of muscle cells. To unravel the molecular mechanisms underlying the target recognition process, we conducted a comprehensive biophysical investigation of the interaction between two calmodulin variants associated with congenital arrhythmia, namely N97I and Q135P, and a highly conserved calmodulin-binding region in RyR1 and RyR2. The structural, thermodynamic, and kinetic properties of protein-peptide interactions were assessed together with an in-depth structural and topological investigation based on molecular dynamics simulations. This integrated approach allowed us to identify amino acids that are crucial in mediating allosteric processes, which enable high selectivity in molecular target recognition. Our results suggest that the ability of calmodulin to discriminate between RyR1 an RyR2 targets depends on kinetic discrimination and robust allosteric communication between Ca²⁺-binding sites (EF1-EF3 and EF3-EF4 pairs), which is perturbed in both N97I and Q135P arrhythmia-associated variants.

Nanomolar affinity of EF-hands in neuronal calcium sensor 1 for bivalent cations Pb2+, Mn2+ and Hg2

Abiogenic metals Pb and Hg are highly toxic since chronic and/or acute exposure often leads to severe neuropathologies. Mn2+ is an essential metal ion but in excess can impair neuronal function. In this study, we address in vitro the interactions between neuronal calcium sensor 1 (NCS1) and divalent cations. Results showed that non-physiological ions (Pb2+, Mn2+ and Hg2+) bind to EF-hands in NCS1 with nanomolar affinity and lower equilibrium dissociation constant than the physiological Ca2+ ion. (Kd,Pb2+ = 7.0±1.0 nM; Kd,Mn2+ = 34.0±6.0 nM; Kd, Hg2+ = 0.5±0.1 nM and 27.0±13.0 nM and Kd,Ca2+ = 96.0±48.0 nM). Native ultra-high resolution mass spectrometry (FT-ICR MS) and trapped ion mobility spectrometry - mass spectrometry (nESI-TIMS-MS) studies provided the NCS1-metal complex compositions - up to four Ca2+ or Mn2+ ions and three Pb2+ ions (M⋅Pb1-3Ca1-3, M⋅Mn1-4Ca1-2, and M⋅Ca1-4) were observed in complex - and similarity across the mobility profiles suggests that the overall native structure is preserved regardless of the number and type of cations. However, the non-physiological metal ions (Pb2+, Mn2+, and Hg2+) binding to NCS1 leads to more efficient quenching of Trp emission and a decrease in W30 and W103 solvent exposure compared to the apo and Ca2+ bound form, although the secondary structural rearrangement and exposure of hydrophobic sites are analogous to those for Ca2+ bound protein. Only Pb2+ and Hg2+ binding to EF-hands leads to the NCS1 dimerization whereas Mn2+ bound NCS1 remains in the monomeric form, suggesting that other factors in addition to metal ion coordination, are required for protein dimerization.

Molecular Properties of Human Guanylate Cyclase-Activating Protein 3 (GCAP3) and Its Possible Association with Retinitis Pigmentosa

The cone-specific guanylate cyclase-activating protein 3 (GCAP3), encoded by the GUCA1C gene, has been shown to regulate the enzymatic activity of membrane-bound guanylate cyclases (GCs) in bovine and teleost fish photoreceptors, to an extent comparable to that of the paralog protein GCAP1. To date, the molecular mechanisms underlying GCAP3 function remain largely unexplored. In this work, we report a thorough characterization of the biochemical and biophysical properties of human GCAP3, moreover, we identified an isolated case of retinitis pigmentosa, in which a patient carried the c.301G>C mutation in GUCA1C, resulting in the substitution of a highly conserved aspartate residue by a histidine (p.(D101H)). We found that myristoylated GCAP3 can activate GC1 with a similar Ca2+-dependent profile, but significantly less efficiently than GCAP1. The non-myristoylated form did not induce appreciable regulation of GC1, nor did the p.D101H variant. GCAP3 forms dimers under physiological conditions, but at odds with its paralogs, it tends to form temperature-dependent aggregates driven by hydrophobic interactions. The peculiar properties of GCAP3 were confirmed by 2 ms molecular dynamics simulations, which for the p.D101H variant highlighted a very high structural flexibility and a clear tendency to lose the binding of a Ca2+ ion to EF3. Overall, our data show that GCAP3 has unusual biochemical properties, which make the protein significantly different from GCAP1 and GCAP2. Moreover, the newly identified point mutation resulting in a substantially unfunctional protein could trigger retinitis pigmentosa through a currently unknown mechanism.

Alterations in calmodulin-cardiac ryanodine receptor molecular recognition in congenital arrhythmias

Calmodulin (CaM), a ubiquitous and highly conserved Ca²⁺-sensor protein involved in the regulation of over 300 molecular targets, has been recently associated with severe forms of lethal arrhythmia. Here, we investigated how arrhythmia-associated mutations in CaM localized at the C-terminal lobe alter the molecular recognition with Ryanodine receptor 2 (RyR2), specifically expressed in cardiomyocytes. We performed an extensive structural, thermodynamic, and kinetic characterization of the variants D95V/H in the EF3 Ca²⁺-binding motif and of the D129V and D131H/E variants in the EF4 motif, and probed their interaction with RyR2. Our results show that the specific structural changes observed for individual CaM variants do not extend to the complex with the RyR2 target. Indeed, some common alterations emerge at the protein–protein interaction level, suggesting the existence of general features shared by the arrhythmia-associated variants. All mutants showed a faster rate of dissociation from the target peptide than wild-type CaM. Integration of spectroscopic data with exhaustive molecular dynamics simulations suggests that, in the presence of Ca²⁺, functional recognition involves allosteric interactions initiated by the N-terminal lobe of CaM, which shows a lower affinity for Ca²⁺ compared to the C-terminal lobe in the isolated protein.

Protein conformational switch discerned via network centrality properties

Network analysis has emerged as a powerful tool for examining structural biology systems. The spatial organization of the components of a biomolecular structure has been rendered as a graph representation and analyses have been performed to deduce the biophysical and mechanistic properties of these components. For proteins, the analysis of protein structure networks (PSNs), especially via network centrality measurements and cluster coefficients, has led to identifying amino acid residues that play key functional roles and classifying amino acid residues in general. Whether these network properties examined in various studies are sensitive to subtle (yet biologically significant) conformational changes remained to be addressed. Here, we focused on four types of network centrality properties (betweenness, closeness, degree, and eigenvector centralities) for conformational changes upon ligand binding of a sensor protein (constitutive androstane receptor) and an allosteric enzyme (ribonucleotide reductase). We found that eigenvector centrality is sensitive and can distinguish salient structural features between protein conformational states while other centrality measures, especially closeness centrality, are less sensitive and rather generic with respect to the structural specificity. We also demonstrated that an ensemble-informed, modified PSN with static edges removed (which we term PSN*) has enhanced sensitivity at discerning structural changes.

Expanding the Clinical and Genetic Spectrum of RAB28-Related Cone-Rod Dystrophy: Pathogenicity of Novel Variants in Italian Families

The small Ras-related GTPase Rab-28 is highly expressed in photoreceptor cells, where it possibly participates in membrane trafficking. To date, six alterations in the RAB28 gene have been associated with autosomal recessive cone-rod dystrophies. Confirmed variants include splicing variants, missense and nonsense mutations. Here, we present a thorough phenotypical and genotypical characterization of five individuals belonging to four Italian families, constituting the largest cohort of RAB28 patients reported in literature to date. All probands displayed similar clinical phenotype consisting of photophobia, decreased visual acuity, central outer retinal thinning, and impaired color vision. By sequencing the four probands, we identified: a novel homozygous splicing variant; two novel nonsense variants in homozygosis; a novel missense variant in compound heterozygous state with a previously reported nonsense variant. Exhaustive molecular dynamics simulations of the missense variant p.(Thr26Asn) in both its active and inactive states revealed an allosteric structural mechanism that impairs the binding of Mg²⁺, thus decreasing the affinity for GTP. The impaired GTP-GDP exchange ultimately locks Rab-28 in a GDP-bound inactive state. The loss-of-function mutation p.(Thr26Asn) was present in a compound heterozygosis with the nonsense variant p.(Arg137*), which does not cause mRNA-mediated decay, but is rather likely degraded due to its incomplete folding. The frameshift p.(Thr26Valfs4*) and nonsense p.(Leu13*) and p.(Trp107*) variants, if translated, would lack several key structural components necessary for the correct functioning of the encoded protein.

Modulation of Guanylate Cyclase Activating Protein 1 (GCAP1) Dimeric Assembly by Ca2+ or Mg2+: Hints to Understand Protein Activity

The guanylyl cyclase-activating protein 1, GCAP1, activates or inhibits retinal guanylyl cyclase (retGC) depending on cellular Ca²⁺ concentrations. Several point mutations of GCAP1 have been associated with impaired calcium sensitivity that eventually triggers progressive retinal degeneration. In this work, we demonstrate that the recombinant human protein presents a highly dynamic monomer-dimer equilibrium, whose dissociation constant is influenced by salt concentration and, more importantly, by protein binding to Ca²⁺ or Mg²⁺. Based on small-angle X-ray scattering data, protein-protein docking, and molecular dynamics simulations we propose two novel three-dimensional models of Ca²⁺-bound GCAP1 dimer. The different propensity of human GCAP1 to dimerize suggests structural differences induced by cation binding potentially involved in the regulation of retGC activity.

Disorder in a two-domain neuronal Ca2+-binding protein regulates domain stability and dynamics using ligand mimicry

Understanding the interplay between sequence, structure and function of proteins has been complicated in recent years by the discovery of intrinsically disordered proteins (IDPs), which perform biological functions in the absence of a well-defined three-dimensional fold. Disordered protein sequences account for roughly 30% of the human proteome and in many proteins, disordered and ordered domains coexist. However, few studies have assessed how either feature affects the properties of the other. In this study, we examine the role of a disordered tail in the overall properties of the two-domain, calcium-sensing protein neuronal calcium sensor 1 (NCS-1). We show that loss of just six of the 190 residues at the flexible C-terminus is sufficient to severely affect stability, dynamics, and folding behavior of both ordered domains. We identify specific hydrophobic contacts mediated by the disordered tail that may be responsible for stabilizing the distal N-terminal domain. Moreover, sequence analyses indicate the presence of an LSL-motif in the tail that acts as a mimic of native ligands critical to the observed order-disorder communication. Removing the disordered tail leads to a shorter life-time of the ligand-bound complex likely originating from the observed destabilization. This close relationship between order and disorder may have important implications for how investigations into mixed systems are designed and opens up a novel avenue of drug targeting exploiting this type of behavior.

Aldo-Keto Reductase 1C1 (AKR1C1) as the First Mutated Gene in a Family with Nonsyndromic Primary Lipedema

Lipedema is an often underdiagnosed chronic disorder that affects subcutaneous adipose tissue almost exclusively in women, which leads to disproportionate fat accumulation in the lower and upper body extremities. Common comorbidities include anxiety, depression, and pain. The correlation between mood disorder and subcutaneous fat deposition suggests the involvement of steroids metabolism and neurohormones signaling, however no clear association has been established so far. In this study, we report on a family with three patients affected by sex-limited autosomal dominant nonsyndromic lipedema. They had been screened by whole exome sequencing (WES) which led to the discovery of a missense variant p.(Leu213Gln) in AKR1C1, the gene encoding for an aldo-keto reductase catalyzing the reduction of progesterone to its inactive form, 20-α-hydroxyprogesterone. Comparative molecular dynamics simulations of the wild-type vs. variant enzyme, corroborated by a thorough structural and functional bioinformatic analysis, suggest a partial loss-of-function of the variant. This would result in a slower and less efficient reduction of progesterone to hydroxyprogesterone and an increased subcutaneous fat deposition in variant carriers. Overall, our results suggest that AKR1C1 is the first candidate gene associated with nonsyndromic lipedema.

Missense mutations affecting Ca2+-coordination in GCAP1 lead to cone-rod dystrophies by altering protein structural and functional properties

Guanylate cyclase activating protein 1 (GCAP1) is a neuronal calcium sensor (NCS) involved in the early biochemical steps underlying the phototransduction cascade. By switching from a Ca²⁺-bound form in the dark to a Mg²⁺-bound state following light activation of the cascade, GCAP1 triggers the activation of the retinal guanylate cyclase (GC), thus replenishing the levels of 3',5'-cyclic monophosphate (cGMP) necessary to re-open CNG channels. Here, we investigated the structural and functional effects of three missense mutations in GCAP1 associated with cone-rod dystrophy, which severely perturb the homeostasis of cGMP and Ca²⁺. Substitutions affect residues directly involved in Ca²⁺ coordination in either EF3 (D100G) or EF4 (E155A and E155G) Ca²⁺ binding motifs. We found that all GCAP1 variants form relatively stable dimers showing decreased apparent affinity for Ca²⁺ and blocking the enzyme in a constitutively active state at physiological levels of Ca²⁺. Interestingly, by corroborating spectroscopic experiments with molecular dynamics simulations we show that beside local structural effects, mutation of the bidentate glutamate in an EF-hand calcium binding motif can profoundly perturb the flexibility of the adjacent EF-hand as well, ultimately destabilizing the whole domain. Therefore, while Ca²⁺-binding to GCAP1 per se occurs sequentially, allosteric effects may connect EF hand motifs, which appear to be essential for the integrity of the structural switch mechanism in GCAP1, and perhaps in other NCS proteins.

New age-related hearing loss candidate genes in humans: an ongoing challenge

Age-related hearing loss (ARHL) is the most frequent sensory disorder in the elderly, affecting approximately one-third of people aged more than 65 years. Despite a large number of people affected, ARHL is still an area of unmet clinical needs, and only a few ARHL susceptibility genes have been detected so far. In order to further investigate the genetics of ARHL, we analyzed a series of 46 ARHL candidate genes, selected according to previous Genome Wide Association Studies (GWAS) data, literature updates and animal models, in a large cohort of 464 Italian ARHL patients. We have filtered the variants according to a) pathogenicity prediction, b) allele frequency in public databases, c) allele frequency in an internal cohort of 113 healthy matched controls, and 81 healthy semi-supercentenarians. After data analysis, all the variants of interest have been tested by functional "in silico" or "in vitro" experiments (i.e., molecular dynamics simulations and protein translation analysis) to assess their pathogenic role, and the expression of the mutated genes have been checked in mouse or zebrafish inner ear. This multi-step approach led to the characterization of a series of ultra-rare likely pathogenic variants in DCLK1, SLC28A3, CEP104, and PCDH20 genes, contributing to describe the first association of these genes with ARHL in humans. These results provide essential insights on the understanding of the molecular bases of such a complex, heterogeneous and frequent disorder, unveiling new possible targets for the future development of innovative therapeutic and preventive approaches that could improve the quality of life of the millions of people affected worldwide.

Constitutive Activation of Guanylate Cyclase by the G86R GCAP1 Variant Is Due to "Locking" Cation-π Interactions that Impair the Activator-to-Inhibitor Structural Transition

Guanylate Cyclase activating protein 1 (GCAP1) mediates the Ca²⁺-dependent regulation of the retinal Guanylate Cyclase (GC) in photoreceptors, acting as a target inhibitor at high [Ca²⁺] and as an activator at low [Ca²⁺]. Recently, a novel missense mutation (G86R) was found in GUCA1A, the gene encoding for GCAP1, in patients diagnosed with cone-rod dystrophy. The G86R substitution was found to affect the flexibility of the hinge region connecting the N- and C-domains of GCAP1, resulting in decreased Ca^2+-sensitivity and abnormally enhanced affinity for GC. Based on a structural model of GCAP1, here, we tested the hypothesis of a cation-π interaction between the positively charged R86 and the aromatic W94 as the main mechanism underlying the impaired activator-to-inhibitor conformational change. W94 was mutated to F or L, thus, resulting in the double mutants G86R+W94L/F. The double mutants showed minor structural and stability changes with respect to the single G86R mutant, as well as lower affinity for both Mg²⁺ and Ca²⁺, moreover, substitutions of W94 abolished "phase II" in Ca²⁺-titrations followed by intrinsic fluorescence. Interestingly, the presence of an aromatic residue in position 94 significantly increased the aggregation propensity of Ca²⁺-loaded GCAP1 variants. Finally, atomistic simulations of all GCAP1 variants in the presence of Ca²⁺ supported the presence of two cation-π interactions involving R86, which was found to act as a bridge between W94 and W21, thus, locking the hinge region in an activator-like conformation and resulting in the constitutive activation of the target under physiological conditions.

Oligomeric state, hydrodynamic properties and target recognition of human Calcium and Integrin Binding protein 2 (CIB2)

Calcium- and Integrin-Binding protein 2 (CIB2) is a small and ubiquitously expressed protein with largely unknown biological function but ascertained role in hearing physiology and disease. Recent studies found that CIB2 binds Ca2+ with moderate affinity and dimerizes under conditions mimicking the physiological ones. Here we provided new lines of evidence on CIB2 oligomeric state and the mechanism of interaction with the α7B integrin target. Based on a combination of native mass spectrometry, chemical cross-linking/mass spectrometry, analytical gel filtration, dynamic light scattering and molecular dynamics simulations we conclude that CIB2 is monomeric under all tested conditions and presents uncommon hydrodynamic properties, most likely due to the high content of hydrophobic solvent accessible surface. Surface plasmon resonance shows that the interaction with α7B occurs with relatively low affinity and is limited to the cytosolic region proximal to the membrane, being kinetically favored in the presence of physiological Mg2+ and in the absence of Ca2+. Although CIB2 binds to an α7B peptide in a 1:1 stoichiometry, the formation of the complex might induce binding of another CIB2 molecule.

Effects of Membrane and Biological Target on the Structural and Allosteric Properties of Recoverin: A Computational Approach

Recoverin (Rec) is a prototypical calcium sensor protein primarily expressed in the vertebrate retina. The binding of two Ca²⁺ ions to the functional EF-hand motifs induces the extrusion of a myristoyl group that increases the affinity of Rec for the membrane and leads to the formation of a complex with rhodopsin kinase (GRK1). Here, unbiased all-atom molecular dynamics simulations were performed to monitor the spontaneous insertion of the myristoyl group into a model multicomponent biological membrane for both isolated Rec and for its complex with a peptide from the GRK1 target. It was found that the functional membrane anchoring of the myristoyl group is triggered by persistent electrostatic protein-membrane interactions. In particular, salt bridges between Arg43, Arg46 and polar heads of phosphatidylserine lipids are necessary to enhance the myristoyl hydrophobic packing in the Rec-GRK1 assembly. The long-distance communication between Ca²⁺-binding EF-hands and residues at the interface with GRK1 is significantly influenced by the presence of the membrane, which leads to dramatic changes in the connectivity of amino acids mediating the highest number of persistent interactions (hubs). In conclusion, specific membrane composition and allosteric interactions are both necessary for the correct assembly and dynamics of functional Rec-GRK1 complex.