Exploring the mechanism of Corbrin capsules in the intervention of AKI-COVID-19 based on network pharmacology combined with GEO dataset.

AIM: of the study: This study used network pharmacology and the Gene Expression Omnibus (GEO) database to investigate the therapeutic effects of Corbrin capsules on acute kidney injury (AKI)-COVID-19 (coronavirus disease 2019). MATERIALS AND METHODS: The active constituents and specific molecular targets of Corbrin capsules were obtained from the Traditional Chinese Medicine Systems Pharmacology (TCMSP) and Swiss Target Prediction databases. The targets related to AKI and COVID-19 disease were obtained from the Online Mendelian Inheritance in Man (OMIM), GeneCards, and GEO databases. A protein-protein interaction (PPI) network was constructed by utilizing Cytoscape. To enhance the analysis of pathways associated with the pathogenesis of AKI-COVID-19, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed. Furthermore, immune infiltration analysis was performed by using single-sample gene set enrichment analysis (ssGSEA) and CIBERSORT. Molecular docking was used to assess interactions between differentially expressed genes and active ingredients. Verification was performed by utilizing GEO databases and in vivo assays. RESULTS: This study revealed an overlap of 18 significantly differentially expressed genes between the Corbrin capsules group and the AKI-COVID-19 target group. Analysis of the PPI network identified TP53, JAK2, PIK3CA, PTGS2, KEAP1, and MCL1 as the top six core protein targets with the highest degrees. The results obtained from GO and KEGG analyses demonstrated that the target genes were primarily enriched in the apoptosis and JAK-STAT signaling pathways. Moreover, the analysis of immune infiltration revealed a notable disparity in the percentage of quiescent memory CD4 + T cells. Western blot analyses provided compelling evidence suggesting that the dysregulation of 6 core protein targets could be effectively reversed by Corbrin capsules. CONCLUSION: This study revealed the key components, targets, and pathways involved in treating AKI-related COVID-19 using Corbrin capsules. This study also provided a new understanding of the molecular mechanisms underlying this treatment.

LncRNAs expression in adriamycin-induced rats reveals the potential role of LncRNAs contributing to chronic glomerulonephritis pathogenesis.

BACKGROUND: Chronic glomerulonephritis (CGN) is the most common form of primary glomerular disease with unclear molecular mechanisms. Currently, limited study on long non-coding RNAs (lncRNAs) in CGN is available. Our study aimed to identify potential lncRNAs and genes in the normal and adriamycin-induced CGN rats, which to explore the potential molecular mechanisms of CGN pathogenesis. METHODS: To identify LncRNAs specifically expressed in CGN, the expression of LncRNAs in glomerular tissues of rats from the adriamycin-induced group (n = 3) was compared with that in the control group (n = 3) using RNA-sequencing and real-time polymerase chain reaction (RT-PCR). Identification of differentially expressed lncRNAs and mRNAs were performed between the 2 groups. Gene ontology (GO) and pathway enrichment analyses were performed to analyze the biological functions and pathways for the differentially expressed mRNAs. LncRNA-mRNA co-expression network was constructed to analyses for the genes. The protein-protein interaction (PPI) network was visualized. RESULTS: A total of 114 significantly up-regulated and 86 down-regulated lncRNAs, 1038 up-regulated and 88 down-regulated mRNAs were identified. Additionally, Il6, with the highest connectivity degree in PPI network, was noteworthy enriched in various kinds pathways. Coding-non-coding gene co-expression networks (CNC network) were drawn based on the correlation analysis between the differentially expressed LncRNAs and mRNAs. Ten LncRNAs, NONRATT000964.2, NONRATT018086.2, NONRATT023684.2, NONRATT009530.2, NONRATT006315.2, NONRATT026805.2, MSTRG.9260.1, NONRATT009155.2, MSTRG.7681.1, NONRATT009275.2, were selected to analyze the relationship between LncRNAs and CGN via the CNC network and GO analysis. Real-time PCR result confirmed that the six LncRNAs were specifically expressed in the CGN rats. CONCLUSIONS: The ten LncRNAs were differentially expressed and might play important roles in the development of CGN. Key genes, such as Il6, Ptprc, TOP2a, Fos, Myc, etc., may be crucial biomarkers for CGN.

Characterization of ribonuclease III from Brucella.

Bacterial ribonuclease III (RNase III) is a highly conserved endonuclease, which plays pivotal roles in RNA maturation and decay pathways by cleaving double-stranded structure of RNAs. Here we cloned rncS gene from the genomic DNA of Brucella melitensis, and analyzed the cleavage properties of RNase III from Brucella. We identified Brucella-encoding small RNA (sRNA) by high-throughput sequencing and northern blot, and found that sRNA of Brucella and Homo miRNA precursor (pre-miRNA) can be bound and cleaved by B.melitensis ribonuclease III (Bm-RNase III). Cleavage activity of Bm-RNase III is bivalent metal cations- and alkaline buffer-dependent. We constructed several point mutations in Bm-RNase III, whose cleavage activity indicated that the 133th Glutamic acid residue was required for catalytic activity. Western blot revealed that Bm-RNase III was differently expressed in Brucella virulence strain 027 and vaccine strain M5-90. Collectively, our data suggest that Brucella RNase III can efficiently bind and cleave stem-loop structure of small RNA, and might participate in regulation of virulence in Brucella.

Immunogenic response to a recombinant pseudorabies virus carrying bp26 gene of Brucella melitensis in mice.

Brucellae are facultative intracellular bacterial pathogens of a zoonotic disease called brucellosis. Live attenuated vaccines are utilized for prophylaxis of brucellosis; however, they retain residual virulence to human and/or animals, as well as interfere with diagnosis. In this study, recombinant virus PRV ΔTK/ΔgE/bp26 was screened and purified. One-step growth curve assay showed that the titer of recombinant virus was comparable to the parent strain. Mice experiments showed the recombinant virus elicited high titer of humoral antibodies against Brucella detected by enzyme-linked immunosorbent assay and against PRV by serum neutralization test. The recombinant virus induced high level of Brucella-specific lymphocyte proliferation response and production of interferon gamma. Collectively, these data suggest that the bivalent virus was capable of inducing both humoral and cellular immunity, and had the potential to be a vaccine candidate to prevent Brucella and/or pseudorabies virus infections.

Detection of persistent organic pollutants binding modes with androgen receptor ligand binding domain by docking and molecular dynamics.

BACKGROUND: Persistent organic pollutants (POPs) are persistent in the environment after release from industrial compounds, combustion productions or pesticides. The exposure of POPs has been related to various reproductive disturbances, such as reduced semen quality, testicular cancer, and imbalanced sex ratio. Among POPs, dichlorodiphenyldichloroethylene (4,4'-DDE) and polychlorinated biphenyls (PCBs) are the most widespread and well-studied compounds. Recent studies have revealed that 4,4'-DDE is an antagonist of androgen receptor (AR). However, the mechanism of the inhibition remains elusive. CB-153 is the most common congener of PCBs, while the action of CB-153 on AR is still under debate. RESULTS: Molecular docking and molecular dynamics (MD) approaches have been employed to study binding modes and inhibition mechanism of 4,4'-DDE and CB-153 against AR ligand binding domain (LBD). Several potential binding sites have been detected and analyzed. One possible binding site is the same binding site of AR natural ligand androgen 5α-dihydrotestosterone (DHT). Another one is on the ligand-dependent transcriptional activation function (AF2) region, which is crucial for the co-activators recruitment. Besides, a novel possible binding site was observed for POPs with low binding free energy with the receptor. Detailed interactions between ligands and the receptor have been represented. The disrupting mechanism of POPs against AR has also been discussed. CONCLUSIONS: POPs disrupt the function of AR through binding to three possible biding sites on AR/LBD. One of them shares the same binding site of natural ligand of AR. Another one is on AF2 region. The third one is in a cleft near N-terminal of the receptor. Significantly, values of binding free energy of POPs with AR/LBD are comparable to that of natural ligand androgen DHT.

Identification of functionally key residues in AMPA receptor with a thermodynamic method.

AMPA receptor mediates the fast excitatory synaptic transmission in the central nervous system, and it is activated by the binding of glutamate that results in the opening of the transmembrane ion channel. In the present work, the thermodynamic method developed by our group was improved and then applied to identify the functionally key residues that regulate the glutamate-binding affinity of AMPA receptor. In our method, the key residues are identified as those whose perturbation largely changes the ligand binding free energy of the protein. It is found that besides the ligand binding sites, other residues distant from the binding cleft can also influence the glutamate binding affinity through a long-range allosteric regulation. These allosteric sites include the hinge region of the ligand binding cleft, the dimer interface of the ligand binding domain, the linkers between the ligand binding domain and the transmembrane domain, and the interface between the N-terminal domain and the ligand binding domain. Our calculation results are consistent with the available experimental data. The results are helpful for our understanding of the mechanism of long-range allosteric communication in the AMPA receptor and the mechanism of channel opening triggered by glutamate binding.

The interactions and recognition of cyclic peptide mimetics of Tat with HIV-1 TAR RNA: a molecular dynamics simulation study.

The interaction of HIV-1 trans-activator protein Tat with its cognate trans-activation response element (TAR) RNA is critical for viral transcription and replication. Therefore, it has long been considered as an attractive target for the development of antiviral compounds. Recently, the conformationally constrained cyclic peptide mimetics of Tat have been tested to be a promising family of lead peptides. Here, we focused on two representative cyclic peptides termed as L-22 and KP-Z-41, both of which exhibit excellent inhibitory potency against Tat and TAR interaction. By means of molecular dynamics simulations, we obtained a detailed picture of the interactions between them and HIV-1 TAR RNA. In results, it is found that the binding modes of the two cyclic peptides to TAR RNA are almost identical at or near the bulge regions, whereas the binding interfaces at the apical loop exhibit large conformational heterogeneity. In addition, it is revealed that electrostatic interaction energy contributes much more to KP-Z-41 complex formation than to L-22 complex, which is the main source of energy that results in a higher binding affinity of KP-Z-41 over-22 for TAR RNA. Furthermore, we identified a conserved motif RRK (Arg-Arg-Lys) that is shown to be essential for specific binding of this class of cyclic peptides to TAR RNA. This work can provide a useful insight into the design and modification of cyclic peptide inhibitors targeting the association of HIV-1 Tat and TAR RNA.

MicroRNA expression profile in RAW264.7 cells in response to Brucella melitensis infection.

MicroRNA (miRNA) is small non-coding RNA with approximate 22 nt in length. Recent studies indicate that miRNAs play significant roles in pathogen-host interactions. Brucella organisms are Gram-negative facultative intracellular bacteria that cause Brucellosis. Brucella strains infect macrophages and establish chronic infection by altering host life activities including apoptosis and autophagy. Here, we report a comprehensive analysis of miRNA expression profiles in mock- and Brucella-infected RAW264.7 cells using high-throughput sequencing approach. In total, 344 unique miRNAs were co-expressed in the two libraries, in which 57 miRNAs were differentially expressed. Eight differentially expressed miRNAs with high abundance were subjected to further analysis. The GO enrichment analysis suggests that the putative target genes of these differentially expressed miRNAs are involved in apoptosis, autophagy and immune response. In particular, a total of 25 target genes are involved in regulating apoptosis and autophagy, indicating that these miRNAs may play important regulatory roles in the Brucella-host interactions. Furthermore, the interactions of miR-1981 and its target genes, Bcl-2 and Bid, were validated by luciferase assay. The results show that miR-1981 mimic up-regulated the luciferase activity of psiCHECK-2 Bcl-2 3' UTR, but the luciferase activity of psiCHECK-2 Bid 3' UTR was not changed significantly. Taken together, these data provide valuable framework on Brucella induced miRNA expression in RAW264.7 cells, and suggest that Brucella may establish chronic infection by regulating miRNA expression profile.

Identification of key residues for protein conformational transition using elastic network model.

Proteins usually undergo conformational transitions between structurally disparate states to fulfill their functions. The large-scale allosteric conformational transitions are believed to involve some key residues that mediate the conformational movements between different regions of the protein. In the present work, a thermodynamic method based on the elastic network model is proposed to predict the key residues involved in protein conformational transitions. In our method, the key functional sites are identified as the residues whose perturbations largely influence the free energy difference between the protein states before and after transition. Two proteins, nucleotide binding domain of the heat shock protein 70 and human/rat DNA polymerase ß, are used as case studies to identify the critical residues responsible for their open-closed conformational transitions. The results show that the functionally important residues mainly locate at the following regions for these two proteins: (1) the bridging point at the interface between the subdomains that control the opening and closure of the binding cleft; (2) the hinge region between different subdomains, which mediates the cooperative motions between the corresponding subdomains; and (3) the substrate binding sites. The similarity in the positions of the key residues for these two proteins may indicate a common mechanism in their conformational transitions.

An analysis of the influence of protein intrinsic dynamical properties on its thermal unfolding behavior.

The influence of the protein topology-encoded dynamical properties on its thermal unfolding motions was studied in the present work. The intrinsic dynamics of protein topology was obtained by the anisotropic network model (ANM). The ANM has been largely used to investigate protein collective functional motions, but it is not well elucidated if this model can also reveal the preferred large-scale motions during protein unfolding. A small protein barnase is used as a typical case study to explore the relationship between protein topology-encoded dynamics and its unfolding motions. Three thermal unfolding simulations at 500 K were performed for barnase and the entire unfolding trajectories were sampled and partitioned into several windows. For each window, the preferred unfolding motions were investigated by essential dynamics analysis, and then associated with the intrinsic dynamical properties of the starting conformation in this window, which is detected by ANM. The results show that only a few slow normal modes imposed by protein structure are sufficient to give a significant overlap with the preferred unfolding motions. Especially, the large amplitude unfolding movements, which imply that the protein jumps out of a local energy basin, can be well described by a single or several ANM slow modes. Besides the global motions, it is also found that the local residual fluctuations encoded in protein structure are highly correlated with those in the protein unfolding process. Furthermore, we also investigated the relationship between protein intrinsic flexibility and its unfolding events. The results show that the intrinsic flexible regions tend to unfold early. Several early unfolding events can be predicted by analysis of protein structural flexibility. These results imply that protein structure-encoded dynamical properties have significant influences on protein unfolding motions.