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

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.

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.

Myristoylation and its effects on the human Golgi Reassembly and Stacking Protein 55

GRASP55 is a myristoylated protein localized in the medial/trans-Golgi faces and involved in the Golgi structure maintenance and the regulation of unconventional secretion pathways. It is believed that GRASP55 achieves its main functionalities in the Golgi organization by acting as a tethering factor. When bound to the lipid bilayer, its orientation relative to the membrane surface is restricted to determine its proper trans-oligomerization. Despite the paramount role of myristoylation in GRASP function, the impact of such protein modification on the membrane-anchoring properties and the structural organization of GRASP remains elusive. Here, an optimized protocol for the myristoylation in E. coli of the membrane-anchoring domain of GRASP55 is presented. The biophysical properties of the myristoylated/non-myristoylated GRASP55 GRASP domain were characterized in a membrane-mimicking micellar environment. Although myristoylation did not cause any impact on the protein's secondary structure, according to our circular dichroism data, it had a significant impact on the protein's thermal stability and solubility. Electrophoresis of negatively charged liposomes incubated with the two GRASP55 constructions showed different electrophoretic mobility for the myristoylated anchored protein only, thus demonstrating that myristoylation is essential for the biological membrane anchoring. Molecular dynamics simulations were used to further explore the anchoring process in determining the restricted orientation of GRASPs in the membrane.

Decrease in DHA and other fatty acids correlates with photoreceptor degeneration in retinitis pigmentosa

Fatty acids, and especially docosahexaenoic acid (DHA), are essential for photoreceptor cell integrity and are involved in the phototransduction cascade. In this study, we analyzed the changes in the fatty acid profile in the retina of the rd10 mouse, model of retinitis pigmentosa, in order to identify potential risk factors for retinal degeneration and possible therapeutic approaches. Fatty acids from C57BL/6J and rd10 mouse retinas were extracted with Folch's method and analyzed by gas chromatography/mass spectrometry. Changes in retinal morphology were evaluated by immunohistochemistry. The rd10 mouse retina showed a decreased number of photoreceptor rows and alterations in photoreceptor morphology compared to C57BL/6J mice. The total amount of fatty acids dropped by 29.4% in the dystrophic retinas compared to C57BL/6J retinas. A positive correlation was found between the retinal content of specific fatty acids and the number of photoreceptor rows. We found that the amount of several short-chain and long-chain saturated fatty acids, as well as monounsaturated fatty acids, decreased in the retina of rd10 mice. Moreover, the content of the n-6 polyunsaturated fatty acid arachidonic acid and the n-3 polyunsaturated DHA decreased markedly in the dystrophic retina. The fall of DHA was more pronounced, hence the n-6/n-3 ratio was significantly increased in the diseased retina. The content of specific fatty acids in the retina decreased with photoreceptor degeneration in retinitis pigmentosa mice, with a remarkable reduction in DHA and other saturated and unsaturated fatty acids. These fatty acids could be essential for photoreceptor cell viability, and they should be evaluated for the design of therapeutical strategies and nutritional supplements.

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.

New Approach for Untangling the Role of Uncommon Calcium-Binding Proteins in the Central Nervous System

Although Ca²⁺ ion plays an essential role in cellular physiology, calcium-binding proteins (CaBPs) were long used for mainly as immunohistochemical markers of specific cell types in different regions of the central nervous system. They are a heterogeneous and wide-ranging group of proteins. Their function was studied intensively in the last two decades and a tremendous amount of information was gathered about them. Girard et al. compiled a comprehensive list of the gene-expression profiles of the entire EF-hand gene superfamily in the murine brain. We selected from this database those CaBPs which are related to information processing and/or neuronal signalling, have a Ca²⁺-buffer activity, Ca²⁺-sensor activity, modulator of Ca²⁺-channel activity, or a yet unknown function. In this way we created a gene function-based selection of the CaBPs. We cross-referenced these findings with publicly available, high-quality RNA-sequencing and in situ hybridization databases (Human Protein Atlas (HPA), Brain RNA-seq database and Allen Brain Atlas integrated into the HPA) and created gene expression heat maps of the regional and cell type-specific expression levels of the selected CaBPs. This represents a useful tool to predict and investigate different expression patterns and functions of the less-known CaBPs of the central nervous system.

New Approach for Untangling the Role of Uncommon Calcium-Binding Proteins in the Central Nervous System

Although Ca²⁺ ion plays an essential role in cellular physiology, calcium-binding proteins (CaBPs) were long used for mainly as immunohistochemical markers of specific cell types in different regions of the central nervous system. They are a heterogeneous and wide-ranging group of proteins. Their function was studied intensively in the last two decades and a tremendous amount of information was gathered about them. Girard et al. compiled a comprehensive list of the gene-expression profiles of the entire EF-hand gene superfamily in the murine brain. We selected from this database those CaBPs which are related to information processing and/or neuronal signalling, have a Ca²⁺-buffer activity, Ca²⁺-sensor activity, modulator of Ca²⁺-channel activity, or a yet unknown function. In this way we created a gene function-based selection of the CaBPs. We cross-referenced these findings with publicly available, high-quality RNA-sequencing and in situ hybridization databases (Human Protein Atlas (HPA), Brain RNA-seq database and Allen Brain Atlas integrated into the HPA) and created gene expression heat maps of the regional and cell type-specific expression levels of the selected CaBPs. This represents a useful tool to predict and investigate different expression patterns and functions of the less-known CaBPs of the central nervous system.

Bringing the Ca2+ sensitivity of myristoylated recoverin into the physiological range

The prototypical Ca²⁺-sensor protein recoverin (Rec) is thought to regulate the activity of rhodopsin kinase (GRK1) in photoreceptors by switching from a relaxed (R) disc membrane-bound conformation in the dark to a more compact, cytosol-diffusing tense (T) conformation upon cell illumination. However, the apparent affinity for Ca²⁺ of its physiologically relevant form (myristoylated recoverin) is almost two orders of magnitude too low to support this mechanism in vivo. In this work, we compared the individual and synergistic roles of the myristic moiety, the GRK1 target and the disc membrane in modulating the calcium sensitivity of Rec. We show that the sole presence of the target or the disc membrane alone are not sufficient to achieve a physiological response to changes in intracellular [Ca²⁺]. Instead, the simultaneous presence of GRK1 and membrane allows the T to R transition to occur in a physiological range of [Ca²⁺] with high cooperativity via a conformational selection mechanism that drives the structural transitions of Rec in the presence of multiple ligands. Our conclusions may apply to other sensory transduction systems involving protein complexes and biological membranes.

Affinity of rhodopsin to raft enables the aligned oligomer formation from dimers: Coarse-grained molecular dynamics simulation of disk membranes

The visual photopigment protein rhodopsin (Rh) is a typical G protein-coupled receptor (GPCR) that initiates the phototransduction cascade in retinal disk membrane of rod-photoreceptor cells. Rh molecule has a tendency to form dimer, and the dimer tends to form rows, which is suggested to heighten phototransduction efficiency in single-photon regime. In addition, the dimerization confers Rh an affinity for lipid raft, i.e. raftophilicity. However, the mechanism by which Rh-dimer raftophilicity contributes to the organization of the higher order structure remains unknown. In this study, we performed coarse-grained molecular dynamics simulations of a disk membrane model containing unsaturated lipids, saturated lipids with cholesterol, and Rh-dimers. We described the Rh-dimers by two-dimensional particle populations where the palmitoyl moieties of each Rh exhibits raftophilicity. We simulated the structuring of Rh in a disk for two types of Rh-dimer, i.e., the most and second most stable Rh dimers, which exposes the raftophilic regions at the dimerization-interface (H1/H8 dimer) and two edges away from the interface (H4/H5 dimer), respectively. Our simulations revealed that only the H1/H8 dimer could form a row structure. A small number of raftophilic lipids recruited to and intercalated in a narrow space between H1/H8 dimers stabilize the side-by-side interaction between dimers in a row. Our results implicate that the nano-sized lipid raft domains act as a “glue” to organize the long row structures of Rh-dimers.