Structural and functional insights into GSDMB isoforms complex roles in pathogenesis

SHADSGasdermins (GSDMs) have garnered significant scientific interest due to their protective and detrimental roles in innate immunity, host defense, inflammation, and cancer alongside with other pathologies. While GSDMs are mostly recognized as key effectors of a lytic type of pro-inflammatory cell death known as pyroptosis, they do also take part in other cell death processes (NETosis, secondary necrosis, or apoptosis) and exhibit cell-death independent functions depending on the cellular context. Among GSDMs, Gasdermin B (GSDMB) pyroptotic capacity has been a subject of conflicting findings in scientific literature even when its processing, and subsequent activation, by Granzyme A (GZMA) was decoded. Nevertheless, recent groundbreaking publications have shed light on the crucial role of alternative splicing in determining the pyroptotic capacity of GSDMB isoforms, which depends on the presence of exon 6-derived elements. This comprehensive review pays attention to the relevant structural differences among recently crystalized GSDMB isoforms. As a novelty, the structural aspects governing GSDMB isoform susceptibility to GZMA-mediated activation have been investigated. By elucidating the complex roles of GSDMB isoforms, this review aims to deepen the understanding of this multifunctional player and its potential implications in disease pathogenesis and therapeutic interventions. [Figure: see text].

Biallelic variants in ZNF142 lead to a syndromic neurodevelopmental disorder

Biallelic variants of the gene encoding for the zinc-finger protein 142 (ZNF142) have recently been associated with intellectual disability (ID), speech impairment, seizures, and movement disorders in nine individuals from five families. In this study, we obtained phenotype and genotype information of 26 further individuals from 16 families. Among the 27 different ZNF142 variants identified in the total of 35 individuals only four were missense. Missense variants may give a milder phenotype by changing the local structure of ZF motifs as suggested by protein modelling; but this correlation should be validated in larger cohorts and pathogenicity of the missense variants should be investigated with functional studies. Clinical features of the 35 individuals suggest that biallelic ZNF142 variants lead to a syndromic neurodevelopmental disorder with mild to moderate ID, varying degrees of delay in language and gross motor development, early onset seizures, hypotonia, behavioral features, movement disorders, and facial dysmorphism. The differences in symptom frequencies observed in the unpublished individuals compared to those of published, and recognition of previously underemphasized facial features are likely to be due to the small sizes of the previous cohorts, which underlines the importance of larger cohorts for the phenotype descriptions of rare genetic disorders. This article is protected by copyright. All rights reserved.

The clinical and molecular spectrum of QRICH1 associated neurodevelopmental disorder

De novo variants in QRICH1 (Glutamine-rich protein 1) has recently been reported in 11 individuals with intellectual disability (ID). The function of QRICH1 is largely unknown but it is likely to play a key role in the unfolded response of endoplasmic reticulum stress through transcriptional control of proteostasis. In this study, we present 27 additional individuals and delineate the clinical and molecular spectrum of the individuals (n = 38) with QRICH1 variants. The main clinical features were mild to moderate developmental delay/ID (71%), nonspecific facial dysmorphism (92%) and hypotonia (39%). Additional findings included poor weight gain (29%), short stature (29%), autism spectrum disorder (29%), seizures (24%) and scoliosis (18%). Minor structural brain abnormalities were reported in 52% of the individuals with brain imaging. Truncating or splice variants were found in 28 individuals and 10 had missense variants. Four variants were inherited from mildly affected parents. This study confirms that heterozygous QRICH1 variants cause a neurodevelopmental disorder including short stature and expands the phenotypic spectrum to include poor weight gain, scoliosis, hypotonia, minor structural brain anomalies, and seizures. Inherited variants from mildly affected parents are reported for the first time, suggesting variable expressivity.

Dysfunctional Homozygous VRK1-D263G Variant Impairs the Assembly of Cajal Bodies and DNA Damage Response in Hereditary Spastic Paraplegia

Background and Objectives

To conduct a genetic and molecular functional study of a family with members affected of hereditary spastic paraplegia (HSP) of unknown origin and carrying a novel pathogenic vaccinia-related kinase 1 (VRK1) variant.

Methods

Whole-exome sequencing was performed in 2 patients, and their parents diagnosed with HSP. The novel VRK1 variant was detected by whole-exome sequencing, molecularly modeled and biochemically characterized in kinase assays. Functionally, we studied the role of this VRK1 variant in DNA damage response and its effect on the assembly of Cajal bodies (CBs).

Results

We have identified a very rare homozygous variant VRK1-D263G with a neurologic phenotype associated with HSP and moderate intellectual disability. The molecular modeling of this VRK1 variant protein predicted an alteration in the folding of a loop that interferes with the access to the kinase catalytic site. The VRK1-D263G variant is kinase inactive and does not phosphorylate histones H2AX and H3, transcription factors activating transcription factor 2 and p53, coilin needed for assembly of CBs, and p53 binding protein 1, a DNA repair protein. Functionally, this VRK1 variant protein impairs CB formation and the DNA damage response.

Discussion

This report expands the neurologic spectrum of neuromotor syndromes associated with a new and rare VRK1 variant, representing a novel pathogenic participant in complicated HSP and demonstrates that CBs and the DNA damage response are impaired in these patients.

Comparing molecular dynamics force fields in the essential subspace

The continued development and utility of molecular dynamics simulations requires improvements in both the physical models used (force fields) and in our ability to sample the Boltzmann distribution of these models. Recent developments in both areas have made available multi-microsecond simulations of two proteins, ubiquitin and Protein G, using a number of different force fields. Although these force fields mostly share a common mathematical form, they differ in their parameters and in the philosophy by which these were derived, and previous analyses showed varying levels of agreement with experimental NMR data. To complement the comparison to experiments, we have performed a structural analysis of and comparison between these simulations, thereby providing insight into the relationship between force-field parameterization, the resulting ensemble of conformations and the agreement with experiments. In particular, our results show that, at a coarse level, many of the motional properties are preserved across several, though not all, force fields. At a finer level of detail, however, there are distinct differences in both the structure and dynamics of the two proteins, which can, together with comparison with experimental data, help to select force fields for simulations of proteins. A noteworthy observation is that force fields that have been reparameterized and improved to provide a more accurate energetic description of the balance between helical and coil structures are difficult to distinguish from their "unbalanced" counterparts in these simulations. This observation implies that simulations of stable, folded proteins, even those reaching 10 microseconds in length, may provide relatively little information that can be used to modify torsion parameters to achieve an accurate balance between different secondary structural elements.

Determination of key residues for catalysis and RNA cleavage specificity: one mutation turns RNase II into a "SUPER-ENZYME"

RNase II is the prototype of a ubiquitous family of enzymes that are crucial for RNA metabolism. In Escherichia coli this protein is a single-stranded-specific 3'-exoribonuclease with a modular organization of four functional domains. In eukaryotes, the RNase II homologue Rrp44 (also known as Dis3) is the catalytic subunit of the exosome, an exoribonuclease complex essential for RNA processing and decay. In this work we have performed a functional characterization of several highly conserved residues located in the RNase II catalytic domain to address their precise role in the RNase II activity. We have constructed a number of RNase II mutants and compared their activity and RNA binding to the wild type using different single- or double-stranded substrates. The results presented in this study substantially improve the RNase II model for RNA degradation. We have identified the residues that are responsible for the discrimination of cleavage of RNA versus DNA. We also show that the Arg-500 residue present in the RNase II active site is crucial for activity but not for RNA binding. The most prominent finding presented is the extraordinary catalysis observed in the E542A mutant that turns RNase II into a "super-enzyme."

Structural properties of the human respiratory syncytial virus P protein: evidence for an elongated homotetrameric molecule that is the smallest orthologue within the family of paramyxovirus polymerase cofactors

The oligomeric state and the hydrodynamic properties of human respiratory syncytial virus (HRSV) phosphoprotein (P), a known cofactor of the viral RNA-dependent RNA polymerase (L), and a trypsin-resistant fragment (X) that includes its oligomerization domain were analyzed by sedimentation equilibrium and velocity using analytical ultracentrifugation. The results obtained demonstrate that both P and fragment X are homotetrameric with elongated shapes, consistent with electron micrographs of the purified P protein in which thin rod-like molecules of approximately 12.5 +/- 1.0 nm in length were observed. A new chymotrypsin resistant fragment (Y*) included in fragment X has been identified and purified by gel filtration chromatography. Fragment Y* may represent a minimal version of the P oligomerization domain. Thermal denaturation curves based on circular dichroism data of P protein showed a complex behavior. In contrast, melting data generated for fragments X and particularly fragment Y* showed more homogeneous transitions indicative of simpler structures. A three-dimensional model of X and Y* fragments was built based on the atomic structure of the P oligomerization domain of the related Sendai virus, which is in good agreement with the experimental data. This model will be an useful tool to make rational mutations and test the role of specific amino acids in the oligomerization and functional properties of the HRSV P protein.

New insights into the mechanism of RNA degradation by ribonuclease II: identification of the residue responsible for setting the RNase II end product

RNase II is a key exoribonuclease involved in the maturation, turnover, and quality control of RNA. RNase II homologues are components of the exosome, a complex of exoribonucleases. The structure of RNase II unraveled crucial aspects of the mechanism of RNA degradation. Here we show that mutations in highly conserved residues at the active site affect the activity of the enzyme. Moreover, we have identified the residue that is responsible for setting the end product of RNase II. In addition, we present for the first time the models of two members of the RNase II family, RNase R from Escherichia coli and human Rrp44, also called Dis3. Our findings improve the present model for RNA degradation by the RNase II family of enzymes.

The role of the S1 domain in exoribonucleolytic activity: substrate specificity and multimerization

RNase II is a 3'-5' exoribonuclease that processively hydrolyzes single-stranded RNA generating 5' mononucleotides. This enzyme contains a catalytic core that is surrounded by three RNA-binding domains. At its C terminus, there is a typical S1 domain that has been shown to be critical for RNA binding. The S1 domain is also present in the other major 3'-5' exoribonucleases from Escherichia coli: RNase R and polynucleotide phosphorylase (PNPase). In this report, we examined the involvement of the S1 domain in the different abilities of these three enzymes to overcome RNA secondary structures during degradation. Hybrid proteins were constructed by replacing the S1 domain of RNase II for the S1 from RNase R and PNPase, and their exonucleolytic activity and RNA-binding ability were examined. The results revealed that both the S1 domains of RNase R and PNPase are able to partially reverse the drop of RNA-binding ability and exonucleolytic activity resulting from removal of the S1 domain of RNase II. Moreover, the S1 domains investigated are not equivalent. Furthermore, we demonstrate that S1 is neither responsible for the ability to overcome secondary structures during RNA degradation, nor is it related to the size of the final product generated by each enzyme. In addition, we show that the S1 domain from PNPase is able to induce the trimerization of the RNaseII-PNP hybrid protein, indicating that this domain can have a role in the biogenesis of multimers.

Ca2+ Activation kinetics of the two aspartate-glutamate mitochondrial carriers, aralar and citrin: role in the heart malate-aspartate NADH shuttle

Ca(2+) regulation of the Ca(2+) binding mitochondrial carriers for aspartate/glutamate (AGCs) is provided by their N-terminal extensions, which face the intermembrane space. The two mammalian AGCs, aralar and citrin, are members of the malate-aspartate NADH shuttle. We report that their N-terminal extensions contain up to four pairs of EF-hand motifs plus a single vestigial EF-hand, and have no known homolog. Aralar and citrin contain one fully canonical EF-hand pair and aralar two additional half-pairs, in which a single EF-hand is predicted to bind Ca(2+). Shuttle activity in brain or skeletal muscle mitochondria, which contain aralar as the major AGC, is activated by Ca(2+) with S(0.5) values of 280-350 nm; higher than those obtained in liver mitochondria (100-150 nm) that contain citrin as the major AGC. We have used aralar- and citrin-deficient mice to study the role of the two isoforms in heart, which expresses both AGCs. The S(0.5) for Ca(2+) activation of the shuttle in heart mitochondria is about 300 nm, and it remains essentially unchanged in citrin-deficient mice, although it undergoes a drastic reduction to about 100 nm in aralar-deficient mice. Therefore, aralar and citrin, when expressed as single isoforms in heart, confer differences in Ca(2+) activation of shuttle activity, probably associated with their structural differences. In addition, the results reveal that the two AGCs fully account for shuttle activity in mouse heart mitochondria and that no other glutamate transporter can replace the AGCs in this pathway.

A single-residue mutation, G203E, causes 3-hydroxy-3-methylglutaric aciduria by occluding the substrate channel in the 3D structural model of HMG-CoA lyase

3-Hydroxy-3-methylglutaric aciduria is a rare autosomal recessive genetic disorder that affects ketogenesis and leucine metabolism. The disease is caused by mutations in the gene coding for 3-hydroxy-3-methylglutaryl-coenzyme A lyase (HL). To date 26 different mutations have been described. A (betaalpha)(8) TIM barrel structure has been proposed for the protein, and almost all missense mutations identified so far localize in the beta sheets that define the inside cavity. We report an Italian patient who bears homozygously a novel HL mutation, c.608G > A (p. G203E) in beta sheet six. A structural model of the mutated protein suggests that glutamic acid 203 impedes catalysis by blocking the entrance to the inner cavity of the enzyme. Loss of functionality has been confirmed in expression studies in E. coli, which demonstrate that the G203E mutation completely abolishes enzyme activity. Beta sheet six and beta sheet two are the two protein regions that accumulate most missense mutations, indicating their importance in enzyme functionality. A model for the mechanism of enzyme function is proposed.

Computational Approaches to Model Ligand Selectivity in Drug Design

To be effective, a designed drug must discriminate successfully the macromolecular target from alternative structures present in the organism. The last few years have witnessed the emergence of different computational tools aimed to the understanding and modeling of this process at molecular level. Although still rudimentary, these methods are shaping a coherent approach to help in the design of molecules with high affinity and specificity, both in lead discovery and in lead optimization. It is the purpose of this review to illustrate the array of computational tools available to consider selectivity in the design process, to summarize the most relevant applications, and to sketch the challenges ahead.

TRAF family proteins link PKR with NF-kappa B activation

The double-stranded RNA (dsRNA)-dependent protein kinase PKR activates NF-kappa B via the I kappa B kinase (IKK) complex, but little is known about additional molecules that may be involved in this pathway. Analysis of the PKR sequence enabled us to identify two putative TRAF-interacting motifs. The viability of such an interaction was further suggested by computer modeling. Here, we present evidence of the colocalization and physical interaction between PKR and TRAF family proteins in vivo, as shown by immunoprecipitation and confocal microscopy experiments. This interaction is induced upon PKR dimerization. Most importantly, we show that the binding between PKR and TRAFs is functionally relevant, as observed by the absence of NF-kappa B activity upon PKR expression in cells genetically deficient in TRAF2 and TRAF5 or after expression of TRAF dominant negative molecules. On the basis of sequence information and mutational and computer docking analyses, we favored a TRAF-PKR interaction model in which the C-terminal domain of TRAF binds to a predicted TRAF interaction motif present in the PKR kinase domain. Altogether, our data suggest that TRAF family proteins are key components located downstream of PKR that have an important role in mediating activation of NF-kappa B by the dsRNA-dependent protein kinase.

A comparative sequence analysis to revise the current taxonomy of the family Coronaviridae

The Coronaviridae family, comprising the Coronavirus and Torovirus genera, is part of the Nidovirales order that also includes two other families, Arteriviridae and Roniviridae. Based on genetic and serological relationships, groups 1, 2 and 3 were previously recognized in the Coronavirus genus. In this report we present results of comparative sequence analysis of the spike (S), envelope (E), membrane (M), and nucleoprotein (N) structural proteins, and the two most conserved replicase domains, putative RNA-dependent RNA polymerase (RdRp) and RNA helicase (HEL), aimed at a revision of the Coronaviridae taxonomy. The results of pairwise comparisons involving structural and replicase proteins of the Coronavirus genus were consistent and produced percentages of sequence identities that were distributed in discontinuous clusters. Inter-group pairwise scores formed a single cluster in the lowest percentile. No homologs of the N and E proteins have been found outside coronaviruses, and the only (very) distant homologs of S and M proteins were identified in toroviruses. Intragroup sequence conservation was higher, although for some pairs, especially those from the most diverse group 1, scores were close or even overlapped with those from the intergroup comparisons. Phylogenetic analysis of six proteins using a neighbor-joining algorithm confirmed three coronavirus groups. Comparative sequence analysis of RdRp and HEL domains were extended to include arterivirus and ronivirus homologs. The pairwise scores between sequences of the genera Coronavirus and Torovirus (22–25% and 21–25%) were found to be very close to or overlapped with the value ranges (12 to 22% and 17 to 25%) obtained for interfamily pairwise comparisons, but were much smaller than values derived from pairwise comparisons within the Coronavirus genus (63–71% and 59–67%). Phylogenetic analysis confirmed toroviruses and coronaviruses to be separated by a large distance that is comparable to those between established nidovirus families. Based on comparison of these scores with those derived from analysis of separate ranks of several multi-genera virus families, like the Picornaviridae, a revision of the Coronaviridae taxonomy is proposed. We suggest the Coronavirus and Torovirus genera to be re-defined as two subfamilies within the Coronavirdae or two families within Nidovirales, and the current three informal coronavirus groups to be converted into three genera within the Coronaviridae.