Добавил:
Sekretar
kiopkiopkiop18@yandex.ru
t.me/Prokururor I Вовсе не секретарь, но почту проверяю
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:Ординатура / Хирургия / Библиотека им академика М.И. Перельмана / Книга_5529_Библиотеки_им_академика_М_И_Перельмана
.pdf
144 Bioinformatics of Autoimmune Diseases
One key application is sequence alignment and homology analysis, where researchers compare
autoimmune-related proteins with sequences from other species or related proteins within humans.
By aligning these sequences using tools like BLAST or Clustal Omega, scientists can identify conserved regions and mutations that may contribute to disease susceptibility. For example, in SLE,
analyzing the sequence variations in HLA proteins can provide insights into disease-associated
genetic predispositions.
FASTA sequences are also critical in epitope prediction and vaccine development for autoimmune diseases. Using computational immunology tools such as IEDB (Immune Epitope Database)
and NetMHC, researchers can analyze protein sequences to identify antigenic regions that trigger
immune responses. This is particularly useful in diseases like type 1 diabetes, where autoantigens
like GAD65 and insulin peptide fragments can be mapped to understand how T cells recognize and
attack self-proteins.
In structural modeling and drug discovery, protein sequences extracted from FASTA les are
used to predict 3D structures of autoimmune-related proteins through computational techniques like
homology modeling (SWISS-MODEL) or molecular docking (AutoDock, AlphaFold). This helps in
designing small molecules or biologics that can inhibit or modulate the function of immune-related
proteins, offering potential treatments.
Another application is phylogenetic analysis, where researchers study evolutionary relationships
among autoimmune disease-associated proteins. By constructing phylogenetic trees from sequence
alignments, scientists can trace the evolution of immune system proteins and understand how
genetic variations might have contributed to disease susceptibility in different populations.
Additionally, FASTA sequences are used in machine learning and articial intelligence (AI)–
based biomarker discovery, where large-scale sequence data is analyzed to identify patterns that
correlate with disease progression, severity, or treatment response. For example, in myasthenia gravis, analyzing sequence variations in the acetylcholine receptor (AChR) and muscle-specic kinase
(MuSK) proteins can help develop predictive models for patient outcomes.
By providing a structured way to access protein sequences, this script enables researchers to
quickly retrieve, analyze, and utilize protein data in various computational and experimental workows for autoimmune disease research.
This Python script “fetch _ uniprot _ seqs _ by _ disease.py” retrieves protein
sequences associated with a specic disease from the UniProt database and saves them in FASTA
format. The user provides a valid disease name, and the script queries the UniProt REST API for
reviewed protein entries, extracting key information such as UniProt accession ID, protein name,
gene name, sequence length, and the amino acid sequence itself. Each protein entry is formatted
into the FASTA standard, ensuring compatibility with bioinformatics tools for sequence alignment,
homology analysis, and structural modeling.
The script efciently handles nested JSON structures to accurately extract sequence data.
A key improvement is the correction of the sequence retrieval process, ensuring that the script
accesses the correct eld within UniProt’s API response. By implementing robust error handling, it lters out incomplete records, preventing missing or malformed sequences from being
included. If a valid sequence is found, the script prints debug messages conrming the extracted
details, such as the UniProt ID and sequence length, providing transparency in the data retrieval
process.
The output (Figure 4.6) is a properly formatted FASTA le where each sequence entry begins
with a deine containing the UniProt ID, protein name, gene name, and sequence length, separated
by pipes (|). This structure allows for seamless integration into computational pipelines. An example
of the output includes protein entries such as CTLA4 and HLA class II histocompatibility antigen
(HLA-DRB1), essential components in immune system regulation. The amino acid sequences are
stored in a format that can be used directly in BLAST searches or multiple sequence alignment
tools. This script provides researchers with an automated, reliable, and structured way to fetch protein sequence data for disease-related studies.
145 Bioinformatics Databases
FIGURE 4.6 Protein FASTA with deines showing ID, name, and length.
4.2.2.3.2.4 Mapping a UniProt ID to PDB IDs Mapping a UniProt ID to PDB IDs is essential for understanding the structural aspects of proteins and their implications in various diseases,
including autoimmune disorders. The UniProt database provides comprehensive information about
proteins, including their sequence, function, and structural data. To retrieve PDB IDs associated with
a UniProt entry, one can use the UniProt API or cross-references in the UniProt web interface. The
PDB contains experimentally determined 3D structures of proteins, which are critical for studying
protein interactions, mutations, and drug design. In the context of autoimmune diseases, structural
insights from PDB can help researchers understand how misfolded proteins, antigen presentation, or
autoantibody interactions contribute to disease pathology. For instance, in diseases like rheumatoid
arthritis or multiple sclerosis, analyzing PDB structures of immune system proteins, such as MHC
molecules, can reveal how genetic variations impact antigen binding and immune response. This
knowledge facilitates the design of targeted therapies, monoclonal antibodies, and small-molecule
inhibitors to modulate immune activity and treat autoimmune conditions effectively.
The following Python script “fetch _ u niprot _ PDB.p y ” retrieves PDB IDs associated
with a given UniProt ID by querying the UniProt API. It constructs a request URL using the provided UniProt ID and sends an HTTP GET request to fetch the corresponding protein data in JSON
format. If the response is successful, the script parses the JSON output to extract relevant structural
data. Specically, it searches for cross-references to the PDB database, which stores experimentally
determined 3D protein structures. These cross-references are listed in the “uniProtKBCross-
References” section of the JSON response, where each entry species the external database and
the corresponding ID. The script iterates through these entries, identifying those linked to PDB, and
collects the corresponding PDB IDs.
Once the PDB IDs are extracted, the script formats them into a structured output and saves them
to a text le. The le is named after the UniProt ID, ensuring easy identication and organization of results. If no PDB entries are found for a given UniProt ID, the script provides appropriate
feedback, indicating that no structural data is available. This automated process simplies the task
of retrieving protein structure information, which is crucial for researchers studying protein function, interactions, and mutations. By leveraging the UniProt API, the script enables rapid access to
protein structural data, which can be further analyzed using molecular modeling tools or structural
bioinformatics techniques.
import requests
def map_uniprot_to_pdb(uniprot_id):
url = f"https://rest.uniprot.org/uniprotkb/{uniprot_id}.json"
response = requests.get(url)
if response.status_code != 200:
print(f"Error: Unable to retrieve data for UniProt ID
{uniprot_id}")
return None
146 Bioinformatics of Autoimmune Diseases
data = response.json()
# Extract PDB information if available
pdb_entries = []
if 'uniProtKBCrossReferences' in data:
for entry in data['uniProtKBCrossReferences']:
if entry['database'] == 'PDB':
pdb_entries.append(entry['id'])
if pdb_entries:
pdb_list = '\n'.join(pdb_entries)
filename = f"{uniprot_id}.txt"
with open(filename, "w") as file:
file.write(f"PDB IDs for UniProt ID {uniprot_id}:\n{pdb_list}")
print(f"PDB IDs saved to {filename}")
else:
print(f"No PDB entries found for UniProt ID {uniprot_id}")
uniprot_id = "P69905" # Example: Hemoglobin subunit alpha (HBA1)
map_uniprot_to_pdb(uniprot_id)
4.2.2.3.2.5 Mapping UniProt Proteins to NCBI Coding Regions Mapping a UniProt protein
to its corresponding CDS in the NCBI RefSeq nucleotide database is a crucial step in understanding the genetic and molecular mechanisms underlying autoimmune diseases and in leveraging deep
learning for PTM site prediction. Autoimmune disorders often result from genetic mutations or
dysregulation in immune-related proteins. By linking a protein sequence to its CDS, researchers can
analyze genetic variations, alternative splicing events, and regulatory elements that may contribute
to disease susceptibility. This mapping enables the identication of pathogenic mutations, facilitates
personalized medicine strategies, and aids in the development of targeted therapies such as gene
editing and RNA-based treatments.
In addition to disease research, CDS data is essential for deep learning-based prediction of PTM
sites, which play a critical role in protein function and immune signaling. PTMs, such as phosphorylation, glycosylation, and ubiquitination, can alter protein activity and are often implicated in autoimmune disorders. By training machine learning models on both genomic (CDS) and proteomic
(PTM) data, researchers can uncover novel disease-associated PTM patterns, improve biomarker
discovery, and rene drug target validation. Recent work, such as CaLMPhosKAN, demonstrates
how codon-aware and amino acid–aware embeddings, fused with advanced neural architectures,
can effectively predict phosphorylation sites using CDS-derived features (Pratyush et al., 2025).
The integration of genomic, proteomic, and AI-driven approaches enhances our understanding of
how genetic variations inuence protein modications, ultimately leading to more precise diagnostics and treatment strategies for autoimmune diseases.
The provided code “fetc h _ u niprot _ c ds.p y ” enables the mapping of a UniProt ID to
its corresponding RefSeq nucleotide sequence (NM_ accession) in the NCBI database and extracts
the CDS. This process is crucial in linking protein sequences to their genetic origins, which can
provide valuable insights into genetic mutations, disease mechanisms, and protein function. The
script begins by retrieving the gene name associated with a given UniProt ID using UniProt’s REST
API. Since UniProt contains extensive annotations for proteins, including their associated genes, the
function rst makes an API request and extracts the gene name from the returned JSON response.
This step is essential because NCBI’s nucleotide database organizes genetic information by gene
name rather than protein ID.
Once the gene name is obtained, the script proceeds to search the NCBI Nucleotide database
for its RefSeq mRNA accession (NM_), which represents the transcribed mRNA from which
the protein is translated. To achieve this, an Entrez search query is constructed using the gene
name and the organism, ensuring that only sequences relevant to H. sapiens (or another specied
species) are retrieved. The search specically lters results to include only RefSeq sequences,
147 Bioinformatics Databases
ensuring that manually curated and high-condence reference sequences are used. The script
fetches up to ve results and checks whether the returned records contain an NM_ accession,
which is necessary for further analysis. The rst valid NM_ sequence found is then selected for
downstream processing.
After identifying the correct RefSeq nucleotide sequence, the script retrieves its GenBank record
using the Entrez Efetch API and parses the retrieved le to extract the CDS. In GenBank annotations, the CDS region is specically marked within the feature table, indicating the exact segment
of the nucleotide sequence that corresponds to the protein-coding region. By iterating through the
feature table and extracting sequences labeled “CDS”, the script isolates the precise nucleotide
sequence that encodes the protein. This is crucial for downstream applications, such as studying
genetic variations, alternative splicing, or deep learning-based predictions of PTMs.
The nal step of the script ensures that the extracted CDS sequences are stored in a FASTA
format for easy reference and computational analysis. The FASTA le name is dynamically generated based on the gene name, ensuring clarity and organization of results. The header (deine) of
the FASTA sequence contains essential information, including the UniProt ID, gene name, RefSeq
accession, and CDS number, making it easy to trace the sequence back to its original database
sources. The script utilizes Biopython’s SeqIO module to format and write the extracted CDS
sequences into a FASTA le. This structured approach ensures that researchers can readily use the
data for further genomic, proteomic, or AI-based studies, enhancing the ability to analyze genetic
inuences on protein function and disease pathology.
4.2.3 KEGG PATHWAY DATABASE
The KEGG Pathway database is a comprehensive resource that systematically organizes information on molecular interaction networks, cellular processes, and human diseases. Developed as part
of the KEGG, it serves as a platform for understanding high-level biological functions and utilities
derived from genomic and molecular datasets. KEGG pathway maps provide graphical representations of biochemical pathways, including metabolism, genetic and environmental information
processing, and various cellular processes. By integrating genes, proteins, metabolites, and other
biological entities into a unied framework, the KEGG Pathway database facilitates the study of
complex biological systems and their functional interactions.
The database is structured into several pathway categories, each covering a different aspect of
biological function. Metabolic pathways form a signicant portion, detailing essential biochemical
reactions such as glycolysis, the citric acid cycle, and amino acid metabolism. These maps illustrate
the ow of metabolites through enzymatic reactions, linking genes and enzymes to their respective reactions and allowing researchers to explore metabolic regulation and adaptations in various
organisms. In addition to metabolism, KEGG also includes pathways related to genetic information processing, which encompasses transcription, translation, replication, and repair mechanisms.
These pathways highlight the coordination of molecular machinery in the synthesis and maintenance of genetic material, revealing how genes encode functional proteins and how cellular processes ensure genomic stability.
Another crucial component of KEGG pathways is its representation of signaling and regulatory networks that govern cellular functions. Signal transduction pathways, such as the MAPK,
PI3K-Akt, and JAK–STAT signaling cascades, illustrate how cells respond to external stimuli
and transmit signals through a cascade of molecular interactions. These pathways are essential
for understanding cellular communication, differentiation, apoptosis, and immune responses. The
database also includes regulatory pathways such as cell cycle control and apoptosis, which help
researchers investigate mechanisms underlying cell proliferation and programmed cell death. By
providing detailed pathway maps, KEGG enables the study of how molecular signals are relayed
and modulated to maintain cellular homeostasis.
148 Bioinformatics of Autoimmune Diseases
One of the most signicant applications of the KEGG Pathway database is in the study of human
diseases. Disease pathways in KEGG link genetic mutations and molecular dysfunctions to pathological conditions, offering insights into the mechanisms of diseases such as cancer, neurodegenerative disorders, cardiovascular diseases, and autoimmune conditions. For example, rheumatoid
arthritis is represented as a pathway that maps the interactions of cytokines, immune cells, and
signaling molecules involved in inammation and joint destruction. Cancer-related pathways integrate genetic alterations with disrupted signaling networks, illustrating how mutations in oncogenes
and tumor suppressor genes drive malignant transformation. By associating genetic variants with
disease mechanisms, KEGG supports translational research in developing therapeutic interventions
and personalized medicine.
KEGG pathways also facilitate the integration of molecular and drug interaction data, making
them valuable for pharmacological studies. The database includes pharmacokinetics and pharmacodynamics pathways that describe drug metabolism, absorption, and excretion, allowing researchers
to explore drug–target interactions. By linking small molecules to their corresponding metabolic
or signaling pathways, KEGG helps in identifying potential drug candidates and understanding
adverse drug reactions. This integration of pathway data with chemical and genomic information
enhances drug discovery efforts and the development of targeted therapies for various diseases.
As a highly interactive and continuously updated resource, KEGG pathway maps support computational biology and bioinformatics applications. Researchers can access KEGG pathways through the
KEGG REST API, enabling programmatic retrieval of pathway data, genes, enzymes, and molecular
interactions. The database is widely used for pathway enrichment analysis, helping researchers identify
signicantly enriched pathways in gene expression studies. Through KEGG’s integration with omics
data, such as transcriptomics and proteomics, scientists can interpret large-scale biological datasets
in the context of well-dened pathways. This capability makes KEGG an essential tool for systems
biology, enabling a holistic understanding of biological networks and their regulatory mechanisms.
By providing a structured framework for understanding biochemical and regulatory interactions,
the KEGG Pathway database serves as a cornerstone for research in genomics, molecular biology,
and medicine. Its extensive collection of manually curated and computationally inferred pathways
offers a knowledge base that supports diverse scientic inquiries, from basic biological research to
applied biomedical studies. Through its dynamic representation of molecular interactions and functional hierarchies, KEGG continues to be an indispensable resource for elucidating the complexities
of life at the molecular level.
KEGG human disease pathways provide a structured framework for understanding the molecular mechanisms underlying autoimmune diseases. By mapping genes, proteins, and signaling molecules involved in immune dysregulation, these pathways help researchers uncover how genetic
and environmental factors contribute to disease progression. Autoimmune diseases result from the
immune system attacking the body’s own tissues. KEGG pathways allow scientists to visualize the
key molecular interactions that drive inammation, autoantibody production, and tissue damage,
providing critical insights for diagnostics and therapeutic development.
The utility of KEGG in autoimmune disease research extends to the identication of key signaling pathways implicated in immune dysfunction. For example, the KEGG rheumatoid arthritis
pathway highlights the involvement of cytokines such as tumor necrosis factor (TNF), interleukins (IL-6, IL-17), and immune-regulatory molecules that mediate chronic inammation and joint
destruction. By integrating genomic and transcriptomic data into these pathways, researchers have
been able to identify novel therapeutic targets. This has contributed to the success of biologic drugs
such as TNF inhibitors, IL-6 blockers, and JAK inhibitors, which have signicantly improved treatment outcomes for rheumatoid arthritis patients.
KEGG pathways also facilitate the study of shared molecular mechanisms across different
autoimmune diseases, revealing common immunological signatures. The Toll-like receptor signaling pathway, for example, is a major component in autoimmune disease pathogenesis, playing
a key role in innate immune activation. Dysregulation of this pathway has been linked to SLE
149 Bioinformatics Databases
and inammatory bowel disease, where excessive immune activation leads to chronic inammation. KEGG pathway–based studies have enabled the discovery of targeted therapies, such as TLR
inhibitors, which have shown promise in preclinical and clinical trials for autoimmune conditions.
One of the major successes of KEGG in autoimmune disease research has been its role in
advancing personalized medicine approaches. By integrating pathway analysis with patient-specic
genetic and transcriptomic data, researchers have been able to stratify patients based on their molecular proles. This has led to precision medicine initiatives in diseases like multiple sclerosis, where
pathway-based biomarker discovery has helped predict treatment responses. For example, studies
using KEGG pathways have identied differential expression of immune-related genes that correlate with response to interferon-beta therapy, allowing clinicians to personalize treatment regimens
for better outcomes.
The use of KEGG human disease pathways has also accelerated drug repurposing efforts for
autoimmune diseases. By analyzing pathway interactions between existing drugs and diseaserelated targets, researchers have successfully identied potential off-label uses of medications. One
such example is the discovery that Janus kinase (JAK) inhibitors, originally developed for myelobrosis, can effectively modulate immune pathways involved in rheumatoid arthritis and ulcerative
colitis. KEGG pathway analysis has played a crucial role in these discoveries by highlighting the
molecular overlaps between diseases, facilitating a more efcient approach to drug repurposing.
Beyond drug development, KEGG pathways have been instrumental in uncovering environmental and microbial inuences on autoimmune disease pathogenesis. The integration of KEGG
metabolic pathways with gut microbiome data has revealed that microbial metabolites can regulate
immune function and contribute to diseases such as type 1 diabetes and multiple sclerosis. By linking dysbiosis to specic metabolic and immune pathways, KEGG has provided a valuable resource
for understanding how environmental factors interact with host immunity. This has led to emerging
strategies that explore microbiome-targeted therapies, such as probiotics and dietary interventions,
for autoimmune disease management.
The success of KEGG human disease pathways in autoimmune disease research highlights its
role as an indispensable tool for unraveling complex immune mechanisms and guiding therapeutic
innovations. By continuously updating and rening pathway models, KEGG ensures that researchers have access to the latest insights into disease-associated molecular networks. The integration
of KEGG pathways with high-throughput omics data, AI, and network-based drug discovery is
expected to further enhance our ability to treat and prevent autoimmune diseases in the future.
4.2.3.1 Retrieving KEGG Disease and Pathway List
Retrieving the KEGG disease list, which includes both the KEGG disease ID and disease name,
enhances database searching by providing a controlled vocabulary that ensures accuracy and consistency in queries. By using the exact KEGG disease name or its corresponding ID, researchers
can avoid issues related to synonyms, variations in terminology, or ambiguous disease classications. This controlled approach streamlines data retrieval, improves reproducibility, and facilitates
cross-referencing with other KEGG resources such as pathway maps, gene associations, and drug
interactions. Additionally, it allows for automated and large-scale searches, making it easier to analyze disease-related biological networks and conduct comparative studies across different datasets.
We may need to install bioservices, which is a Python package designed to interact with various biological databases, including KEGG, UniProt, Ensembl, NCBI, and more. It allows researchers to retrieve, analyze, and manipulate biological data programmatically.
from bioservices import KEGG
def fetch_Kegg_dis_list():
kegg = KEGG()
disease_list = kegg.list("disease")
with open("kegg_diseases.txt", "w") as f:
150 Bioinformatics of Autoimmune Diseases
f.write(disease_list)
print("Disease list saved to kegg_diseases.txt")
fetch_Kegg_dis_list()
The above Python script “kegg _ dis _ list.py” uses the bioservices library to retrieve
the list of diseases from the KEGG database and saves the information to a text le named kegg _
diseases.txt (Figure 4.7).
It initializes a connection to KEGG using the K EG G() class and then fetches the disease list by
calling kegg.list("disease"), which returns the data in a tab-separated format containing
KEGG disease IDs and their corresponding disease names. The script then writes this raw data
directly to a le without additional processing or formatting. Once the le is created, a message is
printed to conrm that the disease list has been successfully saved. This script provides a simple
and direct way to store KEGG disease information locally for reference or further processing.
You can also search the KEGG database for a specic category of diseases using a keyword.
In the following Python script “kegg _ find _ dis.py”, we search for diseases in the KEGG
database whose names include the word “autoimmune”.
from bioservices import KEGG
def find_kegg_dis(keyword):
kegg = KEGG()
diseases = kegg.find("disease", keyword)
print(diseases)
find_kegg_dis("autoimmune")
The following Python script “kegg _ pathway _ list.py” retrieves KEGG pathway list
and save them in a text le:
from bioservices import KEGG
def kegg_pathway_list():
kegg = KEGG()
pathways = kegg.list("pathway")
pathways_list = [line for line in pathways.split("\n")]
with open("kegg_pathways.txt", "w") as f:
f.write("\n".join(pathways_list))
FIGURE 4.7 KEGG disease IDs and their corresponding disease names.
151 Bioinformatics Databases
print("\n".join(pathways_list))
print("Pathway list saved to kegg_pathways.txt")
kegg_pathway_list()
4.2.3.2 Retrieving Autoimmune Disease Pathway
A KEGG pathway record provides detailed information about a specic biological pathway, including its description, involved genes, enzymes, compounds, and interactions. Each pathway is assigned
a unique ID (e.g., hsa05323 for rheumatoid arthritis) and includes annotations about molecular networks, regulatory relationships, and disease associations.
The following Python script “kegg _ pathway _ id.py” uses the bioservices package
to retrieve pathway information from KEGG based on a given pathway ID. It initializes an instance
of the KEGG class, sends a request to fetch the pathway record, and prints the retrieved data.
from bioservices import KEGG
def fetch_pathway_by_id(pathway_id):
kegg = KEGG()
pathway_info = kegg.get(pathway_id)
print(pathway_info)
pathway_id = "hsa05323" #Example pathway for rheumatoid arthritis
fetch_pathway_by_id(pathway_id)
The Python script “kegg _ pathway _ dsease.py” is designed to retrieve and organize
information from the KEGG database regarding human disease pathways, with a focus on autoimmune diseases. It systematically extracts the KEGG disease ID and description, retrieves the
associated pathways, and identies genes involved in those pathways. The script is structured into
multiple steps, beginning with querying the KEGG API to nd the disease of interest. By using
the /find/disease/{disease _ name} endpoint, it identies the KEGG-specic disease ID,
which serves as a reference for subsequent data retrieval. The next step involves obtaining KEGG
pathway IDs linked to the disease using /link/pathway/{disease _ id} API, which returns
pathways associated with the identied disease.
Once the pathways are retrieved, the script proceeds to extract genes involved in these pathways.
Instead of using a direct disease-to-gene approach, which is not supported by KEGG, it links genes
to pathways using the /link/hsa/{pathway _ id} API. This ensures that only human genes are
retrieved and prevents the inclusion of genes from other species. The script stores gene IDs in a set to
eliminate duplicates before proceeding to extract the ofcial gene names. The /get/hsa:{gene _
id} API is used for each gene to obtain its corresponding name, ensuring that the results are meaningful and human-readable. To handle missing data, the script includes an enhanced regular expression
that searches not only for the NAME eld but also for the DESCRIPTION eld if the name is not found.
To provide a structured and accessible output, the script saves the extracted data in two different
formats. The rst output is a JSON le, named after the disease, which contains the KEGG disease
information, pathways, and gene details in a structured format. This le is useful for computational
processing and further data analysis. The second output is a formatted text le, which organizes the
extracted data in a human-readable manner. Instead of listing multiple entries as separate lines, the
text le consolidates pathways and gene data using a pipe (|) separator. The disease ID and description appear at the top, followed by pathway IDs and a list of genes, where each gene entry includes
both the gene ID and its corresponding name.
The output of the script ensures clarity and usability for both computational and manual analysis.
The JSON le provides a structured dataset that can be easily parsed for further processing, while
the text le offers a compact summary of key information. For example, if the script is executed
for rheumatoid arthritis, the JSON le will contain the KEGG disease ID H00019, the description
“Rheumatoid arthritis”, and a pathway list that includes hsa05323, the primary KEGG pathway
for rheumatoid arthritis. The gene list will include key inammatory genes such as TNF, IL6, and
152 Bioinformatics of Autoimmune Diseases
CTLA4, ensuring that researchers and clinicians can quickly identify relevant molecular players in
the disease. The corresponding text le will present this information in a condensed form, listing
the disease details, pathways, and genes in a way that is easy to interpret.
The Python script efciently automates data retrieval from KEGG, avoiding manual searches
and data entry. The integration of rate-limiting measures through ti m e.slee p(1) ensures compliance with KEGG’s API restrictions while preventing potential request failures due to excessive
calls. Additionally, robust error handling prevents the script from crashing if KEGG data is incomplete or if an API request fails. By structuring the output into both JSON and text formats, the script
provides exibility for different types of users, whether they need structured data for computational
analysis or a summarized view for direct interpretation. Overall, the script serves as a powerful
tool for researchers studying autoimmune diseases, allowing them to explore disease mechanisms,
molecular interactions, and potential therapeutic targets in a systematic and automated manner.
4.2.3.3 Automated KEGG Pathway ID Retrieval for Autoimmune Diseases
The Python script “kegg _ multiple _ diseases.py” is designed to automate the retrieval of
KEGG pathway IDs for multiple autoimmune diseases, ensuring that key autoimmune diseases are
included in the dataset. The script takes a more reliable approach by directly querying the KEGG
database for all diseases and ltering those that match a predened list of known autoimmune conditions. This method eliminates the risk of missing relevant diseases due to inconsistencies in KEGG’s
categorization and ensures that only well-established autoimmune diseases are included in the output.
The script works by rst fetching all disease entries from KEGG using the /list/disease end-
point. It then cross-references these results with a predened list of autoimmune diseases to identify
their corresponding KEGG disease IDs. By matching known disease names exactly, the script avoids
missing important conditions due to minor variations in naming conventions. Once the relevant disease IDs are collected, they are stored in structured formats for further analysis. The output includes
a JSON le that provides a structured dataset for computational research and a human-readable text
le for easy reference. Each output le contains the KEGG disease ID along with the full disease
name, allowing researchers to quickly identify and link diseases to their corresponding pathways.
The JSON output can be seamlessly integrated into bioinformatics pipelines, enabling researchers to programmatically retrieve disease-related pathways, genes, and molecular interactions. By
using this dataset, scientists studying autoimmune disorders can explore pathway-based disease
mechanisms, identify potential drug targets, and conduct comparative analyses across multiple autoimmune conditions. The text le provides a simplied reference that can be used for manual inspection, making it useful for researchers who need a quick overview of KEGG disease associations.
This script is particularly useful in systems biology and biomedical research, where pathwaylevel analysis is essential for understanding the genetic and molecular basis of autoimmune diseases. It can be used in combination with pathway enrichment tools to determine which biological
pathways are overrepresented in specic disease conditions. Additionally, it can assist in identifying
common pathways among different autoimmune diseases, providing insights into shared immune
dysregulation mechanisms. Researchers studying drug repurposing can leverage the KEGG pathway data to nd potential therapeutic targets that overlap between diseases, ultimately contributing
to more precise treatment strategies.
By automating the retrieval of KEGG pathway IDs for autoimmune diseases, this script streamlines the process of accessing high-quality biological data. It eliminates the need for manual
searches, reduces errors in disease classication, and ensures consistency across studies. This
makes it a valuable tool for computational biologists, bioinformaticians, and medical researchers
investigating the molecular pathways underlying autoimmune disorders.
4.2.3.4 Retrieving Autoimmune Disease Drugs from the KEGG Drug Database
The KEGG Drug database is a comprehensive resource that provides information on approved and
experimental drugs, their chemical properties, molecular interactions, and therapeutic targets. It
integrates data from various sources to offer a detailed view of drug–target relationships, making
153 Bioinformatics Databases
it a valuable tool for researchers in pharmacology and biomedical sciences. The database is structured to link drugs with diseases, pathways, and molecular targets, enabling users to analyze drug
mechanisms and interactions. By leveraging KEGG Drug, researchers can explore how different
drugs interact with specic proteins, enzymes, and receptors, which is essential for drug repurposing, precision medicine, and understanding the pharmacodynamics of therapeutic agents.
In the context of autoimmune diseases, KEGG Drug can be used to identify drugs that modulate immune system activity, either by inhibiting pro-inammatory pathways or enhancing regulatory mechanisms. Autoimmune diseases often involve complex molecular interactions that make
drug discovery and treatment optimization challenging. By accessing KEGG Drug, researchers can
retrieve information on existing medications and their targets, assess potential new treatment strategies, and compare drugs based on their mechanisms of action. This information is particularly valuable for developing personalized treatment plans and for identifying alternative drugs for patients
who may not respond well to conventional therapies.
The provided Python script “kegg _ drugs _ disease.py” automates the retrieval of drug–
target interactions for a given autoimmune disease using KEGG’s REST API. It rst searches for
disease-related entries, extracts drug information linked to the disease, and then identies the molecular targets of these drugs. The script processes the data in a structured manner, ensuring that target
hierarchies are preserved to maintain clarity. The results are saved in both JSON and tab-separated
values (TSV) format, making them accessible for further analysis. The JSON format allows for easy
integration with other bioinformatics tools, while the TSV format provides a structured, humanreadable table suitable for data visualization and statistical analysis. Researchers can use these outputs
to study drug mechanisms, compare treatment options, and identify potential candidates for further
experimental validation, ultimately facilitating the advancement of autoimmune disease research.
In the output JSON le as shown in Figure 4.8, the drug target represents the biological molecules that a drug interacts with to produce its therapeutic effect. These targets are typically proteins,
FIGURE 4.8 KEGG drug–target interactions for autoimmune disease.
Соседние файлы в папке Библиотека им академика М.И. Перельмана