Z-DNA and Z-RNA are functionally important left-handed structures of nucleic acids, which play a significant role in several molecular and biological processes including DNA replication, gene expression regulation and viral nucleic acid sensing. Most proteins that have been proven to interact with Z-DNA/Z-RNA contain the so-called Zα domain, which is structurally well conserved. To date, only eight proteins with Zα domain have been described within a few organisms (including human, mouse, Danio rerio, Trypanosoma brucei and some viruses). Therefore, this paper aimed to search for new Z-DNA/Z-RNA binding proteins in the complete PDB structures database and from the AlphaFold2 protein models. A structure-based similarity search found 14 proteins with highly similar Zα domain structure in experimentally-defined proteins and 185 proteins with a putative Zα domain using the AlphaFold2 models. Structure-based alignment and molecular docking confirmed high functional conservation of amino acids involved in Z-DNA/Z-RNA, suggesting that Z-DNA/Z-RNA recognition may play an important role in a variety of cellular processes.

SARS-CoV-2 is a novel positive-sense single-stranded RNA virus from the Coronaviridae family (genus Betacoronavirus), which has been established as causing the COVID-19 pandemic. The genome of SARS-CoV-2 is one of the largest among known RNA viruses, comprising of at least 26 known protein-coding loci. Studies thus far have outlined the coding capacity of the positive-sense strand of the SARS-CoV-2 genome, which can be used directly for protein translation. However, it has been recently shown that transcribed negative-sense viral RNA intermediates that arise during viral genome replication from positive-sense viruses can also code for proteins. No studies have yet explored the potential for negative-sense SARS-CoV-2 RNA intermediates to contain protein-coding loci. Thus, using sequence and structure-based bioinformatics methodologies, we have investigated the presence and validity of putative negative-sense ORFs (nsORFs) in the SARS-CoV-2 genome. Nine nsORFs were discovered to contain strong eukaryotic translation initiation signals and high codon adaptability scores, and several of the nsORFs were predicted to interact with RNA-binding proteins. Evolutionary conservation analyses indicated that some of the nsORFs are deeply conserved among related coronaviruses. Three-dimensional protein modeling revealed the presence of higher order folding among all putative SARS-CoV-2 nsORFs, and subsequent structural mimicry analyses suggest similarity of the nsORFs to DNA/RNA-binding proteins and proteins involved in immune signaling pathways. Altogether, these results suggest the potential existence of still undescribed SARS-CoV-2 proteins, which may play an important role in the viral lifecycle and COVID-19 pathogenesis.

Voltage-gated sodium channel α-subunits (NaV1.1–1.9) initiate and propagate action potentials in neurons and myocytes. The NaV β-subunits (β1–4) have been shown to modulate α-subunit properties. Homo-oligomerization of β-subunits on neighboring or opposing plasma membranes has been suggested to facilitate cis or trans interactions, respectively. The interactions between several NaV channel isoforms and β-subunits have been determined using cryogenic electron microscopy (cryo-EM). Interestingly, the NaV cryo-EM structures reveal the presence of N-linked glycosylation sites. However, only the first glycan moieties are typically resolved at each site due to the flexibility of mature glycan trees. Thus, existing cryo-EM structures may risk de-emphasizing the structural implications of glycans on the NaV channels. Herein, molecular modeling and all-atom molecular dynamics simulations were applied to investigate the conformational landscape of N-linked glycans on NaV channel surfaces. The simulations revealed that negatively charged sialic acid residues of two glycan sites may interact with voltage-sensing domains. Notably, two NaV1.5 isoform-specific glycans extensively cover the α-subunit region that, in other NaV channel α-subunit isoforms, corresponds to the binding site for the β1- (and likely β3-) subunit immunoglobulin (Ig) domain. NaV1.8 contains a unique N-linked glycosylation site that likely prevents its interaction with the β2 and β4-subunit Ig-domain. These isoform-specific glycans may have evolved to facilitate specific functional interactions, for example, by redirecting β-subunit Ig-domains outward to permit cis or trans supraclustering within specialized cellular compartments such as the cardiomyocyte perinexal space. Further experimental work is necessary to validate these predictions.

Glycine receptors (GlyRs) are glycine-gated inhibitory pentameric ligand-gated ion channels composed of α or α + β subunits. A number of structures of these proteins have been reported, but to date, these have only revealed details of the extracellular and transmembrane domains, with the intracellular domain (ICD) remaining uncharacterised due to its high flexibility. The ICD is a region that can modulate function in addition to being critical for receptor localisation and clustering via proteins such as gephyrin. Here, we use modelling and molecular dynamics (MD) to reveal details of the ICDs of both homomeric and heteromeric GlyR. At their N and C ends, both the α and β subunit ICDs have short helices, which are major sites of stabilising interactions; there is a large flexible loop between them capable of forming transient secondary structures. The α subunit can affect the β subunit ICD structure, which is more flexible in a 4α₂:1β than in a 4α₁:1β GlyR. We also explore the effects of gephyrin binding by creating GlyR models bound to the gephyrin E domain; MD simulations suggest these are more stable than the unbound forms, and again there are α subunit-dependent differences, despite the fact the gephyrin binds to the β subunit. The bound models also suggest that gephyrin causes compaction of the ICD. Overall, the data expand our knowledge of this important receptor protein and in particular clarify features of the underexplored ICD.

The RGD motif on the SARS-CoV-2 spike protein has been suggested to interact with RGD-binding integrins αVβ3 and α5β1 to enhance viral cell entry and alter downstream signaling cascades. The D405N mutation on the Omicron subvariant spike proteins, resulting in an RGN motif, has recently been shown to inhibit binding to integrin αVβ3. Deamidation of asparagines in protein ligand RGN motifs has been demonstrated to generate RGD and RGisoD motifs that permit binding to RGD-binding integrins. Two asparagines, N481 and N501, on the Wild-type spike receptor-binding domain have been previously shown to have deamidation half-lives of 16.5 and 123 days, respectively, which may occur during the viral life cycle. Deamidation of Omicron subvariant N405 may recover the ability to interact with RGD-binding integrins. Thus, herein, all-atom molecular dynamics simulations of the Wild-type and Omicron subvariant spike protein receptor-binding domains were conducted to investigate the potential for asparagines, the Omicron subvariant N405 in particular, to assume the optimized geometry for deamidation to occur. In summary, the Omicron subvariant N405 was primarily found to be stabilized in a state unfavourable for deamidation after hydrogen bonding with downstream E406. Nevertheless, a small number of RGD or RGisoD motifs on the Omicron subvariant spike proteins may restore the ability to interact with RGD-binding integrins. The simulations also provided structural clarification regarding the deamidation rates of Wild-type N481 and N501 and highlighted the utility of tertiary structure dynamics information in predicting asparagine deamidation. Further work is needed to characterize the effects of deamidation on spike-integrin interactions.

Bioelectrical signalling is fundamental for regulating biological processes in all forms of life. Ion channels and transporters generate and propagate electrical currents by selectively allowing ions to flow across membranes in response to voltage changes. Although recent breakthroughs in structural determination methods, such as cryogenic electron microscopy, have provided novel insights into the structure–function relationships of ion channels and scaffolding proteins, their precise roles in bioelectrical signal generation and propagation within and across different cells and tissues remain unresolved. This article examines the biochemical and ultrastructural features of the three most studied modes of bioelectrical conduction in human tissues—electrotonic, saltatory and ephaptic conduction—and how biophysical constraints set by membranes and proteins give rise to bioelectricity. Notably, ion channel clustering and scaffolding proteins that define intermembrane distances are common key features among all forms of bioelectrical signalling. Techniques like cryogenic electron tomography offer promising avenues for exploring ion channels and their regulatory protein interactions in situ. The central question is: ‘How does the spatial organization of ions, molecules and tissues give rise to bioelectricity?’ These insights may inform novel therapeutic approaches for various diseases, while also potentially offering new perspectives on life, evolution and consciousness.