Allosteric communication pathways routed by Ca2+/Mg2+ exchange in GCAP1 selectively switch target regulation modes

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.

Recombinant protein delivery enables modulation of the phototransduction cascade in mouse retina

Inherited retinal dystrophies are often associated with mutations in the genes involved in the phototransduction cascade in photoreceptors, a paradigmatic signaling pathway mediated by G protein-coupled receptors. Photoreceptor viability is strictly dependent on the levels of the second messengers cGMP and Ca2+. Here we explored the possibility of modulating the phototransduction cascade in mouse rods using direct or liposome-mediated administration of a recombinant protein crucial for regulating the interplay of the second messengers in photoreceptor outer segments. The effects of administration of the free and liposome-encapsulated human guanylate cyclase-activating protein 1 (GCAP1) were compared in biological systems of increasing complexity (in cyto, ex vivo, and in vivo). The analysis of protein biodistribution and the direct measurement of functional alteration in rod photoresponses show that the exogenous GCAP1 protein is fully incorporated into the mouse retina and photoreceptor outer segments. Furthermore, only in the presence of a point mutation associated with cone-rod dystrophy in humans p.(E111V), protein delivery induces a disease-like electrophysiological phenotype, consistent with constitutive activation of the retinal guanylate cyclase. Our study demonstrates that both direct and liposome-mediated protein delivery are powerful complementary tools for targeting signaling cascades in neuronal cells, which could be particularly important for the treatment of autosomal dominant genetic 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.

Binding and Functional Folding (BFF): A Physiological Framework for Studying Biomolecular Interactions and Allostery

EF-hand Ca2+-binding proteins (CBPs), such as S100 proteins (S100s) and calmodulin (CaM), are signaling proteins that undergo conformational changes upon increasing intracellular Ca2+. Upon binding Ca2+, S100 proteins and CaM interact with protein targets and induce important biological responses. The Ca2+-binding affinity of CaM and most S100s in the absence of target is weak (CaKD > 1 μM). However, upon effector protein binding, the Ca2+ affinity of these proteins increases via heterotropic allostery (CaKD < 1 μM). Because of the high number and micromolar concentrations of EF-hand CBPs in a cell, at any given time, allostery is required physiologically, allowing for (i) proper Ca2+ homeostasis and (ii) strict maintenance of Ca2+-signaling within a narrow dynamic range of free Ca2+ ion concentrations, [Ca2+]free. In this review, mechanisms of allostery are coalesced into an empirical "binding and functional folding (BFF)" physiological framework. At the molecular level, folding (F), binding and folding (BF), and BFF events include all atoms in the biomolecular complex under study. The BFF framework is introduced with two straightforward BFF types for proteins (type 1, concerted; type 2, stepwise) and considers how homologous and nonhomologous amino acid residues of CBPs and their effector protein(s) evolved to provide allosteric tightening of Ca2+ and simultaneously determine how specific and relatively promiscuous CBP-target complexes form as both are needed for proper cellular function.

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.

A Novel GUCA1A Variant Associated with Cone Dystrophy Alters cGMP Signaling in Photoreceptors by Strongly Interacting with and Hyperactivating Retinal Guanylate Cyclase

Guanylate cyclase-activating protein 1 (GCAP1), encoded by the GUCA1A gene, is a neuronal calcium sensor protein involved in shaping the photoresponse kinetics in cones and rods. GCAP1 accelerates or slows the cGMP synthesis operated by retinal guanylate cyclase (GC) based on the light-dependent levels of intracellular Ca²⁺, thereby ensuring a timely regulation of the phototransduction cascade. We found a novel variant of GUCA1A in a patient affected by autosomal dominant cone dystrophy (adCOD), leading to the Asn104His (N104H) amino acid substitution at the protein level. While biochemical analysis of the recombinant protein showed impaired Ca²⁺ sensitivity of the variant, structural properties investigated by circular dichroism and limited proteolysis excluded major structural rearrangements induced by the mutation. Analytical gel filtration profiles and dynamic light scattering were compatible with a dimeric protein both in the presence of Mg²⁺ alone and Mg²⁺ and Ca²⁺. Enzymatic assays showed that N104H-GCAP1 strongly interacts with the GC, with an affinity that doubles that of the WT. The doubled IC₅₀ value of the novel variant (520 nM for N104H vs. 260 nM for the WT) is compatible with a constitutive activity of GC at physiological levels of Ca²⁺. The structural region at the interface with the GC may acquire enhanced flexibility under high Ca²⁺ conditions, as suggested by 2 μs molecular dynamics simulations. The altered interaction with GC would cause hyper-activity of the enzyme at both low and high Ca²⁺ levels, which would ultimately lead to toxic accumulation of cGMP and Ca²⁺ in the photoreceptor outer segment, thus triggering cell death.

Structural Insights into Retinal Guanylate Cyclase Activator Proteins (GCAPs)

Retinal guanylate cyclases (RetGCs) promote the Ca²⁺-dependent synthesis of cGMP that coordinates the recovery phase of visual phototransduction in retinal rods and cones. The Ca²⁺-sensitive activation of RetGCs is controlled by a family of photoreceptor Ca²⁺ binding proteins known as guanylate cyclase activator proteins (GCAPs). The Mg²⁺-bound/Ca²⁺-free GCAPs bind to RetGCs and activate cGMP synthesis (cyclase activity) at low cytosolic Ca²⁺ levels in light-activated photoreceptors. By contrast, Ca²⁺-bound GCAPs bind to RetGCs and inactivate cyclase activity at high cytosolic Ca²⁺ levels found in dark-adapted photoreceptors. Mutations in both RetGCs and GCAPs that disrupt the Ca²⁺-dependent cyclase activity are genetically linked to various retinal diseases known as cone-rod dystrophies. In this review, I will provide an overview of the known atomic-level structures of various GCAP proteins to understand how protein dimerization and Ca²⁺-dependent conformational changes in GCAPs control the cyclase activity of RetGCs. This review will also summarize recent structural studies on a GCAP homolog from zebrafish (GCAP5) that binds to Fe²⁺ and may serve as a Fe²⁺ sensor in photoreceptors. The GCAP structures reveal an exposed hydrophobic surface that controls both GCAP1 dimerization and RetGC binding. This exposed site could be targeted by therapeutics designed to inhibit the GCAP1 disease mutants, which may serve to mitigate the onset of retinal cone-rod dystrophies.

Binding mode characterization of 13b in the monomeric and dimeric states of SARS-CoV-2 main protease using molecular dynamics simulations

The main protease, M^pro/3CL^pro, plays an essential role in processing polyproteins translated from viral RNA to produce functional viral proteins and therefore serve as an attractive target for discovering COVID-19 therapeutics. The availability of both monomer and dimer crystal bound with a common ligand, '13b' (α-ketoamide inhibitor), opened up opportunities to understand the M^pro mechanism of action. A comparative analysis of both forms of M^pro was carried out to elucidate the binding site architectural differences in the presence and absence of '13b'. Molecular dynamics simulations suggest that the presence of '13b' enhances the stability of M^pro than the unbound APO form. The N- and C- terminals of both the protomers stabilize each other, and making it's interface essential for the active form of M^pro. In comparison to monomer, the relatively high affinity of '13b' is gained in dimer pocket due to the high stability of the pocket by the interaction of S1 residue of chain B with residues F140, E166 and H172 of chain A, which is absent in monomer. The comprehensive essential dynamics, protein structure network analysis and thermodynamic profiling highlight the hot-spots, pivotal in molecular recognition process at protein-ligand and protein-protein interaction levels, cross-validated through computational alanine scanning study. A comparative description of '13b' binding mechanism in both forms illustrates valuable insights into the inhibition mechanism and the selection of critical residues suitable for the structure-based approaches for the identification of more potent M^pro inhibitors.Communicated by Ramaswamy H. Sarma.

Impaired Ca2+ Sensitivity of a Novel GCAP1 Variant Causes Cone Dystrophy and Leads to Abnormal Synaptic Transmission Between Photoreceptors and Bipolar Cells

Guanylate cyclase-activating protein 1 (GCAP1) is involved in the shutdown of the phototransduction cascade by regulating the enzymatic activity of retinal guanylate cyclase via a Ca²⁺/cGMP negative feedback. While the phototransduction-associated role of GCAP1 in the photoreceptor outer segment is widely established, its implication in synaptic transmission to downstream neurons remains to be clarified. Here, we present clinical and biochemical data on a novel isolate GCAP1 variant leading to a double amino acid substitution (p.N104K and p.G105R) and associated with cone dystrophy (COD) with an unusual phenotype. Severe alterations of the electroretinogram were observed under both scotopic and photopic conditions, with a negative pattern and abnormally attenuated b-wave component. The biochemical and biophysical analysis of the heterologously expressed N104K-G105R variant corroborated by molecular dynamics simulations highlighted a severely compromised Ca²⁺-sensitivity, accompanied by minor structural and stability alterations. Such differences reflected on the dysregulation of both guanylate cyclase isoforms (RetGC1 and RetGC2), resulting in the constitutive activation of both enzymes at physiological levels of Ca²⁺. As observed with other GCAP1-associated COD, perturbation of the homeostasis of Ca²⁺ and cGMP may lead to the toxic accumulation of second messengers, ultimately triggering cell death. However, the abnormal electroretinogram recorded in this patient also suggested that the dysregulation of the GCAP1–cyclase complex further propagates to the synaptic terminal, thereby altering the ON-pathway related to the b-wave generation. In conclusion, the pathological phenotype may rise from a combination of second messengers’ accumulation and dysfunctional synaptic communication with bipolar cells, whose molecular mechanisms remain to be clarified.

Bcl-xL Dynamics under the Lens of Protein Structure Networks

Understanding the finely orchestrated interactions leading to or preventing programmed cell death (apoptosis) is of utmost importance in cancer research because the failure of these systems could eventually lead to the onset of the disease. In this regard, the maintenance of a delicate balance between the promoters and inhibitors of mitochondrial apoptosis is crucial, as demonstrated by the interplay among the Bcl-2 family members. In particular, B-cell lymphoma extra-large (Bcl-x_L) is a target of interest due to the forefront role of its dysfunctions in cancer development. Bcl-x_L prevents apoptosis by binding both the pro-apoptotic BH3-only proteins, like PUMA, and the noncanonical partners, such as p53, at different sites. An allosteric communication between the BH3-only protein binding pocket and the p53 binding site, mediating the release of p53 from Bcl-x_L upon PUMA binding, has been postulated and supported by nuclear magnetic resonance and other biophysical data. The molecular details of this mechanism, especially at the residue level, remain unclear. In this work, we investigated the distal communication between these two sites in Bcl-x_L in its free state and when bound to PUMA. We also evaluated how missense mutations of Bcl-x_L found in cancer samples might impair this communication and therefore the allosteric mechanism. We employed all-atom explicit solvent microsecond molecular dynamics simulations, analyzed through a Protein Structure Network approach and integrated with calculations of changes in free energies upon cancer-related mutations identified by genomics studies. We found a subset of candidate residues responsible for both maintaining protein stability and for conveying structural information between the two binding sites and hypothesized possible communication routes between specific residues at both sites.

Molecular properties of human guanylate cyclase-activating protein 2 (GCAP2) and its retinal dystrophy-associated variant G157R

In murine and bovine photoreceptors, guanylate cyclase-activating protein 2 (GCAP2) activates retinal guanylate cyclases (GCs) at low Ca²⁺ levels, thus contributing to the Ca²⁺/cGMP negative feedback on the cyclase together with its paralog guanylate cyclase-activating protein 1, which has the same function but different Ca²⁺ sensitivity. In humans, a GCAP2 missense mutation (G157R) has been associated with inherited retinal degeneration (IRD) via an unknown molecular mechanism. Here, we characterized the biochemical properties of human GCAP2 and the G157R variant, focusing on its dimerization and the Ca²⁺/Mg²⁺-binding processes in the presence or absence of N-terminal myristoylation. We found that human GCAP2 and its bovine/murine orthologs significantly differ in terms of oligomeric properties, cation binding, and GC regulation. Myristoylated GCAP2 endothermically binds up to 3 Mg²⁺ with high affinity and forms a compact dimer that may reversibly dissociate in the presence of Ca²⁺. Conversely, nonmyristoylated GCAP2 does not bind Mg²⁺ over the physiological range and remains as a monomer in the absence of Ca²⁺. Both myristoylated and nonmyristoylated GCAP2 bind Ca²⁺ with high affinity. At odds with guanylate cyclase-activating protein 1 and independently of myristoylation, human GCAP2 does not significantly activate retinal GC1 in a Ca²⁺-dependent fashion. The IRD-associated G157R variant is characterized by a partly misfolded, molten globule-like conformation with reduced affinity for cations and prone to form aggregates, likely mediated by hydrophobic interactions. Our findings suggest that GCAP2 might be mostly implicated in processes other than phototransduction in human photoreceptors and suggest a possible molecular mechanism for G157R-associated IRD.

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.

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.

Neuronal Calcium Sensor GCAP1 Encoded by GUCA1A Exhibits Heterogeneous Functional Properties in Two Cases of Retinitis Pigmentosa

Genetic heterogeneity leading to retinal disorders impairs biological processes by causing, for example, severe disorder of signal transduction in photoreceptor outer segments. A normal balance of the second messenger homeostasis in photoreceptor cells seems to be a crucial factor for healthy and normal photoreceptor function. Genes like GUCY2D coding for guanylate cyclase GC-E and GUCA1A coding for the Ca²⁺-sensor guanylate cyclase-activating protein GCAP1 are critical for a precisely controlled synthesis of the second messenger cGMP. Mutations in GUCA1A frequently correlate in patients with cone dystrophy and cone-rod dystrophy. Here, we report two mutations in the GUCA1A gene that were found in patients diagnosed with retinitis pigmentosa, a phenotype that was rarely detected among previous cases of GUCA1A related retinopathies. One patient was heterozygous for the missense variant c.55C > T (p.H19Y), while the other patient was heterozygous for the missense variant c.479T > G (p.V160G). Using heterologous expression and cell culture systems, we examined the functional and molecular consequences of these point mutations. Both variants showed a dysregulation of guanylate cyclase activity, either a profound shift in Ca²⁺-sensitivity (H19Y) or a nearly complete loss of activating potency (V160G). Functional heterogeneity became also apparent in Ca²⁺/Mg²⁺-binding properties and protein conformational dynamics. A faster progression of retinal dystrophy in the patient carrying the V160G mutation seems to correlate with the more severe impairment of this variant.

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.

Mutations in PLS1, encoding fimbrin, cause autosomal dominant nonsyndromic hearing loss

Nonsyndromic hearing loss (NSHL), a common sensory disorder, is characterized by high clinical and genetic heterogeneity (i.e., approximately 115 genes and 170 loci so far identified). Nevertheless, almost half of patients submitted for genetic testing fail to receive a conclusive molecular diagnosis. We used next-generation sequencing to identify causal variants in PLS1 (c.805G>A, p.[E269K]; c.713G>T, p.[L238R], and c.383T>C, p.[F128S]) in three unrelated families of European ancestry with autosomal dominant NSHL. PLS1 encodes Plastin 1 (also called fimbrin), one of the most abundant actin-bundling proteins of the stereocilia. In silico protein modeling suggests that all variants destabilize the structure of the actin-binding domain 1, likely reducing the protein's ability to bind F actin. The role of PLS1 gene in hearing function is further supported by the recent demonstration that Pls1^{-/ -} mice show a hearing loss phenotype similar to that of our patients. In summary, we report PLS1 as a novel gene for autosomal dominant NSHL, suggesting that this gene is required for normal hearing in humans and mice.

A novel p.(Glu111Val) missense mutation in GUCA1A associated with cone-rod dystrophy leads to impaired calcium sensing and perturbed second messenger homeostasis in photoreceptors

Guanylate Cyclase-Activating Protein 1 (GCAP1) regulates the enzymatic activity of the photoreceptor guanylate cyclases (GC), leading to inhibition or activation of the cyclic guanosine monophosphate (cGMP) synthesis depending on its Ca2+- or Mg2+-loaded state. By genetically screening a family of patients diagnosed with cone-rod dystrophy, we identified a novel missense mutation with autosomal dominant inheritance pattern (c.332A>T; p.(Glu111Val); E111V from now on) in the GUCA1A gene coding for GCAP1. We performed a thorough biochemical and biophysical investigation of wild type (WT) and E111V human GCAP1 by heterologous expression and purification of the recombinant proteins. The E111V substitution disrupts the coordination of the Ca2+ ion in the high-affinity site (EF-hand 3, EF3), thus significantly decreasing the ability of GCAP1 to sense Ca2+ (∼80-fold higher Kdapp compared to WT). Both WT and E111V GCAP1 form dimers independently on the presence of cations, but the E111V Mg2+-bound form is prone to severe aggregation over time. Molecular dynamics simulations suggest a significantly increased flexibility of both the EF3 and EF4 cation binding loops for the Ca2+-bound form of E111V GCAP1, in line with the decreased affinity for Ca2+. In contrast, a more rigid backbone conformation is observed in the Mg2+-bound state compared to the WT, which results in higher thermal stability. Functional assays confirm that E111V GCAP1 interacts with the target GC with a similar apparent affinity (EC50); however, the mutant shifts the GC inhibition out of the physiological [Ca2+] (IC50E111V ∼10 μM), thereby leading to the aberrant constitutive synthesis of cGMP under conditions of dark-adapted photoreceptors.

Chokepoints in Mechanical Coupling Associated with Allosteric Proteins: The Pyruvate Kinase Example

Although the critical role of allostery in controlling enzymatic processes is well appreciated, there is a current dearth in our understanding of its underlying mechanisms, including communication between binding sites. One potential key aspect of intersite communication is the mechanical coupling between residues in a protein. Here, we introduce a graph-based computational approach to investigate the mechanical coupling between distant parts of a protein, highlighting effective pathways via which protein motion can transfer energy between sites. In this method, each residue is treated as a node on a weighted, undirected graph, in which the edges are defined by locally correlated motions of those residues and weighted by the strength of the correlation. The method was validated against experimental data on allosteric regulation in the human liver pyruvate kinase as obtained from full-protein alanine-scanning mutagenesis (systematic mutation) studies, as well as computational data on two G-protein-coupled receptors. The method provides semiquantitative information on the regulatory importance of specific structural elements. It is shown that these elements are key for the mechanical coupling between distant parts of the protein by providing effective pathways for energy transfer. It is also shown that, although there are a multitude of energy transfer pathways between distant parts of a protein, these pathways share a few common nodes that represent effective "chokepoints" for the communication.

Evolutionary-Conserved Allosteric Properties of Three Neuronal Calcium Sensor Proteins

Neuronal Calcium Sensors (NCS) are highly conserved proteins specifically expressed in neurons. Calcium (Ca²⁺)-binding to their EF-hand motifs results in a conformational change, which is crucial for the recognition of a specific target and the downstream biological process. Here we present a comprehensive analysis of the allosteric communication between Ca²⁺-binding sites and the target interfaces of three NCS, namely NCS1, recoverin (Rec), and GCAP1. In particular, Rec was investigated in different Ca²⁺-loading states and in complex with a peptide from the Rhodopsin Kinase (GRK1) while NCS1 was studied in a Ca²⁺-loaded state in complex with either the same GRK1 target or a peptide from the D₂ Dopamine receptor. A Protein Structure Network (PSN) accounting for persistent non-covalent interactions between amino acids was built for each protein state based on exhaustive Molecular Dynamics simulations. Structural network analysis helped unveiling the role of key amino acids in allosteric mechanisms and their evolutionary conservation among homologous proteins. Results for NCS1 highlighted allosteric inter-domain interactions between Ca²⁺-binding motifs and residues involved in target recognition. Robust long range, allosteric protein-target interactions were found also in Rec, in particular originating from the EF3 motif. Interestingly, Tyr 86, involved the hydrophobic packing of the N-terminal domain, was found to be a key residue for both intra- and inter-molecular communication with EF3, regardless of the presence of target or Ca²⁺ ions. Finally, based on a comprehensive topological PSN analysis for Rec, NCS1, and GCAP1 and multiple sequence alignments with homolog proteins, we propose that an evolution-driven correlation may exist between the amino acids mediating the highest number of persistent interactions (high-degree hubs) and their conservation. Such conservation is apparently fundamental for the specific structural dynamics required in signaling events.

Preferential Binding of Mg2+ Over Ca2+ to CIB2 Triggers an Allosteric Switch Impaired in Usher Syndrome Type 1J

Calcium and integrin binding protein 2 (CIB2) shares with the other members of the CIB family the ability to bind Ca²⁺ and Mg²⁺ via two functional EF-hand motifs, namely EF3 and EF4. As a cation sensor, CIB2 is able to switch to a conformation likely associated with specific biological functions yet to be clarified. Recent findings demonstrate the involvement of CIB2 in hearing physiology and a single, conservative point mutation (p.E64D) has been related to Usher Syndrome type 1J (USH1J) and non-syndromic hearing loss. We present an exhaustive biochemical and biophysical characterization of human wild type (WT) and E64D CIB2. We found that CIB2 does not possibly work as a calcium sensor under physiological conditions, its affinity for Ca²⁺ (K_d^app = 0.5 mM) being too low for detecting normal intracellular levels. Instead, CIB2 displays a significantly high affinity for Mg²⁺ (K_d^app = 290 μM), and it is probably Mg²⁺ -bound under physiological conditions. At odds with the homologous protein CIB1, CIB2 forms a non-covalent dimer under conditions that mimic the physiological ones, and as such it interacts with its physiological target α7B integrin. NMR spectroscopy revealed a long-range allosteric communication between the residue E64, located at the N-terminal domain, and the metal cation binding site EF3, located at the C-terminal domain. The conservative E64D mutation breaks up such inter-domain communication resulting in the impaired ability of CIB2 to switch to its Mg²⁺-bound form. The ability to bind the target integrin peptide was substantially conserved for E64D CIB2, thus suggesting that the molecular defect associated with USH1J resides in its inability to sense Mg²⁺ and adopt the required conformation.

Retinal guanylyl cyclase activating protein 1 forms a functional dimer

Retinal guanylyl cyclases (RetGCs) in vertebrate photoreceptors are regulated by the guanylyl cyclase activator proteins (GCAP1 and GCAP2). Here, we report EPR double electron-electron resonance (DEER) studies on the most ubiquitous GCAP isoform, GCAP1 and site-directed mutagenesis analysis to determine an atomic resolution structural model of a GCAP1 dimer. Nitroxide spin-label probes were introduced at individual GCAP1 residues: T29C, E57C, E133C, and E154C. The intermolecular distance of each spin-label probe (measured by DEER) defined restraints for calculating the GCAP1 dimeric structure by molecular docking. The DEER-derived structural model of the GCAP1 dimer was similar within the experimental error for both the Mg2+-bound activator and Ca2+-bound inhibitor states (RMSD < 2.0 Å). The GCAP1 dimer possesses intermolecular hydrophobic contacts involving the side chain atoms of H19, Y22, F73 and V77. The structural model of the dimer was validated by GCAP1 mutations (H19R, Y22D, F73E, and V77E) at the dimer interface that each abolished protein dimerization. Previous studies have shown that each of these mutants either diminished or completely suppressed the ability of GCAP1 to activate the cyclase. These results suggest that GCAP1 dimerization may affect compartmentalization of GCAP1 in the photoreceptors and/or affect regulation of the cyclase activity.

An optimal distance cutoff for contact-based Protein Structure Networks using side-chain centers of mass

Proteins are highly dynamic entities attaining a myriad of different conformations. Protein side chains change their states during dynamics, causing clashes that are propagated at distal sites. A convenient formalism to analyze protein dynamics is based on network theory using Protein Structure Networks (PSNs). Despite their broad applicability, few efforts have been devoted to benchmarking PSN methods and to provide the community with best practices. In many applications, it is convenient to use the centers of mass of the side chains as nodes. It becomes thus critical to evaluate the minimal distance cutoff between the centers of mass which will provide stable network properties. Moreover, when the PSN is derived from a structural ensemble collected by molecular dynamics (MD), the impact of the MD force field has to be evaluated. We selected a dataset of proteins with different fold and size and assessed the two fundamental properties of the PSN, i.e. hubs and connected components. We identified an optimal cutoff of 5 Å that is robust to changes in the force field and the proteins. Our study builds solid foundations for the harmonization and standardization of the PSN approach.