ZDOCK server: interactive docking prediction of protein-protein complexes and symmetric multimers

Minimal Collagen-Binding Epitope of Glycoprotein VI in Human and Mouse Platelets

Glycoprotein VI (GPVI) is a platelet-specific receptor for collagen and fibrin, regulating important platelet functions such as platelet adhesion and thrombus growth. Although the blockade of GPVI function is widely recognized as a potent anti-thrombotic approach, there are limited studies focused on site-specific targeting of GPVI. Using computational modeling and bioinformatics, we analyzed collagen- and CRP-binding surfaces of GPVI monomers and dimers, and compared the interacting surfaces with other mammalian GPVI isoforms. We could predict a minimal collagen-binding epitope of GPVI dimer and designed an EA-20 antibody that recognizes a linear epitope of this surface. Using platelets and whole blood samples donated from wild-type and humanized GPVI transgenic mice and also humans, our experimental results show that the EA-20 antibody inhibits platelet adhesion and aggregation in response to collagen and CRP, but not to fibrin. The EA-20 antibody also prevents thrombus formation in whole blood, on the collagen-coated surface, in arterial flow conditions. We also show that EA-20 does not influence GPVI clustering or receptor shedding. Therefore, we propose that blockade of this minimal collagen-binding epitope of GPVI with the EA-20 antibody could represent a new anti-thrombotic approach by inhibiting specific interactions between GPVI and the collagen matrix.

The armadillo-repeat domain of Plakophilin 1 binds to human enzyme PADI4

Plakophilin 1 (PKP1), a member of the armadillo repeat family of proteins, is a key structural component of cell-cell adhesion scaffolds, although it can also be found in other cell locations, including the cytoplasm and the nucleus. PADI4 (peptidyl-arginine deiminase 4) is one of the human isoforms of a family of enzymes engaged in the conversion of arginine to citrulline, and is present in monocytes, macrophages, granulocytes, and in several types of cancer cells. It is the only family member observed both within the nucleus and the cytoplasm under ordinary conditions. We studied the binding of the armadillo domain of PKP1 (ARM-PKP1) with PADI4, by using several biophysical methods, namely fluorescence, far-ultraviolet (far-UV) circular dichroism (CD), isothermal titration calorimetry (ITC), and molecular simulations; furthermore, binding was also tested by Western-blot (WB) analyses. Our results show that there was binding between the two proteins, with a dissociation constant in the low micromolar range (∼ 1 μM). Molecular modelling provided additional information on the possible structure of the binding complex, and especially on the binding hot-spot predicted for PADI4. This is the first time that the interaction between these two proteins has been described and studied. Our findings could be of importance to understand the development of tumors, where PKP1 and PADI4 are involved. Moreover, our findings pave the way to describe the formation of neutrophil extracellular traps (NETs), whose construction is modulated by PADI4, and which mediate the proteolysis of cell-cell junctions where PKP1 intervenes.

Design and Prediction of Aptamers Assisted by In Silico Methods

An aptamer is a single-stranded DNA or RNA that binds to a specific target with high binding affinity. Aptamers are developed through the process of systematic evolution of ligands by exponential enrichment (SELEX), which is repeated to increase the binding power and specificity. However, the SELEX process is time-consuming, and the characterization of aptamer candidates selected through it requires additional effort. Here, we describe in silico methods in order to suggest the most efficient way to develop aptamers and minimize the laborious effort required to screen and optimise aptamers. We investigated several methods for the estimation of aptamer-target molecule binding through conformational structure prediction, molecular docking, and molecular dynamic simulation. In addition, examples of machine learning and deep learning technologies used to predict the binding of targets and ligands in the development of new drugs are introduced. This review will be helpful in the development and application of in silico aptamer screening and characterization.

FK506-Binding Protein 2 Participates in Proinsulin Folding

Apart from chaperoning, disulfide bond formation, and downstream processing, the molecular sequence of proinsulin folding is not completely understood. Proinsulin requires proline isomerization for correct folding. Since FK506-binding protein 2 (FKBP2) is an ER-resident proline isomerase, we hypothesized that FKBP2 contributes to proinsulin folding. We found that FKBP2 co-immunoprecipitated with proinsulin and its chaperone GRP94 and that inhibition of FKBP2 expression increased proinsulin turnover with reduced intracellular proinsulin and insulin levels. This phenotype was accompanied by an increased proinsulin secretion and the formation of proinsulin high-molecular-weight complexes, a sign of proinsulin misfolding. FKBP2 knockout in pancreatic β-cells increased apoptosis without detectable up-regulation of ER stress response genes. Interestingly, FKBP2 mRNA was overexpressed in β-cells from pancreatic islets of T2D patients. Based on molecular modeling and an in vitro enzymatic assay, we suggest that proline at position 28 of the proinsulin B-chain (P28) is the substrate of FKBP2's isomerization activity. We propose that this isomerization step catalyzed by FKBP2 is an essential sequence required for correct proinsulin folding.

Mutation in the D1 domain of CD4 receptor modulates the binding affinity to HIV-1 gp120

The gp120 surface subunit of HIV-1 envelope lycoprotein (Env) is the key component for the viral entry process through interaction with the CD4 binding site (CD4bs) of the primary receptor CD4. The point mutant was introduced into SD1, a CD4 D1 variant, to enhance the interaction with HIV-1 gp120.The three-dimensional structures of gp120 and SD1 were determined using homology modeling based on the results previously determined by X-ray crystallography. The binding models were carried out via protein-protein docking tools. The 5 best docking solutions were retained according to the docking scores and were used for structural assessment. Our results demonstrated the consistency between the 3D models of gp120 and SD1 predicted by molecular docking calculations and the co-crystallized data available. We first discovered that most residues in SD1 that interacted with gp120 were located within the region 6-94 of the first N-terminal D1 domain of CD4. SD1 bound to gp120 stably at which 15 residues formed 20 hydrogen bonds with 16 residues of gp120. Five pairs of electrostatic interactions between positively and negatively charged side chains of amino acids were identified in the SD1-gp120 interface, which showed an increased number of electrostatic interactions with gp120. The mutant in the D1 domain of human CD4 receptor could strengthen binding affinity with HIV-1 gp120 and might improve the interaction pattern of the neighboring residues. The sequence analysis of gp120 suggested that Asp186, Asn189, Arg191, Glu293, Phe318 and Tyr319 were located in the variable regions of gp120, which may be HIV-1 AE strain-specific amino acid residues. Together, the results presented in this study contributed to a better understanding of the changes in the interaction between the gp120 protein and the human host CD4 receptor associated with point mutation in the D1 domain. The stabilized derivative of human CD4 D1 should serve as a promising target for therapeutics development in HIV-1 vaccine and viral entry inhibitor and may warrant further investigation.

KLF14 regulates the growth of hepatocellular carcinoma cells via its modulation of iron homeostasis through the repression of iron-responsive element-binding protein 2

BACKGROUND

Hepatocellular carcinoma (HCC) is a multifactor-driven malignant tumor with rapid progression, which causes the difficulty to substantially improve the prognosis of HCC. Limited understanding of the mechanisms in HCC impedes the development of efficacious therapies. Despite Krüpple-Like factors (KLFs) were reported to be participated in HCC pathogenesis, the function of KLF14 in HCC remains largely unexplored.

METHODS

We generated KLF14 overexpressed and silenced liver cancer cells, and nude mouse xenograft models for the in vitro and in vivo study. Luciferase reporter assay, ChIP-qPCR, Co-IP, immunofluorescence were performed for mechanism research. The expression of KLF14 in HCC samples was analyzed by quantitative RT-PCR, Western blotting, and immunohistochemistry (IHC) analysis.

RESULTS

KLF14 was significantly downregulated in human HCC tissues, which was highly correlated with poor prognosis. Inhibition of KLF14 promoted liver cancer cells proliferation and overexpression of KLF14 suppressed cells growth. KLF14 exerts its anti-tumor function by inhibiting Iron-responsive element-binding protein 2 (IRP2), which then causes transferrin receptor-1(TfR1) downregulation and ferritin upregulation on the basis of IRP-IREs system. This then leading to cellular iron deficiency and HCC cells growth suppression in vitro and in vivo. Interestingly, KLF14 suppressed the transcription of IRP2 via recruiting SIRT1 to reduce the histone acetylation of the IRP2 promoter, resulting in iron depletion and cell growth suppression. More important, we found fluphenazine is an activator of KLF14, inhibiting HCC cells growth through inducing iron deficiency.

CONCLUSION

KLF14 acts as a tumor suppressor which inhibits the proliferation of HCC cells by modulating cellular iron metabolism via the repression of IRP2. We identified Fluphenazine, as an activator of KLF14, could be a potential compound for HCC therapy. Our findings therefore provide an innovative insight into the pathogenesis of HCC and a promising therapeutic target.

TLR4 is one of the receptors for Chikungunya virus envelope protein E2 and regulates virus induced pro-inflammatory responses in host macrophages

Toll like receptor 4 (TLR4), a pathogen-associated molecular pattern (PAMP) receptor, is known to exert inflammation in various cases of microbial infection, cancer and autoimmune disorders. However, any such involvement of TLR4 in Chikungunya virus (CHIKV) infection is yet to be explored. Accordingly, the role of TLR4 was investigated towards CHIKV infection and modulation of host immune responses in the current study using mice macrophage cell line RAW264.7, primary macrophage cells of different origins and in vivo mice model. The findings suggest that TLR4 inhibition using TAK-242 (a specific pharmacological inhibitor) reduces viral copy number as well as reduces the CHIKV-E2 protein level significantly using p38 and JNK-MAPK pathways. Moreover, this led to reduced expression of macrophage activation markers like CD14, CD86, MHC-II and pro-inflammatory cytokines (TNF, IL-6, MCP-1) significantly in both the mouse primary macrophages and RAW264.7 cell line, in vitro. Additionally, TAK-242-directed TLR4 inhibition demonstrated a significant reduction of percent E2-positive cells, viral titre and TNF expression in hPBMC-derived macrophages, in vitro. These observations were further validated in TLR4-knockout (KO) RAW cells. Furthermore, the interaction between CHIKV-E2 and TLR4 was demonstrated by immuno-precipitation studies, in vitro and supported by molecular docking analysis, in silico. TLR4-dependent viral entry was further validated by an anti-TLR4 antibody-mediated blocking experiment. It was noticed that TLR4 is necessary for the early events of viral infection, especially during the attachment and entry stages. Interestingly, it was also observed that TLR4 is not involved in the post-entry stages of CHIKV infection in host macrophages. The administration of TAK-242 decreased CHIKV infection significantly by reducing disease manifestations, improving survivability (around 75%) and reducing inflammation in mice model. Collectively, for the first time, this study reports TLR4 as one of the novel receptors to facilitate the attachment and entry of CHIKV in host macrophages, the TLR4-CHIKV-E2 interactions are essential for efficient viral entry and modulation of infection-induced pro-inflammatory responses in host macrophages, which might have translational implication for designing future therapeutics to regulate the CHIKV infection.

Paragraph-antibody paratope prediction using graph neural networks with minimal feature vectors

SUMMARY

The development of new vaccines and antibody therapeutics typically takes several years and requires over $1bn in investment. Accurate knowledge of the paratope (antibody binding site) can speed up and reduce the cost of this process by improving our understanding of antibody-antigen binding. We present Paragraph, a structure-based paratope prediction tool that outperforms current state-of-the-art tools using simpler feature vectors and no antigen information.

AVAILABILITY AND IMPLEMENTATION

Source code is freely available at www.github.com/oxpig/Paragraph.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Soluble domains of cytochrome c-556 and Rieske iron-sulfur protein from Chlorobaculum tepidum: Crystal structures and interaction analysis

In photosynthetic green sulfur bacteria, the electron transfer reaction from menaquinol:cytochrome c oxidoreductase to the P840 reaction center (RC) complex occurs directly without any involvement of soluble electron carrier protein(s). X-ray crystallography has determined the three-dimensional structures of the soluble domains of the CT0073 gene product and Rieske iron-sulfur protein (ISP). The former is a mono-heme cytochrome c with an α-absorption peak at 556 nm. The overall fold of the soluble domain of cytochrome c-556 (designated as cyt c-556_sol) consists of four α-helices and is very similar to that of water-soluble cyt c-554 that independently functions as an electron donor to the P840 RC complex. However, the latter's remarkably long and flexible loop between the α3 and α4 helices seems to make it impossible to be a substitute for the former. The structure of the soluble domain of the Rieske ISP (Rieske_sol protein) shows a typical β-sheets-dominated fold with a small cluster-binding and a large subdomain. The architecture of the Rieske_sol protein is bilobal and belongs to those of b₆f-type Rieske ISPs. Nuclear magnetic resonance (NMR) measurements revealed weak non-polar but specific interaction sites on Rieske_sol protein when mixed with cyt c-556_sol. Therefore, menaquinol:cytochrome c oxidoreductase in green sulfur bacteria features a Rieske/cytb complex tightly associated with membrane-anchored cyt c-556.

Human small-intestinal gluten-degrading bacteria and its potential implication in celiac disease

Celiac disease (CeD) is an immune-mediated chronic disorder triggered by the ingestion of wheat gluten in genetically predisposed individuals. Gluten is a major food ingredient, infamously containing proline and glutamine-rich domains that are highly resistant to digestion by mammalian proteolytic enzymes. Thus, adhering to a gluten-free diet (GFD) is the only known treatment for CeD, albeit with many complications. Therefore, any therapy that eliminates the gluten immunogenic part before it reaches the small intestine is highly desirable. Probiotic therapy containing gluten-degrading bacteria (GDB) and their protease enzymes are possibly new approaches to treating CeD. Our study aimed to identify novel GDB from the duodenal biopsy of the first-degree relative (FDR) subjects (relatives of diseased individuals who are healthy but susceptible to celiac disease) with the potential to reduce gluten immunogenicity. Using the gluten agar plate technique, bacterial strains Brevibacterium casei NAB46 and Staphylococcus arlettae R2AA77 displaying glutenase activity were screened, identified, and characterized. Whole-genome sequencing found gluten-degrading prolyl endopeptidase (PEP) in the B. casei NAB46 genome and glutamyl endopeptidase (GEP) in the S. arlettae R2AA77 genome. Partially purified PEP has a specific activity of 1.15 U/mg, while GEP has a specific activity of 0.84 U/mg, which are, respectively, 6- and 9-fold times higher after concentrating the enzymes. Our results showed that these enzymes could hydrolyse immunotoxic gliadin peptides recognized in western blot using an anti-gliadin antibody. Additionally, a docking model was proposed for representative gliadin peptide PQPQLPYPQPQLP in the active site of the enzymes, where the residues of the N-terminal peptide extensively interact with the catalytic domain of the enzymes. These bacteria and their associated glutenase enzymes efficiently neutralize gliadin immunogenic epitopes, opening possibilities for their application as a dietary supplement in treating CeD patients.

A novel shark single-domain antibody targeting OGT as a tool for detection and intracellular localization

Introduction

O-GlcNAcylation is a type of reversible post-translational modification on Ser/Thr residues of intracellular proteins in eukaryotic cells, which is generated by the sole O-GlcNAc transferase (OGT) and removed by O-GlcNAcase (OGA). Thousands of proteins, that are involved in various physiological and pathological processes, have been found to be O-GlcNAcylated. However, due to the lack of favorable tools, studies of the O-GlcNAcylation and OGT were impeded. Immunoglobulin new antigen receptor (IgNAR) derived from shark is attractive to research tools, diagnosis and therapeutics. The variable domain of IgNARs (VNARs) have several advantages, such as small size, good stability, low-cost manufacture, and peculiar paratope structure.

Methods

We obtained shark single domain antibodies targeting OGT by shark immunization, phage display library construction and panning. ELISA and BIACORE were used to assess the affinity of the antibodies to the antigen and three shark single-domain antibodies with high affinity were successfully screened. The three antibodies were assessed for intracellular function by flow cytometry and immunofluorescence co-localization.

Results

In this study, three anti-OGT VNARs (2D9, 3F7 and 4G2) were obtained by phage display panning. The affinity values were measured using surface plasmon resonance (SPR) that 2D9, 3F7 and 4G2 bound to OGT with KD values of 35.5 nM, 53.4 nM and 89.7 nM, respectively. Then, the VNARs were biotinylated and used for the detection and localization of OGT by ELISA, flow cytometry and immunofluorescence. 2D9, 3F7 and 4G2 were exhibited the EC50 values of 102.1 nM, 40.75 nM and 120.7 nM respectively. VNAR 3F7 was predicted to bind the amino acid residues of Ser375, Phe377, Cys379 and Tyr 380 on OGT.

Discussion

Our results show that shark single-domain antibodies targeting OGT can be used for in vitro detection and intracellular co-localization of OGT, providing a powerful tool for the study of OGT and O-GlcNAcylation.

A scoring function for the prediction of protein complex interfaces based on the neighborhood preferences of amino acids

Proteins often assemble into functional complexes, the structures of which are more difficult to obtain than those of the individual protein molecules. Given the structures of the subunits, it is possible to predict plausible complex models via computational methods such as molecular docking. Assessing the quality of the predicted models is crucial to obtain correct complex structures. Here, an energy-scoring function was developed based on the interfacial residues of structures in the Protein Data Bank. The statistically derived energy function (Nepre) imitates the neighborhood preferences of amino acids, including the types and relative positions of neighboring residues. Based on the preference statistics, a program iNepre was implemented and its performance was evaluated with several benchmarking decoy data sets. The results show that iNepre scores are powerful in model ranking to select the best protein complex structures.

Modeling Protein Complexes and Molecular Assemblies Using Computational Methods

Many biological molecules are assembled into supramolecular complexes that are necessary to perform functions in the cell. Better understanding and characterization of these molecular assemblies are thus essential to further elucidate molecular mechanisms and key protein-protein interactions that could be targeted to modulate the protein binding affinity or develop new binders. Experimental access to structural information on these supramolecular assemblies is often hampered by the size of these systems that make their recombinant production and characterization rather difficult. Computational methods combining both structural data, molecular modeling techniques, and sequence coevolution information can thus offer a good alternative to gain access to the structural organization of protein complexes and assemblies. Herein, we present some computational methods to predict structural models of the protein partners, to search for interacting regions using coevolution information, and to build molecular assemblies. The approach is exemplified using a case study to model the succinate-quinone oxidoreductase heterocomplex.

Genomic distribution of signal transducer and activator of transcription (STAT) family in colorectal cancer

JAK/STAT pathway has been widely acknowledged in the development of human cancers. However, the role of different phosphorylated STAT proteins translocating into nucleus in transcription activation of target genes is not fully understood. In present research, ChIP-seq was carried on to investigate the genome-wide distribution of the activated STAT1, STAT2, STAT3, STAT5 and STAT6 in colorectal cancer HCT-116 cells. Our observations indicated that the homodimers rather than heterodimers of STAT protein predominantly occupied on genomic DNA. STAT3 accounted for the largest proportion among all STAT proteins HCT-116 cells. Furthermore, the biased binding motif targeted by different STAT homodimers suggested the distinct biological functions. Here, we noticed that NR5A2 was a specific co-activator of STAT3 by DNA motif analysis. Co-IP assay determined that NR5A2 indeed interacted with STAT3 homodimer rather than other homodimers or heterodimers. NR5A2 knockdown resulted in a reduced binding affinity of STAT3 homodimer in the original regions. Taken together, we characterize the genome-wide landscape of activated STAT proteins, and reveal the differences of binding patterns as well as the target genes and associated functions between homodimer and heterodimer of STAT proteins in HCT-116 cells. We also present some new findings and possible mechanisms regarding the role of NR5A2 on STAT3 in CRC. Our findings may provide new insights into the design of STAT inhibitors to treat CRC and other diseases.

Staphylococcal virulence factor HlgB targets the endoplasmic-reticulum-resident E3 ubiquitin ligase AMFR to promote pneumonia

Staphylococcus aureus invades cells and persists intracellularly, causing persistent inflammation that is notoriously difficult to treat. Here we investigated host-pathogen interactions underlying intracellular S. aureus infection in macrophages and discovered that the endoplasmic reticulum (ER) is an important cellular compartment for intracellular S. aureus infection. Using CRISPR-Cas9 guide RNA library screening, we determined that the autocrine motility factor receptor (AMFR), an ER-resident E3 ubiquitin ligase, played an essential role in mediating intracellular S. aureus-induced inflammation. AMFR directly interacted with TAK1-binding protein 3 (TAB3) in the ER, inducing K27-linked polyubiquitination of TAB3 on lysine 649 and promoting TAK1 activation. Moreover, the virulence factor γ-haemolysin B (HIgB) of S. aureus bound to the AMFR and regulated TAB3. Our findings highlight an unknown role of AMFR in intracellular S. aureus infection-induced pneumonia and suggest that pharmacological interruption of AMFR-mediated TAB3 signalling cascades and HIgB targeting may prevent invasive staphylococci-mediated pneumonia.

"From molecular to clinic": The pivotal role of CDC42 in pathophysiology of human papilloma virus related cancers and a correlated sensitivity of afatinib

Background

Human papilloma virus (HPV)-related cancers are global health challenge. Insufficient comprehension of these cancers has impeded the development of novel therapeutic interventions. Bioinformatics empowered us to investigate these cancers from new entry points.

Methods

DNA methylation data of cervical squamous cell carcinoma (CESC) and anal squamous cell carcinoma (ASCC) were analyzed to identify the significantly altered pathways. Through analyses integrated with RNA sequencing data of genes in these pathways, genes with strongest correlation to the TNM staging of CESC was identified and their correlations with overall survival in patients were assessed. To find a potential promising drug, correlation analysis of gene expression levels and compound sensitivity was performed. In vitro experiments were conducted to validate these findings. We further performed molecular docking experiments to explain our findings.

Results

Significantly altered pathways included immune, HPV infection, oxidative stress, ferroptosis and necroptosis. 10 hub genes in these pathways (PSMD11, RB1, SAE1, TAF15, TFDP1, CORO1C, JOSD1, CDC42, KPNA2 and NUP62) were identified, in which only CDC42 high expression was statistically significantly correlated with overall survival (Hazard Ratio: 1.6, P = 0.045). Afatinib was then screened out to be tested. In vitro experiments exhibited that the expression level of CDC42 was upregulated in HaCaT/A431 cells transfected with HPV E6 and E7, and the inhibitory effect of afatinib on proliferation was enhanced after transfection. CDC42-GTPase-effector interface-EGFR-afatinib was found to be a stable complex with a highest ZDOCK score of 1264.017.

Conclusion

We identified CDC42 as a pivotal gene in the pathophysiology of HPV-related cancers. The upregulation of CDC42 could be a signal for afatinib treatment and the mechanism in which may be an increased affinity of EGFR to afatinib, inferred from a high stability in the quaternary complex of CDC42-GTPase-effector interface-EGFR-afatinib.

Bovine lactoferrin inhibits SARS-CoV-2 and SARS-CoV-1 by targeting the RdRp complex and alleviates viral infection in the hamster model

Breast milk has been found to inhibit coronavirus infection, while the key components and mechanisms are unknown. We aimed to determine the components that contribute to the antiviral effects of breastmilk and explore their potential mechanism. Lactoferrin (Lf) and milk fat globule membrane inhibit severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-related coronavirus GX_P2V and transcription- and replication-competent SARS-CoV-2 virus-like particles in vitro and block viral entry into cells. We confirmed that bovine Lf (bLf) blocked the binding between human angiotensin-converting enzyme 2 and SARS-CoV-2 spike protein by combining receptor-binding domain (RBD). Importantly, bLf inhibited RNA-dependent RNA polymerase (RdRp) activity of both SARS-CoV-2 and SARS-CoV in vitro in the nanomolar range. So far, no biological macromolecules have been reported to inhibit coronavirus RdRp. Our result indicated that bLf plays a major role in inhibiting viral replication. bLf treatment reduced viral load in lungs and tracheae and alleviated pathological damage. Our study provides evidence that bLf prevents SARS-CoV-2 infection by combining SARS-CoV-2 spike protein RBD and inhibiting coronaviruses' RdRp activity, and may be a promising candidate for the treatment of coronavirus disease 2019.

Pairwise and Multi-chain Protein Docking Enhanced Using LZerD Web Server

Interactions of proteins with other macromolecules have important structural and functional roles in the basic processes of living cells. To understand and elucidate the mechanisms of interactions, it is important to know the 3D structures of the complexes. Proteomes contain numerous protein-protein complexes, for which experimentally determined structures often do not exist. Computational techniques can be a practical alternative to obtain useful complex structure models. Here, we present a web server that provides access to the LZerD and Multi-LZerD protein docking tools, which can perform both pairwise and multi-chain docking. The web server is user-friendly, with options to visualize the distribution and structures of binding poses of top-scoring models. The LZerD web server is available at https://lzerd.kiharalab.org . This chapter dictates the algorithm and step-by-step procedure to model the monomeric structures with AttentiveDist, and also provides the detail of pairwise LZerD docking, and multi-LZerD. This also provided case studies for each of the three modules.

Potential antiviral peptides against the nucleoprotein of SARS-CoV-2

Nucleoprotein is a conserved structural protein of SARS-CoV-2, which is involved in several functions, including replication, packaging, and transcription. In this research, 21 antiviral peptides that are known to have inhibitory function against nucleoprotein in several other viruses, were screened computationally against the nucleoprotein of SARS-CoV-2. The complexes of five best performing peptides (AVP1142, AVP1145, AVP1148, AVP1150, AVP1155) with nucleoprotein were selected for subsequent screening via 5 ns molecular dynamics (MD) simulation. Two peptides, namely AVP1145 and AVP1155, came out as promising candidates and hence were selected for 200 ns MD simulation for further validation, incorporating a DMPC-based membrane environment. In the long MD simulation, both AVP1155 and AVP1145 utilized multiple residues-mainly aromatic, acidic, and nonpolar residues-as interacting points to remain in contact with the nucleoprotein and formed predominantly hydrogen bonds along with hydrophobic and electrostatic interactions. However, AVP1155 proved to be superior to AVP1145 when its complex with nucleoprotein was analyzed in terms of root-mean-square deviation, root-mean-square fluctuation, radius of gyration, solvent accessible surface area and free energy landscape. In a nutshell, the findings of this research may guide future studies in the development of selective peptide inhibitors of SARS-CoV-2 nucleoprotein.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11696-022-02514-4.

Modulation of interaction of BRCA1-RAD51 and BRCA1-AURKA protein complexes by natural metabolites using as possible therapeutic intervention toward cardiotoxic effects of cancer drugs: an in-silico approach

Breast cancer type 1 susceptibility protein (BRCA1) plays an important role in maintaining genome stability and is known to interact with several proteins involved in cellular pathways, gene transcription regulation and DNA damage response. More than 40% of inherited breast cancer cases are due to BRCA1 mutation. It is also a prognostic marker in non-small cell lung cancer patients as well as a gatekeeper of cardiac function. Interaction of mutant BRCA1 with other proteins is known to disrupt the tumor suppression mechanism. Two directly interacting proteins with BRCA1 namely, DNA repair protein RAD51 (RAD51) and Aurora kinase A (AURKA), known to regulate homologous recombination (HR) and G/M cell cycle transition, respectively, form protein complex with both wild and mutant BRCA1. To analyze the interactions, protein-protein complexes were generated for each pair of proteins. In order to combat the cardiotoxic effects of cancer drugs, pharmacokinetically screened natural metabolites derived from plant, marine and bacterial sources and along with FDA-approved cancer drugs as control, were subjected to molecular docking. Piperoleine B and dihydrocircumin were the best docked natural metabolites in both RAD51 and AURKA complexes, respectively. Molecular dynamics simulation (MDS) analysis and binding free energy calculations for the best docked natural metabolite and drug for both the mutant BRCA1 complexes suggested better stability for the natural metabolites piperolein B and dihydrocurcumin as compared to drug. Thus, both natural metabolites could be further analyzed for their role against the cardiotoxic effects of cancer drugs through wet lab experiments.Communicated by Ramaswamy H. Sarma.

Recent developments of sequence-based prediction of protein-protein interactions

The identification of protein-protein interactions (PPIs) can lead to a better understanding of cellular functions and biological processes of proteins and contribute to the design of drugs to target disease-causing PPIs. In addition, targeting host-pathogen PPIs is useful for elucidating infection mechanisms. Although several experimental methods have been used to identify PPIs, these methods can yet to draw complete PPI networks. Hence, computational techniques are increasingly required for the prediction of potential PPIs, which have never been seen experimentally. Recent high-performance sequence-based methods have contributed to the construction of PPI networks and the elucidation of pathogenetic mechanisms in specific diseases. However, the usefulness of these methods depends on the quality and quantity of training data of PPIs. In this brief review, we introduce currently available PPI databases and recent sequence-based methods for predicting PPIs. Also, we discuss key issues in this field and present future perspectives of the sequence-based PPI predictions.

Understanding the Allosteric Modulation of PTH1R by a Negative Allosteric Modulator

The parathyroid hormone type 1 receptor (PTH1R) acts as a canonical class B G protein-coupled receptor, regulating crucial functions including calcium homeostasis and bone formation. The identification and development of PTH1R non-peptide allosteric modulators have obtained widespread attention. It has been found that a negative allosteric modulator (NAM) could inhibit the activation of PTH1R, but the implied mechanism remains unclear. Herein, extensive molecular dynamics simulations together with multiple analytical approaches are utilized to unravel the mechanism of PTH1R allosteric inhibition. The results suggest that the binding of NAM destabilizes the structure of the PTH1R-PTH-spep/qpep (the C terminus of Gs/Gq proteins) complexes. Moreover, the presence of NAM weakens the binding of PTH/peps (spep and qpep) and PTH1R. The intra- and inter-molecular couplings are also weakened in PTH1R upon NAM binding. Interestingly, compared with our previous study of the positive allosteric effects induced by extracellular Ca²⁺, the enhanced correlation between the PTH and G-protein binding sites is significantly reduced by the replacement of this negative allosteric regulator. Our findings might contribute to the development of new therapeutic agents for diseases caused by the abnormal activation of PTH1R.

PTMint database of experimentally verified PTM regulation on protein-protein interaction

MOTIVATION

Post-translational modification (PTM) is an important biochemical process. which includes six most well-studied types: phosphorylation, acetylation, methylation, sumoylation, ubiquitylation and glycosylation. PTM is involved in various cell signaling pathways and biological processes. Abnormal PTM status is closely associated with severe diseases (such as cancer and neurologic diseases) by regulating protein functions, such as protein-protein interactions (PPIs). A set of databases was constructed separately for PTM sites and PPI; however, the resource of regulation for PTM on PPI is still unsolved.

RESULTS

Here, we firstly constructed a public accessible database of PTMint (PTMs that are associated with PPIs) (https://ptmint.sjtu.edu.cn/) that contains manually curated complete experimental evidence of the PTM regulation on PPIs in multiple organisms, including Homo sapiens, Arabidopsis thaliana, Caenorhabditis elegans, Drosophila melanogaster, Saccharomyces cerevisiae and Schizosaccharomyces pombe. Currently, the first version of PTMint encompassed 2477 non-redundant PTM sites in 1169 proteins affecting 2371 protein-protein pairs involving 357 diseases. Various annotations were systematically integrated, such as protein sequence, structure properties and protein complex analysis. PTMint database can help to insight into disease mechanism, disease diagnosis and drug discovery associated with PTM and PPI.

AVAILABILITY AND IMPLEMENTATION

PTMint is freely available at: https://ptmint.sjtu.edu.cn/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Protein complex prediction using Rosetta, AlphaFold, and mass spectrometry covalent labeling

Covalent labeling (CL) in combination with mass spectrometry can be used as an analytical tool to study and determine structural properties of protein-protein complexes. However, data from these experiments is sparse and does not unambiguously elucidate protein structure. Thus, computational algorithms are needed to deduce structure from the CL data. In this work, we present a hybrid method that combines models of protein complex subunits generated with AlphaFold with differential CL data via a CL-guided protein-protein docking in Rosetta. In a benchmark set, the RMSD (root-mean-square deviation) of the best-scoring models was below 3.6 Å for 5/5 complexes with inclusion of CL data, whereas the same quality was only achieved for 1/5 complexes without CL data. This study suggests that our integrated approach can successfully use data obtained from CL experiments to distinguish between nativelike and non-nativelike models.

Protein-protein interaction prediction methods: from docking-based to AI-based approaches

Protein-protein interactions (PPIs), such as protein-protein inhibitor, antibody-antigen complex, and supercomplexes play diverse and important roles in cells. Recent advances in structural analysis methods, including cryo-EM, for the determination of protein complex structures are remarkable. Nevertheless, much room remains for improvement and utilization of computational methods to predict PPIs because of the large number and great diversity of unresolved complex structures. This review introduces a wide array of computational methods, including our own, for estimating PPIs including antibody-antigen interactions, offering both historical and forward-looking perspectives.

In Silico Study of the Interactions of Anle138b Isomer, an Inhibitor of Amyloid Aggregation, with Partner Proteins

Herein, we aimed to highlight current "gaps" in the understanding of the potential interactions between the Anle138b isomer ligand, a promising agent for clinical research, and the intrinsically disordered alpha-synuclein protein. The presence of extensive unstructured areas in alpha-synuclein determines its existence in the cell of partner proteins, including the cyclophilin A chaperone, which prevents the aggregation of alpha-synuclein molecules that are destructive to cell life. Using flexible and cascaded molecular docking techniques, we aimed to expand our understanding of the molecular architecture of the protein complex between alpha-synuclein, cyclophilin A and the Anle138b isomer ligand. We demonstrated the possibility of intricate complex formation under cellular conditions and revealed that the main interactions that stabilize the complex are hydrophobic and involve hydrogen.

Interactive Flexible-Receptor Molecular Docking in Virtual Reality Using DockIT

Interactive docking enables the user to guide and control the docking of two biomolecules into a binding pose. It is of particular use when the binding site is known and is thought to be applicable to structure-based drug design (SBDD) and educating students about biomolecular interactions. For SBDD, it enables expertise and intuition to be brought to bear in the drug design process. In education, it can teach students about the most basic level of biomolecular function. Here, we introduce DockIT for virtual reality (VR) that uses a VR headset and hand-held controllers. Using the method of linear response on explicit solvent molecular dynamics simulations, DockIT can model both global and local conformational changes within the receptor due to forces of interaction with the ligand. It has real-time flexible molecular surface rendering and can show the real-time formation and breaking of hydrogen bonds, both between the ligand and receptor and within the receptor itself as it smoothly changes conformation.

Cystatin B increases autophagic flux by sustaining proteolytic activity of cathepsin B and fuels glycolysis in pancreatic cancer: CSTB orchestrates autophagy and glycolysis in PDAC

BACKGROUND

Both autophagy and glycolysis are essential for pancreatic ductal adenocarcinoma (PDAC) survival due to desmoplasia. We investigated whether targeting a hub gene which participates in both processes could be an efficient strategy for PDAC treatment.

METHODS

The expression pattern of glycolysis signatures (GS) and autophagy signatures (AS) and their correlation with cystatin B (CSTB) in PDAC were analysed. It was discovered how CSTB affected the growth, glycolysis, and autophagy of PDAC cells. We assessed competitive binding to cathepsin B (CTSB) between CSTB and cystatin C (CSTC) via immunoprecipitation (IP) and immunofluorescence (IF). Chromatin immunoprecipitation quantitative polymerase chain reaction (ChIP-qPCR) and luciferase reporter gene assays were used to unveil the mechanism underlying CSTB upregulation. The expression pattern of CSTB was examined in clinical samples and KrasG12D/+, Trp53R172H/+, Pdx1-Cre (KPC) mice.

RESULTS

GS and AS were enriched and closely associated in PDAC tissues. CSTB increased autophagic flux and provided substrates for glycolysis. CSTB knockdown attenuated the proliferation of PDAC cells and patient-derived xenografts. The liquid chromatography-tandem mass spectrometry assay indicated CSTB interacted with CTSB and contributed to the proteolytic activity of CTSB in lysosomes. IF and IP assays demonstrated that CSTB competed with CSTC to bind to CTSB. Mutation of the key sites of CSTB abolished the interaction between CSTB and CTSB. CSTB was highly expressed in PDAC due to H3K27acetylation and SP1 expression. High expression of CSTB in PDAC was observed in tissue microarray and patients' serum samples.

CONCLUSIONS

Our work demonstrated the tumorigenic roles of autophagy and glycolysis in PDAC. CSTB is a key role in orchestrating these processes to ensure energy supply of PDAC cells.

HBD-2 variants and SARS-CoV-2: New insights into inter-individual susceptibility

Background

A deep understanding of the causes of liability to SARS-CoV-2 is essential to develop new diagnostic tests and therapeutics against this serious virus in order to overcome this pandemic completely. In the light of the discovered role of antimicrobial peptides [such as human b-defensin-2 (hBD-2) and cathelicidin LL-37] in the defense against SARS-CoV-2, it became important to identify the damaging missense mutations in the genes of these molecules and study their role in the pathogenesis of COVID-19.

Methods

We conducted a comprehensive analysis with multiple in silico approaches to identify the damaging missense SNPs for hBD-2 and LL-37; moreover, we applied docking methods and molecular dynamics analysis to study the impact of the filtered mutations.

Results

The comprehensive analysis reveals the presence of three damaging SNPs in hBD-2; these SNPs were predicted to decrease the stability of hBD-2 with a damaging impact on hBD-2 structure as well. G51D and C53G mutations were located in highly conserved positions and were associated with differences in the secondary structures of hBD-2. Docking-coupled molecular dynamics simulation analysis revealed compromised binding affinity for hBD-2 SNPs towards the SARS-CoV-2 spike domain. Different protein-protein binding profiles for hBD-2 SNPs, in relation to their native form, were guided through residue-wise levels and differential adopted conformation/orientation.

Conclusions

The presented model paves the way for identifying patients prone to COVID-19 in a way that would guide the personalization of both the diagnostic and management protocols for this serious disease.

Multiple HA substitutions in highly pathogenic avian influenza H5Nx viruses contributed to the change in the NA subtype preference

Novel highly pathogenic avian influenza (HPAI) H5Nx viruses are predominantly circulating worldwide, with an increasing potential threat of an outbreak in humans. It remains largely unknown how the stably maintained HPAI H5N1 suddenly altered its neuraminidase (NA) to other NA subtypes, which resulted in the emergence and evolution of H5Nx viruses. Here, we found that a combination of four specific amino acid (AA) substitutions (S123P-T156A-D183N- S223 R) in the hemagglutinin (HA) protein consistently observed in the H5Nx markedly altered the NA preference of H5N1 viruses. These molecular changes in H5N1 impaired its fitness, particularly viral growth and the functional activities of the HA and NA proteins. Among the AA substitutions identified, the T156A substitution, which contributed to the NA shift, also dramatically altered the antigenicity of H5N1 viruses, suggesting an occurrence of antigenic drift triggered by selective pressure. Our study shows the importance of how HA and NA complement each other and that antigenic drift in HA can potentially cause a shift in the NA protein in influenza A virus evolution.

A novel protein encoded by ZCRB1-induced circHEATR5B suppresses aerobic glycolysis of GBM through phosphorylation of JMJD5

Background

RNA-binding proteins (RBPs) and circular RNAs (circRNAs) play important roles in glioblastoma multiforme (GBM). Aerobic glycolysis is a metabolic characteristic of GBM. However, the roles of RBPs and circRNAs in aerobic glycolysis in GBM remain unclear. The aim of this study is to explore the mechanisms by which RBPs and circRNAs regulate aerobic glycolysis in GBM cells.

Methods

RNA sequencing and circRNA microarray analysis were performed to identify RBPs and circRNAs for further study. Mass spectrometry validated the encoded protein and its interacting proteins. Quantitative reverse transcription PCR and western blot assays were used to determine the mRNA and protein expression, respectively. Furthermore, immunofluorescence and fluorescence in situ hybridization assays were used to determine the protein and RNA localization, respectively. Glucose and lactate measurement assays, Seahorse XF glycolysis stress assays and cell viability assays were conducted to investigate the effects on glycolysis and proliferation in GBM cells.

Results

We selected zinc finger CCHC-type and RNA-binding motif 1 (ZCRB1) and circRNA HEAT repeat containing 5B (circHEATR5B) as candidates for this study. These genes were expressed at low levels in GBM tissues and cells. Both ZCRB1 and circHEATR5B overexpression suppressed aerobic glycolysis and proliferation in GBM cells. ZCRB1 overexpression promoted the Alu element-mediated formation of circHEATR5B. In addition, circHEATR5B encoded a novel protein HEATR5B-881aa which interacted directly with Jumonji C-domain-containing 5 (JMJD5) and reduced its stability by phosphorylating S361. JMJD5 knockdown increased pyruvate kinase M2 (PKM2) enzymatic activity and suppressed glycolysis and proliferation in GBM cells. Finally, ZCRB1, circHEATR5B and HEATR5B-881aa overexpression inhibited GBM xenograft growth and prolonged the survival time of nude mice.

Conclusions

This study reveals a novel mechanism of regulating aerobic glycolysis and proliferation in GBM cells through the ZCRB1/circHEATR5B/HEATR5B-881aa/JMJD5/PKM2 pathway, which can provide novel strategies and potential targets for GBM therapy.

A PARP1 PROTAC as a novel strategy against PARP inhibitor resistance via promotion of ferroptosis in p53-positive breast cancer

Therapeutic targeting of the nuclear enzyme poly (ADP-ribose) polymerase 1 (PARP1) with PARP inhibitors (PARPis) in patients with a homologous recombination (HR)- deficient phenotype based on the mechanism of synthetic lethality has been shown tremendous success in cancer therapy. With the clinical use of various PARPis, emerging evidence has shown that some PARPis offer hope for breakthroughs in triple-negative breast cancer (TNBC) therapy, regardless of HR status. However, similar to other conventional cytotoxic drugs, PARPis are also subject to the intractable problem of drug resistance. Notably, acquired resistance to PARPis caused by point mutations in the PARP1 protein is hard to overcome with current strategies. To explore modalities to overcome resistance and identify patients who are most likely to benefit from PARP1-targeted therapy, we developed a proteolysis-targeted chimaera (PROTAC) to degrade mutant PARP1 in TNBC. Here, we investigated a PARP1 PROTAC termed "NN3″, which triggered ubiquitination and proteasome-mediated degradation of PARP1. Moreover, NN3 degraded PARP1 with resistance-related mutations. Interestingly, compared with other reported PARP1 degraders, NN3 exhibited a unique antitumor mechanism in p53-positive breast cancer cells that effectively promoted ferroptosis by downregulating the SLC7A11 pathway. Furthermore, NN3 showed potent activity and low toxicity in vivo. In conclusion, we propose PROTAC-mediated degradation of PARP1 as a novel strategy against mutation-related PARPi resistance and a paradigm for targeting breast cancer with functional p53 via ferroptosis induction.

Targeting the lncMST-EPRS/HSP90AB1 complex as novel therapeutic strategy for T-2 toxin-induced growth retardation

Growth retardation is a global public health problem that is highly prevalent especially in low-and middle-income countries, which is closely related to the consumption of grains contaminated with T-2 toxin, a risk for human and animal health. However, the possible targets that can relieve T-2 toxin-induced growth retardation still need to be explored. In the present study, T-2 toxin was used as an environmental exposure factor to induce growth retardation and further explore the regulatory role of lncRNA in growth retardation. The present study systematically characterised the expression profiles of lncRNAs and identified a lncRNA lncMST that is related to growth retardation in T-2 toxin-administered rats. Functionally, lncMST could alleviate cell cycle arrest and apoptosis in T-2 toxin-treated GH3 cells. Mechanistically, lncMST, serve as an inducible chaperone RNA, involved in the paradigm "Chemical-induced stress related growth retardation", through recruiting the EPRS/HSP90AB1 complex to increase HDAC6 expression, thus further alleviating T-2 toxin-induced growth retardation. These findings for the first time demonstrate that the probable therapeutic relationship between lncMST and growth retardation, providing an explanation and therapeutic targets for the pathogenesis of growth retardation.

Reverse engineering approach: a step towards a new era of vaccinology with special reference to Salmonella

INTRODUCTION

Salmonella is responsible for causing enteric fever, septicemia, and gastroenteritis in humans. Due to high disease burden and emergence of multi- and extensively drug-resistant Salmonella strains, it is becoming difficult to treat the infection with existing battery of antibiotics as we are not able to discover newer antibiotics at the same pace at which the pathogens are acquiring resistance. Though vaccines against Salmonella are available commercially, they have limited efficacy. Advancements in genome sequencing technologies and immunoinformatics approaches have solved the problem significantly by giving rise to a new era of vaccine designing, i.e. 'Reverse engineering.' Reverse engineering/vaccinology has expedited the vaccine identification process. Using this approach, multiple potential proteins/epitopes can be identified and constructed as a single entity to tackle enteric fever.

AREAS COVERED

This review provides details of reverse engineering approach and discusses various protein and epitope-based vaccine candidates identified using this approach against typhoidal Salmonella.

EXPERT OPINION

Reverse engineering approach holds great promise for developing strategies to tackle the pathogen(s) by overcoming the limitations posed by existing vaccines. Progressive advancements in the arena of reverse vaccinology, structural biology, and systems biology combined with an improved understanding of host-pathogen interactions are essential components to design new-generation vaccines.

The endocytic TPLATE complex internalizes ubiquitinated plasma membrane cargo

Endocytosis controls the perception of stimuli by modulating protein abundance at the plasma membrane. In plants, clathrin-mediated endocytosis is the most prominent internalization pathway and relies on two multimeric adaptor complexes, the AP-2 and the TPLATE complex (TPC). Ubiquitination is a well-established modification triggering endocytosis of cargo proteins, but how this modification is recognized to initiate the endocytic event remains elusive. Here we show that TASH3, one of the large subunits of TPC, recognizes ubiquitinated cargo at the plasma membrane via its SH3 domain-containing appendage. TASH3 lacking this evolutionary specific appendage modification allows TPC formation but the plants show severely reduced endocytic densities, which correlates with reduced endocytic flux. Moreover, comparative plasma membrane proteomics identified differential accumulation of multiple ubiquitinated cargo proteins for which we confirm altered trafficking. Our findings position TPC as a key player for ubiquitinated cargo internalization, allowing future identification of target proteins under specific stress conditions.

Could artificial intelligence revolutionize the development of nanovectors for gene therapy and mRNA vaccines?

Gene therapy enables the introduction of nucleic acids like DNA and RNA into host cells, and is expected to revolutionize the treatment of a wide range of diseases. This growth has been further accelerated by the discovery of CRISPR/Cas technology, which allows accurate genomic editing in a broad range of cells and organisms in vitro and in vivo. Despite many advances in gene delivery and the development of various viral and non-viral gene delivery vectors, the lack of highly efficient non-viral systems with low cellular toxicity remains a challenge. The application of cutting-edge technologies such as artificial intelligence (AI) has great potential to find new paradigms to solve this issue. Herein, we review AI and its major subfields including machine learning (ML), neural networks (NNs), expert systems, deep learning (DL), computer vision and robotics. We discuss the potential of AI-based models and algorithms in the design of targeted gene delivery vehicles capable of crossing extracellular and intracellular barriers by viral mimicry strategies. We finally discuss the role of AI in improving the function of CRISPR/Cas systems, developing novel nanobots, and mRNA vaccine carriers.

Successes and Failures of Static Aptamer-Target 3D Docking Models

While Molecular Dynamics simulation programs are probably superior for predicting the binding and affinity of aptamers and their cognate ligands, such molecular dynamics programs require more computing power and analysis time than static docking programs that are more widely accessible to the scientific community on the internet. Static docking programs can be used to investigate the geometric fit of rigid DNA or RNA aptamer 3D structures and their ligands to aid in predicting the relative affinities and cross-reactivity of various potential ligands. Herein, the author describes when such static 3D docking analysis has worked well to produce useful predictions or confirmation of high-affinity aptamer interactions or successful aptamer beacon behavior and when it has not worked well. The analysis of why failures may occur with static 3D computer models is also discussed.

Tryptophan self-assembly yields cytotoxic nanofibers containing amyloid-mimicking and cross-seeding competent conformers

Dietary consumption of Trp via protein-based foods is essential for the maintenance of crucial metabolic processes including the synthesis of proteins and several vital metabolites such as serotonin, melatonin, acetyl CoA, and NADP. However, the abnormal build-up of Trp is known to cause familial hypertryptophanemia and several brain-related medical complications. The molecular mechanism of the onset of such Trp-driven health issues is largely unknown. Here, we show that Trp, under the physiologically mimicked conditions of temperature and buffer, undergoes a concentration driven self-assembly process, yielding amyloid-mimicking nanofibers. Viable H-bonds, π-π interactions and hydrophobic contacts between optimally coordinated Trp molecules become important factors for the formation of a Trp nanoassembly that displays a hydrophobic exterior and a hydrophilic interior. Importantly, Trp nanofibers were found to possess high affinity for native proteins, and they act as cross-seeding competent conformers capable of nucleating amyloid formation in globular proteins including whey protein β-lactoglobulin and type II diabetes linked insulin hormone. Moreover, these amyloid mimicking Trp nanostructures showed toxic effects on neuroblastoma cells. Since the key symptoms in hypertryptophanemia such as behavioural defects and brain-damaging oxidative stress are also observed in amyloid related disorders, our findings on amyloid-like Trp-nanofibers may help in the mechanistic understanding of Trp-related complications and these findings are equally important for innovation in applied nanomaterials design and strategies.

The human RNA polymerase I structure reveals an HMG-like docking domain specific to metazoans

We characterize the human RNA polymerase I by evolutionary biochemistry and cryo-EM revealing a built-in structural domain that apparently serves as transcription factor–binding platform in metazoans.

Transcription of the ribosomal RNA precursor by RNA polymerase (Pol) I is a major determinant of cellular growth, and dysregulation is observed in many cancer types. Here, we present the purification of human Pol I from cells carrying a genomic GFP fusion on the largest subunit allowing the structural and functional analysis of the enzyme across species. In contrast to yeast, human Pol I carries a single-subunit stalk, and in vitro transcription indicates a reduced proofreading activity. Determination of the human Pol I cryo-EM reconstruction in a close-to-native state rationalizes the effects of disease-associated mutations and uncovers an additional domain that is built into the sequence of Pol I subunit RPA1. This “dock II” domain resembles a truncated HMG box incapable of DNA binding which may serve as a downstream transcription factor–binding platform in metazoans. Biochemical analysis, in situ modelling, and ChIP data indicate that Topoisomerase 2a can be recruited to Pol I via the domain and cooperates with the HMG box domain–containing factor UBF. These adaptations of the metazoan Pol I transcription system may allow efficient release of positive DNA supercoils accumulating downstream of the transcription bubble.

Development of CD163 receptor-based enzyme-linked immunosorbent assay for diagnosis of porcine reproductive and respiratory syndrome virus

Porcine reproductive and respiratory syndrome (PRRS) is an important economical disease in the global swine industry. The accurate detection of the PRRS virus (PRRSV) antigen is essential for the disease control and prevention programme. In this study, an indirect enzyme-linked immunosorbent test (PRRSVCD163-iELISA) was developed for the detection of the PRRSV antigen in samples of post-mortem swine tissue using the recombinant pig CD163 receptor protein as the capture ligand. The test was found to be specific for PRRSV, with no cross-reactions with other prevalent pig viral pathogens. The assay was validated by testing 217 post-mortem porcine tissue samples and the results were found to be satisfactory with a relative accuracy of 88.88%. Our assay is also quite precise, with intra- and inter-assay CVs of 6% and 10%, respectively. These findings imply that the PRRSVCD163-iELISA developed is capable of detecting the PRRSV antigen in swine post-mortem tissue samples. This research showed that porcine CD163, the PRRSV cellular receptor, can be exploited to build a diagnostic technique for the detection of PRRSV antigen.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03376-z.

Computationally Designed Anti-LuxP DNA Aptamer Suppressed Flagellar Assembly- and Quorum Sensing-Related Gene Expression in Vibrio parahaemolyticus

(1) Background: Quorum sensing (QS) is the chemical communication between bacteria that sense chemical signals in the bacterial population to control phenotypic changes through the regulation of gene expression. The inhibition of QS has various potential applications, particularly in the prevention of bacterial infection. QS can be inhibited by targeting the LuxP, a periplasmic receptor protein that is involved in the sensing of the QS signaling molecule known as the autoinducer 2 (AI-2). The sensing of AI-2 by LuxP transduces the chemical information through the inner membrane sensor kinase LuxQ protein and activates the QS cascade. (2) Methods: An in silico approach was applied to design DNA aptamers against LuxP in this study. A method combining molecular docking and molecular dynamics simulations was used to select the oligonucleotides that bind to LuxP, which were then further characterized using isothermal titration calorimetry. Subsequently, the bioactivity of the selected aptamer was examined through comparative transcriptome analysis. (3) Results: Two aptamer candidates were identified from the ITC, which have the lowest dissociation constants (Kd) of 0.2 and 0.5 micromolar. The aptamer with the lowest Kd demonstrated QS suppression and down-regulated the flagellar-assembly-related gene expression. (4) Conclusions: This study developed an in silico approach to design an aptamer that possesses anti-QS properties.

HDOCK update for modeling protein-RNA/DNA complex structures

Protein-nucleic acid interactions are involved in various cellular processes. Therefore, determining the structures of protein-nucleic acid complexes can provide insights into the mechanisms of the interactions and thus guide the rational drug design to modulate these interactions. Due to the high cost and technical difficulties of solving complex structures experimentally, computational modeling such as molecular docking has been playing an important role in the study of molecular interactions. In order to make it easier for researchers to obtain biomolecular complex structures through molecular docking, we developed the HDOCK server for protein-protein and protein-RNA/DNA docking (accessed at http://hdock.phys.hust.edu.cn/). Since its first release in 2017, HDOCK has been widely used in the scientific community. As nucleic acids may include single-stranded (ss) RNA/DNA and double-stranded (ds) RNA/DNA, we now present an updated version of HDOCK, which offers new options for structural modeling of ssRNA, ssDNA, dsRNA, and dsDNA. We hope this update will better help the scientific community solve important biological problems, thereby advancing the field. In this article, we describe the general protocol of HDOCK with emphasis on the new functions on RNA/DNA modeling. Several application examples are also given to illustrate the usage of the new functions.

TssI2-TsiI2 of Vibrio fluvialis VflT6SS2 delivers pesticin domain-containing periplasmic toxin and cognate immunity that modulates bacterial competitiveness

Vibrio fluvialis is a halophilic Gram-negative bacterium regarded as an emerging unusual enteric pathogen of increasing public health concern. Our previous work has identified two type VI secretion systems (T6SSs) in V. fluvialis, VflT6SS1, and VflT6SS2, and the latter is functional in mediating interbacterial competitiveness. However, its antibacterial effectors remain to be clarified. In this work, we focused on a new potential effector/immunity pair TssI2/TsiI2. Bioinformatics analysis revealed that the C-terminal domain of TssI2 belongs to a widespread family of pesticin, and its antibacterial toxicity and corresponding protection by TsiI2 were proved via bacterial killing assays, and their action sites were localized to the periplasm of bacterial cells. The interaction of TssI2 and TsiI2 was demonstrated by the bacterial adenylate cyclase two-hybrid, protein pull-down and isothermal titration calorimetry assays. Site-directed mutagenesis demonstrated that, in addition to Glu-844, Thr-863, and Asp-869, which correspond to three reported residues in pesticin of Yersinia pestis, additional residues including Phe-837, Gly-845, Tyr-851, Gly-867, Gln-963, Trp-975, and Arg-1000 were also proved to be crucial to the bactericidal activity of TssI2. Muramidase/lysozyme-related peptidoglycan (PG) hydrolase activities of TssI2 and its variants were validated with permeabilized Escherichia coli cells and purified PG substrate. Based on sequence homologies at C-terminals in various V. fluvialis isolates, TssI2 was subdivided into five clusters (12-22% identity among them), and the antibacterial activities of representative effectors from other four Clusters were also confirmed through periplasmic over-expression in E. coli host. Two selected cognate immunities were proved to confer protection against the toxicities of their effectors. Additionally, TsiI2, which belongs to Cluster I, exhibited cross-protection to effector from Cluster V. Together, current findings expand our knowledge of the diversity and consistency of evolved VgrG effectors in V. fluvialis and on how VflT6SS2 mediates a competitive advantage to gain a better survival.

CD82 attenuates TGF-β1-mediated epithelial-mesenchymal transition by blocking smad-dependent signaling in ARPE-19 cells

In retinal pigment epithelial (RPE) cells, transforming growth factor-beta (TGF-β) plays a critical role in epithelial-mesenchymal transition (EMT), which contributes to various fibrotic retinal disorders. In the present study, we investigated the effect of recombinant human cluster of differentiation 82 (rhCD82), a tumor metastasis suppressor, on TGF-β-induced EMT in the human RPE cell line APRE-19. The results show that TGF-β1 significantly enhanced cell migration, invasion and the expression of EMT-mediate factors in ARPE-19 cells. However, rhCD82 markedly inhibited cell mobility and the expression of epithelial marker, zonula occludens-1, as well as increased the expression of mesenchymal markers, such as vimentin and α-smooth muscle actin in TGF-β1-treated APRE-19 cells. In addition, TGF-β1 upregulated the phosphorylation of Smad, extracellular signal regulated kinase (ERK) and glycogen synthase kinase-3β (GSK-3β), but only phosphorylation of Smad was suppressed by rhCD82. Noteworthy, rhCD82 greatly suppressed the expression of TGF-β receptor I (TGFRI), TGFRII and integrins in TGF-β1-treated APRE-19 cells. In particular, the result of molecular docking analysis and structural modeling show that rhCD82 partially interacts with the TGF-β1 binding sites of TGFRI, TGFRII, integrin β1 and integrin αv. Taken together, this finding suggested that rhCD82 suppressed TGF-β1-induced EMT of RPE by blocking of Smad-dependent pathway, which is caused by rhCD82 interaction with TGFRs and integrins, suggesting new insight into CD82 as a potential therapeutic strategy in fibrotic retinal disorders.

The production of the first functional antibody mimetic in higher plants: the chloroplast makes the DARPin G3 for HER2 imaging in oncology

Background

Designed mimetic molecules are attractive tools in biopharmaceuticals and synthetic biology. They require mass and functional production for the assessment of upcoming challenges in the near future. The DARPin family is considered a mimetic pharmaceutical peptide group with high affinity binding to specific targets. DARPin G3 is designed to bind to the HER2 (human epidermal growth factor receptor 2) tyrosine kinase receptor. Overexpression of HER2 is common in some cancers, including breast cancer, and can be used as a prognostic and predictive tool for cancer. The chloroplasts are cost-effective alternatives, equal to, and sometimes better than, bacterial, yeast, or mammalian expression systems. This research examined the possibility of the production of the first antibody mimetic, DARPin G3, in tobacco chloroplasts for HER2 imaging in oncology.

Results

The chloroplast specific DARPin G3 expression cassette was constructed and transformed into N. tabacum chloroplasts. PCR and Southern blot analysis confirmed integration of transgenes as well as chloroplastic and cellular homoplasmy. The Western blot analysis and ELISA confirmed the production of DARPin G3 at the commercial scale and high dose with the rate of 20.2% in leaf TSP and 33.7% in chloroplast TSP. The functional analysis by ELISA confirmed the binding of IMAC purified chloroplast-made DARPin G3 to the extracellular domain of the HER2 receptor with highly effective picomolar affinities. The carcinoma cellular studies by flow cytometry and immunofluorescence microscopy confirmed the correct functioning by the specific binding of the chloroplast-made DARPin G3 to the HER2 receptor on the surface of HER2-positive cancer cell lines.

Conclusion

The efficient functional bioactive production of DARPin G3 in chloroplasts led us to introduce plant chloroplasts as the site of efficient production of the first antibody mimetic molecules. This report, as the first case of the cost-effective production of mimetic molecules, enables researchers in pharmaceuticals, synthetic biology, and bio-molecular engineering to develop tool boxes by producing new molecular substitutes for diverse purposes.

A Novel ssDNA Aptamer Targeting Carcinoembryonic Antigen: Selection and Characterization

One of the major causes of a drastically shorter life expectancy and one of the most prevalent diseases in the world today is cancer. Given the data on the rise in cancer cases throughout the world, it is obvious that, despite the diagnostic techniques currently being used, there is a pressing need to develop precise and sensitive techniques for early diagnosis of the disease. A high degree of affinity and specificity towards particular targets is maintained by the short nucleic acid molecules known as aptamers. Aptamers outperform antibodies due to their unique benefits, such as their simplicity in synthesis and modification, lack of toxicity, and long-term stability. Utilizing an accurate recognition element and a robust signal transduction mechanism, molecular diagnostics can be extremely sensitive and specific. In this study, development of new single-stranded DNA aptamers against CEA for use in cancer diagnostics was accomplished using SELEX and NGS methods. As a result of 12 iterative SELEX rounds, nine aptamer candidates against CEA were developed. NGS comparative analysis revealed that round twelve had an enriched number of aptamers that were specifically bound, as opposed to round eight. Among the selected nine sequences characterized by bioinformatics analysis and ELONA, an aptamer sequence with the highest specificity and affinity for the target protein was identified and further examined. Aptamer sequence (6) was screened in a concentration-dependent assay, specificity analysis was performed, and its potential secondary and tertiary structures were predicted, which enabled us to test one of the possible putative interactions with CEA. Finally, aptamer sequence (6) labelled with a Cy5 fluorescent tag was used in confocal microscopy to observe its binding towards the CEA expressed in HT-29 human colon adenocarcinoma cell line.