COVID-19 liver and gastroenterology findings: An in silico analysis of SARS-CoV-2 interactions with liver molecules

BACKGROUND

Coronavirus disease 19 (COVID-19) has not only been shown to affect the respiratory system, but has also demonstrated variable clinical presentations including gastrointestinal tract disorders. In addition, abnormalities in liver enzymes have been reported indicating hepatic injury. It is known that severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) might infect cells via the viral receptor angiotensin-converting enzyme 2 (ACE2) which is expressed in several organs including the liver. The viral Spike glycoprotein binds to ACE2 and must be cleaved by Furin and Type 2 Serine Protease to enter the cells. After that, the Akt/mTOR signaling pathway is activated and several COVID-19 changes are triggered.

AIM

To analyze liver and gastrointestinal symptoms and cell signaling pathways triggered by SARS-CoV-2 infection due to virus-liver interactions in silico.

METHODS

In this in silico study, the three-dimensional structures of the Akt, mTORC1 and Furin (receptors) were selected from the Protein Data Bank (PDB) and the structures of inhibitors (ligands) MK-2206, CC-223 and Naphthofluorescein were selected from PubChem and ZINC databases. Ligand files were downloaded as 2D structures and converted to optimized 3D structures using ViewerLite 4.2 software. Marvin Sketch^® software was used to calculate prediction of the protonated form of inhibitors in a physiological environment (pH 7.4). AutoDock Tools (ADT) software was used to calculate and delimit the Grid box used in the molecular docking of each structure selected in the PDB. In addition, protonated ligands were prepared for molecular docking using ADT software. Molecular docking was performed using ADT software tools connected to Vina software. Analysis of the amino acid residues involved in ligand interactions, as well as ligand twists, the atoms involved in interactions, bond type and strength of interactions were performed using PyMol^® and Discovery Studio^® (BIOVIA) software.

RESULTS

Molecular docking analysis showed that the mTORC1/CC-223 complex had affinity energy between the receptor and ligand of -7.7 kcal/moL with interactions ranging from 2.7 to 4.99 Å. There were four significant chemical bonds which involved two of five polypeptide chains that formed the FKBP12–Rapamycin-Binding (FRB) domain. The strongest was a hydrogen bond, the only polar interaction, and Van der Waals interactions shown to be present in 12 residues of mTORC1’s FRB domain. With regard to the Akt/MK-2206 complex there were three Van der Waals interactions and 12 chemical bonds in which seven residues of Akt were involved with all five rings of the MK-2206 structure. In this way, both ASP 388 and GLN 391 bind to the same MK-2206 ring, the smaller one. However, LYS 386 had four chemical bonds with the inhibitor, one with each structure ring, while LYS 387 binds two distinct rings. One of the MK-2206 inhibitor's rings which binds to LYS 387 also binds simultaneously to ILE 367 and LEU 385 residues, and the fifth ring of the structure was involved in a bond with the ALA 382 residue. The hydrogen bonds were the shortest bonds in the complex (2.61 and 3.08 Å) and all interactions had an affinity energy of -8.8 kcal/moL. The affinity energy in the Furin/Naphhofluorescein complex was -9.8 kcal/moL and involved six interactions ranging from 2.57 to 4.98 Å. Among them, two were polar and the others were non-polar, in addition to twelve more Van der Waals interactions. Two distinct hydrogen bonds were formed between Furin and its inhibitor involving GLN 388 and ALA 532 residues. ALA 532 also binds to two distinct rings of Naphthofluorescein, while TRP 531 residue has two simultaneous bonds with the inhibitor.

CONCLUSION

Liver infection and signaling pathways altered by SARS-CoV-2 can be modulated by inhibitors that demonstrate significant interaction affinity with human proteins, which could prevent the development of infection and symptoms.

Molecular docking of DS-3032B, a mouse double minute 2 enzyme antagonist with potential for oncology treatment development

BACKGROUND

It is known that p53 suppression is an important marker of poor prognosis of cancers, especially in solid tumors of the breast, lung, stomach, and esophagus; liposarcomas, glioblastomas, and leukemias. Because p53 has mouse double minute 2 (MDM2) as its primary negative regulator, this molecular docking study seeks to answer the following hypotheses: Is the interaction between DS-3032B and MDM2 stable enough for this drug to be considered as a promising neoplastic inhibitor?

AIM

To analyze, in silico, the chemical bonds between the antagonist DS-3032B and its binding site in MDM2.

METHODS

For molecular docking simulations, the file containing structures of MDM2 (receptor) and the drug DS-3032B (ligand) were selected. The three-dimensional structure of MDM2 was obtained from Protein Data Bank, and the one for DS-3032B was obtained from PubChem database. The location and dimensions of the Grid box was determined using AutoDock Tools software. In this case, the dimensions of the Grid encompassed the entire receptor. The ligand DS-3032B interacts with the MDM2 receptor in a physiological environment with pH 7.4; thus, to simulate more reliably, its interaction was made with the calculation for the prediction of its protonation state using the MarvinSketch^® software. Both ligands, with and without the protonation, were prepared for molecular docking using the AutoDock Tools software. This software detects the torsion points of the drug and calculates the angle of the torsions. Molecular docking simulations were performed using the tools of the AutoDock platform connected to the Vina software. The analyses of the amino acid residues involved in the interactions between the receptor and the ligand as well as the twists of the ligand, atoms involved in the interactions, and type, strength, and length of the interactions were performed using the PyMol software (pymol.org/2) and Discovery Studio from BIOVIA^®.

RESULTS

The global alignment indicated crystal structure 5SWK was more suitable for docking simulations by presenting the p53 binding site. The three-dimensional structure 5SWK for MDM2 was selected from Protein Data Bank and the three-dimensional structure of DS-3032B was selected from PubChem (Compound CID: 73297272; Milademetan). After molecular docking simulations, the most stable conformer was selected for both protonated and non-protonated DS-3032B. The interaction between MDM2 and DS-3032B occurs with high affinity; no significant difference was observed in the affinity energies between the MDM2/pronated DS-3032B (-9.9 kcal/mol) and MDM2/non-protonated DS-3032B conformers (-10.0 kcal/mol). Sixteen amino acid residues of MDM2 are involved in chemical bonds with the protonated DS-3032B; these 16 residues of MDM2 belong to the p53 biding site region and provide high affinity to interaction and stability to drug-protein complex.

CONCLUSION

Molecular docking indicated that DS-3032B antagonist binds to the same region of the p53 binding site in the MDM2 with high affinity and stability, and this suggests therapeutic efficiency.

Epitope mapping from Mycobacterium leprae proteins: Convergent data from in silico and in vitro approaches for serodiagnosis of leprosy

Knowledge of immunodominant B-cell epitopes is essential to design powerful diagnostic strategies aiming for antibody detection. Outstanding progress in computational prediction has achieved a significant contribution to the biomedical fields, including immunodiagnosis. In silico analysis may have an even more important role when information concerning antigens from etiologic agents of neglected diseases, such as leprosy, is scarce. The aim of this study was to provide mapping of B-cell epitopes from two Mycobacterium leprae-derived antigens (Ag85B and ML2055), confirm their antigenicity, and to assess the ability of in silico immunoinformatics tools to accurately predict them. Linear B-cell epitopes predicted by ABCpred and SVMTrip servers were compared to antigenic regions of synthetic overlapping peptides that exhibited reactivity to antibodies from patients with leprosy. Our in vitro results identified several immunodominant regions that had also been indicated by in silico prediction, providing agreement between experimental and simulated data. After chemical synthesis, we used enzyme-linked immunosorbent assays to determine the effectiveness of the first identified sequence (GTNVPAEFLENFVHG) which had 72 % sensitivity and 78 % specificity (AUC = 0.79) while the second one (PVSSEAQPGDPNAPS) had 72 % sensitivity and 93.8 % specificity (AUC = 0.85). Using dot blotting, an easy-to-read visual test, both peptides could distinguish sera from patients with leprosy from those with tuberculosis and from sera of healthy volunteers. Our findings suggest that these synthetic peptides, with some refinement, may be useful as serological diagnostic antigens for leprosy. In addition, it was displayed that immunoinformatics provides reliable information for mapping potential B-cell epitopes for development of peptide-based diagnostic assays for neglected diseases.

Análise da hemoglobina como modelo de estudo de interação droga-proteína in vitro e in silico

Modelo do estudo: Trata-se de um estudo experimental in vitro com abordagem computacional. Objetivo: Analisar a existência de interação entre as drogas hidrofóbicas bezafibrato e hidroclorotiazida com a hemoglobina a fim de prever alterações na biodisponibilidade das drogas, bem como na função proteica. Metodologia: Testes de interação in vitro entre a hemoglobina bovina e bezafibrato ou hidroclorotiazida foram realizados por espectrofotometria; análises dos sítios de interação e extrapolações para a hemoglobina humana foram feitas por técnicas de bioinformática. Resultados: Os testes in vitro demonstraram diminuição de absorbância (k) em 405 nm igual a 8,75 x 10-4 min-1 para o bezafibrato e 6,25 x 10-4 min--1 para a hidroclorotiazida. A diminuição sugere interação das drogas com a hemoglobina, sendo que o bezafibrato parece interagir com afinidade ligeiramente maior. As análises in silico mostraram que as drogas se ligam à porção proteica da hemoglobina. A constante de afinidade de ligação obtida por ancoragem molecular para o bezafibrato com a hemoglobina bovina (-8,3 kcal/mol) corrobora com o valor experimental de k e com o maior número de interações observadas, em relação à hidroclorotiazida (-6,6 kcal/mol). O mesmo padrão é observado para a interação do bezafibrato (-7,6 kcal/mol) e da hidroclorotiazida (-6,7 kcal/mol) com a hemoglobina humana. Conclusão: As técnicas de espectrofotometria e bioinformática utilizadas sugerem a possibilidade de interação da hemoglobina com drogas de natureza hidrofóbica, como bezafibrato e hidroclorotiazida, sendo que essa interação pode afetar a função normal da hemoglobina e alterar a farmacodinâmica e farmacocinética das drogas prejudicando sua eficiência terapêutica.

Biomolecular engineering of antidehydroepiandrosterone antibodies: a new perspective in cancer diagnosis and treatment using single-chain antibody variable fragment

Aim: To develop a monoclonal antibody against dehydroepiandrosterone (DHEA) and miniaturize it, generating a single-chain antibody variable fragment (scFv) against DHEA as an adrenocortical carcinoma (ACC) marker. Material & methods: DHEA conjugated to keyhole limpet hemocyanin was used as an immunogen to obtain anti-DHEA hybridomas. Variable fragments were cloned from hybridoma 5B7 total RNA, and used to detect DHEA in normal adrenal tissue and ACC cells. Results: IgM monoclonal antibody was highly specific, and the recombinant scFv preserved parental antibody characteristics, allowing tissue localization of DHEA. Conclusion: Undefined small lesions are challenges for clinicians and impact clinical adrenocortical tumor management. Generating an anti-DHEA scFv facilitates development of imaging tests for early diagnosis of pediatric ACC.

STAT3 polymorphism and Helicobacter pylori CagA strains with higher number of EPIYA-C segments independently increase the risk of gastric cancer

Background

Because to date there is no available study on STAT3 polymorphism and gastric cancer in Western populations and taking into account that Helicobacter pylori CagA EPIYA-C segment deregulates SHP-2/ERK-JAK/STAT3 pathways, we evaluated whether the two variables are independently associated with gastric cancer.

Methods

We included 1048 subjects: H. pylori-positive patients with gastric carcinoma (n = 232) and with gastritis (n = 275) and 541 blood donors. Data were analyzed using logistic regression model.

Results

The rs744166 polymorphic G allele (p = 0.01; OR = 1.76; 95 % CI = 1.44-2.70), and CagA-positive (OR = 12.80; 95 % CI = 5.58-19.86) status were independently associated with gastric cancer in comparison with blood donors. The rs744166 polymorphism (p = 0.001; OR = 1.64; 95 % CI = 1.16-2.31) and infection with H. pylori CagA-positive strains possessing higher number of EPIYA-C segments (p = 0.001; OR = 2.28; 95 % CI = 1.41-3.68) were independently associated with gastric cancer in comparison with gastritis. The association was stronger when host and bacterium genotypes were combined (p < 0.001; OR = 3.01; 95 % CI = 2.29-3.98). When stimulated with LPS (lipopolysaccharide) or Pam3Cys, peripheral mononuclear cells of healthy carriers of the rs744166 GG and AG genotypes expressed higher levels of STAT3 mRNA than those carrying AA genotype (p = 0.04 for both). The nuclear expression of phosphorylated p-STAT3 protein was significantly higher in the antral gastric tissue of carriers of rs744166 GG genotype than in carriers of AG and AA genotypes.

Conclusions

Our study provides evidence that STAT3 rs744166 G allele and infection with CagA-positive H. pylori with higher number of EPIYA-C segments are independent risk factors for gastric cancer. The odds ratio of having gastric cancer was greater when bacterium and host high risk genotypes were combined.

Th1 immune response to H. pylori infection varies according to the age of the patients and influences the gastric inflammatory patterns

To compare children and adults in respect to the effect of H. pylori infection on the gastric concentrations of cytokines linked to innate and Th1 immune response, as well as to investigate the changes in the gastric concentrations of the studied cytokines according to the age. We studied 245 children (142 H. pylori-negative and 103 H. pylori-positive) and 140 adults (40 H. pylori-negative and 100 H. pylori-positive). The gastric concentrations of cytokines representative of the innate and Th1 response were higher in the H. pylori-positive than in the -negative children and adults. The gastric concentrations of IL-1α and TNF-α were significantly higher, while those of IL-2, IL-12p70 and IFN-γ were lower in the infected children than in the infected adults. In the infected children, the gastric concentration of IL-1α, IL-2, IL-12p70 and IFN-γ increased, whereas in adults, the gastric concentrations of IFN-γ and IL-12p70 decreased with the aging. Increased gastric concentration of Th1 associated cytokines correlated with increased degree of gastritis that is the background lesion for the development of the H. pylori associated severe diseases. Concluding, Th1 response to H. pylori infection varies according to the age and seems to have determinant implication in the H. pylori infection outcomes.

Increased gastric IL-1β concentration and iron deficiency parameters in H. pylori infected children

Association between H. pylori infection, iron deficiency and iron deficiency anaemia has been described, but the mechanisms involved have not been established. We hypothesized that in H. pylori infected children increased gastric concentrations of IL-1β and/or TNF-α, both potent inhibitors of gastric acid secretion that is essential for iron absorption, are predictors for low blood concentrations of ferritin and haemoglobin, markers of early depletion of iron stores and anaemia, respectively. We evaluated 125 children undergoing endoscopy to clarify the origin of gastrointestinal symptoms. Gastric specimens were obtained for H. pylori status and cytokine evaluation and blood samples for determination of iron deficiency/iron deficiency anaemia parameters and IL1 cluster and TNFA polymorphisms that are associated with increased cytokine secretions. Higher IL-1β and TNF-α gastric concentrations were observed in H. pylori-positive (n = 47) than in -negative (n = 78) children. Multiple linear regression models revealed gastric IL-1β, but not TNF-α, as a significant predictor of low ferritin and haemoglobin concentrations; results were reproduced in young children in whom IL1RN polymorphic genotypes associated with higher gastric IL-1β expression and lower blood ferritin and haemoglobin concentrations. In conclusion, high gastric levels of IL-1β can be the link between H. pylori infection and iron deficiency/iron deficiency anaemia in childhood.

Size polymorphism in alleles of the myoglobin gene from biomphalaria mollusks

Introns are common among all eukaryotes, while only a limited number of introns are found in prokaryotes. Globin and globin-like proteins are widely distributed in nature, being found even in prokaryotes and a wide range of patterns of intron-exon have been reported in several eukaryotic globin genes. Globin genes in invertebrates show considerable variation in the positions of introns; globins can be found without introns, with only one intron or with three introns in different positions. In this work we analyzed the introns in the myoglobin gene from Biomphalaria glabrata, B. straminea and B. tenagophila. In the Biomphalaria genus, the myoglobin gene has three introns; these were amplified by PCR and analyzed by PCR-RFLP. Results showed that the size (number or nucleotides) and the nucleotide sequence of the coding gene of the myoglobin are variable in the three species. We observed the presence of size polymorphisms in intron 2 and 3; this characterizes a homozygous/heterozygous profile and it indicates the existence of two alleles which are different in size in each species of Biomphalaria. This polymorphism could be explored for specific identification of Biomphalaria individuals.

Analysis of the oxidase activity induced by CCl(4) and H(2)O(2) in different recombinant myoglobins

Hemoproteins may present several functions due to their prosthetic groups. After a long time, well-studied proteins such as myoglobin have surprised us with new functions. Myoglobin is a hemoprotein which has some well described and unexpected functions within the organism. Oxidase activity in standard myoglobins has been described and this activity was attributed to a covalent linkage between heme and some amino acid residues such as histidine, when myoglobins are treated with alkyl halides, and tyrosine, and when myoglobins are treated with H(2)O(2). We have found that the oxidase activity, due to H(2)O(2) treatment, can appear in different myoglobins, which presents no key residue, such as Tyr 103, for the oxidase activity previously described in the literature.

Immune responses induced by the Leishmania (Leishmania) donovani A2 antigen, but not by the LACK antigen, are protective against experimental Leishmania (Leishmania) amazonensis infection

Leishmania amazonensis is one of the major etiologic agents of a broad spectrum of clinical forms of leishmaniasis and has a wide geographical distribution in the Americas, which overlaps with the areas of transmission of many other Leishmania species. The LACK and A2 antigens are shared by various Leishmania species. A2 was previously shown to induce a potent Th1 immune response and protection against L. donovani infection in BALB/c mice. LACK is effective against L. major infection, but no significant protection against L. donovani infection was observed, in spite of the induction of a potent Th1 immune response. In an attempt to select candidate antigens for an American leishmaniasis vaccine, we investigated the protective effect of these recombinant antigens (rLACK and rA2) and recombinant interleukin-12 (rIL-12) against L. amazonensis infection in BALB/c mice. As expected, immunization with either rA2-rIL-12 or rLACK-rIL-12 induced a robust Th1 response prior to infection. However, only the BALB/c mice immunized with rA2-rIL-12 were protected against infection. Sustained gamma interferon (IFN-gamma) production, high levels of anti-A2 antibodies, and low levels of parasite-specific antibodies were detected in these mice after infection. In contrast, mice immunized with rLACK-rIL-12 displayed decreased levels of IFN-gamma and high levels of both anti-LACK and parasite-specific antibodies. Curiously, the association between rA2 and rLACK antigens in the same vaccine completely inhibited the rA2-specific IFN-gamma and humoral responses and, consequently, the protective effect of the rA2 antigen against L. amazonensis infection. We concluded that A2, but not LACK, fits the requirements for a safe vaccine against American leishmaniasis.