MCCS: a novel recognition pattern-based method for fast track discovery of anti-SARS-CoV-2 drugs

Acarbose reduces Pseudomonas aeruginosa respiratory tract infection in type 2 diabetic mice

BACKGROUND

Type 2 diabetes mellitus (T2DM) is widely prevalent worldwide, and respiratory tract infections (RTIs) have become the primary cause of death for T2DM patients who develop concurrent infections. Among these, Pseudomonas aeruginosa infection has been found to exhibit a high mortality rate and poor prognosis and is frequently observed in bacterial infections that are concurrent with COVID-19. Studies have suggested that acarbose can be used to treat T2DM and reduce inflammation. Our objective was to explore the effect of acarbose on P. aeruginosa RTI in T2DM individuals and elucidate its underlying mechanism.

METHODS

High-fat diet (HFD) induction and P. aeruginosa inhalation were used to establish a RTI model in T2DM mice. The effect and mechanism of acarbose administered by gavage on P. aeruginosa RTI were investigated in T2DM and nondiabetic mice using survival curves, pathological examination, and transcriptomics.

RESULTS

We found that P. aeruginosa RTI was more severe in T2DM mice than in nondiabetic individuals, which could be attributed to the activation of the NF-κB and TREM-1 signaling pathways. When acarbose alleviated P. aeruginosa RTI in T2DM mice, both HIF-1α and NF-κB signaling pathways were inhibited. Furthermore, inhibition of the calcium ion signaling pathway and NF-κB signaling pathway contributed to the attenuation of P. aeruginosa RTI by acarbose in nondiabetic mice.

CONCLUSIONS

This study confirmed the attenuating effect of acarbose on P. aeruginosa RTIs in T2DM and nondiabetic mice and investigated its mechanism, providing novel support for its clinical application in related diseases.

In Silico Models for Anti-COVID-19 Drug Discovery: A Systematic Review

The coronavirus disease 2019 (COVID-19) is a severe worldwide pandemic. Due to the emergence of various SARS-CoV-2 variants and the presence of only one Food and Drug Administration (FDA) approved anti-COVID-19 drug (remdesivir), the disease remains a mindboggling global public health problem. Developing anti-COVID-19 drug candidates that are effective against SARS-CoV-2 and its various variants is a pressing need that should be satisfied. This systematic review assesses the existing literature that used in silico models during the discovery procedure of anti-COVID-19 drugs. Cochrane Library, Science Direct, Google Scholar, and PubMed were used to conduct a literature search to find the relevant articles utilizing the search terms "In silico model," "COVID-19," "Anti-COVID-19 drug," "Drug discovery," "Computational drug designing," and "Computer-aided drug design." Studies published in English between 2019 and December 2022 were included in the systematic review. From the 1120 articles retrieved from the databases and reference lists, only 33 were included in the review after the removal of duplicates, screening, and eligibility assessment. Most of the articles are studies that use SARS-CoV-2 proteins as drug targets. Both ligand-based and structure-based methods were utilized to obtain lead anti-COVID-19 drug candidates. Sixteen articles also assessed absorption, distribution, metabolism, excretion, toxicity (ADMET), and drug-likeness properties. Confirmation of the inhibitory ability of the candidate leads by in vivo or in vitro assays was reported in only five articles. Virtual screening, molecular docking (MD), and molecular dynamics simulation (MDS) emerged as the most commonly utilized in silico models for anti-COVID-19 drug discovery.

Interactome-Based Machine Learning Predicts Potential Therapeutics for COVID-19

COVID-19, the disease caused by SARS-CoV-2, has been disrupting our lives for more than two years now. SARS-CoV-2 interacts with human proteins to pave its way into the human body, thereby wreaking havoc. Moreover, the mutating variants of the virus that take place in the SARS-CoV-2 genome are also a cause of concern among the masses. Thus, it is very important to understand human-spike protein-protein interactions (PPIs) in order to predict new PPIs and consequently propose drugs for the human proteins in order to fight the virus and its different mutated variants, with the mutations occurring in the spike protein. This fact motivated us to develop a complete pipeline where PPIs and drug-protein interactions can be predicted for human-SARS-CoV-2 interactions. In this regard, initially interacting data sets are collected from the literature, and noninteracting data sets are subsequently created for human-SARS-CoV-2 by considering only spike glycoprotein. On the other hand, for drug-protein interactions both interacting and noninteracting data sets are considered from DrugBank and ChEMBL databases. Thereafter, a model based on a sequence-based feature is used to code the protein sequences of human and spike proteins using the well-known Moran autocorrelation technique, while the drugs are coded using another well-known technique, viz., PaDEL descriptors, to predict new human-spike PPIs and eventually new drug-protein interactions for the top 20 predicted human proteins interacting with the original spike protein and its different mutated variants like Alpha, Beta, Delta, Gamma, and Omicron. Such predictions are carried out by random forest as it is found to perform better than other predictors, providing an accuracy of 90.53% for human-spike PPI and 96.15% for drug-protein interactions. Finally, 40 unique drugs like eicosapentaenoic acid, doxercalciferol, ciclesonide, dexamethasone, methylprednisolone, etc. are identified that target 32 human proteins like ACACA, DST, DYNC1H1, etc.

Bioinformatics-based investigation on the genetic influence between SARS-CoV-2 infections and idiopathic pulmonary fibrosis (IPF) diseases, and drug repurposing

Some recent studies showed that severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections and idiopathic pulmonary fibrosis (IPF) disease might stimulate each other through the shared genes. Therefore, in this study, an attempt was made to explore common genomic biomarkers for SARS-CoV-2 infections and IPF disease highlighting their functions, pathways, regulators and associated drug molecules. At first, we identified 32 statistically significant common differentially expressed genes (cDEGs) between disease (SARS-CoV-2 and IPF) and control samples of RNA-Seq profiles by using a statistical r-package (edgeR). Then we detected 10 cDEGs (CXCR4, TNFAIP3, VCAM1, NLRP3, TNFAIP6, SELE, MX2, IRF4, UBD and CH25H) out of 32 as the common hub genes (cHubGs) by the protein-protein interaction (PPI) network analysis. The cHubGs regulatory network analysis detected few key TFs-proteins and miRNAs as the transcriptional and post-transcriptional regulators of cHubGs. The cDEGs-set enrichment analysis identified some crucial SARS-CoV-2 and IPF causing common molecular mechanisms including biological processes, molecular functions, cellular components and signaling pathways. Then, we suggested the cHubGs-guided top-ranked 10 candidate drug molecules (Tegobuvir, Nilotinib, Digoxin, Proscillaridin, Simeprevir, Sorafenib, Torin 2, Rapamycin, Vancomycin and Hesperidin) for the treatment against SARS-CoV-2 infections with IFP diseases as comorbidity. Finally, we investigated the resistance performance of our proposed drug molecules compare to the already published molecules, against the state-of-the-art alternatives publicly available top-ranked independent receptors by molecular docking analysis. Molecular docking results suggested that our proposed drug molecules would be more effective compare to the already published drug molecules. Thus, the findings of this study might be played a vital role for diagnosis and therapies of SARS-CoV-2 infections with IPF disease as comorbidity risk.

Identification of host genomic biomarkers from multiple transcriptomics datasets for diagnosis and therapies of SARS-CoV-2 infections

The pandemic of COVID-19 is a severe threat to human life and the global economy. Despite the success of vaccination efforts in reducing the spread of the virus, the situation remains largely uncontrolled due to the random mutation in the RNA sequence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which demands different variants of effective drugs. Disease-causing gene-mediated proteins are usually used as receptors to explore effective drug molecules. In this study, we analyzed two different RNA-Seq and one microarray gene expression profile datasets by integrating EdgeR, LIMMA, weighted gene co-expression network and robust rank aggregation approaches, which revealed SARS-CoV-2 infection causing eight hub-genes (HubGs) including HubGs; REL, AURKA, AURKB, FBXL3, OAS1, STAT4, MMP2 and IL6 as the host genomic biomarkers. Gene Ontology and pathway enrichment analyses of HubGs significantly enriched some crucial biological processes, molecular functions, cellular components and signaling pathways that are associated with the mechanisms of SARS-CoV-2 infections. Regulatory network analysis identified top-ranked 5 TFs (SRF, PBX1, MEIS1, ESR1 and MYC) and 5 miRNAs (hsa-miR-106b-5p, hsa-miR-20b-5p, hsa-miR-93-5p, hsa-miR-106a-5p and hsa-miR-20a-5p) as the key transcriptional and post-transcriptional regulators of HubGs. Then, we conducted a molecular docking analysis to determine potential drug candidates that could interact with HubGs-mediated receptors. This analysis resulted in the identification of top-ranked ten drug agents, including Nilotinib, Tegobuvir, Digoxin, Proscillaridin, Olysio, Simeprevir, Hesperidin, Oleanolic Acid, Naltrindole and Danoprevir. Finally, we investigated the binding stability of the top-ranked three drug molecules Nilotinib, Tegobuvir and Proscillaridin with the three top-ranked proposed receptors (AURKA, AURKB, OAS1) by using 100 ns MD-based MM-PBSA simulations and observed their stable performance. Therefore, the findings of this study might be useful resources for diagnosis and therapies of SARS-CoV-2 infections.

Molecular Modeling Study of a Receptor-Orthosteric Ligand-Allosteric Modulator Signaling Complex

Allosteric modulators (AMs) are considered as a perpetual hotspot in research for their higher selectivity and various effects on orthosteric ligands (OL). They are classified in terms of their functionalities as positive, negative, or silent allosteric modulators (PAM, NAM, or SAM, respectively). In the present work, 11 pairs of three-dimensional (3D) structures of receptor-orthosteric ligand and receptor-orthosteric ligand-allosteric modulator complexes have been collected for the studies, including three different systems: GPCR, enzyme, and ion channel. Molecular dynamics (MD) simulations are applied to quantify the dynamic interactions in both the orthosteric and allosteric binding pockets and the structural fluctuation of the involved proteins. Our results showed that MD simulations of moderately large molecules or peptides undergo insignificant changes compared to crystal structure results. Furthermore, we also studied the conformational changes of receptors that bound with PAM and NAM, as well as the different allosteric binding sites in a receptor. There should be no preference for the position of the allosteric binding pocket after comparing the allosteric binding pockets of these three systems. Finally, we aligned four distinct β2 adrenoceptor structures and three N-methyl-d-aspartate receptor (NMDAR) structures to investigate conformational changes. In the β2 adrenoceptor systems, the aligned results revealed that transmembrane (TM) helices 1, 5, and 6 gradually increased outward movement from an enhanced inactive state to an improved active state. TM6 endured the most significant conformational changes (around 11 Å). For NMDAR, the bottom section of NMDAR's ligand-binding domain (LBD) experienced an upward and outward shift during the gradually activating process. In conclusion, our research provides insight into receptor-orthosteric ligand-allosteric modulator studies and the design and development of allosteric modulator drugs using MD simulation.

Discovery of Novel Chinese Medicine Compounds Targeting 3CL Protease by Virtual Screening and Molecular Dynamics Simulation

The transmission and infectivity of COVID-19 have caused a pandemic that has lasted for several years. This is due to the constantly changing variants and subvariants that have evolved rapidly from SARS-CoV-2. To discover drugs with therapeutic potential for COVID-19, we focused on the 3CL protease (3CL^pro) of SARS-CoV-2, which has been proven to be an important target for COVID-19 infection. Computational prediction techniques are quick and accurate enough to facilitate the discovery of drugs against the 3CL^pro of SARS-CoV-2. In this paper, we used both ligand-based virtual screening and structure-based virtual screening to screen the traditional Chinese medicine small molecules that have the potential to target the 3CL^pro of SARS-CoV-2. MD simulations were used to confirm these results for future in vitro testing. MCCS was then used to calculate the normalized free energy of each ligand and the residue energy contribution. As a result, we found ZINC15676170, ZINC09033700, and ZINC12530139 to be the most promising antiviral therapies against the 3CL^pro of SARS-CoV-2.

Elucidation of Binding Features and Dissociation Pathways of Inhibitors and Modulators in SARS-CoV-2 Main Protease by Multiple Molecular Dynamics Simulations

COVID-19 can cause different neurological symptoms in some people, including smell, inability to taste, dizziness, confusion, delirium, seizures, stroke, etc. Owing to the issue of vaccine effectiveness, update and coverage, we still need one or more diversified strategies as the backstop to manage illness. Characterizing the structural basis of ligand recognition in the main protease (Mpro) of SARS-CoV-2 will facilitate its rational design and development of potential drug candidates with high affinity and selectivity against COVID-19. Up to date, covalent-, non-covalent inhibitors and allosteric modulators have been reported to bind to different active sites of Mpro. In the present work, we applied the molecular dynamics (MD) simulations to systematically characterize the potential binding features of catalytic active site and allosteric binding sites in Mpro using a dataset of 163 3D structures of Mpro-inhibitor complexes, in which our results are consistent with the current studies. In addition, umbrella sampling (US) simulations were used to explore the dissociation processes of substrate pathway and allosteric pathway. All the information provided new insights into the protein features of Mpro and will facilitate its rational drug design for COVID-19.

How Do Modulators Affect the Orthosteric and Allosteric Binding Pockets?

Allosteric modulators (AMs) that bind allosteric sites can exhibit greater selectivity than the orthosteric ligands and can either enhance agonist-induced receptor activity (termed positive allosteric modulator or PAM), inhibit agonist-induced activity (negative AM or NAM), or have no effect on activity (silent AM or SAM). Until now, it is not clear what the exact effects of AMs are on the orthosteric active site or the allosteric binding pocket(s). In the present work, we collected both the three-dimensional (3D) structures of receptor-orthosteric ligand and receptor-orthosteric ligand-AM complexes of a specific target protein. Using our novel algorithm toolset, molecular complex characterizing system (MCCS), we were able to quantify the key residues in both the orthosteric and allosteric binding sites along with potential changes of the binding pockets. After analyzing 21 pairs of 3D crystal or cryo-electron microscopy (cryo-EM) complexes, including 4 pairs of GPCRs, 5 pairs of ion channels, 11 pairs of enzymes, and 1 pair of transcription factors, we found that the binding of AMs had little impact on both the orthosteric and allosteric binding pockets. In return, given the accurately predicted allosteric binding pocket(s) of a drug target of medicinal interest, we can confidently conduct the virtual screening or lead optimization without concern that the huge conformational change of the pocket could lead to the low accuracy of virtual screening.

A new glimpse on the active site of SARS-CoV-2 3CLpro, coupled with drug repurposing study

Coronavirus disease 2019 (COVID-19) is caused by novel severe acute respiratory syndrome coronavirus (SARS-CoV-2). Its main protease, 3C-like protease (3CLpro), is an attractive target for drug design, due to its importance in virus replication. The analysis of the radial distribution function of 159 3CLpro structures reveals a high similarity index. A study of the catalytic pocket of 3CLpro with bound inhibitors reveals that the influence of the inhibitors is local, perturbing dominantly only residues in the active pocket. A machine learning based model with high predictive ability against SARS-CoV-2 3CLpro is designed and validated. The model is used to perform a drug-repurposing study, with the main aim to identify existing drugs with the highest 3CLpro inhibition power. Among antiviral agents, lopinavir, idoxuridine, paritaprevir, and favipiravir showed the highest inhibition potential. Enzyme - ligand interactions as a key ingredient for successful drug design.

COVID-19 and Diabetes: A Comprehensive Review of Angiotensin Converting Enzyme 2, Mutual Effects and Pharmacotherapy

The potential relationship between diabetes and COVID-19 has been evaluated. However, new knowledge is rapidly emerging. In this study, we systematically reviewed the relationship between viral cell surface receptors (ACE2, AXL, CD147, DC-SIGN, L-SIGN and DPP4) and SARS-CoV-2 infection risk, and emphasized the implications of ACE2 on SARS-CoV-2 infection and COVID-19 pathogenesis. Besides, we updated on the two-way interactions between diabetes and COVID-19, as well as the treatment options for COVID-19 comorbid patients from the perspective of ACE2. The efficacies of various clinical chemotherapeutic options, including anti-diabetic drugs, renin-angiotensin-aldosterone system inhibitors, lipid-lowering drugs, anticoagulants, and glucocorticoids for COVID-19 positive diabetic patients were discussed. Moreover, we reviewed the significance of two different forms of ACE2 (mACE2 and sACE2) and gender on COVID-19 susceptibility and severity. This review summarizes COVID-19 pathophysiology and the best strategies for clinical management of diabetes patients with COVID-19.

SYNBIP: synthetic binding proteins for research, diagnosis and therapy

The success of protein engineering and design has extensively expanded the protein space, which presents a promising strategy for creating next-generation proteins of diverse functions. Among these proteins, the synthetic binding proteins (SBPs) are smaller, more stable, less immunogenic, and better of tissue penetration than others, which make the SBP-related data attracting extensive interest from worldwide scientists. However, no database has been developed to systematically provide the valuable information of SBPs yet. In this study, a database named ‘Synthetic Binding Proteins for Research, Diagnosis, and Therapy (SYNBIP)’ was thus introduced. This database is unique in (a) comprehensively describing thousands of SBPs from the perspectives of scaffolds, biophysical & functional properties, etc.; (b) panoramically illustrating the binding targets & the broad application of each SBP and (c) enabling a similarity search against the sequences of all SBPs and their binding targets. Since SBP is a human-made protein that has not been found in nature, the discovery of novel SBPs relied heavily on experimental protein engineering and could be greatly facilitated by in-silico studies (such as AI and computational modeling). Thus, the data provided in SYNBIP could lay a solid foundation for the future development of novel SBPs. The SYNBIP is accessible without login requirement at both official (https://idrblab.org/synbip/) and mirror (http://synbip.idrblab.net/) sites.

Structure-Based Discovery of Novel Nonpeptide Inhibitors Targeting SARS-CoV-2 Mpro

The continual spread of novel coronavirus disease 2019 (COVID-19) is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), posing a severe threat to the health worldwide. The main protease (Mpro, alias 3CLpro) of SARS-CoV-2 is a crucial enzyme for the maturation of viral particles and is a very attractive target for designing drugs to treat COVID-19. Here, we propose a multiple conformation-based virtual screening strategy to discover inhibitors that can target SARS-CoV-2 Mpro. Based on this strategy, nine Mpro structures and a protein mimetics library with 8960 commercially available compounds were prepared to carry out ensemble docking for the first time. Five of the nine structures are apo forms presented in different conformations, whereas the other four structures are holo forms complexed with different ligands. The surface plasmon resonance assay revealed that 6 out of 49 compounds had the ability to bind to SARS-CoV-2 Mpro. The fluorescence resonance energy transfer experiment showed that the biochemical half-maximal inhibitory concentration (IC50) values of the six compounds could hamper Mpro activities ranged from 0.69 ± 0.05 to 2.05 ± 0.92 μM. Evaluation of antiviral activity using the cell-based assay indicated that two compounds (Z1244904919 and Z1759961356) could strongly inhibit the cytopathic effect and reduce replication of the living virus in Vero E6 cells with the half-maximal effective concentrations (EC50) of 4.98 ± 1.83 and 8.52 ± 0.92 μM, respectively. The mechanism of the action for the two inhibitors were further elucidated at the molecular level by molecular dynamics simulation and subsequent binding free energy analysis. As a result, the discovered noncovalent reversible inhibitors with novel scaffolds are promising antiviral drug candidates, which may be used to develop the treatment of COVID-19.

Machine Learning Models Identify Inhibitors of SARS-CoV-2

With the rapidly evolving SARS-CoV-2 variants of concern, there is an urgent need for the discovery of further treatments for the coronavirus disease (COVID-19). Drug repurposing is one of the most rapid strategies for addressing this need, and numerous compounds have already been selected for in vitro testing by several groups. These have led to a growing database of molecules with in vitro activity against the virus. Machine learning models can assist drug discovery through prediction of the best compounds based on previously published data. Herein, we have implemented several machine learning methods to develop predictive models from recent SARS-CoV-2 in vitro inhibition data and used them to prioritize additional FDA-approved compounds for in vitro testing selected from our in-house compound library. From the compounds predicted with a Bayesian machine learning model, lumefantrine, an antimalarial was selected for testing and showed limited antiviral activity in cell-based assays while demonstrating binding (K_d 259 nM) to the spike protein using microscale thermophoresis. Several other compounds which we prioritized have since been tested by others and were also found to be active in vitro. This combined machine learning and in vitro testing approach can be expanded to virtually screen available molecules with predicted activity against SARS-CoV-2 reference WIV04 strain and circulating variants of concern. In the process of this work, we have created multiple iterations of machine learning models that can be used as a prioritization tool for SARS-CoV-2 antiviral drug discovery programs. The very latest model for SARS-CoV-2 with over 500 compounds is now freely available at www.assaycentral.org.

Tenofovir, Another Inexpensive, Well-Known and Widely Available Old Drug Repurposed for SARS-COV-2 Infection

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection is spreading worldwide with different clinical manifestations. Age and comorbidities may explain severity in critical cases and people living with human immunodeficiency virus (HIV) might be at particularly high risk for severe progression. Nonetheless, current data, although sometimes contradictory, do not confirm higher morbidity, risk of more severe COVID-19 or higher mortality in HIV-infected people with complete access to antiretroviral therapy (ART). A possible protective role of ART has been hypothesized to explain these observations. Anti-viral drugs used to treat HIV infection have been repurposed for COVID-19 treatment; this is also based on previous studies on severe acute respiratory syndrome virus (SARS-CoV) and Middle East respiratory syndrome virus (MERS-CoV). Among them, lopinavir/ritonavir, an inhibitor of viral protease, was extensively used early in the pandemic but it was soon abandoned due to lack of effectiveness in clinical trials. However, remdesivir, a nucleotide analog that acts as reverse-transcriptase inhibitor, which was tested early during the pandemic because of its wide range of antiviral activity against several RNA viruses and its safety profile, is currently the only antiviral medication approved for COVID-19. Tenofovir, another nucleotide analog used extensively for HIV treatment and pre-exposure prophylaxis (PrEP), has also been hypothesized as effective in COVID-19. No data on tenofovir's efficacy in coronavirus infections other than COVID-19 are currently available, although information relating to SARS-CoV-2 infection is starting to come out. Here, we review the currently available evidence on tenofovir's efficacy against SARS-CoV-2.

Binding Characterization of Agonists and Antagonists by MCCS: A Case Study from Adenosine A2A Receptor

Characterizing the structural basis of ligand recognition of adenosine A_2A receptor (AA_2AR) will facilitate its rational design and development of small molecules with high affinity and selectivity, as well as optimal therapeutic effects for pain, cancers, drug abuse disorders, etc. In the present work, we applied our reported algorithm, molecular complex characterizing system (MCCS), to characterize the binding features of AA_2AR based on its reported 3D structures of protein-ligand complexes. First, we compared the binding score to the reported experimental binding affinities of each compound. Then, we calculated an output example of residue energy contribution using MCCS and compared the results with data obtained from MM/GBSA. The consistency in results indicated that MCCS is a powerful, fast, and accurate method. Sequentially, using a receptor-ligand data set of 57 crystallized structures of AA_2ARs, we characterized the binding features of the binding pockets in AA_2AR, summarized the key residues that distinguish antagonist from agonist, produced heatmaps of residue energy contribution for clustering various statuses of AA_2ARs, explored the selectivity between AA_2AR and AA_1AR, etc. All the information provided new insights into the protein features of AA_2AR and will facilitate its rational drug design.

IsAb: a computational protocol for antibody design

The design of therapeutic antibodies has attracted a large amount of attention over the years. Antibodies are widely used to treat many diseases due to their high efficiency and low risk of adverse events. However, the experimental methods of antibody design are time-consuming and expensive. Although computational antibody design techniques have had significant advances in the past years, there are still some challenges that need to be solved, such as the flexibility of antigen structure, the lack of antibody structural data and the absence of standard antibody design protocol. In the present work, we elaborated on an in silico antibody design protocol for users to easily perform computer-aided antibody design. First, the Rosetta web server will be applied to generate the 3D structure of query antibodies if there is no structural information available. Then, two-step docking will be used to identify the binding pose of an antibody-antigen complex when the binding information is unknown. ClusPro is the first method to be used to conduct the global docking, and SnugDock is applied for the local docking. Sequentially, based on the predicted binding poses, in silico alanine scanning will be used to predict the potential hotspots (or key residues). Finally, computational affinity maturation protocol will be used to modify the structure of antibodies to theoretically increase their affinity and stability, which will be further validated by the bioassays in the future. As a proof of concept, we redesigned antibody D44.1 and compared it with previously reported data in order to validate IsAb protocol. To further illustrate our proposed protocol, we used cemiplimab antibody, a PD-1 checkpoint inhibitor, as an example to showcase a step-by-step tutorial.