Molecular modeling of antibodies for the treatment of TNFα-related immunological diseases

Combined in silico/in vitro approaches for identifying modulators of the activity of the p.Tyr110Cys Carnitine O-Acetyltransferase (CRAT) variant associated to an early onset case of Leigh syndrome

Carnitine O-acetyltransferase (CRAT) is a crucial enzyme involved in mitochondrial energy metabolism. Alterations in CRAT activity have emerged as significant contributors to the pathogenesis of Leigh syndrome and related mitochondrial disorders. In this study we employed an integrated approach combining in silico docking analysis and virtual screening of chemical libraries with subsequent in vitro validation to identify small molecule modulators of the activity of the wild type (WT) CRAT and the p.Tyr110Cys (Y110C) variant associated to an early onset case of Leigh syndrome. Through 3D molecular modeling, docking simulations, and virtual screening of chemical libraries, potential CRAT modulators were prioritized based on their predicted binding affinities and interactions with the 3D models of the WT-CRAT and of the p.Tyr110Cys-CRAT mutant. The performed in silico analyses were validated through in vitro assays on the purified recombinant CRAT proteins and cell-lysates from control fibroblasts and the fibroblasts of a patient with genetic diagnosis of CRAT-deficiency, carrying the compound heterozygous missense variants in the CRAT gene, namely p.Tyr110Cys and p.Val569Met. Based on the above screening by applying the indicated filtering strategy and mentioned criteria, 3 commercially available approved drugs (also known for their possible interactions with mitochondria) namely glimepiride, artemisinin and dorzolamide, as well as suramin (already known for its ability to interact with mitochondrial proteins) were tested in in vitro assays. We found that suramin (1-1000 μM) dose-dependently inhibited the activity of both WT-CRAT and p.Tyr110Cys_CRAT variant. Artemisinin (0.1-200 μM) dose-dependently stimulated the activity of the recombinant p.Tyr110Cys CRAT mutant, whereas glimepiride and dorzolamide did not change the activity of these proteins towards acetyl-CoA. This study demonstrates the effectiveness of this combined approach in identifying novel compounds for modulating CRAT enzyme activity, providing valuable insights for potential therapeutic interventions targeting CRAT-related disorders.

Insights into Transient Dimerisation of Carnitine/Acylcarnitine Carrier (SLC25A20) from Sarkosyl/PAGE, Cross-Linking Reagents, and Comparative Modelling Analysis

The carnitine/acylcarnitine carrier (CAC) is a crucial protein for cellular energy metabolism, facilitating the exchange of acylcarnitines and free carnitine across the mitochondrial membrane, thereby enabling fatty acid β-oxidation and oxidative phosphorylation (OXPHOS). Although CAC has not been crystallised, structural insights are derived from the mitochondrial ADP/ATP carrier (AAC) structures in both cytosolic and matrix conformations. These structures underpin a single binding centre-gated pore mechanism, a common feature among mitochondrial carrier (MC) family members. The functional implications of this mechanism are well-supported, yet the structural organization of the CAC, particularly the formation of dimeric or oligomeric assemblies, remains contentious. Recent investigations employing biochemical techniques on purified and reconstituted CAC, alongside molecular modelling based on crystallographic AAC dimeric structures, suggest that CAC can indeed form dimers. Importantly, this dimerization does not alter the transport mechanism, a phenomenon observed in various other membrane transporters across different protein families. This observation aligns with the ping-pong kinetic model, where the dimeric form potentially facilitates efficient substrate translocation without necessitating mechanistic alterations. The presented findings thus contribute to a deeper understanding of CAC's functional dynamics and its structural parallels with other MC family members.

The Community Pharmacy as a Study Center for the Epidemiological Analysis of the Population Vaccination against SARS-CoV-2: Evaluation of Vaccine Safety and Pharmaceutical Service

We conducted a monocentric observational study aimed at evaluating the vaccine safety and the pharmaceutical service provided at a community pharmacy (C.PHARM) in the Puglia Region in the period from 29 December 2021 to 12 March 2022 using data from 550 patients of various ages and sexes and with concomitant diseases. We collected anamnestic data, the number of hospitalizations, and any post-vaccination adverse reactions. Interviews using the integrated EQ5 method were also performed to evaluate the quality of the service offered and any therapy preference. As expected, the vaccines were reactogenic after the first dose in the patients with mild-moderate reactions, with younger age and female gender as risk factors. Immune-allergic reactions of a moderate-severe degree were observed in adult females. In the elderly, the vaccination was well tolerated. Comirnaty® showed a favorable O.R. < 1 vs. other vaccines. No cardiovascular events or hospitalizations were observed up to May 2023. Regional data indicate that all treatments during May 2023 were correlated with the viremia. PaxlovidTM was prescribed in 3% of the patients in our center and in 1.46% in the region, and distributed/dispensed on behalf of third parties in accordance with a novel distribution/dispensation protocol of the C.PHARM that resulted in a safe vaccination center providing appropriate patient inclusion during vaccination.

Proline/Glycine residues of the PG-levels guide conformational changes along the transport cycle in the mitochondrial carnitine/acylcarnitine carrier (SLC25A20)

Mitochondrial carnitine/acylcarnitine carrier (CAC) is a member of the mitochondrial carrier (MC) family and imports acylcarnitine into the mitochondrial matrix in exchange for carnitine, playing a pivotal role in carnitine shuttle, crucial for fatty acid oxidation. The crystallized structure of CAC has not been solved yet, however, the availability of several in vitro/in silico studies, also based on the crystallized structures of the ADP/ATP carrier in the cytosolic-conformation and in the matrix-conformation, has made possible to confirm the hypothesis of the single-binding centered-gated pore mechanism for all the members of the MC family. In addition, our recent bioinformatics analyses allowed quantifying in silico the importance of protein residues of MC substrate binding region, of those involved in the formation of the matrix and cytosolic gates, and of those belonging to the Pro/Gly (PG) levels, proposed to be crucial for the tilting/kinking/bending of the six MC transmembrane helices, funneling the substrate translocation pathway. Here we present a combined in silico/in vitro analysis employed for investigating the role played by a group of 6 proline residues and 6 glycine residues, highly conserved in CAC, belonging to MC PG-levels. Residues of the PG-levels surround the similarly located MC common substrate binding region, and were proposed to lead conformational changes and substrate translocation, following substrate binding. For our analysis, we employed 3D molecular modeling approaches, alanine scanning site-directed mutagenesis and in vitro transport assays. Our analysis reveals that P130 (H3), G268 (H6) and G220 (H5), mutated in alanine, affect severely CAC transport activity (mutant catalytic efficiency lower than 5 % compared to the wild type CAC), most likely due to their major role in triggering CAC conformational changes, following carnitine binding. Notably, P30A (H1) and G121A (H3) CAC mutants, increase the carnitine uptake up to 217 % and 112 %, respectively, compared to the wild type CAC.

3D structures inferred from cDNA clones identify the CD1D-Restricted γδ T cell receptor in dromedaries

The Camelidae species occupy an important immunological niche within the humoral as well as cell mediated immune response. Although recent studies have highlighted that the somatic hypermutation (SHM) shapes the T cell receptor gamma (TRG) and delta (TRD) repertoire in Camelus dromedarius, it is still unclear how γδ T cells use the TRG/TRD receptors and their respective variable V-GAMMA and V-DELTA domains to recognize antigen in an antibody-like fashion. Here we report about 3D structural analyses of the human and dromedary γδ T cell receptor. First, we have estimated the interaction energies at the interface within the human crystallized paired TRG/TRD chains and quantified interaction energies within the same human TRG/TRD chains in complex with the CD1D, an RPI-MH1-LIKE antigen presenting glycoprotein. Then, we used the human TRG/TRD-CD1D complex as template for the 3D structure of the dromedary TRG/TRD-CD1D complex and for guiding the 3D human/dromedary comparative analysis. The choice of mutated TRG alternatively combined with mutated TRD cDNA clones originating from the spleen of one single dromedary was crucial to quantify the strength of the interactions at the protein-protein interface between the paired C. dromedarius TRG and TRD V-domains and between the C. dromedarius TRG/TRD V-domains and CD1D G-domains. Interacting amino acids located in the V-domain Complementarity Determining Regions (CDR) and Framework Regions (FR) according to the IMGT unique numbering for V-domains were identified. The resulting 3D dromedary TRG V-GAMMA combined with TRD V-DELTA protein complexes allowed to deduce the most stable gamma/delta chains pairings and to propose a candidate CD1D-restricted γδ T cell receptor complex.

BRAFV600E;K601Q metastatic melanoma patient-derived organoids and docking analysis to predict the response to targeted therapy

The V600E mutation in BRAF is associated with increased phosphorylation of Erk1/2 and high sensitivity to BRAFi/MEKi combination in metastatic melanoma. In very few patients, a tandem mutation in BRAF, V600 and K601, causes a different response to BRAFi/MEKi combination. BRAF^V600E;K601Q patient-derived organoids (PDOs) were generated to investigate targeted therapy efficacy and docking analysis was used to assess BRAF^V600E;K601Q interactions with Vemurafenib. PDOs were not sensitive to Vemurafenib and Cobimetinib given alone and sensitive to their combination, although not as responsive as BRAF^V600E PDOs. The docking analysis justified such a result showing that the tandem mutation in BRAF reduced the affinity for Vemurafenib. Tumor analysis showed that BRAF^V600E;K601Q displayed both increased phosphorylation of Erk1/2 at cytoplasmic level and activation of Notch resistance signaling. This prompted us to inhibit Notch signaling with Nirogacestat, achieving a greater antitumor response and providing PDOs-based evaluation of treatment efficacy in such rare metastatic melanoma.

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.

The Non-Gastric H+/K+ ATPase (ATP12A) Is Expressed in Mammalian Spermatozoa

H⁺/K⁺ ATPase Type 2 is an heteromeric membrane protein involved in cation transmembrane transport and consists of two subunits: a specific α subunit (ATP12A) and a non-specific β subunit. The aim of this study was to demonstrate the presence and establish the localization of ATP12A in spermatozoa from Bubalus bubalis, Bos taurus and Ovis aries. Immunoblotting revealed, in all three species, a major band (100 kDa) corresponding to the expected molecular mass. The ATP12A immunolocalization pattern showed, consistently in the three species, a strong signal at the acrosome. These results, described here for the first time in spermatozoa, are consistent with those observed for the β₁ subunit of Na⁺/K⁺ ATPase, suggesting that the latter may assemble with the α subunit to produce a functional ATP12A dimer in sperm cells. The above scenario appeared to be nicely supported by 3D comparative modeling and interaction energy calculations. The expression of ATP12A during different stages of bovine sperm maturation progressively increased, moving from epididymis to deferent ducts. Based on overall results, we hypothesize that ATP12A may play a role in acrosome reactions. Further studies will be required in order to address the functional role of this target protein in sperm physiology.

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.

Overview of Humira® Biosimilars: Current European Landscape and Future Implications

Humira® (adalimumab) by AbbVie has been the top-selling biologic drug product for the last few years - reaching nearly $20 billion in annual sales in 2018. Upon the October 2018 release of four adalimumab biosimilars into the European market, those sales began to shrink. By the end of 2019, the annual sales of Humira®, albeit still high, dipped closer to $19 billion as nearly 35% of European patients had been switched from Humira® to a biosimilar. Diminishing sales are expected to continue as the adoption of adalimumab biosimilars increases in Europe and Humira®'s patent protection is lost in the United States come 2023. In this review we discuss how impactful the availability of biosimilars has been to the European adalimumab market approximately two years after their release. We further analyze the marketed biosimilars with regards to differences in their formulation, delivery devices, biological activity, physicochemical properties, clinical trials data, and current financial foothold. More importantly, though, we highlight how "similar" these biosimilars are to Humira®. In doing so, we seek to educate the public on what they may be able to expect once adalimumab biosimilars enter the United States market in 2023.

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.

Targeting Penicillium expansum GMC Oxidoreductase with High Affinity Small Molecules for Reducing Patulin Production

Simple Summary

With the urgent necessity of potential treatments for limiting mycotoxin production and postharvest fungal rots, we propose a combined in silico/in vitro/in vivo strategy for the rapid and effective identification of bioactive small molecules, chosen among a chemical library hosting approved drugs and phytochemicals, to be used after harvest. The molecular target of our analysis was the GMC oxidoreductase from Penicillium expansum involved in the biosynthesis of patulin, a mycotoxin that can contaminate many foods, especially fruits and fruit-based products. The employed in silico/in vitro/in vivo assays described in our study proved the effectiveness of our strategy and in particular of two small molecules, 6-hydroxycoumarin (structurally related to umbelliferon, an already characterized patulin synthase inhibitor) and meticrane (an already approved drug) in reducing patulin accumulation. Our findings highly recommend the mentioned ligands to be subjected to further analysis for being used in the next future in place of other more toxic compounds, in postharvest treatments based on dipping or drenching methods.

Abstract

Flavine adenine dinucleotide (FAD) dependent glucose methanol choline oxidoreductase (GMC oxidoreductase) is the terminal key enzyme of the patulin biosynthetic pathway. GMC oxidoreductase catalyzes the oxidative ring closure of (E)-ascladiol to patulin. Currently, no protein involved in the patulin biosynthesis in Penicillium expansum has been experimentally characterized or solved by X-ray diffraction. Consequently, nothing is known about P. expansum GMC oxidoreductase substrate-binding site and mode of action. In the present investigation, a 3D comparative model for P. expansum GMC oxidoreductase has been described. Furthermore, a multistep computational approach was used to identify P. expansum GMC oxidoreductase residues involved in the FAD binding and in substrate recognition. Notably, the obtained 3D comparative model of P. expansum GMC oxidoreductase was used for performing a virtual screening of a chemical/drug library, which allowed to predict new GMC oxidoreductase high affinity ligands to be tested in in vitro/in vivo assays. In vitro assays performed in presence of 6-hydroxycoumarin and meticrane, among the highly affinity predicted binders, confirmed a dose-dependent inhibition (17–81%) of patulin production by 6-hydroxycoumarin (10 µM–1 mM concentration range), whereas the approved drug meticrane inhibited patulin production by 43% already at 10 µM. Furthermore, 6-hydroxycoumarin and meticrane caused a 60 and 41% reduction of patulin production, respectively, in vivo on apples at 100 µg/wound.

Quaternary variations in the structural assembly of N-acetylglucosamine-6-phosphate deacetylase from Pasteurella multocida

N-acetylglucosamine 6-phosphate deacetylase (NagA) catalyzes the conversion of N-acetylglucosamine-6-phosphate to glucosamine-6-phosphate in amino sugar catabolism. This conversion is an essential step in the catabolism of sialic acid in several pathogenic bacteria, including Pasteurella multocida, and thus NagA is identified as a potential drug target. Here, we report the unique structural features of NagA from P. multocida (PmNagA) resolved to 1.95 Å. PmNagA displays an altered quaternary architecture with unique interface interactions compared to its close homolog, the Escherichia coli NagA (EcNagA). We confirmed that the altered quaternary structure is not a crystallographic artifact using single particle electron cryo-microscopy. Analysis of the determined crystal structure reveals a set of hot-spot residues involved in novel interactions at the dimer-dimer interface. PmNagA binds to one Zn2+ ion in the active site and demonstrates kinetic parameters comparable to other bacterial homologs. Kinetic studies reveal that at high substrate concentrations (~10-fold the KM ), the tetrameric PmNagA displays hysteresis similar to its distant neighbor, the dimeric Staphylococcus aureus NagA (SaNagA). Our findings provide key information on structural and functional properties of NagA in P. multocida that could be utilized to design novel antibacterials.

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.

FAD/NADH Dependent Oxidoreductases: From Different Amino Acid Sequences to Similar Protein Shapes for Playing an Ancient Function

Flavoprotein oxidoreductases are members of a large protein family of specialized dehydrogenases, which include type II NADH dehydrogenase, pyridine nucleotide-disulphide oxidoreductases, ferredoxin-NAD+ reductases, NADH oxidases, and NADH peroxidases, playing a crucial role in the metabolism of several prokaryotes and eukaryotes. Although several studies have been performed on single members or protein subgroups of flavoprotein oxidoreductases, a comprehensive analysis on structure-function relationships among the different members and subgroups of this great dehydrogenase family is still missing. Here, we present a structural comparative analysis showing that the investigated flavoprotein oxidoreductases have a highly similar overall structure, although the investigated dehydrogenases are quite different in functional annotations and global amino acid composition. The different functional annotation is ascribed to their participation in species-specific metabolic pathways based on the same biochemical reaction, i.e., the oxidation of specific cofactors, like NADH and FADH₂. Notably, the performed comparative analysis sheds light on conserved sequence features that reflect very similar oxidation mechanisms, conserved among flavoprotein oxidoreductases belonging to phylogenetically distant species, as the bacterial type II NADH dehydrogenases and the mammalian apoptosis-inducing factor protein, until now retained as unique protein entities in Bacteria/Fungi or Animals, respectively. Furthermore, the presented computational analyses will allow consideration of FAD/NADH oxidoreductases as a possible target of new small molecules to be used as modulators of mitochondrial respiration for patients affected by rare diseases or cancer showing mitochondrial dysfunction, or antibiotics for treating bacterial/fungal/protista infections.

Monoclonal Antibodies Targeting Ion Channels and Their Therapeutic Potential

Monoclonal antibodies (mAbs) represent a rapidly growing pharmaceutical class of protein drugs that becomes an important part of the precision therapy. mAbs are characterized by their high specificity and affinity for the target antigen, which is mostly present on the cell surface. Ion channels are a large family of transmembrane proteins that control ion transport across the cell membrane. They are involved in almost all biological processes in both health and disease and are widely considered as prospective targets. However, no antibody-based drug targeting ion channel has been developed so far that has progressed to clinical use. Thus, we provide a comprehensive review of the elaborated mAbs against ion channels, describe their mechanisms of action, and discuss their therapeutic potential.

Effect of cariporide on ram sperm pH regulation and motility: possible role of NHE1

Sperm motility, a feature essential for in vivo fertilization, is influenced by intracellular pH (pH_i) homeostasis. Several mechanisms are involved in pH_i regulation, among which sodium-hydrogen exchangers (NHEs), a family of integral transmembrane proteins that catalyze the exchange of Na⁺ for H⁺ across lipid bilayers. A preliminary characterization of NHE activity and kinetic parameters, followed by analysis of the expression and localization of the protein in ram spermatozoa was performed. NHE activity showed an apparent K_m for external Na⁺ of 17.61 mM. Immunoblotting revealed a molecular mass of 85 kDa. Immunolocalization pattern showed some species-specific aspects, such as positive labeling at the equatorial region of the sperm head. Cariporide, a selective NHE1 inhibitor, significantly reduced pH_i recovery (85%). Similarly, exposure to cariporide significantly inhibited different motility parameters, including those related to sperm capacitation. In vitro fertilization (IVF) was not affected by cariporide, possibly due to the non-dramatic, although significant, drop in motility and velocity parameters or due to prolonged exposure during IVF, which may have caused progressive loss of its inhibitory effect. In conclusion, this is the first study documenting, in a large animal model (sheep) of well-known translational relevance, a direct functional role of NHE on sperm pH_i and motility. The postulated specificity of cariporide toward isoform 1 of the Na⁺/H⁺ exchanger seems to suggest that NHE1 may contribute to the observed effects on sperm cell functionality.

The Immunogenicity of Branded and Biosimilar Infliximab in Rheumatoid Arthritis According to Th9-Related Responses

Our objective was to evaluate the immunogenicity of branded and biosimilar infliximab by detecting changes in T-helper-9 (Th9) percentages induced by an in vitro stimulation test.

METHODS

Peripheral blood mononuclear cells collected from 55 consecutive rheumatoid arthritis (RA) outpatients (15 drug free, 20 successfully treated with branded infliximab, 20 branded infliximab inadequate responders) and 10 healthy controls were cultured, with or without 50 μg/mL of infliximab originator (Remicade^®) or 50 μg/mL of infliximab biosimilar (Remsima^®) for 18 h. Th9 lymphocytes were identified by means of flow cytometry as PU.1 and IRF4-expressing, IL-9-secreting CD4⁺ T cells. Furthermore, the markers CCR7 and CD45RA were used to distinguish naïve from memory IL-9 producer cells.

RESULTS

Under unstimulated conditions, the drug-free RA patients had the highest percentages of Th9 lymphocytes. Following stimulation with branded infliximab, the percentages of PU.1 and IRF4-expressing Th9 cells, CCR7⁺, CD45RA^- (central memory) and CCR7^-, CD45RA^- (effector memory) cells significantly increased in the group of inadequate responders, but no significant variation was observed after exposure to the biosimilar of infliximab.

CONCLUSIONS

Th9 cells seem to be involved in the immune response to the epitopes of branded, but not biosimilar, infliximab, and this may depend on the recall and stimulation of both central and effector memory cells.

The process defines the product: what really matters in biosimilar design and production?

Biologic drugs are highly complex molecules produced by living cells through a multistep manufacturing process. The key characteristics of these molecules, known as critical quality attributes (CQAs), can vary based on post-translational modifications that occur in the cellular environment or during the manufacturing process. The extent of the variation in each of the CQAs must be characterized for the originator molecule and systematically matched as closely as possible by the biosimilar developer to ensure bio-similarity. The close matching of the originator fingerprint is the foundation of the biosimilarity exercise, as the analytical tools designed to measure differences at the molecular level are far more sensitive and specific than tools available to physicians during clinical trials. Biosimilar development, therefore, has a greater focus on preclinical attributes compared with the development of an original biological agent. As changes in CQAs can occur at different stages of the manufacturing process, even small modifications to the process can alter biosimilar attributes beyond the point of similarity and impact clinical effectiveness and safety. The manufacturer's ability to provide consistent production and quality control will greatly influence the acceptance of biosimilars. To this end, preventing drift from the required specifications over time and avoiding the various implications brought by product shortage will enhance biosimilar integration into daily practice. As most prescribers are not familiar with this new drug development paradigm, educational programmes will be needed so that prescribers see biosimilars as fully equivalent, efficacious and safe medicines when compared with originator products.

Insights on the use of biosimilars in the treatment of inflammatory bowel disease

Biologic therapy, such as those that target tumor necrosis factor (TNF) signaling, has proven to be an efficacious method of treatment for patients with inflammatory bowel disease (IBD) with regards to symptom management and mucosal healing. However, the rising prevalence of IBD worldwide and the ever-increasing burden of biologic pharmaceuticals in the health care industry is alarming for insurance companies, clinicians, and patients. The impending patent expiry and the relatively high costs of biologics, particularly anti-TNF agents, have paved the way for biosimilar development for IBD. The United States Food and Drug Administration defines a biosimilar as a biological product that is highly similar to its reference medicinal product, with no clinically meaningful differences in terms of safety, purity, and potency. The hope with biosimilars is that their entry into the market will be able to drive competition between pharmaceutical companies to reduce prices like that of the generic market, and that access to appropriate biologic treatments for IBD patients is increased in the long-term. Yet, there are challenging issues such as indication extrapolation and interchangeability that are still being debated in the field of IBD and must be addressed in future issued guidance. This review will discuss the issues and implications concerning the use of biosimilar therapy for IBD.

Molecular modeling of antibodies for the treatment of TNFα-related immunological diseases

Therapeutic monoclonal antibodies (mAbs) have high efficacy in treating TNF α-related immunological diseases. Other than neutralizing TNF α, these IgG1 antibodies exert Fc receptor-mediated effector functions such as the complement-dependent cytotoxicity (CDC) and antibody-dependent cell cytotoxicity (ADCC). The crystallizable fragment (Fc) of these IgG1 contains a single glycosylation site at Asn 297/300 that is essential for the CDC and ADCC. Glycosylated antibodies lacking core fucosylation showed an improved ADCC. However, no structural data are available concerning the ligand-binding interaction of these mAbs used in TNF α-related diseases and the role of the fucosylation. We therefore used comparative modeling for generating complete 3D mAb models that include the antigen-binding fragment (Fab) portions of infliximab, complexed with TNF α (4G3Y.pdb), the Fc region of the human IGHG1 fucosylated (3SGJ) and afucosylated (3SGK) complexed with the Fc receptor subtype Fcγ RIIIA, and the Fc region of a murine immunoglobulin (1IGT). After few thousand steps of energy minimization on the resulting 3D mAb models, minimized final models were used to quantify interactions occurring between Fcγ RIIIA and the fucosylated/afucosylated Fc fragments. While fucosylation does not affect Fab-TNF α interactions, we found that in the absence of fucosylation the Fc-mAb domain and Fcγ RIIIA are closer and new strong interactions are established between G129 of the receptor and S301 of the Chimera 2 Fc mAb; new polar interactions are also established between the Chimera 2 Fc residues Y299, N300, and S301 and the Fcγ RIIIA residues K128, G129, R130, and R155. These data help to explain the reduced ADCC observed in the fucosylated mAbs suggesting the specific AA residues involved in binding interactions.