Bioinformatics in protein analysis

Modeling SARS-CoV-2 spike/ACE2 protein-protein interactions for predicting the binding affinity of new spike variants for ACE2, and novel ACE2 structurally related human protein targets, for COVID-19 handling in the 3PM context

Aims

The rapid spread of new SARS-CoV-2 variants has highlighted the crucial role played in the infection by mutations occurring at the SARS-CoV-2 spike receptor binding domain (RBD) in the interactions with the human ACE2 receptor. In this context, it urgently needs to develop new rapid tools for quickly predicting the affinity of ACE2 for the SARS-CoV-2 spike RBD protein variants to be used with the ongoing SARS-CoV-2 genomic sequencing activities in the clinics, aiming to gain clues about the transmissibility and virulence of new variants, to prevent new outbreaks and to quickly estimate the severity of the disease in the context of the 3PM.

Methods

In our study, we used a computational pipeline for calculating the interaction energies at the SARS-CoV-2 spike RBD/ACE2 protein-protein interface for a selected group of characterized infectious variants of concern/interest (VoC/VoI). By using our pipeline, we built 3D comparative models of the SARS-CoV-2 spike RBD/ACE2 protein complexes for the VoC B.1.1.7-United Kingdom (carrying the mutations of concern/interest N501Y, S494P, E484K at the RBD), P.1-Japan/Brazil (RBD mutations: K417T, E484K, N501Y), B.1.351-South Africa (RBD mutations: K417N, E484K, N501Y), B.1.427/B.1.429-California (RBD mutations: L452R), the B.1.141 (RBD mutations: N439K), and the recent B.1.617.1-India (RBD mutations: L452R; E484Q) and the B.1.620 (RBD mutations: S477N; E484K). Then, we used the obtained 3D comparative models of the SARS-CoV-2 spike RBD/ACE2 protein complexes for predicting the interaction energies at the protein-protein interface.

Results

Along SARS-CoV-2 mutation database screening and mutation localization analysis, it was ascertained that the most dangerous mutations at VoC/VoI spike proteins are located mainly at three regions of the SARS-CoV-2 spike "boat-shaped" receptor binding motif, on the RBD domain. Notably, the P.1 Japan/Brazil variant present three mutations, K417T, E484K, N501Y, located along the entire receptor binding motif, which apparently determines the highest interaction energy at the SARS-CoV-2 spike RBD/ACE2 protein-protein interface, among those calculated. Conversely, it was also observed that the replacement of a single acidic/hydrophilic residue with a basic residue (E484K or N439K) at the "stern" or "bow" regions, of the boat-shaped receptor binding motif on the RBD, appears to determine an interaction energy with ACE2 receptor higher than that observed with single mutations occurring at the "hull" region or with other multiple mutants. In addition, our pipeline allowed searching for ACE2 structurally related proteins, i.e., THOP1 and NLN, which deserve to be investigated for their possible involvement in interactions with the SARS-CoV-2 spike protein, in those tissues showing a low expression of ACE2, or as a novel receptor for future spike variants. A freely available web-tool for the in silico calculation of the interaction energy at the SARS-CoV-2 spike RBD/ACE2 protein-protein interface, starting from the sequences of the investigated spike and/or ACE2 variants, was made available for the scientific community at: https://www.mitoairm.it/covid19affinities.

Conclusion

In the context of the PPPM/3PM, the employment of the described pipeline through the provided webservice, together with the ongoing SARS-CoV-2 genomic sequencing, would help to predict the transmissibility of new variants sequenced from future patients, depending on SARS-CoV-2 genomic sequencing activities and on the specific amino acid replacement and/or on its location on the SARS-CoV-2 spike RBD, to put in play all the possible counteractions for preventing the most deleterious scenarios of new outbreaks, taking into consideration that a greater transmissibility has not to be necessarily related to a more severe manifestation of the disease.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13167-021-00267-w.

Zoledronic Acid as a Novel Dual Blocker of KIR6.1/2-SUR2 Subunits of ATP-Sensitive K+ Channels: Role in the Adverse Drug Reactions

Zoledronic acid (ZOL) is used as a bone-specific antiresorptive drug with antimyeloma effects. Adverse drug reactions (A.D.R.) are associated with ZOL-therapy, whose mechanics are unknown. ZOL is a nitrogen-containing molecule whose structure shows similarities with nucleotides, ligands of ATP-sensitive K⁺ (KATP) channels. We investigated the action of ZOL by performing in vitro patch-clamp experiments on native KATP channels in murine skeletal muscle fibers, bone cells, and recombinant subunits in cell lines, and by in silico docking the nucleotide site on KIR and SUR, as well as the glibenclamide site. ZOL fully inhibited the KATP currents recorded in excised macro-patches from Extensor digitorum longus (EDL) and Soleus (SOL) muscle fibers with an IC₅₀ of 1.2 ± 1.4 × 10⁻⁶ and 2.1 ± 3.7 × 10⁻¹⁰ M, respectively, and the KATP currents recorded in cell-attached patches from primary long bone cells with an IC₅₀ of 1.6 ± 2.8 × 10⁻¹⁰ M. ZOL fully inhibited a whole-cell KATP channel current of recombinant KIR6.1-SUR2B and KIR6.2-SUR2A subunits expressed in HEK293 cells with an IC₅₀ of 3.9 ± 2.7 × 10⁻¹⁰ M and 7.1 ± 3.1 × 10⁻⁶ M, respectively. The rank order of potency in inhibiting the KATP currents was: KIR6.1-SUR2B/SOL-KATP/osteoblast-KATP > KIR6.2-SUR2A/EDL-KATP >>> KIR6.2-SUR1 and KIR6.1-SUR1. Docking investigation revealed that the drug binds to the ADP/ATP sites on KIR6.1/2 and SUR2A/B and on the sulfonylureas site showing low binding energy <6 Kcal/mol for the KIR6.1/2-SUR2 subunits vs. the <4 Kcal/mol for the KIR6.2-SUR1. The IC₅₀ of ZOL to inhibit the KIR6.1/2-SUR2A/B channels were correlated with its musculoskeletal and cardiovascular risks. We first showed that ZOL blocks at subnanomolar concentration musculoskeletal KATP channels and cardiac and vascular KIR6.2/1-SUR2 channels.

The T Cell Receptor (TRB) Locus in Tursiops truncatus: From Sequence to Structure of the Alpha/Beta Heterodimer in the Human/Dolphin Comparison

The bottlenose dolphin (Tursiops truncatus) belongs to the Cetartiodactyla and, similarly to other cetaceans, represents the most successful mammalian colonization of the aquatic environment. Here we report a genomic, evolutionary, and expression study of T. truncatus T cell receptor beta (TRB) genes. Although the organization of the dolphin TRB locus is similar to that of the other artiodactyl species, with three in tandem D-J-C clusters located at its 3' end, its uniqueness is given by the reduction of the total length due essentially to the absence of duplications and to the deletions that have drastically reduced the number of the germline TRBV genes. We have analyzed the relevant mature transcripts from two subjects. The simultaneous availability of rearranged T cell receptor α (TRA) and TRB cDNA from the peripheral blood of one of the two specimens, and the human/dolphin amino acids multi-sequence alignments, allowed us to calculate the most likely interactions at the protein interface between the alpha/beta heterodimer in complex with major histocompatibility class I (MH1) protein. Interacting amino acids located in the complementarity-determining region according to IMGT numbering (CDR-IMGT) of the dolphin variable V-alpha and beta domains were identified. According to comparative modelization, the atom pair contact sites analysis between the human MH1 grove (G) domains and the T cell receptor (TR) V domains confirms conservation of the structure of the dolphin TR/pMH.

Protein structure analysis of the interactions between SARS-CoV-2 spike protein and the human ACE2 receptor: from conformational changes to novel neutralizing antibodies

The recent severe acute respiratory syndrome, known as Coronavirus Disease 2019 (COVID-19) has spread so much rapidly and severely to induce World Health Organization (WHO) to declare a state of emergency over the new coronavirus SARS-CoV-2 pandemic. While several countries have chosen the almost complete lock-down for slowing down SARS-CoV-2 spread, the scientific community is called to respond to the devastating outbreak by identifying new tools for diagnosis and treatment of the dangerous COVID-19. With this aim, we performed an in silico comparative modeling analysis, which allows gaining new insights into the main conformational changes occurring in the SARS-CoV-2 spike protein, at the level of the receptor-binding domain (RBD), along interactions with human cells angiotensin-converting enzyme 2 (ACE2) receptor, that favor human cell invasion. Furthermore, our analysis provides (1) an ideal pipeline to identify already characterized antibodies that might target SARS-CoV-2 spike RBD, aiming to prevent interactions with the human ACE2, and (2) instructions for building new possible neutralizing antibodies, according to chemical/physical space restraints and complementary determining regions (CDR) mutagenesis of the identified existing antibodies. The proposed antibodies show in silico high affinity for SARS-CoV-2 spike RBD and can be used as reference antibodies also for building new high-affinity antibodies against present and future coronaviruses able to invade human cells through interactions of their spike proteins with the human ACE2. More in general, our analysis provides indications for the set-up of the right biological molecular context for investigating spike RBD–ACE2 interactions for the development of new vaccines, diagnostic kits, and other treatments based on the targeting of SARS-CoV-2 spike protein.

Using Phylogenetic Analysis to Investigate Eukaryotic Gene Origin

Phylogenetic analysis uses nucleotide or amino acid sequences or other parameters, such as domain sequences and three-dimensional structure, to construct a tree to show the evolutionary relationship among different taxa (classification units) at the molecular level. Phylogenetic analysis can also be used to investigate domain relationships within an individual taxon, particularly for organisms that have undergone substantial change in morphology and physiology, but for which researchers lack fossil evidence due to the organisms' long evolutionary history or scarcity of fossilization. In this text, a detailed protocol is described for using the phylogenetic method, including amino acid sequence alignment using Clustal Omega, and subsequent phylogenetic tree construction using both Maximum Likelihood (ML) of Molecular Evolutionary Genetics Analysis (MEGA) and Bayesian Inference via MrBayes. To investigate the origin of eukaryotic Sugars Will Eventually be Exported Transporters (SWEET) genes, 228 SWEETs including 35 SWEET proteins from unicellular eukaryotes and 57 SemiSWEET proteins from prokaryotes were analyzed. Interestingly, SemiSWEETs were found in prokaryotes, but SWEETs were found in eukaryotes. Two phylogenetic trees constructed using theoretically distinct methods have consistently suggested that the first eukaryotic SWEET gene might stem from the fusion of a bacterial SemiSWEET gene and an archaeal SemiSWEET gene. It is worth noting that one should be cautious to draw a conclusion based only on phylogenetic analysis, although it is useful to explain the underlying relationship between different taxa, which is difficult or even impossible to discern through experimental means.

Human neuroimaging as a "Big Data" science

The maturation of in vivo neuroimaging has led to incredible quantities of digital information about the human brain. While much is made of the data deluge in science, neuroimaging represents the leading edge of this onslaught of "big data". A range of neuroimaging databasing approaches has streamlined the transmission, storage, and dissemination of data from such brain imaging studies. Yet few, if any, common solutions exist to support the science of neuroimaging. In this article, we discuss how modern neuroimaging research represents a multifactorial and broad ranging data challenge, involving the growing size of the data being acquired; sociological and logistical sharing issues; infrastructural challenges for multi-site, multi-datatype archiving; and the means by which to explore and mine these data. As neuroimaging advances further, e.g. aging, genetics, and age-related disease, new vision is needed to manage and process this information while marshalling of these resources into novel results. Thus, "big data" can become "big" brain science.

Mutations in the Mitochondrial Citrate Carrier SLC25A1 are Associated with Impaired Neuromuscular Transmission

Background and Objective

Congenital myasthenic syndromes are rare inherited disorders characterized by fatigable weakness caused by malfunction of the neuromuscular junction. We performed whole exome sequencing to unravel the genetic aetiology in an English sib pair with clinical features suggestive of congenital myasthenia.

Methods

We used homozygosity mapping and whole exome sequencing to identify the candidate gene variants. Mutant protein expression and function were assessed in vitro and a knockdown zebrafish model was generated to assess neuromuscular junction development.

Results

We identified a novel homozygous missense mutation in the SLC25A1 gene, encoding the mitochondrial citrate carrier. Mutant SLC25A1 showed abnormal carrier function. SLC25A1 has recently been linked to a severe, often lethal clinical phenotype. Our patients had a milder phenotype presenting primarily as a neuromuscular (NMJ) junction defect. Of note, a previously reported patient with different compound heterozygous missense mutations of SLC25A1 has since been shown to suffer from a neuromuscular transmission defect. Using knockdown of SLC25A1 expression in zebrafish, we were able to mirror the human disease in terms of variable brain, eye and cardiac involvement. Importantly, we show clear abnormalities in the neuromuscular junction, regardless of the severity of the phenotype.

Conclusions

Based on the axonal outgrowth defects seen in SLC25A1 knockdown zebrafish, we hypothesize that the neuromuscular junction impairment may be related to pre-synaptic nerve terminal abnormalities. Our findings highlight the complex machinery required to ensure efficient neuromuscular function, beyond the proteomes exclusive to the neuromuscular synapse.

Characterization of mitochondrial dicarboxylate/tricarboxylate transporters from grape berries

Grape berries (Vitis vinifera L fruit) exhibit a double-sigmoid pattern of development that results from two successive periods of vacuolar swelling during which the nature of accumulated solutes changes significantly. Throughout the first period, called green or herbaceous stage, berries accumulate high levels of organic acids, mainly malate and tartrate. At the cellular level fruit acidity comprises both metabolism and vacuolar storage. Malic acid compartmentation is critical for optimal functioning of cytosolic enzymes. Therefore, the identification and characterization of the carriers involved in malate transport across sub-cellular compartments is of great importance. The decrease in acid content during grape berry ripening has been mainly associated to mitochondrial malate oxidation. However, no Vitis vinifera mitochondrial carrier involved in malate transport has been reported to date. Here we describe the identification of three V. vinifera mitochondrial dicarboxylate/tricarboxylate carriers (VvDTC1-3) putatively involved in mitochondrial malate, citrate and other di/tricarboxylates transport. The three VvDTCs are very similar, sharing a percentage of identical residues of at least 83 %. Expression analysis of the encoding VvDTC genes in grape berries shows that they are differentially regulated exhibiting a developmental pattern of expression. The simultaneous high expression of both VvDTC2 and VvDTC3 in grape berry mesocarp close to the onset of ripening suggests that these carriers might be involved in the transport of malate into mitochondria.

Bartonella genotypes in fleas (insecta: siphonaptera) collected from rodents in the negev desert, Israel

Fleas collected from rodents in the Negev Desert in southern Israel were molecularly screened for Bartonella species. A total of 1,148 fleas, collected from 122 rodents belonging to six species, were pooled in 245 pools based on flea species, sex, and rodent host species. Two Bartonella gene fragments, corresponding to RNA polymerase B (rpoB) and citrate synthase (gltA), were targeted, and 94 and 74 flea pools were found positive by PCR, respectively. The Bartonella 16S-23S internal transcribed spacer (ITS) region was also targeted, and 66 flea pools were found to be positive by PCR. Sixteen different Bartonella gltA genotypes were detected in 94 positive flea pools collected from 5 different rodent species, indicating that fleas collected from each rodent species can harbor several Bartonella genotypes. Based on gltA analysis, identified Bartonella genotypes were highly similar or identical to strains previously detected in rodent species from different parts of the world. A gltA fragment 100% similar to Bartonella henselae was detected in one flea pool. Another 2 flea pools contained gltA fragments that were closely related to B. henselae (98% similarity). The high sequence similarities to the zoonotic pathogen B. henselae warrant further investigation.

Multiple genotypes of influenza B viruses cocirculated in Taiwan in 2004 and 2005

An influenza B outbreak occurred in Taiwan in 2004 and 2005, during which both Victoria (Vic) and Yamagata (Ya) lineages cocirculated. This study examined 36 influenza B viral genomes isolated during the outbreak to reveal their reassortment patterns. According to the isolate groupings in phylogenetic analysis, we were able to categorize those 36 isolates as being of either the Victoria or Yamagata lineage for all eight influenza B virus genomic segments, except for the NS gene, in which clades A and B existed. Based on these groupings, three genome patterns clearly emerged, namely, pattern I (Vic+Vic+Ya+Vic+Ya+Ya+Ya+A, from segments 1 to 8), pattern II (Ya+Ya+Ya+Ya+Ya+Ya+Ya+B), and pattern III (Ya+Ya+Ya+Ya+Ya+Ya+Ya+A). According to the timeline of those isolates under investigation, it appears that pattern I and II viruses could have generated pattern III via reassortment of the NS gene. On the other hand, a genomewide comparison of all six pattern III Taiwanese viruses with 37 international influenza B viral genomes showed that two international strains, B/Oslo/71/04 and B/England/23/04, were consistently clustered with the pattern III viruses isolated in Taiwan in 2004 and 2005, suggesting that Taiwanese pattern III viruses might also have been imported due to their matching genomic composition.

Interactions between Nef and AIP1 proliferate multivesicular bodies and facilitate egress of HIV-1

Background

Nef is an accessory protein of primate lentiviruses, HIV-1, HIV-2 and SIV. Besides removing CD4 and MHC class I from the surface and activating cellular signaling cascades, Nef also binds GagPol during late stages of the viral replicative cycle. In this report, we investigated further the ability of Nef to facilitate the replication of HIV-1.

Results

To this end, first the release of new viral particles was much lower in the absence of Nef in a T cell line. Since the same results were obtained in the absence of the viral envelope using pseudo-typed viruses, this phenomenon was independent of CD4 and enhanced infectivity. Next, we found that Nef not only possesses a consensus motif for but also binds AIP1 in vitro and in vivo. AIP1 is the critical intermediate in the formation of multivesicular bodies (MVBs), which play an important role in the budding and release of viruses from infected cells. Indeed, Nef proliferated MVBs in cells, but only when its AIP1-binding site was intact. Finally, these functions of Nef were reproduced in primary macrophages, where the wild type but not mutant Nef proteins led to increased release of new viral particles from infected cells.

Conclusion

We conclude that by binding GagPol and AIP1, Nef not only proliferates MVBs but also contributes to the egress of viral particles from infected cells.

Laboratory-based surveillance and molecular epidemiology of influenza virus in Taiwan

A laboratory-based surveillance network of 11 clinical virological laboratories for influenza viruses was established in Taiwan under the coordination of the Center for Disease Control and Prevention (CDC), Taiwan. From October 2000 to March 2004, 3,244 influenza viruses were isolated, including 1,969 influenza A and 1,275 influenza B viruses. The influenza infections usually occurred frequently in winter in the northern hemisphere. However, the influenza seasonality in Taiwan was not clear during the four seasons under investigation. For example, the influenza A viruses peaked during the winters of 2001, 2002, and 2003. However, some isolated peaks were also found in the summer and fall (June to November) of 2001 and 2002. An unusual peak of influenza B also occurred in the summer of 2002 (June to August). Phylogenetic analysis shows that influenza A isolates from the same year were often grouped together. However, influenza B isolates from the year 2002 clustered into different groups, and the data indicate that both B/Victoria/2/87-like and B/Yamagata/16/88-like lineages of influenza B viruses were cocirculating. Sequence comparison of epidemic strains versus vaccine strains shows that many vaccine-like Taiwanese strains were circulating at least 2 years before the vaccine strains were introduced. No clear seasonality of influenza reports in Taiwan occurred in contrast to other more continental regions.

Influenza A virus PB1-F2 gene in recent Taiwanese isolates

Influenza A virus contains eight RNA segments and encodes 10 viral proteins. However, an 11th protein, called PB1-F2, was found in A/Puerto Rico/8/34 (H1N1). This novel protein is translated from an alternative open reading frame (ORF) in the PB1 gene. We analyzed the PB1 gene of 42 recent influenza A isolates in Taiwan, including 24 H1N1 and 18 H3N2 strains. One H1N1 isolate and 17 H3N2 isolates contained the entire PB1-F2 ORF of 90 residues, three amino acids (aa) longer than the PB1-F2 of A/Puerto Rico/8/34 at the C terminal. The one remaining H3N2 isolate encoded a truncated PB1-F2 with 79 residues. The other 23 H1N1 isolates contained a truncated PB1-F2 of 57 aa. Phylogenetic analysis of both the HA and the PB1 genes showed that they shared similar clustering of these Taiwanese isolates, suggesting that no obvious reassortment occurred between the two genomic segments.

The SmtB/ArsR family of metalloregulatory transcriptional repressors: Structural insights into prokaryotic metal resistance

The SmtB/ArsR family of prokaryotic metalloregulatory transcriptional repressors represses the expression of operons linked to stress-inducing concentrations of di- and multivalent heavy metal ions. Derepression results from direct binding of metal ions by these homodimeric "metal sensor" proteins. An evolutionary analysis, coupled with comparative structural and spectroscopic studies of six SmtB/ArsR family members, suggests a unifying "theme and variations" model, in which individual members have evolved distinct metal selectivity profiles by alteration of one or both of two structurally distinct metal coordination sites. These two metal sites are designated alpha3N (or alpha3) and alpha5 (or alpha5C), named for the location of the metal binding ligands within the known or predicted secondary structure of individual family members. The alpha3N/alpha3 sensors, represented by Staphylococcus aureus pI258 CadC, Listeria monocytogenes CadC and Escherichia coli ArsR, form cysteine thiolate-rich coordination complexes (S(3) or S(4)) with thiophilic heavy metal pollutants including Cd(II), Pb(II), Bi(III) and As(III) via inter-subunit coordination by ligands derived from the alpha3 helix and the N-terminal "arm" (CadCs) or from the alpha3 helix only (ArsRs). The alpha5/alpha5C sensors Synechococcus SmtB, Synechocystis ZiaR, S. aureus CzrA, and Mycobacterium tuberculosis NmtR form metal complexes with biologically required metal ions Zn(II), Co(II) and Ni(II) characterized by four or more coordination bonds to a mixture of histidine and carboxylate ligands derived from the C-terminal alpha5 helices on opposite subunits. Direct binding of metal ions to either the alpha3N or alpha5 sites leads to strong, negative allosteric regulation of repressor operator/promoter binding affinity, consistent with a simple model for derepression. We hypothesize that distinct allosteric pathways for metal sensing have co-evolved with metal specificities of distinct alpha3N and alpha5 coordination complexes.